1
|
Loureiro RJS, Faísca PFN. The Early Phase of β2-Microglobulin Aggregation: Perspectives From Molecular Simulations. Front Mol Biosci 2020; 7:578433. [PMID: 33134317 PMCID: PMC7550760 DOI: 10.3389/fmolb.2020.578433] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2020] [Accepted: 09/08/2020] [Indexed: 11/24/2022] Open
Abstract
Protein β2-microglobulin is the causing agent of two amyloidosis, dialysis related amyloidosis (DRA), affecting the bones and cartilages of individuals with chronic renal failure undergoing long-term hemodialysis, and a systemic amyloidosis, found in one French family, which impairs visceral organs. The protein’s small size and its biomedical significance attracted the attention of theoretical scientists, and there are now several studies addressing its aggregation mechanism in the context of molecular simulations. Here, we review the early phase of β2-microglobulin aggregation, by focusing on the identification and structural characterization of monomers with the ability to trigger aggregation, and initial small oligomers (dimers, tetramers, hexamers etc.) formed in the so-called nucleation phase. We focus our analysis on results from molecular simulations and integrate our views with those coming from in vitro experiments to provide a broader perspective of this interesting field of research. We also outline directions for future computer simulation studies.
Collapse
Affiliation(s)
- Rui J S Loureiro
- Faculty of Sciences, BioISI - Biosystems & Integrative Sciences Institute, University of Lisboa, Lisbon, Portugal
| | - Patrícia F N Faísca
- Faculty of Sciences, BioISI - Biosystems & Integrative Sciences Institute, University of Lisboa, Lisbon, Portugal.,Department of Physics, Faculty of Sciences, University of Lisboa, Lisbon, Portugal
| |
Collapse
|
2
|
Narang SS, Shuaib S, Goyal D, Goyal B. In silico-guided identification of potential inhibitors against β2m aggregation in dialysis-related amyloidosis. J Biomol Struct Dyn 2019; 38:3927-3941. [DOI: 10.1080/07391102.2019.1668852] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Simranjeet Singh Narang
- Department of Chemistry, Faculty of Basic and Applied Sciences, Sri Guru Granth Sahib World University, Fatehgarh Sahib, Punjab, India
| | - Suniba Shuaib
- Department of Chemistry, Faculty of Basic and Applied Sciences, Sri Guru Granth Sahib World University, Fatehgarh Sahib, Punjab, India
| | - Deepti Goyal
- Department of Chemistry, Faculty of Basic and Applied Sciences, Sri Guru Granth Sahib World University, Fatehgarh Sahib, Punjab, India
| | - Bhupesh Goyal
- School of Chemistry and Biochemistry, Thapar Institute of Engineering & Technology, Patiala, Punjab, India
| |
Collapse
|
3
|
Narang SS, Shuaib S, Goyal B. Molecular insights into the inhibitory mechanism of rifamycin SV against β 2 –microglobulin aggregation: A molecular dynamics simulation study. Int J Biol Macromol 2017; 102:1025-1034. [DOI: 10.1016/j.ijbiomac.2017.04.086] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2017] [Revised: 04/24/2017] [Accepted: 04/24/2017] [Indexed: 01/30/2023]
|
4
|
Han W, Wan CK, Wu YD. PACE Force Field for Protein Simulations. 2. Folding Simulations of Peptides. J Chem Theory Comput 2015; 6:3390-402. [PMID: 26617093 DOI: 10.1021/ct100313a] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
We present the application of our recently developed PACE force field to the folding of peptides. These peptides include α-helical (AK17 and Fs), β-sheet (GB1m2 and Trpzip2), and mixed helical/coil (Trp-cage) peptides. With replica exchange molecular dynamics (REMD), our force field can fold the five peptides into their native structures while maintaining their stabilities reasonably well. Our force field is also able to capture important thermodynamic features of the five peptides that have been observed in previous experimental and computational studies, such as different preferences for a helix-turn-helix topology for AK17 and Fs, the relative contribution of four hydrophobic side chains of GB1p to the stability of β-hairpin, and the distinct role of a hydrogen bond involving Trp-Hε and a D9/R16 salt bridge in stabilizing the Trp-cage native structure. Furthermore, multiple folding and unfolding events are observed in our microsecond-long normal MD simulations of AK17, Trpzip2, and Trp-cage. These simulations provide mechanistic information such as a "zip-out" pathway of the folding mechanism of Trpzip2 and the folding times of AK17 and Trp-cage, which are estimated to be about 51 ± 43 ns and 270 ± 110 ns, respectively. A 600 ns simulation of the peptides can be completed within one day. These features of our force field are potentially applicable to the study of thermodynamics and kinetics of real protein systems.
Collapse
Affiliation(s)
- Wei Han
- Department of Chemistry, The Hong Kong University of Science & Technology, Clear Water Bay, Kowloon, Hong Kong, China, School of Chemical Biology and Biotechnology, Laboratory of Chemical Genomics, Peking University Shenzhen Graduate School, Shenzhen, China, and College of Chemistry, Peking University, Beijing, China
| | - Cheuk-Kin Wan
- Department of Chemistry, The Hong Kong University of Science & Technology, Clear Water Bay, Kowloon, Hong Kong, China, School of Chemical Biology and Biotechnology, Laboratory of Chemical Genomics, Peking University Shenzhen Graduate School, Shenzhen, China, and College of Chemistry, Peking University, Beijing, China
| | - Yun-Dong Wu
- Department of Chemistry, The Hong Kong University of Science & Technology, Clear Water Bay, Kowloon, Hong Kong, China, School of Chemical Biology and Biotechnology, Laboratory of Chemical Genomics, Peking University Shenzhen Graduate School, Shenzhen, China, and College of Chemistry, Peking University, Beijing, China
| |
Collapse
|
5
|
Zhou P, Deng L, Wang Y, Lu JR, Xu H. Different nanostructures caused by competition of intra- and inter-β-sheet interactions in hierarchical self-assembly of short peptides. J Colloid Interface Sci 2015; 464:219-28. [PMID: 26619132 DOI: 10.1016/j.jcis.2015.11.030] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2015] [Revised: 11/12/2015] [Accepted: 11/13/2015] [Indexed: 10/22/2022]
Abstract
To understand how molecular interactions lead to the self-assembly of twisted, helical and flat nanoribbons, we have compared the hierarchical self-assembly processes of three selected octapeptides with the same amino acid composition but different sequences by both experiments and molecular dynamics (MD) simulations. KE-F8 (NH2-KEFFFFKE-CONH2) and EK-F8 (NH2-KEFFFFEK-CONH2) have the same distribution of hydrophobic residues and only differ by swapping the positive and negative charged residues at their C-terminals, while KFE-8 (NH2-KFEFKFEF-CONH2) differs from KE-F8 and EK-F8 by having all hydrophobic and charged residues evenly distributed. MD simulations indicated that the competition between electrostatic and hydrophobic interactions at the molecular level results in different initial packing modes: KE-F8 monomers form completely matched anti-parallel β-sheets, EK-F8 monomers align with one residue shifting, and KFE-8 monomers pack β-sheets with two heterogeneous surfaces, consistent with previously suggested models. Driven by inter-strand and inter-sheet interactions, further growth of these molecular templates leads to larger oligomers with different twisting and stacking degrees, which are structurally consistent with the experimentally observed self-assembled morphologies. Further MD simulations showed that the competition between intra-β-sheet and inter-β-sheet interactions is responsible for the different twisting and stacking degrees of β-sheets and the subsequent formation of different nanostructures (twisted ribbons for KE-F8, helical ribbons/tubes for EK-F8 and flat ribbons for KFE-8). This study thus provided an important mechanistic insight into the fine tuning of molecular packing and interactions via peptide sequence variation leading to controllable self-assembly of twisted, helical and flat nanostructures.
Collapse
Affiliation(s)
- Peng Zhou
- State Key Laboratory of Heavy Oil Processing and Centre for Bioengineering and Biotechnology, China University of Petroleum (East China), 66 Changjiang West Road, Qingdao 266580, China
| | - Li Deng
- State Key Laboratory of Heavy Oil Processing and Centre for Bioengineering and Biotechnology, China University of Petroleum (East China), 66 Changjiang West Road, Qingdao 266580, China
| | - Yanting Wang
- State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, 55 East Zhongguancun Road, P.O. Box 2735, Beijing 100190, China.
| | - Jian R Lu
- Biological Physics Group, School of Physics and Astronomy, The University of Manchester, Manchester M13 9PL, United Kingdom.
| | - Hai Xu
- State Key Laboratory of Heavy Oil Processing and Centre for Bioengineering and Biotechnology, China University of Petroleum (East China), 66 Changjiang West Road, Qingdao 266580, China.
| |
Collapse
|
6
|
Nasica-Labouze J, Nguyen PH, Sterpone F, Berthoumieu O, Buchete NV, Coté S, De Simone A, Doig AJ, Faller P, Garcia A, Laio A, Li MS, Melchionna S, Mousseau N, Mu Y, Paravastu A, Pasquali S, Rosenman DJ, Strodel B, Tarus B, Viles JH, Zhang T, Wang C, Derreumaux P. Amyloid β Protein and Alzheimer's Disease: When Computer Simulations Complement Experimental Studies. Chem Rev 2015; 115:3518-63. [PMID: 25789869 DOI: 10.1021/cr500638n] [Citation(s) in RCA: 478] [Impact Index Per Article: 53.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Jessica Nasica-Labouze
- †Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique (IBPC), UPR9080 CNRS, Université Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Phuong H Nguyen
- †Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique (IBPC), UPR9080 CNRS, Université Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Fabio Sterpone
- †Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique (IBPC), UPR9080 CNRS, Université Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Olivia Berthoumieu
- ‡LCC (Laboratoire de Chimie de Coordination), CNRS, Université de Toulouse, Université Paul Sabatier (UPS), Institut National Polytechnique de Toulouse (INPT), 205 route de Narbonne, BP 44099, Toulouse F-31077 Cedex 4, France
| | | | - Sébastien Coté
- ∥Département de Physique and Groupe de recherche sur les protéines membranaires (GEPROM), Université de Montréal, C.P. 6128, succursale Centre-ville, Montréal, Québec H3C 3T5, Canada
| | - Alfonso De Simone
- ⊥Department of Life Sciences, Imperial College London, London SW7 2AZ, United Kingdom
| | - Andrew J Doig
- #Manchester Institute of Biotechnology, University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom
| | - Peter Faller
- ‡LCC (Laboratoire de Chimie de Coordination), CNRS, Université de Toulouse, Université Paul Sabatier (UPS), Institut National Polytechnique de Toulouse (INPT), 205 route de Narbonne, BP 44099, Toulouse F-31077 Cedex 4, France
| | | | - Alessandro Laio
- ○The International School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy
| | - Mai Suan Li
- ◆Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-668 Warsaw, Poland.,¶Institute for Computational Science and Technology, SBI Building, Quang Trung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam
| | - Simone Melchionna
- ⬠Instituto Processi Chimico-Fisici, CNR-IPCF, Consiglio Nazionale delle Ricerche, 00185 Roma, Italy
| | | | - Yuguang Mu
- ▲School of Biological Sciences, Nanyang Technological University, 60 Nanyang Drive, 637551 Singapore
| | - Anant Paravastu
- ⊕National High Magnetic Field Laboratory, 1800 East Paul Dirac Drive, Tallahassee, Florida 32310, United States
| | - Samuela Pasquali
- †Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique (IBPC), UPR9080 CNRS, Université Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | | | - Birgit Strodel
- △Institute of Complex Systems: Structural Biochemistry (ICS-6), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Bogdan Tarus
- †Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique (IBPC), UPR9080 CNRS, Université Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - John H Viles
- ▼School of Biological and Chemical Sciences, Queen Mary University of London, London E1 4NS, United Kingdom
| | - Tong Zhang
- †Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique (IBPC), UPR9080 CNRS, Université Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France.,▲School of Biological Sciences, Nanyang Technological University, 60 Nanyang Drive, 637551 Singapore
| | | | - Philippe Derreumaux
- †Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique (IBPC), UPR9080 CNRS, Université Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France.,□Institut Universitaire de France, 75005 Paris, France
| |
Collapse
|
7
|
Nguyen PH, Tarus B, Derreumaux P. Familial Alzheimer A2 V Mutation Reduces the Intrinsic Disorder and Completely Changes the Free Energy Landscape of the Aβ1–28 Monomer. J Phys Chem B 2014; 118:501-10. [DOI: 10.1021/jp4115404] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Phuong H. Nguyen
- Laboratoire
de Biochimie Théorique, UPR 9080 CNRS, Université Paris Diderot, Sorbonne Paris Cité, IBPC, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Bogdan Tarus
- Laboratoire
de Biochimie Théorique, UPR 9080 CNRS, Université Paris Diderot, Sorbonne Paris Cité, IBPC, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Philippe Derreumaux
- Laboratoire
de Biochimie Théorique, UPR 9080 CNRS, Université Paris Diderot, Sorbonne Paris Cité, IBPC, 13 rue Pierre et Marie Curie, 75005 Paris, France
- Institut Universitaire de France, IUF, 103 Boulevard Saint-Michel, 75005 Paris, France
| |
Collapse
|
8
|
Nguyen PH, Derreumaux P. Conformational ensemble and polymorphism of the all-atom Alzheimer's Aβ(37-42) amyloid peptide oligomers. J Phys Chem B 2013; 117:5831-40. [PMID: 23581814 DOI: 10.1021/jp401563n] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Although the Aβ(37-42) peptide has two opposite terminal charges, counterintuitively its current fibril amyloid structure reveals in register parallel β-strands, as formed by the full length Aβ peptide. In this study, we carried out a replica exchange molecular dynamics simulation of 16 all-atom Aβ(37-42) peptides in explicit water starting from randomized and dispersed chains. The extensive conformational sampling (48 replicas, 460 ns/replica) with a total simulation time of 23 μs allows us to obtain a full picture on the equilibrium conformational distribution of oligomers and β-sheet sizes and gain some insights into the oligomerization process at 300 K. At the peptide concentration of 12 mM, self-assembly is described by the condensation-polymerization mechanism with conversion from micelle-like to high β-sheet structures. At equilibrium, the oligomer distribution consists of large aggregates and free monomers, representing 70% and 25% of all species, respectively. Though the formation/dissociation of β-strand is high, the population of 4-5 fully parallel β-strands, consistent with the arrangement in the current fibril, is marginal and that of 4-5 fully antiparallel β-strands, consistent with amyloid polymorphism, is non-negligible. However, the system adopts essentially mixed parallel/antiparallel β-strands. This indicates that a system of 16 Aβ(37-42) chains in explicit solvent still does not form more stable species that will irreversibly grow to a fibril, independently of polymorphism. Our results also suggest that the Aβ(37-42) fibril may display packing polymorphism with antiparallel β-strands, in addition to the experimentally observed in register parallel β-strands.
Collapse
Affiliation(s)
- Phuong H Nguyen
- Laboratoire de Biochimie Théorique, UPR 9080, CNRS, Université Denis Diderot, Sorbonne Paris Cité, IBPC, 13 rue Pierre et Marie Curie, 75005 Paris, France.
| | | |
Collapse
|
9
|
Fogolari F, Corazza A, Varini N, Rotter M, Gumral D, Codutti L, Rennella E, Viglino P, Bellotti V, Esposito G. Molecular dynamics simulation of β₂-microglobulin in denaturing and stabilizing conditions. Proteins 2010; 79:986-1001. [PMID: 21287627 DOI: 10.1002/prot.22940] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2010] [Revised: 10/22/2010] [Accepted: 11/02/2010] [Indexed: 11/11/2022]
Abstract
β₂-Microglobulin has been a model system for the study of fibril formation for 20 years. The experimental study of β₂-microglobulin structure, dynamics, and thermodynamics in solution, at atomic detail, along the pathway leading to fibril formation is difficult because the onset of disorder and aggregation prevents signal resolution in Nuclear Magnetic Resonance experiments. Moreover, it is difficult to characterize conformers in exchange equilibrium. To gain insight (at atomic level) on processes for which experimental information is available at molecular or supramolecular level, molecular dynamics simulations have been widely used in the last decade. Here, we use molecular dynamics to address three key aspects of β₂-microglobulin, which are known to be relevant to amyloid formation: (1) 60 ns molecular dynamics simulations of β₂-microglobulin in trifluoroethanol and in conditions mimicking low pH are used to study the behavior of the protein in environmental conditions that are able to trigger amyloid formation; (2) adaptive biasing force molecular dynamics simulation is used to force cis-trans isomerization at Proline 32 and to calculate the relative free energy in the folded and unfolded state. The native-like trans-conformer (known as intermediate 2 and determining the slow phase of refolding), is simulated for 10 ns, detailing the possible link between cis-trans isomerization and conformational disorder; (3) molecular dynamics simulation of highly concentrated doxycycline (a molecule able to suppress fibril formation) in the presence of β₂-microglobulin provides details of the binding modes of the drug and a rationale for its effect.
Collapse
Affiliation(s)
- Federico Fogolari
- Dipartimento di Scienze e Tecnologie Biomediche, Universita' di Udine, Piazzale Kolbe 4, 33100 Udine, Italy.
| | | | | | | | | | | | | | | | | | | |
Collapse
|
10
|
Knecht V. Model Amyloid Peptide B18 Monomer and Dimer Studied by Replica Exchange Molecular Dynamics Simulations. J Phys Chem B 2010; 114:12701-7. [DOI: 10.1021/jp1048698] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Volker Knecht
- Max Planck Institute of Colloids and Interfaces, Science Park Golm, 14424 Potsdam, Germany
| |
Collapse
|
11
|
Molecular dynamics simulations of the minor ampullate spidroin modular amino acid sequence from Parawixia bistriatra: insights into silk tertiary structure and fibre formation. J Mol Model 2010; 17:1183-9. [DOI: 10.1007/s00894-010-0823-4] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2010] [Accepted: 07/26/2010] [Indexed: 10/19/2022]
|
12
|
Thukral L, Smith JC, Daidone I. Common folding mechanism of a beta-hairpin peptide via non-native turn formation revealed by unbiased molecular dynamics simulations. J Am Chem Soc 2010; 131:18147-52. [PMID: 19919102 DOI: 10.1021/ja9064365] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The folding of a 15-residue beta-hairpin peptide (Peptide 1) is characterized using multiple unbiased, atomistic molecular dynamics (MD) simulations. Fifteen independent MD trajectories, each 2.5 micros-long for a total of 37.5 micros, are performed of the peptide in explicit solvent, at room temperature, and without the use of enhanced sampling techniques. The computed folding time of 1-1.5 micros obtained from the simulations is in good agreement with experiment [Xu, Y.; et al. J. Am. Chem. Soc. 2003, 125, 15388-15394]. A common folding mechanism is observed, in which the turn is always found to be the major determinant in initiating the folding process, followed by cooperative formation of the interstrand hydrogen bonds and the side-chain packing. Furthermore, direct transition to the folded state from fully unstructured conformations does not take place. Instead, the peptide is always observed to form partially structured conformations involving a non-native (ESYI) turn from which the native (NPDG) turn forms, triggering the folding to the beta-hairpin.
Collapse
Affiliation(s)
- Lipi Thukral
- Interdisciplinary Center for Scientific Computing, University of Heidelberg, Im Neuenheimer Feld 368, 69120 Heidelberg, Germany
| | | | | |
Collapse
|
13
|
Mustata M, Capone R, Jang H, Arce FT, Ramachandran S, Lal R, Nussinov R. K3 fragment of amyloidogenic beta(2)-microglobulin forms ion channels: implication for dialysis related amyloidosis. J Am Chem Soc 2010; 131:14938-45. [PMID: 19824733 DOI: 10.1021/ja9049299] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Beta(2)-microglobulin (beta(2)m) amyloid deposits are linked to dialysis-related amyloidosis (DRA) in hemodialysis patients. The mechanism by which beta(2)m causes DRA is not understood. It is also unclear whether only the full-length beta(2)m induces pathophysiology or if proteolytic fragments are sufficient for inducing this effect. Ser20-Lys41 (K3) is a digestion fragment of full-length beta(2)m. Solid state NMR (ssNMR) combined with X-ray diffraction and atomic force microscopy (AFM) revealed the characteristic oligomeric amyloid conformation of the U-turn beta-strand-turn-beta-strand motif stacked in parallel and stabilized by intermolecular interactions also shown by Abeta(9-40)/Abeta(17-42) and the CA150 WW domain. Here we use the K3 U-turn atomic coordinates and molecular dynamic (MD) simulations to model K3 channels in the membrane. Consistent with previous AFM imaging of other amyloids that show channel-like structures in the membrane, in the simulations K3 also forms ion channels with 3-6 loosely attached mobile subunits. We carry out AFM, single channel electrical recording, and fluorescence imaging experiments. AFM images display 3D ion channel topography with shapes, morphologies, and dimensions consistent with the theoretical model. Electrical conductance measurements indicate multiple single channel conductances, suggesting that various K3 oligomer sizes can constitute the channel structure. Fluorescence measurements in kidney cells show channel-mediated cell calcium uptake. These results suggest that the beta(2)m-induced DRA can be mediated by ion channels formed by its K3 fragment. Because the beta-strand-turn-beta-strand motif appears to be a universal amyloid feature, its ability to form ion channels further suggests that the motif may play a generic role in toxicity.
Collapse
Affiliation(s)
- Mirela Mustata
- Center for Nanomedicine and Department of Medicine, University of Chicago, Chicago, Illinois 60637, USA
| | | | | | | | | | | | | |
Collapse
|
14
|
Chebaro Y, Dong X, Laghaei R, Derreumaux P, Mousseau N. Replica exchange molecular dynamics simulations of coarse-grained proteins in implicit solvent. J Phys Chem B 2009; 113:267-74. [PMID: 19067549 DOI: 10.1021/jp805309e] [Citation(s) in RCA: 66] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Current approaches aimed at determining the free energy surface of all-atom medium-size proteins in explicit solvent are slow and are not sufficient to converge to equilibrium properties. To ensure a proper sampling of the configurational space, it is preferable to use reduced representations such as implicit solvent and/or coarse-grained protein models, which are much lighter computationally. Each model must be verified, however, to ensure that it can recover experimental structures and thermodynamics. Here we test the coarse-grained implicit solvent OPEP model with replica exchange molecular dynamics (REMD) on six peptides ranging in length from 10 to 28 residues: two alanine-based peptides, the second beta-hairpin from protein G, the Trp-cage and zinc-finger motif, and a dimer of a coiled coil peptide. We show that REMD-OPEP recovers the proper thermodynamics of the systems studied, with accurate structural description of the beta-hairpin and Trp-cage peptides (within 1-2 A from experiments). The light computational burden of REMD-OPEP, which enables us to generate many hundred nanoseconds at each temperature and fully assess convergence to equilibrium ensemble, opens the door to the determination of the free energy surface of larger proteins and assemblies.
Collapse
Affiliation(s)
- Yassmine Chebaro
- Institut de Biologie Physico-Chimique et Universite Paris 7 Denis Diderot, 75005 Paris, France
| | | | | | | | | |
Collapse
|
15
|
Dong X, Chen W, Mousseau N, Derreumaux P. Energy landscapes of the monomer and dimer of the Alzheimer's peptide Abeta(1-28). J Chem Phys 2008; 128:125108. [PMID: 18376983 DOI: 10.1063/1.2890033] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
The cytotoxicity of Alzheimer's disease has been linked to the self-assembly of the 4042 amino acid of the amyloid-beta (Abeta) peptide into oligomers. To understand the assembly process, it is important to characterize the very first steps of aggregation at an atomic level of detail. Here, we focus on the N-terminal fragment 1-28, known to form fibrils in vitro. Circular dichroism and NMR experiments indicate that the monomer of Abeta(1-28) is alpha-helical in a membranelike environment and random coil in aqueous solution. Using the activation-relaxation technique coupled with the OPEP coarse grained force field, we determine the structures of the monomer and of the dimer of Abeta(1-28). In agreement with experiments, we find that the monomer is predominantly random coil in character, but displays a non-negligible beta-strand probability in the N-terminal region. Dimerization impacts the structure of each chain and leads to an ensemble of intertwined conformations with little beta-strand content in the region Leu17-Ala21. All these structural characteristics are inconsistent with the amyloid fibril structure and indicate that the dimer has to undergo significant rearrangement en route to fibril formation.
Collapse
Affiliation(s)
- Xiao Dong
- Département de Physique and Regroupement Québécois sur les Matériaux de Pointe, Université de Montréal, C.P. 6128, Succursale Centre-Ville, Montréal, Québec H3C 3J7, Canada
| | | | | | | |
Collapse
|
16
|
The beta-strand-loop-beta-strand conformation is marginally populated in beta2-microglobulin (20-41) peptide in solution as revealed by replica exchange molecular dynamics simulations. Biophys J 2008; 95:510-7. [PMID: 18408040 DOI: 10.1529/biophysj.107.125054] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Solid-state NMR study shows that the 22-residue K3 peptide (Ser(20)-Lys(41)) from beta(2)-microglobulin (beta(2)m) adopts a beta-strand-loop-beta-strand conformation in its fibril state. Residue Pro(32) has a trans conformation in the fibril state of the peptide, while it adopts a cis conformation in the native state of full-length beta(2)m. To get insights into the structural properties of the K3 peptide, and determine whether the strand-loop-strand conformation is encoded at the monomeric level, we run all-atom explicit solvent replica exchange molecular dynamics on both the cis and trans variants. Our simulations show that the conformational space of the trans- and cis-K3 peptides is very different, with 1% of the sampled conformations in common at room temperature. In addition, both variants display only 0.3-0.5% of the conformations with beta-strand-loop-beta-strand character. This finding, compared to results on the Alzheimer's Abeta peptide, suggests that the biases toward aggregation leading to the beta-strand-loop-beta-strand conformation in fibrils are peptide-dependent.
Collapse
|