Effects of force fields on the conformational and dynamic properties of amyloid β(1-40) dimer explored by replica exchange molecular dynamics simulations

Conformation-Dependent Influences of Hydrophobic Amino Acids in Two In-Register Parallel β-Sheet Amyloids, an α-Synuclein Amyloid and a Local Structural Model of PrPSc

Prions are unconventional pathogens that encode the pathogenic information in conformations of the constituent abnormal isoform of prion protein (PrP^Sc), independently of the nucleotide genome. Therefore, conformational diversity of PrP^Sc underlies the existence of many prion strains and species barriers of prions, although the conformational information is extremely limited. Interestingly, differences between polymorphic or species-specific residues responsible for the species/strain barriers are often caused by conservative replacements between hydrophobic amino acids. This implies that subtle differences among hydrophobic amino acids are significant for PrP^Sc structures. Here we analyzed the influence of different hydrophobic residues on the structures of an in-register parallel β-sheet amyloid of α-synuclein (αSyn) using molecular dynamics (MD) simulation and applied the knowledge from the αSyn amyloid to modeling a local structure of human PrP^Sc encompassing residues 107-143. We found that mutations equivalent to polymorphisms that cause transmission barriers substantially affect the stabilities of the local structures; for example, the G127V mutation, which makes the host resistant to various human prion diseases, greatly destabilized the local structure of the model amyloid. Our study indicates that subtle differences among hydrophobic side chains can considerably affect the interaction network, including hydrogen bonds, and demonstrates specifically how and in what structures hydrophobic residues can exert unique effects on in-register parallel β-sheet amyloids.

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.

Forcefield evaluation and accelerated molecular dynamics simulation of Zn(II) binding to N-terminus of amyloid-β

We report conventional and accelerated molecular dynamics simulation of Zn(II) bound to the N-terminus of amyloid-β. By comparison against NMR data for the experimentally determined binding mode, we find that certain combinations of forcefield and solvent model perform acceptably in describing the size, shape and secondary structure, and that there is no appreciable difference between implicit and explicit solvent models. We therefore used the combination of ff14SB forcefield and GBSA solvent model to compare the result of different binding modes of Zn(II) to the same peptide, using accelerated MD to enhance sampling and comparing the free peptide simulated in the same way. We show that Zn(II) imparts significant rigidity to the peptide, disrupts the secondary structure and pattern of salt bridges seen in the free peptide, and induces closer contact between residues. Free energy surfaces in 1 or 2 dimensions further highlight the effect of metal coordination on peptide's spatial extent. We also provide evidence that accelerated MD provides improved sampling over conventional MD by visiting as many or more configurations in much shorter simulation times.

The Pathogenesis Mechanism, Structure Properties, Potential Drugs and Therapeutic Nanoparticles against the Small Oligomers of Amyloid-β

Alzheimer's Disease (AD) is a devastating neurodegenerative disease that affects millions of people in the world. The abnormal aggregation of amyloid β protein (Aβ) is regarded as the key event in AD onset. Meanwhile, the Aβ oligomers are believed to be the most toxic species of Aβ. Recent studies show that the Aβ dimers, which are the smallest form of Aβ oligomers, also have the neurotoxicity in the absence of other oligomers in physiological conditions. In this review, we focus on the pathogenesis, structure and potential therapeutic molecules against small Aβ oligomers, as well as the nanoparticles (NPs) in the treatment of AD. In this review, we firstly focus on the pathogenic mechanism of Aβ oligomers, especially the Aβ dimers. The toxicity of Aβ dimer or oligomers, which attributes to the interactions with various receptors and the disruption of membrane or intracellular environments, were introduced. Then the structure properties of Aβ dimers and oligomers are summarized. Although some structural information such as the secondary structure content is characterized by experimental technologies, detailed structures are still absent. Following that, the small molecules targeting Aβ dimers or oligomers are collected; nevertheless, all of these ligands have failed to come into the market due to the rising controversy of the Aβ-related "amyloid cascade hypothesis". At last, the recent progress about the nanoparticles as the potential drugs or the drug delivery for the Aβ oligomers are present.

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.

Folding and self-assembly of short intrinsically disordered peptides and protein regions

Proteins and peptide fragments are highly relevant building blocks in self-assembly for nanostructures with plenty of applications. Intrinsically disordered proteins (IDPs) and protein regions (IDRs) are defined by the absence of a well-defined secondary structure, yet IDPs/IDRs show a significant biological activity. Experimental techniques and computational modelling procedures for the characterization of IDPs/IDRs are discussed. Directed self-assembly of IDPs/IDRs allows reaching a large variety of nanostructures. Hybrid materials based on the derivatives of IDPs/IDRs show a promising performance as alternative biocides and nanodrugs. Cell mimicking, in vivo compartmentalization, and bone regeneration are demonstrated for IDPs/IDRs in biotechnological applications. The exciting possibilities of IDPs/IDRs in nanotechnology with relevant biological applications are shown.

Activation of G-protein-coupled receptors is thermodynamically linked to lipid solvation

Preferential lipid solvation of the G-protein-coupled A_2A adenosine receptor (A_2AR) is evaluated from 35 μs of all-atom molecular dynamics simulation. A coarse-grained transition matrix algorithm is developed to overcome slow equilibration of the first solvation shell, obtaining estimates of the free energy of solvation by different lipids for the receptor in different activation states. Results indicate preference for solvation by unsaturated chains, which favors the active receptor. A model for lipid-dependent G-protein-coupled receptor activity is proposed in which the chemical potential of lipids in the bulk membrane modulates receptor activity. The entropies associated with moving saturated and unsaturated lipids from bulk to A_2AR's first solvation shell are evaluated. Overall, the acyl chains are more disordered (i.e., obtain a favorable entropic contribution) when partitioning to the receptor surface, and this effect is augmented for the saturated chains, which are relatively more ordered in bulk.

Evaluating the Multifunctionality of a New Modulator of Zinc-Induced Aβ Aggregation Using a Novel Computational Approach

The high concentration of zinc metal ions in Aβ aggregations is one of the most cited hallmarks of Alzheimer's disease (AD), and several substantial pieces of evidence emphasize the key role of zinc metal ions in the pathogenesis of AD. In this study, while designing a multifunctional peptide for simultaneous targeting Aβ aggregation and chelating the zinc metal ion, a novel and comprehensive approach is introduced for evaluating the multifunctionality of a multifunctional drugs based on computational methods. The multifunctional peptide consists of inhibitor and chelator domains, which are included in the C-terminal hydrophobic region of Aβ, and the first four amino acids of human albumin. The ability of the multifunctional peptide in zinc ion chelation has been investigated using molecular dynamics (MD) simulations of the peptide-zinc interaction for 300 ns, and Bennett's acceptance ratio (BAR) method has been used to accurately calculate the chelation free energy. Data analysis demonstrates that the peptide chelating domain can be stably linked to the zinc ion. Besides, the introduced method used for evaluating chelation and calculating the free energy of peptide binding to zinc ions was successfully validated by comparison with previous experimental and theoretical published data. The results indicate that the multifunctional peptide, coordinating with the zinc metal ion, can be effective in Aβ inhibition by preserving the native helical structure of the Aβ₄₂ monomer as well as disrupting the β-sheet structure of Aβ₄₂ aggregates. Detailed assessments of the Aβ₄₂-peptide interactions elucidate that the inhibition of Aβ is achieved by considerable hydrophobic interactions and hydrogen bonding between the multifunctional peptide and the hydrophobic Aβ regions, along with interfering in stable bridges formed inside the Aβ aggregate.

On the Conformational Dynamics of β-Amyloid Forming Peptides: A Computational Perspective

Understanding the conformational dynamics of proteins and peptides involved in important functions is still a difficult task in computational structural biology. Because such conformational transitions in β-amyloid (Aβ) forming peptides play a crucial role in many neurological disorders, researchers from different scientific fields have been trying to address issues related to the folding of Aβ forming peptides together. Many theoretical models have been proposed in the recent years for studying Aβ peptides using mathematical, physicochemical, and molecular dynamics simulation, and machine learning approaches. In this article, we have comprehensively reviewed the developmental advances in the theoretical models for Aβ peptide folding and interactions, particularly in the context of neurological disorders. Furthermore, we have extensively reviewed the advances in molecular dynamics simulation as a tool used for studying the conversions between polymorphic amyloid forms and applications of using machine learning approaches in predicting Aβ peptides and aggregation-prone regions in proteins. We have also provided details on the theoretical advances in the study of Aβ peptides, which would enhance our understanding of these peptides at the molecular level and eventually lead to the development of targeted therapies for certain acute neurological disorders such as Alzheimer's disease in the future.

Combining molecular dynamics simulations and experimental analyses in protein misfolding

The fold of a protein determines its function and its misfolding can result in loss-of-function defects. In addition, for certain proteins their misfolding can lead to gain-of-function toxicities resulting in protein misfolding diseases such as Alzheimer's, Parkinson's, or the prion diseases. In all of these diseases one or more proteins misfold and aggregate into disease-specific assemblies, often in the form of fibrillar amyloid deposits. Most, if not all, protein misfolding diseases share a fundamental molecular mechanism that governs the misfolding and subsequent aggregation. A wide variety of experimental methods have contributed to our knowledge about misfolded protein aggregates, some of which are briefly described in this review. The misfolding mechanism itself is difficult to investigate, as the necessary timescale and resolution of the misfolding events often lie outside of the observable parameter space. Molecular dynamics simulations fill this gap by virtue of their intrinsic, molecular perspective and the step-by-step iterative process that forms the basis of the simulations. This review focuses on molecular dynamics simulations and how they combine with experimental analyses to provide detailed insights into protein misfolding and the ensuing diseases.

Effects of in vivo conditions on amyloid aggregation

One of the grand challenges of biophysical chemistry is to understand the principles that govern protein misfolding and aggregation, which is a highly complex process that is sensitive to initial conditions, operates on a huge range of length- and timescales, and has products that range from protein dimers to macroscopic amyloid fibrils. Aberrant aggregation is associated with more than 25 diseases, which include Alzheimer's, Parkinson's, Huntington's, and type II diabetes. Amyloid aggregation has been extensively studied in the test tube, therefore under conditions that are far from physiological relevance. Hence, there is dire need to extend these investigations to in vivo conditions where amyloid formation is affected by a myriad of biochemical interactions. As a hallmark of neurodegenerative diseases, these interactions need to be understood in detail to develop novel therapeutic interventions, as millions of people globally suffer from neurodegenerative disorders and type II diabetes. The aim of this review is to document the progress in the research on amyloid formation from a physicochemical perspective with a special focus on the physiological factors influencing the aggregation of the amyloid-β peptide, the islet amyloid polypeptide, α-synuclein, and the hungingtin protein.

Transition Metal Ion Interactions with Disordered Amyloid-β Peptides in the Pathogenesis of Alzheimer's Disease: Insights from Computational Chemistry Studies

Monomers and oligomers of the amyloid-β peptide aggregate to form the fibrils found in the brains of Alzheimer's disease patients. These monomers and oligomers are largely disordered and can interact with transition metal ions, affecting the mechanism and kinetics of amyloid-β aggregation. Due to the disordered nature of amyloid-β, its rapid aggregation, as well as solvent and paramagnetic effects, experimental studies face challenges in the characterization of transition metal ions bound to amyloid-β monomers and oligomers. The details of the coordination chemistry between transition metals and amyloid-β obtained from experiments remain debated. Furthermore, the impact of transition metal ion binding on the monomeric or oligomeric amyloid-β structures and dynamics are still poorly understood. Computational chemistry studies can serve as an important complement to experimental studies and can provide additional knowledge on the binding between amyloid-β and transition metal ions. Many research groups conducted first-principles calculations, ab initio molecular dynamics simulations, quantum mechanics/classical mechanics simulations, and classical molecular dynamics simulations for studying the interplay between transition metal ions and amyloid-β monomers and oligomers. This review summarizes the current understanding of transition metal interactions with amyloid-β obtained from computational chemistry studies. We also emphasize the current view of the coordination chemistry between transition metal ions and amyloid-β. This information represents an important foundation for future metal ion chelator and drug design studies aiming to combat Alzheimer's disease.

Norepinephrine Inhibits Alzheimer's Amyloid-β Peptide Aggregation and Destabilizes Amyloid-β Protofibrils: A Molecular Dynamics Simulation Study

The abnormal self-assembly of amyloid-β (Aβ) peptides into toxic fibrillar aggregates is associated with the pathogenesis of Alzheimer's disease (AD). The inhibition of β-sheet-rich oligomer formation is considered as the primary therapeutic strategy for AD. Previous experimental studies reported that norepinephrine (NE), one of the neurotransmitters, is able to inhibit Aβ aggregation and disaggregate the preformed fibrils. Moreover, exercise can markedly increase the level of NE. However, the underlying inhibitory and disruptive mechanisms remain elusive. In this work, we performed extensive replica-exchange molecular dynamic (REMD) simulations to investigate the conformational ensemble of Aβ_1-42 dimer with and without NE molecules. Our results show that without NE molecules, Aβ_1-42 dimer transiently adopts a β-hairpin-containing structure, and the β-strand regions of this β-hairpin (residues 15QKLVFFA21 and 33GLMVGGVV40) strongly resemble those of the Aβ fibril structure (residues 15QKLVFFA21 and 30AIIGLMVG37) reported in an electron paramagnetic resonance spectroscopy study. NE molecules greatly reduce the interpeptide β-sheet content and suppress the formation of the above-mentioned β-hairpin, leading to a more disordered coil-rich Aβ dimer. Five dominant binding sites are identified, and the central hydrophobic core 16KLVFFA21 site and C-terminal 31IIGLMV36 hydrophobic site are the two most favorable ones. Our data reveal that hydrophobic, aromatic stacking, hydrogen-bonding and cation-π interactions synergistically contribute to the binding of NE molecules to Aβ peptides. MD simulations of Aβ_1-42 protofibril show that NE molecules destabilize Aβ protofibril by forming H-bonds with residues D1, A2, D23, and A42. This work reveals the molecular mechanism by which NE molecules inhibit Aβ_1-42 aggregation and disaggregate Aβ protofibrils, providing valuable information for developing new drug candidates and exercise therapy against AD.

Effects of All-Atom Molecular Mechanics Force Fields on Amyloid Peptide Assembly: The Case of Aβ16-22 Dimer

We investigated the effects of 17 widely used atomistic molecular mechanics force fields (MMFFs) on the structures and kinetics of amyloid peptide assembly. To this end, we performed large-scale all-atom molecular dynamics simulations in explicit water on the dimer of the seven-residue fragment of the Alzheimer's amyloid-β peptide, Aβ16-22, for a total time of 0.34 ms. We compared the effects of these MMFFs by analyzing various global reaction coordinates, secondary structure contents, the fibril population, the in-register and out-of-register architectures, and the fibril formation time at 310 K. While the AMBER94, AMBER99, and AMBER12SB force fields do not predict any β-sheets, the seven force fields, AMBER96, GROMOS45a3, GROMOS53a5, GROMOS53a6, GROMOS43a1, GROMOS43a2, and GROMOS54a7, form β-sheets rapidly. In contrast, the following five force fields, AMBER99-ILDN, AMBER14SB, CHARMM22*, CHARMM36, and CHARMM36m, are the best candidates for studying amyloid peptide assembly, as they provide good balances in terms of structures and kinetics. We also investigated the assembly mechanisms of dimeric Aβ16-22 and found that the fibril formation rate is predominantly controlled by the total β-strand content.

Molecular Dynamics Simulation for Mechanism Elucidation of Food Processing and Safety: State of the Art

Molecular dynamics (MD) simulation is a useful technique to study the interaction between molecules and how they are affected by various processes and processing conditions. This review summarizes the application of MD simulations in food processing and safety, with an emphasis on the effects that emerging nonthermal technologies (for example, high hydrostatic pressure, pulsed electric field) have on the molecular and structural characteristics of foods and biomaterials. The advances and potential projection of MD simulations in the science and engineering aspects of food materials are discussed and focused on research work conducted to study the effects of emerging technologies on food components. It is expected by showing key case studies that it will stir novel developments as a valuable tool to study the effects of emerging food technologies on biomaterials. This review is useful to food researchers and the food industry, as well as researchers and practitioners working on flavor and nutraceutical encapsulations, dietary carbohydrate product developments, modified starches, protein engineering, and other novel food applications.