Solvent effects on the conformational transition of a model polyalanine peptide

Effect of the geometry of confining media on the stability and folding rate of α -helix proteins

Protein folding in confined media has attracted wide attention over the past 15 years due to its importance to both in vivo and in vitro applications. It is generally believed that protein stability increases by decreasing the size of the confining medium, if the medium's walls are repulsive, and that the maximum folding temperature in confinement is in a pore whose size D ₀ is only slightly larger than the smallest dimension of a protein's folded state. Until recently, the stability of proteins in pores with a size very close to that of the folded state has not received the attention it deserves. In a previous paper [L. Javidpour and M. Sahimi, J. Chem. Phys. 135, 125101 (2011)], we showed that, contrary to the current theoretical predictions, the maximum folding temperature occurs in larger pores for smaller α-helices. Moreover, in very tight pores, the free energy surface becomes rough, giving rise to a new barrier for protein folding close to the unfolded state. In contrast to unbounded domains, in small nanopores proteins with an α-helical native state that contain the β structures are entropically stabilized implying that folding rates decrease notably and that the free energy surface becomes rougher. In view of the potential significance of such results to interpretation of many sets of experimental data that could not be explained by the current theories, particularly the reported anomalously low rates of folding and the importance of entropic effects on proteins' misfolded states in highly confined environments, we address the following question in the present paper: To what extent the geometry of a confined medium affects the stability and folding rates of proteins? Using millisecond-long molecular dynamics simulations, we study the problem in three types of confining media, namely, cylindrical and slit pores and spherical cavities. Most importantly, we find that the prediction of the previous theories that the dependence of the maximum folding temperature T _f on the size D of a confined medium occurs in larger media for larger proteins is correct only in spherical geometry, whereas the opposite is true in the two other geometries that we study. Also studied is the effect of the strength of the interaction between the confined media's walls and the proteins. If the walls are only weakly or moderately attractive, a complex behavior emerges that depends on the size of the confining medium.

Dynamics of proteins aggregation. II. Dynamic scaling in confined media

In this paper, the second in a series devoted to molecular modeling of protein aggregation, a mesoscale model of proteins together with extensive discontinuous molecular dynamics simulation is used to study the phenomenon in a confined medium. The medium, as a model of a crowded cellular environment, is represented by a spherical cavity, as well as cylindrical tubes with two aspect ratios. The aggregation process leads to the formation of β sheets and eventually fibrils, whose deposition on biological tissues is believed to be a major factor contributing to many neuro-degenerative diseases, such as Alzheimer's, Parkinson's, and amyotrophic lateral sclerosis diseases. Several important properties of the aggregation process, including dynamic evolution of the total number of the aggregates, the mean aggregate size, and the number of peptides that contribute to the formation of the β sheets, have been computed. We show, similar to the unconfined media studied in Paper I [S. Zheng et al., J. Chem. Phys. 145, 134306 (2016)], that the computed properties follow dynamic scaling, characterized by power laws. The existence of such dynamic scaling in unconfined media was recently confirmed by experiments. The exponents that characterize the power-law dependence on time of the properties of the aggregation process in spherical cavities are shown to agree with those in unbounded fluids at the same protein density, while the exponents for aggregation in the cylindrical tubes exhibit sensitivity to the geometry of the system. The effects of the number of amino acids in the protein, as well as the size of the confined media, have also been studied. Similarities and differences between aggregation in confined and unconfined media are described, including the possibility of no fibril formation, if confinement is severe.

Dynamics of proteins aggregation. I. Universal scaling in unbounded media

It is well understood that in some cases proteins do not fold correctly and, depending on their environment, even properly-folded proteins change their conformation spontaneously, taking on a misfolded state that leads to protein aggregation and formation of large aggregates. An important factor that contributes to the aggregation is the interactions between the misfolded proteins. Depending on the aggregation environment, the aggregates may take on various shapes forming larger structures, such as protein plaques that are often toxic. Their deposition in tissues is a major contributing factor to many neuro-degenerative diseases, such as Alzheimer's, Parkinson's, amyotrophic lateral sclerosis, and prion. This paper represents the first part in a series devoted to molecular simulation of protein aggregation. We use the PRIME, a meso-scale model of proteins, together with extensive discontinuous molecular dynamics simulation to study the aggregation process in an unbounded fluid system, as the first step toward MD simulation of the same phenomenon in crowded cellular environments. Various properties of the aggregates have been computed, including dynamic evolution of aggregate-size distribution, mean aggregate size, number of peptides that contribute to the formation of β sheets, number of various types of hydrogen bonds formed in the system, radius of gyration of the aggregates, and the aggregates' diffusivity. We show that many of such quantities follow dynamic scaling, similar to those for aggregation of colloidal clusters. In particular, at long times the mean aggregate size S(t) grows with time as, S(t) ∼ t^z, where z is the dynamic exponent. To our knowledge, this is the first time that the qualitative similarity between aggregation of proteins and colloidal aggregates has been pointed out.

Polymorphism of fibrillar structures depending on the size of assembled Aβ17-42 peptides

The size of assembled Aβ_17-42 peptides can determine polymorphism during oligomerization and fibrillization, but the mechanism of this effect is unknown. Starting from separate random monomers, various fibrillar oligomers with distinct structural characteristics were identified using discontinuous molecular dynamics simulations based on a coarse-grained protein model. From the structures observed in the simulations, two characteristic oligomer sizes emerged, trimer and paranuclei, which generated distinct structural patterns during fibrillization. A majority of the simulations for trimers and tetramers formed non-fibrillar oligomers, which primarily progress to off-pathway oligomers. Pentamers and hexamers were significantly converted into U-shape fibrillar structures, meaning that these oligomers, called paranuclei, might be potent on-pathway intermediates in fibril formation. Fibrillar oligomers larger than hexamers generated substantial polymorphism in which hybrid structures were readily formed and homogeneous fibrillar structures appeared infrequently.

Structural Conversion of Aβ17-42 Peptides from Disordered Oligomers to U-Shape Protofilaments via Multiple Kinetic Pathways

Discovering the mechanisms by which proteins aggregate into fibrils is an essential first step in understanding the molecular level processes underlying neurodegenerative diseases such as Alzheimer's and Parkinson's. The goal of this work is to provide insights into the structural changes that characterize the kinetic pathways by which amyloid-β peptides convert from monomers to oligomers to fibrils. By applying discontinuous molecular dynamics simulations to PRIME20, a force field designed to capture the chemical and physical aspects of protein aggregation, we have been able to trace out the entire aggregation process for a system containing 8 Aβ17-42 peptides. We uncovered two fibrillization mechanisms that govern the structural conversion of Aβ17-42 peptides from disordered oligomers into protofilaments. The first mechanism is monomeric conversion templated by a U-shape oligomeric nucleus into U-shape protofilament. The second mechanism involves a long-lived and on-pathway metastable oligomer with S-shape chains, having a C-terminal turn, en route to the final U-shape protofilament. Oligomers with this C-terminal turn have been regarded in recent experiments as a major contributing element to cell toxicity in Alzheimer's disease. The internal structures of the U-shape protofilaments from our PRIME20/DMD simulation agree well with those from solid state NMR experiments. The approach presented here offers a simple molecular-level framework to describe protein aggregation in general and to visualize the kinetic evolution of a putative toxic element in Alzheimer's disease in particular.

Conformational behavior of polyalanine peptides with and without protecting groups of varying chain lengths: population of PP-II structure!

Oculopharyngeal muscular dystrophy (OPMD), a polyalanine myopathy, occurs due to expansion of homo-polyalanine stretch in normal polyadenylating binding protein nuclear 1 (PABPN1) protein from Ala10 to Ala11-17. Therefore, the conformational behavior of polyalanine peptides with n = 10-17, with and without terminal protecting groups, have been investigated with different starting geometries in water by molecular dynamics simulation studies. Alanine peptides are shown to give rise to unordered structure irrespective of starting geometry and not more than two residues at a stretch have the same/similar set of φ, ψ values. However, the final structure with terminal protecting groups look like β-strand. Unprotected poly-Ala peptides adopt twisted β-hairpin/multi hairpin like structure with increasing chain length. The number of residues having φ, ψ values in collagen region is found to be less in peptides with unprotected termini as compared to peptides with protected termini of same chain length. The results have been supported by recent synchrotron radiation circular dichroism spectroscopy of polyproline II and unordered secondary structures. Opening of the helical structure in poly-Ala peptides with protecting groups has been shown to take place from C-terminal and in peptides without protecting groups opening of helix starts from both terminals. Further, opening of helix takes more time in poly-Ala peptides without terminal protecting groups. The deviations in amide bond planarity have been discussed and compared with available experimental and computational results.

Role of side-chain interactions on the formation of α-helices in model peptides

The role played by side-chain interactions on the formation of α-helices is studied using extensive all-atom molecular dynamics simulations of polyalanine-like peptides in explicit TIP4P water. The peptide is described by the OPLS-AA force field except for the Lennard-Jones interaction between Cβ-Cβ atoms, which is modified systematically. We identify values of the Lennard-Jones parameter that promote α-helix formation. To rationalize these results, potentials of mean force (PMF) between methane-like molecules that mimic side chains in our polyalanine-like peptides are computed. These PMF exhibit a complex distance dependence where global and local minima are separated by an energy barrier. We show that α-helix propensity correlates with values of these PMF at distances corresponding to Cβ-Cβ of i-i+3 and other nearest neighbors in the α-helix. In particular, the set of Lennard-Jones parameters that promote α-helices is characterized by PMF that exhibit a global minimum at distances corresponding to i-i+3 neighbors in α-helices. Implications of these results are discussed.

Coarse-grained simulation study of sequence effects on DNA hybridization in a concentrated environment

A novel coarse-grained model is developed to elucidate thermodynamics and kinetic mechanisms of DNA self-assembly. It accounts for sequence and solvent conditions to capture key experimental results such as sequence-dependent thermal property and salt-dependent persistence length of ssDNA and dsDNA. Moreover, constant-temperature simulations on two single strands of a homogeneous sequence show two main mechanisms of hybridization: a slow slithering mechanism and a one-order faster zippering mechanism. Furthermore, large-scale simulations at a high DNA strand concentration demonstrate that DNA self-assembly is a robust and enthalpically driven process in which the formation of double helices is deciphered to occur via multiple self-assembly pathways including the strand displacement mechanism. However, sequence plays an important role in shifting the majority of one pathway over the others and controlling size distribution of self-assembled aggregates. This study yields a complex picture on the role of sequence on programmable self-assembly and demonstrates a promising simulation tool that is suitable for studies in DNA nanotechnology.

Solvent effects on kinetic mechanisms of self-assembly by peptide amphiphiles via molecular dynamics simulations

Peptide amphiphiles are known to form a variety of distinctive self-assembled nanostructures (including cylindrical nanofibers in hydrogels) dependent upon the solvent conditions. Using a novel coarse-grained model, large-scale molecular dynamics simulations are performed on a system of 800 peptide amphiphiles (sequence, palmitoyl-Val3Ala3Glu3) to elucidate kinetic mechanisms of molecular assembly as a function of the solvent conditions. The assembly process is found to occur via a multistep process with transient intermediates that ultimately leads to the stabilized nanostructures including open networks of β-sheets, cylindrical nanofibers, and elongated micelles. Different kinetic mechanisms are compared in terms of peptide secondary structures, solvent-accessible surface area, radius of gyration, relative shape anisotropy, intra/intermolecular interactions, and aggregate size dynamics to provide insightful information for the design of functional biomaterials.

Mechanism of the pH-Controlled Self-Assembly of Nanofibers from Peptide Amphiphiles

Stimuli-responsive, self-assembling nanomaterials hold a great promise to revolutionize medicine and technology. However, current discovery is slow and often serendipitous. Here we report a multiscale modeling study to elucidate the pH-controlled self-assembly of nanofibers from the peptide amphiphiles, palmitoyl-I-A₃E₄-NH₂. The coarse-grained simulations revealed the formation of random-coil based spherical micelles at strong electrostatic repulsion. However, at weak or no electrostatic repulsion, the micelles merge into a nanofiber driven by the β-sheet formation between the peptide segments. The all-atom constant pH molecular dynamics revealed a cooperative transition between random coil and β-sheet in the pH range 6-7, matching the CD data. Interestingly, although the bulk pK_a is more than one unit below the transition pH, consistent with the titration data, the highest pK_a's coincide with the transition pH, suggesting that the latter may be tuned by modulating the pK_a's of a few solvent-buried Glu side chains. Together, these data offer, to our best knowledge, the first multiresolution and quantitative view of the pH-dependent self-assembly of nanofibers. The novel protocols and insights gained are expected to advance the computer-aided design and discovery of pH-responsive nanomaterials.

Role of hydrophobicity on self-assembly by peptide amphiphiles via molecular dynamics simulations

Using a novel coarse-grained model, large-scale molecular dynamics simulations were performed to examine self-assembly of 800 peptide amphiphiles (sequence palmitoyl-V3A3E3). Under suitable physiological conditions, these molecules readily assemble into nanofibers leading to hydrogel construction as observed in experiments. Our simulations capture this spontaneous self-assembly process, including formation of secondary structure, to identify morphological transitions of distinctive nanostructures. As the hydrophobic interaction is increased, progression from open networks of secondary structures toward closed cylindrical nanostructures containing either β-sheets or random coils are observed. Moreover, temperature effects are also determined to play an important role in regulating formation of secondary structures within those nanostructures. These understandings of the molecular interactions involved and the role of environmental factors on hydrogel formation provide useful insight for development of innovative smart biomaterials for biomedical applications.

Impact of sequence on the molecular assembly of short amyloid peptides

The goal of this work is to understand how the sequence of a protein affects the likelihood that it will form an amyloid fibril and the kinetics along the fibrillization pathway. The focus is on very short fragments of amyloid proteins since these play a role in the fibrillization of the parent protein and can form fibrils themselves. Discontinuous molecular dynamics simulations using the PRIME20 force field were performed of the aggregation of 48-peptide systems containing SNQNNF (PrP (170-175)), SSTSAA (RNaseA(15-20)), MVGGVV (Aβ(35-40)), GGVVIA (Aβ(37-42)), and MVGGVVIA (Aβ(35-42)). In our simulations SNQQNF, SSTTSAA, and MVGGVV form large numbers of fibrillar structures spontaneously (as in experiment). GGVVIA forms β-sheets that do not stack into fibrils (unlike experiment). The combination sequence MVGGVVIA forms less fibrils than MVGGVV, hindered by the presence of the hydrophobic residues at the C-terminal. Analysis of the simulation kinetics and energetics reveals why MVGGVV forms fibrils and GGVVIA does not, and why adding I and A to MVGGVVIA reduces fibrillization and enhances amorphous aggregation into oligomeric structures. The latter helps explain why Aβ(1-42) assembles into more complex oligomers than Aβ(1-40), a consequence of which is that it is more strongly associated with Alzheimer's disease.

The role of electrostatics and temperature on morphological transitions of hydrogel nanostructures self-assembled by peptide amphiphiles via molecular dynamics simulations

Smart biomaterials that are self-assembled from peptide amphiphiles (PA) are known to undergo morphological transitions in response to specific physiological stimuli. The design of such customizable hydrogels is of significant interest due to their potential applications in tissue engineering, biomedical imaging, and drug delivery. Using a novel coarse-grained peptide/polymer model, which has been validated by comparison of equilibrium conformations from atomistic simulations, large-scale molecular dynamics simulations are performed to examine the spontaneous self-assembly process. Starting from initial random configurations, these simulations result in the formation of nanostructures of various sizes and shapes as a function of the electrostatics and temperature. At optimal conditions, the self-assembly mechanism for the formation of cylindrical nanofibers is deciphered involving a series of steps: (1) PA molecules quickly undergo micellization whose driving force is the hydrophobic interactions between alkyl tails; (2) neighboring peptide residues within a micelle engage in a slow ordering process that leads to the formation of β-sheets exposing the hydrophobic core; (3) spherical micelles merge together through an end-to-end mechanism to form cylindrical nanofibers that exhibit high structural fidelity to the proposed structure based on experimental data. As the temperature and electrostatics vary, PA molecules undergo alternative kinetic mechanisms, resulting in the formation of a wide spectrum of nanostructures. A phase diagram in the electrostatics-temperature plane is constructed delineating regions of morphological transitions in response to external stimuli.

Evolution of the genetic code by incorporation of amino acids that improved or changed protein function

Fifty years have passed since the genetic code was deciphered, but how the genetic code came into being has not been satisfactorily addressed. It is now widely accepted that the earliest genetic code did not encode all 20 amino acids found in the universal genetic code as some amino acids have complex biosynthetic pathways and likely were not available from the environment. Therefore, the genetic code evolved as pathways for synthesis of new amino acids became available. One hypothesis proposes that early in the evolution of the genetic code four amino acids-valine, alanine, aspartic acid, and glycine-were coded by GNC codons (N = any base) with the remaining codons being nonsense codons. The other sixteen amino acids were subsequently added to the genetic code by changing nonsense codons into sense codons for these amino acids. Improvement in protein function is presumed to be the driving force behind the evolution of the code, but how improved function was achieved by adding amino acids has not been examined. Based on an analysis of amino acid function in proteins, an evolutionary mechanism for expansion of the genetic code is described in which individual coded amino acids were replaced by new amino acids that used nonsense codons differing by one base change from the sense codons previously used. The improved or altered protein function afforded by the changes in amino acid function provided the selective advantage underlying the expansion of the genetic code. Analysis of amino acid properties and functions explains why amino acids are found in their respective positions in the genetic code.

Influence of temperature on formation of perfect tau fragment fibrils using PRIME20/DMD simulations

We investigate the fibrillization process for amyloid tau fragment peptides (VQIVYK) by applying the discontinuous molecular dynamics method to a system of 48 VQIVYK peptides modeled using a new protein model/force field, PRIME20. The aim of the article is to ascertain which factors are most important in determining whether or not a peptide system forms perfect coherent fibrillar structures. Two different directional criteria are used to determine when a hydrogen bond occurs: the original H-bond constraints and a parallel preference H-bond constraint that imparts a slight bias towards the formation of parallel versus antiparallel strands in a β-sheet. Under the original H-bond constraints, the resulting fibrillar structures contain a mixture of parallel and antiparallel pairs of strands within each β-sheet over the whole fibrillization temperature range. Under the parallel preference H-bond constraints, the β-sheets within the fibrillar structures are more likely to be parallel and indeed become perfectly parallel, consistent with X-ray crystallography, at a high temperature slightly below the fibrillization temperature. The high temperature environment encourages the formation of perfect fibril structures by providing enough time and space for peptides to rearrange during the aggregation process. There are two different kinetic mechanisms, template assembly with monomer addition at high temperature and merging/rearrangement without monomer addition at low temperature, which lead to significant differences in the final fibrillar structure. This suggests that the diverse fibril morphologies generally observed in vitro depend more on environmental conditions than has heretofore been appreciated.

Structural transitions and oligomerization along polyalanine fibril formation pathways from computer simulations

The results of a computer simulation study of the aggregation kinetics of a large system of model peptides with particular focus on the formation of intermediates are presented. Discontinuous molecular dynamic simulations were used in combination with our intermediate-resolution protein model, PRIME, to simulate the aggregation of a system of 192 polyalanine (KA(14) K) peptides at a concentration of 5 mM and a reduced temperature of T* = 0.13 starting from a random configuration and ending in the assembly of a fibrillar structure. The population of various structures, including free monomers, beta sheets, amorphous aggregates, hybrid aggregates, and fibrils, and the transitions between the structures were tracked over the course of 30 independent simulations and averaged together. The aggregation pathway for this system starts with the association of free monomers into small amorphous aggregates that then grow to moderate size by incorporating other free monomers or merging with other small amorphous aggregates. These then rearrange into either small beta sheets or hybrid aggregates formed by association between unstructured chains and beta sheets, both of which grow in size by adding free monomer chains or other small aggregates, one at a time. Fibrillar structures are formed initially either by the stacking of beta sheets, rearrangement of hybrid aggregates or association between beta sheets and hybrid aggregates. They grow by the addition of beta sheets, hybrid aggregates, and other small fibrillar structures. The rearrangement of amorphous aggregates into beta sheets is a critical and necessary step in the fibril formation pathway.

Confinement in nanopores can destabilize α-helix folding proteins and stabilize the β structures

Protein folding in confined media has attracted wide attention over the past decade due to its importance in both in vivo and in vitro applications. Currently, it is generally believed that protein stability increases by decreasing the size of the confining medium, if its interaction with the confining walls is repulsive, and that the maximum folding temperature in confinement occurs for a pore size only slightly larger than the smallest dimension of the folded state of a protein. Protein stability in pore sizes, very close to the size of the folded state, has not however received the attention that it deserves. Using detailed, 0.3-ms-long molecular dynamics simulations, we show that proteins with an α-helix native state can have an optimal folding temperature in pore sizes that do not affect the folded-state structure. In contradiction to the current theoretical explanations, we find that the maximum folding temperature occurs in larger pores for smaller α-helices. In highly confined pores the free energy surface becomes rough, and a new barrier for protein folding may appear close to the unfolded state. In addition, in small nanopores the protein states that contain the β structures are entropically stabilized, in contrast to the bulk. As a consequence, folding rates decrease notably and the free energy surface becomes rougher. The results shed light on many recent experimental observations that cannot be explained by the current theories, and demonstrate the importance of entropic effects on proteins' misfolded states in highly confined environments. They also support the concept of passive effect of chaperonin GroEL on protein folding by preventing it from aggregation in crowded environment of biological cells, and provide deeper clues to the α → β conformational transition, believed to contribute to Alzheimer's and Parkinson's diseases. The strategy of protein and enzyme stabilization in confined media may also have to be revisited in the case of tight confinement. For in silico studies of protein folding in confined media, use of non-Go potentials may be more appropriate.

Gay-Berne and electrostatic multipole based coarse-grain potential in implicit solvent

A general, transferable coarse-grain (CG) framework based on the Gay-Berne potential and electrostatic point multipole expansion is presented for polypeptide simulations. The solvent effect is described by the Generalized Kirkwood theory. The CG model is calibrated using the results of all-atom simulations of model compounds in solution. Instead of matching the overall effective forces produced by atomic models, the fundamental intermolecular forces such as electrostatic, repulsion-dispersion, and solvation are represented explicitly at a CG level. We demonstrate that the CG alanine dipeptide model is able to reproduce quantitatively the conformational energy of all-atom force fields in both gas and solution phases, including the electrostatic and solvation components. Replica exchange molecular dynamics and microsecond dynamic simulations of polyalanine of 5 and 12 residues reveal that the CG polyalanines fold into "alpha helix" and "beta sheet" structures. The 5-residue polyalanine displays a substantial increase in the "beta strand" fraction relative to the 12-residue polyalanine. The detailed conformational distribution is compared with those reported from recent all-atom simulations and experiments. The results suggest that the new coarse-graining approach presented in this study has the potential to offer both accuracy and efficiency for biomolecular modeling.

Fibrillization propensity for short designed hexapeptides predicted by computer simulation

Assembly of normally soluble proteins into ordered aggregates, known as amyloid fibrils, is a cause or associated symptom of numerous human disorders, including Alzheimer's and the prion diseases. Here, we test the ability of discontinuous molecular dynamics (DMD) simulations based on PRIME20, a new intermediate-resolution protein force field, to predict which designed hexapeptide sequences will form fibrils, which will not, and how this depends on temperature and concentration. Simulations were performed on 48-peptide systems containing STVIIE, STVIFE, STVIVE, STAIIE, STVIAE, STVIGE, and STVIEE starting from random-coil configurations. By the end of the simulations, STVIIE and STVIFE (which form fibrils in vitro) form fibrils over a range of temperatures, STVIEE (which does not form fibrils in vitro) does not form fibrils, and STVIVE, STAIIE, STVIAE, and STVIGE (which do not form fibrils in vitro) form fibrils at lower temperatures but stop forming fibrils at higher temperatures. At the highest temperatures simulated, the results on the fibrillization propensity of the seven short de novo designed peptides all agree with the experiments of López de la Paz and Serrano. Our results suggest that the fibrillization temperature (temperature above which fibrils cease to form) is a measure of fibril stability and that by rank ordering the fibrillization temperatures of various sequences, PRIME20/DMD simulations could be used to ascertain their relative fibrillization propensities. A phase diagram showing regions in the temperature-concentration plane where fibrils are formed in our simulations is presented.

Hydrophobic interactions in the formation of secondary structures in small peptides

Effects of the attractive and repulsive parts of hydrophobic interactions on α helices and β sheets in small peptides are investigated using a simple atomic potential. Typically, a physical spatial range of attraction tends to favor β sheets, but α helices would be favored if the attractive range were more extended. We also found that desolvation barriers favor β sheets in collapsed conformations of polyalanine, polyvaline, polyleucine, and three fragments of amyloid peptides tested in this study. Our results provide insight into the multifaceted role of hydrophobicity in secondary structure formation, including the α to β transitions in certain amyloid peptides.

Effects of surface hydrophobicity on the conformational changes of polypeptides of different length

We studied the effects of surface hydrophobicity on the conformational changes of different length polypeptides by calculating the free energy difference between peptide structures using the bias-potential Monte Carlo technique and the probability ratio method. It was found that the hydrophobic surface plays an important role in the stability of secondary structures of the polypeptides with hydrophobic side chains. For short GAAAAG peptides, the hydrophobic surface destabilizes the α helix but stabilizes the β hairpin in the entire temperature region considered in our study. Interestingly, when the surface hydrophobic strength ε(hpsf)≥ε(hp), the most stable structure in the low temperature region changes from α helix to β hairpin, and the corresponding phase transition temperature increases slightly. For longer GAAAAAAAAAAG peptides, the effects of the relatively weak hydrophobic surface (ε(hpsf) < ε(hp)) on α-helical structures may be neglected, while the relatively strongly hydrophobic surface (ε(hpsf)≥ε(hp)) leads to the obvious partial helicity loss. In contrast, the stability of β structures can be enhanced significantly by the hydrophobic surface, especially by the strongly hydrophobic surface, at low and intermediate temperatures. At high temperatures, in addition to thermal fluctuations, the strongly hydrophobic surface (ε(hpsf)>ε(hp)) may further disturb the formation of both α-helical and β structures. Moreover, the phase transition temperature between α-helical structures and random coils significantly decreases due to the helicity loss when ε(hpsf)>ε(hp). Our findings provide a basic and quantitative picture for understanding the effects of a hydrophobic surface on the conformational changes of the polypeptides with hydrophobic side chains. From an application viewpoint, the present study is helpful in developing alternative strategies of producing high-quality biological fibrillar materials and functional nanoscale devices by the self-assembly of the polypeptides on hydrophobic surfaces.

Computer simulation study of amyloid fibril formation by palindromic sequences in prion peptides

We simulate the aggregation of large systems containing palindromic peptides from the Syrian hamster prion protein SHaPrP 113-120 (AGAAAAGA) and the mouse prion protein MoPrP 111-120 (VAGAAAAGAV) and eight sequence variations: GAAAAAAG, (AG)(4) , A8, GAAAGAAA, A10, V10, GAVAAAAVAG, and VAVAAAAVAV The first two peptides are thought to act as the Velcro that holds the parent prion proteins together in amyloid structures and can form fibrils themselves. Kinetic events along the fibrillization pathway influence the types of structures that occur and variations in the sequence affect aggregation kinetics and fibrillar structure. Discontinuous molecular dynamics simulations using the PRIME20 force field are performed on systems containing 48 peptides starting from a random coil configuration. Depending on the sequence, fibrillar structures form spontaneously over a range of temperatures, below which amorphous aggregates form and above which no aggregation occurs. AGAAAAGA forms well organized fibrillar structures whereas VAGAAAAGAV forms less well organized structures that are partially fibrillar and partially amorphous. The degree of order in the fibrillar structure stems in part from the types of kinetic events leading up to its formation, with AGAAAAGA forming less amorphous structures early in the simulation than VAGAAAAGAV. The ability to form fibrils increases as the chain length and the length of the stretch of hydrophobic residues increase. However as the hydrophobicity of the sequence increases, the ability to form well-ordered structures decreases. Thus, longer hydrophobic sequences form slightly disordered aggregates that are partially fibrillar and partially amorphous. Subtle changes in sequence result in slightly different fibril structures.

Molecular level studies on binding modes of labeling molecules with polyalanine peptides

In this work, the binding modes of typical labeling molecules (thioflavin T (ThT), Congo red (CR) and copper(II) phthalocyanine tetrasulfonic acid tetrasodium salt (PcCu(SO(3)Na)(4))) on pentaalanine, which is a model peptide segment of amyloid peptides, have been resolved at the molecular level by using scanning tunneling microscopy (STM). In the STM images, ThT molecules are predominantly adsorbed parallel to the peptide strands and two binding modes could be identified. It was found that ThT molecules are preferentially binding on top of the peptide strand, and the mode of intercalated between neighboring peptides also exists. The parallel binding mode of CR molecules can be observed with pentaalanine peptides. Besides the binding modes of labeling molecules, the CR and PcCu(SO(3)Na)(4) display different adsorption affinity with the pentaalanine peptides. The results could be beneficial for obtaining molecular level insight of the interactions between labeling molecules and peptides.

Self-assembly of polypeptides into left-handedly twisted fibril-like structures

In this paper, we investigated the spontaneous formation of aggregation structures of amyloid-forming peptide (GGVVIA) using a coarse-grained model and Monte Carlo simulations. The effects of concentration and temperature on the formation of different aggregation structures were studied. Three types of aggregation structures, single-layer beta sheet, amorphous beta-sheet aggregate, and fibril-like structures, were observed in our simulations. The fibril-like structures obtained in simulations have a common cross-beta spine structure in which beta sheets twist in a left-handed fashion. The averaged twisting angle of the beta sheet in the fibril-like structures is 12 degrees +/-2 degrees. Moreover, it was found that the peptides in the same beta sheets prefer to arrange in a parallel way, which is consistent with the corresponding GGVVIA crystalline structure. On the other hand, it was found that there is a rich family of beta-sheet stacking patterns in the fibril-like structures suggesting that the fibril structures are more complex than the corresponding crystalline structure and there exist many local free-energy minima rather than a distinct global minimum.

Structure and stability of helices in square-well homopolymers

Recently, it has been demonstrated [Magee, Phys. Rev. Lett. 96, 207802 (2006)] that isolated square-well homopolymers can spontaneously break chiral symmetry and "freeze" into helical structures at sufficiently low temperatures. This behavior is interesting because the square-well homopolymer is itself achiral. In this work, we use event-driven molecular dynamics combined with an optimized parallel tempering scheme to study this polymer model over a wide range of parameters. We examine the conditions where the helix structure is stable and determine how the interaction parameters of the polymer govern the details of the helix structure. The width of the square well (proportional to lambda) is found to control the radius of the helix, which decreases with increasing well width until the polymer forms a coiled sphere for sufficiently large wells. The helices are found to be stable for only a "window" of molecular weights. If the polymer is too short, the helix will not form. If the polymer is too long, the helix is no longer the minimum energy structure, and other folded structures will form. The size of this window is governed by the chain stiffness, which in this model is a function of the ratio of the monomer size to the bond length. Outside this window, the polymer still freezes into a locked structure at low temperature; however, unless the chain is sufficiently stiff, this structure will not be unique and is similar to a glassy state.

Molecular simulation of protein dynamics in nanopores. II. Diffusion

A novel combination of discontinuous molecular dynamics and the Langevin equation, together with an intermediate-resolution model of proteins, is used to carry out long (several microsecond) simulations in order to study transport of proteins in nanopores. We simulated single-domain proteins with the alpha-helical native structure. Both attractive and repulsive interaction potentials between the proteins and the pores' walls are considered. The diffusivity D of the proteins is computed not only under the bulk conditions but also as a function of their "length" (the number of the amino-acid groups), temperature T, pore size, and interaction potentials with the walls. Compared with the experimental data, the computed diffusivities under the bulk conditions are of the correct order of magnitude. The diffusivities both in the bulk and in the pores follow a power law in the length [script-l] of the proteins and are larger in pores with repulsive walls. D(+)/D(-), the ratio of the diffusivities in pores with attractive and repulsive walls, exhibits two local maxima in its dependence on the pore size h, which are attributed to the pore sizes and protein configurations that induce long-lasting simultaneous interactions with both walls of the pores. Far from the folding temperature T(f), D increases about linearly with T, but due to the thermal fluctuations and their effect on the proteins' structure near T(f), the dependence of D on T in this region is nonlinear. We propose a novel and general "phase diagram," consisting of four regions, that describes qualitatively the effect of h, T, and interaction potentials with the walls on the diffusivity D of a protein.

Activation of immobilized lipase in non-aqueous systems by hydrophobic poly-DL-tryptophan tethers

Many industrially important reactions use immobilized enzymes in non-aqueous, organic systems, particularly for the production of chiral compounds such as pharmaceutical precursors. The addition of a spacer molecule ("tether") between a supporting surface and enzyme often substantially improves the activity and stability of enzymes in aqueous solution. Most "long" linkers (e.g., polyethylene oxide derivatives) are relatively hydrophilic, improving the solubility of the linker-enzyme conjugate in polar environments, but this provides little benefit in non-polar environments such as organic solvents. We present a novel method for the covalent immobilization of enzymes on solid surfaces using a long, hydrophobic polytryptophan tether. Candida antarctica lipase B (CALB) was covalently immobilized on non-porous, functionalized 1-microm silica microspheres, with and without an intervening hydrophobic poly-DL-tryptophan tether (n approximately 78). The polytryptophan-tethered enzyme exhibited 35 times greater esterification of n-propanol with lauric acid in the organic phase and five times the hydrolytic activity against p-nitrophenol palmitate, compared to the activity of the same enzyme immobilized without tethers. In addition, the hydrophobic tethers caused the silica microspheres to disperse more readily in the organic phase, while the surface-immobilized control treatment was less lipophilic and quickly settled out of the organic phase when the suspensions were not vigorously mixed.

Folding of gas-phase polyalanines in a static electric field: alignment, deformations, and polarization effects

Monte Carlo simulations of the temperature-induced unfolding of small gas-phase polyalanines in a static, homogeneous electric field are reported, based on the AMBER ff96 force field. The peptides exhibit a structural transition from the native alpha-helix state to entropically favored beta-sheet conformations, before eventually turning to extended coil at higher temperatures. Upon switching the electric field, the molecules undergo preferential alignment of their dipole moment vector toward the field axis and a shift of the alpha-beta transition to higher temperatures. At higher field strengths (>10(8) V/m) the molecules stretch and the alpha-beta and beta-coil transitions merge. A simple three-state model is shown to account for the observed behavior. Under even higher fields, density functional theory calculations and a polarizable force field both show that electronic rearrangements tend to further increase the dipole moment, polarization effects being approximately half in magnitude with respect to stretching effect. Finally a tentative (temperature, field-strength) phase diagram is sketched.

Molecular simulation of protein dynamics in nanopores. I. Stability and folding

Discontinuous molecular dynamics simulations, together with the protein intermediate resolution model, an intermediate-resolution model of proteins, are used to carry out several microsecond-long simulations and study folding transition and stability of alpha-de novo-designed proteins in slit nanopores. Both attractive and repulsive interaction potentials between the proteins and the pore walls are considered. Near the folding temperature T(f) and in the presence of the attractive potential, the proteins undergo a repeating sequence of folding/partially folding/unfolding transitions, with T(f) decreasing with decreasing pore sizes. The unfolded states may even be completely adsorbed on the pore's walls with a negative potential energy. In such pores the energetic effects dominate the entropic effects. As a result, the unfolded state is stabilized, with a folding temperature T(f) which is lower than its value in the bulk and that, compared with the bulk, the folding rate decreases. The opposite is true in the presence of a repulsive interaction potential between the proteins and the walls. Moreover, for short proteins in very tight pores with attractive walls, there exists an unfolded state with only one alpha-helical hydrogen bond and an energy nearly equal to that of the folded state. The proteins have, however, high entropies, implying that they cannot fold onto their native structure, whereas in the presence of repulsive walls the proteins do attain their native structure. There is a pronounced asymmetry between the two termini of the protein with respect to their interaction with the pore walls. The effect of a variety of factors, including the pore size and the proteins' length, as well as the temperature, is studied in detail.

Modulation of self-association and subsequent fibril formation in an alanine-rich helical polypeptide

Thermal unfolding, reversible self-association, and irreversible aggregation were investigated for an alanine-rich helical polypeptide, 17-H-6, with sequence [AAAQEAAAAQAAAQAEAAQAAQ] 6. Dynamic light scattering, transmission electron microscopy, and thermal unfolding measurements indicate that 17-H-6 spontaneously and reversibly self-associates at acidic pH and low temperature. The resulting multimers have a compact, globular morphology with an average hydrodynamic radius approximately 10-20 nm and reversibly dissociate to monomers upon an increase to pH 7.4. Both free monomer and 17-H-6 chains within the multimers are alpha-helical and folded at low temperature. Reversible unfolding of the monomer occurs upon heating of solutions at pH 7.4. At pH 2.3, heating first causes incomplete dissociation and unfolding of the constituent chains. Further incubation at elevated temperature induces additional structural and morphological changes and results in fibrils with a beta-sheet 2 degrees structure and a characteristic diameter of 5-10 nm (7 nm mean). The ability to modulate association and aggregation suggests opportunities for this class of polypeptides in nanotechnology and biomedical applications.