A multiscale approach to characterize the early aggregation steps of the amyloid-forming peptide GNNQQNY from the yeast prion sup-35

Energy gap of conformational transition related with temperature for the NACore of α-synuclein

Pathological aggregation of α-synuclein (α-syn) into amyloid fibrils is a major feature of Parkinson's disease (PD). The self-assembly of α-syn is mainly governed by a non-amyloid-β component core (NACore). However, the effects of concentrations and temperatures on their conformational transition remain unclear. To answer this question, we investigated the aggregation kinetics of NACore oligomers in silico by performing several independent all-atom molecular dynamics simulations. The simulation results show that tetramers are more prone to form β-sheets at 300 K than dimers and octamers. We also found that the NACore oligomers had higher β-sheet and β-barrel contents at 310 K. The inter-chain hydrophobic interactions, the backbone hydrogen bonding, the residue-residue interactions between V70-V77 as well as V77-V77 play important roles in the aggregation tendency of NACore octamers at 310 K. Interestingly, the energy gap analysis revealed that the conformational transition of NACore oligomers from intermediate states (β-barrel conformation) to stable structures (β-sheet layers) was dependent on the temperatures. In short, our study provides insight into the kinetic and thermodynamic mechanisms of the conformational transition of NACore at different concentrations and temperatures, contributing to a better understanding of the aggregation process of α-syn in Parkinson's disease.

Structures and Dynamics of β-Rich Oligomers of ATTR (105-115) Assembly

Transthyretin (TTR) is a tetrameric homologous protein that can dissociate into monomers. Misfolding and aggregation of TTR can lead to amyloid transthyretin amyloidosis (ATTR), which can cause many diseases (e.g., senile systemic amyloidosis, familial amyloid cardiomyopathy, and familial amyloid polyneuropathy). Despite growing evidence indicating that small oligomers play a critical role in regulating cytotoxicity, the structures of these oligomeric intermediates and their conformational transformations are still unclear, impeding our understanding of neurodegenerative mechanisms and the development of therapeutics targeting early aggregation species. The TTR monomer protein consists of various fragments prone to self-aggregation, including the residue 105-115 sequence. Therefore, our study investigated the assembly progress of ATTR (105-115) peptides using all-atom molecular dynamics simulations. The findings indicate that the probability of β-sheet content increases with increasing numbers of peptides. Additionally, interactions between hydrophobic residues L110 and L111 are crucial for the formation of a β-rich oligomer formation. These β-rich oligomers may adopt β-barrel conformations, potentially toxic oligomer species. Free-energy analysis reveals that β-barrel conformations serve as intermediates for these β-rich oligomers. Our insights into the structural ensemble dynamics of ATTR (105-115) contribute to understanding the physical mechanisms underlying the β-barrel oligomers of ATTR. These findings may shed light on the pathological role of ATTR in neurodegenerative diseases and offer potential therapeutic targets.

Molecular Insights into the Dynamics of Amyloid Fibril Growth: Elongation and Lateral Assembly of GNNQQNY Protofibrils

The self-assembly of peptides and proteins into β-sheet rich amyloid fibrils is linked to both functional and pathological states. In this study, the growth of fibrillar structures of the short peptide GNNQQNY, a fragment from the yeast prion Sup35 protein, was examined. Molecular dynamics simulations were used to study alternative mechanisms of fibril growth, including elongation through binding of monomers as well as fibril self-assembly into larger, more mature structures. It was found that after binding, monomers diffused along preformed fibrils toward the ends, supporting the mechanism of fibril growth via elongation. Lateral assembly of protofibrils was found to occur readily, suggesting that this could be the key to transitioning from isolated fibrils to mature multilayer structures. Overall, the work provides mechanistic insights into the competitive pathways that govern amyloid fibril growth.

Mapping the configurational landscape and aggregation phase behavior of the tau protein fragment PHF6

The PHF6 (Val-Gln-Ile-Val-Tyr-Lys) motif, found in all isoforms of the microtubule-associated protein tau, forms an integral part of ordered cores of amyloid fibrils formed in tauopathies and is thought to play a fundamental role in tau aggregation. Because PHF6 as an isolated hexapeptide assembles into ordered fibrils on its own, it is investigated as a minimal model for insight into the initial stages of aggregation of larger tau fragments. Even for this small peptide, however, the large length and time scales associated with fibrillization pose challenges for simulation studies of its dynamic assembly, equilibrium configurational landscape, and phase behavior. Here, we develop an accurate, bottom-up coarse-grained model of PHF6 for large-scale simulations of its aggregation, which we use to uncover molecular interactions and thermodynamic driving forces governing its assembly. The model, not trained on any explicit information about fibrillar structure, predicts coexistence of formed fibrils with monomers in solution, and we calculate a putative equilibrium phase diagram in concentration-temperature space. We also characterize the configurational and free energetic landscape of PHF6 oligomers. Importantly, we demonstrate with a model of heparin that this widely studied cofactor enhances the aggregation propensity of PHF6 by ordering monomers during nucleation and remaining associated with growing fibrils, consistent with experimentally characterized heparin-tau interactions. Overall, this effort provides detailed molecular insight into PHF6 aggregation thermodynamics and pathways and, furthermore, demonstrates the potential of modern multiscale modeling techniques to produce predictive models of amyloidogenic peptides simultaneously capturing sequence-specific effects and emergent aggregate structures.

Multistep molecular mechanisms of Aβ16-22 fibril formation revealed by lattice Monte Carlo simulations

As a model of self-assembly from disordered monomers to fibrils, the amyloid-β fragment Aβ16-22 was subject to past numerous experimental and computational studies. Because dynamics information between milliseconds and seconds cannot be assessed by both studies, we lack a full understanding of its oligomerization. Lattice simulations are particularly well suited to capture pathways to fibrils. In this study, we explored the aggregation of 10 Aβ16-22 peptides using 65 lattice Monte Carlo simulations, each simulation consisting of 3 × 109 steps. Based on a total of 24 and 41 simulations that converge and do not converge to the fibril state, respectively, we are able to reveal the diversity of the pathways leading to fibril structure and the conformational traps slowing down the fibril formation.

Evolution of large Aβ16-22 aggregates at atomic details and potential of mean force associated to peptide unbinding and fragmentation events

Atomic characterization of large nonfibrillar aggregates of amyloid polypeptides cannot be determined by experimental means. Starting from β-rich aggregates of Y and elongated topologies predicted by coarse-grained simulations and consisting of more than 100 Aβ16-22 peptides, we performed atomistic molecular dynamics (MD), replica exchange with solute scaling (REST2), and umbrella sampling simulations using the CHARMM36m force field in explicit solvent. Here, we explored the dynamics within 3 μs, the free energy landscape, and the potential of mean force associated with either the unbinding of one single peptide in different configurations within the aggregate or fragmentation events of a large number of peptides. Within the time scale of MD and REST2, we find that the aggregates experience slow global conformational plasticity, and remain essentially random coil though we observe slow beta-strand structuring with a dominance of antiparallel beta-sheets over parallel beta-sheets. Enhanced REST2 simulation is able to capture fragmentation events, and the free energy of fragmentation of a large block of peptides is found to be similar to the free energy associated with fibril depolymerization by one chain for longer Aβ sequences.

Dissecting how ALS-associated D290V mutation enhances pathogenic aggregation of hnRNPA2286-291 peptides: Dynamics and conformational ensembles

The aggregation of RNA binding proteins, including hnRNPA1/2, TDP-43 and FUS, is heavily implicated in causing or increasing disease risk for a series of neurodegenerative diseases such as amyotrophic lateral sclerosis (ALS). A recent experimental study demonstrated that an ALS-related D290V mutation in the low complexity domain (LCD) of hnRNPA2 can enhance the aggregation propensity of wild type (WT) hnRNPA2_286-291 peptide. However, the underlying molecular mechanisms remain elusive. Herein, we investigated effects of D290V mutation on aggregation dynamics of hnRNPA2_286-291 peptide and the conformational ensemble of hnRNPA2_286-291 oligomers by performing all-atom molecular dynamic and replica-exchange molecular dynamic simulations. Our simulations demonstrate that D290V mutation greatly reduces the dynamics of hnRNPA2_286-291 peptide and that D290V oligomers possess higher compactness and β-sheet content than WT, indicative of mutation-enhanced aggregation capability. Specifically, D290V mutation strengthens inter-peptide hydrophobic, main-chain hydrogen bonding and side-chain aromatic stacking interactions. Those interactions collectively lead to the enhancement of aggregation capability of hnRNPA2_286-291 peptides. Overall, our study provides insights into the dynamics and thermodynamic mechanisms underlying D290V-induced disease-causing aggregation of hnRNPA2_286-291, which could contribute to better understanding of the transitions from reversible condensates to irreversible pathogenic aggregates of hnRNPA2 LCD in ALS-related diseases.

Disassembly of Amyloid Fibril with Infrared Free Electron Laser

Amyloid fibril causes serious amyloidosis such as neurodegenerative diseases. The structure is composed of rigid β-sheet stacking conformation which makes it hard to disassemble the fibril state without denaturants. Infrared free electron laser (IR-FEL) is an intense picosecond pulsed laser that is oscillated through a linear accelerator, and the oscillation wavelengths are tunable from 3 μm to 100 μm. Many biological and organic compounds can be structurally altered by the mode-selective vibrational excitations due to the wavelength variability and the high-power oscillation energy (10-50 mJ/cm²). We have found that several different kinds of amyloid fibrils in amino acid sequences were commonly disassembled by the irradiation tuned to amide I (6.1-6.2 μm) where the abundance of β-sheet decreased while that of α-helix increased by the vibrational excitation of amide bonds. In this review, we would like to introduce the IR-FEL oscillation system briefly and describe combination studies of experiments and molecular dynamics simulations on disassembling amyloid fibrils of a short peptide (GNNQQNY) from yeast prion and 11-residue peptide (NFLNCYVSGFH) from β2-microglobulin as representative models. Finally, possible applications of IR-FEL for amyloid research can be proposed as a future outlook.

Computer Simulations Aimed at Exploring Protein Aggregation and Dissociation

Protein aggregation can lead to well-defined structures that are functional, but is also the cause of the death of neuron cells in many neurodegenerative diseases. The complexity of the molecular events involved in the aggregation kinetics of amyloid proteins and the transient and heterogeneous characters of all oligomers prevent high-resolution structural experiments. As a result, computer simulations have been used to determine the atomic structures of amyloid proteins at different association stages as well as to understand fibril dissociation. In this chapter, we first review the current computer simulation methods used for aggregation with some atomistic and coarse-grained results aimed at better characterizing the early formed oligomers and amyloid fibril formation. Then we present the applications of non-equilibrium molecular dynamics simulations to comprehend the dissociation of protein assemblies.

Dynamics of Amyloid Formation from Simplified Representation to Atomistic Simulations

Amyloid fibril formation is an intrinsic property of short peptides, non-disease proteins, and proteins associated with neurodegenerative diseases. Aggregates of the Aβ and tau proteins, the α-synuclein protein, and the prion protein are observed in the brain of Alzheimer's, Parkinson's, and prion disease patients, respectively. Due to the transient short-range and long-range interactions of all species and their high aggregation propensities, the conformational ensemble of these devastating proteins, the exception being for the monomeric prion protein, remains elusive by standard structural biology methods in bulk solution and in lipid membranes. To overcome these limitations, an increasing number of simulations using different sampling methods and protein models have been performed. In this chapter, we first review our main contributions to the field of amyloid protein simulations aimed at understanding the early aggregation steps of short linear amyloid peptides, the conformational ensemble of the Aβ40/42 dimers in bulk solution, and the stability of Aβ aggregates in lipid membrane models. Then we focus on our studies on the interactions of amyloid peptides/inhibitors to prevent aggregation, and long amyloid sequences, including new results on a monomeric tau construct.

Nucleation-dependent Aggregation Kinetics of Yeast Sup35 Fragment GNNQQNY

An N-terminal hepta-peptide sequence of yeast prion protein Sup35 with the sequence GNNQQNY is widely used as a model system for amyloid fibril formation. In this study, we used a reproducible solubilisation protocol that allows the generation of a homogenous monomeric solution of GNNQQNY to uncover the molecular details of its self-assembly mechanism. The aggregation kinetics data show that the GNNQQNY sequence follows nucleation-dependent aggregation kinetics with a critical nucleus of size ~7 monomers and that the efficiency of nucleation were found to be inversely related to the reaction temperature. The nucleus reduces the thermodynamic energy barrier by acting as a template for further self-assembly and results in highly ordered amyloid fibrils. The fibers grown at different temperatures showed similar Thioflavin T fluorescence, Congo-red binding and β-sheet rich structures displaying a characteristic cross-β diffraction pattern. These aggregates also share morphological and structural identity with those reported earlier. The mature GNNQQNY fibers did not exert significant oxidative stress or cytotoxicity upon incubating with differentiated SHSY5Y cells. To our knowledge, this is the first study to experimentally validate previous nucleus size predictions based on theoretical and molecular dynamics simulations. These findings provide the basis for understanding the kinetics and thermodynamics of amyloid nucleation and elongation of amyloidogenic proteins/peptides associated with many systemic and neurodegenerative diseases.

Peptide-based vesicles and droplets: a review

Peptide assembly is an increasingly important field of study due to the versatility, tunability and vast design space of amino acid based biomolecular assemblies. Peptides can be precisely engineered to possess various useful properties such as the ability to form supramolecular assemblies, desired response to pH, or thermal stability. These peptide supramolecular assemblies have diverse morphologies including vesicles, nanotubes, nanorods and ribbons. Of specific interest is the domain of engineering peptides that aggregate into spherical nanostructures due to their encapsulation properties: the ability to hold, transport and release chemical payloads in a controllable manner. This is invaluable to the fields of nanomedicine and targeted drug delivery. In this review, the state of the art in the domain of peptide-based vesicles and nanospheres is summarized. Specifically, an overview of the assembly of peptides into nanovesicles and nanospheres is provided. Both aromatic as well as aliphatic side chain amino acids are discussed. The domain of aromatic side chained amino acid residues is largely dominated by phenylalanine based peptides and variants thereof. Tyrosine also demonstrates similar aggregation properties. Both experimentally and computationally driven approaches are discussed. The domain of aliphatic amino acid residues based vesicles and droplets is broader, and details multiple amino acid residues such as alanine, valine, lysine, glycine, proline, and aspartic acid. Finally, a discussion on potential future directions is provided.

Sampling of the conformational landscape of small proteins with Monte Carlo methods

Computer simulation provides an increasingly realistic picture of large-scale conformational change of proteins, but investigations remain fundamentally constrained by the femtosecond timestep of molecular dynamics simulations. For this reason, many biologically interesting questions cannot be addressed using accessible state-of-the-art computational resources. Here, we report the development of an all-atom Monte Carlo approach that permits the modelling of the large-scale conformational change of proteins using standard off-the-shelf computational hardware and standard all-atom force fields. We demonstrate extensive thermodynamic characterization of the folding process of the α-helical Trp-cage, the Villin headpiece and the β-sheet WW-domain. We fully characterize the free energy landscape, transition states, energy barriers between different states, and the per-residue stability of individual amino acids over a wide temperature range. We demonstrate that a state-of-the-art intramolecular force field can be combined with an implicit solvent model to obtain a high quality of the folded structures and also discuss limitations that still remain.

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.

Infrared Laser-Induced Amyloid Fibril Dissociation: A Joint Experimental/Theoretical Study on the GNNQQNY Peptide

Neurodegenerative diseases are usually characterized by plaques made of well-ordered aggregates of distinct amyloid proteins. Dissociating these very stable amyloid plaques is a critical clinical issue. In this study, we present a joint mid-infrared free electron laser experiment/nonequilibrium molecular dynamics simulation to understand the dissociation process of a representative example GNNQQNY fibril. By tuning the laser frequency to the amide I band of the fibril, the resonance takes place and dissociation is occurred. With the calculated and observed wide-angle X-ray scattering profiles and secondary structures before and after laser irradiation being identical, we can propose a dissociation mechanism with high confidence from our simulations. We find that dissociation starts in the core of the fibrils by fragmenting the intermolecular hydrogen bonds and separating the peptides and then propagates to the fibril extremities leading to the formation of unstructured expanded oligomers. We suggest that this should be a generic mechanism of the laser-induced dissociation of amyloid fibrils.

Structures of the intrinsically disordered Aβ, tau and α-synuclein proteins in aqueous solution from computer simulations

Intrinsically disordered proteins (IDPs) play many biological roles in the human proteome ranging from vesicular transport, signal transduction to neurodegenerative diseases. The Aβ and tau proteins, and the α-synuclein protein, key players in Alzheimer's and Parkinson's diseases, respectively are fully disordered at the monomer level. The structural heterogeneity of the monomeric and oligomeric states and the high self-assembly propensity of these three IDPs have precluded experimental structural determination. Simulations have been used to determine the atomic structures of these IDPs. In this article, we review recent computer models to capture the equilibrium ensemble of Aβ, tau and α-synuclein proteins at different association steps in aqueous solution and present new results of the PEP-FOLD framework on α-synuclein monomer.

Theoretical and computational advances in protein misfolding

Misfolded proteins escape the cellular quality control mechanism and fail to fold properly or remain correctly folded leading to a loss in their functional specificity. Thus misfolding of proteins cause a large number of very different diseases ranging from errors in metabolism to various types of complex neurodegenerative diseases. A theoretical and computational perspective of protein misfolding is presented with a special emphasis on its salient features, mechanism and consequences. These insights quantitatively analyze different determinants of misfolding, that may be applied to design disease specific molecular targets.

Aggregation of disease-related peptides

Protein misfolding and aggregation of amyloid proteins is the fundamental cause of more than 20 diseases. Molecular mechanisms of the self-assembly and the formation of the toxic aggregates are still elusive. Computer simulations have been intensively used to study the aggregation of amyloid peptides of various amino acid lengths related to neurodegenerative diseases. We review atomistic and coarse-grained simulations of short amyloid peptides aimed at determining their transient oligomeric structures and the early and late aggregation steps.

A computational approach to get insights into multiple faces of additives in modulation of protein aggregation pathways

An enormous population worldwide is presently confronted with debilitating neurodegenerative diseases. The etiology of the disease is connected to protein aggregation and the events involved therein. Thus, a complete understanding of an inhibitor at different stages in the process is imperative for the formulation of a drug molecule. This review presents a detailed summary of the current status of different cosolvents. It further develops how the complex aggregation pathway can be simplified into three steps common to all proteins and the way computer simulations can be exploited to gain insights into the ways by which known inhibitors can affect all these stages. Computation of theoretical parameters in this regard and their correlation with experimental techniques is accentuated. In addition to providing an outline of the scope of different additives, this review showcases the way by which the problem of analyzing an effect of an additive can be addressed effectively via MD simulations.

Molecular Model for the Self-Assembly of the Cyclic Lipodepsipeptide Pseudodesmin A

Self-assembly of peptides into supramolecular structures represents an active field of research with potential applications ranging from material science to medicine. Their study typically involves the application of a large toolbox of spectroscopic and imaging techniques. However, quite often, the structural aspects remain underexposed. Besides, molecular modeling of the self-assembly process is usually difficult to handle, since a vast conformational space has to be sampled. Here, we have used an approach that combines short molecular dynamics simulations for peptide dimerization and NMR restraints to build a model of the supramolecular structure from the dimeric units. Experimental NMR data notably provide crucial information about the conformation of the monomeric units, the supramolecular assembly dimensions, and the orientation of the individual peptides within the assembly. This in silico/in vitro mixed approach enables us to define accurate atomistic models of supramolecular structures of the bacterial cyclic lipodepsipeptide pseudodesmin A.

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.

Amyloid-β(29-42) Dimeric Conformations in Membranes Rich in Omega-3 and Omega-6 Polyunsaturated Fatty Acids

The omega-3 and omega-6 polyunsaturated fatty acids are two important components of cell membranes in human brains. When incorporated into phospholipids, omega-3 slows the progression of Alzheimer's disease (AD), whereas omega-6 is linked to increased risk of AD. Little is known on the amyloid-β (Aβ) conformations in membranes rich in omega-3 and omega-6 phospholipids. Herein, the structural properties of the Aβ_29-42 dimer embedded in both fatty acid membranes were comparatively studied to a 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphocholine (POPC) bilayer using all-atom molecular dynamics (MD) simulations. Starting from α-helix, both omega-6 and omega-3 membranes promote new orientations and conformations of the dimer, in agreement with the observed dependence of Aβ production upon addition of these two fatty acids. This conformational result is corroborated by atomistic MD simulations of the dimer of the 99 amino acid C-terminal fragment of amyloid precursor protein spanning the residues 15-55. Starting from β-sheet, omega-6 membrane promotes helical and disordered structures of Aβ_29-42 dimer, whereas omega-3 membrane preserves the β-sheet structures differing however from those observed in POPC. Remarkably, the mixture of the two fatty acids and POPC depicts another conformational ensemble of the Aβ_29-42 dimer. This finding demonstrates that variation in the abundance of the molecular phospholipids, which changes with age, modulates membrane-embedded Aβ oligomerization.

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.

Multi-scale simulations of biological systems using the OPEP coarse-grained model

Biomolecules are complex machines that are optimized by evolution to properly fulfill or contribute to a variety of biochemical tasks in the cellular environment. Computer simulations based on quantum mechanics and atomistic force fields have been proven to be a powerful microscope for obtaining valuable insights into many biological, physical, and chemical processes. Many interesting phenomena involve, however, a time scale and a number of degrees of freedom, notably if crowding is considered, that cannot be explored at an atomistic resolution. To bridge the gap between reality and simulation, many different advanced computational techniques and coarse-grained (CG) models have been developed. Here, we report some applications of the CG OPEP protein model to amyloid fibril formation, the response of catch-bond proteins to two types of fluid flow, and interactive simulations to fold peptides with well-defined 3D structures or with intrinsic disorder.

Mechanism of Initiation, Association, and Formation of Amyloid Fibrils Modeled with the N-Terminal Peptide Fragment, IKYLEFIS, of Myoglobin G-Helix

Some peptides and proteins undergo self-aggregation under certain conditions, leading to amyloid fibrils formation, which is related to many disease conditions. It is important to understand such amyloid fibrils formation to provide mechanistic detail that governs the process. A predominantly α-helical myoglobin has been reported recently to readily form amyloid fibrils at a higher temperature, similar to its G-helix segment. Here, we have investigated the mechanism of amyloid fibrils formation by performing multiple long molecular dynamics simulations (27 μs) on the N-terminal segment of the G-helix of myoglobin. These simulations resulted in the formation of a single-layered tetrameric β-sheet with mixed parallel and antiparallel β-strands and this is the most common event irrespective of many different starting structures. Formation of the single-layered tetrameric β-sheet takes place following three distinctive pathways. The process of fibril initiation is dependent on temperature. Further, this study provides mechanistic insights into the formation of multilayered fibrilar structure, which could be applicable to a wider variety of peptides or proteins to understand the amyloidogenesis.

Lattice model for amyloid peptides: OPEP force field parametrization and applications to the nucleus size of Alzheimer's peptides

Coarse-grained protein lattice models approximate atomistic details and keep the essential interactions. They are, therefore, suitable for capturing generic features of protein folding and amyloid formation at low computational cost. As our aim is to study the critical nucleus sizes of two experimentally well-characterized peptide fragments Aβ16-22 and Aβ37-42 of the full length Aβ1-42 Alzheimer's peptide, it is important that simulations with the lattice model reproduce all-atom simulations. In this study, we present a comprehensive force field parameterization based on the OPEP (Optimized Potential for Efficient protein structure Prediction) force field for an on-lattice protein model, which incorporates explicitly the formation of hydrogen bonds and directions of side-chains. Our bottom-up approach starts with the determination of the best lattice force parameters for the Aβ16-22 dimer by fitting its equilibrium parallel and anti-parallel β-sheet populations to all-atom simulation results. Surprisingly, the calibrated force field is transferable to the trimer of Aβ16-22 and the dimer and trimer of Aβ37-42. Encouraged by this finding, we characterized the free energy landscapes of the two decamers. The dominant structure of the Aβ16-22 decamer matches the microcrystal structure. Pushing the simulations for aggregates between 4-mer and 12-mer suggests a nucleus size for fibril formation of 10 chains. In contrast, the Aβ37-42 decamer is largely disordered with mixed by parallel and antiparallel chains, suggesting that the nucleus size is >10 peptides. Our refined force field coupled to this on-lattice model should provide useful insights into the critical nucleation number associated with neurodegenerative diseases.

Hydrodynamic effects on β-amyloid (16-22) peptide aggregation

Computer simulations based on simplified representations are routinely used to explore the early steps of amyloid aggregation. However, when protein models with implicit solvent are employed, these simulations miss the effect of solvent induced correlations on the aggregation kinetics and lifetimes of metastable states. In this work, we apply the multi-scale Lattice Boltzmann Molecular Dynamics technique (LBMD) to investigate the initial aggregation phases of the amyloid Aβ16-22 peptide. LBMD includes naturally hydrodynamic interactions (HIs) via a kinetic on-lattice representation of the fluid kinetics. The peptides are represented by the flexible OPEP coarse-grained force field. First, we have tuned the essential parameters that control the coupling between the molecular and fluid evolutions in order to reproduce the experimental diffusivity of elementary species. The method is then deployed to investigate the effect of HIs on the aggregation of 100 and 1000 Aβ16-22 peptides. We show that HIs clearly impact the aggregation process and the fluctuations of the oligomer sizes by favouring the fusion and exchange dynamics of oligomers between aggregates. HIs also guide the growth of the leading largest cluster. For the 100 Aβ16-22 peptide system, the simulation of ∼300 ns allowed us to observe the transition from ellipsoidal assemblies to an elongated and slightly twisted aggregate involving almost the totality of the peptides. For the 1000 Aβ16-22 peptides, a system of unprecedented size at quasi-atomistic resolution, we were able to explore a branched disordered fibril-like structure that has never been described by other computer simulations, but has been observed experimentally.

Inhibitory effect of hydrophobic fullerenes on the β-sheet-rich oligomers of a hydrophilic GNNQQNY peptide revealed by atomistic simulations

Fullerenes suppress fibril-like β-sheet oligomers by interacting strongly with the nonpolar aliphatic groups of polar residues of GNNQQNY peptide, thus inhibit peptide aggregation.

Assembly of Triblock Amphiphilic Peptides into One-Dimensional Aggregates and Network Formation

Peptide assembly plays a key role in both neurological diseases and development of novel biomaterials with well-defined nanostructures. Synthetic model peptides provide a unique platform to explore the role of intermolecular interactions in the assembly process. A triblock peptide architecture designed by the Hartgerink group is a versatile system which relies on Coulomb interactions, hydrogen bonding, and hydrophobicity to guide these peptides' assembly at three different length scales: β-sheets, double-wall ribbon-like aggregates, and finally a highly porous network structure which can support gels with ≤1% by weight peptide concentration. In this study, by using molecular dynamics simulations of a structure based implicit solvent coarse grained model, we analyzed this hierarchical assembly process. Parametrization of our CG model is based on multiple-state points from atomistic simulations, which enables this model to represent the conformational adaptability of the triblock peptide molecule based on the surrounding medium. Our results indicate that emergence of the double-wall β-sheet packing mechanism, proposed in light of the experimental evidence, strongly depends on the subtle balance of the intermolecular forces. We demonstrate that, even though backbone hydrogen bonding dominates the early nucleation stages, depending on the strength of the hydrophobic and Coulomb forces, alternative structures such as zero-dimensional aggregates with two β-sheets oriented orthogonally (which we refer to as a cross-packed structure) and β-sheets with misoriented hydrophobic side chains are also feasible. We discuss the implications of these competing structures for the three different length scales of assembly by systematically investigating the influence of density, counterion valency, and hydrophobicity.

An Atomistic View of Amyloidogenic Self-assembly: Structure and Dynamics of Heterogeneous Conformational States in the Pre-nucleation Phase

The formation of well-defined filamentous amyloid structures involves a polydisperse collection of oligomeric states for which relatively little is known in terms of structural organization. Here we use extensive, unbiased explicit solvent molecular dynamics (MD) simulations to investigate the structural and dynamical features of oligomeric aggregates formed by a number of highly amyloidogenic peptides at atomistic resolution on the μs time scale. A consensus approach has been adopted to analyse the simulations in multiple force fields, yielding an in-depth characterization of pre-fibrillar oligomers and their global and local structure properties. A collision cross section analysis revealed structurally heterogeneous aggregate ensembles for the individual oligomeric states that lack a single defined quaternary structure during the pre-nucleation phase. To gain insight into the conformational space sampled in early aggregates, we probed their substructure and found emerging β-sheet subunit layers and a multitude of ordered intermolecular β-structure motifs with growing aggregate size. Among those, anti-parallel out-of-register β-strands compatible with toxic β-barrel oligomers were particularly prevalent already in smaller aggregates and formed prior to ordered fibrillar structure elements. Notably, also distinct fibril-like conformations emerged in the oligomeric state and underscore the notion that pre-nucleated oligomers serve as a critical intermediate step on-pathway to fibrils.

Understanding Amyloid-β Oligomerization at the Molecular Level: The Role of the Fibril Surface

The aggregation of the amyloid β-peptide into fibrils is a complex process that involves mechanisms such as primary and secondary nucleation, fibril elongation and fibril fragmentation. Some of these processes generate neurotoxic Aβ oligomers, which are involved in the development of Alzheimer's disease. Recent experimental studies have emphasized the role of the fibril as a catalytic surface for the production of highly toxic oligomers during secondary nucleation. By using molecular dynamics simulations, we show that it is the hydrophobic fibril region that causes the structural changes required for catalyzing the formation of β-sheet-rich Aβ1-42 oligomers on the fibril surface. These results reveal, for the first time, the molecular basis of the secondary nucleation pathway.

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.

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.

Protein Simulations in Fluids: Coupling the OPEP Coarse-Grained Force Field with Hydrodynamics

A novel simulation framework that integrates the OPEP coarse-grained (CG) model for proteins with the Lattice Boltzmann (LB) methodology to account for the fluid solvent at mesoscale level is presented. OPEP is a very efficient, water-free and electrostatic-free force field that reproduces at quasi-atomistic detail processes like peptide folding, structural rearrangements, and aggregation dynamics. The LB method is based on the kinetic description of the solvent in order to solve the fluid mechanics under a wide range of conditions, with the further advantage of being highly scalable on parallel architectures. The capabilities of the approach are presented, and it is shown that the strategy is effective in exploring the role of hydrodynamics on protein relaxation and peptide aggregation. The end result is a strategy for modeling systems of thousands of proteins, such as in the case of dense protein suspensions. The future perspectives of the multiscale approach are also discussed.

The architecture of amyloid-like peptide fibrils revealed by X-ray scattering, diffraction and electron microscopy

Structural analysis of protein fibrillation is inherently challenging. Given the crucial role of fibrils in amyloid diseases, method advancement is urgently needed. A hybrid modelling approach is presented enabling detailed analysis of a highly ordered and hierarchically organized fibril of the GNNQQNY peptide fragment of a yeast prion protein. Data from small-angle X-ray solution scattering, fibre diffraction and electron microscopy are combined with existing high-resolution X-ray crystallographic structures to investigate the fibrillation process and the hierarchical fibril structure of the peptide fragment. The elongation of these fibrils proceeds without the accumulation of any detectable amount of intermediate oligomeric species, as is otherwise reported for, for example, glucagon, insulin and α-synuclein. Ribbons constituted of linearly arranged protofilaments are formed. An additional hierarchical layer is generated via the pairing of ribbons during fibril maturation. Based on the complementary data, a quasi-atomic resolution model of the protofilament peptide arrangement is suggested. The peptide structure appears in a β-sheet arrangement reminiscent of the β-zipper structures evident from high-resolution crystal structures, with specific differences in the relative peptide orientation. The complexity of protein fibrillation and structure emphasizes the need to use multiple complementary methods.

Role of monomer arrangement in the amyloid self-assembly

Assembly of amyloid proteins into aggregates requires the ordering of the monomers in oligomers and especially in such highly organized structures as fibrils. This ordering is accompanied by structural transitions leading to the formation of ordered β-structural motifs in proteins and peptides lacking secondary structures. To characterize the effect of the monomer arrangements on the aggregation process at various stages, we performed comparative studies of the yeast prion protein Sup35 heptapeptide (GNNQQNY) along with its dimeric form CGNNQQNY-(d-Pro)-G-GNNQQNY. The (d-Pro)-G linker in this construct is capable of adopting a β-turn, facilitating the assembly of the dimer into the dimeric antiparallel hairpin structure (AP-hairpin). We applied Atomic Force Microscopy (AFM) techniques to follow peptide-peptide interactions at the single molecule level, to visualize the morphology of aggregates formed by both constructs, thioflavin T (ThT) fluorescence to follow the aggregation kinetics, and circular dichroism (CD) spectroscopy to characterize the secondary structure of the constructs. The ThT fluorescence data showed that the AP-hairpin aggregation kinetics is insensitive to the external environment such as ionic strength and pH contrary to the monomers the kinetics of which depends dramatically on the ionic strength and pH. The AFM topographic imaging revealed that AP-hairpins primarily assemble into globular aggregates, whereas linear fibrils are primary assemblies of the monomers suggesting that both constructs follow different aggregation pathways during the self-assembly. These morphological differences are in line with the AFM force spectroscopy experiments and CD spectroscopy measurements, suggesting that the AP-hairpin is structurally rigid regardless of changes of environmental factors.

C-terminal domain swapping of SSB changes the size of the ssDNA binding site

Single-stranded DNA-binding protein (SSB) plays an important role in DNA metabolism, including DNA replication, repair, and recombination, and is therefore essential for cell survival. Bacterial SSB consists of an N-terminal ssDNA-binding/oligomerization domain and a flexible C-terminal protein-protein interaction domain. We characterized the ssDNA-binding properties of Klebsiella pneumoniae SSB (KpSSB), Salmonella enterica Serovar Typhimurium LT2 SSB (StSSB), Pseudomonas aeruginosa PAO1 SSB (PaSSB), and two chimeric KpSSB proteins, namely, KpSSBnStSSBc and KpSSBnPaSSBc. The C-terminal domain of StSSB or PaSSB was exchanged with that of KpSSB through protein chimeragenesis. By using the electrophoretic mobility shift assay, we characterized the stoichiometry of KpSSB, StSSB, PaSSB, KpSSBnStSSBc, and KpSSBnPaSSBc, complexed with a series of ssDNA homopolymers. The binding site sizes were determined to be 26 ± 2, 21 ± 2, 29 ± 2, 21 ± 2, and 29 ± 2 nucleotides (nt), respectively. Comparison of the binding site sizes of KpSSB, KpSSBnStSSBc, and KpSSBnPaSSBc showed that the C-terminal domain swapping of SSB changes the size of the binding site. Our observations suggest that not only the conserved N-terminal domain but also the C-terminal domain of SSB is an important determinant for ssDNA binding.

The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems

The OPEP coarse-grained protein model has been applied to a wide range of applications since its first release 15 years ago. The model, which combines energetic and structural accuracy and chemical specificity, allows the study of single protein properties, DNA-RNA complexes, amyloid fibril formation and protein suspensions in a crowded environment. Here we first review the current state of the model and the most exciting applications using advanced conformational sampling methods. We then present the current limitations and a perspective on the ongoing developments.