Elongation of ordered peptide aggregate of an amyloidogenic hexapeptide NFGAIL observed in molecular dynamics simulations with explicit solvent

The mechanisms by which amyloidogenic peptides and proteins form soluble toxic oligomers remain elusive. We have studied the formation of partially ordered tetramers and well-ordered octamers of an amyloidogenic hexapeptide NFGAIL (residues 22-27 of the human islet amyloid polypeptide) in our previous work. Continuing the effort, we here probe the beta-sheet elongation process by a combined total of 2.0 micros molecular dynamics simulations with explicit solvent. In a set of 10 simulations with the peptides restrained to the extended conformation, we observed that the main growth mode was elongation along the beta-sheet hydrogen bonds through primarily a two-stage process. Driven by hydrophobic forces, the peptides initially attached to the surface of the ordered oligomer, moved quickly to the beta-sheet edges, and formed stable beta-sheet hydrogen bonds. Addition of peptides to the existing oligomer notably improved the order of the peptide aggregate in which labile outer layer beta-sheets were stabilized, which provides good templates for further elongation. These simulations suggested that elongation along the beta-sheet hydrogen bonds occurs at the intermediate stage when low-weight oligomers start to form. We did not observe significant preference toward either parallel or antiparallel beta-sheets at the elongation stage for this peptide. In another set of 10 unrestrained simulations, the dominant growth mode was disordered aggregation. Taken together, these results offered a glimpse at the molecular events leading to the formation of ordered and disordered low-weight oligomers.     

Mechanism of unassisted ion transport across membrane bilayers

To establish how charged species move from water to the nonpolar membrane interior and to determine the energetic and structural effects accompanying this process, we performed molecular dynamics simulations of the transport of Na+ and Cl- across a lipid bilayer located between two water lamellae. The total length of molecular dynamics trajectories generated for each ion was 10 ns. Our simulations demonstrate that permeation of ions into the membrane is accompanied by the formation of deep, asymmetric thinning defects in the bilayer, whereby polar lipid head groups and water penetrate the nonpolar membrane interior. Once the ion crosses the midplane of the bilayer the deformation "switches sides"; the initial defect slowly relaxes, and a defect forms in the outgoing side of the bilayer. As a result, the ion remains well solvated during the process; the total number of oxygen atoms from water and lipid head groups in the first solvation shell remains constant. A similar membrane deformation is formed when the ion is instantaneously inserted into the interior of the bilayer. The formation of defects considerably lowers the free energy barrier to transfer of the ion across the bilayer and, consequently, increases the permeabilities of the membrane to ions, compared to the rigid, planar structure, by approximately 14 orders of magnitude. Our results have implications for drug delivery using liposomes and peptide insertion into membranes.     

Molecular insights into the surface morphology, layering structure, and aggregation kinetics of surfactant-stabilized graphene dispersions

The production of graphene with open band gaps for the manufacturing of graphene-based electronic and optical devices requires synthesis methods to either control the number of layers to enrich AB-stacked bilayer or trilayer graphene or control the extent of functionalization of monolayer graphene. Solution-phase dispersion of graphene is promising for both methods to create printable electronics and nanocomposites. However, both methods face common challenges, including controlling the surface morphology, reducing the turbostratic layering, and enhancing the dispersion stability. To address these challenges at the molecular level, we successfully combined molecular simulations, theoretical modeling, and experimental measurements. First, we probed the surface structure and electrostatic potential of monolayer graphene dispersed in a sodium cholate (SC) surfactant aqueous solution, which exhibits 2D sheets partially covered with a monolayer of negatively charged cholate ions. Similar to the case of carbon nanotube functionalization, one may regulate the binding affinity of charged reactants for graphene functionalization by manipulating the surface morphology. Subsequently, we quantified the interactions between two graphene-surfactant assemblies by calculating the potential of mean force (PMF) between two surfactant-covered graphene sheets, which confirmed the existence of a metastable bilayer graphene structure due to the steric hindrance of the confined surfactant molecules. The traditional Derjaguin-Landau-Verwey-Overbeek (DLVO) theory was found to be adequate to explain the long-range electrostatic repulsions between the ionic surfactant-covered graphene sheets but was unable to account for the dominant, short-range steric hindrance imparted by the confined surfactant molecules. Interestingly, one faces a dilemma when using surfactants to disperse and stabilize graphene in aqueous solution: on the one hand, surfactants can stabilize graphene aqueous dispersions, but on the other hand, they prevent the formation of new AB-stacked bilayer and trilayer graphene resulting from the reaggregation process. Finally, the lifetime and time-dependent distribution of various graphene layer types were predicted using a kinetic model of colloid aggregation, and each graphene layer type was further decomposed into subtypes, including the AB-stacked species and various turbostratic species. The kinetic model of colloid aggregation developed here can serve as a useful tool to evaluate the quality of graphene dispersions for subsequent substrate-transferring or functionalization processes.     

Coarse-grained protein molecular dynamics simulations

A limiting factor in biological science is the time-scale gap between experimental and computational trajectories. At this point, all-atom explicit solvent molecular dynamics (MD) are clearly too expensive to explore long-range protein motions and extract accurate thermodynamics of proteins in isolated or multimeric forms. To reach the appropriate time scale, we must then resort to coarse graining. Here we couple the coarse-grained OPEP model, which has already been used with activated methods, to MD simulations. Two test cases are studied: the stability of three proteins around their experimental structures and the aggregation mechanisms of the Alzheimer's Abeta16-22 peptides. We find that coarse-grained isolated proteins are stable at room temperature within 50 ns time scale. Based on two 220 ns trajectories starting from disordered chains, we find that four Abeta16-22 peptides can form a three-stranded beta sheet. We also demonstrate that the reptation move of one chain over the others, first observed using the activation-relaxation technique, is a kinetically important mechanism during aggregation. These results show that MD-OPEP is a particularly appropriate tool to study qualitatively the dynamics of long biological processes and the thermodynamics of molecular assemblies.     

Self-assembly of amylin(20-29) amide-bond derivatives into helical ribbons and peptide nanotubes rather than fibrils

Uncontrolled aggregation of proteins or polypeptides can be detrimental for normal cellular processes in healthy organisms. Proteins or polypeptides that form these amyloid deposits differ in their primary sequence but share a common structural motif: the (anti)parallel beta sheet. A well-accepted approach for interfering with beta-sheet formation is the design of soluble beta-sheet peptides to disrupt the hydrogen-bonding network; this ultimately leads to the disassembly of the aggregates or fibrils. Here, we describe the synthesis, spectroscopic analysis, and aggregation behavior, imaged by electron microscopy, of several backbone-modified amylin(20-29) derivatives. It was found that these amylin derivatives were not able to form fibrils and to some extent were able to inhibit fibril growth of native amylin(20-29). However, two of the amylin peptides were able to form large supramolecular assemblies, like helical ribbons and peptide nanotubes, in which beta-sheet formation was clearly absent. This was quite unexpected since these peptides have been designed as soluble beta-sheet breakers for disrupting the characteristic hydrogen-bonding network of (anti)parallel beta sheets. The increased hydrophobicity and the presence of essential amino acid side chains in the newly designed amylin(20-29) derivatives were found to be the driving force for self-assembly into helical ribbons and peptide nanotubes. This example of controlled and desired peptide aggregation may be a strong impetus for research on bionanomaterials in which special shapes and assemblies are the focus of interest.     

A molecular dynamics simulation study of C60 fullerenes inside a dimyristoylphosphatidylcholine lipid bilayer

We have carried out atomistic molecular dynamics simulations of C60 fullerenes inside a dimyristoylphosphatidylcholine lipid bilayer and an alkane melt. Simulations reveal that the preferred position of a single C60 fullerene is about 6-7 A off of the center plane, allowing the fullerene to take advantage of strong dispersion interactions with denser regions of the bilayer. Further displacement (>8 A) of the fullerene away from the center plane results in a rapid increase in free energy likely due to distortion of the lipid head group layer. The effective interaction between fullerenes (direct interaction plus environment (bilayer)-induced interaction), measured as the potential of mean force (POMF) between two fullerenes as a function of their separation, was found to be significantly less attractive in the lipid bilayer than in an alkane melt of the same molecular weight as the lipid tails. Only part of this difference can be accounted for by the more favorable interaction of the fullerene with the relatively denser bilayer. Additionally, our POMF studies indicate that the bilayer is less able to accommodate the larger aggregated fullerene pair than isolated single fullerenes, again likely due to distortion of the bilayer structure. The implications of these effects on aggregation of fullerenes within lipid bilayer are considered.     

Engineering a beta-sheet protein toward the folding speed limit

Recent experimental studies have shown that alpha-helical proteins can approach the folding "speed limit", where folding switches from an activated to a downhill process in free energy. beta-sheet proteins are generally thought to fold more slowly than helix bundles. However, based on studies of hairpins, folding should still be able to approach the microsecond time scale. Here we demonstrate how the hPin1 WW domain, a triple-stranded beta-sheet protein with a sharp thermodynamic melting transition, can be engineered toward the folding "speed limit" without a significant loss in thermal denaturation cooperativity.     

Beta-hairpin folding by a model amyloid peptide in solution and at an interface

The development of specific agents against amyloidoses requires an understanding of the conformational distribution of fibrillogenic peptides at a microscopic level. Here, I present molecular dynamics simulations of the model amyloid peptide LSFD with sequence LSFDNSGAITIG-NH2 in explicit water and at a water/vapor interface for a total time scale of approximately 1.8 micros. An extended structure was used as initial peptide configuration. At approximately 290 K, solvated LSFD was kinetically trapped in diverse misfolded beta-sheet/coil conformations. At 350 K, in contrast, the same type II' beta-hairpin in equilibrium with less ordered but also U-shaped conformations was observed for the core residues DNSGAITI in solution and at the interface in multiple independent simulations. The most stable structural unit of the beta-hairpin was the two residue turn (GA). The core residues exhibited a well-defined folded state in which the beta-hairpin was stabilized by a hydrogen bond between the side chain of Asn-385 and the main chain carbonyl group of Gly-387. My results suggest that beta-sheet conformations indicated from previous Fourier-transform infrared spectroscopy measurements immediately after preparation of the peptide solution may not arise from protofilaments as speculated by others but are a property of LSFD monomers. In addition, combined with previous results from X-ray scattering, my findings suggest that interfacial aggregation of LSFD implies a transition from U-shaped to extended peptide conformations. This work including the first simulations of reversible beta-hairpin folding at an interface is an essential step toward a microscopic understanding of interfacial peptide folding and self-assembly. Knowledge of the main conformation of the peptide core may facilitate the design of possible inhibitors of LSFD aggregation as a test ground for future computational therapeutic strategies against amyloid diseases.     

Aggregation of lysozyme and of poly(ethylene glycol)-modified lysozyme after adsorption to silica

Surface-induced aggregation is a common instability during protein storage, delivery and purification. This aggregation can lead to the formation of fibrils rich in intermolecular beta-sheet structure. Techniques to probe surface-clustering are limited. Here we use protein intrinsic fluorescence and thioflavin T probe fluorescence in a total internal reflection fluorescence (TIRF) sampling geometry to simultaneously monitor the kinetics of adsorption and aggregation for chicken egg lysozyme on a silica surface. We observe a slow surface-induced aggregation process that continues well after the lysozyme adsorption kinetics have plateaued. The rate of surface-induced aggregation is independent of the lysozyme concentration in solution. Consistent with the clustering observed via thioflavin T fluorescence, infrared amide I band spectra also show a 1.5-fold increase in intermolecular beta-sheet content upon lysozyme adsorption. Tryptophan emission spectra show no evidence for any tertiary structural change upon adsorption. Furthermore, we observe that the covalent modification of lysozyme with a single poly(ethylene glycol) (PEG) grafted chain does not inhibit aggregation on the surface, but a second PEG graft significantly inhibits the intermolecular beta-sheet formation.     

Dissecting Cooperative Communications in a Protein with a High‐Throughput Single‐Molecule Scalpel

Miscued communication often leads to misfolding and aggregation of the proteins involved in many diseases. Owing to the ensemble average property of conventional techniques, detailed communication diagrams are difficult to obtain. Mechanical unfolding affords an unprecedented perspective on cooperative transitions by observing a protein along a trajectory defined by two mutated cysteine residues. Nevertheless, this approach requires tedious sample preparation at the risk of altering native protein conformations. To address these issues, we applied click chemistry to tether a protein to the two dsDNA handles through primary amines in lysine residues as well as at the N terminus. As a proof of concept, we used laser tweezers to mechanically unfold and refold calmodulin along 36 trajectories, maximally allowed by this strategy in a single batch of protein preparation. Without a priori knowledge of the particular residues to which the double‐stranded DNA handles attach, we used hierarchical cluster analysis to identify 20 major trajectories, according to the size and the pattern of unfolding transitions. We dissected the cooperativity into all‐or‐none and partially cooperative events, which represent strong and weak high‐order interactions in proteins, respectively. Although the overall cooperativity is higher within the N or C lobe than that between the lobes, the all‐or‐none cooperativity is anisotropic among different the unfolding trajectories and becomes relatively more predominant when the size of the protein segments increases. The average cooperativity for all‐or‐none transitions falls within the expected range observed by ensemble techniques, which supports the hypothesis that unfolding of a free protein can be reconstituted from individual trajectories.     

ThX – a next-generation probe for the early detection of amyloid aggregates

Neurodegenerative diseases such as Alzheimer''s and Parkinson''s are associated with protein misfolding and aggregation. Recent studies suggest that the small, rare and heterogeneous oligomeric species, formed early on in the aggregation process, may be a source of cytotoxicity. Thioflavin T (ThT) is currently the gold-standard fluorescent probe for the study of amyloid proteins and aggregation processes. However, the poor photophysical and binding properties of ThT impairs the study of oligomers. To overcome this challenge, we have designed Thioflavin X, (ThX), a next-generation fluorescent probe which displays superior properties; including a 5-fold increase in brightness and 7-fold increase in binding affinity to amyloidogenic proteins. As an extrinsic dye, this can be used to study unique structural amyloid features both in bulk and on a single-aggregate level. Furthermore, ThX can be used as a super-resolution imaging probe in single-molecule localisation microscopy. Finally, the improved optical properties (extinction coefficient, quantum yield and brightness) of ThX can be used to monitor structural differences in oligomeric species, not observed via traditional ThT imaging.

Introducing ThX, a next-generation ThT derivative that allows for the early detection of amyloid aggregates at the bulk and single-aggregate levels.     

Inhibition of amyloid fibril formation of human amylin by N-alkylated amino acid and alpha-hydroxy acid residue containing peptides

Amyloid deposits are formed as a result of uncontrolled aggregation of (poly)peptides or proteins. Today several diseases are known, for example Alzheimer's disease, Creutzfeldt-Jakob disease, mad cow disease, in which amyloid formation is involved. Amyloid fibrils are large aggregates of beta-pleated sheets and here a general method is described to introduce molecular mutations in order to achieve disruption of beta-sheet formation. Eight backbone-modified amylin derivatives, an amyloidogenic peptide involved in maturity onset diabetes, were synthesized. Their beta-sheet forming properties were studied by IR spectroscopy and electron microscopy. Modification of a crucial amide NH by an alkyl chain led to a complete loss of the beta-sheet forming capacity of amylin. The resulting molecular mutated amylin derivative could be used to break the beta-sheet thus retarding beta-sheet formation of unmodified amylin. Moreover, it was found that the replacement of this amide bond by an ester moiety suppressed fibrillogenesis significantly. Introduction of N-alkylated amino acids and/or ester functionalities-leading to depsipeptides-into amyloidogenic peptides opens new avenues towards novel peptidic beta-sheet breakers for inhibition of beta-amyloid aggregation.     

Combination of COSMOmic and molecular dynamics simulations for the calculation of membrane–water partition coefficients

The importance of membrane‐water partition coefficients led to the recent extension of the conductor‐like screening model for realistic solvation (COSMO‐RS) to micelles and biomembranes termed COSMOmic. Compared to COSMO‐RS, this new approach needs structural information to account for the anisotropy of colloidal systems. This information can be obtained from molecular dynamics (MD) simulations. In this work, we show that this combination of molecular methods can efficiently be used to predict partition coefficients with good agreement to experimental data and enables screening studies. However, there is a discrepancy between the amount of data generated by MD simulations and the structural information needed for COSMOmic. Therefore, a new scheme is presented to extract data from MD trajectories for COSMOmic calculations. In particular, we show how to calculate the system structure from MD, the influence of lipid conformers, the relation to the COSMOmic layer size, and the water/lipid ratio impact. For a 1,2‐dimyristoyl‐sn‐glycero‐3‐phosphocholine (DMPC) bilayer, 66 partition coefficients for various solutes were calculated. Further, 52 partition coefficients for a 1‐palmitoyl‐2‐oleoyl‐sn‐glycero‐3‐phosphocholine (POPC) bilayer system were calculated. All these calculations were compared to experimental data. © 2013 Wiley Periodicals, Inc.     

Amyloid fibril polymorphism is under kinetic control

Self-assembly of proteins into amyloid aggregates displays a broad diversity of morphologies, both at the protofibrillar and final fibrillar species. These polymorphic species can coexist at fixed experimental conditions, and their relative abundance can be controlled by changing the solvent composition, or stirring the solution. However, the extent to which external conditions regulate the equilibrium of morphologically distinct species is still unknown. Here we investigate the nucleation of distinct fibril morphologies using computer simulations of a simplified model of an amyloid polypeptide. Counterintuitively, the energetically less favorable fibril morphologies nucleate more frequently than the morphologies of higher stability for models with low aggregation propensity. The free-energy profiles of the aggregation process indicate that the nucleation barrier determines the population fractions of different fibril morphologies, i.e., amyloid polymorphism is under kinetic control.     

Probing amyloid fibril formation of the NFGAIL peptide by computer simulations

Amyloid fibril formation, as observed in Alzheimer's disease and type II diabetes, is currently described by a nucleation-condensation mechanism, but the details of the process preceding the formation of the nucleus are still lacking. In this study, using an activation-relaxation technique coupled to a generic energy model, we explore the aggregation pathways of 12 chains of the hexapeptide NFGAIL. The simulations show, starting from a preformed parallel dimer and ten disordered chains, that the peptides form essentially amorphous oligomers or more rarely ordered beta-sheet structures where the peptides adopt a parallel orientation within the sheets. Comparison between the simulations indicates that a dimer is not a sufficient seed for avoiding amorphous aggregates and that there is a critical threshold in the number of connections between the chains above which exploration of amorphous aggregates is preferred.     

On the extended beta-conformation propensity of polypeptides at high temperature

At room temperature, natural polypeptides exposed to high concentrations of a strong denaturant nearly attain the circular dichroism spectra characteristic of random coils. As temperature is increased, the spectra begin to show the signature of a substantial fraction of extended chain, the structure common in beta-sheets. This structural propensity at high temperature is not altered by concentration changes over a greater than 1000-fold range, so it is not caused by aggregation. Four proteins with different folds and varying amounts of alpha-helical and beta-sheet secondary structure, in the presence or absence of denaturant, all were subject to extended chain formation upon heating. This effect arises naturally from the steric constraints associated with polypeptides and is probably counteracted, not enhanced, by hydrophobic interactions. Molecular dynamics simulations in the 298-1000 K range reveal an attractive potential of mean force in the extended chain region of the Ramachandran diagram, which broadens as the temperature is raised. We also demonstrate a direct correlation between extended structure content and the rate of aggregation kinetics. Thus pre-existing extended structure could funnel proteins into aggregates.     

Improving Internal Peptide Dynamics in the Coarse-Grained MARTINI Model: Toward Large-Scale Simulations of Amyloid- and Elastin-like Peptides

We present an extension of the coarse-grained MARTINI model for proteins and apply this extension to amyloid- and elastin-like peptides. Atomistic simulations of tetrapeptides, octapeptides, and longer peptides in solution are used as a reference to parametrize a set of pseudodihedral potentials that describe the internal flexibility of MARTINI peptides. We assess the performance of the resulting model in reproducing various structural properties computed from atomistic trajectories of peptides in water. The addition of new dihedral angle potentials improves agreement with the contact maps computed from atomistic simulations significantly. We also address the question of which parameters derived from atomistic trajectories are transferable between different lengths of peptides. The modified coarse-grained model shows reasonable transferability of parameters for the amyloid- and elastin-like peptides. In addition, the improved coarse-grained model is also applied to investigate the self-assembly of β-sheet forming peptides on the microsecond time scale. The octapeptides SNNFGAIL and (GV)(4) are used to examine peptide aggregation in different environments, in water, and at the water-octane interface. At the interface, peptide adsorption occurs rapidly, and peptides spontaneously aggregate in favor of stretched conformers resembling β-strands.     

Folding of small proteins using a single continuous potential

Extensive Monte Carlo folding simulations for four proteins of various structural classes are carried out, using a single continuous potential (united-residue force field). In all cases, collapse occurs at a very early stage, and proteins fold into their nativelike conformations at appropriate temperatures. We also observe that glassy transitions occur at low temperatures. The simulation results demonstrate that the folding mechanism is controlled not only by thermodynamic factors but also by kinetic factors: The way a protein folds into its native structure is also determined by the convergence point of early folding trajectories, which cannot be obtained by the free energy surface.     

Nanosecond field alignment of head group and water dipoles in electroporating phospholipid bilayers

To investigate the mechanism of biological cell membrane electroporation at the nanosecond and nanometer scale, we tracked pore-forming lipids and water in molecular dynamics simulations of a palmitoyloleoylphosphatidylcholine bilayer in a minimum porating electric field. Although the field-generated torque tilts the mean head group dipole a few degrees away from its equilibrium, zero-field position relative to the bilayer plane, this change in conformation does not appear to contribute directly to the development of the pore-initiating aggregation of lipid head groups and water that leads to the formation of a membrane-spanning hydrophilic pore. Field-directed rotation of the head group dipoles in the plane of the incipient pore wall, in combination with water dipole and solvation interactions at the aqueous-lipid interface, is one component in the coordinated ensemble of electroporation events.