Transformation between α-helix and β-sheet structures of one and two polyglutamine peptides in explicit water molecules by replica-exchange molecular dynamics simulations

Dissociation process of polyalanine aggregates by free electron laser irradiation

Polyalanine (polyA) disease-causative proteins with an expansion of alanine repeats can be aggregated. Although curative treatments for polyA diseases have not been explored, the dissociation of polyA aggregates likely reduces the cytotoxicity of polyA. Mid-infrared free electron laser (FEL) successfully dissociated multiple aggregates. However, whether the FEL dissociates polyA aggregates like other aggregates has not been tested. Here, we show that FEL at 6.1 μm experimentally weakened the extent of aggregation of a peptide with 13 alanine repeats (13A), and the irradiated 13A exerted lesser cytotoxicity to neuron-like cells than non-irradiated 13A. Then, we applied molecular dynamics (MD) simulation to follow the dissociation process by FEL. We successfully observed how the intermolecular β-sheet of polyA aggregates was dissociated and separated into monomers with helix structures upon FEL irradiation. After the dissociation by FEL, water molecules inhibited the reformation of polyA aggregates. We recently verified the same dissociation process using FEL-treated amyloid-β aggregates. Thus, a common mechanism underlies the dissociation of different protein aggregates that cause different diseases, polyA disease and Alzheimer's disease. However, MD simulation indicated that polyA aggregates are less easily dissociated than amyloid-β aggregates and require longer laser irradiation due to hydrophobic alanine repeats.

The Double-Layered Structure of Amyloid-β Assemblage on GM1-Containing Membranes Catalytically Promotes Fibrillization

Alzheimer's disease (AD) is associated with progressive accumulation of amyloid-β (Aβ) cross-β fibrils in the brain. Aβ species tightly associated with GM1 ganglioside, a glycosphingolipid abundant in neuronal membranes, promote amyloid fibril formation; therefore, they could be attractive clinical targets. However, the active conformational state of Aβ in GM1-containing lipid membranes is still unknown. The present solid-state nuclear magnetic resonance study revealed a nonfibrillar Aβ assemblage characterized by a double-layered antiparallel β-structure specifically formed on GM1 ganglioside clusters. Our data show that this unique assemblage was not transformed into fibrils on GM1-containing membranes but could promote conversion of monomeric Aβ into fibrils, suggesting that a solvent-exposed hydrophobic layer provides a catalytic surface evoking Aβ fibril formation. Our findings offer structural clues for designing drugs targeting catalytically active Aβ conformational species for the development of anti-AD therapeutics.

Ingenuity in performing replica permutation: How to order the state labels for improving sampling efficiency

In the replica-permutation method, an advanced version of the replica-exchange method, all combinations of replicas and parameters are considered for parameter permutation, and a list of all the combinations is prepared. Here, we report that the temperature transition probability depends on how the list is created, especially in replica permutation with solute tempering (RPST). We found that the transition probabilities decrease at large replica indices when the combinations are sequentially assigned to the state labels as in the originally proposed list. To solve this problem, we propose to modify the list by randomly assigning the combinations to the state labels. We performed molecular dynamics simulations of amyloid-β(16-22) peptides using RPST with the "randomly assigned" list (RPST-RA) and RPST with the "sequentially assigned" list (RPST-SA). The results show the decreases in the transition probabilities in RPST-SA are eliminated, and the sampling efficiency is improved in RPST-RA.

Key Residue for Aggregation of Amyloid-β Peptides

It is known that oligomers of amyloid-β (Aβ) peptide are associated with Alzheimer's disease. Aβ has two isoforms: Aβ40 and Aβ42. Although the difference between Aβ40 and Aβ42 is only two additional C-terminal residues, Aβ42 aggregates much faster than Aβ40. It is unknown what role the C-terminal two residues play in accelerating aggregation. Since Aβ42 is more toxic than Aβ40, its oligomerization process needs to be clarified. Moreover, clarifying the differences between the oligomerization processes of Aβ40 and Aβ42 is essential to elucidate the key factors of oligomerization. Therefore, to investigate the dimerization process, which is the early oligomerization process, Hamiltonian replica-permutation molecular dynamics simulations were performed for Aβ40 and Aβ42. We identified a key residue, Arg5, for the Aβ42 dimerization. The two additional residues in Aβ42 allow the C-terminus to form contact with Arg5 because of the electrostatic attraction between them, and this contact stabilizes the β-hairpin. This β-hairpin promotes dimer formation through the intermolecular β-bridges. Thus, we examined the effects of amino acid substitutions of Arg5, thereby confirming that the mutations remarkably suppressed the aggregation of Aβ42. Moreover, the mutations of Arg5 suppressed the Aβ40 aggregation. It was found by analyzing the simulations that Arg5 is important for Aβ40 to form intermolecular contacts. Thus, it was clarified that the role of Arg5 in the oligomerization process varies due to the two additional C-terminal residues.

State-of-the-Art Molecular Dynamics Simulation Studies of RNA-Dependent RNA Polymerase of SARS-CoV-2

Molecular dynamics (MD) simulations are powerful theoretical methods that can reveal biomolecular properties, such as structure, fluctuations, and ligand binding, at the level of atomic detail. In this review article, recent MD simulation studies on these biomolecular properties of the RNA-dependent RNA polymerase (RdRp), which is a multidomain protein, of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are presented. Although the tertiary structures of RdRps in SARS-CoV-2 and SARS-CoV are almost identical, the RNA synthesis activity of RdRp of SARS-CoV is higher than SARS-CoV-2. Recent MD simulations observed a difference in the dynamic properties of the two RdRps, which may cause activity differences. RdRp is also a drug target for Coronavirus disease 2019 (COVID-19). Nucleotide analogs, such as remdesivir and favipiravir, are considered to be taken up by RdRp and inhibit RNA replication. Recent MD simulations revealed the recognition mechanism of RdRp for these drug molecules and adenosine triphosphate (ATP). The ligand-recognition ability of RdRp decreases in the order of remdesivir, favipiravir, and ATP. As a typical recognition process, it was found that several lysine residues of RdRp transfer these ligand molecules to the binding site such as a “bucket brigade.” This finding will contribute to understanding the mechanism of the efficient ligand recognition by RdRp. In addition, various simulation studies on the complexes of SARS-CoV-2 RdRp with several nucleotide analogs are reviewed, and the molecular mechanisms by which these compounds inhibit the function of RdRp are discussed. The simulation studies presented in this review will provide useful insights into how nucleotide analogs are recognized by RdRp and inhibit the RNA replication.

Molecular Dynamics Simulation Studies on the Aggregation of Amyloid-β Peptides and Their Disaggregation by Ultrasonic Wave and Infrared Laser Irradiation

Alzheimer’s disease is understood to be caused by amyloid fibrils and oligomers formed by aggregated amyloid-β (Aβ) peptides. This review article presents molecular dynamics (MD) simulation studies of Aβ peptides and Aβ fragments on their aggregation, aggregation inhibition, amyloid fibril conformations in equilibrium, and disruption of the amyloid fibril by ultrasonic wave and infrared laser irradiation. In the aggregation of Aβ, a β-hairpin structure promotes the formation of intermolecular β-sheet structures. Aβ peptides tend to exist at hydrophilic/hydrophobic interfaces and form more β-hairpin structures than in bulk water. These facts are the reasons why the aggregation is accelerated at the interface. We also explain how polyphenols, which are attracting attention as aggregation inhibitors of Aβ peptides, interact with Aβ. An MD simulation study of the Aβ amyloid fibrils in equilibrium is also presented: the Aβ amyloid fibril has a different structure at one end from that at the other end. The amyloid fibrils can be destroyed by ultrasonic wave and infrared laser irradiation. The molecular mechanisms of these amyloid fibril disruptions are also explained, particularly focusing on the function of water molecules. Finally, we discuss the prospects for developing treatments for Alzheimer’s disease using MD simulations.

Tardigrade Secretory-Abundant Heat-Soluble Protein Varies Entrance Propensity Depending on the Amino-Acid Sequence

Secretory-abundant heat-soluble (SAHS) proteins, which constitute a protein family unique to tardigrades, are thought to be essential for anhydrobiosis. Our previous study has revealed that one of the SAHS proteins of Ramazzottius varieornatus (RvSAHS1) has a more flexible entrance than a mammalian fatty-acid-binding protein, which has a crystal structure similar to that of RvSAHS1. Recently, SAHS paralogs that are expressed abundantly and specifically in the early embryos of this tardigrade and Hypsibius exemplaris have been identified. Comparing these amino-acid sequences with that of RvSAHS1, we have found characteristic differences as I113F and D146T. In this study, we investigate I113F and D146T mutants' properties of RvSAHS1 using molecular dynamics simulations and compare the structures and fluctuations of their entrances with those of the wild type. The two mutants exhibit different properties at the entrance of the β-barrel structure. The I113F mutant tends to close the entrance more than the wild type due to the enhanced hydrophobic network inside the cavity. The D146T mutant, in contrast to the I113F mutant, tends to open the entrance. The mechanism by which this mutation opens the entrance is also discussed. Even though only a single mutation located far from the entrance is added to the wild type, there is a clear difference in the tendency to open and close the β-barrel entrance. It indicates that the entrance properties of the SAHS protein are sensitive to the amino-acid sequence.

Replica permutation with solute tempering for molecular dynamics simulation and its application to the dimerization of amyloid-β fragments

We propose the replica permutation with solute tempering (RPST) by combining the replica-permutation method (RPM) and the replica exchange with solute tempering (REST). Temperature permutations are performed among more than two replicas in RPM, whereas temperature exchanges are performed between two replicas in the replica-exchange method (REM). The temperature transition in RPM occurs more efficiently than in REM. In REST, only the temperatures of the solute region, the solute temperatures, are exchanged to reduce the number of replicas compared to REM. Therefore, RPST is expected to be an improved method taking advantage of these methods. For comparison, we applied RPST, REST, RPM, and REM to two amyloid-β(16–22) peptides in explicit water. We calculated the transition ratio and the number of tunneling events in the temperature space and the number of dimerization events of amyloid-β(16–22) peptides. The results indicate that, in RPST, the number of replicas necessary for frequent random walks in the temperature and conformational spaces is reduced compared to the other three methods. In addition, we focused on the dimerization process of amyloid-β(16–22) peptides. The RPST simulation with a relatively small number of replicas shows that the two amyloid-β(16–22) peptides form the intermolecular antiparallel β-bridges due to the hydrophilic side-chain contact between Lys and Glu and hydrophobic side-chain contact between Leu, Val, and Phe, which stabilizes the dimer of the peptides.

All-Atom Molecular Dynamics Simulation Methods for the Aggregation of Protein and Peptides: Replica Exchange/Permutation and Nonequilibrium Simulations

Protein aggregates are associated with more than 40 serious human diseases. To understand the formation mechanism of protein aggregates at atomic level, all-atom molecular dynamics (MD) simulation is a powerful computational tool. In this chapter, we review the all-atom MD simulation methods that are useful for study on the protein aggregation. We first explain conventional MD simulation methods in physical statistical ensembles, such as the canonical and isothermal-isobaric ensembles. We then describe the generalized-ensemble algorithms such as replica-exchange and replica-permutation MD methods. These methods can overcome a difficulty, in which simulations tend to get trapped in local-minimum free-energy states. Finally we explain the nonequilibrium MD method. Some simulation results based on these methods are also presented.

"Bucket brigade" using lysine residues in RNA-dependent RNA polymerase of SARS-CoV-2

The RNA-dependent RNA polymerase (RdRp) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a promising drug target for coronavirus disease 2019 (COVID-19) because it plays the most important role in the replication of the RNA genome. Nucleotide analogs such as remdesivir and favipiravir are thought to interfere with the RNA replication by RdRp. More specifically, they are expected to compete with nucleoside triphosphates, such as ATP. However, the process in which these drug molecules and nucleoside triphosphates are taken up by RdRp remains unknown. In this study, we performed all-atom molecular dynamics simulations to clarify the recognition mechanism of RdRp for these drug molecules and ATP that were at a distance. The ligand recognition ability of RdRp decreased in the order of remdesivir, favipiravir, and ATP. We also identified six recognition paths. Three of them were commonly found in all ligands, and the remaining three paths were ligand-dependent ones. In the common two paths, it was observed that the multiple lysine residues of RdRp carried the ligands to the binding site like a "bucket brigade." In the remaining common path, the ligands directly reached the binding site. Our findings contribute to the understanding of the efficient ligand recognition by RdRp at the atomic level.

Dynamic properties of SARS-CoV and SARS-CoV-2 RNA-dependent RNA polymerases studied by molecular dynamics simulations

One of the promising drug targets against COVID-19 is an RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2. The tertiary structures of the SARS-CoV-2 and SARS-CoV RdRps are almost the same. However, the RNA-synthesizing activity of the SARS-CoV RdRp is higher than that of the SARS-CoV-2 RdRp. We performed molecular dynamics simulations and found differences in their dynamic properties. In the SARS-CoV RdRp, motifs A-G, which form the active site, are up to 63% closer to each other. We also observed cooperative domain motion in the SARS-CoV RdRp. Such dynamic differences may cause the activity differences between the two RdRps.

Tardigrade Secretory-Abundant Heat-Soluble Protein Has a Flexible β-Barrel Structure in Solution and Keeps This Structure in Dehydration

Secretory-abundant heat-soluble (SAHS) proteins are unique heat-soluble proteins of Tardigrada and are believed to play an essential role in anhydrobiosis, a latent state of life induced by desiccation. To investigate the dynamic properties, molecular dynamics (MD) simulations of a SAHS protein, RvSAHS1, were performed in solution and under dehydrating conditions. For comparison purposes, MD simulations of a human liver-type fatty-acid binding protein (LFABP) were performed in solution. Furthermore, high-speed atomic force microscopy observations were conducted to ascertain the results of the MD simulations. Three properties of RvSAHS1 were found as follows. (1) The entrance region of RvSAHS1 is more flexible and can be more extensive in solutions compared with that of a human LFABP because there is no salt bridge between the βD and βE strands. (2) The intrinsically disordered domain in the N-terminal region significantly fluctuates and can form an amphiphilic α-helix. (3) The size of the entrance region gets smaller along with dehydration, keeping the β-barrel structure. Overall, the obtained results provide atomic-level dynamics of SAHS proteins.

Molecular dynamics simulations of amyloid-β(16-22) peptide aggregation at air-water interfaces

Oligomers of amyloid-β (Aβ) peptides are known to be related to Alzheimer's disease, and their formation is accelerated at hydrophilic-hydrophobic interfaces, such as the cell membrane surface and air-water interface. Here, we report molecular dynamics simulations of aggregation of Aβ(16-22) peptides at air-water interfaces. First, 100 randomly distributed Aβ(16-22) peptides moved to the interface. The high concentration of peptides then accelerated their aggregation and formation of antiparallel β-sheets. Two layers of oligomers were observed near the interface. In the first layer from the interface, the oligomer with less β-bridges exposed the hydrophobic residues to the air. The second layer consisted of oligomers with more β-bridges that protruded into water. They are more soluble in water because the hydrophobic residues are covered by N- and C-terminal hydrophilic residues that are aligned well along the oligomer edge. These results indicate that amyloid protofibril formation mainly occurs in the second layer.

Role of Water Molecules and Helix Structure Stabilization in the Laser-Induced Disruption of Amyloid Fibrils Observed by Nonequilibrium Molecular Dynamics Simulations

Water plays a crucial role in the formation and destruction of biomolecular structures. The mechanism for destroying biomolecular structures was thought to be an active breaking of hydrogen bonds by water molecules. However, using nonequilibrium molecular dynamics simulations, in which an amyloid-β amyloid fibril was destroyed via infrared free-electron laser (IR-FEL) irradiation, we discovered a new mechanism, in which water molecules disrupt protein aggregates. The intermolecular hydrogen bonds formed by C═O and N-H in the fibril are broken at each pulse of laser irradiation. These bonds spontaneously re-form after the irradiation in many cases. However, when a water molecule happens to enter the gap between C═O and N-H, it inhibits the re-formation of the hydrogen bonds. Such sites become defects in the regularly aligned hydrogen bonds, from which all hydrogen bonds in the intermolecular β-sheet are broken as the fraying spreads. This role of water molecules is entirely different from other known mechanisms. This new mechanism can explain the recent experiments showing that the amyloid fibrils are not destroyed by laser irradiation under dry conditions. Additionally, we found that helix structures form more after the amyloid disruption; this is because the resonance frequency is different in a helix structure. Our findings provide a theoretical basis for the application of IR-FEL to the future treatment of amyloidosis.

Dimerization of α-Synuclein Fragments Studied by Isothermal-Isobaric Replica-Permutation Molecular Dynamics Simulation

Aggregates and fibrils of intrinsically disordered α-synuclein are associated with Parkinson's disease. Within a non-amyloid β component (NAC) spanning from the 61st to the 95th residue of α-synuclein, an 11-residue segment called NACore (₆₈GAVVTGVTAVA₇₈) is an essential region for both fibril formation and cytotoxicity. Although NACore peptides alone are known to form aggregates and amyloid fibrils, the mechanisms of aggregation and fibrillation remain unknown. This study investigated the dimerization process of NACore peptides as the initial stage of the aggregation and fibrillation processes. We performed an isothermal-isobaric replica-permutation molecular dynamics simulation, which is one of the efficient sampling methods, for the two NACore peptides in explicit water over 96 μs. The simulation succeeded in sampling a variety of dimer structures. An analysis of secondary structure revealed that most of the NACore dimers form intermolecular β-bridges. In particular, more antiparallel β-bridges were observed than parallel β-bridges. We also found that intramolecular secondary structures such as α-helix and antiparallel β-bridge are stabilized in the pre-dimer state. However, we identified that the intermolecular β-bridges tend to form directly between residues with no specific structure rather than via the intramolecular β-bridges. This is because the NACore peptides still have a low propensity to form the intramolecular secondary structures even though they are stabilized in the pre-dimer state.

Promotion and Inhibition of Amyloid-β Peptide Aggregation: Molecular Dynamics Studies

Aggregates of amyloid-β (Aβ) peptides are known to be related to Alzheimer’s disease. Their aggregation is enhanced at hydrophilic–hydrophobic interfaces, such as a cell membrane surface and air-water interface, and is inhibited by polyphenols, such as myricetin and rosmarinic acid. We review molecular dynamics (MD) simulation approaches of a full-length Aβ peptide, Aβ40, and Aβ(16–22) fragments in these environments. Since these peptides have both hydrophilic and hydrophobic amino acid residues, they tend to exist at the interfaces. The high concentration of the peptides accelerates the aggregation there. In addition, Aβ40 forms a β-hairpin structure, and this structure accelerates the aggregation. We also describe the inhibition mechanism of the Aβ(16–22) aggregation by polyphenols. The aggregation of Aβ(16–22) fragments is caused mainly by the electrostatic attraction between charged amino acid residues known as Lys16 and Glu22. Since polyphenols form hydrogen bonds between their hydroxy and carboxyl groups and these charged amino acid residues, they inhibit the aggregation.

A multifunnel energy landscape encodes the competing α-helix and β-hairpin conformations for a designed peptide

Depending on the amino acid sequence, as well as the local environment, some peptides have the capability to fold into multiple secondary structures. Conformational switching between such structures is a key element of protein folding and aggregation. Specifically, understanding the molecular mechanism underlying the transition from an α-helix to a β-hairpin is critical because it is thought to be a harbinger of amyloid assembly. In this study, we explore the energy landscape for an 18-residue peptide (DP5), designed by Araki and Tamura to exhibit equal propensities for the α-helical and β-hairpin forms. We find that the degeneracy is encoded in the multifunnel nature of the underlying free energy landscape. In agreement with experiment, we also observe that mutation of tyrosine at position 12 to serine shifts the equilibrium in favor of the α-helix conformation, by altering the landscape topography. The transition from the α-helix to the β-hairpin is a complex stepwise process, and occurs via collapsed coil-like intermediates. Our findings suggest that even a single mutation can tune the emergent features of the landscape, providing an efficient route to protein design. Interestingly, the transition pathways for the conformational switch seem to be minimally perturbed upon mutation, suggesting that there could be universal microscopic features that are conserved among different switch-competent protein sequences.

G37V mutation of Aβ42 induces a nontoxic ellipse-like aggregate: An in vitro and in silico study

The glycine zipper motif at the C-terminus of the β-amyloid (Aβ) peptide have been shown to strongly influence the formation of neurotoxic aggregates. A previous study showed that the G37L mutation dramatically reduces the Aβ toxicity in vivo and in vitro. However, the primary cause and mechanism of the glycine zipper motif on Aβ properties remain unknown. To gain molecular insights into the impact of glycine zipper on Aβ properties, we substituted the residue 37 of Glycine by Valine and studied the structural and biochemical properties of G37V mutation, Aβ42(37V), by using in vitro and in silico approaches. Unlike G37L mutation, the G37V mutation reduced toxicity substantially but did not significantly accelerate the aggregation rate or change the content of secondary structures. Further TEM analyses showed that the G37V mutation formed an ellipse-like aggregate rather than a network-like fibril as wild type or G37L mutation of Aβ42 form. This different aggregation morphology may be highly linked with the reduction of toxicity. To gain the insight for the different properties of Aβ42(37V), we studied the structure of Aβ42 and G37V mutation using the replica exchange molecular dynamics simulation. Our results demonstrate that although the overall secondary structure population is similar with Aβ42 and Aβ42(G37V), Aβ42(G37V) shows an increase in the β-turn and β-hairpin at residues 36-37 and the flexibility of the Asp23-Lys28 salt bridge. These unique structural features may be the possible reason to account for the ellipse-like morphology.

Structural Analysis of a Trimer of β2-Microgloblin Fragment by Molecular Dynamics Simulations

A peptide β₂-m_21-31, which is a fragment from residue 21 to residue 31 of β₂-microgloblin, is experimentally known to self-assemble and form amyloid fibrils. In order to understand the mechanism of amyloid fibril formations, we applied the replica-exchange molecular dynamics method to the system consisting of three fragments of β₂-m_21-31. From the analyses on the temperature dependence, we found that there is a clear phase transition temperature in which the peptides aggregate with each other. Moreover, we found by the free energy analyses that there are two major stable states: One of them is like amyloid fibrils and the other is amorphous aggregates.

Structural insights into the aggregation mechanism of huntingtin exon 1 protein fragment with different polyQ-lengths

Huntington disease is a neurodegenerative disorder caused by the expansion of polyglutamine (polyQ) at the N-terminal of the huntingtin exon 1 protein. The detailed structure and the mechanism behind this aggregation remain unclear and it is assumed that the polyQ undergoes a conformational transition to the β-sheet structure when it aggregates. Investigating the misfolding of polyQ facilitates the determination of the molecular mechanism of aggregation and can potentially help in developing a novel approach to inhibit polyQ aggregation. Moreover, the flanking sequences of the polyQ region play a vital role in structural changes and the aggregation mechanism. We performed all-atom molecular dynamics simulations to gain structural insights into the aggregation mechanism using eight different models with glutamine repeat lengths Q₂₇ , Q₂₇ P₁₁ , Q₃₄ , Q₃₅ , Q₃₆ , Q₄₀ , Q₅₀ , and Q₅₀ P₁₁ . In the models without flanking polyPs, we noticed that the transformation of a random coil to β-sheet occurs when the number of Q increases. We also found that the flanking polyPs prevent aggregation by decreasing the probability of forming a β-sheet structure. When polyQ length increases, the 17 N-terminal flanking residues are more likely to adopt a β-sheet conformation from α-helix and coil. From our simulations, we suggest that at least 34 glutamines are required for initiating aggregation and 40 residues length is critical for the aggregation of huntingtin exon 1 protein for disease onset. This study provides structural insights into misfolding and the role of flanking sequences in huntingtin aggregation which will further help in developing therapeutic strategies for Huntington's disease.

Effects of a Hydrophilic/Hydrophobic Interface on Amyloid-β Peptides Studied by Molecular Dynamics Simulations and NMR Experiments

Oligomer formation of amyloid-β peptides (Aβ) is accelerated at a hydrophilic/hydrophobic interface. However, details of the acceleration mechanism have not been elucidated. To understand the effects of the interface on oligomerization at the atomic level, we performed molecular dynamics simulations for an Aβ40 monomer in the presence and absence of the hydrophilic/hydrophobic interface. Nuclear magnetic resonance experiments of Aβ40 peptides with gangliosidic micelles were also carried out. In the simulations and experiments, the hydrophobic residues of Aβ40 bound to the interface stably. Moreover, we found that Aβ40 formed a hairpin structure at the interface more readily than in bulk water. From these results, we discussed the acceleration mechanism of the oligomer formation at the interface.

Theoretical approaches for dynamical ordering of biomolecular systems

BACKGROUND

Living systems are characterized by the dynamic assembly and disassembly of biomolecules. The dynamical ordering mechanism of these biomolecules has been investigated both experimentally and theoretically. The main theoretical approaches include quantum mechanical (QM) calculation, all-atom (AA) modeling, and coarse-grained (CG) modeling. The selected approach depends on the size of the target system (which differs among electrons, atoms, molecules, and molecular assemblies). These hierarchal approaches can be combined with molecular dynamics (MD) simulation and/or integral equation theories for liquids, which cover all size hierarchies.

SCOPE OF REVIEW

We review the framework of quantum mechanical/molecular mechanical (QM/MM) calculations, AA MD simulations, CG modeling, and integral equation theories. Applications of these methods to the dynamical ordering of biomolecular systems are also exemplified.

MAJOR CONCLUSIONS

The QM/MM calculation enables the study of chemical reactions. The AA MD simulation, which omits the QM calculation, can follow longer time-scale phenomena. By reducing the number of degrees of freedom and the computational cost, CG modeling can follow much longer time-scale phenomena than AA modeling. Integral equation theories for liquids elucidate the liquid structure, for example, whether the liquid follows a radial distribution function.

GENERAL SIGNIFICANCE

These theoretical approaches can analyze the dynamic behaviors of biomolecular systems. They also provide useful tools for exploring the dynamic ordering systems of biomolecules, such as self-assembly. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

Molecular dynamics analysis of the aggregation propensity of polyglutamine segments

Protein misfolding and aggregation is a pathogenic feature shared among at least ten polyglutamine (polyQ) neurodegenerative diseases. While solvent-solution interaction is a key factor driving protein folding and aggregation, the solvation properties of expanded polyQ tracts are not well understood. By using GPU-enabled all-atom molecular dynamics simulations of polyQ monomers in an explicit solvent environment, this study shows that solvent-polyQ interaction propensity decreases as the lengths of polyQ tract increases. This study finds a predominance in long-distance interactions between residues far apart in polyQ sequences with longer polyQ segments, that leads to significant conformational differences. This study also indicates that large loops, comprised of parallel β-structures, appear in long polyQ tracts and present new aggregation building blocks with aggregation driven by long-distance intra-polyQ interactions. Finally, consistent with previous observations using coarse-grain simulations, this study demonstrates that there is a gain in the aggregation propensity with increased polyQ length, and that this gain is correlated with decreasing ability of solvent-polyQ interaction. These results suggest the modulation of solvent-polyQ interactions as a possible therapeutic strategy for treating polyQ diseases.

Structural and fluctuational difference between two ends of Aβ amyloid fibril: MD simulations predict only one end has open conformations

Aβ amyloid fibrils, which are related to Alzheimer’s disease, have a cross-β structure consisting of two β-sheets: β1 and β2. The Aβ peptides are thought to be serially arranged in the same molecular conformation along the fibril axis. However, to understand the amyloid extension mechanism, we must understand the amyloid fibril structure and fluctuation at the fibril end, which has not been revealed to date. Here, we reveal these features by all-atom molecular dynamics (MD) simulations of Aβ42 and Aβ40 fibrils in explicit water. The structure and fluctuation were observed to differ between the two ends. At the even end, the Aβ peptide always took a closed form wherein β1 and β2 were closely spaced. The Aβ peptide fluctuated more at the odd end and took an open form wherein the two β-sheets were well separated. The differences are attributed to the stronger β-sheet formation by the β1 exposed at the even end than the β2 exposed at the odd end. Along with the small fluctuations at the even end, these results explain why the fibril extends from one end only, as observed in experiments. Our MD results agree well with recent observations by high-speed atomic force microscopy.

Simulated tempering based on global balance or detailed balance conditions: Suwa-Todo, heat bath, and Metropolis algorithms

Simulated tempering (ST) is a useful method to enhance sampling of molecular simulations. When ST is used, the Metropolis algorithm, which satisfies the detailed balance condition, is usually applied to calculate the transition probability. Recently, an alternative method that satisfies the global balance condition instead of the detailed balance condition has been proposed by Suwa and Todo. In this study, ST method with the Suwa-Todo algorithm is proposed. Molecular dynamics simulations with ST are performed with three algorithms (the Metropolis, heat bath, and Suwa-Todo algorithms) to calculate the transition probability. Among the three algorithms, the Suwa-Todo algorithm yields the highest acceptance ratio and the shortest autocorrelation time. These suggest that sampling by a ST simulation with the Suwa-Todo algorithm is most efficient. In addition, because the acceptance ratio of the Suwa-Todo algorithm is higher than that of the Metropolis algorithm, the number of temperature states can be reduced by 25% for the Suwa-Todo algorithm when compared with the Metropolis algorithm.

Exploring biomolecular dynamics and interactions using advanced sampling methods

Molecular dynamics (MD) and Monte Carlo (MC) simulations have emerged as a valuable tool to investigate statistical mechanics and kinetics of biomolecules and synthetic soft matter materials. However, major limitations for routine applications are due to the accuracy of the molecular mechanics force field and due to the maximum simulation time that can be achieved in current simulations studies. For improving the sampling a number of advanced sampling approaches have been designed in recent years. In particular, variants of the parallel tempering replica-exchange methodology are widely used in many simulation studies. Recent methodological advancements and a discussion of specific aims and advantages are given. This includes improved free energy simulation approaches and conformational search applications.

Effect of pH on the Aggregation of α-syn12 Dimer in Explicit Water by Replica-Exchange Molecular Dynamics Simulation

The dimeric structure of the N-terminal 12 residues drives the interaction of α-synuclein protein with membranes. Moreover, experimental studies indicated that the aggregation of α-synuclein is faster at low pH than neutral pH. Nevertheless, the effects of different pH on the structural characteristics of the α-syn12 dimer remain poorly understood. We performed 500 ns temperature replica exchange molecular dynamics (T-REMD) simulations of two α-syn12 peptides in explicit solvent. The free energy surfaces contain ten highly populated regions at physiological pH, while there are only three highly populated regions contained at acidic pH. The anti-parallel β-sheet conformations were found as the lowest free energy state. Additionally, these states are nearly flat with a very small barrier which indicates that these states can easily transit between themselves. The dimer undergoes a disorder to order transition from physiological pH to acidic pH and the α-syn12 dimer at acidic pH involves a faster dimerization process. Further, the Lys6–Asp2 contact may prevent the dimerization.

Dimerization process of amyloid-β(29-42) studied by the Hamiltonian replica-permutation molecular dynamics simulations

The amyloid-β peptides form amyloid fibrils which are associated with Alzheimer's disease. Amyloid-β(29-42) is its C-terminal fragment and a critical determinant of the amyloid formation rate. This fragment forms the amyloid fibril by itself. However, the fragment conformation in the fibril has yet to be determined. The oligomerization process including the dimerization process is also still unknown. The dimerization process corresponds to an early process of the amyloidogenesis. In order to investigate the dimerization process and conformations, we applied the Hamiltonian replica-permutation method, which is a better alternative to the Hamiltonian replica-exchange method, to two amyloid-β(29-42) molecules in explicit water solvent. At the first step of the dimerization process, two amyloid-β(29-42) molecules came close to each other and had intermolecular side chain contacts. When two molecules had the intermolecular side chain contacts, the amyloid-β(29-42) tended to have intramolecular secondary structures, especially β-hairpin structures. The two molecules had intermolecular β-bridge structures by coming much closer at the second step of the dimerization process. Formation of these intermolecular β-bridge structures was induced by the β-hairpin structures. The intermolecular β-sheet structures elongated at the final step. Structures of the amyloid-β(29-42) in the monomer and dimer states are also shown with the free-energy landscapes, which were obtained by performing efficient sampling in the conformational space in our simulations.