Low molecular weight oligomers of amyloid peptides display β-barrel conformations: A replica exchange molecular dynamics study in explicit solvent

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.

Self-Assembly of Amyloid-Beta (Aβ) Peptides from Solution to Near In Vivo Conditions

Understanding the atomistic resolution changes during the self-assembly of amyloid peptides or proteins is important to develop compounds or conditions to alter the aggregation pathways and suppress the toxicity and potentially aid in the development of drugs. However, the complexity of protein aggregation and the transient order/disorder of oligomers along the pathways to fibril are very challenging. In this Perspective, we discuss computational studies of amyloid polypeptides carried out under various conditions, including conditions closely mimicking in vivo and point out the challenges in obtaining physiologically relevant results, focusing mainly on the amyloid-beta Aβ peptides.

Computational Models for the Study of Protein Aggregation

Protein aggregation has been studied by many groups around the world for many years because it can be the cause of a number of neurodegenerative diseases that have no effective treatment. Obtaining the structure of related fibrils and toxic oligomers, as well as describing the pathways and main factors that govern the self-organization process, is of paramount importance, but it is also very difficult. To solve this problem, experimental and computational methods are often combined to get the most out of each method. The effectiveness of the computational approach largely depends on the construction of a reasonable molecular model. Here we discussed different versions of the four most popular all-atom force fields AMBER, CHARMM, GROMOS, and OPLS, which have been developed for folded and intrinsically disordered proteins, or both. Continuous and discrete coarse-grained models, which were mainly used to study the kinetics of aggregation, are also summarized.

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.

Spontaneous Formation of β-sheet Nano-barrels during the Early Aggregation of Alzheimer's Amyloid Beta

Soluble low-molecular-weight oligomers formed during the early aggregation of amyloid peptides have been hypothesized as a major toxic species of amyloidogenesis. Herein, we performed the first synergic in silico, in vitro and in vivo validations of the structure, dynamics and toxicity of Aβ42 oligomers. Aβ peptides readily assembled into β-rich oligomers comprised of extended β-hairpins and β-strands. Nanosized β-barrels were observed with certainty with simulations, transmission electron microscopy and Fourier transform infrared spectroscopy, corroborated by immunohistochemistry, cell viability, apoptosis, inflammation, autophagy and animal behavior assays. Secondary and tertiary structural proprieties of these oligomers, such as the sequence regions with high β-sheet propensities and inter-residue contact frequency patterns, were similar to the properties known for Aβ fibrils. The unambiguous spontaneous formation of β-barrels in the early aggregation of Aβ42 supports their roles as the common toxic intermediates in Alzheimer's pathobiology and a target for Alzheimer's therapeutics.

Defining the Neuropathological Aggresome across in Silico, in Vitro, and ex Vivo Experiments

The loss of proteostasis over the life course is associated with a wide range of debilitating degenerative diseases and is a central hallmark of human aging. When left unchecked, proteins that are intrinsically disordered can pathologically aggregate into highly ordered fibrils, plaques, and tangles (termed amyloids), which are associated with countless disorders such as Alzheimer's disease, Parkinson's disease, type II diabetes, cancer, and even certain viral infections. However, despite significant advances in protein folding and solution biophysics techniques, determining the molecular cause of these conditions in humans has remained elusive. This has been due, in part, to recent discoveries showing that soluble protein oligomers, not insoluble fibrils or plaques, drive the majority of pathological processes. This has subsequently led researchers to focus instead on heterogeneous and often promiscuous protein oligomers. Unfortunately, significant gaps remain in how to prepare, model, experimentally corroborate, and extract amyloid oligomers relevant to human disease in a systematic manner. This Review will report on each of these techniques and their successes and shortcomings in an attempt to standardize comparisons between protein oligomers across disciplines, especially in the context of neurodegeneration. By standardizing multiple techniques and identifying their common overlap, a clearer picture of the soluble neuropathological aggresome can be constructed and used as a baseline for studying human disease and aging.

Direct Observation of β-Barrel Intermediates in the Self-Assembly of Toxic SOD128-38 and Absence in Nontoxic Glycine Mutants

Soluble low-molecular-weight oligomers formed during the early stage of amyloid aggregation are considered the major toxic species in amyloidosis. The structure-function relationship between oligomeric assemblies and the cytotoxicity in amyloid diseases are still elusive due to the heterogeneous and transient nature of these aggregation intermediates. To uncover the structural characteristics of toxic oligomeric intermediates, we compared the self-assembly dynamics and structures of SOD1_28-38, a cytotoxic fragment of the superoxide dismutase 1 (SOD1) associated with the amyotrophic lateral sclerosis, with its two nontoxic mutants G33V and G33W using molecular dynamics simulations. Single-point glycine substitutions in SOD1_28-38 have been reported to abolish the amyloid toxicity. Our simulation results showed that the toxic SOD1_28-38 and its nontoxic mutants followed different aggregation pathways featuring distinct aggregation intermediates. Specifically, wild-type SOD1_28-38 initially self-assembled into random-coil-rich oligomers, among which fibrillar aggregates composed of well-defined curved single-layer β-sheets were nucleated via coil-to-sheet conversions and the formation of β-barrels as intermediates. In contrast, the nontoxic G33V/G33W mutants readily assembled into small β-sheet-rich oligomers and then coagulated with each other into cross-β fibrils formed by two-layer β-sheets without forming β-barrels as the intermediates. The direct observation of β-barrel oligomers during the assembly of toxic SOD1_28-38 fragments but not the nontoxic glycine-substitution mutants strongly supports β-barrels as the toxic oligomers in amyloidosis, probably via interactions with the cell membrane and forming amyloid pores. With well-defined structures, the β-barrel might serve as a novel therapeutic target against amyloid-related diseases.

The Early Phase of β2-Microglobulin Aggregation: Perspectives From Molecular Simulations

Protein β2-microglobulin is the causing agent of two amyloidosis, dialysis related amyloidosis (DRA), affecting the bones and cartilages of individuals with chronic renal failure undergoing long-term hemodialysis, and a systemic amyloidosis, found in one French family, which impairs visceral organs. The protein’s small size and its biomedical significance attracted the attention of theoretical scientists, and there are now several studies addressing its aggregation mechanism in the context of molecular simulations. Here, we review the early phase of β2-microglobulin aggregation, by focusing on the identification and structural characterization of monomers with the ability to trigger aggregation, and initial small oligomers (dimers, tetramers, hexamers etc.) formed in the so-called nucleation phase. We focus our analysis on results from molecular simulations and integrate our views with those coming from in vitro experiments to provide a broader perspective of this interesting field of research. We also outline directions for future computer simulation studies.

Stability of Aβ11-40 Trimers with Parallel and Antiparallel β-Sheet Organizations in a Membrane-Mimicking Environment by Replica Exchange Molecular Dynamics Simulation

The aggregation of the amyloid (Aβ) peptide of 39-43 amino acids into plaques is observed in the brain of Alzheimer's disease (AD) patients, but the mechanisms underlying the neurotoxicity of Aβ oligomers are still elusive. One suggested initial mechanism is related to the implications of amyloid membrane interactions, but characterization of these assemblies is challenging by experimental means. In this study, we have explored the stability of a trimer of Aβ11-40 in parallel and antiparallel β-sheet structures for the wild-type sequence and its F20W mutant in a dipalmitoylphosphatidylcholine membrane using atomistic replica exchange molecular dynamic simulations. We show that both the U-shape organization and the assembly of β-hairpins are maintained in the membrane and are resistant to the mutation F20W. In contrast the models are destabilized by the F19P mutation. Overall, our results indicate that these two assemblies represent minimal seeds or nuclei for the formation of either amyloid fibrils, a variety of β-barrel pores, or various aggregates for many Aβ sequences in a membrane-mimicking environment.

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.

ATP Controls the Aggregation of Aβ16-22 Peptides

The oligomerization of Aβ_16-22 peptide, which is the hydrophobic core region of full-length Aβ_1-42, causes Alzheimer's disease (AD). This progressive neurodegenerative disease affects over 44 million people worldwide. However, very few synthesized drug molecules are available to inhibit the aggregation of Aβ. Recently, experimental studies have shown that the biological ATP molecule prevents Aβ fibrillation at the millimolar scale; however, the significance of ATP molecules on Aβ fibrillation and the mechanism behind it remain elusive. We have carried out a total of 7.5 μs extensive all-atom molecular dynamics and 8.82 μs of umbrella sampling in explicit water using AMBER14SB, AMBER99SB-ILDN, and AMBER-FB15 force fields for Aβ_16-22 peptide, to investigate the role of ATP on the disruption of Aβ_16-22 prefibrils. From various analyses, such as secondary structure analysis, residue-wise contact map, SASA, and interaction energies, we have observed that, in the presence of ATP, the aggregation of Aβ_16-22 peptide is very unfavorable. Moreover, the biological molecule ATP interacts with the Aβ_16-22 peptide via hydrogen bonding, π-π stacking, and NH-π interactions which, ultimately, prevent the aggregation of Aβ_16-22 peptide. Hence, we assume that the deficiency of ATP may cause Alzheimer's disease (AD).

Protein structure predictions by enhanced conformational sampling methods

In this Special Festschrift Issue for the celebration of Professor Nobuhiro Gō's 80th birthday, we review enhanced conformational sampling methods for protein structure predictions. We present several generalized-ensemble algorithms such as multicanonical algorithm, replica-exchange method, etc. and parallel Monte Carlo or molecular dynamics method with genetic crossover. Examples of the results of these methods applied to the predictions of protein tertiary structures are also presented.

Amyloid Self-Assembly of hIAPP8-20 via the Accumulation of Helical Oligomers, α-Helix to β-Sheet Transition, and Formation of β-Barrel Intermediates

The self-assembly of human islet amyloid polypeptide (hIAPP) into β-sheet-rich nanofibrils is associated with the pathogeny of type 2 diabetes. Soluble hIAPP is intrinsically disordered with N-terminal residues 8-17 as α-helices. To understand the contribution of the N-terminal helix to the aggregation of full-length hIAPP, here the oligomerization dynamics of the hIAPP fragment 8-20 (hIAPP8-20) are investigated with combined computational and experimental approaches. hIAPP8-20 forms cross-β nanofibrils in silico from isolated helical monomers via the helical oligomers and α-helices to β-sheets transition, as confirmed by transmission electron microscopy, atomic force microscopy, circular dichroism spectroscopy, Fourier transform infrared spectroscopy, and reversed-phase high performance liquid chromatography. The computational results also suggest that the critical nucleus of aggregation corresponds to hexamers, consistent with a recent mass-spectroscopy study of hIAPP8-20 aggregation. hIAPP8-20 oligomers smaller than hexamers are helical and unstable, while the α-to-β transition starts from the hexamers. Converted β-sheet-rich oligomers first form β-barrel structures as intermediates before aggregating into cross-β nanofibrils. This study uncovers a complete picture of hIAPP8-20 peptide oligomerization, aggregation nucleation via conformational conversion, formation of β-barrel intermediates, and assembly of cross-β protofibrils, thereby shedding light on the aggregation of full-length hIAPP, a hallmark of pancreatic beta-cell degeneration.

Multiscale Aggregation of the Amyloid Aβ16-22 Peptide: From Disordered Coagulation and Lateral Branching to Amorphous Prefibrils

In this work we investigate the multiscale dynamics of the aggregation process of an amyloid peptide, Aβ_16-22. By performing massive coarse-grained simulations at the quasi-atomistic resolution and including hydrodynamic effects, we followed the formation and growth of a large elongated aggregate and its slow structuring. The elongation proceeds via a two-step nucleation mechanism with disordered aggregates formed initially and subsequently fusing to elongate the amorphous prefibril. A variety of coagulation events coexist, including lateral growth. The latter mechanism, sustained by long-range hydrodynamics correlations, actually can create a large branched structure spanning a few tens of nanometers. Our findings confirm the experimental hypothesis of a critical contribution of lateral growth to the amyloid aggregation kinetics and the capability of our model to sample critical structures like prefibril hosting annular pores.

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

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

Nucleation of β-rich oligomers and β-barrels in the early aggregation of human islet amyloid polypeptide

The self-assembly of human islet amyloid polypeptide (hIAPP) into β-sheet rich amyloid aggregates is associated with pancreatic β-cell death in type 2 diabetes (T2D). Prior experimental studies of hIAPP aggregation reported the early accumulation of α-helical intermediates before the rapid conversion into β-sheet rich amyloid fibrils, as also corroborated by our experimental characterizations with transmission electron microscopy and Fourier transform infrared spectroscopy. Although increasing evidence suggests that small oligomers populating early hIAPP aggregation play crucial roles in cytotoxicity, structures of these oligomer intermediates and their conformational conversions remain unknown, hindering our understanding of T2D disease mechanism and therapeutic design targeting these early aggregation species. We further applied large-scale discrete molecule dynamics simulations to investigate the oligomerization of full-length hIAPP, employing multiple molecular systems of increasing number of peptides. We found that the oligomerization process was dynamic, involving frequent inter-oligomeric exchanges. On average, oligomers had more α-helices than β-sheets, consistent with ensemble-based experimental measurements. However, in ~4-6% independent simulations, β-rich oligomers expected as the fibrillization intermediates were observed, especially in the pentamer and hexamer simulations. These β-rich oligomers could adopt β-barrel conformations, recently postulated to be the toxic oligomer species but only observed computationally in the aggregates of short amyloid protein fragments. Free-energy analysis revealed high energies of these β-rich oligomers, supporting the nucleated conformational changes of oligomers in amyloid aggregation. β-barrel oligomers of full-length hIAPP with well-defined three-dimensional structures may play an important pathological role in T2D etiology and may be a therapeutic target for the disease.

On the Applicability of Force Fields To Study the Aggregation of Amyloidogenic Peptides Using Molecular Dynamics Simulations

Molecular dynamics simulations play an essential role in understanding biomolecular processes such as protein aggregation at temporal and spatial resolutions which are not attainable by experimental methods. For a correct modeling of protein aggregation, force fields must accurately represent molecular interactions. Here, we study the effect of five different force fields on the oligomer formation of Alzheimer's Aβ_16-22 peptide and two of its mutants: Aβ_16-22(F19V,F20V), which does not form fibrils, and Aβ_16-22(F19L) which forms fibrils faster than the wild type. We observe that while oligomer formation kinetics depends strongly on the force field, structural properties, such as the most relevant protein-protein contacts, are similar between them. The oligomer formation kinetics obtained with different force fields differ more from each other than the kinetics between aggregating and nonaggregating peptides simulated with a single force field. We discuss the difficulties in comparing atomistic simulations of amyloid oligomer formation with experimental observables.

β-barrel Oligomers as Common Intermediates of Peptides Self-Assembling into Cross-β Aggregates

Oligomers populated during the early amyloid aggregation process are more toxic than mature fibrils, but pinpointing the exact toxic species among highly dynamic and heterogeneous aggregation intermediates remains a major challenge. β-barrel oligomers, structurally-determined recently for a slow-aggregating peptide derived from αB crystallin, are attractive candidates for exerting amyloid toxicity due to their well-defined structures as therapeutic targets and compatibility to the "amyloid-pore" hypothesis of toxicity. To assess whether β-barrel oligomers are common intermediates to amyloid peptides - a necessary step toward associating β-barrel oligomers with general amyloid cytotoxicity, we computationally studied the oligomerization and fibrillization dynamics of seven well-studied fragments of amyloidogenic proteins with different experimentally-determined aggregation morphologies and cytotoxicity. In our molecular dynamics simulations, β-barrel oligomers were only observed in five peptides self-assembling into the characteristic cross-β aggregates, but not the other two that formed polymorphic β-rich aggregates as reported experimentally. Interestingly, the latter two peptides were previously found nontoxic. Hence, the observed correlation between β-barrel oligomers formation and cytotoxicity supports the hypothesis of β-barrel oligomers as the common toxic intermediates of amyloid aggregation.

Mechanistic insights into the inhibition and size effects of graphene oxide nanosheets on the aggregation of an amyloid-β peptide fragment

The aggregation of amyloid-β (Aβ), which involves the formation of small oligomers and mature fibrils, has received considerable attention in the past few decades due to its close link with Alzheimer's disease (AD). The inhibition of β-sheet formation has been considered as the primary therapeutic strategy for AD. In this respect, graphene oxide (GO) has gained significant attention because of its high solubility, good biocompatibility and inhibitory effect on the aggregation of Aβ and the 33-42 fragment (Aβ33-42). However, the inhibitory mechanism at the atomic level remains elusive. Herein, we investigated the oligomerization of Aβ33-42 by performing replica exchange molecular dynamics simulations on four Aβ33-42 peptide chains in the absence and presence of two different sizes of GO. Our simulations show that isolated Aβ33-42 can form fibril-prone extended β-sheets and barrel-like structures, whereas they are suppressed in the presence of GO nanosheets. Our data reveal that GO inhibits Aβ33-42 oligomerization by making Aβ33-42 peptides separate from each other through strong interactions with M35. With the same total number of atoms, GO120 displays better inhibitory effect than GO60 by providing a larger effective contact surface area. This study provides the molecular mechanism of GO in inhibiting the aggregation of Aβ33-42, which might offer a theoretical insight into the design of drugs against AD at the atomic level.

Structures and dynamics of β-barrel oligomer intermediates of amyloid-beta16-22 aggregation

Accumulating evidence suggests that soluble oligomers are more toxic than final fibrils of amyloid aggregations. Among the mixture of inter-converting intermediates with continuous distribution of sizes and secondary structures, oligomers in the β-barrel conformation - a common class of protein folds with a closed β-sheet - have been postulated as the toxic species with well-defined three-dimensional structures to perform pathological functions. A common mechanism for amyloid toxicity, therefore, implies that all amyloid peptides should be able to form β-barrel oligomers as the aggregation intermediates. Here, we applied all-atom discrete molecular dynamics (DMD) simulations to evaluate the formation of β-barrel oligomers and characterize their structures and dynamics in the aggregation of a seven-residue amyloid peptide, corresponding to the amyloid core of amyloid-β with a sequence of ¹⁶KLVFFAE²² (Aβ16-22). We carried out aggregation simulations with various numbers of peptides to study the size dependence of aggregation dynamics and assembly structures. Consistent with previous computational studies, we observed the formation of β-barrel oligomers in all-atom DMD simulations. Using a network-based approach to automatically identify β-barrel conformations, we systematically characterized β-barrels of various sizes. Our simulations revealed the conformational inter-conversion between β-barrels and double-layer β-sheets due to increased structural strains upon forming a closed β-barrel while maximizing backbone hydrogen bonds. The potential of mean force analysis further characterized the free energy barriers between these two states. The obtained structural and dynamic insights of β-barrel oligomers may help better understand the molecular mechanism of oligomer toxicities and design novel therapeutics targeting the toxic β-barrel oligomers. This article is part of a Special Issue entitled: Protein Aggregation and Misfolding at the Cell Membrane Interface edited by Ayyalusamy Ramamoorthy.

Resolving the Atomistic Modes of Anle138b Inhibitory Action on Peptide Oligomer Formation

The diphenyl-pyrazole compound anle138b is a known inhibitor of oligomeric aggregate formation in vitro and in vivo. Therefore, anle138b is considered a promising drug candidate to beneficially interfere with neurodegenerative processes causing devastating pathologies in humans. The atomistic details of the aggregation inhibition mechanism, however, are to date unknown since the ensemble of small nonfibrillar aggregates is structurally heterogeneous and inaccessible to direct structural characterization. Here, we set out to elucidate anle138b's mode of action using all-atom molecular dynamics simulations on the multi-microsecond time scale. By comparing simulations of dimeric to tetrameric aggregates from fragments of four amyloidogenic proteins (Aβ, hTau40, hIAPP, and Sup35N) in the presence and absence of anle138b, we show that the compound reduces the overall number of intermolecular hydrogen bonds, disfavors the sampling of the aggregated state, and remodels the conformational distributions within the small oligomeric peptide aggregates. Most notably, anle138b preferentially interacts with the disordered structure ensemble via its pyrazole moiety, thereby effectively blocking interpeptide main chain interactions and impeding the spontaneous formation of ordered β-sheet structures, in particular those with out-of-register antiparallel β-strands. The structurally very similar compound anle234b was previously identified as inactive by in vitro experiments. Here, we show that anle234b has no significant effect on the aggregation process in terms of reducing the β-structure content. Moreover, we demonstrate that the hydrogen bonding capabilities are autoinhibited due to steric effects imposed by the molecular geometry of anle234b and thereby indirectly confirm the proposed inhibitory mechanism of anle138b. We anticipate that the prominent binding of anle138b to partially disordered and dynamical aggregate structures is a generic basis for anle138b's ability to suppress toxic oligomer formation in a wide range of amyloidogenic peptides and proteins.

How accurate are your simulations? Effects of confined aqueous volume and AMBER FF99SB and CHARMM22/CMAP force field parameters on structural ensembles of intrinsically disordered proteins: Amyloid-β42 in water

Amyloid-β₄₂ (Aβ₄₂) is an intrinsically disordered peptide intimately related to the pathogenesis of several neurodegenerative diseases. Molecular dynamics (MD) simulations are extensively utilized in the characterization of the structures and conformational dynamics of intrinsically disordered proteins (IDPs) including Aβ₄₂, with AMBER and CHARMM parameters being commonly used in these studies. Recently, comparison of the effects of force field parameters on the Aβ₄₂ structures has started to gain significant attention. In this study, the structures of Aβ₄₂ are simulated using AMBER FF99SB and CHARMM22/CMAP parameters via replica exchange MD simulations utilizing a widely used clustering algorithm. These analyses show that the structural properties (extent and positioning of the elements of secondary and tertiary structure), radius of gyration values, number and position of salt bridges are extremely dependent on the chosen force field parameters notably with the usage of clustering algorithms. For example, predicted secondary structure elements, which are of the great importance for better understanding of the molecular mechanisms of neurodegenerative diseases, deviate enormously in models generated using currently available force field parameters for proteins. Based on the derived models, chemical shift values are calculated and compared to the experimentally determined data. This comparison revealed that although both force field parameters yield results in agreement with experiments, the obtained structural properties were rather different using a clustering algorithm. In other words, these results show that the predicted structures depend heavily on the force field parameters. Importantly, since none of the force field parameters currently utilized in MD studies were developed specifically taking into account the disordered nature of IDPs, these findings clearly indicate that new force field parameters have to be developed for IDPs considering their rapid flexibility and dynamics with high amplitude. Furthermore, molecular simulations of IDPs are typically conducted using one water volume. We show that the confined aqueous volume impacts the predicted structural properties of Aβ₄₂ in water. Although up to date, confined aqueous volume effects have been ignored in the MD simulations of IDPs in water, our data indicate that these effects have to be taken into account in predicting the structural and thermodynamic properties of disordered proteins in solution.

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.

Probing oligomerization of amyloid beta peptide in silico

Aggregation of amyloid β (Aβ) peptide is implicated in fatal Alzheimer's disease, for which no cure is available. Understanding the mechanisms responsible for this aggregation is required in order for therapies to be developed. In an effort to better understand the molecular mechanisms involved in spontaneous aggregation of Aβ peptide, extensive molecular dynamics simulations are reported, and the results are analyzed through a combination of structural biology tools and a novel essential collective dynamics method. Several model systems composed of ten or twelve Aβ_17-42 chains in water are investigated, and the influence of metal ions is probed. The results suggest that Aβ monomers tend to aggregate into stable globular-like oligomers with 13-23% of β-sheet content. Two stages of oligomer formation have been identified: quick collapse within the first 40 ns of the simulation, characterized by a decrease in inter-chain separation and build-up of β-sheets, and the subsequent slow relaxation of the oligomer structure. The resulting oligomers comprise a stable, coherently moving sub-aggregate of 6-9 strongly inter-correlated chains. Cu²⁺ and Fe²⁺ ions have been found to develop coordination bonds with carboxylate groups of E22, D23 and A42, which remain stable during 200 ns simulations. The presence of Fe²⁺, and particularly Cu²⁺ ions, in negatively charged cavities has been found to cause significant changes in the structure and dynamics of the oligomers. The results indicate, in particular, that formation of non-fibrillar oligomers might be involved in early template-free aggregation of Aβ_17-42 monomers, with charged species such as Cu²⁺ or Fe²⁺ ions playing an important role.

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.

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.

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.

Peptides Composed of Alternating L- and D-Amino Acids Inhibit Amyloidogenesis in Three Distinct Amyloid Systems Independent of Sequence

There is now substantial evidence that soluble oligomers are primary toxic agents in amyloid diseases. The development of an antibody recognizing the toxic soluble oligomeric forms of different and unrelated amyloid species suggests a common conformational intermediate during amyloidogenesis. We previously observed a common occurrence of a novel secondary structure element, which we call α-sheet, in molecular dynamics (MD) simulations of various amyloidogenic proteins, and we hypothesized that the toxic conformer is composed of α-sheet structure. As such, α-sheet may represent a conformational signature of the misfolded intermediates of amyloidogenesis and a potential unique binding target for peptide inhibitors. Recently, we reported the design and characterization of a novel hairpin peptide (α1 or AP90) that adopts stable α-sheet structure and inhibits the aggregation of the β-Amyloid Peptide Aβ42 and transthyretin. AP90 is a 23-residue hairpin peptide featuring alternating D- and L-amino acids with favorable conformational propensities for α-sheet formation, and a designed turn. For this study, we reverse engineered AP90 to identify which of its design features is most responsible for conferring α-sheet stability and inhibitory activity. We present experimental characterization (CD and FTIR) of seven peptides designed to accomplish this. In addition, we measured their ability to inhibit aggregation in three unrelated amyloid species: Aβ42, transthyretin, and human islet amylin polypeptide. We found that a hairpin peptide featuring alternating L- and D-amino acids, independent of sequence, is sufficient for conferring α-sheet structure and inhibition of aggregation. Additionally, we show a correlation between α-sheet structural stability and inhibitory activity.

Critical Nucleus Structure and Aggregation Mechanism of the C-terminal Fragment of Copper-Zinc Superoxide Dismutase Protein

The aggregation of the copper-zinc superoxide dismutase (SOD1) protein is linked to familial amyotrophic lateral sclerosis, a progressive neurodegenerative disease. A recent experimental study has shown that the (147)GVIGIAQ(153) SOD1 C-terminal segment not only forms amyloid fibrils in isolation but also accelerates the aggregation of full-length SOD1, while substitution of isoleucine at site 149 by proline blocks its fibril formation. Amyloid formation is a nucleation-polymerization process. In this study, we investigated the oligomerization and the nucleus structure of this heptapeptide. By performing extensive replica-exchange molecular dynamics (REMD) simulations and conventional MD simulations, we found that the GVIGIAQ hexamers can adopt highly ordered bilayer β-sheets and β-barrels. In contrast, substitution of I149 by proline significantly reduces the β-sheet probability and results in the disappearance of bilayer β-sheet structures and the increase of disordered hexamers. We identified mixed parallel-antiparallel bilayer β-sheets in both REMD and conventional MD simulations and provided the conformational transition from the experimentally observed parallel bilayer sheets to the mixed parallel-antiparallel bilayer β-sheets. Our simulations suggest that the critical nucleus consists of six peptide chains and two additional peptide chains strongly stabilize this critical nucleus. The stabilized octamer is able to recruit additional random peptides into the β-sheet. Therefore, our simulations provide insights into the critical nucleus formation and the smallest stable nucleus of the (147)GVIGIAQ(153) peptide.

Amyloid β-Protein C-Terminal Fragments: Formation of Cylindrins and β-Barrels

In order to evaluate potential therapeutic targets for treatment of amyloidoses such as Alzheimer's disease (AD), it is essential to determine the structures of toxic amyloid oligomers. However, for the amyloid β-protein peptide (Aβ), thought to be the seminal neuropathogenetic agent in AD, its fast aggregation kinetics and the rapid equilibrium dynamics among oligomers of different size pose significant experimental challenges. Here we use ion-mobility mass spectrometry, in combination with electron microscopy, atomic force microscopy, and computational modeling, to test the hypothesis that Aβ peptides can form oligomeric structures resembling cylindrins and β-barrels. These structures are hypothesized to cause neuronal injury and death through perturbation of plasma membrane integrity. We show that hexamers of C-terminal Aβ fragments, including Aβ(24-34), Aβ(25-35) and Aβ(26-36), have collision cross sections similar to those of cylindrins. We also show that linking two identical fragments head-to-tail using diglycine increases the proportion of cylindrin-sized oligomers. In addition, we find that larger oligomers of these fragments may adopt β-barrel structures and that β-barrels can be formed by folding an out-of-register β-sheet, a common type of structure found in amyloid proteins.

Oligomer Formation of Toxic and Functional Amyloid Peptides Studied with Atomistic Simulations

Amyloids are associated with diseases, including Alzheimer's, as well as functional roles such as storage of peptide hormones. It is still unclear what differences exist between aberrant and functional amyloids. However, it is known that soluble oligomers formed during amyloid aggregation are more toxic than the final fibrils. Here, we perform molecular dynamics simulations to study the aggregation of the amyloid-β peptide Aβ25-35, associated with Alzheimer's disease, and two functional amyloid-forming tachykinin peptides: kassinin and neuromedin K. Although the three peptides have similar primary sequences, tachykinin peptides, in contrast to Aβ25-35, form nontoxic amyloids. Our simulations reveal that the charge of the C-terminus is essential to controlling the aggregation process. In particular, when the kassinin C-terminus is not amidated, the aggregation kinetics decreases considerably. In addition, we observe that the monomeric peptides in extended conformations aggregate faster than those in collapsed hairpin-like conformations.

Energetic contributions of residues to the formation of early amyloid-β oligomers

Low-weight amyloid-β (Aβ) oligomers formed at early stages of oligomerization rather than fibril assemblies seem to be the toxic components that drive neurodegeneration in Alzheimer's disease. Unfortunately, detailed knowledge of the structure of these early oligomers at the residue level is not yet available. In this study, we performed all-atom explicit solvent molecular dynamics simulations to examine the oligomerization process of Aβ10-35 monomers when forming dimers, trimers, tetramers and octamers, with four independent simulations of a total simulated time of 3 μs for each oligomer system. The decomposition of the stability free energy by MM-GBSA methodology allowed us to unravel the network of energetic interactions that stabilize such oligomers. The contribution of the intermonomeric van der Waals term is the most significant energy feature of the oligomerization process, consistent with the so-called hydrophobic effect. Furthermore, the decomposition of the stability free energy into residues and residue-pairwise terms revealed that it is mainly apolar interactions between the three specific hydrophobic fragments 31-35 (C-terminal region), 17-20 (central hydrophobic core) and 12-14 (N-terminal region) that are responsible for such a favourable effect. The conformation in which the hydrophobic cthr-chc interaction is oriented perpendicularly is particularly important. We propose three other model substructures that favour the oligomerization process and can thus be considered as molecular targets for future inhibitors. Understanding Aβ oligomerization at the residue level could lead to more efficient design of inhibitors of this process.

Stability of amyloid oligomers

Molecular simulations are now commonly used to complement experimental techniques in investigating amyloids and their role in human diseases. In this chapter, we will summarize techniques and approaches often used in amyloid simulations and will present recent success stories. Our examples will be focused on lessons learned from molecular dynamics simulations in aqueous environments that start from preformed aggregates. These studies explore the limitations that arise from the choice of force field, the role of mutations in the growth of amyloid aggregates, segmental polymorphism, and the importance of cross-seeding. Furthermore, they give evidence for potential toxicity mechanisms. We finally discuss the role of molecular simulations in the search for aggregation inhibitors.

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.

Acidic pH retards the fibrillization of human Islet Amyloid Polypeptide due to electrostatic repulsion of histidines

The human Islet Amyloid Polypeptide (hIAPP) is the major constituent of amyloid deposits in pancreatic islets of type-II diabetes. IAPP is secreted together with insulin from the acidic secretory granules at a low pH of approximately 5.5 to the extracellular environment at a neutral pH. The increased accumulation of extracellular hIAPP in diabetes indicates that changes in pH may promote amyloid formation. To gain insights and underlying mechanisms of the pH effect on hIAPP fibrillogenesis, all-atom molecular dynamics simulations in explicit solvent model were performed to study the structural properties of five hIAPP protofibrillar oligomers, under acidic and neutral pH, respectively. In consistent with experimental findings, simulation results show that acidic pH is not conducive to the structural stability of these oligomers. This provides a direct evidence for a recent experiment [L. Khemtemourian, E. Domenech, J. P. F. Doux, M. C. Koorengevel, and J. A. Killian, J. Am. Chem. Soc. 133, 15598 (2011)], which suggests that acidic pH inhibits the fibril formation of hIAPP. In addition, a complementary coarse-grained simulation shows the repulsive electrostatic interactions among charged His18 residues slow down the dimerization process of hIAPP by twofold. Besides, our all-atom simulations reveal acidic pH mainly affects the local structure around residue His18 by destroying the surrounding hydrogen-bonding network, due to the repulsive interactions between protonated interchain His18 residues at acidic pH. It is also disclosed that the local interactions nearby His18 operating between adjacent β-strands trigger the structural transition, which gives hints to the experimental findings that the rate of hIAPP fibril formation and the morphologies of the fibrillar structures are strongly pH-dependent.

Understanding amyloid fibril nucleation and aβ oligomer/drug interactions from computer simulations

Evolution has fine-tuned proteins to accomplish a variety of tasks. Yet, with aging, some proteins assemble into harmful amyloid aggregates associated with neurodegenerative diseases, such as Alzheimer's disease (AD), which presents a complex and costly challenge to our society. Thus, far, drug after drug has failed to slow the progression of AD, characterized by the self-assembly of the 39-43 amino acid β-amyloid (Aβ) protein into extracellular senile plaques that form a cross-β structure. While there is experimental evidence that the Aβ small oligomers are the primary toxic species, standard tools of biology have failed to provide structures of these transient, inhomogeneous assemblies. Despite extensive experimental studies, researchers have not successfully characterized the nucleus ensemble, the starting point for rapid fibril formation. Similarly scientists do not have atomic data to show how the compounds that reduce both fibril formation and toxicity in cells bind to Aβ42 oligomers. In this context, computer simulations are important tools for gaining insights into the self-assembly of amyloid peptides and the molecular mechanism of inhibitors. This Account reviews what analytical models and simulations at different time and length scales tell us about the dynamics, kinetics, and thermodynamics of amyloid fibril formation and, notably, the nucleation process. Though coarse-grained and mesoscopic protein models approximate atomistic details by averaging out unimportant degrees of freedom, they provide generic features of amyloid formation and insights into mechanistic details of the self-assembly process. The thermodynamics and kinetics vary from linear peptides adopting straight β-strands in fibrils to longer peptides adopting in parallel U shaped conformations in fibrils. In addition, these properties change with the balance between electrostatic and hydrophobic interactions and the intrinsic disorder of the system. However, simulations suggest that the critical nucleus size might be on the order of 20 chains under physiological conditions. The transition state might be characterized by a simultaneous change from mixed antiparallel/parallel β-strands with random side-chain packing to the final antiparallel or parallel states with the steric zipper packing of the side chains. Second, we review our current computer-based knowledge of the 3D structures of inhibitors with Aβ42 monomer and oligomers, a prerequisite for developing new drugs against AD. Recent extensive all-atom simulations of Aβ42 dimers with known inhibitors such as the green tea compound epigallocatechin-3-gallate and 1,4-naphthoquinon-2-yl-l-tryptophan provide a spectrum of initial Aβ42/inhibitor structures useful for screening and drug design. We conclude by discussing future directions that may offer opportunities to fully understand nucleation and further AD drug development.

Identification of small-molecule binding pockets in the soluble monomeric form of the Aβ42 peptide

The aggregation of intrinsically disordered peptides and proteins is associated with a wide range of highly debilitating neurological and systemic disorders. In this work we explored the potential of a structure-based drug discovery procedure to target one such system, the soluble monomeric form of the Aβ42 peptide. We utilised for this purpose a set of structures of the Aβ42 peptide selected from clusters of conformations within an ensemble generated by molecular dynamics simulations. Using these structures we carried out fragment mapping calculations to identify binding "hot spots" on the monomeric form of the Aβ42 peptide. This procedure provided a set of hot spots with ligand efficiencies comparable to those observed for structured proteins, and clustered into binding pockets. Such binding pockets exhibited a propensity to bind small molecules known to interact with the Aβ42 peptide. Taken together these results provide an initial indication that fragment-based drug discovery may represent a potential therapeutic strategy for diseases associated with the aggregation of intrinsically disordered proteins.

Conformational ensemble and polymorphism of the all-atom Alzheimer's Aβ(37-42) amyloid peptide oligomers

Although the Aβ(37-42) peptide has two opposite terminal charges, counterintuitively its current fibril amyloid structure reveals in register parallel β-strands, as formed by the full length Aβ peptide. In this study, we carried out a replica exchange molecular dynamics simulation of 16 all-atom Aβ(37-42) peptides in explicit water starting from randomized and dispersed chains. The extensive conformational sampling (48 replicas, 460 ns/replica) with a total simulation time of 23 μs allows us to obtain a full picture on the equilibrium conformational distribution of oligomers and β-sheet sizes and gain some insights into the oligomerization process at 300 K. At the peptide concentration of 12 mM, self-assembly is described by the condensation-polymerization mechanism with conversion from micelle-like to high β-sheet structures. At equilibrium, the oligomer distribution consists of large aggregates and free monomers, representing 70% and 25% of all species, respectively. Though the formation/dissociation of β-strand is high, the population of 4-5 fully parallel β-strands, consistent with the arrangement in the current fibril, is marginal and that of 4-5 fully antiparallel β-strands, consistent with amyloid polymorphism, is non-negligible. However, the system adopts essentially mixed parallel/antiparallel β-strands. This indicates that a system of 16 Aβ(37-42) chains in explicit solvent still does not form more stable species that will irreversibly grow to a fibril, independently of polymorphism. Our results also suggest that the Aβ(37-42) fibril may display packing polymorphism with antiparallel β-strands, in addition to the experimentally observed in register parallel β-strands.

Enhanced sampling algorithms

In biomolecular systems (especially all-atom models) with many degrees of freedom such as proteins and nucleic acids, there exist an astronomically large number of local-minimum-energy states. Conventional simulations in the canonical ensemble are of little use, because they tend to get trapped in states of these energy local minima. Enhanced conformational sampling techniques are thus in great demand. A simulation in generalized ensemble performs a random walk in potential energy space and can overcome this difficulty. From only one simulation run, one can obtain canonical-ensemble averages of physical quantities as functions of temperature by the single-histogram and/or multiple-histogram reweighting techniques. In this article we review uses of the generalized-ensemble algorithms in biomolecular systems. Three well-known methods, namely, multicanonical algorithm, simulated tempering, and replica-exchange method, are described first. Both Monte Carlo and molecular dynamics versions of the algorithms are given. We then present various extensions of these three generalized-ensemble algorithms. The effectiveness of the methods is tested with short peptide and protein systems.

Coarse-grained models for protein folding and aggregation

Coarse-grained models for protein folding and aggregation are used to explore large dimension scales and timescales that are inaccessible to all-atom models in explicit aqueous solution. Combined with enhanced configuration search methods, these simplified models with various levels of granularity offer the possibility to determine equilibrium structures, compare folding kinetics and thermodynamics with experiments for single proteins and understand the dynamic assembly of amyloid proteins leading to neurodegenerative diseases. I shall describe recent progress in developing such models, and discuss their potentials and limitations in probing the folding and misfolding of proteins with computer simulations.

Out-of-register β-sheets suggest a pathway to toxic amyloid aggregates

Although aberrant protein aggregation has been conclusively linked to dozens of devastating amyloid diseases, scientists remain puzzled about the molecular features that render amyloid fibrils or small oligomers toxic. Here, we report a previously unobserved type of amyloid fibril that tests as cytotoxic: one in which the strands of the contributing β-sheets are out of register. In all amyloid fibrils previously characterized at the molecular level, only in-register β-sheets have been observed, in which each strand makes its full complement of hydrogen bonds with the strands above and below it in the fibril. In out-of-register sheets, strands are sheared relative to one another, leaving dangling hydrogen bonds. Based on this finding, we designed out-of-register β-sheet amyloid mimics, which form both cylindrin-like oligomers and fibrils, and these mimics are cytotoxic. Structural and energetic considerations suggest that out-of-register fibrils can readily convert to toxic cylindrins. We propose that out-of-register β-sheets and their related cylindrins are part of a toxic amyloid pathway, which is distinct from the more energetically favored in-register amyloid pathway.

Structural, thermodynamical, and dynamical properties of oligomers formed by the amyloid NNQQ peptide: insights from coarse-grained simulations

Characterizing the early formed oligomeric intermediates of amyloid peptides is of particular interest due to their links with neurodegenerative diseases. Here we study the NNQQ peptide, known to display parallel β-strands in amyloid fibrils by x-ray microcrystallography, and investigate the structural, thermodynamical, and dynamical properties of 20 NNQQ peptides using molecular dynamics and replica exchange molecular dynamics simulations coupled to a coarse-grained force field. All simulations are initiated from randomized and fully dispersed monomeric conformations. Our simulations reveal that the phase transition is characterized by a change in the oligomer and β-sheet size distributions and the percentage of mixed parallel/antiparallel β-strands when the sheets are formed. At all temperatures, however, the fraction of parallel β-strands remains low, though there are many association/fragmentation events. This work and a growing body of computational studies provide strong evidence that the critical nucleus goes beyond 20 chains and reordering of the β-strands occurs in larger oligomers.

Nanomaterials design and tests for neural tissue engineering

Nanostructured scaffolds recently showed great promise in tissue engineering: nanomaterials can be tailored at the molecular level and scaffold morphology may more closely resemble features of extracellular matrix components in terms of porosity, framing and biofunctionalities. As a consequence, both biomechanical properties of scaffold microenvironments and biomaterial-protein interactions can be tuned, allowing for improved transplanted cell engraftment and better controlled diffusion of drugs. Easier said than done, a nanotech-based regenerative approach encompasses different fields of know-how, ranging from in silico simulations, nanomaterial synthesis and characterization at the nano-, micro- and mesoscales to random library screening methods (e.g. phage display), in vitro cellular-based experiments and validation in animal models of the target injury. All of these steps of the "assembly line" of nanostructured scaffolds are tightly interconnected both in their standard analysis techniques and in their most recent breakthroughs: indeed their efforts have to jointly provide the deepest possible analyses of the diverse facets of the challenging field of neural tissue engineering. The purpose of this review is therefore to provide a critical overview of the recent advances in and drawbacks and potential of each mentioned field, contributing to the realization of effective nanotech-based therapies for the regeneration of peripheral nerve transections, spinal cord injuries and brain traumatic injuries. Far from being the ultimate overview of such a number of topics, the reader will acknowledge the intrinsic complexity of the goal of nanotech tissue engineering for a conscious approach to the development of a regenerative therapy and, by deciphering the thread connecting all steps of the research, will gain the necessary view of its tremendous potential if each piece of stone is correctly placed to work synergically in this impressive mosaic.