Molecular simulation of protein aggregation

Computer simulation offers unique possibilities for investigating molecular-level phenomena difficult to probe experimentally. Drawing from a wealth of studies concerning protein folding, computational studies of protein aggregation are emerging. These studies have been successful in capturing aspects of aggregation known from experiment and are being used to refine experimental methods aimed at abating aggregation. Here we review molecular-simulation studies of protein aggregation conducted in our laboratory. Specific attention is devoted to issues with implications for biotechnology.     

Lattice simulations of aggregation funnels for protein folding.

A computer model of protein aggregation competing with productive folding is proposed. Our model adapts techniques from lattice Monte Carlo studies of protein folding to the problem of aggregation. However, rather than starting with a single string of residues, we allow independently folding strings to undergo collisions and consider their interactions in different orientations. We first present some background into the nature and significance of protein aggregation and the use of lattice Monte Carlo simulations in understanding other aspects of protein folding. The results of a series of simulation experiments involving simple versions of the model illustrate the importance of considering aggregation in simulations of protein folding and provide some preliminary understanding of the characteristics of the model. Finally, we discuss the value of the model in general and of our particular design decisions and experiments. We conclude that computer simulation techniques developed to study protein folding can provide insights into protein aggregation, and that a better understanding of aggregation may in turn provide new insights into and constraints on the more general protein folding problem.     

Solvent sensitivity of protein aggregation in Cu,Zn superoxide dismutase: a molecular dynamics simulation study

Misfolding and aggregation of Cu, Zn Superoxide Dismutase (SOD1) is often found in amyotrophic lateral sclerosis (ALS) patients. The central apo SOD1 barrel was involved in protein maturation and pathological aggregation in ALS. In this work, we employed atomistic molecular dynamics (MD) simulations to study the conformational dynamics of SOD1^barrel monomer in different concentrations of trifluoroethanol (TFE). We find concentration dependence unusual structural and dynamical features, characterized by the local unfolding of SOD1^barrel. This partially unfolded structure is characterized by the exposure of hydrophobic core, is highly dynamic in nature, and is the precursor of aggregation seen in SOD1^barrel. Our computational studies supports the hypothesis of the formation of aggregation ‘building blocks’ by means of local unfolding of apo monomer as the mechanism of SOD1 fibrillar aggregation. The non-monotonic TFE concentration dependence of protein conformational changes was explored through simulation studies. Our results suggest that altered protein conformation and dynamics within its structure may underlie the aggregation of SOD1 in ALS.     

Competition between protein folding and aggregation with molecular chaperones in crowded solutions: insight from mesoscopic simulations

The living cell is inherently crowded with proteins and macromolecules. To avoid aggregation of denatured proteins in the living cell, molecular chaperones play important roles. Here we introduce a simple model to describe crowded protein solutions with chaperone-like species based on a dynamic density functional theory. As predicted by others, our simulations show that macromolecular crowding enhances the association of proteins and chaperones. However, when the intrinsic folding rate of the protein is slow, it is possible that crowding also enhances aggregation of proteins. The results of simulation suggest that, when the concentration of the crowding agent is as high as that in the cell, the association of the protein and unbound chaperone becomes correlated with the aggregation process, and that the protein-bound chaperones efficiently destroy the potential nuclei of aggregates and thus prevent the aggregation.     

Amyloid Aggregation on Lipid Bilayers and Its Impact on Membrane Permeability

Fibrillar protein aggregates (amyloids) are involved in several common pathologies, e.g., Alzheimer's disease and type II diabetes. Accumulating evidence suggests that toxicity in amyloid-related diseases originates from the deposition of protein aggregates on the cell membrane, which results in bilayer disruption and cell leakage. The molecular mechanism of damage to the membrane, however, is still obscure. To shed light on it we have performed coarse-grained molecular dynamics simulations of fibril-forming amphipathic peptides in the presence of lipid vesicles. The simulation results show that highly amyloidogenic peptides fibrillate on the surface of the vesicle, damaging the bilayer and promoting leakage. In contrast, the ordered aggregation of peptides with low amyloidogenicity is hindered by the vesicles. Remarkably, leakage from the vesicle is caused by growing aggregates, but not mature fibrils. The simulation results provide a basis for understanding the range of aggregation behavior that is observed in experiments with fibril-forming (poly)peptides.     

Receptor-dependent compartmentalization of PPIP5K1, a kinase with a cryptic polyphosphoinositide binding domain

Amyloid-related diseases are a group of illnesses in which an abnormal accumulation of proteins into fibrillar structures is evident. Results from a wide range of studies, ranging from identification of amyloid-β dimers in the brain to biophysical characterization of the interactions between amyloidogenic peptides and lipid membranes during fibril growth shed light on the initial events which take place during amyloid aggregation. Accounts of fibril disaggregation and formation of globular aggregates due to interactions with lipids or fatty acids further demonstrate the complexity of the aggregation process and the difficulty to treat amyloid-related diseases. There is an inherent difficulty in generalizing from studies of aggregation in vitro, but the involvement of too many cellular components limits the ability to follow amyloid aggregation in a cellular (or extracellular) context. Fortunately, the development of experimental methods to generate stable globular aggregates suggests new means of studying the molecular events associated with amyloid aggregation. Furthermore, simulation studies enable deeper understanding of the experimental results and provide useful predictions that can be tested in the laboratory. Computer simulations can nowadays provide molecular or even atomistic details that are experimentally not available or very difficult to obtain. In the present review, recent developments on modelling and experiments of amyloid aggregation are reviewed, and an integrative account on how isolated interactions (as observed in vitro and in silico) combine during the course of amyloid-related diseases is presented. Finally, it is argued that an integrative approach is necessary to get a better understanding of the protein aggregation process.     

High‐throughput screening of optimal solution conditions for structural biological studies by fluorescence correlation spectroscopy

Protein aggregation is an essential molecular event in a wide variety of biological situations, and is a causal factor in several degenerative diseases. The aggregation of proteins also frequently hampers structural biological analyses, such as solution NMR studies. Therefore, precise detection and characterization of protein aggregation are of crucial importance for various research fields. In this study, we demonstrate that fluorescence correlation spectroscopy (FCS) using a single‐molecule fluorescence detection system enables the detection of otherwise invisible aggregation of proteins at higher protein concentrations, which are suitable for structural biological experiments, and consumes relatively small amounts of protein over a short measurement time. Furthermore, utilizing FCS, we established a method for high‐throughput screening of protein aggregation and optimal solution conditions for structural biological experiments.     

Thermodynamic analysis of structural transitions during GNNQQNY aggregation

Amyloid protein aggregation characterizes many neurodegenerative disorders, including Alzheimer's, Parkinson's, and Creutzfeldt‐Jakob disease. Evidence suggests that amyloid aggregates may share similar aggregation pathways, implying simulation of full‐length amyloid proteins is not necessary for understanding amyloid formation. In this study, we simulate GNNQQNY, the N‐terminal prion‐determining domain of the yeast protein Sup35 to investigate the thermodynamics of structural transitions during aggregation. Utilizing a coarse‐grained model permits equilibration on relevant time scales. Replica‐exchange molecular dynamics is used to gather simulation statistics at multiple temperatures and clear energy traps that would aversely impact results. Investigating the association of 3‐, 6‐, and 12‐chain GNNQQNY systems by calculating thermodynamic quantities and orientational order parameters, we determine the aggregation pathway by studying aggregation states of GNNQQNY. We find that the aggregation of the hydrophilic GNNQQNY sequence is mainly driven by H‐bond formation, leading to the formation of β‐sheets from the very beginning of the assembly process. Condensation (aggregation) and ordering take place simultaneously, which is underpinned by the occurrence of a single heat capacity peak. Proteins 2013; 81:1141–1155. © 2013 Wiley Periodicals, Inc.     

Chaperone-like activity and hydrophobicity of alpha-crystallin

alpha-Crystallin, a prominent member of small heat shock protein (sHsp) family and a major structural protein of the eye lens is a large polydisperse oligomer of two isoforms, alphaA- and alphaB-crystallins. Numerous studies have demonstrated that alpha-crystallin functions like a molecular chaperone in preventing the aggregation of various proteins under a wide range of stress conditions. The molecular chaperone function of alpha-crystallin is thus considered to be vital in the maintenance of lens transparency and in cataract prevention. alpha-Crystallin selectively interacts with non-native proteins thereby preventing them from aggregation and helps maintain them in a folding competent state. It has been proposed and generally accepted that alpha-crystallin suppresses the aggregation of other proteins through the interaction between hydrophobic patches on its surface and exposed hydrophobic sites of partially unfolded substrate protein. However, a quantifiable relationship between hydrophobicity and chaperone-like activity remains a matter to be concerned about. On an attentive review of studies on alpha-crystallin chaperone-like activity, particularly the studies that have direct or indirect implications to hydrophobicity and chaperone-like activity, we found several instances wherein the correlation between hydrophobicity and its chaperone-like activity is paradoxical. We thus attempted to provide an overview on the role of hydrophobicity in chaperone-like activity of alpha-crystallin, the kind of evaluation done for the first time.     

Charge-Rich Regions Modulate the Anti-Aggregation Activity of Hsp90

Protein aggregation can have dramatic effects on cellular function and plays a causative role in many human diseases. In all cells, molecular chaperones bind to aggregation-prone proteins and hinder aggregation. The ability of a protein to resist aggregation and remain soluble in aqueous solution is linked to the physical properties of the protein. Numerous physical studies demonstrate that charged atoms favor solubility. We note that many molecular chaperones possess a substantial negative charge that may allow them to impart solubility on aggregation-prone proteins. Hsp90 is one such negatively charged molecular chaperone. The charge on Hsp90 is largely concentrated in two highly acidic regions. To investigate the relationship between chaperone charge and protein solubility, we deleted these charge-rich regions and analyzed the resulting Hsp90 constructs for anti-aggregation activity. We found that deletion of both charge-rich regions dramatically impaired Hsp90 anti-aggregation activity. The anti-aggregation role of the deleted charge-rich regions could be due to net charge or sequence-specific features. To distinguish these possibilities, we attached an acid-rich region with a distinct amino acid sequence to our double-deleted Hsp90 construct. This charge rescue construct displayed effective anti-aggregation activity indicating that the net charge of Hsp90 contributes to its anti-aggregation activity.     

Structure modeling and dynamics driven mutation and phosphorylation analysis of Beta-amyloid peptides

The most common characteristics of diverse age-related neurodegenerative diseases are aggregation and accumulation of the misfolded protein in the brain. Alzheimer׳s disease (AD) is one of these protein conformational diseases. Extracellular accumulation of amyloid β (Aβ) is one the neuropathological hallmarks of Alzheimer disease. Various studies have shown that mutation in specific hydrophobic region of Aβ protein inhibit the formation of β sheet, thus aggregation of this protein is stalled. The identification of such mutation in Aβ protein can help us in elucidating the etiology of sporadic Aβ. In our study we have selected three positions: 19ILU, 21ALA and 41ILU in Aβ protein based on their hydrophobic nature and substituted them with PRO ( βSheet breaker). The effects of the substitutions were analysed using molecular dynamics simulation studies. The results validated that the mutations in the specified regions change the hydrophobicity of the protein and the βsheet formation was declined to zero per cent.     

How native proteins aggregate in solution: a dynamic Monte Carlo simulation

Aggregation of native proteins in solution is of fundamental importance with regard to both the processing and the utilization of proteins. In the present work, a dynamic Monte Carlo simulation has been performed to give a molecular insight into the way in which native proteins aggregate in solution and to explore means of suppressing aggregation, using two proteins of different compositions and conformations represented by a two-dimensional (2D) lattice model (HP model). It is shown that the native HP protein with accessible hydrophobic beads on its surface is prone to aggregation. The aggregation of this protein is intensified when the solution conditions favor the partially unfolded conformation as opposed to either the native or fully unfolded conformations. In this case, the partially unfolded proteins form the cores of aggregates, which may also encapsulate the native protein. One way to inhibit protein aggregation is to introduce polymers of appropriate hydrophobicity and chain length into the solution, such that these polymer molecules wrap around the hydrophobic regions of both the unfolded and folded proteins, thereby segregating the protein molecules. Our simulation is consistent with experimental observations reported elsewhere and provides a molecular basis for the behavior of proteins in liquid environments.     

DNA Facilitates Oligomerization and Prevents Aggregation via DNA Networks

Previous studies have shown that nucleic acids can nucleate protein aggregation in disease-related proteins, but in other cases, they can act as molecular chaperones that prevent protein aggregation, even under extreme conditions. In this study, we describe the link between these two behaviors through a combination of electron microscopy and aggregation kinetics. We find that two different proteins become soluble under harsh conditions through oligomerization with DNA. These DNA/protein oligomers form “networks,” which increase the speed of oligomerization. The cases of DNA both increasing and preventing protein aggregation are observed to stem from this enhanced oligomerization. This observation raises interesting questions about the role of nucleic acids in aggregate formation in disease states.     

分子伴侣HdeA与HdeB的作用机制

分子伴侣能够与其他蛋白质的不稳定构象相结合并使其稳定．它的功能之一是能够帮助蛋白质进行正确的折叠与组装．最新研究发现,在肠道致病菌的周质空间中存在着酸性条件下能帮助周质蛋白复性的分子伴侣HdeA和HdeB．HdeA在极端酸性的胃部环境中由二聚体迅速解离成具有伴侣活性的单体,HdeA单体能够和变性的底物蛋白结合防止它们酸诱导聚集,从而保护肠道致病菌安全到达肠道．本文对肠道致病菌的耐酸机制进行了总结,最后对 HdeA和HdeB作用机制的研究近况进行综述,最后对HdeA和HdeB以后的研究方向进行了展望．     

What the use of disease-unrelated model proteins can tell us about the molecular basis of amyloid aggregation and toxicity

Recent advances in the studies on protein aggregation have led to a reappraisal of the concepts underlying this process. The data reported in the last few years showing that protein aggregation into assemblies of amyloid type can be considered a generic property of the polypeptide chains suggest that protein aggregation in cells can be a more common phenomenon than previously believed. Furthermore, the findings that aggregates of disease-unrelated proteins display the same cytotoxicity as those formed by proteins and peptides associated with disease suggest that toxicity is a consequence of the common structure of aggregates and that, at least in most cases, it proceeds by impairing common cellular parameters such as free Ca2+ and ROS levels. The new view that aggregation of polypeptide chains and aggregate toxicity are not linked to specific amino acid sequences rises dramatically the number of sequences one can investigate to assess the molecular features underlying protein aggregation and the molecular basis of aggregate toxicity. In addition, it rises intriguing considerations on protein and cell evolution as well as on amyloid disease pathogenesis.     

Dynamical Phase Transitions Reveal Amyloid-like States on Protein Folding Landscapes

Developing an understanding of protein misfolding processes presents a crucial challenge for unlocking the mysteries of human disease. In this article, we present our observations of β-sheet-rich misfolded states on a number of protein dynamical landscapes investigated through molecular dynamics simulation and Markov state models. We employ a nonequilibrium statistical mechanical theory to identify the glassy states in a protein’s dynamics, and we discuss the nonnative, β-sheet-rich states that play a distinct role in the slowest dynamics within seven protein folding systems. We highlight the fundamental similarity between these states and the amyloid structures responsible for many neurodegenerative diseases, and we discuss potential consequences for mechanisms of protein aggregation and intermolecular amyloid formation.     

Effects of hesperidin,a flavanone glycoside interaction on the conformation,stability, and aggregation of lysozyme: multispectroscopic and molecular dynamic simulation studies?

Hesperidin (HESP), a flavanone glycoside, shows high antioxidant properties and posses ability to go through the blood–brain barrier. Therefore, it could be a potential drug molecule against aggregation based diseases such as Alzheimer’s, Parkinson’s, and systemic amyloidoses. In this work, we investigated the potential of HESP to interact with hen egg-white lysozyme (HEWL) monomer and prevent its aggregation. The HESP–HEWL binding studies were performed using a fluorescence quenching technique, molecular docking and molecular dynamics simulations. We found a strong interaction of HESP with the lysozyme monomer (K_a, ~ 5 × 10⁴ M^?1) mainly through hydrogen bonding, water bridges, and hydrophobic interactions. We showed that HESP molecule spanned the highly aggregation prone region (amino acid residues 48-101) of HEWL and prevented its fibrillar aggregation. Further, we found that HESP binding completely inhibited amorphous aggregation of the protein induced by disulfide-reducing agent tries-(2-carboxyethyl) phosphine. Conformational and stability studies as followed by various tertiary and secondary structure probes revealed that HESP binding only marginally affected the lysozyme monomer conformation and increased both stability and reversibility of the protein against thermal denaturation. Future studies should investigate detail effects of HESP on solvent dynamics, structure, and toxicity of various aggregates. The answers to these questions will not only target the basic sciences, but also have application in biomedical and biotechnological sciences.     

Increased Aggregation Is More Frequently Associated to Human Disease-Associated Mutations Than to Neutral Polymorphisms

Protein aggregation is a hallmark of over 30 human pathologies. In these diseases, the aggregation of one or a few specific proteins is often toxic, leading to cellular degeneration and/or organ disruption in addition to the loss-of-function resulting from protein misfolding. Although the pathophysiological consequences of these diseases are overt, the molecular dysregulations leading to aggregate toxicity are still unclear and appear to be diverse and multifactorial. The molecular mechanisms of protein aggregation and therefore the biophysical parameters favoring protein aggregation are better understood. Here we perform an in silico survey of the impact of human sequence variation on the aggregation propensity of human proteins. We find that disease-associated variations are statistically significantly enriched in mutations that increase the aggregation potential of human proteins when compared to neutral sequence variations. These findings suggest that protein aggregation might have a broader impact on human disease than generally assumed and that beyond loss-of-function, the aggregation of mutant proteins involved in cancer, immune disorders or inflammation could potentially further contribute to disease by additional burden on cellular protein homeostasis.     

A molecular dynamics study of acylphosphatase in aggregation-promoting conditions: the influence of trifluoroethanol/water solvent

The 98-residue protein acylphosphatase exhibits a high propensity for aggregation under certain conditions. Aggregates formed from wild-type acylphosphatase in the presence of 2,2,2-trifluoroethanol and from highly destabilized mutants are essentially identical in structure. Furthermore, it has been shown by mutational studies that different regions of the protein are important for aggregation and folding. In the present molecular dynamics study, we compare the behavior of the protein in aqueous solution and in a 25% (v/v) 2,2,2-trifluoroethanol/water environment mimicking the experimental conditions. The 2,2,2-trifluoroethanol surrounding affects the structure of the protein mostly in the regions important for aggregation, in good agreement with experimental data. This suggests that the early step of (partly) unfolding, which precedes the aggregation process, has been observed.     

A coarse-grain MD (molecular dynamic) simulation of PCL–PEG and PLA–PEG aggregation as a computational model for prediction of the drug-loading efficacy of doxorubicin

Abstract

Formulating a hydrophobic drug that is water-soluble is a pharmaceutical challenge. One way is to incorporate the drug in an amphiphilic micelle composed from an aggregation of block copolymers. Design of a good nano-micelle requires many trial-and-error attempts. In this article, we developed a computational model based on a coarse-grained molecular dynamic (MD) simulation and correlated outputs with previous studies. A good correlation shows that this model reliably simulates poly-lactic acid–poly-ethylene glycol (PLA–PEG) and poly-caprolactone (PCL)–PEG aggregation in water with and without the presence of doxorubicin.

Communicated by Ramaswamy H. Sarma