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

Inhibition of α-Synuclein Misfolding into β-Sheet Domains on Medium-Sized Gold Nanoclusters: Evidence from Enhanced Sampling MD Simulations

Targeting Parkinson's disease (PD) related protein, α-synuclein (αS), via gold nanoclusters (AuNCs) has received considerable attention in PD treatments, but its molecular basis on the initial interactions between αS and AuNCs remains elusive due to the absence of a unique secondary structure of αS chains. Here, at the single-cluster level, we incorporate well-tempered metadynamics simulations to explore the structural and thermodynamic characteristics of the full length αS adsorbed on different-sized AuNCs (Aun, n = 25, 36, 44, 68, 102) with modeled thiolated ligands (Aun@Lig). The conformational landscapes of αS indicate that uncharged Aun@SCH2OH chaperones the native intrinsically disordered conformations of αS, while negatively and positively charged AuNCs greatly increase the likelihood of forming intramolecular β-sheet domains, which are necessary for αS fibrillation and are a hallmark of PD. The binding details further demonstrate the significant inhibitory effect of the medium-sized Au36@SCH2OH on αS misfolding into β-sheet domains. This provides a valuable guideline for customizing AuNCs to precisely manipulate protein folding and misfolding behaviors, with potential implications for disease treatments.

Why Is Arginine the Only Amino Acid That Inhibits Polyglutamine Monomers from Taking on Toxic Conformations?

Polyglutamine (polyQ) diseases are devastating neurodegenerative disorders characterized by abnormal expansion of glutamine repeats within specific proteins. The aggregation of polyQ proteins is a critical pathological hallmark of these diseases. Arginine was identified as a promising inhibitory compound because it prevents polyQ-protein monomers from forming intra- and intermolecular β-sheet structures and hinders polyQ proteins from aggregating to form oligomers. Such an aggregation inhibitory effect was not observed in other amino acids. However, the underlying molecular mechanism of the aggregation inhibition and the factors that differentiate arginine from other amino acids, in terms of the inhibition of the polyQ-protein aggregation, remain poorly understood. Here, we performed replica-permutation molecular dynamics simulations to elucidate the molecular mechanism by which arginine inhibits the formation of the intramolecular β-sheet structure of a polyQ monomer. We found that the intramolecular β-sheet structure with more than four β-bridges of the polyQ monomer with arginine is more unstable than without any ligand and with lysine. We also found that arginine has 1.6-2.1 times more contact with polyQ than lysine. In addition, we revealed that arginine forms more hydrogen bonds with the main chain of the polyQ monomer than lysine. More hydrogen bonds formed between arginine and polyQ inhibit polyQ from forming the long intramolecular β-sheet structure. It is known that intramolecular β-sheet structure enhances intermolecular β-sheet structure between proteins. These effects are thought to be the reason for the inhibition of polyQ aggregation. This study provides insights into the molecular events underlying arginine's inhibition of polyQ-protein aggregation.

Antimicrobial drugs for Parkinson's disease: Existing therapeutic strategies and novel drugs exploration

Parkinson's disease (PD), the second most common neurodegenerative disorder, is characterized by loss of dopaminergic neurons in the substantia nigra, as well as the abnormal accumulation of misfolded α-synuclein. Clinically, PD is featured by typical motor symptoms and some non-motor symptoms. Up to now, although considerable progress has been made in understanding the pathogenesis of PD, there is still no effective therapeutic treatment for the disease. Thus, exploring new therapeutic strategies has been a topic that needs to be addressed urgently. Noteworthy, with the proposal of the microbiota-gut-brain axis theory, antimicrobial drugs have received significant attention due to their effects on regulating the intestinal microbiota. Nowadays, there is growing evidence showing that some antimicrobial drugs may be promising drugs for the treatment of PD. Data from pre-clinical and clinical studies have shown that some antimicrobial drugs may play neuroprotective roles in PD by modulating multiple biochemical and molecular pathways, including reducing α-synuclein aggregation, inhibiting neuroinflammation, regulating mitochondrial structure and function, as well as suppressing oxidative stress. In this paper, we summarized the effects of some antimicrobial drugs on PD treatment from recent pre-clinical and clinical studies. Then, we further discussed the potential of a few antimicrobial drugs for treating PD based on molecular docking and molecular dynamics simulation. Importantly, we highlighted the potential of clorobiocin as the therapeutic strategy for PD owing to its ability to inhibit α-synuclein aggregation. These results will help us to better understand the potential of antimicrobial drugs in treating PD and how antimicrobial drugs may alleviate or reverse the pathological symptoms of PD.

Hysteresis Elimination for an Anisotropic Liquid-Crystal Model via Molecule Design and Replica-Exchange Optimization

The phenomenon of hysteresis in simulations, in which a system's current state is correlated to previous states and inhibits the transition to a more stable phase, may often lead to misleading results in physical chemistry. In this study, in addition to the replica exchange method (REM), a novel approach was taken by combining an evolution strategy based on the evolutionary principles of nature to predict phase transitions for the Hess-Su liquid-crystal model. In this model, an anisotropy term is added to the simple 6-12 Lennard-Jones model to intuitively reproduce the behavior of liquid crystals. We first applied the pressure-temperature REM to the Hess-Su model and optimized the replica spacing for the energy distribution to gain the maximum advantage from the REM. We then used the same approach as for the Hamiltonian REM, seeking to optimize the replica spacing in the same way. Based on both results, we attempted to predict this coarse-grained liquid-crystal model's exact phase transition point. In the Hamiltonian REM, replicas were prepared with different molecular aspect ratios corresponding to the values of the anisotropy terms in the potential function. The Hess-Su liquid-crystal model, which undergoes a direct transition from the nematic to the solid phase without going through a smectic phase, is a challenging research target for understanding phase transitions. Despite the tremendous computational difficulty in overcoming the strong hysteresis present in this system, our method could predict the phase transition point clearly and significantly reduce the extent of hysteresis. Our approach is beneficial when simulating more complex systems and, above all, shows great potential for more accurate and efficient phase transition predictions in the field of molecular simulation in the future.

Interactions between Lipid Vesicle Membranes and Single Amino Acid Fibrils: Probable Origin of Specific Neurological Disorders

Amyloid fibrils are known to be responsible for several neurological disorders, like Alzheimer's disease (AD), Parkinson's disease (PD), etc. For decades, mostly proteins and peptide-based amyloid fibrils have been focused on, and the topic has acknowledged the rise, development, understanding of, and controversy, as well. However, the single amino acid based amyloid fibrils, responsible for several disorders, such as phenylketonuria, tyrosenimia type II, hypermethioninemia, etc., have gotten scientific attention lately. To understand the molecular level pathogenesis of such disorders originated due to the accumulation of single amino acid-based amyloid fibrils, interaction of these fibrils with phospholipid vesicle membranes is found to be an excellent cell-free in vitro setup. Based on such an in vitro setup, these fibrils show a generic mechanism of membrane insertion driven by electrostatic and hydrophobic effects inside the membrane that reduces the integral rigidity of the membrane. Alteration of such fundamental properties of the membrane, therefore, might be referred to as one of the prime pathological factors for the development of these neurological disorders. Hence, such interactions must be investigated in cellular and intracellular compartments to design suitable therapeutic modulators against fibrils.

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.

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.

Insights into aggregation dynamics of NACore peptides from coarse-grained simulations

Alpha(α)-synuclein is closely related to the pathogenesis of Parkinson's disease (PD). The NACore, a fragment of α-synuclein, is considered to be the key region of α-synuclein that causes PD. The aggregation dynamics of NACores are studied via coarse-grained molecular dynamics simulations. We find that NACores can self-assemble into a large cluster at high concentrations. The aggregation dynamics can be divided into three stages. The growth kinetics for the first and second stages follows the power law, S_max  ~ t^γ , with the second stage faster than the first one. The characteristic lifetime for the high concentration is 40 times larger than that for the low concentration, implying the low fluidity. Understanding the aggregation dynamics of NACores is helpful to develop drugs for therapeutic prevention and intervention.

Wild-Type α-Synuclein and Variants Occur in Different Disordered Dimers and Pre-Fibrillar Conformations in Early Stage of Aggregation

α-Synuclein is a 140 amino-acid intrinsically disordered protein mainly found in the brain. Toxic α-synuclein aggregates are the molecular hallmarks of Parkinson’s disease. In vitro studies showed that α-synuclein aggregates in oligomeric structures of several 10th of monomers and into cylindrical structures (fibrils), comprising hundred to thousands of proteins, with polymorphic cross-β-sheet conformations. Oligomeric species, formed at the early stage of aggregation remain, however, poorly understood and are hypothezised to be the most toxic aggregates. Here, we studied the formation of wild-type (WT) and mutant (A30P, A53T, and E46K) dimers of α-synuclein using coarse-grained molecular dynamics. We identified two principal segments of the sequence with a higher propensity to aggregate in the early stage of dimerization: residues 36–55 and residues 66–95. The transient α-helices (residues 53–65 and 73–82) of α-synuclein monomers are destabilized by A53T and E46K mutations, which favors the formation of fibril native contacts in the N-terminal region, whereas the helix 53–65 prevents the propagation of fibril native contacts along the sequence for the WT in the early stages of dimerization. The present results indicate that dimers do not adopt the Greek key motif of the monomer fold in fibrils but form a majority of disordered aggregates and a minority (9–15%) of pre-fibrillar dimers both with intra-molecular and intermolecular β-sheets. The percentage of residues in parallel β-sheets is by increasing order monomer < disordered dimers < pre-fibrillar dimers. Native fibril contacts between the two monomers are present in the NAC domain for WT, A30P, and A53T and in the N-domain for A53T and E46K. Structural properties of pre-fibrillar dimers agree with rupture-force atomic force microscopy and single-molecule Förster resonance energy transfer available data. This suggests that the pre-fibrillar dimers might correspond to the smallest type B toxic oligomers. The probability density of the dimer gyration radius is multi-peaks with an average radius that is 10 Å larger than the one of the monomers for all proteins. The present results indicate that even the elementary α-synuclein aggregation step, the dimerization, is a complicated phenomenon that does not only involve the NAC region.

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.

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.

Strong inhibition of peptide amyloid formation by a fatty acid

The aggregation of peptides into amyloid fibrils is associated with several diseases, including Alzheimer’s and Parkinson’s disease. Because hydrophobic interactions often play an important role in amyloid formation, the presence of various hydrophobic or amphiphilic molecules, such as lipids, may influence the aggregation process. We have studied the effect of a fatty acid, linoleic acid, on the fibrillation process of the amyloid-forming model peptide NACore (GAVVTGVTAVA). NACore is a peptide fragment spanning residue 68–78 of the protein α-synuclein involved in Parkinson’s disease. Based primarily on circular dichroism measurements, we found that even a very small amount of linoleic acid can substantially inhibit the fibrillation of NACore. This inhibitory effect manifests itself through a prolongation of the lag phase of the peptide fibrillation. The effect is greatest when the fatty acid is present from the beginning of the process together with the monomeric peptide. Cryogenic transmission electron microscopy revealed the presence of nonfibrillar clusters among NACore fibrils formed in the presence of linoleic acid. We argue that the observed inhibitory effect on fibrillation is due to co-association of peptide oligomers and fatty acid aggregates at the early stage of the process. An important aspect of this mechanism is that it is nonmonomeric peptide structures that associate with the fatty acid aggregates. Similar mechanisms of action could be relevant in amyloid formation occurring in vivo, where the aggregation takes place in a lipid-rich environment.

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.