Replica exchange molecular dynamics simulations of amyloid peptide aggregation

Aggregation of Amyloidogenic Peptide Uperin-Molecular Dynamics Simulations

Uperin 3.5 is a remarkable natural peptide obtained from the skin of toadlets comprised of 17 amino acids which exhibits both antimicrobial and amyloidogenic properties. Molecular dynamics simulations were performed to study the β-aggregation process of uperin 3.5 as well as two of its mutants, in which the positively charged residues Arg7 and Lys8 have been replaced by alanine. All three peptides rapidly underwent spontaneous aggregation and conformational transition from random coils to beta-rich structures. The simulations reveal that the initial and essential step of the aggregation process involves peptide dimerization and the formation of small beta-sheets. A decrease in positive charge and an increase in the number of hydrophobic residues in the mutant peptides lead to an increase in the rate of their aggregation.

Influence of ALS-linked M337V mutation on the conformational ensembles of TDP-43321-340 peptide monomer and dimer

The transactive response (TAR) DNA/RNA-binding protein 43 (TDP-43) can self-assemble into both functional stress granules via liquid-liquid phase separation (LLPS) and pathogenic amyloid fibrillary aggregates that are closely linked to amyotrophic lateral sclerosis. Previous experimental studies reported that the low complexity domain (LCD) of TDP-43 plays an essential role in the LLPS and aggregation of the full-length protein, and it alone can also undergo LLPS to form liquid droplets mainly via intermolecular interactions in the 321-340 region. And the ALS-associated M337V mutation impairs LCD's LLPS and facilitates liquid-solid phase transition. However, the underlying atomistic mechanism is not well understood. Herein, as a first step to understand the M337V-caused LLPS disruption of TDP-43 LCD mediated by the 321-340 region and the fibrillization enhancement, we investigated the conformational properties of monomer/dimer of TDP-43321-340 peptide and its M337V mutant by performing extensive all-atom explicit-solvent replica exchange molecular dynamic simulations. Our simulations demonstrate that M337V mutation alters the residue regions with high helix/β-structure propensities and thus affects the conformational ensembles of both monomer and dimer. M337V mutation inhibits helix formation in the N-terminal Ala-rich region and the C-terminal mutation site region, while facilitating their long β-sheet formation, albeit with a minor impact on the average probability of both helix structure and β-structure. Further analysis of dimer system shows that M337V mutation disrupts inter-molecular helix-helix interactions and W334-W334 π-π stacking interactions which were reported to be important for the LLPS of TDP-43 LCD, whereas enhances the overall peptide residue-residue interactions and weakens peptide-water interactions, which is conducive to peptide fibrillization. This study provides mechanistic insights into the M337V-mutation-induced impairment of phase separation and facilitation of fibril formation of TDP-43 LCD.

Early-Stage Oligomerization of Prion-like Polypeptides Reveals the Molecular Mechanism of Amyloid-Disrupting Capacity by Proline Residues

Proline cis/trans isomerization governs protein local conformational changes via its local mechanical rigidity. The amyloid-disrupting capacity of proline is widely acknowledged; however, the molecular mechanism is still not clear. To understand how proline residues in polypeptide chains influence amyloid propensity, we study several truncated sequences of the TDP-43 C-terminal region (287-322) and their triple proline variants (308PPP310). We use coarse-grained molecular simulation to study the time evolution of the process of aggregation in the early stages in an effective high-concentration condition (∼25 mM). This ensures the long time scales for protein association at laboratory concentrations. We use several experimentally determined structure templates as initial structures of monomer conformations. We carry out oligomer size analysis and cluster analysis, along with several structural measures, to characterize the size distributions of oligomers and their morphological/structural properties. We show that average oligomer size is not a good indicator of amyloid propensity. Structural order and/or morphological properties are better alternatives. We show that proline variants can efficiently maintain the formation of large "ordered" oligomers of shorter truncated sequences, i.e., 307-322. This "order" maintenance is weakened when using longer truncated sequences (i.e., 287-322), leading to the formation of "disordered" oligomers. From an energy trade-off perspective, if the entropic effect is weak (short sequence length), the shape-complementarity of proline variants effectively guides the oligomerization process to form "ordered" oligomer intermediates. This leads to a distinct aggregation pathway that promotes amyloid formation (on-pathway). Strong entropic effects (long sequence length), however, would cause the formation of "disordered" oligomers. This in turn will suppress amyloid formation (off-pathway). The proline shape-complementary effects provide a guided morphological restraint to facilitate the pathways of amyloid formation. Our study supports the importance of structure-based kinetic heterogeneity of prion-like sequence fragments in driving different aggregation pathways. This work sheds light on the role of morphological and structural order of early-stage oligomeric species in regulating amyloid-disrupting capacity by prolines.

Virtual Quasi-2D Intermediates as Building Blocks for Plausible Structural Models of Amyloid Fibrils from Proteins with Complex Topologies: A Case Study of Insulin

Conformational transitions of globular proteins into amyloid fibrils are complex multistage processes exceedingly challenging to simulate using molecular dynamics (MD). Slow monomer diffusion rates and rugged free energy landscapes disfavor swift self-assembly of orderly amyloid architectures within timescales accessible to all-atom MD. Here, we conduct a multiscale MD study of the amyloidogenic self-assembly of insulin: a small protein with a complex topology defined by two polypeptide chains interlinked by three disulfide bonds. To avoid kinetic traps, unconventional preplanarized insulin conformations are used as amyloid building blocks. These starting conformers generated through uniaxial compression of the native monomer in various spatial directions represent 6 distinct (out of 16 conceivable) two-dimensional (2D) topological classes varying in N-/C-terminal segments of insulin's A- and B-chains being placed inside or outside of the central loop constituted by the middle sections of both chains and Cys7A-Cys7B/Cys19B-Cys20A disulfide bonds. Simulations of the fibrillar self-assembly are initiated through a biased in-register alignment of two, three, or four layers of flat conformers belonging to a single topological class. The various starting topologies are conserved throughout the self-assembly process resulting in polymorphic amyloid fibrils varying in structural features such as helical twist, presence of cavities, and overall stability. Some of the protofilament structures obtained in this work are highly compatible with the earlier biophysical studies on insulin amyloid and high-resolution studies on insulin-derived amyloidogenic peptide models postulating the presence of steric zippers. Our approach provides in silico means to study amyloidogenic tendencies and viable amyloid architectures of larger disulfide-constrained proteins with complex topologies.

Influencing factors and characterization methods of nanoparticles regulating amyloid aggregation

Human disorders associated with amyloid aggregation, such as Alzheimer's disease and Parkinson's disease, afflict the lives of millions worldwide. When peptides and proteins in the body are converted to amyloids, which have a tendency to aggregate, the toxic oligomers produced during the aggregation process can trigger a range of diseases. Nanoparticles (NPs) have been found to possess surface effects that can modulate the amyloid aggregation process and they have potential application value in the treatment of diseases related to amyloid aggregation and fibrillary tangles. In this review, we discuss recent progress relating to studies of nanoparticles that regulate amyloid aggregation. The review focuses on the factors influencing this regulation, which are important as guidelines for the future design of NPs for the treatment of amyloid aggregation. We describe the characterization methods that have been utilized so far in such studies. This review provides research information and characterization methods for the rational design of NPs, which should result in therapeutic strategies for amyloid diseases.

Fibril Surface-Dependent Amyloid Precursors Revealed by Coarse-Grained Molecular Dynamics Simulation

Amyloid peptides are known to self-assemble into larger aggregates that are linked to the pathogenesis of many neurodegenerative disorders. In contrast to primary nucleation, recent experimental and theoretical studies have shown that many toxic oligomeric species are generated through secondary processes on a pre-existing fibrillar surface. Nucleation, for example, can also occur along the surface of a pre-existing fibril—secondary nucleation—as opposed to the primary one. However, explicit pathways are still not clear. In this study, we use molecular dynamics simulation to explore the free energy landscape of a free Abeta monomer binding to an existing fibrillar surface. We specifically look into several potential Abeta structural precursors that might precede some secondary events, including elongation and secondary nucleation. We find that the overall process of surface-dependent events can be described at least by the following three stages: 1. Free diffusion 2. Downhill guiding 3. Dock and lock. And we show that the outcome of adding a new monomer onto a pre-existing fibril is pathway-dependent, which leads to different secondary processes. To understand structural details, we have identified several monomeric amyloid precursors over the fibrillar surfaces and characterize their heterogeneity using a probability contact map analysis. Using the frustration analysis (a bioinformatics tool), we show that surface heterogeneity correlates with the energy frustration of specific local residues that form binding sites on the fibrillar structure. We further investigate the helical twisting of protofilaments of different sizes and observe a length dependence on the filament twisting. This work presents a comprehensive survey over the properties of fibril growth using a combination of several openMM-based platforms, including the GPU-enabled openAWSEM package for coarse-grained modeling, MDTraj for trajectory analysis, and pyEMMA for free energy calculation. This combined approach makes long-timescale simulation for aggregation systems as well as all-in-one analysis feasible. We show that this protocol allows us to explore fibril stability, surface binding affinity/heterogeneity, as well as fibrillar twisting. All these properties are important for understanding the molecular mechanism of surface-catalyzed secondary processes of fibril growth.

On the Emergence of Orientational Order in Folded Proteins with Implications for Allostery

The beautiful structures of single- and multi-domain proteins are clearly ordered in some fashion but cannot be readily classified using group theory methods that are successfully used to describe periodic crystals. For this reason, protein structures are considered to be aperiodic, and may have evolved this way for functional purposes, especially in instances that require a combination of softness and rigidity within the same molecule. By analyzing the solved protein structures, we show that orientational symmetry is broken in the aperiodic arrangement of the secondary structure elements (SSEs), which we deduce by calculating the nematic order parameter, P2. We find that the folded structures are nematic droplets with a broad distribution of P2. We argue that a non-zero value of P2, leads to an arrangement of the SSEs that can resist external forces, which is a requirement for allosteric proteins. Such proteins, which resist mechanical forces in some regions while being flexible in others, transmit signals from one region of the protein to another (action at a distance) in response to binding of ligands (oxygen, ATP, or other small molecules).

Disaggregation mechanism of prion amyloid for tweezer inhibitor

The aggregation of amyloid has been an important event in the pathology of amyloidogenicity. A number of small molecules have been designed for Amyloidosis treatment. Molecular tweezer CLR01, a potential drug for misfolded β-amyloids inhibition, was reportedly bind directly to Lysine residues and interrupt oligomerization. However, the disaggregation mechanism of amyloid for this inhibitor is unclear. Here we used long timescale of molecular dynamic simulation to reveal the mechanism of disaggregation for pentamer prion amyloid. Molecular docking and molecular dynamics simulation demonstrate that CLR01 is attached with Lysine₂₂₂ nitrogen by π-cation interaction of its nine aromatic rings and formation of salt bridge/hydrogen bond of one of the two rotatable peripheral anionic phosphate groups. Upon CLR01 binding, we found a major shifting occurs in initial conformation of the oligomer and stretch out the N-terminal chain A from the rest of the amyloid which seems to be the first stage of disaggregated the fibrils slowly yet efficiently. Moreover, the CLR01 remodelled the pentamer Prion_220-272 into a compact structure which might be the resistant conformation for further oligomerization. Our work will contribute to better understand the interaction and deterioration mechanism of molecular tweezer for prions and similar amyloids, and offer significant insights into therapeutic development for Amyloidosis treatment.

Clustering and Fibril Formation during GNNQQNY Aggregation: A Molecular Dynamics Study

The precise kinetic pathways of peptide clustering and fibril formation are not fully understood. Here we study the initial clustering kinetics and transient cluster morphologies during aggregation of the heptapeptide fragment GNNQQNY from the yeast prion protein Sup35. We use a mid-resolution coarse-grained molecular dynamics model of Bereau and Deserno to explore the aggregation pathways from the initial random distribution of free monomers to the formation of large clusters. By increasing the system size to 72 peptides we could follow directly the molecular events leading to the formation of stable fibril-like structures. To quantify those structures we developed a new cluster helicity parameter. We found that the formation of fibril-like structures is a cooperative processes that requires a critical number of monomers, M⋆≈25, in a cluster. The terminal tyrosine residue is the structural determinant in the formation of helical fibril-like structures. This work supports and quantifies the two-step aggregation model where the initially formed amorphous clusters grow and, when they are large enough, rearrange into mature twisted structures. However, in addition to the nucleated fibrillation, growing aggregates undergo further internal reorganization, which leads to more compact structures of large aggregates.

Latent Models of Molecular Dynamics Data: Automatic Order Parameter Generation for Peptide Fibrillization

Variational autoencoders are artificial neural networks with the capability to reduce highly dimensional sets of data to smaller dimensional, latent representations. In this work, these models are applied to molecular dynamics simulations of the self-assembly of coarse-grained peptides to obtain a singled-valued order parameter for amyloid aggregation. This automatically learned order parameter is constructed by time-averaging the latent parametrizations of internal coordinate representations and compared to the nematic order parameter which is commonly used to study ordering of similar systems in literature. It is found that the latent space value provides more tailored insight into the aggregation mechanism's details, correctly identifying fibril formation in instances where the nematic order parameter fails to do so. A means is provided by which the latent space value can be analyzed so that the major contributing internal coordinates are identified, allowing for a direct interpretation of the latent space order parameter in terms of the behavior of the system. The latent model is found to be an effective and convenient way of representing the data from the dynamic ensemble and provides a means of reducing the dimensionality of a system whose scale exceeds molecular systems so-far considered with similar tools. This bypasses a need for researcher speculation on what elements of a system best contribute to summarizing major transitions and suggests latent models are effective and insightful when applied to large systems with a diversity of complex behaviors.

Designer Amyloid Cell-Penetrating Peptides for Potential Use as Gene Transfer Vehicles

Cell-penetrating peptides are used extensively to deliver molecules into cells due to their unique characteristics such as rapid internalization, charge, and non-cytotoxicity. Amyloid fibril biomaterials were reported as gene transfer or retroviral infection enhancers; no cell internalization of the peptides themselves is reported so far. In this study, we focus on two rationally and computationally designed peptides comprised of β-sheet cores derived from naturally occurring protein sequences and designed positively charged and aromatic residues exposed at key residue positions. The β-sheet cores bestow the designed peptides with the ability to self-assemble into amyloid fibrils. The introduction of positively charged and aromatic residues additionally promotes DNA condensation and cell internalization by the self-assembled material formed by the designed peptides. Our results demonstrate that these designer peptide fibrils can efficiently enter mammalian cells while carrying packaged luciferase-encoding plasmid DNA, and they can act as a protein expression enhancer. Interestingly, the peptides additionally exhibited strong antimicrobial activity against the enterobacterium Escherichia coli.

Binding Process and Free Energy Characteristics of Cellulose Chain into the Catalytic Domain of Cellobiohydrolase TrCel7A

It was observed in experiments that the catalytic domain (CD) of Trichoderma reesei Cel7A (TrCel7A) hydrolyzes crystalline cellulose in a processive manner, but the underlying binding mechanism is still unknown. Here, through replica-exchange molecular dynamics simulations, we find that the loading and sucking-in process of the cellulose chain into CD is entropy-driven and enthalpy-unfavorable, which firmly relate to the desolvation of the binding channel of CD. During the loading process, hydrophobic interactions play a dominant role because several aromatic residues have been identified to guide the cellulose chain processing. At the active site, a transition from enthalpy- to entropy-driven is detected for the driving force. Such a finding reveals the indispensability of the catalytic reaction of the glycosidic bond to provide the energy to drive the movements of the cellulose chain. Our study reveals the interaction pictures between the cellulose chain and TrCel7A at the atomic level, which helps better understand the catalytic mechanism of TrCel7A.

Molecular Dynamics Simulations Reveal the Inhibitory Mechanism of Dopamine against Human Islet Amyloid Polypeptide (hIAPP) Aggregation and Its Destabilization Effect on hIAPP Protofibrils

The aberrant self-assembly of human islet amyloid polypeptide (hIAPP) into toxic oligomers, protofibrils, and mature fibrils is associated with the pathogenesis of type 2 diabetes (T2D). Inhibition of hIAPP aggregation and destabilization of preformed hIAPP fibrils are considered as two major therapeutic strategies for treating T2D. Previous experimental studies reported that dopamine prevented the formation of hIAPP oligomers and fibrils. However, the underlying inhibitory mechanism at the atomic level remains elusive. Herein we investigated the conformational ensembles of hIAPP dimer with and without dopamine using replica-exchange molecular dynamics simulations. The simulations demonstrated that dopamine preferentially bound to R11, L12, F15, H18, F23, I26, L27, and Y37 residues, inhibited the formation of β-sheets in the amyloidogenic regions spanning residues ¹¹RLANFLVH¹⁸, ²²NFGAIL²⁷, and ³⁰TNVGSNT³⁶, and resulted in more disordered hIAPP dimers, thus hindering the amyloid formation of hIAPP. Protonated and deprotonated dopamine molecules displayed distinct binding capabilities but bound to similar residue sites on hIAPP. Additional microsecond molecular dynamics simulations showed that dopamine mainly bound to the β1 and turn regions of hIAPP protofibril and destabilized the protofibril structure. This study not only revealed the molecular mechanism of dopamine toward the inhibition of hIAPP aggregation but also demonstrated the protofibril-destabilizing effects of dopamine, which may be helpful for the design of drug candidates to treat T2D.

In Silico Prediction of the Thermodynamic Equilibrium of Solute Partition in Multiphase Complex Fluids: A Case Study of Oil-Water Microemulsion

Multiphase complex fluids such as micelles, microemulsions, and dispersions are ubiquitous in product formulations of foods, pharmaceuticals, cosmetics, and fine chemicals. Quantifying how active solutes partition in the microstructure of such multiphase fluids is necessary for designing formulations that can optimally deliver the benefits of functional actives. In this paper, we at first predict the structure of a heptane/butanol/sodium dodecyl sulfate droplet in water that self-assembled to form a microemulsion through the molecular dynamics (MD) simulation and subsequently investigate the thermodynamic equilibrium of solute partitioning using COSMOmic. To our knowledge, this is the first time that the MD/COSMOmic approach is used for predicting solute partitioning in a microemulsion. The predicted partition coefficients are compared to experimental values derived from retention measurements of the same microemulsion. We show that the experimental data of droplet-water partition coefficients (K_droplet/w) can be reliably predicted by the method that combines MD simulations with COSMOmic.

Aggregation kinetics of short peptides: All-atom and coarse-grained molecular dynamics study

Peptides can aggregate into ordered structures with different morphologies. The aggregation mechanism and evolving structures are the subject of intense research. In this paper we have used molecular dynamics to examine the sequence-dependence of aggregation kinetics for three short peptides: octaalanine (Ala8), octaasparagine (Asn8), and the heptapeptide GNNQQNY (abbreviated as GNN). First, we compared the aggregation of 20 randomly distributed peptides using the coarse-grained MARTINI force field and the atomistic OPLS-AA force field. We found that the MARTINI and OPLS-AA aggregation kinetics are similar for Ala8, Asn8, and GNN. Second, we used the MARTINI force field to study the early stages of aggregation kinetics for a larger system with 72 peptides. In the initial stage of aggregation small clusters grow by monomer addition. In the second stage, when the free monomers are depleted, the dominant cluster growth path is cluster-cluster coalescence. We quantified the aggregation kinetics in terms of rate equations. Our study shows that the initial aggregation kinetics are similar for Ala8, Asn8, and GNN but the molecular details can be different, especially for MARTINI Ala8. We hypothesize that peptide aggregation proceed in two steps. In the first step amorphous aggregates are formed, and then, in the second step, they reorganize into ordered structures. We conclude that sequence-specific differences show up in the second step of aggregation.

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.

A Biomphalaria glabrata peptide that stimulates significant behaviour modifications in aquatic free-living Schistosoma mansoni miracidia

The human disease schistosomiasis (or bilharzia) is caused by the helminth blood fluke parasite Schistosoma mansoni, which requires an intermediate host, the freshwater gastropod snail Biomphalaria glabrata (the most common intermediate host). The free-swimming parasite miracidia utilise an excellent chemosensory sense to detect and locate an appropriate host. This study investigated the biomolecules released by the snail that stimulate changes in the behaviour of the aquatic S. mansoni miracidia. To achieve this, we have performed an integrated analysis of the snail-conditioned water, through chromatography and bioassay-guided behaviour observations, followed by mass spectrometry. A single fraction containing multiple putative peptides could stimulate extreme swimming behaviour modifications (e.g. velocity, angular variation) similar to those observed in response to crude snail mucus. One peptide (P12;—R-DITSGLDPEVADD-KR—) could replicate the stimulation of miracidia behaviour changes. P12 is derived from a larger precursor protein with a signal peptide and multiple dibasic cleavage sites, which is synthesised in various tissues of the snail, including the central nervous system and foot. P12 consists of an alpha helix secondary structure as indicated by circular dichroism spectroscopy. This information will be helpful for the development of approaches to manipulate this parasites life cycle, and opens up new avenues for exploring other parasitic diseases which have an aquatic phase using methods detailed in this investigation.

In aquatic environments, where the vast majority of animals live in darkness, key relationships are often formed and maintained by chemical communication (including smell and taste). Parasites with an aquatic life phase rely on an exquisite sense of chemosensation to detect host biomolecules (kairomones), allowing them to locate and infect their host. Our study identifies the first kairomone released by the freshwater gastropod snail Biomphalaria glabrata, an intermediate host for the helminth blood fluke parasite Schistosoma mansoni. This is a key aspect of the S. mansoni life-cycle that ultimately leads to human infection, causing the disease schistosomiasis (or bilharzia), which is considered the most devastating human helminth infection in terms of global morbidity and mortality. The kairomone we identify is a peptide that does not appear to share any similarity with any other known animal peptide. This information will be helpful as we explore methods to interrupt parasite infection, and therefore break the cycle of infection that causes a major human disease.

Entropic forces drive clustering and spatial localization of influenza A M2 during viral budding

The influenza A matrix 2 (M2) transmembrane protein facilitates virion release from the infected host cell. In particular, M2 plays a role in the induction of membrane curvature and/or in the scission process whereby the envelope is cut upon virion release. Here we show using coarse-grained computer simulations that various M2 assembly geometries emerge due to an entropic driving force, resulting in compact clusters or linearly extended aggregates as a direct consequence of the lateral membrane stresses. Conditions under which these protein assemblies will cause the lipid membrane to curve are explored, and we predict that a critical cluster size is required for this to happen. We go on to demonstrate that under the stress conditions taking place in the cellular membrane as it undergoes large-scale membrane remodeling, the M2 protein will, in principle, be able to both contribute to curvature induction and sense curvature to line up in manifolds where local membrane line tension is high. M2 is found to exhibit linactant behavior in liquid-disordered-liquid-ordered phase-separated lipid mixtures and to be excluded from the liquid-ordered phase, in near-quantitative agreement with experimental observations. Our findings support a role for M2 in membrane remodeling during influenza viral budding both as an inducer and a sensor of membrane curvature, and they suggest a mechanism by which localization of M2 can occur as the virion assembles and releases from the host cell, independent of how the membrane curvature is produced.

A novel amyloid designable scaffold and potential inhibitor inspired by GAIIG of amyloid beta and the HIV-1 V3 loop

The GAIIG sequence, common to the amyloid beta peptide (residues 29-33) and to the HIV-1 gp120 (residues 24-28 in a typical V3 loop), self-assembles into amyloid fibrils, as suggested by theory and the experiments presented here. The longer YATGAIIGNII sequence from the V3 loop also self-assembles into amyloid fibrils, of which the first three and the last two residues are outside the amyloid GAIIG core. We postulate that this sequence, with suitably selected modifications at the flexible positions, can serve as a designable scaffold for novel amyloid-based materials. Moreover, we report the single crystal X-ray structure of the beta-breaker peptide GAIPIG at 1.05 Å resolution. The structural information provided in this study could serve as the basis for structure-based design of potential inhibitors of amyloid formation.

Inhibition of GNNQQNY prion peptide aggregation by trehalose: a mechanistic view

Deposition of amyloid fibrils is the seminal event in the pathogenesis of numerous neurodegenerative diseases. The formation of this amyloid assembly is the manifestation of a cascade of structural transitions including toxic oligomer formation in the early stages of aggregation. Thus a viable therapeutic strategy involves the use of small molecular ligands to interfere with this assembly. In this perspective, we have explored the kinetics of aggregate formation of the fibril forming GNNQQNY peptide fragment from the yeast prion protein SUP35 using multiple all atom MD simulations with explicit solvent and provided mechanistic insights into the way trehalose, an experimentally known aggregation inhibitor, modulates the aggregation pathway. The results suggest that the assimilation process is impeded by different barriers at smaller and larger oligomeric sizes: the initial one being easily surpassed at higher temperatures and peptide concentrations. The kinetic profile demonstrates that trehalose delays the aggregation process by increasing both these activation barriers, specifically the latter one. It increases the sampling of small-sized aggregates that lack the beta sheet conformation. Analysis reveals that the barrier in the growth of larger stable oligomers causes the formation of multiple stable small oligomers which then fuse together bimolecularly. The PCA of 26 properties was carried out to deconvolute the events within the temporary lag phases, which suggested dynamism in lags involving an increase in interchain contacts and burial of SASA. The predominant growth route is monomer addition, which changes to condensation on account of a large number of depolymerisation events in the presence of trehalose. The favourable interaction of trehalose specifically with the sidechain of the peptide promotes crowding of trehalose molecules in its vicinity - the combination of both these factors imparts the observed behaviour. Furthermore, increasing trehalose concentration leads to faster expulsion of water molecules than interpeptide interactions. These expelled water molecules have larger translational movement, suggesting an entropy factor to favor the assembly process. Different conformations observed under this condition suggest the role of water molecules in guiding the morphology of the aggregates as well. A similar scenario exists on increasing peptide concentration.

Insights into the Molecular Mechanisms of Alzheimer's and Parkinson's Diseases with Molecular Simulations: Understanding the Roles of Artificial and Pathological Missense Mutations in Intrinsically Disordered Proteins Related to Pathology

Amyloid-β and α-synuclein are intrinsically disordered proteins (IDPs), which are at the center of Alzheimer's and Parkinson's disease pathologies, respectively. These IDPs are extremely flexible and do not adopt stable structures. Furthermore, both amyloid-β and α-synuclein can form toxic oligomers, amyloid fibrils and other type of aggregates in Alzheimer's and Parkinson's diseases. Experimentalists face challenges in investigating the structures and thermodynamic properties of these IDPs in their monomeric and oligomeric forms due to the rapid conformational changes, fast aggregation processes and strong solvent effects. Classical molecular dynamics simulations complement experiments and provide structural information at the atomic level with dynamics without facing the same experimental limitations. Artificial missense mutations are employed experimentally and computationally for providing insights into the structure-function relationships of amyloid-β and α-synuclein in relation to the pathologies of Alzheimer's and Parkinson's diseases. Furthermore, there are several natural genetic variations that play a role in the pathogenesis of familial cases of Alzheimer's and Parkinson's diseases, which are related to specific genetic defects inherited in dominant or recessive patterns. The present review summarizes the current understanding of monomeric and oligomeric forms of amyloid-β and α-synuclein, as well as the impacts of artificial and pathological missense mutations on the structural ensembles of these IDPs using molecular dynamics simulations. We also emphasize the recent investigations on residual secondary structure formation in dynamic conformational ensembles of amyloid-β and α-synuclein, such as β-structure linked to the oligomerization and fibrillation mechanisms related to the pathologies of Alzheimer's and Parkinson's diseases. This information represents an important foundation for the successful and efficient drug design studies.

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

The conformational space and structural ensembles of amyloid beta (Aβ) peptides and their oligomers in solution are inherently disordered and proven to be challenging to study. Optimum force field selection for molecular dynamics (MD) simulations and the biophysical relevance of results are still unknown. We compared the conformational space of the Aβ(1-40) dimers by 300 ns replica exchange MD simulations at physiological temperature (310 K) using: the AMBER-ff99sb-ILDN, AMBER-ff99sb*-ILDN, AMBER-ff99sb-NMR, and CHARMM22* force fields. Statistical comparisons of simulation results to experimental data and previously published simulations utilizing the CHARMM22* and CHARMM36 force fields were performed. All force fields yield sampled ensembles of conformations with collision cross sectional areas for the dimer that are statistically significantly larger than experimental results. All force fields, with the exception of AMBER-ff99sb-ILDN (8.8 ± 6.4%) and CHARMM36 (2.7 ± 4.2%), tend to overestimate the α-helical content compared to experimental CD (5.3 ± 5.2%). Using the AMBER-ff99sb-NMR force field resulted in the greatest degree of variance (41.3 ± 12.9%). Except for the AMBER-ff99sb-NMR force field, the others tended to under estimate the expected amount of β-sheet and over estimate the amount of turn/bend/random coil conformations. All force fields, with the exception AMBER-ff99sb-NMR, reproduce a theoretically expected β-sheet-turn-β-sheet conformational motif, however, only the CHARMM22* and CHARMM36 force fields yield results compatible with collapse of the central and C-terminal hydrophobic cores from residues 17-21 and 30-36. Although analyses of essential subspace sampling showed only minor variations between force fields, secondary structures of lowest energy conformers are different.

Structural properties of amyloid β(1-40) dimer explored by replica exchange molecular dynamics simulations

Replica exchange molecular dynamics simulations (300 ns) were used to study the dimerization of amyloid β(1-40) (Aβ(1-40)) polypeptide. Configurational entropy calculations revealed that at physiological temperature (310 K, 37°C) dynamic dimers are formed by randomly docked monomers. Free energy of binding of the two chains to each other was -93.56 ± 6.341 kJ mol^-1 . Prevalence of random coil conformations was found for both chains with the exceptions of increased β-sheet content from residues 16-21 and 29-32 of chain A and residues 15-21 and 30-33 of chain B with β-turn/β-bend conformations in both chains from residues 1-16, 21-29 of chain A, 1-16, and 21-29 of chain B. There is a mixed β-turn/β-sheet region from residues 33-38 of both chains. Analysis of intra- and interchain residue distances shows that, although the individual chains are highly flexible, the dimer system stays in a loosely packed antiparallel β-sheet configuration with contacts between residues 17-21 of chain A with residues 17-21 and 31-36 of chain B as well as residues 31-36 of chain A with residues 17-21 and 31-36 of chain B. Based on dihedral principal component analysis, the antiparallel β-sheet-loop-β-sheet conformational motif is favored for many low energy sampled conformations. Our results show that Aβ(1-40) can form dynamic dimers in aqueous solution that have significant conformational flexibility and are stabilized by collapse of the central and C-terminal hydrophobic cores with the expected β-sheet-loop-β-sheet conformational motif. Proteins 2017; 85:1024-1045. © 2017 Wiley Periodicals, Inc.

An Acrobatic Substrate Metamorphosis Reveals a Requirement for Substrate Conformational Dynamics in Trypsin Proteolysis

The molecular basis of enzyme catalytic power and specificity derives from dynamic interactions between enzyme and substrate during catalysis. Although considerable effort has been devoted to understanding how conformational dynamics within enzymes affect catalysis, the role of conformational dynamics within protein substrates has not been addressed. Here, we examine the importance of substrate dynamics in the cleavage of Kunitz-bovine pancreatic trypsin inhibitor protease inhibitors by mesotrypsin, finding that the varied conformational dynamics of structurally similar substrates can profoundly impact the rate of catalysis. A 1.4-Å crystal structure of a mesotrypsin-product complex formed with a rapidly cleaved substrate reveals a dramatic conformational change in the substrate upon proteolysis. By using long all-atom molecular dynamics simulations of acyl-enzyme intermediates with proteolysis rates spanning 3 orders of magnitude, we identify global and local dynamic features of substrates on the nanosecond-microsecond time scale that correlate with enzymatic rates and explain differential susceptibility to proteolysis. By integrating multiple enhanced sampling methods for molecular dynamics, we model a viable conformational pathway between substrate-like and product-like states, linking substrate dynamics on the nanosecond-microsecond time scale with large collective substrate motions on the much slower time scale of catalysis. Our findings implicate substrate flexibility as a critical determinant of catalysis.

Prediction of a stable associated liquid of short amyloidogenic peptides

Amyloid fibril formation is believed to be a nucleation-controlled process. Depending on the nature of peptide sequence, fibril nucleation can occur in one step, straight from a dilute solution, or in multiple steps via oligomers or disordered aggregates. What determines this process is poorly understood. Since the fibril formation kinetics is driven by thermodynamic forces, knowledge of the phase behavior is crucial. Here, we investigated the phase behavior of three short peptide sequences of varying side-chain hydrophobicity. Replica exchange molecular dynamics simulations of a mid-resolution model indicate that the weakly hydrophobic peptide forms fibrils directly from solution, whereas the most hydrophobic peptide forms a dense liquid phase before crystallizing into ordered fibrils at low temperatures. For the medium hydrophobic peptide we found evidence of a novel additional transition to a liquid phase consisting of clusters of aligned peptides, implying a three-step nucleation process. We tested the robustness of this prediction by applying Wertheim's theory and statistical associating fluid theory to a hard-sphere model dressed with isotropic and anisotropic attractions. We found that the ratio of interaction strengths strongly affects the phase behavior, and under certain conditions indeed gives rise to a stable polymerized liquid phase. The peptide clusters in the associated liquid tend to be slow and long-lived, which may give the oligomer droplet more time to act as a toxic oligomer, before turning into a fibril.

Preformed template fluctuations promote fibril formation: insights from lattice and all-atom models

Fibril formation resulting from protein misfolding and aggregation is a hallmark of several neurodegenerative diseases such as Alzheimer's and Parkinson's diseases. Despite the fact that the fibril formation process is very slow and thus poses a significant challenge for theoretical and experimental studies, a number of alternative pictures of molecular mechanisms of amyloid fibril formation have been recently proposed. What seems to be common for the majority of the proposed models is that fibril elongation involves the formation of pre-nucleus seeds prior to the creation of a critical nucleus. Once the size of the pre-nucleus seed reaches the critical nucleus size, its thermal fluctuations are expected to be small and the resulting nucleus provides a template for sequential (one-by-one) accommodation of added monomers. The effect of template fluctuations on fibril formation rates has not been explored either experimentally or theoretically so far. In this paper, we make the first attempt at solving this problem by two sets of simulations. To mimic small template fluctuations, in one set, monomers of the preformed template are kept fixed, while in the other set they are allowed to fluctuate. The kinetics of addition of a new peptide onto the template is explored using all-atom simulations with explicit water and the GROMOS96 43a1 force field and simple lattice models. Our result demonstrates that preformed template fluctuations can modulate protein aggregation rates and pathways. The association of a nascent monomer with the template obeys the kinetics partitioning mechanism where the intermediate state occurs in a fraction of routes to the protofibril. It was shown that template immobility greatly increases the time of incorporating a new peptide into the preformed template compared to the fluctuating template case. This observation has also been confirmed by simulation using lattice models and may be invoked to understand the role of template fluctuations in slowing down fibril elongation in vivo.

Protective V127 prion variant prevents prion disease by interrupting the formation of dimer and fibril from molecular dynamics simulations

Recent studies uncovered a novel protective prion protein variant: V127 variant, which was reported intrinsically resistant to prion conversion and propagation. However, the structural basis of its protective effect is still unknown. To uncover the origin of the protective role of V127 variant, molecular dynamics simulations were performed to explore the influence of G127V mutation on two key processes of prion propagation: dimerization and fibril formation. The simulation results indicate V127 variant is unfavorable to form dimer by reducing the main-chain H-bond interactions. The simulations of formed fibrils consisting of β1 strand prove V127 variant will make the formed fibril become unstable and disorder. The weaker interaction energies between layers and reduced H-bonds number for V127 variant reveal this mutation is unfavorable to the formation of stable fibril. Consequently, we find V127 variant is not only unfavorable to the formation of dimer but also unfavorable to the formation of stable core and fibril, which can explain the mechanism on the protective role of V127 variant from the molecular level. Our findings can deepen the understanding of prion disease and may guide the design of peptide mimetics or small molecule to mimic the protective effect of V127 variant.

A minimal conformational switching-dependent model for amyloid self-assembly

Amyloid formation is associated with various pathophysiological conditions like Alzheimer’s and Parkinson’s diseases as well as many useful functions. The hallmark of amyloid assemblies is a conformational transition of the constituent proteins into a β - sheet rich filament. Accounting for this conformational transition in amyloidogenic proteins, we develop an analytically solvable model that can probe the dynamics of an ensemble of single filaments. Using the theory and Monte Carlo simulations, we show the presence of two kinetic regimes for the growth of a self-assembling filament – switching-dependent and –independent growth regimes. We observe a saturation in fibril elongation velocities at higher concentrations in the first regime, providing a novel explanation to the concentration-independence of growth velocities observed experimentally. We also compute the length fluctuation of the filaments to characterize aggregate heterogeneity. From the early velocities and length fluctuation, we propose a novel way of estimating the conformational switching rate. Our theory predicts a kinetic phase diagram that has three distinct phases – short oligomers/monomers, disordered aggregates and β -rich filaments. The model also predicts the force generation potential and the intermittent growth of amyloid fibrils evident from single molecular experiments. Our model could contribute significantly to the physical understanding of amyloid aggregation.

Dynamics of Seeded Aβ40-Fibril Growth from Atomistic Molecular Dynamics Simulations: Kinetic Trapping and Reduced Water Mobility in the Locking Step

Filamentous β-amyloid aggregates are crucial for the pathology of Alzheimer's disease. Despite the tremendous biomedical importance, the molecular pathway of growth propagation is not completely understood and remains challenging to investigate by simulations due to the long time scales involved. Here, we apply extensive all-atom molecular dynamics simulations in explicit water to obtain free energy profiles and kinetic information from position-dependent diffusion profiles for three different Aβ9-40-growth processes: fibril elongation by single monomers at the structurally unequal filament tips and association of larger filament fragments. Our approach provides insight into the molecular steps of the kinetic pathway and allows close agreement with experimental binding free energies and macroscopic growth rates. Water plays a decisive role, and solvent entropy is identified as the main driving force for assembly. Fibril growth is disfavored energetically due to cancellation of direct peptide-peptide interactions and solvation effects. The kinetics of growth is consistent with the characteristic dock/lock mechanism, and docking is at least 2 orders of magnitude faster. During initial docking, interactions are mediated by transient non-native hydrogen bonds, which efficiently catch the incoming monomer or fragment already at separations of about 3 nm. In subsequent locking, the dynamics is much slower due to formation of kinetically trapped conformations caused by long-lived non-native hydrogen bonds. Fibril growth additionally requires collective motion of water molecules to create a dry binding interface. Fibril growth is further retarded due to reduced mobility of the involved hydration water, evident from a 2-fold reduction of the diffusion coefficient.

Interplay between the hydrophobic effect and dipole interactions in peptide aggregation at interfaces

Amphipathic octapeptide aggregation at hydrophobic–hydrophilic interfaces is largely driven by backbone dipole interactions in peptide aggregation at interfaces.

The structure of chromophore-grafted amyloid-β12–28 dimers in the gas-phase: FRET-experiment guided modelling

Theoretical modelling, ion mobility spectrometry and action-FRET experiments are combined to an experiment guided approach and used to elucidate the structure of chromophore-grafted amyloid-β_12–28 dimers in the gas-phase.

Simulated Solute Tempering

For the enhanced conformational sampling in molecular dynamics (MD) simulations, we present "simulated solute tempering" (SST) which is an easy to implement variant of simulated tempering. SST extends conventional simulated tempering (CST) by key concepts of "replica exchange with solute tempering" (REST, Liu et al. Proc. Natl. Acad. Sci. U.S.A. 2005, 102, 13749). We have applied SST, CST, and REST to molecular dynamics (MD) simulations of an alanine octapeptide in explicit water. The weight parameters required for CST and SST are determined by two different formulas whose performance is compared. For SST only one of them yields a uniform sampling of the temperature space. Compared to CST and REST, SST provides the highest exchange probabilities between neighboring rungs in the temperature ladder. Concomitantly, SST leads to the fastest diffusion of the simulation system through the temperature space, in particular, if the "even-odd" exchange scheme is employed in SST. As a result, SST exhibits the highest sampling speed of the investigated tempering methods.

Membrane Interaction of the Factor VIIIa Discoidin Domains in Atomistic Detail

A recently developed membrane-mimetic model was applied to study membrane interaction and binding of the two anchoring C2-like discoidin domains of human coagulation factor VIIIa (FVIIIa), the C1 and C2 domains. Both individual domains, FVIII C1 and FVIII C2, were observed to bind the phospholipid membrane by partial or full insertion of their extruding loops (the spikes). However, the two domains adopted different molecular orientations in their membrane-bound states; FVIII C2 roughly was positioned normal to the membrane plane, while FVIII C1 displayed a multitude of tilted orientations. The results indicate that FVIII C1 may be important in modulating the orientation of the FVIIIa molecule to optimize the interaction with FIXa, which is anchored to the membrane via its γ-carboxyglutamic acid-rich (Gla) domain. Additionally, a structural change was observed in FVIII C1 in the coiled main chain leading the first spike. A tight interaction with one lipid per domain, similar to what has been suggested for the homologous FVa C2, is characterized. Finally, we rationalize known FVIII antibody epitopes and the scarcity of documented hemophilic missense mutations related to improper membrane binding of FVIIIa, based on the prevalent nonspecificity of ionic interactions in the simulated membrane-bound states of FVIII C1 and FVIII C2.

Greedy replica exchange algorithm for heterogeneous computing grids

Replica exchange molecular dynamics (REMD) has become a valuable tool in studying complex biomolecular systems. However, its application on distributed computing grids is limited by the heterogeneity of this environment. In this study, we propose a REMD implementation referred to as greedy REMD (gREMD) suitable for computations on heterogeneous grids. To decentralize replica management, gREMD utilizes a precomputed schedule of exchange attempts between temperatures. Our comparison of gREMD against standard REMD suggests four main conclusions. First, gREMD accelerates grid REMD simulations by as much as 40 %. Second, gREMD increases CPU utilization rates in grid REMD by up to 60 %. Third, we argue that gREMD is expected to maintain approximately constant CPU utilization rates and simulation wall-clock times with the increase in the number of replicas. Finally, we show that gREMD correctly implements the REMD algorithm and reproduces the conformational ensemble of a short peptide sampled in our previous standard REMD simulations. We believe that gREMD can find its place in large-scale REMD simulations on heterogeneous computing grids.

Flexible receptor docking for drug discovery

INTRODUCTION

Molecular docking has become a popular method for virtual screening. Docking small molecules to a rigid biological receptor is fast but could produce many false negatives and identify less diverse compounds. Flexible receptor docking has alleviated this problem.

AREAS COVERED

This article focuses on reviewing ensemble docking as an approximate but inexpensive method to incorporate receptor flexibility in molecular docking. It outlines key features and recent advances of this method and points out problem areas that need to be addressed to make it even more useful in drug discovery.

EXPERT OPINION

Among the different methods introduced for flexible receptor docking, ensemble docking represents one of the most popular approaches, especially for high-throughput virtual screening. One can generate structural ensembles by using experimental structures, by structural modeling and by various types of molecular simulations. In building a structural ensemble, a judicious choice of the structures to be included can improve performance. Furthermore, reducing the size of the structural ensemble can cut computational costs, and removing the structures that can bind few ligands well could enrich the number of true actives identified by ensemble docking. The ability of ensemble docking to identify more true positives at the top of a rank-ordered list also depends on the choice of the methods to score and rank compounds, an area that needs further research.

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.

Exploring the role of hydration and confinement in the aggregation of amyloidogenic peptides Aβ(16-22) and Sup35(7-13) in AOT reverse micelles

Knowledge of how intermolecular interactions of amyloid-forming proteins cause protein aggregation and how those interactions are affected by sequence and solution conditions is essential to our understanding of the onset of many degenerative diseases. Of particular interest is the aggregation of the amyloid-β (Aβ) peptide, linked to Alzheimer's disease, and the aggregation of the Sup35 yeast prion peptide, which resembles the mammalian prion protein linked to spongiform encephalopathies. To facilitate the study of these important peptides, experimentalists have identified small peptide congeners of the full-length proteins that exhibit amyloidogenic behavior, including the KLVFFAE sub-sequence, Aβ16-22, and the GNNQQNY subsequence, Sup357-13. In this study, molecular dynamics simulations were used to examine these peptide fragments encapsulated in reverse micelles (RMs) in order to identify the fundamental principles that govern how sequence and solution environment influence peptide aggregation. Aβ16-22 and Sup357-13 are observed to organize into anti-parallel and parallel β-sheet arrangements. Confinement in the sodium bis(2-ethylhexyl) sulfosuccinate (AOT) reverse micelles is shown to stabilize extended peptide conformations and enhance peptide aggregation. Substantial fluctuations in the reverse micelle shape are observed, in agreement with earlier studies. Shape fluctuations are found to facilitate peptide solvation through interactions between the peptide and AOT surfactant, including direct interaction between non-polar peptide residues and the aliphatic surfactant tails. Computed amide I IR spectra are compared with experimental spectra and found to reflect changes in the peptide structures induced by confinement in the RM environment. Furthermore, examination of the rotational anisotropy decay of water in the RM demonstrates that the water dynamics are sensitive to the presence of peptide as well as the peptide sequence. Overall, our results demonstrate that the RM is a complex confining environment where substantial direct interaction between the surfactant and peptides plays an important role in determining the resulting ensemble of peptide conformations. By extension the results suggest that similarly complex sequence-dependent interactions may determine conformational ensembles of amyloid-forming peptides in a cellular environment.