Solvent environment conducive to protein aggregation

Trehalose Effect on the Aggregation of Model Proteins into Amyloid Fibrils

Protein aggregation into amyloid fibrils is a phenomenon that attracts attention from a wide and composite part of the scientific community. Indeed, the presence of mature fibrils is associated with several neurodegenerative diseases, and in addition these supramolecular aggregates are considered promising self-assembling nanomaterials. In this framework, investigation on the effect of cosolutes on protein propensity to aggregate into fibrils is receiving growing interest, and new insights on this aspect might represent valuable steps towards comprehension of highly complex biological processes. In this work we studied the influence exerted by the osmolyte trehalose on fibrillation of two model proteins, that is, lysozyme and insulin, investigated during concomitant variation of the solution ionic strength due to NaCl. In order to monitor both secondary structures and the overall tridimensional conformations, we have performed UV spectroscopy measurements with Congo Red, Circular Dichroism, and synchrotron Small Angle X-ray Scattering. For both proteins we describe the effect of trehalose in changing the fibrillation pattern and, as main result, we observe that ionic strength in solution is a key factor in determining trehalose efficiency in slowing down or blocking protein fibrillation. Ionic strength reveals to be a competitive element with respect to trehalose, being able to counteract its inhibiting effects toward amyloidogenesis. Reported data highlight the importance of combining studies carried out on cosolutes with valuation of other physiological parameters that may affect the aggregation process. Also, the obtained experimental results allow to hypothesize a plausible mechanism adopted by the osmolyte to preserve protein surface and prevent protein fibrillation.

Probing protein aggregation using discrete molecular dynamics

Understanding the role of biomolecular dynamics in cellular processes leading to human diseases and the ability to rationally manipulate these processes is of fundamental importance in scientific research. The last decade has witnessed significant progress in probing biophysical behavior of proteins. However, we are still limited in understanding how changes in protein dynamics and inter-protein interactions occurring in short length- and time-scales lead to aberrations in their biological function. Bridging this gap in biology probed using computer simulations marks a challenging frontier in computational biology. Here we examine hypothesis-driven simplified protein models in conjunction with discrete molecular dynamics in the study of protein aggregation, implicated in series of neurodegenerative diseases, such as Alzheimer's and Huntington's diseases. Discrete molecular dynamics simulations of simplified protein models have emerged as a powerful methodology with its ability to bridge the gap in time and length scales from protein dynamics to aggregation, and provide an indispensable tool for probing protein aggregation.

Effect of single-point sequence alterations on the aggregation propensity of a model protein

Sequences of contemporary proteins are believed to have evolved through a process that optimized their overall fitness, including their resistance to deleterious aggregation. Biotechnological processing may expose therapeutic proteins to conditions that are much more conducive to aggregation than those encountered in a cellular environment. An important task of protein engineering is to identify alternative sequences that would protect proteins when processed at high concentrations without altering their native structure associated with specific biological function. Our computational studies exploit parallel tempering simulations of coarse-grained model proteins to demonstrate that isolated amino acid residue substitutions can result in significant changes in the aggregation resistance of the protein in a crowded environment while retaining protein structure in isolation. A thermodynamic analysis of protein clusters subject to competing processes of folding and association shows that moderate mutations can produce effects similar to those caused by changes in system conditions, including temperature, concentration, and solvent composition, that affect the aggregation propensity. The range of conditions where a protein can resist aggregation can therefore be tuned by sequence alterations, although the protein generally may retain its generic ability for aggregation.

Spontaneous fibril formation by polyalanines; discontinuous molecular dynamics simulations

Fibrillary protein aggregates rich in beta-sheet structure have been implicated in the pathology of several neurodegenerative diseases. In this work, we investigate the formation of fibrils by performing discontinuous molecular dynamics simulations on systems containing 12 to 96 model Ac-KA(14)K-NH(2) peptides using our newly developed off-lattice, implicit-solvent, intermediate-resolution model, PRIME. We find that, at a low concentration, random-coil peptides assemble into alpha-helices at low temperatures. At intermediate concentrations, random-coil peptides assemble into alpha-helices at low temperatures and large beta-sheet structures at high temperatures. At high concentrations, the system forms beta-sheets over a wide range of temperatures. These assemble into fibrils above a critical temperature which decreases with concentration and exceeds the isolated peptide's folding temperature. At very high temperatures and all concentrations, the system is in a random-coil state. All of these results are in good qualitative agreement with those by Blondelle and co-workers on Ac-KA(14)K-NH(2) peptides. The fibrils observed in our simulations mimic the structural characteristics observed in experiments in terms of the number of sheets formed, the values of the intra- and intersheet separations, and the parallel peptide arrangement within each beta-sheet. Finally, we find that when the strength of the hydrophobic interaction between nonpolar side chains is high compared to the strength of hydrogen bonding, amorphous aggregates, rather than fibrillar aggregates, are formed.

Kinetics of fibril formation by polyalanine peptides

Ordered beta-sheet complexes, termed amyloid fibrils, are the underlying structural components of the intra- and extracellular fibrillar protein deposits that are associated with a variety of human diseases, including Alzheimer's, Parkinson's, and the prion diseases. In this work, we investigated the kinetics of fibril formation using our newly developed off-lattice intermediate resolution model, PRIME. The model is simple enough to allow the treatment of large multichain systems while maintaining a fairly realistic description of protein dynamics without built-in bias toward any conformation when used in conjunction with constant temperature discontinuous molecular dynamics, a fast alternative to conventional molecular dynamics. Simulations were performed on systems containing 48-96 model Ac-KA14K-NH2 peptides. We found that fibril formation for polyalanines incorporate features that are characteristic of three models, the templated assembly, nucleated polymerization, and nucleated conformational conversion models, but that none of them gave a completely satisfactory description of the simulation kinetics. Fibril formation was nucleation-dependent, occurring after a lag time that decreased with increasing peptide concentration and increased with increasing temperature. Fibril formation appeared to be a conformational conversion process in which small amorphous aggregates --> beta-sheets --> ordered nucleus --> subsequent rapid growth of a small stable fibril or protofilament. Fibril growth in our simulations involved both beta-sheet elongation, in which the fibril grew by adding individual peptides to the end of each beta-sheet, and lateral addition, in which the fibril grew by adding already formed beta-sheets to its side. The initial rate of fibril formation increased with increasing concentration and decreased with increasing temperature.

Replica exchange molecular dynamics simulations of amyloid peptide aggregation

The replica exchange molecular dynamics (REMD) approach is applied to four oligomeric peptide systems. At physiologically relevant temperature values REMD samples conformation space and aggregation transitions more efficiently than constant temperature molecular dynamics (CTMD). During the aggregation process the energetic and structural properties are essentially the same in REMD and CTMD. A condensation stage toward disordered aggregates precedes the beta-sheet formation. Two order parameters, borrowed from anisotropic fluid analysis, are used to monitor the aggregation process. The order parameters do not depend on the peptide sequence and length and therefore allow to compare the amyloidogenic propensity of different peptides

Molecular dynamics simulations of spontaneous fibril formation by random-coil peptides

Assembly of normally soluble proteins into amyloid fibrils is a cause or associated symptom of numerous human disorders, including Alzheimer's and the prion diseases. We report molecular-level simulation of spontaneous fibril formation. Systems containing 12-96 model polyalanine peptides form fibrils at temperatures greater than a critical temperature that decreases with peptide concentration and exceeds the peptide's folding temperature, consistent with experimental findings. Formation of small amorphous aggregates precedes ordered nucleus formation and subsequent rapid fibril growth through addition of beta-sheets laterally and monomeric peptides at fibril ends. The fibril's structure is similar to that observed experimentally.

Phase diagrams describing fibrillization by polyalanine peptides

Amyloid fibrils are the structural components underlying the intra- and extracellular protein deposits that are associated with a variety of human diseases, including Alzheimer's, Parkinson's, and the prion diseases. In this work, we examine the thermodynamics of fibril formation using our newly-developed off-lattice intermediate-resolution protein model, PRIME. The model is simple enough to allow the treatment of large multichain systems while maintaining a fairly realistic description of protein dynamics when used in conjunction with constant-temperature discontinuous molecular dynamics, a fast alternative to conventional molecular dynamics. We conduct equilibrium simulations on systems containing 96 Ac-KA14K-NH2 peptides over a wide range of temperatures and peptide concentrations using the replica-exchange method. Based on measured values of the heat capacity, radius of gyration, and percentage of peptides that form the various structures, a phase diagram in the temperature-concentration plane is constructed delineating the regions where each structure is stable. There are four distinct single-phase regions: alpha-helices, fibrils, nonfibrillar beta-sheets, and random coils; and four two-phase regions: random coils/nonfibrillar beta-sheets, random coils/fibrils, fibrils/nonfibrillar beta-sheets, and alpha-helices/nonfibrillar beta-sheets. The alpha-helical region is at low temperature and low concentration. The nonfibrillar beta-sheet region is at intermediate temperatures and low concentrations and expands to higher temperatures as concentration is increased. The fibril region occurs at intermediate temperatures and intermediate concentrations and expands to lower as the peptide concentration is increased. The random-coil region is at high temperatures and all concentrations; this region shifts to higher temperatures as the concentration is increased.

Dehydron: a structurally encoded signal for protein interaction

We introduce a quantifiable structural motif, called dehydron, that is shown to be central to protein-protein interactions. A dehydron is a defectively packed backbone hydrogen bond suggesting preformed monomeric structure whose Coulomb energy is highly sensitive to binding-induced water exclusion. Such preformed hydrogen bonds are effectively adhesive, since water removal from their vicinity contributes to their stability. At the structural level, a significant correlation is established between dehydrons and sites for protein complexation, with the HIV-1 capsid protein P24 complexed with antibody light-chain FAB25.3 providing the most dramatic correlation. Furthermore, the number of dehydrons in homologous similar-fold proteins from different species is shown to be a signature of proteomic complexity. The techniques are then applied to higher levels of organization: The formation of the capsid and its organization in picornaviruses correlates strongly with the distribution of dehydrons on the rim of the virus unit. Furthermore, antibody contacts and crystal contacts may be assigned to dehydrons still prevalent after the capsid has been assembled. The implications of the dehydron as an encoded signal in proteomics, bioinformatics, and inhibitor drug design are emphasized.

The role of side-chain interactions in the early steps of aggregation: Molecular dynamics simulations of an amyloid-forming peptide from the yeast prion Sup35

Understanding the early steps of aggregation at atomic detail might be crucial for the rational design of therapeutics preventing diseases associated with amyloid deposits. In this paper, aggregation of the heptapeptide GNNQQNY, from the N-terminal prion-determining domain of the yeast protein Sup35, was studied by 20 molecular dynamics runs for a total simulation time of 20 micros. The simulations generate in-register parallel packing of GNNQQNY beta-strands that is consistent with x-ray diffraction and Fourier transform infrared data. The statistically preferred aggregation pathway does not correspond to a purely downhill profile of the energy surface because of the presence of enthalpic barriers that originate from out-of-register interactions. The parallel beta-sheet arrangement is favored over the antiparallel because of side-chain contacts; in particular, stacking interactions of the tyrosine rings and hydrogen bonds between amide groups. No ordered aggregation was found in control simulations with the mutant sequence SQNGNQQRG in accord with experimental data and the strong sequence dependence of aggregation.

DSC studies of a family of natively disordered fragments from Escherichia coli thioredoxin: surface burial in intrinsic coils

The accumulating data from proteome analysis indicates that numerous proteins have segments and/or domains, involved in regulatory functions of the eukaryotic cell, which are entirely unstructured under physiological conditions, challenging the structure-function paradigm. Although many such natively unfolded proteins have been structurally analyzed by NMR spectroscopy, little is known about solvent inaccessible surfaces in premolten globules and intrinsic coils. Recent DSC studies of two protein fragments have shown a promising way to estimate the predominantly hydrophobic buried surfaces [Georgescu, R. E., García-Mira, M. M., Tasayco, M. L., and Sánchez-Ruiz, J. M. (2001) Eur. J. Biochem. 268, 1-10]. Here we report a systematic heat capacity analysis of a family of natively disordered complementary fragments of oxidized Escherichia coli thioredoxin (1-31/32-108, 1-37/38-108, 1-50/51-108, and 38-73) which provides insights into the local and nonlocal interactions contributing to the burial of predominantly hydrophobic surface in intrinsic coils.