Interconversion between Serum Amyloid A Native and Fibril Conformations
- Fatih Yasar*Fatih Yasar*Email: [email protected]Department of Chemistry & Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United StatesMore by Fatih Yasar
- ,
- Miranda S. Sheridan*Miranda S. Sheridan*Email: [email protected]Department of Chemistry & Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United StatesMore by Miranda S. Sheridan
- , and
- Ulrich H. E. Hansmann*Ulrich H. E. Hansmann*Email: [email protected]Department of Chemistry & Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United StatesMore by Ulrich H. E. Hansmann
Abstract

Overexpression of serum amyloid A (SAA) can lead to a form of amyloidosis where the fibrils are made of SAA fragments, most often SAA1–76. Using Replica Exchange with Tunneling, we study the conversion of a SAA1–76 chain between the folded conformation and a fibril conformation. We find that the basins in the free energy landscape corresponding to the two motifs are separated by barriers of only about 2–3 kBT. Crucial for the assembly into the fibril structure is the salt bridge 26E–34K that provides a scaffold for forming the fibril conformation.
This publication is licensed under
License Summary*
You are free to share (copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share (copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share (copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share (copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share (copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
1. Introduction
Figure 1

2. Materials and Methods
2.1. Replica Exchange with Tunneling


2.2. Simulation Setup
2.3. Analysis Tools and Protocols
3. Results and Discussion
Figure 2

Figure 2. (a) A typical example of a replica walking through λ space starting from a replica where the physical model is initially biased toward the folded SAA structure. While the system walks between a replica with bias toward the folded structure (upper half) and a replica with bias toward the fibril structure (lower half), its configuration changes accordingly. This can be seen in (b), where we show the corresponding time evolution of the RMSD to the native structure (in magenta) and the fibril structure (in green).
Figure 3

Figure 3. Free energy landscape as obtained from RET simulations, with data taken at λ = 0, i.e., where the physical models are not biased by any Go-term. Energies are listed in units of kT. The prospective transition pathway is drawn in black, and the five regions crossed by this path are marked in capital letters.
Figure 4

Figure 4. Characteristic conformations of of SAA1–76 as seen in each of the five regions (labeled A, B, C, D, and E) identified on the proposed transition pathway. The N-terminus of the chains is colored in blue. Unlike for the transition regions B, C, and D, these conformations superimposed on the respective reference structures for region A (dominated by folded-conformations) and region E (where fibril-like conformations are dominant).
Figure 5

Figure 5. (a) Number of contacts (normalized to one) nNS that are shared with the folded structure as measured in each of the five regions A to E of the transition pathway. The subset of long-range contacts nLR, again normalized to one, is drawn separately. Shown are also the number nFS of contacts also found in the fibril reference structure. In (b), we show the relative frequency with which one of the three characteristic helices of the folded structure, or the two main β-strands of the fibril structure, is observed.
Figure 6

Figure 6. Residue–residue map of the average minimal distance between heavy atoms in a pair of residues, shown for each of the five regions A, B, C, D, and E. Unlike in the transition regions B, C, and D, folded conformations dominate in region A, and fibril-like ones dominate in region E. Residue pairs whose average contact distance is more than 10 Å are excluded. Numbers on the X and Y axis mark residues, and the average distance in Å is given by the color coding shown in the middle row.
4. Conclusions
Acknowledgments
The simulations in this work were done using the SCHOONER cluster of the University of Oklahoma, XSEDE resources allocated under Grant MCB160005, and TACC resources allocated under Grant MCB20016 (both National Science Foundation). We acknowledge financial support from the National Institutes of Health under Grants GM120634 and GM120578. F.Y. also wants to thank the Scientific and Technological Research Council of Turkey (TUBITAK) under the BIDEB programs and the Department of Chemistry and Biochemistry for kind hospitality during the his sabbatical stay at the University of Oklahoma.
References
This article references 45 other publications.
- 1Chiti, F.; Dobson, C. M. Protein misfolding, amyloid formation, and human disease: A summary of progress over the last decade. Annu. Rev. Biochem. 2017, 86, 27– 68, DOI: 10.1146/annurev-biochem-061516-045115[Crossref], [PubMed], [CAS], Google Scholar1https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXns1Wjsr0%253D&md5=e1802f20bcaacdd5097f5ba13a42e4eaProtein Misfolding, Amyloid Formation, and Human Disease: A Summary of Progress Over the Last DecadeChiti, Fabrizio; Dobson, Christopher M.Annual Review of Biochemistry (2017), 86 (), 27-68CODEN: ARBOAW; ISSN:0066-4154. (Annual Reviews)Peptides and proteins have been found to possess an inherent tendency to convert from their native functional states into intractable amyloid aggregates. This phenomenon is assocd. with a range of increasingly common human disorders, including Alzheimer and Parkinson diseases, type II diabetes, and a no. of systemic amyloidoses. In this review, we describe this field of science with particular ref. to the advances that have been made over the last decade in our understanding of its fundamental nature and consequences. We list the proteins that are known to be deposited as amyloid or other types of aggregates in human tissues and the disorders with which they are assocd., as well as the proteins that exploit the amyloid motif to play specific functional roles in humans. In addn., we summarize the genetic factors that have provided insight into the mechanisms of disease onset. We describe recent advances in our knowledge of the structures of amyloid fibrils and their oligomeric precursors and of the mechanisms by which they are formed and proliferate to generate cellular dysfunction. We show evidence that a complex proteostasis network actively combats protein aggregation and that such an efficient system can fail in some circumstances and give rise to disease. Finally, we anticipate the development of novel therapeutic strategies with which to prevent or treat these highly debilitating and currently incurable conditions.
- 2Sipe, J. D.; Benson, M. D.; Buxbaum, J. N.; Ikeda, S. I.; Merlini, G.; Saraiva, M. J.; Westermark, P. Amyloid fibril proteins and amyloidosis: chemical identification and clinical classification International Society of Amyloidosis 2016 Nomenclature Guidelines. Amyloid 2016, 23, 209– 213, DOI: 10.1080/13506129.2016.1257986[Crossref], [PubMed], [CAS], Google Scholar2https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhvF2msbzI&md5=ffcab18f122fb48c87aba322d03ed0f0Amyloid fibril proteins and amyloidosis: chemical identification and clinical classification International Society of Amyloidosis 2016 Nomenclature GuidelinesSipe, Jean D.; Benson, Merrill D.; Buxbaum, Joel N.; Ikeda, Shu-ichi; Merlini, Giampaolo; Saraiva, Maria J. M.; Westermark, PerAmyloid (2016), 23 (4), 209-213CODEN: AIJIET; ISSN:1350-6129. (Taylor & Francis Ltd.)The Nomenclature Committee of the International Society of Amyloidosis (ISA) met during the XVth Symposium of the Society, 3 July-7 July 2016, Uppsala, Sweden, to assess and formulate recommendations for nomenclature for amyloid fibril proteins and the clin. classification of the amyloidoses. An amyloid fibril must exhibit affinity for Congo red and with green, yellow or orange birefringence when the Congo red-stained deposits are viewed with polarized light. While congophilia and birefringence remain the gold std. for demonstration of amyloid deposits, new staining and imaging techniques are proving useful. To be included in the nomenclature list, in addn. to congophilia and birefringence, the chem. identity of the protein must be unambiguously characterized by protein sequence anal. when possible. In general, it is insufficient to identify a mutation in the gene of a candidate amyloid protein without confirming the variant changes in the amyloid fibril protein. Each distinct form of amyloidosis is uniquely characterized by the chem. identity of the amyloid fibril protein that deposits in the extracellular spaces of tissues and organs and gives rise to the disease syndrome. The fibril proteins are designated as protein A followed by a suffix that is an abbreviation of the parent or precursor protein name. To date, there are 36 known extracellular fibril proteins in humans, 2 of which are iatrogenic in nature and 9 of which have also been identified in animals. Two newly recognized fibril proteins, AApoCII derived from apolipoprotein CII and AApoCIII derived from apolipoprotein CIII, have been added. AApoCII amyloidosis and AApoCIII amyloidosis are hereditary systemic amyloidoses. Intracellular protein inclusions displaying some of the properties of amyloid, "intracellular amyloid" have been reported. Two proteins which were previously characterized as intracellular inclusions, tau and α-synuclein, are now recognized to form extracellular deposits upon cell death and thus have been included in Table 1 as ATau and AαSyn.
- 3Obici, L.; Merlini, G. AA amyloidosis: Basic knowledge, unmet needs and future treatments. Swiss Med. Wkly. 2012, 142, w13580, DOI: 10.4414/smw.2012.13580
- 4Röcken, C.; Shakespeare, A. Pathology, diagnosis and pathogenesis of AA amyloidosis. Virchows Arch 2002, 440, 111– 122, DOI: 10.1007/s00428-001-0582-9[Crossref], [PubMed], [CAS], Google Scholar4https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BD383is12ktw%253D%253D&md5=578f51b866526454b68445925df76deePathology, diagnosis and pathogenesis of AA amyloidosisRocken Christoph; Shakespeare AnnVirchows Archiv : an international journal of pathology (2002), 440 (2), 111-122 ISSN:0945-6317.Amyloid is defined as a proteinaceous tissue deposit that shows a typical green birefringence in polarised light after staining with Congo red, the presence of non-branching linear fibrils of indefinite length with an approximate diameter of 10-12 nm and a distinct X-ray diffraction pattern consistent with Pauling's model of a cross-beta fibril. Approximately 45% of generalised amyloidoses are secondary or reactive (AA) amyloidosis. Among the causes of AA amyloidosis are rheumatic diseases, idiopathic diseases, inherited diseases, infectious diseases and malignant tumours. Recent decades have provided significant advances in our understanding of the pathology and pathogenesis of AA amyloidosis. Its pathogenesis is multifactorial involving many variables such as primary structure of the precursor protein, acute phase response, the presence of non-fibril proteins (e.g. amyloid P component, apolipoprotein E, glycosaminoglycans, proteoglycans and basement membrane proteins), receptors, lipid metabolism and proteases. Study of the pathogenesis of AA amyloidosis has provided many insights into the nature of conformational diseases, which may help in the understanding of other members of this particularly heterogeneous group of diseases, such as Alzheimer's disease and transmissible spongiform encephalopathies.
- 5Real de Asua, D.; Galvan, J. M.; Filigghedu, M. T.; Trujillo, D.; Costa, R.; Cadinanos, J. Systemic AA amyloidosis: epidemiology, diagnosis, and management. Clin. Epidemiol. 2014, 6, 369, DOI: 10.2147/CLEP.S39981[Crossref], [PubMed], [CAS], Google Scholar5https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BC2M3ms1eisg%253D%253D&md5=fcef6da5b8c6d8f0eda36e6e9fab8d28Systemic AA amyloidosis: epidemiology, diagnosis, and managementReal de Asua Diego; Costa Ramon; Galvan Jose Maria; Filigheddu Maria Teresa; Trujillo Davinia; Cadinanos JulenClinical epidemiology (2014), 6 (), 369-77 ISSN:1179-1349.The term "amyloidosis" encompasses the heterogeneous group of diseases caused by the extracellular deposition of autologous fibrillar proteins. The global incidence of amyloidosis is estimated at five to nine cases per million patient-years. While amyloid light-chain (AL) amyloidosis is more frequent in developed countries, amyloid A (AA) amyloidosis is more common in some European regions and in developing countries. The spectrum of AA amyloidosis has changed in recent decades owing to: an increase in the median age at diagnosis; a percent increase in the frequency of primary AL amyloidosis with respect to the AA type; and a substantial change in the epidemiology of the underlying diseases. Diagnosis of amyloidosis is based on clinical organ involvement and histological evidence of amyloid deposits. Among the many tinctorial characteristics of amyloid deposits, avidity for Congo red and metachromatic birefringence under unidirectional polarized light remain the gold standard. Once the initial diagnosis has been made, the amyloid subtype must be identified and systemic organ involvement evaluated. In this sense, the (123)I-labeled serum amyloid P component scintigraphy is a safe and noninvasive technique that has revolutionized the diagnosis and monitoring of treatment in systemic amyloidosis. It can successfully identify anatomical patterns of amyloid deposition throughout the body and enables not only an initial estimation of prognosis, but also the monitoring of the course of the disease and the response to treatment. Given the etiologic diversity of AA amyloidosis, common therapeutic strategies are scarce. All treatment options should be based upon a greater control of the underlying disease, adequate organ support, and treatment of symptoms. Nevertheless, novel therapeutic strategies targeting the formation of amyloid fibrils and amyloid deposition may generate new expectations for patients with AA amyloidosis.
- 6Nguyen, P. H.; Ramamoorthy, A.; Sahoo, B. R.; Zheng, J.; Faller, P.; Straub, J. E.; Dominguez, L.; Shea, J.-E.; Dokholyan, N. V.; De Simone, A. Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer’s Disease, Parkinson’s Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis. Chem. Rev. 2021, 121, 2545– 2647, DOI: 10.1021/acs.chemrev.0c01122[ACS Full Text
], [CAS], Google Scholar
6https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXjtVGmtb0%253D&md5=242576b480f2ce99fa61e6a2bb2dfcaeAmyloid Oligomers: A joint experimental/computational perspective on Alzheimer's disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral SclerosisNguyen, Phuong H.; Ramamoorthy, Ayyalusamy; Sahoo, Bikash R.; Zheng, Jie; Faller, Peter; Straub, John E.; Dominguez, Laura; Shea, Joan-Emma; Dokholyan, Nikolay V.; De Simone, Alfonso; Ma, Buyong; Nussinov, Ruth; Najafi, Saeed; Ngo, Son Tung; Loquet, Antoine; Chiricotto, Mara; Ganguly, Pritam; McCarty, James; Li, Mai Suan; Hall, Carol; Wang, Yiming; Miller, Yifat; Melchionna, Simone; Habenstein, Birgit; Timr, Stepan; Chen, Jiaxing; Hnath, Brianna; Strodel, Birgit; Kayed, Rakez; Lesne, Sylvain; Wei, Guanghong; Sterpone, Fabio; Doig, Andrew J.; Derreumaux, PhilippeChemical Reviews (Washington, DC, United States) (2021), 121 (4), 2545-2647CODEN: CHREAY; ISSN:0009-2665. (American Chemical Society)A review. Protein misfolding and aggregation is obsd. in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by exptl. and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacol. expts. tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, resp., for many years. - 7Cai, L.; De Beer, M. C.; De Beer, F. C.; Van Der Westhuyzen, D. R. Serum amyloid a is a ligand for scavenger receptor class B type I and inhibits high density lipoprotein binding and selective lipid uptake. J. Biol. Chem. 2005, 280, 2954– 2961, DOI: 10.1074/jbc.M411555200[Crossref], [PubMed], [CAS], Google Scholar7https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXmt1yqtA%253D%253D&md5=296babba7b8568a797fb9fe485a5843bSerum Amyloid A Is a Ligand for Scavenger Receptor Class B Type I and Inhibits High Density Lipoprotein Binding and Selective Lipid UptakeCai, Lei; de Beer, Maria C.; de Beer, Frederick C.; van der Westhuyzen, Deneys R.Journal of Biological Chemistry (2005), 280 (4), 2954-2961CODEN: JBCHA3; ISSN:0021-9258. (American Society for Biochemistry and Molecular Biology)Serum amyloid A is an acute phase protein that is carried in the plasma largely as an apolipoprotein of high d. lipoprotein (HDL). In this study we investigated whether SAA is a ligand for the HDL receptor, scavenger receptor class B type I (SR-BI), and how SAA may influence SR-BI-mediated HDL binding and selective cholesteryl ester uptake. Studies using Chinese hamster ovary cells expressing SR-BI showed that 125I-labeled SAA, both in lipid-free form and in reconstituted HDL particles, functions as a high affinity ligand for SR-BI. SAA also bound with high affinity to the hepatocyte cell line, HepG2. Alexa-labeled SAA was shown by fluorescence confocal microscopy to be internalized by cells in a SR-BI-dependent manner. To assess how SAA assocn. with HDL influences HDL interaction with SR-BI, SAA-contg. HDL was isolated from mice overexpressing SAA through adenoviral gene transfer. SAA presence on HDL had little effect on HDL binding to SR-BI but decreased (30-50%) selective cholesteryl ester uptake. Lipid-free SAA, unlike lipid-free apoA-I, was an effective inhibitor of both SR-BI-dependent binding and selective cholesteryl ester uptake of HDL. We have concluded that SR-BI plays a key role in SAA metab. through its ability to interact with and internalize SAA and, further, that SAA influences HDL cholesterol metab. through its inhibitory effects on SR-BI-mediated selective lipid uptake.
- 8Wang, W.; Khatua, P.; Hansmann, U. H. Cleavage, Downregulation, and Aggregation of Serum Amyloid A. J. Phys. Chem. B 2020, 124, 1009– 1019, DOI: 10.1021/acs.jpcb.9b10843[ACS Full Text
], [CAS], Google Scholar
8https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhtlWnsLs%253D&md5=9e32a66ffbdd75e1b8d3b4a8cc94d18aCleavage, downregulation, and aggregation of serum amyloid AWang, Wenhua; Khatua, Prabir; Hansmann, Ulrich H. E.Journal of Physical Chemistry B (2020), 124 (6), 1009-1019CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)Various diseases cause overexpression of the serum amyloid A (SAA) protein, which in some cases, but not in all cases, leads to amyloidosis as a secondary disease. Response to the overexpression involves dissocn. of the SAA hexamer and subsequent cleavage of the released monomers, most commonly yielding fragments SAA1-76 of the full-sized SAA1-104. We report results from mol. dynamic simulations that probe the role of this cleavage for downregulating the activity and concn. of SAA. We propose a mechanism that relies on two elements. First, the probability to assemble into hexamers is lower for the fragments than it is for the full-sized protein. Second, unlike other fragments, SAA1-76 can switch between two distinct configurations. The first kind is easy to proteolyze (allowing a fast redn. of the SAA concn.) but prone to aggregation, whereas the situation is opposite for the second kind. If the time scale for amyloid formation is longer than the one for proteolysis, the aggregation-prone species dominates. However, if environmental conditions such as low pH increases the risk of amyloid formation, the ensemble shifts toward the more protected form. We speculate that SAA amyloidosis is a failure of this switching mechanism leading to accumulation of the aggregation-prone species and subsequent amyloid formation. - 9Yaşar, F.; Bernhardt, N. A.; Hansmann, U. H. Replica-exchange-with-tunneling for fast exploration of protein landscapes. J. Chem. Phys. 2015, 143, 224102, DOI: 10.1063/1.4936968[Crossref], [PubMed], [CAS], Google Scholar9https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXhvF2ksbjF&md5=f496bc29c5375597099d9b54c87378f4Replica-exchange-with-tunneling for fast exploration of protein landscapesYasar, Fatih; Bernhardt, Nathan A.; Hansmann, Ulrich H. E.Journal of Chemical Physics (2015), 143 (22), 224102/1-224102/6CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)While the use of replica-exchange mol. dynamics in protein simulations has become ubiquitous, its utility is limited in many practical applications. The authors propose to overcome some shortcomings that hold back its use in settings such as multi-scale or explicit solvent simulations by integrating ideas of hybrid MC/MD into the replica-exchange protocol. This Replica-Exchange-with-Tunneling method is tested by simulating the Trp-cage protein, a system often used in mol. biophysics for testing sampling techniques. (c) 2015 American Institute of Physics.
- 10Bernhardt, N. A.; Xi, W.; Wang, W.; Hansmann, U. H. Simulating Protein Fold Switching by Replica Exchange with Tunneling. J. Chem. Theory Comput. 2016, 12, 5656– 5666, DOI: 10.1021/acs.jctc.6b00826[ACS Full Text
], [CAS], Google Scholar
10https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhslWjsLjN&md5=06db8e4a33672ac6a9e26208037590c1Simulating Protein Fold Switching by Replica Exchange with TunnelingBernhardt, Nathan A.; Xi, Wenhui; Wang, Wei; Hansmann, Ulrich H. E.Journal of Chemical Theory and Computation (2016), 12 (11), 5656-5666CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)Recent expts. suggest that an amino acid sequence encodes not only the native fold of a protein but also other forms essential for its function, or important during folding or assocn. These various forms populate a multi-funnel folding and assocn. landscape where mutations, changes in environment, or interaction with other mols. switch between the encoded folds. Here, the authors introduced replica-exchange-with-tunneling as a way to simulate efficiently switching between distinct folds of proteins and protein aggregates. The correctness and efficiency of this approach was demonstrated in a series of simulations covering a wide range of proteins, from a small 11-residue large designed peptide up to 2 56-residue large mutants of the A and B domain of streptococcal protein G. - 11Zhang, H.; Xi, W.; Hansmann, U. H.; Wei, Y. Fibril-Barrel Transitions in Cylindrin Amyloids. J. Chem. Theory Comput. 2017, 13, 3936– 3944, DOI: 10.1021/acs.jctc.7b00383[ACS Full Text
], [CAS], Google Scholar
11https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtFSjtrrO&md5=d144e2f7757c674122513919ccff1640Fibril-Barrel Transitions in Cylindrin AmyloidsZhang, Huiling; Xi, Wenhui; Hansmann, Ulrich H. E.; Wei, YanjieJournal of Chemical Theory and Computation (2017), 13 (8), 3936-3944CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We introduce Replica-Exchange-with-Tunneling (RET) simulations as a tool for studies of the conversion between polymorph amyloids. For the eleven-residue amyloid- forming cylindrin peptide we show that this technique allows for a more efficient sampling of the formation and interconversion between fibril-like and barrel-like assemblies. We describe a protocol for optimized anal. of RET simulations that allows us to propose a mechanism for formation and interconversion between various cylindrin assemblies. Esp., we show that an interchain salt bridge between residues K3 and D7 is crucial for formation of the barrel structure. - 12Bernhardt, N. A.; Hansmann, U. H. Multifunnel Landscape of the Fold-Switching Protein RfaH-CTD. J. Phys. Chem. B 2018, 122, 1600– 1607, DOI: 10.1021/acs.jpcb.7b11352[ACS Full Text
], [CAS], Google Scholar
12https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXnvVCnsQ%253D%253D&md5=2bf1bd93a92bbe0444c16e73cdcd2dc4Multifunnel Landscape of the Fold-Switching Protein RfaH-CTDBernhardt, Nathan A.; Hansmann, Ulrich H. E.Journal of Physical Chemistry B (2018), 122 (5), 1600-1607CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)Proteins such as transcription factor RfaH can change biol. function by switching between distinct 3-dimensional folds. RfaH regulates transcription if the C-terminal domain folds into a double helix bundle, and promotes translation when this domain assumes a β-barrel form. This fold-switch has been also obsd. for the isolated C-terminal domain (RfaH-CTD) and was studied here with a variant of the RET (Replica-Exchange-with-Tunneling) approach recently introduced by us. We used the enhanced sampling properties of this technique to map the free energy landscape of RfaH-CTD and to propose a mechanism for the conversion process. - 13Lu, J.; Yu, Y.; Zhu, I.; Cheng, Y.; Sun, P. D. Structural mechanism of serum amyloid A-mediated inflammatory amyloidosis. Proc. Natl. Acad. Sci. U. S. A. 2014, 111, 5189– 5194, DOI: 10.1073/pnas.1322357111[Crossref], [PubMed], [CAS], Google Scholar13https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXkslGqt7w%253D&md5=cec09402517df4fee69ac5d78acaf3f8Structural mechanism of serum amyloid A-mediated inflammatory amyloidosisLu, Jinghua; Yu, Yadong; Zhu, Iowis; Cheng, Yifan; Sun, Peter D.Proceedings of the National Academy of Sciences of the United States of America (2014), 111 (14), 5189-5194CODEN: PNASA6; ISSN:0027-8424. (National Academy of Sciences)Serum amyloid A (SAA) represents an evolutionarily conserved family of inflammatory acute-phase proteins. It is also a major constituent of secondary amyloidosis. Here, to understand its function and structural transition to amyloid, the authors detd. a structure of human SAA isoform 1.1 (SAA1.1) in 2 crystal forms, representing a prototypic member of the family. Native SAA1.1 exists as a hexamer, with subunits displaying a unique 4-helix bundle fold stabilized by its long C-terminal tail. Structure-based mutational studies revealed 2 pos.-charge clusters, near the center and apex of the hexamer, that were involved in SAA assocn. with heparin. The binding of high-d. lipoprotein involved only the apex region of SAA and could be inhibited by heparin. Peptide amyloid formation assays identified N-terminal helixes 1 and 3 as amyloidogenic peptides of SAA1.1. Both peptides were secluded in the hexameric structure of SAA1.1, suggesting that native SAA is non-pathogenic. Furthermore, dissocn. of the SAA hexamer appeared insufficient to initiate the amyloidogenic transition, and proteolytic cleavage or removal of the C-terminal tail of SAA resulted in the formation of various-sized structural aggregates contg. ∼5-nm regular repeating protofibril-like units. The combined structural and functional studies provided mechanistic insights into the pathogenic contribution of glycosaminoglycan in SAA1.1-mediated AA amyloid formation.
- 14Liberta, F.; Loerch, S.; Rennegarbe, M.; Schierhorn, A.; Westermark, P.; Westermark, G. T.; Hazenberg, B. P.; Grigorieff, N.; Fändrich, M.; Schmidt, M. Cryo-EM fibril structures from systemic AA amyloidosis reveal the species complementarity of pathological amyloids. Nature Commun. 2019, 10, 1104, DOI: 10.1038/s41467-019-09033-z[Crossref], [PubMed], [CAS], Google Scholar14https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BB3cbhs12nug%253D%253D&md5=af5ff643e7e3d8503a088d93c58df447Cryo-EM fibril structures from systemic AA amyloidosis reveal the species complementarity of pathological amyloidsLiberta Falk; Rennegarbe Matthies; Fandrich Marcus; Schmidt Matthias; Loerch Sarah; Grigorieff Nikolaus; Schierhorn Angelika; Westermark Per; Westermark Gunilla T; Hazenberg Bouke P CNature communications (2019), 10 (1), 1104 ISSN:.Systemic AA amyloidosis is a worldwide occurring protein misfolding disease of humans and animals. It arises from the formation of amyloid fibrils from the acute phase protein serum amyloid A. Here, we report the purification and electron cryo-microscopy analysis of amyloid fibrils from a mouse and a human patient with systemic AA amyloidosis. The obtained resolutions are 3.0 ÅA and 2.7 ÅA for the murine and human fibril, respectively. The two fibrils differ in fundamental properties, such as presence of right-hand or left-hand twisted cross-β sheets and overall fold of the fibril proteins. Yet, both proteins adopt highly similar β-arch conformations within the N-terminal ~21 residues. Our data demonstrate the importance of the fibril protein N-terminus for the stability of the analyzed amyloid fibril morphologies and suggest strategies of combating this disease by interfering with specific fibril polymorphs.
- 15Bansal, A.; Schmidt, M.; Rennegarbe, M.; Haupt, C.; Liberta, F.; Stecher, S.; Puscalau-Girtu, I.; Biedermann, A.; Fändrich, M. AA amyloid fibrils from diseased tissue are structurally different from in vitro formed SAA fibrils. Nature Commun. 2021, 12, 1013, DOI: 10.1038/s41467-021-21129-z[Crossref], [PubMed], [CAS], Google Scholar15https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXksV2kt7Y%253D&md5=d38267d13a31b0be0bce42df1582fe08AA amyloid fibrils from diseased tissue are structurally different from in vitro formed SAA fibrilsBansal, Akanksha; Schmidt, Matthias; Rennegarbe, Matthies; Haupt, Christian; Liberta, Falk; Stecher, Sabrina; Puscalau-Girtu, Ioana; Biedermann, Alexander; Faendrich, MarcusNature Communications (2021), 12 (1), 1013CODEN: NCAOBW; ISSN:2041-1723. (Nature Research)Systemic AA amyloidosis is a world-wide occurring protein misfolding disease of humans and animals. It arises from the formation of amyloid fibrils from serum amyloid A (SAA) protein. Using cryo electron microscopy we here show that amyloid fibrils which were purified from AA amyloidotic mice are structurally different from fibrils formed from recombinant SAA protein in vitro. Ex vivo amyloid fibrils consist of fibril proteins that contain more residues within their ordered parts and possess a higher β-sheet content than in vitro fibril proteins. They are also more resistant to proteolysis than their in vitro formed counterparts. These data suggest that pathogenic amyloid fibrils may originate from proteolytic selection, allowing specific fibril morphologies to proliferate and to cause damage to the surrounding tissue.
- 16Fukunishi, H.; Watanabe, O.; Takada, S. On the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure prediction. J. Chem. Phys. 2002, 116, 9058– 9067, DOI: 10.1063/1.1472510[Crossref], [CAS], Google Scholar16https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38XjsFKmsLo%253D&md5=7ac571a5afdd63b0b4b29cfdec06f53bOn the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure predictionFukunishi, Hiroaki; Watanabe, Osamu; Takada, ShojiJournal of Chemical Physics (2002), 116 (20), 9058-9067CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Motivated by the protein structure prediction problem, we develop two variants of the Hamiltonian replica exchange methods (REMs) for efficient configuration sampling, (1) the scaled hydrophobicity REM and (2) the phantom chain REM, and compare their performance with the ordinary REM. We first point out that the ordinary REM has a shortage for the application to large systems such as biomols. and that the Hamiltonian REM, an alternative formulation of the REM, can give a remedy for it. We then propose two examples of the Hamiltonian REM that are suitable for a coarse-grained protein model. (1) The scaled hydrophobicity REM preps. replicas that are characterized by various strengths of hydrophobic interaction. The strongest interaction that mimics aq. soln. environment makes proteins folding, while weakened hydrophobicity unfolds proteins as in org. solvent. Exchange between these environments enables proteins to escape from misfolded traps and accelerate conformational search. This resembles the roles of mol. chaperone that assist proteins to fold in vivo. (2) The phantom chain REM uses replicas that allow various degrees of at. overlaps. By allowing at. overlap in some of replicas, the peptide chain can cross over itself, which can accelerate conformation sampling. Using a coarse-gained model we developed, we compute equil. probability distributions for poly-alanine 16-mer and for a small protein by these REMs and compare the accuracy of the results. We see that the scaled hydrophobicity REM is the most efficient method among the three REMs studied.
- 17Kwak, W.; Hansmann, U. H. Efficient sampling of protein structures by model hopping. Phys. Rev. Lett. 2005, 95, 138102, DOI: 10.1103/PhysRevLett.95.138102[Crossref], [PubMed], [CAS], Google Scholar17https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXhtVelsL7K&md5=8d057f6efdfce47c6ed7a245a4c0da77Efficient Sampling of Protein Structures by Model HoppingKwak, Wooseop; Hansmann, Ulrich H. E.Physical Review Letters (2005), 95 (13), 138102/1-138102/4CODEN: PRLTAO; ISSN:0031-9007. (American Physical Society)A review. We introduce a novel simulation method, model hopping, that enhances sampling of low-energy configurations in complex systems. The approach is illustrated for a protein-folding problem. Thermodn. quantities of proteins with up to 46 residues are evaluated from all-atom simulations with this method.
- 18Hansmann, U. H. Parallel tempering algorithm for conformational studies of biological molecules. Chem. Phys. Lett. 1997, 281, 140– 150, DOI: 10.1016/S0009-2614(97)01198-6[Crossref], [CAS], Google Scholar18https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2sXotFWktr4%253D&md5=ddfed8309757ca88a834a1cf3c80cb08Parallel tempering algorithm for conformational studies of biological moleculesHansmann, Ulrich H. E.Chemical Physics Letters (1997), 281 (1,2,3), 140-150CODEN: CHPLBC; ISSN:0009-2614. (Elsevier Science B.V.)The effectiveness of a new algorithm, parallel tempering, is studied for numerical simulations of biol. mols. These mols. suffer from a rough energy landscape. The resulting slowing down in numerical simulations is overcome by the new method. This is demonstrated by performing simulations with high statistics for one of the simplest peptides, Met-enkephalin. The numerical effectiveness of the new technique was found to be much better than traditional methods and is comparable to sophisticated methods like generalized ensemble techniques.
- 19Best, R. B.; Zhu, X.; Shim, J.; Lopes, P. E.; Mittal, J.; Feig, M.; MacKerell, A. D. Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone ϕ, ψ and side-chain χ1 and χ2 Dihedral Angles. J. Chem. Theory Comput. 2012, 8, 3257– 3273, DOI: 10.1021/ct300400x[ACS Full Text
], [CAS], Google Scholar
19https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38XhtVKqurfP&md5=9a48a0c5770fb1e887c3bb34d45b1354Optimization of the Additive CHARMM All-Atom Protein Force Field Targeting Improved Sampling of the Backbone .vphi., ψ and Side-Chain χ1 and χ2 Dihedral AnglesBest, Robert B.; Zhu, Xiao; Shim, Jihyun; Lopes, Pedro E. M.; Mittal, Jeetain; Feig, Michael; MacKerell, Alexander D.Journal of Chemical Theory and Computation (2012), 8 (9), 3257-3273CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)While the quality of the current CHARMM22/CMAP additive force field for proteins has been demonstrated in a large no. of applications, limitations in the model with respect to the equil. between the sampling of helical and extended conformations in folding simulations have been noted. To overcome this, as well as make other improvements in the model, we present a combination of refinements that should result in enhanced accuracy in simulations of proteins. The common (non-Gly, -Pro) backbone CMAP potential has been refined against exptl. soln. NMR data for weakly structured peptides, resulting in a rebalancing of the energies of the α-helix and extended regions of the Ramachandran map, correcting the α-helical bias of CHARMM22/CMAP. The Gly and Pro CMAPs have been refitted to more accurate quantum-mech. energy surfaces. Side-chain torsion parameters have been optimized by fitting to backbone-dependent quantum-mech. energy surfaces, followed by addnl. empirical optimization targeting NMR scalar couplings for unfolded proteins. A comprehensive validation of the revised force field was then performed against a collection of exptl. data: (i) comparison of simulations of eight proteins in their crystal environments with crystal structures; (ii) comparison with backbone scalar couplings for weakly structured peptides; (iii) comparison with NMR residual dipolar couplings and scalar couplings for both backbone and side-chains in folded proteins; (iv) equil. folding of mini-proteins. The results indicate that the revised CHARMM 36 parameters represent an improved model for modeling and simulation studies of proteins, including studies of protein folding, assembly, and functionally relevant conformational changes. - 20Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 1983, 79, 926– 935, DOI: 10.1063/1.445869[Crossref], [CAS], Google Scholar20https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL3sXksF2htL4%253D&md5=a1161334e381746be8c9b15a5e56f704Comparison of simple potential functions for simulating liquid waterJorgensen, William L.; Chandrasekhar, Jayaraman; Madura, Jeffry D.; Impey, Roger W.; Klein, Michael L.Journal of Chemical Physics (1983), 79 (2), 926-35CODEN: JCPSA6; ISSN:0021-9606.Classical Monte Carlo simulations were carried out for liq. H2O in the NPT ensemble at 25° and 1 atm using 6 of the simpler intermol. potential functions for the dimer. Comparisons were made with exptl. thermodn. and structural data including the neutron diffraction results of Thiessen and Narten (1982). The computed densities and potential energies agree with expt. except for the original Bernal-Fowler model, which yields an 18% overest. of the d. and poor structural results. The discrepancy may be due to the correction terms needed in processing the neutron data or to an effect uniformly neglected in the computations. Comparisons were made for the self-diffusion coeffs. obtained from mol. dynamics simulations.
- 21Pettersen, E. F.; Goddard, T. D.; Huang, C. C.; Couch, G. S.; Greenblatt, D. M.; Meng, E. C.; Ferrin, T. E. UCSF Chimera - A visualization system for exploratory research and analysis. J. Comput. Chem. 2004, 25, 1605– 1612, DOI: 10.1002/jcc.20084[Crossref], [PubMed], [CAS], Google Scholar21https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXmvVOhsbs%253D&md5=944b175f440c1ff323705987cf937ee7UCSF Chimera-A visualization system for exploratory research and analysisPettersen, Eric F.; Goddard, Thomas D.; Huang, Conrad C.; Couch, Gregory S.; Greenblatt, Daniel M.; Meng, Elaine C.; Ferrin, Thomas E.Journal of Computational Chemistry (2004), 25 (13), 1605-1612CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)The design, implementation, and capabilities of an extensible visualization system, UCSF Chimera, are discussed. Chimera is segmented into a core that provides basic services and visualization, and extensions that provide most higher level functionality. This architecture ensures that the extension mechanism satisfies the demands of outside developers who wish to incorporate new features. Two unusual extensions are presented: Multiscale, which adds the ability to visualize large-scale mol. assemblies such as viral coats, and Collab., which allows researchers to share a Chimera session interactively despite being at sep. locales. Other extensions include Multalign Viewer, for showing multiple sequence alignments and assocd. structures; ViewDock, for screening docked ligand orientations; Movie, for replaying mol. dynamics trajectories; and Vol. Viewer, for display and anal. of volumetric data. A discussion of the usage of Chimera in real-world situations is given, along with anticipated future directions. Chimera includes full user documentation, is free to academic and nonprofit users, and is available for Microsoft Windows, Linux, Apple Mac OS X, SGI IRIX, and HP Tru64 Unix from http://www.cgl.ucsf.edu/chimera/.
- 22Noel, J. K.; Whitford, P. C.; Sanbonmatsu, K. Y.; Onuchic, J. N. SMOG@ctbp: Simplified deployment of structure-based models in GROMACS. Nucleic Acids Res. 2010, 38, W657, DOI: 10.1093/nar/gkq498[Crossref], [PubMed], [CAS], Google Scholar22https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXotVSqt7s%253D&md5=937e9d055c9a72ef70e2deb3002f3813SMOG@ctbp: simplified deployment of structure-based models in GROMACSNoel, Jeffrey K.; Whitford, Paul C.; Sanbonmatsu, Karissa Y.; Onuchic, Jose N.Nucleic Acids Research (2010), 38 (Web Server), W657-W661CODEN: NARHAD; ISSN:0305-1048. (Oxford University Press)Mol. dynamics simulations with coarse-grained and/or simplified Hamiltonians are an effective means of capturing the functionally important long-time and large-length scale motions of proteins and RNAs. Structure-based Hamiltonians, simplified models developed from the energy landscape theory of protein folding, have become a std. tool for investigating biomol. dynamics. SMOG@ctbp is an effort to simplify the use of structure-based models. The purpose of the web server is two fold. First, the web tool simplifies the process of implementing a well-characterized structure-based model on a state-of-the-art, open source, mol. dynamics package, GROMACS. Second, the tutorial-like format helps speed the learning curve of those unfamiliar with mol. dynamics. A web tool user is able to upload any multi-chain biomol. system consisting of std. RNA, DNA and amino acids in PDB format and receive as output all files necessary to implement the model in GROMACS. Both Cα and all-atom versions of the model are available. SMOG@ctbp resides at http://smog.ucsd.edu.
- 23Zhang, W.; Chen, J. Accelerate sampling in atomistic energy landscapes using topology-based coarse-grained models. J. Chem. Theory Comput. 2014, 10, 918– 923, DOI: 10.1021/ct500031v[ACS Full Text
], [CAS], Google Scholar
23https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhsVOjtro%253D&md5=e8368763a3870dac4a993d744da02ba8Accelerate Sampling in Atomistic Energy Landscapes Using Topology-Based Coarse-Grained ModelsZhang, Weihong; Chen, JianhanJournal of Chemical Theory and Computation (2014), 10 (3), 918-923CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We describe a multiscale enhanced sampling (MSES) method where efficient topol.-based coarse-grained models are coupled with all-atom ones to enhance the sampling of atomistic protein energy landscape. The bias from the coupling is removed by Hamiltonian replica exchange, thus allowing one to benefit simultaneously from faster transitions of coarse-grained modeling and accuracy of atomistic force fields. The method is demonstrated by calcg. the conformational equil. of several small but nontrivial β-hairpins with varied stabilities. - 24Moritsugu, K.; Terada, T.; Kidera, A. Scalable free energy calculation of proteins via multiscale essential sampling. J. Chem. Phys. 2010, 133, 224105, DOI: 10.1063/1.3510519[Crossref], [PubMed], [CAS], Google Scholar24https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXhsFGksbrM&md5=1e8de2879382edb9b892f1cac9933f82Scalable free energy calculation of proteins via multiscale essential samplingMoritsugu, Kei; Terada, Tohru; Kidera, AkinoriJournal of Chemical Physics (2010), 133 (22), 224105/1-224105/6CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A multiscale simulation method, multiscale essential sampling (MSES), is proposed for calcg. free energy surface of proteins in a sizable dimensional space with good scalability. In MSES, the configurational sampling of a full-dimensional model is enhanced by coupling with the accelerated dynamics of the essential degrees of freedom. Applying the Hamiltonian exchange method to MSES can remove the biasing potential from the coupling term, deriving the free energy surface of the essential degrees of freedom. The form of the coupling term ensures good scalability in the Hamiltonian exchange. As a test application, the free energy surface of the folding process of a mini-protein, chignolin, was calcd. in the continuum solvent model. Results agreed with the free energy surface derived from the multi-canonical simulation. Significantly improved scalability with the MSES method was clearly shown in the free energy calcn. of chignolin in explicit solvent, which was achieved without increasing the no. of replicas in the Hamiltonian exchange. (c) 2010 American Institute of Physics.
- 25Hess, B.; Kutzner, C.; Van Der Spoel, D.; Lindahl, E. GROMACS 4: Algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput. 2008, 4, 435– 447, DOI: 10.1021/ct700301q[ACS Full Text
], [CAS], Google Scholar
25https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXhsVSqurc%253D&md5=d53c94901386260221792ea30f151c5fGROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular SimulationHess, Berk; Kutzner, Carsten; van der Spoel, David; Lindahl, ErikJournal of Chemical Theory and Computation (2008), 4 (3), 435-447CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)Mol. simulation is an extremely useful, but computationally very expensive tool for studies of chem. and biomol. systems. Here, we present a new implementation of our mol. simulation toolkit GROMACS which now both achieves extremely high performance on single processors from algorithmic optimizations and hand-coded routines and simultaneously scales very well on parallel machines. The code encompasses a minimal-communication domain decompn. algorithm, full dynamic load balancing, a state-of-the-art parallel constraint solver, and efficient virtual site algorithms that allow removal of hydrogen atom degrees of freedom to enable integration time steps up to 5 fs for atomistic simulations also in parallel. To improve the scaling properties of the common particle mesh Ewald electrostatics algorithms, we have in addn. used a Multiple-Program, Multiple-Data approach, with sep. node domains responsible for direct and reciprocal space interactions. Not only does this combination of algorithms enable extremely long simulations of large systems but also it provides that simulation performance on quite modest nos. of std. cluster nodes. - 26Hess, B. P-LINCS: A parallel linear constraint solver for molecular simulation. J. Chem. Theory Comput. 2008, 4, 116– 122, DOI: 10.1021/ct700200b[ACS Full Text
], [CAS], Google Scholar
26https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXhtlKru7zL&md5=476d5ca2eb25574d44b775996fff7b75P-LINCS: A Parallel Linear Constraint Solver for Molecular SimulationHess, BerkJournal of Chemical Theory and Computation (2008), 4 (1), 116-122CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)By removing the fastest degrees of freedom, constraints allow for an increase of the time step in mol. simulations. In the last decade parallel simulations have become commonplace. However, up till now efficient parallel constraint algorithms have not been used with domain decompn. In this paper the parallel linear constraint solver (P-LINCS) is presented, which allows the constraining of all bonds in macromols. Addnl. the energy conservation properties of (P-)LINCS are assessed in view of improvements in the accuracy of uncoupled angle constraints and integration in single precision. - 27Swope, W. C.; Andersen, H. C.; Berens, P. H.; Wilson, K. R. A computer simulation method for the calculation of equilibrium constants for the formation of physical clusters of molecules: Application to small water clusters. J. Chem. Phys. 1982, 76, 637– 649, DOI: 10.1063/1.442716[Crossref], [CAS], Google Scholar27https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL38XlvVWhtg%253D%253D&md5=696100391d72d68cb6eee764c9691d01A computer simulation method for the calculation of equilibrium constants for the formation of physical clusters of molecules: Application to small water clustersSwope, William C.; Andersen, Hans C.; Berens, Peter H.; Wilson, Kent R.Journal of Chemical Physics (1982), 76 (1), 637-49CODEN: JCPSA6; ISSN:0021-9606.A mol. dynamics computer simulation method is described for calcg. equil. consts. for the formation of phys. clusters of mols. The method is based on Hill's formal theory of phys. clusters. A mol. dynamics calcn. is used to calc. the av. potential energy of a cluster of mols. as a function of temp., and the equil. consts. are calcd. from the integral of the energy with respect to reciprocal temp. The method is illustrated by calcns. of the equil. consts. for the formation of clusters of 2-5 water mols. that interact with each other by an intermol. potential devised by Watts. The method is compared with other procedures for calcg. the thermodn. properties of clusters.
- 28Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 2007, 126, 014101, DOI: 10.1063/1.2408420[Crossref], [PubMed], [CAS], Google Scholar28https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXosVCltg%253D%253D&md5=9c182b57bfc65bca6be23c8c76b4be77Canonical sampling through velocity rescalingBussi, Giovanni; Donadio, Davide; Parrinello, MicheleJournal of Chemical Physics (2007), 126 (1), 014101/1-014101/7CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The authors present a new mol. dynamics algorithm for sampling the canonical distribution. In this approach the velocities of all the particles are rescaled by a properly chosen random factor. The algorithm is formally justified and it is shown that, in spite of its stochastic nature, a quantity can still be defined that remains const. during the evolution. In numerical applications this quantity can be used to measure the accuracy of the sampling. The authors illustrate the properties of this new method on Lennard-Jones and TIP4P water models in the solid and liq. phases. Its performance is excellent and largely independent of the thermostat parameter also with regard to the dynamic properties.
- 29McGibbon, R. T.; Beauchamp, K. A.; Harrigan, M. P.; Klein, C.; Swails, J. M.; Hernández, C. X.; Schwantes, C. R.; Wang, L. P.; Lane, T. J.; Pande, V. S. MDTraj: A Modern Open Library for the Analysis of Molecular Dynamics Trajectories. Biophys. J. 2015, 109, 1528– 1532, DOI: 10.1016/j.bpj.2015.08.015[Crossref], [PubMed], [CAS], Google Scholar29https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXhsVKhu7%252FI&md5=9cbd2138ce200e7d63fda8ae756df43fMDTraj: A Modern Open Library for the Analysis of Molecular Dynamics TrajectoriesMcGibbon, Robert T.; Beauchamp, Kyle A.; Harrigan, Matthew P.; Klein, Christoph; Swails, Jason M.; Hernandez, Carlos X.; Schwantes, Christian R.; Wang, Lee-Ping; Lane, Thomas J.; Pande, Vijay S.Biophysical Journal (2015), 109 (8), 1528-1532CODEN: BIOJAU; ISSN:0006-3495. (Cell Press)As mol. dynamics (MD) simulations continue to evolve into powerful computational tools for studying complex biomol. systems, the necessity of flexible and easy-to-use software tools for the anal. of these simulations is growing. We have developed MDTraj, a modern, lightwt., and fast software package for analyzing MD simulations. MDTraj reads and writes trajectory data in a wide variety of commonly used formats. It provides a large no. of trajectory anal. capabilities including minimal root-mean-square-deviation calcns., secondary structure assignment, and the extn. of common order parameters. The package has a strong focus on interoperability with the wider scientific Python ecosystem, bridging the gap between MD data and the rapidly growing collection of industry-std. statistical anal. and visualization tools in Python. MDTraj is a powerful and user-friendly software package that simplifies the anal. of MD data and connects these datasets with the modern interactive data science software ecosystem in Python.
- 30Dijkstra, E. W. A note on two problems in connexion with graphs. Numer. Math. 1959, 1, 269– 271, DOI: 10.1007/BF01386390
- 31Marcos-Alcalde, I.; Setoain, J.; Mendieta-Moreno, J. I.; Mendieta, J.; Gómez-Puertas, P. MEPSA: minimum energy pathway analysis for energy landscapes. Bioinformatics 2015, 31, 3853– 3855, DOI: 10.1093/bioinformatics/btv453[Crossref], [PubMed], [CAS], Google Scholar31https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28Xht1CitLzJ&md5=aea18fdc9e14b7eda6dfacea20145116MEPSA: minimum energy pathway analysis for energy landscapesMarcos-Alcalde, Inigo; Setoain, Javier; Mendieta-Moreno, Jesus I.; Mendieta, Jesus; Gomez-Puertas, PaulinoBioinformatics (2015), 31 (23), 3853-3855CODEN: BOINFP; ISSN:1367-4803. (Oxford University Press)Summary: From conformational studies to atomistic descriptions of enzymic reactions, potential and free energy landscapes can be used to describe biomol. systems in detail. However, extg. the relevant data of complex 3D energy surfaces can sometimes be laborious. In this article, we present MEPSA (Min. Energy Path Surface Anal.), a cross-platform user friendly tool for the anal. of energy landscapes from a transition state theory perspective. Some of its most relevant features are: identification of all the barriers and min. of the landscape at once, description of maxima edge profiles, detection of the lowest energy path connecting two min. and generation of transition state theory diagrams along these paths. In addn. to a built-in plotting system, MEPSA can save most of the generated data into easily parseable text files, allowing more versatile uses of MEPSA's output such as the generation of mol. dynamics restraints from a calcd. path.
- 32Khatua, P.; Ray, A. J.; Hansmann, U. H. Bifurcated Hydrogen Bonds and the Fold Switching of Lymphotactin. J. Phys. Chem. B 2020, 124, 6555– 6564, DOI: 10.1021/acs.jpcb.0c04565[ACS Full Text
], [CAS], Google Scholar
32https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhtlSisL7F&md5=35d96fe846f57c94f6546bc33602fa47Bifurcated Hydrogen Bonds and the Fold Switching of LymphotactinKhatua, Prabir; Ray, Alan J.; Hansmann, Ulrich H. E.Journal of Physical Chemistry B (2020), 124 (30), 6555-6564CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)Lymphotactin (Ltn) exists under physiol. conditions in an equil. between two interconverting structures with distinct biol. functions. Using replica-exchange-with-tunneling, we study the conversion between the 2-folds. Unlike previously proposed, we find that the fold switching does not require unfolding of lymphotactin but proceeds through a series of intermediates that remain partially structured. This process relies on two bifurcated hydrogen bonds that connect the β2 and β3 strands and ease the transition between the hydrogen bond pattern by which the central three-stranded β-sheet in the two forms differs. - 33Bolhuis, P. G.; Chandler, D.; Dellago, C.; Geissler, P. L. Transition Path Sampling: Throwing Ropes over Rough Mountain Passes, in the Dark. Annu. Rev. Phys. Chem. 2002, 53, 291– 318, DOI: 10.1146/annurev.physchem.53.082301.113146[Crossref], [PubMed], [CAS], Google Scholar33https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38Xks1Ois7g%253D&md5=fbf79b52e751dfd87a73c345b9581898Transition path sampling: throwing ropes over rough mountain passes, in the darkBolhuis, Peter G.; Chandler, David; Dellago, Christoph; Geissler, Phillip L.Annual Review of Physical Chemistry (2002), 53 (), 291-318CODEN: ARPLAP; ISSN:0066-426X. (Annual Reviews Inc.)A review is given of the concepts and methods of transition path sampling. These methods allow computational studies of rare events without requiring prior knowledge of mechanisms, reaction coordinates, and transition states. Based upon a statistical mechanics of trajectory space, they provide a perspective with which time dependent phenomena, even for systems driven far from equil., can be examd. with the same types of importance sampling tools that in the past have been applied so successfully to static equil. properties.
- 34Dellago, C.; Bolhuis, P. G.; Csajka, F. S.; Chandler, D. Transition Path Sampling and the Calculation of Rate Constants. J. Chem. Phys. 1998, 108, 1964– 1977, DOI: 10.1063/1.475562[Crossref], [CAS], Google Scholar34https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1cXkvFChsQ%253D%253D&md5=4bb2f7b8316dbb4526be81bd24083814Transition path sampling and the calculation of rate constantsDellago, Christoph; Bolhuis, Peter G.; Csajka, Felix S.; Chandler, DavidJournal of Chemical Physics (1998), 108 (5), 1964-1977CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We have developed a method to study transition pathways for rare events in complex systems. The method can be used to det. rate consts. for transitions between stable states by turning the calcn. of reactive flux correlation functions into the computation of an isomorphic reversible work. In contrast to previous dynamical approaches, the method relies neither on prior knowledge nor on explicit specification of transition states. Rather, it provides an importance sampling from which transition states can be characterized statistically. A simple model is analyzed to illustrate the methodol.
- 35Dellago, C.; Bolhuis, P. G.; Chandler, D. Efficient Transition Path Sampling: Application to Lennard-Jones Cluster Rearrangements. J. Chem. Phys. 1998, 108, 9236– 9245, DOI: 10.1063/1.476378[Crossref], [CAS], Google Scholar35https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1cXjtFSitb8%253D&md5=0d02dabdf6c5811d0bf11847d8949c7eEfficient transition path sampling: Application to Lennard-Jones cluster rearrangementsDellago, Christoph; Bolhuis, Peter G.; Chandler, DavidJournal of Chemical Physics (1998), 108 (22), 9236-9245CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We develop an efficient Monte-Carlo algorithm to sample an ensemble of stochastic transition paths between stable states. In our description, paths are represented by chains of states linked by Markovian transition probabilities. Rate consts. and mechanisms characterizing the transition may be detd. from the path ensemble. We have previously devised several algorithms for sampling the path ensemble. For these algorithms, the numerical effort scales with the square of the path length. In the new simulation scheme, the required computation scales linearly with the length of the transition path. This improved efficiency allows the calcn. of rate consts. in complex mol. systems. As an example, we study rearrangement processes in a cluster consisting of seven Lennard-Jones particles in two dimensions. Using a quenching technique we are able to identify the relevant transition mechanisms and to locate the related transition states. We furthermore calc. transition rate consts. for various isomerization processes.
- 36Weinan, E.; Ren, W.; Vanden-Eijnden, E. String Method for the Study of Rare Events. Phys. Rev. B 2002, 66, 052301, DOI: 10.1103/PhysRevB.66.052301
- 37Maragliano, L.; Fischer, A.; Vanden-Eijnden, E.; Ciccotti, G. String Method in Collective Variables: Minimum Free Energy Paths and Isocommittor Surfaces. J. Chem. Phys. 2006, 125, 024106, DOI: 10.1063/1.2212942[Crossref], [PubMed], [CAS], Google Scholar37https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28Xnt1Khurk%253D&md5=b3a7b0b167df6980ac6c8e7bb36b16fcString method in collective variables: Minimum free energy paths and isocommittor surfacesMaragliano, Luca; Fischer, Alexander; Vanden-Eijnden, Eric; Ciccotti, GiovanniJournal of Chemical Physics (2006), 125 (2), 024106/1-024106/15CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A computational technique is proposed which combines the string method with a sampling technique to det. min. free energy paths. The technique only requires to compute the mean force and another conditional expectation locally along the string, and therefore can be applied even if the no. of collective variables kept in the free energy calcn. is large. This is in contrast with other free energy sampling techniques which aim at mapping the full free energy landscape and whose cost increases exponentially with the no. of collective variables kept in the free energy. Provided that the no. of collective variables is large enough, the new technique captures the mechanism of transition in that it allows to det. the committor function for the reaction and, in particular, the transition state region. The new technique is illustrated on the example of alanine dipeptide, in which we compute the min. free energy path for the isomerization transition using either two or four dihedral angles as collective variables. It is shown that the mechanism of transition can be captured using the four dihedral angles, but it cannot be captured using only two of them.
- 38Unarta, I. C.; Zhu, L.; Tse, C. K. M.; Cheung, P. P.-H.; Yu, J.; Huang, X. Molecular Mechanisms of RNA Polymerase II Transcription Elongation Elucidated by Kinetic Network Models. Curr. Opin. Struct. Biol. 2018, 49, 54– 62, DOI: 10.1016/j.sbi.2018.01.002[Crossref], [PubMed], [CAS], Google Scholar38https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhsVensbg%253D&md5=efc857e4ccb759523af7030715a53d1fMolecular mechanisms of RNA polymerase II transcription elongation elucidated by kinetic network modelsUnarta, Ilona Christy; Zhu, Lizhe; Tse, Carmen Ka Man; Cheung, Peter Pak-Hang; Yu, Jin; Huang, XuhuiCurrent Opinion in Structural Biology (2018), 49 (), 54-62CODEN: COSBEF; ISSN:0959-440X. (Elsevier Ltd.)Transcription elongation cycle (TEC) of RNA polymerase II (Pol II) is a process of adding a nucleoside triphosphate to the growing mRNA chain. Due to the long timescale events in Pol II TEC, an advanced computational technique, such as Markov State Model (MSM), is needed to provide atomistic mechanism and reaction rates. The combination of MSM and exptl. results can be used to build a kinetic network model (KNM) of the whole TEC. This review provides a brief protocol to build MSM and KNM of the whole TEC, along with the latest findings of MSM and other computational studies of Pol II TEC. Lastly, we offer a perspective on potentially using a sequence dependent KNM to predict genome-wide transcription error.
- 39Zhu, L.; Sheong, F. K.; Cao, S.; Liu, S.; Unarta, I. C.; Huang, X. TAPS: A Traveling-Salesman Based Automated Path Searching Method for Functional Conformational Changes of Biological Macromolecules. J. Chem. Phys. 2019, 150, 124105, DOI: 10.1063/1.5082633[Crossref], [PubMed], [CAS], Google Scholar39https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXmtVClt74%253D&md5=18ac3c2a1ef78dff772057ec9ad3aff9TAPS: A traveling-salesman based automated path searching method for functional conformational changes of biological macromoleculesZhu, Lizhe; Sheong, Fu Kit; Cao, Siqin; Liu, Song; Unarta, Ilona C.; Huang, XuhuiJournal of Chemical Physics (2019), 150 (12), 124105/1-124105/10CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Locating the min. free energy paths (MFEPs) between two conformational states is among the most important tasks of biomol. simulations. For example, knowledge of the MFEP is crit. for focusing the effort of unbiased simulations that were used for the construction of Markov state models to the biol. relevant regions of the system. Typically, existing path searching methods perform local sampling around the path nodes in a pre-selected collective variable (CV) space to allow a gradual downhill evolution of the path toward the MFEP. Despite the wide application of such a strategy, the gradual path evolution and the non-trivial a priori choice of CVs are also limiting its overall efficiency and automation. Non-local perpendicular sampling can be pursued to accelerate the search, provided that all nodes are reordered thereafter via a traveling-salesman scheme. Moreover, path-CVs can be computed on-the-fly and used as a coordinate system, minimizing the necessary prior knowledge about the system. The authors' traveling-salesman based automated path searching method achieves a 5-8 times speedup over the string method with swarms-of-trajectories for two peptide systems in vacuum and soln., making it a promising method for obtaining initial pathways when studying functional conformational changes between a pair of structures. (c) 2019 American Institute of Physics.
- 40Schrödinger, L. The PyMOL Molecular Graphics System , Version 1.8; 2015.Google ScholarThere is no corresponding record for this reference.
- 41Westermark, G. T.; Engström, U.; Westermark, P. The N-terminal segment of protein AA determines its fibrillogenic property. Biochem. Biophys. Res. Commun. 1992, 182, 27– 33, DOI: 10.1016/S0006-291X(05)80107-X[Crossref], [PubMed], [CAS], Google Scholar41https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK38Xht1egur4%253D&md5=80fe3a456eac7d41a32920d051c0c60dThe N-terminal segment of protein AA determines its fibrillogenic propertyWestermark, Gunilla T.; Engstroem, Ulla; Westermark, PerBiochemical and Biophysical Research Communications (1992), 182 (1), 27-33CODEN: BBRCA9; ISSN:0006-291X.The amyloid fibril protein AA consists of a varying long N-terminal part of the precursor protein serum AA. By using synthetic peptides corresponding to human and murine protein AA segments and cyanogen bromide fragments of human protein AA, evidence is shown that the amyloidogenic part of the mol. is the first 10-15 amino acid long segment. Amino acid substitutions in this part of the mol. may explain why only one of the two mouse SAA isoforms is amyloidogenic.
- 42Yassine, H. N.; Trenchevska, O.; He, H.; Borges, C. R.; Nedelkov, D.; Mack, W.; Kono, N.; Koska, J.; Reaven, P. D.; Nelson, R. W. Serum Amyloid A Truncations in Type 2 Diabetes Mellitus. PLoS One 2015, 10, e0115320, DOI: 10.1371/journal.pone.0115320
- 43Wang, W.; Hansmann, U. H. Stability of Human Serum Amyloid A Fibrils. J. Phys. Chem. B 2020, 124, 10708– 10717, DOI: 10.1021/acs.jpcb.0c08280[ACS Full Text
], [CAS], Google Scholar
43https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXitlagtLbF&md5=b9dce5f2c90c561d7a9c9b50d1829c39Stability of Human Serum Amyloid A FibrilsWang, Wenhua; Hansmann, Ulrich H. E.Journal of Physical Chemistry B (2020), 124 (47), 10708-10717CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)In systemic amyloidosis, serum amyloid A (SAA) fibril deposits cause widespread damages to tissues and organs that eventually may lead to death. A therapeutically intervention therefore has either to dissolve these fibrils or inhibit their formation. However, only recently has the human SAA fibril structure been resolved at a resoln. that is sufficient for development of drug candidates. Here, we use mol. dynamic simulations to probe the factors that modulate the stability of this fibril model. Our simulations suggest that fibril formation starts with the stacking of two misfolded monomers into metastable dimers, with the stacking depending on the N-terminal amyloidogenic regions of different chains forming anchors. The resulting dimers pack in a second step into a 2-fold two-layer tetramer that is stable enough to nucleate fibril formation. The stability of the initial dimers is enhanced under acidic conditions by a strong salt bridge and side-chain hydrogen bond network in the C-terminal cavity (residues 23-51) but is not affected by the presence of the disordered C-terminal tail. - 44Wang, W.; Xi, W.; Hansmann, U. H. Stability of the N-Terminal Helix and Its Role in Amyloid Formation of Serum Amyloid A. ACS Omega 2018, 3, 16184– 16190, DOI: 10.1021/acsomega.8b02377[ACS Full Text
], [CAS], Google Scholar
44https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXitlaktbbL&md5=59f3032458272f0eff6faa7ca0e70f80Stability of the N-Terminal Helix and Its Role in Amyloid Formation of Serum Amyloid AWang, Wenhua; Xi, Wenhui; Hansmann, Ulrich H. E.ACS Omega (2018), 3 (11), 16184-16190CODEN: ACSODF; ISSN:2470-1343. (American Chemical Society)Colonic amyloidosis is the result of overexpression of the serum amyloid A protein (SAA) in inflammatory bowel disease or colon cancer. Crucial for amyloid formation are the 1st ten N-terminal residues, which in the crystal structure are part of a 27-residue long helix. Here we studied this 27-residue N-terminal region of SAA by a multi-exchange variant of replica exchange mol. dynamics (ME-REMD) simulations. An ensemble of configurations was obsd., dominated by 3 motifs: the single helix of the crystal structure, a helix-turn-helix configuration, and such where residues 14-27 were part of a helix but the 1st 13 residues formed an extended and disordered segment that was prone to aggregation. A single point mutation E9A shifted the equil. to the latter motif indicating the importance of interactions involving this residue for the stability of the SAA protein. - 45Frame, N. M.; Kumanan, M.; Wales, T. E.; Bandara, A.; Fändrich, M.; Straub, J. E.; Engen, J. R.; Gursky, O. Structural Basis for Lipid Binding and Function by an Evolutionarily Conserved Protein, Serum Amyloid A. J. Mol. Bio. 2020, 432, 1978– 1995, DOI: 10.1016/j.jmb.2020.01.029[Crossref], [PubMed], [CAS], Google Scholar45https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXivFyntro%253D&md5=892c2c04204d39ad65267b17c04c3fadStructural Basis for Lipid Binding and Function by an Evolutionarily Conserved Protein, Serum Amyloid AFrame, Nicholas M.; Kumanan, Meera; Wales, Thomas E.; Bandara, Asanga; Fandrich, Marcus; Straub, John E.; Engen, John R.; Gursky, OlgaJournal of Molecular Biology (2020), 432 (7), 1978-1995CODEN: JMOBAK; ISSN:0022-2836. (Elsevier Ltd.)Serum amyloid A (SAA) is a plasma protein that transports lipids during inflammation. To explore SAA soln. conformations and lipid-binding mechanism, we used hydrogen-deuterium exchange mass spectrometry, lipoprotein reconstitution, amino acid sequence anal., and mol. dynamics simulations. Soln. conformations of lipid-bound and lipid-free mSAA1 at pH∼7.4 agreed in details with the crystal structures but also showed important differences. The results revealed that amphipathic α-helixes h1 and h3 comprise a lipid-binding site that is partially pre-formed in soln., is stabilized upon binding lipids, and shows lipid-induced folding of h3. This site sequesters apolar ligands via a concave hydrophobic surface in SAA oligomers. The largely disordered/dynamic C-terminal region is conjectured to mediate the promiscuous binding of other ligands. The h1-h2 linker region is predicted to form an unexpected β-hairpin that may represent an early amyloidogenic intermediate. The results help establish structural underpinnings for understanding SAA interactions with its key functional ligands, its evolutional conservation, and its transition to amyloid.
Cited By
This article is cited by 2 publications.
- Sourav Nandi, Anurup Mukhopadhyay, Pratyush Kiran Nandi, Nanigopal Bera, Ritwik Hazra, Jyotirmoy Chatterjee, Nilmoni Sarkar. Amyloids Formed by Nonaromatic Amino Acid Methionine and Its Cross with Phenylalanine Significantly Affects Phospholipid Vesicle Membrane: An Insight into Hypermethioninemia Disorder. Langmuir 2022, 38 (27) , 8252-8265. https://doi.org/10.1021/acs.langmuir.2c00648
- Asis K. Jana, Chance W. Lander, Andrew D. Chesney, Ulrich H. E. Hansmann. Effect of an Amyloidogenic SARS-COV-2 Protein Fragment on α-Synuclein Monomers and Fibrils. The Journal of Physical Chemistry B 2022, 126 (20) , 3648-3658. https://doi.org/10.1021/acs.jpcb.2c01254
Abstract
Figure 1
Figure 2
Figure 2. (a) A typical example of a replica walking through λ space starting from a replica where the physical model is initially biased toward the folded SAA structure. While the system walks between a replica with bias toward the folded structure (upper half) and a replica with bias toward the fibril structure (lower half), its configuration changes accordingly. This can be seen in (b), where we show the corresponding time evolution of the RMSD to the native structure (in magenta) and the fibril structure (in green).
Figure 3
Figure 3. Free energy landscape as obtained from RET simulations, with data taken at λ = 0, i.e., where the physical models are not biased by any Go-term. Energies are listed in units of kT. The prospective transition pathway is drawn in black, and the five regions crossed by this path are marked in capital letters.
Figure 4
Figure 4. Characteristic conformations of of SAA1–76 as seen in each of the five regions (labeled A, B, C, D, and E) identified on the proposed transition pathway. The N-terminus of the chains is colored in blue. Unlike for the transition regions B, C, and D, these conformations superimposed on the respective reference structures for region A (dominated by folded-conformations) and region E (where fibril-like conformations are dominant).
Figure 5
Figure 5. (a) Number of contacts (normalized to one) nNS that are shared with the folded structure as measured in each of the five regions A to E of the transition pathway. The subset of long-range contacts nLR, again normalized to one, is drawn separately. Shown are also the number nFS of contacts also found in the fibril reference structure. In (b), we show the relative frequency with which one of the three characteristic helices of the folded structure, or the two main β-strands of the fibril structure, is observed.
Figure 6
Figure 6. Residue–residue map of the average minimal distance between heavy atoms in a pair of residues, shown for each of the five regions A, B, C, D, and E. Unlike in the transition regions B, C, and D, folded conformations dominate in region A, and fibril-like ones dominate in region E. Residue pairs whose average contact distance is more than 10 Å are excluded. Numbers on the X and Y axis mark residues, and the average distance in Å is given by the color coding shown in the middle row.
References
ARTICLE SECTIONSThis article references 45 other publications.
- 1Chiti, F.; Dobson, C. M. Protein misfolding, amyloid formation, and human disease: A summary of progress over the last decade. Annu. Rev. Biochem. 2017, 86, 27– 68, DOI: 10.1146/annurev-biochem-061516-045115[Crossref], [PubMed], [CAS], Google Scholar1https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXns1Wjsr0%253D&md5=e1802f20bcaacdd5097f5ba13a42e4eaProtein Misfolding, Amyloid Formation, and Human Disease: A Summary of Progress Over the Last DecadeChiti, Fabrizio; Dobson, Christopher M.Annual Review of Biochemistry (2017), 86 (), 27-68CODEN: ARBOAW; ISSN:0066-4154. (Annual Reviews)Peptides and proteins have been found to possess an inherent tendency to convert from their native functional states into intractable amyloid aggregates. This phenomenon is assocd. with a range of increasingly common human disorders, including Alzheimer and Parkinson diseases, type II diabetes, and a no. of systemic amyloidoses. In this review, we describe this field of science with particular ref. to the advances that have been made over the last decade in our understanding of its fundamental nature and consequences. We list the proteins that are known to be deposited as amyloid or other types of aggregates in human tissues and the disorders with which they are assocd., as well as the proteins that exploit the amyloid motif to play specific functional roles in humans. In addn., we summarize the genetic factors that have provided insight into the mechanisms of disease onset. We describe recent advances in our knowledge of the structures of amyloid fibrils and their oligomeric precursors and of the mechanisms by which they are formed and proliferate to generate cellular dysfunction. We show evidence that a complex proteostasis network actively combats protein aggregation and that such an efficient system can fail in some circumstances and give rise to disease. Finally, we anticipate the development of novel therapeutic strategies with which to prevent or treat these highly debilitating and currently incurable conditions.
- 2Sipe, J. D.; Benson, M. D.; Buxbaum, J. N.; Ikeda, S. I.; Merlini, G.; Saraiva, M. J.; Westermark, P. Amyloid fibril proteins and amyloidosis: chemical identification and clinical classification International Society of Amyloidosis 2016 Nomenclature Guidelines. Amyloid 2016, 23, 209– 213, DOI: 10.1080/13506129.2016.1257986[Crossref], [PubMed], [CAS], Google Scholar2https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhvF2msbzI&md5=ffcab18f122fb48c87aba322d03ed0f0Amyloid fibril proteins and amyloidosis: chemical identification and clinical classification International Society of Amyloidosis 2016 Nomenclature GuidelinesSipe, Jean D.; Benson, Merrill D.; Buxbaum, Joel N.; Ikeda, Shu-ichi; Merlini, Giampaolo; Saraiva, Maria J. M.; Westermark, PerAmyloid (2016), 23 (4), 209-213CODEN: AIJIET; ISSN:1350-6129. (Taylor & Francis Ltd.)The Nomenclature Committee of the International Society of Amyloidosis (ISA) met during the XVth Symposium of the Society, 3 July-7 July 2016, Uppsala, Sweden, to assess and formulate recommendations for nomenclature for amyloid fibril proteins and the clin. classification of the amyloidoses. An amyloid fibril must exhibit affinity for Congo red and with green, yellow or orange birefringence when the Congo red-stained deposits are viewed with polarized light. While congophilia and birefringence remain the gold std. for demonstration of amyloid deposits, new staining and imaging techniques are proving useful. To be included in the nomenclature list, in addn. to congophilia and birefringence, the chem. identity of the protein must be unambiguously characterized by protein sequence anal. when possible. In general, it is insufficient to identify a mutation in the gene of a candidate amyloid protein without confirming the variant changes in the amyloid fibril protein. Each distinct form of amyloidosis is uniquely characterized by the chem. identity of the amyloid fibril protein that deposits in the extracellular spaces of tissues and organs and gives rise to the disease syndrome. The fibril proteins are designated as protein A followed by a suffix that is an abbreviation of the parent or precursor protein name. To date, there are 36 known extracellular fibril proteins in humans, 2 of which are iatrogenic in nature and 9 of which have also been identified in animals. Two newly recognized fibril proteins, AApoCII derived from apolipoprotein CII and AApoCIII derived from apolipoprotein CIII, have been added. AApoCII amyloidosis and AApoCIII amyloidosis are hereditary systemic amyloidoses. Intracellular protein inclusions displaying some of the properties of amyloid, "intracellular amyloid" have been reported. Two proteins which were previously characterized as intracellular inclusions, tau and α-synuclein, are now recognized to form extracellular deposits upon cell death and thus have been included in Table 1 as ATau and AαSyn.
- 3Obici, L.; Merlini, G. AA amyloidosis: Basic knowledge, unmet needs and future treatments. Swiss Med. Wkly. 2012, 142, w13580, DOI: 10.4414/smw.2012.13580
- 4Röcken, C.; Shakespeare, A. Pathology, diagnosis and pathogenesis of AA amyloidosis. Virchows Arch 2002, 440, 111– 122, DOI: 10.1007/s00428-001-0582-9[Crossref], [PubMed], [CAS], Google Scholar4https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BD383is12ktw%253D%253D&md5=578f51b866526454b68445925df76deePathology, diagnosis and pathogenesis of AA amyloidosisRocken Christoph; Shakespeare AnnVirchows Archiv : an international journal of pathology (2002), 440 (2), 111-122 ISSN:0945-6317.Amyloid is defined as a proteinaceous tissue deposit that shows a typical green birefringence in polarised light after staining with Congo red, the presence of non-branching linear fibrils of indefinite length with an approximate diameter of 10-12 nm and a distinct X-ray diffraction pattern consistent with Pauling's model of a cross-beta fibril. Approximately 45% of generalised amyloidoses are secondary or reactive (AA) amyloidosis. Among the causes of AA amyloidosis are rheumatic diseases, idiopathic diseases, inherited diseases, infectious diseases and malignant tumours. Recent decades have provided significant advances in our understanding of the pathology and pathogenesis of AA amyloidosis. Its pathogenesis is multifactorial involving many variables such as primary structure of the precursor protein, acute phase response, the presence of non-fibril proteins (e.g. amyloid P component, apolipoprotein E, glycosaminoglycans, proteoglycans and basement membrane proteins), receptors, lipid metabolism and proteases. Study of the pathogenesis of AA amyloidosis has provided many insights into the nature of conformational diseases, which may help in the understanding of other members of this particularly heterogeneous group of diseases, such as Alzheimer's disease and transmissible spongiform encephalopathies.
- 5Real de Asua, D.; Galvan, J. M.; Filigghedu, M. T.; Trujillo, D.; Costa, R.; Cadinanos, J. Systemic AA amyloidosis: epidemiology, diagnosis, and management. Clin. Epidemiol. 2014, 6, 369, DOI: 10.2147/CLEP.S39981[Crossref], [PubMed], [CAS], Google Scholar5https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BC2M3ms1eisg%253D%253D&md5=fcef6da5b8c6d8f0eda36e6e9fab8d28Systemic AA amyloidosis: epidemiology, diagnosis, and managementReal de Asua Diego; Costa Ramon; Galvan Jose Maria; Filigheddu Maria Teresa; Trujillo Davinia; Cadinanos JulenClinical epidemiology (2014), 6 (), 369-77 ISSN:1179-1349.The term "amyloidosis" encompasses the heterogeneous group of diseases caused by the extracellular deposition of autologous fibrillar proteins. The global incidence of amyloidosis is estimated at five to nine cases per million patient-years. While amyloid light-chain (AL) amyloidosis is more frequent in developed countries, amyloid A (AA) amyloidosis is more common in some European regions and in developing countries. The spectrum of AA amyloidosis has changed in recent decades owing to: an increase in the median age at diagnosis; a percent increase in the frequency of primary AL amyloidosis with respect to the AA type; and a substantial change in the epidemiology of the underlying diseases. Diagnosis of amyloidosis is based on clinical organ involvement and histological evidence of amyloid deposits. Among the many tinctorial characteristics of amyloid deposits, avidity for Congo red and metachromatic birefringence under unidirectional polarized light remain the gold standard. Once the initial diagnosis has been made, the amyloid subtype must be identified and systemic organ involvement evaluated. In this sense, the (123)I-labeled serum amyloid P component scintigraphy is a safe and noninvasive technique that has revolutionized the diagnosis and monitoring of treatment in systemic amyloidosis. It can successfully identify anatomical patterns of amyloid deposition throughout the body and enables not only an initial estimation of prognosis, but also the monitoring of the course of the disease and the response to treatment. Given the etiologic diversity of AA amyloidosis, common therapeutic strategies are scarce. All treatment options should be based upon a greater control of the underlying disease, adequate organ support, and treatment of symptoms. Nevertheless, novel therapeutic strategies targeting the formation of amyloid fibrils and amyloid deposition may generate new expectations for patients with AA amyloidosis.
- 6Nguyen, P. H.; Ramamoorthy, A.; Sahoo, B. R.; Zheng, J.; Faller, P.; Straub, J. E.; Dominguez, L.; Shea, J.-E.; Dokholyan, N. V.; De Simone, A. Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer’s Disease, Parkinson’s Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis. Chem. Rev. 2021, 121, 2545– 2647, DOI: 10.1021/acs.chemrev.0c01122[ACS Full Text
], [CAS], Google Scholar
6https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXjtVGmtb0%253D&md5=242576b480f2ce99fa61e6a2bb2dfcaeAmyloid Oligomers: A joint experimental/computational perspective on Alzheimer's disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral SclerosisNguyen, Phuong H.; Ramamoorthy, Ayyalusamy; Sahoo, Bikash R.; Zheng, Jie; Faller, Peter; Straub, John E.; Dominguez, Laura; Shea, Joan-Emma; Dokholyan, Nikolay V.; De Simone, Alfonso; Ma, Buyong; Nussinov, Ruth; Najafi, Saeed; Ngo, Son Tung; Loquet, Antoine; Chiricotto, Mara; Ganguly, Pritam; McCarty, James; Li, Mai Suan; Hall, Carol; Wang, Yiming; Miller, Yifat; Melchionna, Simone; Habenstein, Birgit; Timr, Stepan; Chen, Jiaxing; Hnath, Brianna; Strodel, Birgit; Kayed, Rakez; Lesne, Sylvain; Wei, Guanghong; Sterpone, Fabio; Doig, Andrew J.; Derreumaux, PhilippeChemical Reviews (Washington, DC, United States) (2021), 121 (4), 2545-2647CODEN: CHREAY; ISSN:0009-2665. (American Chemical Society)A review. Protein misfolding and aggregation is obsd. in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by exptl. and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacol. expts. tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, resp., for many years. - 7Cai, L.; De Beer, M. C.; De Beer, F. C.; Van Der Westhuyzen, D. R. Serum amyloid a is a ligand for scavenger receptor class B type I and inhibits high density lipoprotein binding and selective lipid uptake. J. Biol. Chem. 2005, 280, 2954– 2961, DOI: 10.1074/jbc.M411555200[Crossref], [PubMed], [CAS], Google Scholar7https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXmt1yqtA%253D%253D&md5=296babba7b8568a797fb9fe485a5843bSerum Amyloid A Is a Ligand for Scavenger Receptor Class B Type I and Inhibits High Density Lipoprotein Binding and Selective Lipid UptakeCai, Lei; de Beer, Maria C.; de Beer, Frederick C.; van der Westhuyzen, Deneys R.Journal of Biological Chemistry (2005), 280 (4), 2954-2961CODEN: JBCHA3; ISSN:0021-9258. (American Society for Biochemistry and Molecular Biology)Serum amyloid A is an acute phase protein that is carried in the plasma largely as an apolipoprotein of high d. lipoprotein (HDL). In this study we investigated whether SAA is a ligand for the HDL receptor, scavenger receptor class B type I (SR-BI), and how SAA may influence SR-BI-mediated HDL binding and selective cholesteryl ester uptake. Studies using Chinese hamster ovary cells expressing SR-BI showed that 125I-labeled SAA, both in lipid-free form and in reconstituted HDL particles, functions as a high affinity ligand for SR-BI. SAA also bound with high affinity to the hepatocyte cell line, HepG2. Alexa-labeled SAA was shown by fluorescence confocal microscopy to be internalized by cells in a SR-BI-dependent manner. To assess how SAA assocn. with HDL influences HDL interaction with SR-BI, SAA-contg. HDL was isolated from mice overexpressing SAA through adenoviral gene transfer. SAA presence on HDL had little effect on HDL binding to SR-BI but decreased (30-50%) selective cholesteryl ester uptake. Lipid-free SAA, unlike lipid-free apoA-I, was an effective inhibitor of both SR-BI-dependent binding and selective cholesteryl ester uptake of HDL. We have concluded that SR-BI plays a key role in SAA metab. through its ability to interact with and internalize SAA and, further, that SAA influences HDL cholesterol metab. through its inhibitory effects on SR-BI-mediated selective lipid uptake.
- 8Wang, W.; Khatua, P.; Hansmann, U. H. Cleavage, Downregulation, and Aggregation of Serum Amyloid A. J. Phys. Chem. B 2020, 124, 1009– 1019, DOI: 10.1021/acs.jpcb.9b10843[ACS Full Text
], [CAS], Google Scholar
8https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhtlWnsLs%253D&md5=9e32a66ffbdd75e1b8d3b4a8cc94d18aCleavage, downregulation, and aggregation of serum amyloid AWang, Wenhua; Khatua, Prabir; Hansmann, Ulrich H. E.Journal of Physical Chemistry B (2020), 124 (6), 1009-1019CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)Various diseases cause overexpression of the serum amyloid A (SAA) protein, which in some cases, but not in all cases, leads to amyloidosis as a secondary disease. Response to the overexpression involves dissocn. of the SAA hexamer and subsequent cleavage of the released monomers, most commonly yielding fragments SAA1-76 of the full-sized SAA1-104. We report results from mol. dynamic simulations that probe the role of this cleavage for downregulating the activity and concn. of SAA. We propose a mechanism that relies on two elements. First, the probability to assemble into hexamers is lower for the fragments than it is for the full-sized protein. Second, unlike other fragments, SAA1-76 can switch between two distinct configurations. The first kind is easy to proteolyze (allowing a fast redn. of the SAA concn.) but prone to aggregation, whereas the situation is opposite for the second kind. If the time scale for amyloid formation is longer than the one for proteolysis, the aggregation-prone species dominates. However, if environmental conditions such as low pH increases the risk of amyloid formation, the ensemble shifts toward the more protected form. We speculate that SAA amyloidosis is a failure of this switching mechanism leading to accumulation of the aggregation-prone species and subsequent amyloid formation. - 9Yaşar, F.; Bernhardt, N. A.; Hansmann, U. H. Replica-exchange-with-tunneling for fast exploration of protein landscapes. J. Chem. Phys. 2015, 143, 224102, DOI: 10.1063/1.4936968[Crossref], [PubMed], [CAS], Google Scholar9https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXhvF2ksbjF&md5=f496bc29c5375597099d9b54c87378f4Replica-exchange-with-tunneling for fast exploration of protein landscapesYasar, Fatih; Bernhardt, Nathan A.; Hansmann, Ulrich H. E.Journal of Chemical Physics (2015), 143 (22), 224102/1-224102/6CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)While the use of replica-exchange mol. dynamics in protein simulations has become ubiquitous, its utility is limited in many practical applications. The authors propose to overcome some shortcomings that hold back its use in settings such as multi-scale or explicit solvent simulations by integrating ideas of hybrid MC/MD into the replica-exchange protocol. This Replica-Exchange-with-Tunneling method is tested by simulating the Trp-cage protein, a system often used in mol. biophysics for testing sampling techniques. (c) 2015 American Institute of Physics.
- 10Bernhardt, N. A.; Xi, W.; Wang, W.; Hansmann, U. H. Simulating Protein Fold Switching by Replica Exchange with Tunneling. J. Chem. Theory Comput. 2016, 12, 5656– 5666, DOI: 10.1021/acs.jctc.6b00826[ACS Full Text
], [CAS], Google Scholar
10https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhslWjsLjN&md5=06db8e4a33672ac6a9e26208037590c1Simulating Protein Fold Switching by Replica Exchange with TunnelingBernhardt, Nathan A.; Xi, Wenhui; Wang, Wei; Hansmann, Ulrich H. E.Journal of Chemical Theory and Computation (2016), 12 (11), 5656-5666CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)Recent expts. suggest that an amino acid sequence encodes not only the native fold of a protein but also other forms essential for its function, or important during folding or assocn. These various forms populate a multi-funnel folding and assocn. landscape where mutations, changes in environment, or interaction with other mols. switch between the encoded folds. Here, the authors introduced replica-exchange-with-tunneling as a way to simulate efficiently switching between distinct folds of proteins and protein aggregates. The correctness and efficiency of this approach was demonstrated in a series of simulations covering a wide range of proteins, from a small 11-residue large designed peptide up to 2 56-residue large mutants of the A and B domain of streptococcal protein G. - 11Zhang, H.; Xi, W.; Hansmann, U. H.; Wei, Y. Fibril-Barrel Transitions in Cylindrin Amyloids. J. Chem. Theory Comput. 2017, 13, 3936– 3944, DOI: 10.1021/acs.jctc.7b00383[ACS Full Text
], [CAS], Google Scholar
11https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtFSjtrrO&md5=d144e2f7757c674122513919ccff1640Fibril-Barrel Transitions in Cylindrin AmyloidsZhang, Huiling; Xi, Wenhui; Hansmann, Ulrich H. E.; Wei, YanjieJournal of Chemical Theory and Computation (2017), 13 (8), 3936-3944CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We introduce Replica-Exchange-with-Tunneling (RET) simulations as a tool for studies of the conversion between polymorph amyloids. For the eleven-residue amyloid- forming cylindrin peptide we show that this technique allows for a more efficient sampling of the formation and interconversion between fibril-like and barrel-like assemblies. We describe a protocol for optimized anal. of RET simulations that allows us to propose a mechanism for formation and interconversion between various cylindrin assemblies. Esp., we show that an interchain salt bridge between residues K3 and D7 is crucial for formation of the barrel structure. - 12Bernhardt, N. A.; Hansmann, U. H. Multifunnel Landscape of the Fold-Switching Protein RfaH-CTD. J. Phys. Chem. B 2018, 122, 1600– 1607, DOI: 10.1021/acs.jpcb.7b11352[ACS Full Text
], [CAS], Google Scholar
12https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXnvVCnsQ%253D%253D&md5=2bf1bd93a92bbe0444c16e73cdcd2dc4Multifunnel Landscape of the Fold-Switching Protein RfaH-CTDBernhardt, Nathan A.; Hansmann, Ulrich H. E.Journal of Physical Chemistry B (2018), 122 (5), 1600-1607CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)Proteins such as transcription factor RfaH can change biol. function by switching between distinct 3-dimensional folds. RfaH regulates transcription if the C-terminal domain folds into a double helix bundle, and promotes translation when this domain assumes a β-barrel form. This fold-switch has been also obsd. for the isolated C-terminal domain (RfaH-CTD) and was studied here with a variant of the RET (Replica-Exchange-with-Tunneling) approach recently introduced by us. We used the enhanced sampling properties of this technique to map the free energy landscape of RfaH-CTD and to propose a mechanism for the conversion process. - 13Lu, J.; Yu, Y.; Zhu, I.; Cheng, Y.; Sun, P. D. Structural mechanism of serum amyloid A-mediated inflammatory amyloidosis. Proc. Natl. Acad. Sci. U. S. A. 2014, 111, 5189– 5194, DOI: 10.1073/pnas.1322357111[Crossref], [PubMed], [CAS], Google Scholar13https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXkslGqt7w%253D&md5=cec09402517df4fee69ac5d78acaf3f8Structural mechanism of serum amyloid A-mediated inflammatory amyloidosisLu, Jinghua; Yu, Yadong; Zhu, Iowis; Cheng, Yifan; Sun, Peter D.Proceedings of the National Academy of Sciences of the United States of America (2014), 111 (14), 5189-5194CODEN: PNASA6; ISSN:0027-8424. (National Academy of Sciences)Serum amyloid A (SAA) represents an evolutionarily conserved family of inflammatory acute-phase proteins. It is also a major constituent of secondary amyloidosis. Here, to understand its function and structural transition to amyloid, the authors detd. a structure of human SAA isoform 1.1 (SAA1.1) in 2 crystal forms, representing a prototypic member of the family. Native SAA1.1 exists as a hexamer, with subunits displaying a unique 4-helix bundle fold stabilized by its long C-terminal tail. Structure-based mutational studies revealed 2 pos.-charge clusters, near the center and apex of the hexamer, that were involved in SAA assocn. with heparin. The binding of high-d. lipoprotein involved only the apex region of SAA and could be inhibited by heparin. Peptide amyloid formation assays identified N-terminal helixes 1 and 3 as amyloidogenic peptides of SAA1.1. Both peptides were secluded in the hexameric structure of SAA1.1, suggesting that native SAA is non-pathogenic. Furthermore, dissocn. of the SAA hexamer appeared insufficient to initiate the amyloidogenic transition, and proteolytic cleavage or removal of the C-terminal tail of SAA resulted in the formation of various-sized structural aggregates contg. ∼5-nm regular repeating protofibril-like units. The combined structural and functional studies provided mechanistic insights into the pathogenic contribution of glycosaminoglycan in SAA1.1-mediated AA amyloid formation.
- 14Liberta, F.; Loerch, S.; Rennegarbe, M.; Schierhorn, A.; Westermark, P.; Westermark, G. T.; Hazenberg, B. P.; Grigorieff, N.; Fändrich, M.; Schmidt, M. Cryo-EM fibril structures from systemic AA amyloidosis reveal the species complementarity of pathological amyloids. Nature Commun. 2019, 10, 1104, DOI: 10.1038/s41467-019-09033-z[Crossref], [PubMed], [CAS], Google Scholar14https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BB3cbhs12nug%253D%253D&md5=af5ff643e7e3d8503a088d93c58df447Cryo-EM fibril structures from systemic AA amyloidosis reveal the species complementarity of pathological amyloidsLiberta Falk; Rennegarbe Matthies; Fandrich Marcus; Schmidt Matthias; Loerch Sarah; Grigorieff Nikolaus; Schierhorn Angelika; Westermark Per; Westermark Gunilla T; Hazenberg Bouke P CNature communications (2019), 10 (1), 1104 ISSN:.Systemic AA amyloidosis is a worldwide occurring protein misfolding disease of humans and animals. It arises from the formation of amyloid fibrils from the acute phase protein serum amyloid A. Here, we report the purification and electron cryo-microscopy analysis of amyloid fibrils from a mouse and a human patient with systemic AA amyloidosis. The obtained resolutions are 3.0 ÅA and 2.7 ÅA for the murine and human fibril, respectively. The two fibrils differ in fundamental properties, such as presence of right-hand or left-hand twisted cross-β sheets and overall fold of the fibril proteins. Yet, both proteins adopt highly similar β-arch conformations within the N-terminal ~21 residues. Our data demonstrate the importance of the fibril protein N-terminus for the stability of the analyzed amyloid fibril morphologies and suggest strategies of combating this disease by interfering with specific fibril polymorphs.
- 15Bansal, A.; Schmidt, M.; Rennegarbe, M.; Haupt, C.; Liberta, F.; Stecher, S.; Puscalau-Girtu, I.; Biedermann, A.; Fändrich, M. AA amyloid fibrils from diseased tissue are structurally different from in vitro formed SAA fibrils. Nature Commun. 2021, 12, 1013, DOI: 10.1038/s41467-021-21129-z[Crossref], [PubMed], [CAS], Google Scholar15https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXksV2kt7Y%253D&md5=d38267d13a31b0be0bce42df1582fe08AA amyloid fibrils from diseased tissue are structurally different from in vitro formed SAA fibrilsBansal, Akanksha; Schmidt, Matthias; Rennegarbe, Matthies; Haupt, Christian; Liberta, Falk; Stecher, Sabrina; Puscalau-Girtu, Ioana; Biedermann, Alexander; Faendrich, MarcusNature Communications (2021), 12 (1), 1013CODEN: NCAOBW; ISSN:2041-1723. (Nature Research)Systemic AA amyloidosis is a world-wide occurring protein misfolding disease of humans and animals. It arises from the formation of amyloid fibrils from serum amyloid A (SAA) protein. Using cryo electron microscopy we here show that amyloid fibrils which were purified from AA amyloidotic mice are structurally different from fibrils formed from recombinant SAA protein in vitro. Ex vivo amyloid fibrils consist of fibril proteins that contain more residues within their ordered parts and possess a higher β-sheet content than in vitro fibril proteins. They are also more resistant to proteolysis than their in vitro formed counterparts. These data suggest that pathogenic amyloid fibrils may originate from proteolytic selection, allowing specific fibril morphologies to proliferate and to cause damage to the surrounding tissue.
- 16Fukunishi, H.; Watanabe, O.; Takada, S. On the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure prediction. J. Chem. Phys. 2002, 116, 9058– 9067, DOI: 10.1063/1.1472510[Crossref], [CAS], Google Scholar16https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38XjsFKmsLo%253D&md5=7ac571a5afdd63b0b4b29cfdec06f53bOn the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure predictionFukunishi, Hiroaki; Watanabe, Osamu; Takada, ShojiJournal of Chemical Physics (2002), 116 (20), 9058-9067CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Motivated by the protein structure prediction problem, we develop two variants of the Hamiltonian replica exchange methods (REMs) for efficient configuration sampling, (1) the scaled hydrophobicity REM and (2) the phantom chain REM, and compare their performance with the ordinary REM. We first point out that the ordinary REM has a shortage for the application to large systems such as biomols. and that the Hamiltonian REM, an alternative formulation of the REM, can give a remedy for it. We then propose two examples of the Hamiltonian REM that are suitable for a coarse-grained protein model. (1) The scaled hydrophobicity REM preps. replicas that are characterized by various strengths of hydrophobic interaction. The strongest interaction that mimics aq. soln. environment makes proteins folding, while weakened hydrophobicity unfolds proteins as in org. solvent. Exchange between these environments enables proteins to escape from misfolded traps and accelerate conformational search. This resembles the roles of mol. chaperone that assist proteins to fold in vivo. (2) The phantom chain REM uses replicas that allow various degrees of at. overlaps. By allowing at. overlap in some of replicas, the peptide chain can cross over itself, which can accelerate conformation sampling. Using a coarse-gained model we developed, we compute equil. probability distributions for poly-alanine 16-mer and for a small protein by these REMs and compare the accuracy of the results. We see that the scaled hydrophobicity REM is the most efficient method among the three REMs studied.
- 17Kwak, W.; Hansmann, U. H. Efficient sampling of protein structures by model hopping. Phys. Rev. Lett. 2005, 95, 138102, DOI: 10.1103/PhysRevLett.95.138102[Crossref], [PubMed], [CAS], Google Scholar17https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXhtVelsL7K&md5=8d057f6efdfce47c6ed7a245a4c0da77Efficient Sampling of Protein Structures by Model HoppingKwak, Wooseop; Hansmann, Ulrich H. E.Physical Review Letters (2005), 95 (13), 138102/1-138102/4CODEN: PRLTAO; ISSN:0031-9007. (American Physical Society)A review. We introduce a novel simulation method, model hopping, that enhances sampling of low-energy configurations in complex systems. The approach is illustrated for a protein-folding problem. Thermodn. quantities of proteins with up to 46 residues are evaluated from all-atom simulations with this method.
- 18Hansmann, U. H. Parallel tempering algorithm for conformational studies of biological molecules. Chem. Phys. Lett. 1997, 281, 140– 150, DOI: 10.1016/S0009-2614(97)01198-6[Crossref], [CAS], Google Scholar18https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2sXotFWktr4%253D&md5=ddfed8309757ca88a834a1cf3c80cb08Parallel tempering algorithm for conformational studies of biological moleculesHansmann, Ulrich H. E.Chemical Physics Letters (1997), 281 (1,2,3), 140-150CODEN: CHPLBC; ISSN:0009-2614. (Elsevier Science B.V.)The effectiveness of a new algorithm, parallel tempering, is studied for numerical simulations of biol. mols. These mols. suffer from a rough energy landscape. The resulting slowing down in numerical simulations is overcome by the new method. This is demonstrated by performing simulations with high statistics for one of the simplest peptides, Met-enkephalin. The numerical effectiveness of the new technique was found to be much better than traditional methods and is comparable to sophisticated methods like generalized ensemble techniques.
- 19Best, R. B.; Zhu, X.; Shim, J.; Lopes, P. E.; Mittal, J.; Feig, M.; MacKerell, A. D. Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone ϕ, ψ and side-chain χ1 and χ2 Dihedral Angles. J. Chem. Theory Comput. 2012, 8, 3257– 3273, DOI: 10.1021/ct300400x[ACS Full Text
], [CAS], Google Scholar
19https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38XhtVKqurfP&md5=9a48a0c5770fb1e887c3bb34d45b1354Optimization of the Additive CHARMM All-Atom Protein Force Field Targeting Improved Sampling of the Backbone .vphi., ψ and Side-Chain χ1 and χ2 Dihedral AnglesBest, Robert B.; Zhu, Xiao; Shim, Jihyun; Lopes, Pedro E. M.; Mittal, Jeetain; Feig, Michael; MacKerell, Alexander D.Journal of Chemical Theory and Computation (2012), 8 (9), 3257-3273CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)While the quality of the current CHARMM22/CMAP additive force field for proteins has been demonstrated in a large no. of applications, limitations in the model with respect to the equil. between the sampling of helical and extended conformations in folding simulations have been noted. To overcome this, as well as make other improvements in the model, we present a combination of refinements that should result in enhanced accuracy in simulations of proteins. The common (non-Gly, -Pro) backbone CMAP potential has been refined against exptl. soln. NMR data for weakly structured peptides, resulting in a rebalancing of the energies of the α-helix and extended regions of the Ramachandran map, correcting the α-helical bias of CHARMM22/CMAP. The Gly and Pro CMAPs have been refitted to more accurate quantum-mech. energy surfaces. Side-chain torsion parameters have been optimized by fitting to backbone-dependent quantum-mech. energy surfaces, followed by addnl. empirical optimization targeting NMR scalar couplings for unfolded proteins. A comprehensive validation of the revised force field was then performed against a collection of exptl. data: (i) comparison of simulations of eight proteins in their crystal environments with crystal structures; (ii) comparison with backbone scalar couplings for weakly structured peptides; (iii) comparison with NMR residual dipolar couplings and scalar couplings for both backbone and side-chains in folded proteins; (iv) equil. folding of mini-proteins. The results indicate that the revised CHARMM 36 parameters represent an improved model for modeling and simulation studies of proteins, including studies of protein folding, assembly, and functionally relevant conformational changes. - 20Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 1983, 79, 926– 935, DOI: 10.1063/1.445869[Crossref], [CAS], Google Scholar20https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL3sXksF2htL4%253D&md5=a1161334e381746be8c9b15a5e56f704Comparison of simple potential functions for simulating liquid waterJorgensen, William L.; Chandrasekhar, Jayaraman; Madura, Jeffry D.; Impey, Roger W.; Klein, Michael L.Journal of Chemical Physics (1983), 79 (2), 926-35CODEN: JCPSA6; ISSN:0021-9606.Classical Monte Carlo simulations were carried out for liq. H2O in the NPT ensemble at 25° and 1 atm using 6 of the simpler intermol. potential functions for the dimer. Comparisons were made with exptl. thermodn. and structural data including the neutron diffraction results of Thiessen and Narten (1982). The computed densities and potential energies agree with expt. except for the original Bernal-Fowler model, which yields an 18% overest. of the d. and poor structural results. The discrepancy may be due to the correction terms needed in processing the neutron data or to an effect uniformly neglected in the computations. Comparisons were made for the self-diffusion coeffs. obtained from mol. dynamics simulations.
- 21Pettersen, E. F.; Goddard, T. D.; Huang, C. C.; Couch, G. S.; Greenblatt, D. M.; Meng, E. C.; Ferrin, T. E. UCSF Chimera - A visualization system for exploratory research and analysis. J. Comput. Chem. 2004, 25, 1605– 1612, DOI: 10.1002/jcc.20084[Crossref], [PubMed], [CAS], Google Scholar21https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2cXmvVOhsbs%253D&md5=944b175f440c1ff323705987cf937ee7UCSF Chimera-A visualization system for exploratory research and analysisPettersen, Eric F.; Goddard, Thomas D.; Huang, Conrad C.; Couch, Gregory S.; Greenblatt, Daniel M.; Meng, Elaine C.; Ferrin, Thomas E.Journal of Computational Chemistry (2004), 25 (13), 1605-1612CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)The design, implementation, and capabilities of an extensible visualization system, UCSF Chimera, are discussed. Chimera is segmented into a core that provides basic services and visualization, and extensions that provide most higher level functionality. This architecture ensures that the extension mechanism satisfies the demands of outside developers who wish to incorporate new features. Two unusual extensions are presented: Multiscale, which adds the ability to visualize large-scale mol. assemblies such as viral coats, and Collab., which allows researchers to share a Chimera session interactively despite being at sep. locales. Other extensions include Multalign Viewer, for showing multiple sequence alignments and assocd. structures; ViewDock, for screening docked ligand orientations; Movie, for replaying mol. dynamics trajectories; and Vol. Viewer, for display and anal. of volumetric data. A discussion of the usage of Chimera in real-world situations is given, along with anticipated future directions. Chimera includes full user documentation, is free to academic and nonprofit users, and is available for Microsoft Windows, Linux, Apple Mac OS X, SGI IRIX, and HP Tru64 Unix from http://www.cgl.ucsf.edu/chimera/.
- 22Noel, J. K.; Whitford, P. C.; Sanbonmatsu, K. Y.; Onuchic, J. N. SMOG@ctbp: Simplified deployment of structure-based models in GROMACS. Nucleic Acids Res. 2010, 38, W657, DOI: 10.1093/nar/gkq498[Crossref], [PubMed], [CAS], Google Scholar22https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXotVSqt7s%253D&md5=937e9d055c9a72ef70e2deb3002f3813SMOG@ctbp: simplified deployment of structure-based models in GROMACSNoel, Jeffrey K.; Whitford, Paul C.; Sanbonmatsu, Karissa Y.; Onuchic, Jose N.Nucleic Acids Research (2010), 38 (Web Server), W657-W661CODEN: NARHAD; ISSN:0305-1048. (Oxford University Press)Mol. dynamics simulations with coarse-grained and/or simplified Hamiltonians are an effective means of capturing the functionally important long-time and large-length scale motions of proteins and RNAs. Structure-based Hamiltonians, simplified models developed from the energy landscape theory of protein folding, have become a std. tool for investigating biomol. dynamics. SMOG@ctbp is an effort to simplify the use of structure-based models. The purpose of the web server is two fold. First, the web tool simplifies the process of implementing a well-characterized structure-based model on a state-of-the-art, open source, mol. dynamics package, GROMACS. Second, the tutorial-like format helps speed the learning curve of those unfamiliar with mol. dynamics. A web tool user is able to upload any multi-chain biomol. system consisting of std. RNA, DNA and amino acids in PDB format and receive as output all files necessary to implement the model in GROMACS. Both Cα and all-atom versions of the model are available. SMOG@ctbp resides at http://smog.ucsd.edu.
- 23Zhang, W.; Chen, J. Accelerate sampling in atomistic energy landscapes using topology-based coarse-grained models. J. Chem. Theory Comput. 2014, 10, 918– 923, DOI: 10.1021/ct500031v[ACS Full Text
], [CAS], Google Scholar
23https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhsVOjtro%253D&md5=e8368763a3870dac4a993d744da02ba8Accelerate Sampling in Atomistic Energy Landscapes Using Topology-Based Coarse-Grained ModelsZhang, Weihong; Chen, JianhanJournal of Chemical Theory and Computation (2014), 10 (3), 918-923CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We describe a multiscale enhanced sampling (MSES) method where efficient topol.-based coarse-grained models are coupled with all-atom ones to enhance the sampling of atomistic protein energy landscape. The bias from the coupling is removed by Hamiltonian replica exchange, thus allowing one to benefit simultaneously from faster transitions of coarse-grained modeling and accuracy of atomistic force fields. The method is demonstrated by calcg. the conformational equil. of several small but nontrivial β-hairpins with varied stabilities. - 24Moritsugu, K.; Terada, T.; Kidera, A. Scalable free energy calculation of proteins via multiscale essential sampling. J. Chem. Phys. 2010, 133, 224105, DOI: 10.1063/1.3510519[Crossref], [PubMed], [CAS], Google Scholar24https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXhsFGksbrM&md5=1e8de2879382edb9b892f1cac9933f82Scalable free energy calculation of proteins via multiscale essential samplingMoritsugu, Kei; Terada, Tohru; Kidera, AkinoriJournal of Chemical Physics (2010), 133 (22), 224105/1-224105/6CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A multiscale simulation method, multiscale essential sampling (MSES), is proposed for calcg. free energy surface of proteins in a sizable dimensional space with good scalability. In MSES, the configurational sampling of a full-dimensional model is enhanced by coupling with the accelerated dynamics of the essential degrees of freedom. Applying the Hamiltonian exchange method to MSES can remove the biasing potential from the coupling term, deriving the free energy surface of the essential degrees of freedom. The form of the coupling term ensures good scalability in the Hamiltonian exchange. As a test application, the free energy surface of the folding process of a mini-protein, chignolin, was calcd. in the continuum solvent model. Results agreed with the free energy surface derived from the multi-canonical simulation. Significantly improved scalability with the MSES method was clearly shown in the free energy calcn. of chignolin in explicit solvent, which was achieved without increasing the no. of replicas in the Hamiltonian exchange. (c) 2010 American Institute of Physics.
- 25Hess, B.; Kutzner, C.; Van Der Spoel, D.; Lindahl, E. GROMACS 4: Algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput. 2008, 4, 435– 447, DOI: 10.1021/ct700301q[ACS Full Text
], [CAS], Google Scholar
25https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXhsVSqurc%253D&md5=d53c94901386260221792ea30f151c5fGROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular SimulationHess, Berk; Kutzner, Carsten; van der Spoel, David; Lindahl, ErikJournal of Chemical Theory and Computation (2008), 4 (3), 435-447CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)Mol. simulation is an extremely useful, but computationally very expensive tool for studies of chem. and biomol. systems. Here, we present a new implementation of our mol. simulation toolkit GROMACS which now both achieves extremely high performance on single processors from algorithmic optimizations and hand-coded routines and simultaneously scales very well on parallel machines. The code encompasses a minimal-communication domain decompn. algorithm, full dynamic load balancing, a state-of-the-art parallel constraint solver, and efficient virtual site algorithms that allow removal of hydrogen atom degrees of freedom to enable integration time steps up to 5 fs for atomistic simulations also in parallel. To improve the scaling properties of the common particle mesh Ewald electrostatics algorithms, we have in addn. used a Multiple-Program, Multiple-Data approach, with sep. node domains responsible for direct and reciprocal space interactions. Not only does this combination of algorithms enable extremely long simulations of large systems but also it provides that simulation performance on quite modest nos. of std. cluster nodes. - 26Hess, B. P-LINCS: A parallel linear constraint solver for molecular simulation. J. Chem. Theory Comput. 2008, 4, 116– 122, DOI: 10.1021/ct700200b[ACS Full Text
], [CAS], Google Scholar
26https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXhtlKru7zL&md5=476d5ca2eb25574d44b775996fff7b75P-LINCS: A Parallel Linear Constraint Solver for Molecular SimulationHess, BerkJournal of Chemical Theory and Computation (2008), 4 (1), 116-122CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)By removing the fastest degrees of freedom, constraints allow for an increase of the time step in mol. simulations. In the last decade parallel simulations have become commonplace. However, up till now efficient parallel constraint algorithms have not been used with domain decompn. In this paper the parallel linear constraint solver (P-LINCS) is presented, which allows the constraining of all bonds in macromols. Addnl. the energy conservation properties of (P-)LINCS are assessed in view of improvements in the accuracy of uncoupled angle constraints and integration in single precision. - 27Swope, W. C.; Andersen, H. C.; Berens, P. H.; Wilson, K. R. A computer simulation method for the calculation of equilibrium constants for the formation of physical clusters of molecules: Application to small water clusters. J. Chem. Phys. 1982, 76, 637– 649, DOI: 10.1063/1.442716[Crossref], [CAS], Google Scholar27https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL38XlvVWhtg%253D%253D&md5=696100391d72d68cb6eee764c9691d01A computer simulation method for the calculation of equilibrium constants for the formation of physical clusters of molecules: Application to small water clustersSwope, William C.; Andersen, Hans C.; Berens, Peter H.; Wilson, Kent R.Journal of Chemical Physics (1982), 76 (1), 637-49CODEN: JCPSA6; ISSN:0021-9606.A mol. dynamics computer simulation method is described for calcg. equil. consts. for the formation of phys. clusters of mols. The method is based on Hill's formal theory of phys. clusters. A mol. dynamics calcn. is used to calc. the av. potential energy of a cluster of mols. as a function of temp., and the equil. consts. are calcd. from the integral of the energy with respect to reciprocal temp. The method is illustrated by calcns. of the equil. consts. for the formation of clusters of 2-5 water mols. that interact with each other by an intermol. potential devised by Watts. The method is compared with other procedures for calcg. the thermodn. properties of clusters.
- 28Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 2007, 126, 014101, DOI: 10.1063/1.2408420[Crossref], [PubMed], [CAS], Google Scholar28https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXosVCltg%253D%253D&md5=9c182b57bfc65bca6be23c8c76b4be77Canonical sampling through velocity rescalingBussi, Giovanni; Donadio, Davide; Parrinello, MicheleJournal of Chemical Physics (2007), 126 (1), 014101/1-014101/7CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The authors present a new mol. dynamics algorithm for sampling the canonical distribution. In this approach the velocities of all the particles are rescaled by a properly chosen random factor. The algorithm is formally justified and it is shown that, in spite of its stochastic nature, a quantity can still be defined that remains const. during the evolution. In numerical applications this quantity can be used to measure the accuracy of the sampling. The authors illustrate the properties of this new method on Lennard-Jones and TIP4P water models in the solid and liq. phases. Its performance is excellent and largely independent of the thermostat parameter also with regard to the dynamic properties.
- 29McGibbon, R. T.; Beauchamp, K. A.; Harrigan, M. P.; Klein, C.; Swails, J. M.; Hernández, C. X.; Schwantes, C. R.; Wang, L. P.; Lane, T. J.; Pande, V. S. MDTraj: A Modern Open Library for the Analysis of Molecular Dynamics Trajectories. Biophys. J. 2015, 109, 1528– 1532, DOI: 10.1016/j.bpj.2015.08.015[Crossref], [PubMed], [CAS], Google Scholar29https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXhsVKhu7%252FI&md5=9cbd2138ce200e7d63fda8ae756df43fMDTraj: A Modern Open Library for the Analysis of Molecular Dynamics TrajectoriesMcGibbon, Robert T.; Beauchamp, Kyle A.; Harrigan, Matthew P.; Klein, Christoph; Swails, Jason M.; Hernandez, Carlos X.; Schwantes, Christian R.; Wang, Lee-Ping; Lane, Thomas J.; Pande, Vijay S.Biophysical Journal (2015), 109 (8), 1528-1532CODEN: BIOJAU; ISSN:0006-3495. (Cell Press)As mol. dynamics (MD) simulations continue to evolve into powerful computational tools for studying complex biomol. systems, the necessity of flexible and easy-to-use software tools for the anal. of these simulations is growing. We have developed MDTraj, a modern, lightwt., and fast software package for analyzing MD simulations. MDTraj reads and writes trajectory data in a wide variety of commonly used formats. It provides a large no. of trajectory anal. capabilities including minimal root-mean-square-deviation calcns., secondary structure assignment, and the extn. of common order parameters. The package has a strong focus on interoperability with the wider scientific Python ecosystem, bridging the gap between MD data and the rapidly growing collection of industry-std. statistical anal. and visualization tools in Python. MDTraj is a powerful and user-friendly software package that simplifies the anal. of MD data and connects these datasets with the modern interactive data science software ecosystem in Python.
- 30Dijkstra, E. W. A note on two problems in connexion with graphs. Numer. Math. 1959, 1, 269– 271, DOI: 10.1007/BF01386390
- 31Marcos-Alcalde, I.; Setoain, J.; Mendieta-Moreno, J. I.; Mendieta, J.; Gómez-Puertas, P. MEPSA: minimum energy pathway analysis for energy landscapes. Bioinformatics 2015, 31, 3853– 3855, DOI: 10.1093/bioinformatics/btv453[Crossref], [PubMed], [CAS], Google Scholar31https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28Xht1CitLzJ&md5=aea18fdc9e14b7eda6dfacea20145116MEPSA: minimum energy pathway analysis for energy landscapesMarcos-Alcalde, Inigo; Setoain, Javier; Mendieta-Moreno, Jesus I.; Mendieta, Jesus; Gomez-Puertas, PaulinoBioinformatics (2015), 31 (23), 3853-3855CODEN: BOINFP; ISSN:1367-4803. (Oxford University Press)Summary: From conformational studies to atomistic descriptions of enzymic reactions, potential and free energy landscapes can be used to describe biomol. systems in detail. However, extg. the relevant data of complex 3D energy surfaces can sometimes be laborious. In this article, we present MEPSA (Min. Energy Path Surface Anal.), a cross-platform user friendly tool for the anal. of energy landscapes from a transition state theory perspective. Some of its most relevant features are: identification of all the barriers and min. of the landscape at once, description of maxima edge profiles, detection of the lowest energy path connecting two min. and generation of transition state theory diagrams along these paths. In addn. to a built-in plotting system, MEPSA can save most of the generated data into easily parseable text files, allowing more versatile uses of MEPSA's output such as the generation of mol. dynamics restraints from a calcd. path.
- 32Khatua, P.; Ray, A. J.; Hansmann, U. H. Bifurcated Hydrogen Bonds and the Fold Switching of Lymphotactin. J. Phys. Chem. B 2020, 124, 6555– 6564, DOI: 10.1021/acs.jpcb.0c04565[ACS Full Text
], [CAS], Google Scholar
32https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhtlSisL7F&md5=35d96fe846f57c94f6546bc33602fa47Bifurcated Hydrogen Bonds and the Fold Switching of LymphotactinKhatua, Prabir; Ray, Alan J.; Hansmann, Ulrich H. E.Journal of Physical Chemistry B (2020), 124 (30), 6555-6564CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)Lymphotactin (Ltn) exists under physiol. conditions in an equil. between two interconverting structures with distinct biol. functions. Using replica-exchange-with-tunneling, we study the conversion between the 2-folds. Unlike previously proposed, we find that the fold switching does not require unfolding of lymphotactin but proceeds through a series of intermediates that remain partially structured. This process relies on two bifurcated hydrogen bonds that connect the β2 and β3 strands and ease the transition between the hydrogen bond pattern by which the central three-stranded β-sheet in the two forms differs. - 33Bolhuis, P. G.; Chandler, D.; Dellago, C.; Geissler, P. L. Transition Path Sampling: Throwing Ropes over Rough Mountain Passes, in the Dark. Annu. Rev. Phys. Chem. 2002, 53, 291– 318, DOI: 10.1146/annurev.physchem.53.082301.113146[Crossref], [PubMed], [CAS], Google Scholar33https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38Xks1Ois7g%253D&md5=fbf79b52e751dfd87a73c345b9581898Transition path sampling: throwing ropes over rough mountain passes, in the darkBolhuis, Peter G.; Chandler, David; Dellago, Christoph; Geissler, Phillip L.Annual Review of Physical Chemistry (2002), 53 (), 291-318CODEN: ARPLAP; ISSN:0066-426X. (Annual Reviews Inc.)A review is given of the concepts and methods of transition path sampling. These methods allow computational studies of rare events without requiring prior knowledge of mechanisms, reaction coordinates, and transition states. Based upon a statistical mechanics of trajectory space, they provide a perspective with which time dependent phenomena, even for systems driven far from equil., can be examd. with the same types of importance sampling tools that in the past have been applied so successfully to static equil. properties.
- 34Dellago, C.; Bolhuis, P. G.; Csajka, F. S.; Chandler, D. Transition Path Sampling and the Calculation of Rate Constants. J. Chem. Phys. 1998, 108, 1964– 1977, DOI: 10.1063/1.475562[Crossref], [CAS], Google Scholar34https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1cXkvFChsQ%253D%253D&md5=4bb2f7b8316dbb4526be81bd24083814Transition path sampling and the calculation of rate constantsDellago, Christoph; Bolhuis, Peter G.; Csajka, Felix S.; Chandler, DavidJournal of Chemical Physics (1998), 108 (5), 1964-1977CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We have developed a method to study transition pathways for rare events in complex systems. The method can be used to det. rate consts. for transitions between stable states by turning the calcn. of reactive flux correlation functions into the computation of an isomorphic reversible work. In contrast to previous dynamical approaches, the method relies neither on prior knowledge nor on explicit specification of transition states. Rather, it provides an importance sampling from which transition states can be characterized statistically. A simple model is analyzed to illustrate the methodol.
- 35Dellago, C.; Bolhuis, P. G.; Chandler, D. Efficient Transition Path Sampling: Application to Lennard-Jones Cluster Rearrangements. J. Chem. Phys. 1998, 108, 9236– 9245, DOI: 10.1063/1.476378[Crossref], [CAS], Google Scholar35https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1cXjtFSitb8%253D&md5=0d02dabdf6c5811d0bf11847d8949c7eEfficient transition path sampling: Application to Lennard-Jones cluster rearrangementsDellago, Christoph; Bolhuis, Peter G.; Chandler, DavidJournal of Chemical Physics (1998), 108 (22), 9236-9245CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We develop an efficient Monte-Carlo algorithm to sample an ensemble of stochastic transition paths between stable states. In our description, paths are represented by chains of states linked by Markovian transition probabilities. Rate consts. and mechanisms characterizing the transition may be detd. from the path ensemble. We have previously devised several algorithms for sampling the path ensemble. For these algorithms, the numerical effort scales with the square of the path length. In the new simulation scheme, the required computation scales linearly with the length of the transition path. This improved efficiency allows the calcn. of rate consts. in complex mol. systems. As an example, we study rearrangement processes in a cluster consisting of seven Lennard-Jones particles in two dimensions. Using a quenching technique we are able to identify the relevant transition mechanisms and to locate the related transition states. We furthermore calc. transition rate consts. for various isomerization processes.
- 36Weinan, E.; Ren, W.; Vanden-Eijnden, E. String Method for the Study of Rare Events. Phys. Rev. B 2002, 66, 052301, DOI: 10.1103/PhysRevB.66.052301
- 37Maragliano, L.; Fischer, A.; Vanden-Eijnden, E.; Ciccotti, G. String Method in Collective Variables: Minimum Free Energy Paths and Isocommittor Surfaces. J. Chem. Phys. 2006, 125, 024106, DOI: 10.1063/1.2212942[Crossref], [PubMed], [CAS], Google Scholar37https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28Xnt1Khurk%253D&md5=b3a7b0b167df6980ac6c8e7bb36b16fcString method in collective variables: Minimum free energy paths and isocommittor surfacesMaragliano, Luca; Fischer, Alexander; Vanden-Eijnden, Eric; Ciccotti, GiovanniJournal of Chemical Physics (2006), 125 (2), 024106/1-024106/15CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A computational technique is proposed which combines the string method with a sampling technique to det. min. free energy paths. The technique only requires to compute the mean force and another conditional expectation locally along the string, and therefore can be applied even if the no. of collective variables kept in the free energy calcn. is large. This is in contrast with other free energy sampling techniques which aim at mapping the full free energy landscape and whose cost increases exponentially with the no. of collective variables kept in the free energy. Provided that the no. of collective variables is large enough, the new technique captures the mechanism of transition in that it allows to det. the committor function for the reaction and, in particular, the transition state region. The new technique is illustrated on the example of alanine dipeptide, in which we compute the min. free energy path for the isomerization transition using either two or four dihedral angles as collective variables. It is shown that the mechanism of transition can be captured using the four dihedral angles, but it cannot be captured using only two of them.
- 38Unarta, I. C.; Zhu, L.; Tse, C. K. M.; Cheung, P. P.-H.; Yu, J.; Huang, X. Molecular Mechanisms of RNA Polymerase II Transcription Elongation Elucidated by Kinetic Network Models. Curr. Opin. Struct. Biol. 2018, 49, 54– 62, DOI: 10.1016/j.sbi.2018.01.002[Crossref], [PubMed], [CAS], Google Scholar38https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhsVensbg%253D&md5=efc857e4ccb759523af7030715a53d1fMolecular mechanisms of RNA polymerase II transcription elongation elucidated by kinetic network modelsUnarta, Ilona Christy; Zhu, Lizhe; Tse, Carmen Ka Man; Cheung, Peter Pak-Hang; Yu, Jin; Huang, XuhuiCurrent Opinion in Structural Biology (2018), 49 (), 54-62CODEN: COSBEF; ISSN:0959-440X. (Elsevier Ltd.)Transcription elongation cycle (TEC) of RNA polymerase II (Pol II) is a process of adding a nucleoside triphosphate to the growing mRNA chain. Due to the long timescale events in Pol II TEC, an advanced computational technique, such as Markov State Model (MSM), is needed to provide atomistic mechanism and reaction rates. The combination of MSM and exptl. results can be used to build a kinetic network model (KNM) of the whole TEC. This review provides a brief protocol to build MSM and KNM of the whole TEC, along with the latest findings of MSM and other computational studies of Pol II TEC. Lastly, we offer a perspective on potentially using a sequence dependent KNM to predict genome-wide transcription error.
- 39Zhu, L.; Sheong, F. K.; Cao, S.; Liu, S.; Unarta, I. C.; Huang, X. TAPS: A Traveling-Salesman Based Automated Path Searching Method for Functional Conformational Changes of Biological Macromolecules. J. Chem. Phys. 2019, 150, 124105, DOI: 10.1063/1.5082633[Crossref], [PubMed], [CAS], Google Scholar39https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXmtVClt74%253D&md5=18ac3c2a1ef78dff772057ec9ad3aff9TAPS: A traveling-salesman based automated path searching method for functional conformational changes of biological macromoleculesZhu, Lizhe; Sheong, Fu Kit; Cao, Siqin; Liu, Song; Unarta, Ilona C.; Huang, XuhuiJournal of Chemical Physics (2019), 150 (12), 124105/1-124105/10CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Locating the min. free energy paths (MFEPs) between two conformational states is among the most important tasks of biomol. simulations. For example, knowledge of the MFEP is crit. for focusing the effort of unbiased simulations that were used for the construction of Markov state models to the biol. relevant regions of the system. Typically, existing path searching methods perform local sampling around the path nodes in a pre-selected collective variable (CV) space to allow a gradual downhill evolution of the path toward the MFEP. Despite the wide application of such a strategy, the gradual path evolution and the non-trivial a priori choice of CVs are also limiting its overall efficiency and automation. Non-local perpendicular sampling can be pursued to accelerate the search, provided that all nodes are reordered thereafter via a traveling-salesman scheme. Moreover, path-CVs can be computed on-the-fly and used as a coordinate system, minimizing the necessary prior knowledge about the system. The authors' traveling-salesman based automated path searching method achieves a 5-8 times speedup over the string method with swarms-of-trajectories for two peptide systems in vacuum and soln., making it a promising method for obtaining initial pathways when studying functional conformational changes between a pair of structures. (c) 2019 American Institute of Physics.
- 40Schrödinger, L. The PyMOL Molecular Graphics System , Version 1.8; 2015.Google ScholarThere is no corresponding record for this reference.
- 41Westermark, G. T.; Engström, U.; Westermark, P. The N-terminal segment of protein AA determines its fibrillogenic property. Biochem. Biophys. Res. Commun. 1992, 182, 27– 33, DOI: 10.1016/S0006-291X(05)80107-X[Crossref], [PubMed], [CAS], Google Scholar41https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK38Xht1egur4%253D&md5=80fe3a456eac7d41a32920d051c0c60dThe N-terminal segment of protein AA determines its fibrillogenic propertyWestermark, Gunilla T.; Engstroem, Ulla; Westermark, PerBiochemical and Biophysical Research Communications (1992), 182 (1), 27-33CODEN: BBRCA9; ISSN:0006-291X.The amyloid fibril protein AA consists of a varying long N-terminal part of the precursor protein serum AA. By using synthetic peptides corresponding to human and murine protein AA segments and cyanogen bromide fragments of human protein AA, evidence is shown that the amyloidogenic part of the mol. is the first 10-15 amino acid long segment. Amino acid substitutions in this part of the mol. may explain why only one of the two mouse SAA isoforms is amyloidogenic.
- 42Yassine, H. N.; Trenchevska, O.; He, H.; Borges, C. R.; Nedelkov, D.; Mack, W.; Kono, N.; Koska, J.; Reaven, P. D.; Nelson, R. W. Serum Amyloid A Truncations in Type 2 Diabetes Mellitus. PLoS One 2015, 10, e0115320, DOI: 10.1371/journal.pone.0115320
- 43Wang, W.; Hansmann, U. H. Stability of Human Serum Amyloid A Fibrils. J. Phys. Chem. B 2020, 124, 10708– 10717, DOI: 10.1021/acs.jpcb.0c08280[ACS Full Text
], [CAS], Google Scholar
43https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXitlagtLbF&md5=b9dce5f2c90c561d7a9c9b50d1829c39Stability of Human Serum Amyloid A FibrilsWang, Wenhua; Hansmann, Ulrich H. E.Journal of Physical Chemistry B (2020), 124 (47), 10708-10717CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)In systemic amyloidosis, serum amyloid A (SAA) fibril deposits cause widespread damages to tissues and organs that eventually may lead to death. A therapeutically intervention therefore has either to dissolve these fibrils or inhibit their formation. However, only recently has the human SAA fibril structure been resolved at a resoln. that is sufficient for development of drug candidates. Here, we use mol. dynamic simulations to probe the factors that modulate the stability of this fibril model. Our simulations suggest that fibril formation starts with the stacking of two misfolded monomers into metastable dimers, with the stacking depending on the N-terminal amyloidogenic regions of different chains forming anchors. The resulting dimers pack in a second step into a 2-fold two-layer tetramer that is stable enough to nucleate fibril formation. The stability of the initial dimers is enhanced under acidic conditions by a strong salt bridge and side-chain hydrogen bond network in the C-terminal cavity (residues 23-51) but is not affected by the presence of the disordered C-terminal tail. - 44Wang, W.; Xi, W.; Hansmann, U. H. Stability of the N-Terminal Helix and Its Role in Amyloid Formation of Serum Amyloid A. ACS Omega 2018, 3, 16184– 16190, DOI: 10.1021/acsomega.8b02377[ACS Full Text
], [CAS], Google Scholar
44https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXitlaktbbL&md5=59f3032458272f0eff6faa7ca0e70f80Stability of the N-Terminal Helix and Its Role in Amyloid Formation of Serum Amyloid AWang, Wenhua; Xi, Wenhui; Hansmann, Ulrich H. E.ACS Omega (2018), 3 (11), 16184-16190CODEN: ACSODF; ISSN:2470-1343. (American Chemical Society)Colonic amyloidosis is the result of overexpression of the serum amyloid A protein (SAA) in inflammatory bowel disease or colon cancer. Crucial for amyloid formation are the 1st ten N-terminal residues, which in the crystal structure are part of a 27-residue long helix. Here we studied this 27-residue N-terminal region of SAA by a multi-exchange variant of replica exchange mol. dynamics (ME-REMD) simulations. An ensemble of configurations was obsd., dominated by 3 motifs: the single helix of the crystal structure, a helix-turn-helix configuration, and such where residues 14-27 were part of a helix but the 1st 13 residues formed an extended and disordered segment that was prone to aggregation. A single point mutation E9A shifted the equil. to the latter motif indicating the importance of interactions involving this residue for the stability of the SAA protein. - 45Frame, N. M.; Kumanan, M.; Wales, T. E.; Bandara, A.; Fändrich, M.; Straub, J. E.; Engen, J. R.; Gursky, O. Structural Basis for Lipid Binding and Function by an Evolutionarily Conserved Protein, Serum Amyloid A. J. Mol. Bio. 2020, 432, 1978– 1995, DOI: 10.1016/j.jmb.2020.01.029[Crossref], [PubMed], [CAS], Google Scholar45https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXivFyntro%253D&md5=892c2c04204d39ad65267b17c04c3fadStructural Basis for Lipid Binding and Function by an Evolutionarily Conserved Protein, Serum Amyloid AFrame, Nicholas M.; Kumanan, Meera; Wales, Thomas E.; Bandara, Asanga; Fandrich, Marcus; Straub, John E.; Engen, John R.; Gursky, OlgaJournal of Molecular Biology (2020), 432 (7), 1978-1995CODEN: JMOBAK; ISSN:0022-2836. (Elsevier Ltd.)Serum amyloid A (SAA) is a plasma protein that transports lipids during inflammation. To explore SAA soln. conformations and lipid-binding mechanism, we used hydrogen-deuterium exchange mass spectrometry, lipoprotein reconstitution, amino acid sequence anal., and mol. dynamics simulations. Soln. conformations of lipid-bound and lipid-free mSAA1 at pH∼7.4 agreed in details with the crystal structures but also showed important differences. The results revealed that amphipathic α-helixes h1 and h3 comprise a lipid-binding site that is partially pre-formed in soln., is stabilized upon binding lipids, and shows lipid-induced folding of h3. This site sequesters apolar ligands via a concave hydrophobic surface in SAA oligomers. The largely disordered/dynamic C-terminal region is conjectured to mediate the promiscuous binding of other ligands. The h1-h2 linker region is predicted to form an unexpected β-hairpin that may represent an early amyloidogenic intermediate. The results help establish structural underpinnings for understanding SAA interactions with its key functional ligands, its evolutional conservation, and its transition to amyloid.