Protein flexibility from discrete molecular dynamics simulations using quasi‐physical potentials

RNA molecules with novel functions have revived interest in the accurate prediction of RNA three-dimensional (3D) structure and folding dynamics. However, existing methods are inefficient in automated 3D structure prediction. Here, we report a robust computational approach for rapid folding of RNA molecules. We develop a simplified RNA model for discrete molecular dynamics (DMD) simulations, incorporating base-pairing and base-stacking interactions. We demonstrate correct folding of 150 structurally diverse RNA sequences. The majority of DMD-predicted 3D structures have <4 A deviations from experimental structures. The secondary structures corresponding to the predicted 3D structures consist of 94% native base-pair interactions. Folding thermodynamics and kinetics of tRNA(Phe), pseudoknots, and mRNA fragments in DMD simulations are in agreement with previous experimental findings. Folding of RNA molecules features transient, non-native conformations, suggesting non-hierarchical RNA folding. Our method allows rapid conformational sampling of RNA folding, with computational time increasing linearly with RNA length. We envision this approach as a promising tool for RNA structural and functional analyses.

α-Synuclein (α-syn) is the major component of the intraneuronal inclusions called Lewy bodies, which are the pathological hallmark of Parkinson's disease. α-Syn is capable of self-assembly into many different species, such as soluble oligomers and fibrils. Even though attempts to resolve the structures of the protein have been made, detailed understanding about the structures and their relationship with the different aggregation steps is lacking, which is of interest to provide insights into the pathogenic mechanism of Parkinson's disease. Here we report the structural flexibility of α-syn monomers and dimers in an aqueous solution environment as probed by single-molecule time-lapse high-speed AFM. In addition, we present the molecular basis for the structural transitions using discrete molecular dynamics (DMD) simulations. α-Syn monomers assume a globular conformation, which is capable of forming tail-like protrusions over dozens of seconds. Importantly, a globular monomer can adopt fully extended conformations. Dimers, on the other hand, are less dynamic and show a dumbbell conformation that experiences morphological changes over time. DMD simulations revealed that the α-syn monomer consists of several tightly packed small helices. The tail-like protrusions are also helical with a small β-sheet, acting as a "hinge". Monomers within dimers have a large interfacial interaction area and are stabilized by interactions in the non-amyloid central (NAC) regions. Furthermore, the dimer NAC-region of each α-syn monomer forms a β-rich segment. Moreover, NAC-regions are located in the hydrophobic core of the dimer.

This article reviews heteronuclear NMR techniques that are used to characterize the 3D structures and internal molecular dynamics of proteins in solution. Applications to several proteins that have been studied in our laboratory are presented. Topics discussed include sequential assignments, the use of nuclear Overhauser enhancement spectroscopy (NOESY) and residual dipolar coupling restraints to determine protein structure, and the use of relaxation measurements to obtain information about protein internal motions. Also discussed are recent computational approaches directed at increasing the efficiency of structure determination, as well as new methods for studying large proteins.

A transition in hydraulically induced granular displacement patterns is studied by means of discrete numerical molecular dynamics simulations. During this transition the patterns change from fractures and fingers to finely dispersed bubbles. The dynamics of the displacement patterns are studied in a rectangular Hele-Shaw cell filled with a dense but permeable two-dimensional granular layer. At one side of the cell the pressure of the compressible interstitial gas is increased. At the opposite side from the inlet of the cell a semipermeable boundary is located. This boundary is only permeable towards the gas phase while preventing grains from leaving the cell. The imposed pressure gradient compacts the grains. In the process we can identify and describe a mechanism that controls the transition of the emerging displacement patterns from fractures and fingers to finely dispersed bubbles as a function of the interstitial gas's properties and the characteristics of the granular phase.

Computational insights into the aggregation mechanism of human calcitonin.

The abnormal aggregation of human calcitonin (hCT) hormone peptides impairs their physiological function, leading to harmful immune responses and cytotoxicity, which limits their clinical utility. Interestingly, a representative hCT analog incorporating Y12L and N17H substitutions (DM-hCT) has shown reduced aggregation tendencies while maintaining bioactivity. But the molecular mechanism of Y12L and N17H substitutions on the conformational dynamics of hCT remains unclear. Here, we systematically investigated the folding and self-assembly dynamics of hCT and DM-hCT using atomistic discrete molecular dynamics (DMD) simulations. Our findings revealed that hCT monomers predominantly adopted unstructured conformations with dynamic helices. Oligomerization of hCT resulted in the formation of β-sheet-rich aggregates and β-barrel intermediates. The Y12L and N17H substitutions enhanced helical conformations and suppressed β-sheet formation in both monomers and oligomers. These substitutions stabilized the dynamic helices and disrupted aromatic interactions responsible for β-sheet formation at residue 12. Notably, DM-hCT assemblies still exhibited β-sheets in phenylalanine-rich and C-terminal hydrophobic regions, suggesting that future optimizations should focus on these areas. Our simulations provide insights into the molecular mechanisms underlying hCT aggregation and the amyloid-resistant effects of Y12L and N17H substitutions. These findings have valuable implications for the development of clinical hCT analogs.

Large-scale domain dynamics in proteins are found when flexible linkers or hinges connect domains. The related conformational changes are often related to the function of the protein, for example by arranging the active center after substrate binding or allowing transport and release of products. The adaptation of a specific active structure is referred to as ‘induced fit’ and is challenged by models such as ‘conformational sampling’. Newer models about protein function include some flexibility within the protein structure or even internal dynamics of the protein. As larger domains contribute to the configurational changes, the timescale of the involved motions is slowed down. The role of slow domain dynamics is being increasingly recognized as essential to understanding the function of proteins. Neutron spin echo spectroscopy (NSE) is a technique that is able to access the related timescales from 0.1 up to several hundred nanoseconds and simultaneously covers the length scale relevant for protein domain movements of several nanometers distance between domains. Here we focus on these large-scale domain fluctuations and show how the structure and dynamics of proteins can be assessed by small-angle neutron scattering and NSE.

Multiscale Modeling of Nucleosome Dynamics

Experimental findings suggest that oligomeric forms of the amyloid beta protein (Abeta) play a critical role in Alzheimer's disease. Thus, elucidating their structure and the mechanisms of their formation is critical for developing therapeutic agents. We use discrete molecular dynamics simulations and a four-bead protein model to study oligomerization of two predominant alloforms, Abeta40 and Abeta42, at the atomic level. The four-bead model incorporates backbone hydrogen-bond interactions and amino acid-specific interactions mediated through hydrophobic and hydrophilic elements of the side chains. During the simulations we observe monomer folding and aggregation of monomers into oligomers of variable sizes. Abeta40 forms significantly more dimers than Abeta42, whereas pentamers are significantly more abundant in Abeta42 relative to Abeta40. Structure analysis reveals a turn centered at Gly-37-Gly-38 that is present in a folded Abeta42 monomer but not in a folded Abeta40 monomer and is associated with the first contacts that form during monomer folding. Our results suggest that this turn plays an important role in Abeta42 pentamer formation. Abeta pentamers have a globular structure comprising hydrophobic residues within the pentamer's core and hydrophilic N-terminal residues at the surface of the pentamer. The N termini of Abeta40 pentamers are more spatially restricted than Abeta42 pentamers. Abeta40 pentamers form a beta-strand structure involving Ala-2-Phe-4, which is absent in Abeta42 pentamers. These structural differences imply a different degree of hydrophobic core exposure between pentamers of the two alloforms, with the hydrophobic core of the Abeta42 pentamer being more exposed and thus more prone to form larger oligomers.

The mechanical behavior of nanoparticle assemblies depends strongly on particle size, yet the underlying mechanisms remain insufficiently understood. In present study, we employ a scheme combining discrete element method (DEM) and molecular dynamics (MD) simulations to examine size-dependent strength and deformation in disordered copper nanoparticle assemblies. Granular packings generated by DEM were transformed into atomic models and subjected to uniaxial compression in MD simulations. Assemblies composed of nanoparticles with radius smaller than ~2.5 nm fully densify during relaxation, forming nanopolycrystalline solids, whereas larger particles preserve porous architectures. This structural divergence governs subsequent deformation. Small-particle assemblies deform through grain boundary migration and grain growth, exhibiting an inverse Hall–Petch-type strength dependence. In contrast, large-particle assemblies deform primarily via interparticle contact evolution and densification, with strength conforming to a Gibson–Ashby-type prediction. A scaling law captures the strength variation across size range in this regime. These results establish the competition between surface energy-driven densification and contact-dominated deformation as the controlling factor in the mechanical response of nanoparticle assemblies, providing guidance for designing nanoparticle-based materials with tailored mechanical performance.

Equilibrium constants governing certain reversible changes in structurally simple proteins are formulated by the method of partition functions, adopting as a model for the protein molecule a regular cylinder, capable of cleavage, partial or complete, about a plane parallel to the base. The types of equilibria considered are: (1) the complete separation into fragments, (2) the conversion of the rigid cylinder into one, the two parts of which are capable of independent rotation about the principal axis, and (3) the isomerism which results when the internal rotator couples firmly in two different positions. By means of the method of imaginary frequencies, expressions are derived for the corresponding velocity constants.

Because of unparalleled sensitivity and tolerance to protein size, mass spectrometry (MS) has become a popular method for measuring the solution hydrogen (1H/2H) exchange (HX) of biologically relevant protein states. While incorporated deuterium can be localized to different regions by pepsin proteolysis of the labeled protein, the assignment of deuteriums to individual residues is typically not obtained, thereby limiting a detailed understanding of HX and the dynamics of protein structure. Here we use gas-phase fragmentation of peptic peptides by electron transfer dissociation (ETD) to measure the HX of individual amide linkages in the amyloidogenic protein beta2-microglobulin. A comparison of the deuterium levels of 60 individual backbone amides of beta2-microglobulin measured by HX-ETD-MS analysis to the corresponding values measured by NMR spectroscopy shows an excellent correlation. The deuterium labeling pattern of beta2-microglobulin is retained in the gaseous fragment ions by employing mild declustering conditions for electrospray ionization. A recently developed model peptide is used to arrive at such ion source declustering conditions that prevent the occurrence of intramolecular gas-phase hydrogen (1H/2H) migration (i.e., hydrogen scrambling). This article demonstrates that ETD can be implemented in a mass spectrometric method to monitor the conformational dynamics of proteins in solution at single-residue resolution.

Medin is a principal component of localized amyloid found in the vasculature of individuals over 50 years old. Its amyloid aggregation has been linked to endothelial dysfunction and vascular inflammation, contributing to the pathogenesis of various vascular diseases. Despite its significance, the structures of the medin monomer, oligomer, and fibril remain elusive, and the dynamic processes of medin aggregation are not fully understood. In this study, we comprehensively investigated the medin folding and dimerization dynamics and conformations using atomistic discrete molecular dynamics simulations. Our simulation results suggested that the folding initiation of the medin involved the formation of β-sheets around medin30-41 and medin42-50, with subsequent capping of other segments to their β-sheet edges. Medin monomers typically consisted of three or four β-strands, along with a dynamic N-terminal helix. Two isolated medin peptides readily aggregated into a β-sheet-rich dimer, displaying a strong aggregation propensity. Dimerization of medin not only enhanced the β-sheet conformations but also led to the formation of β-barrel oligomers. The aggregation tendencies of medin1-18 and medin19-29 were relatively weak. However, the segments of medin30-41 and medin42-50 played a crucial role as they primarily formed a β-sheet core and facilitated medin1-18 and medin19-29 to form intra- and interpeptide β-sheets. The findings highlight the critical role of the medin30-41 and medin42-50 regions in stabilizing the monomer structure and driving the medin amyloid aggregation. These regions could potentially serve as promising targets for designing antiamyloid inhibitors against amyloid aggregation of medin. Additionally, our study provides a full picture of the monomer conformations and dimerization dynamics for medin, which will help better understand the pathology of medin aggregation.

Lipid aggregation into fluid bilayers is an essential process for sustaining life. Simplified models of lipid structure, which allow for long time scales or large length scales not obtainable with all-atom simulations, have recently been developed and show promise for describing lipid dynamics in biological systems. Here, we describe two simplified models, a reduced-lipid model and a bola-lipid model for thermophilic bacterial membranes, developed for use with the rapid discrete molecular dynamics simulation method. In the reduced-lipid model, we represent the lipid chain by a series of three beads interacting through pairwise discrete potentials that model hydrophobic attractions between hydrocarbon tails in implicit solvent. Our phase diagram recapitulates those produced by continuous potential models with similar coarse-grained lipid representations. We also find that phase transition temperatures for our reduced-lipid model are dependent upon the flexibility of the lipid chain, giving an insight into archaebacterial membrane stability and prompting development of a bola-lipid model specific for archaebacteria lipids. With both the reduced-lipid and bola-lipid model, we find that the reduced flexibility inherent in archaebacteria lipids yields more stable bilayers as manifested by increased phase transition temperatures. The results of these studies provide a simulation methodology for lipid molecules in biological systems and show that discrete molecular dynamics is applicable to lipid aggregation and dynamics.

The N-methyl-D-aspartate (NMDA)-sensitive glutamate receptor (NMDAR) helps assemble downstream signaling pathways through protein interactions within the postsynaptic density (PSD), which are mediated by its intracellular C-terminal domain (CTD). The most abundant NMDAR subunits in the brain are GluN2A and GluN2B, which are associated with a developmental switch in NMDAR composition. Previously, we used single molecule fluorescence resonance energy transfer (smFRET) to show that the GluN2B CTD contained an intrinsically disordered region with slow, hop-like conformational dynamics. The CTD from GluN2B also undergoes liquid–liquid phase separation (LLPS) with synaptic proteins. Here, we extend these observations to the GluN2A CTD. Sequence analysis showed that both subunits contain a form of intrinsic disorder classified as weak polyampholytes. However, only GluN2B contained matched patterning of arginine and aromatic residues, which are linked to LLPS. To examine the conformational distribution, we used discrete molecular dynamics (DMD), which revealed that GluN2A favors extended disordered states containing secondary structures while GluN2B favors disordered globular states. In contrast to GluN2B, smFRET measurements found that GluN2A lacked slow conformational dynamics. Thus, simulation and experiments found differences in the form of disorder. To understand how this affects protein interactions, we compared the ability of these two NMDAR isoforms to undergo LLPS. We found that GluN2B readily formed condensates with PSD-95 and SynGAP, while GluN2A failed to support LLPS and instead showed a propensity for colloidal aggregation. That GluN2A fails to support this same condensate formation suggests a developmental switch in LLPS propensity.