51
|
Nguyen TH, Nguyen PH, Ngo ST, Derreumaux P. Effect of Cholesterol Molecules on Aβ1-42 Wild-Type and Mutants Trimers. Molecules 2022; 27:molecules27041395. [PMID: 35209177 PMCID: PMC8879133 DOI: 10.3390/molecules27041395] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2022] [Revised: 02/14/2022] [Accepted: 02/15/2022] [Indexed: 12/04/2022] Open
Abstract
Alzheimer’s disease displays aggregates of the amyloid-beta (Aβ) peptide in the brain, and there is increasing evidence that cholesterol may contribute to the pathogenesis of the disease. Though many experimental and theoretical studies have focused on the interactions of Aβ oligomers with membrane models containing cholesterol, an understanding of the effect of free cholesterol on small Aβ42 oligomers is not fully established. To address this question, we report on replica exchange with a solute tempering simulation of an Aβ42 trimer with cholesterol and compare it with a previous replica exchange molecular dynamics simulation. We show that the binding hot spots of cholesterol are rather complex, involving hydrophobic residues L17–F20 and L30–M35 with a non-negligible contribution of loop residues D22–K28 and N-terminus residues. We also examine the effects of cholesterol on the trimers of the disease-causing A21G and disease-protective A2T mutations by molecular dynamics simulations. We show that these two mutations moderately impact cholesterol-binding modes. In our REST2 simulations, we find that cholesterol is rarely inserted into aggregates but rather attached as dimers and trimers at the surface of Aβ42 oligomers. We propose that cholesterol acts as a glue to speed up the formation of larger aggregates; this provides a mechanistic link between cholesterol and Alzheimer’s disease.
Collapse
Affiliation(s)
- Trung Hai Nguyen
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City, Vietnam; (T.H.N.); (S.T.N.)
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Phuong H. Nguyen
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS, Université de Paris, 13 rue Pierre et Marie Curie, 75005 Paris, France;
- Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
| | - Son Tung Ngo
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City, Vietnam; (T.H.N.); (S.T.N.)
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS, Université de Paris, 13 rue Pierre et Marie Curie, 75005 Paris, France;
- Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
- Institut Universitaire de France (IUF), 75005 Paris, France
- Correspondence:
| |
Collapse
|
52
|
Chávez-García C, Hénin J, Karttunen M. Multiscale Computational Study of the Conformation of the Full-Length Intrinsically Disordered Protein MeCP2. J Chem Inf Model 2022; 62:958-970. [PMID: 35130441 DOI: 10.1021/acs.jcim.1c01354] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The malfunction of the methyl-CpG binding protein 2 (MeCP2) is associated with the Rett syndrome, one of the most common causes of cognitive impairment in females. MeCP2 is an intrinsically disordered protein (IDP), making its experimental characterization a challenge. There is currently no structure available for the full-length MeCP2 in any of the databases, and only the structure of its MBD domain has been solved. We used this structure to build a full-length model of MeCP2 by completing the rest of the protein via ab initio modeling. Using a combination of all-atom and coarse-grained simulations, we characterized its structure and dynamics as well as the conformational space sampled by the ID and transcriptional repression domain (TRD) domains in the absence of the rest of the protein. The present work is the first computational study of the full-length protein. Two main conformations were sampled in the coarse-grained simulations: a globular structure similar to the one observed in the all-atom force field and a two-globule conformation. Our all-atom model is in good agreement with the available experimental data, predicting amino acid W104 to be buried, amino acids R111 and R133 to be solvent-accessible, and having a 4.1% α-helix content, compared to the 4% found experimentally. Finally, we compared the model predicted by AlphaFold to our Modeller model. The model was not stable in water and underwent further folding. Together, these simulations provide a detailed (if perhaps incomplete) conformational ensemble of the full-length MeCP2, which is compatible with experimental data and can be the basis of further studies, e.g., on mutants of the protein or its interactions with its biological partners.
Collapse
Affiliation(s)
- Cecilia Chávez-García
- Department of Chemistry, the University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 5B7, Canada.,The Centre of Advanced Materials and Biomaterials Research, the University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 5B7, Canada
| | - Jérôme Hénin
- Laboratoire de Biochimie Théorique UPR 9080, CNRS and Université de Paris, Paris 75005, France
| | - Mikko Karttunen
- Department of Chemistry, the University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 5B7, Canada.,The Centre of Advanced Materials and Biomaterials Research, the University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 5B7, Canada.,Department of Physics and Astronomy, the University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 3K7, Canada
| |
Collapse
|
53
|
Rahman MU, Song K, Da LT, Chen HF. Early aggregation mechanism of Aβ 16-22 revealed by Markov state models. Int J Biol Macromol 2022; 204:606-616. [PMID: 35134456 DOI: 10.1016/j.ijbiomac.2022.02.001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2021] [Revised: 01/24/2022] [Accepted: 02/01/2022] [Indexed: 12/19/2022]
Abstract
Aβ16-22 is believed to have critical role in early aggregation of full length amyloids that are associated with the Alzheimer's disease and can aggregate to form amyloid fibrils. However, the early aggregation mechanism is still unsolved. Here, multiple long-term molecular dynamics simulations combining with Markov state model were used to probe the early oligomerization mechanism of Aβ16-22 peptides. The identified dimeric form adopted either globular random-coil or extended β-strand like conformations. The observed dimers of these variants shared many overall conformational characteristics but differed in several aspects at detailed level. In all cases, the most common type of secondary structure was intermolecular antiparallel β-sheets. The inter-state transitions were very frequent ranges from few to hundred nanoseconds. More strikingly, those states which contain fraction of β secondary structure and significant amount of extended coiled structures, therefore exposed to the solvent, were majorly participated in aggregation. The assembly of low-energy dimers, in which the peptides form antiparallel β sheets, occurred by multiple pathways with the formation of an obligatory intermediates. We proposed that these states might facilitate the Aβ16-22 aggregation through a significant component of the conformational selection mechanism, because they might increase the aggregates population by promoting the inter-chain hydrophobic and the hydrogen bond contacts. The formation of early stage antiparallel β sheet structures is critical for oligomerization, and at the same time provided a flat geometry to seed the ordered β-strand packing of the fibrils. Our findings hint at reorganization of this part of the molecule as a potentially critical step in Aβ aggregation and will insight into early oligomerization for large β amyloids.
Collapse
Affiliation(s)
- Mueed Ur Rahman
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Kaiyuan Song
- Key Laboratory of System Biomedicine (Ministry of Education), Shanghai Center for Systems Biomedicine, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Lin-Tai Da
- Key Laboratory of System Biomedicine (Ministry of Education), Shanghai Center for Systems Biomedicine, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China; Shanghai Center for Bioinformation Technology, Shanghai, 200235, China.
| |
Collapse
|
54
|
Jalali S, Yang Y, Mahmoudinobar F, Singh SM, Nilsson BL, Dias C. Using all-atom simulations in explicit solvent to study aggregation of amphipathic peptides into amyloid-like fibrils. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.118283] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
55
|
Sahoo BR, Souders CL, Watanabe-Nakayama T, Deng Z, Linton H, Suladze S, Ivanova MI, Reif B, Ando T, Martyniuk CJ, Ramamoorthy A. Conformational Tuning of Amylin by Charged Styrene-Maleic-Acid Copolymers. J Mol Biol 2022; 434:167385. [PMID: 34883118 PMCID: PMC8752516 DOI: 10.1016/j.jmb.2021.167385] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2021] [Revised: 11/25/2021] [Accepted: 11/30/2021] [Indexed: 02/01/2023]
Abstract
Human amylin forms structurally heterogeneous amyloids that have been linked to type-2 diabetes. Thus, understanding the molecular interactions governing amylin aggregation can provide mechanistic insights in its pathogenic formation. Here, we demonstrate that fibril formation of amylin is altered by synthetic amphipathic copolymer derivatives of the styrene-maleic-acid (SMAQA and SMAEA). High-speed AFM is used to follow the real-time aggregation of amylin by observing the rapid formation of de novo globular oligomers and arrestment of fibrillation by the positively-charged SMAQA. We also observed an accelerated fibril formation in the presence of the negatively-charged SMAEA. These findings were further validated by fluorescence, SOFAST-HMQC, DOSY and STD NMR experiments. Conformational analysis by CD and FT-IR revealed that the SMA copolymers modulate the conformation of amylin aggregates. While the species formed with SMAQA are α-helical, the ones formed with SMAEA are rich in β-sheet structure. The interacting interfaces between SMAEA or SMAQA and amylin are mapped by NMR and microseconds all-atom MD simulation. SMAEA displayed π-π interaction with Phe23, electrostatic π-cation interaction with His18 and hydrophobic packing with Ala13 and Val17; whereas SMAQA showed a selective interaction with amylin's C terminus (residues 31-37) that belongs to one of the two β-sheet regions (residues 14-19 and 31-36) involved in amylin fibrillation. Toxicity analysis showed both SMA copolymers to be non-toxic in vitro and the amylin species formed with the copolymers showed minimal deformity to zebrafish embryos. Together, this study demonstrates that chemical tools, such as copolymers, can be used to modulate amylin aggregation, alter the conformation of species.
Collapse
Affiliation(s)
- Bikash R Sahoo
- Biophysics Program, The University of Michigan, Ann Arbor, MI 48109, USA; Department of Chemistry, Biomedical Engineering, and Macromolecular Science and Engineering, The University of Michigan, Ann Arbor, MI 48109, USA
| | | | | | - Zhou Deng
- Biophysics Program, The University of Michigan, Ann Arbor, MI 48109, USA
| | - Hunter Linton
- Department of Neurology, University of Michigan, Ann Arbor, MI 48109-2200, USA
| | - Saba Suladze
- Department of Chemistry, Technische Universität München, Garching 85748, Germany
| | - Magdalena I Ivanova
- Biophysics Program, The University of Michigan, Ann Arbor, MI 48109, USA; Department of Neurology, University of Michigan, Ann Arbor, MI 48109-2200, USA
| | - Bernd Reif
- Department of Chemistry, Technische Universität München, Garching 85748, Germany
| | - Toshio Ando
- Bio-AFM Frontier Research Center, Kanazawa University, Kanazawa 9201192, Japan
| | - Christopher J Martyniuk
- Department of Chemistry, Biomedical Engineering, and Macromolecular Science and Engineering, The University of Michigan, Ann Arbor, MI 48109, USA
| | - Ayyalusamy Ramamoorthy
- Biophysics Program, The University of Michigan, Ann Arbor, MI 48109, USA; Department of Chemistry, Biomedical Engineering, and Macromolecular Science and Engineering, The University of Michigan, Ann Arbor, MI 48109, USA.
| |
Collapse
|
56
|
Roy R, Paul S. hIAPP-Amyloid-Core Derived d-Peptide Prevents hIAPP Aggregation and Destabilizes Its Protofibrils. J Phys Chem B 2022; 126:822-839. [DOI: 10.1021/acs.jpcb.1c10395] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Rituparna Roy
- Department of Chemistry, Indian Institute of Technology, Guwahati, Assam, India, 781039
| | - Sandip Paul
- Department of Chemistry, Indian Institute of Technology, Guwahati, Assam, India, 781039
| |
Collapse
|
57
|
Piskorz T, de Vries AH, van Esch JH. How the Choice of Force-Field Affects the Stability and Self-Assembly Process of Supramolecular CTA Fibers. J Chem Theory Comput 2022; 18:431-440. [PMID: 34812627 PMCID: PMC8757428 DOI: 10.1021/acs.jctc.1c00257] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2021] [Indexed: 01/21/2023]
Abstract
In recent years, computational methods have become an essential element of studies focusing on the self-assembly process. Although they provide unique insights, they face challenges, from which two are the most often mentioned in the literature: the temporal and spatial scale of the self-assembly. A less often mentioned issue, but not less important, is the choice of the force-field. The repetitive nature of the supramolecular structure results in many similar interactions. Consequently, even a small deviation in these interactions can lead to significant energy differences in the whole structure. However, studies comparing different force-fields for self-assembling systems are scarce. In this article, we compare molecular dynamics simulations for trifold hydrogen-bonded fibers performed with different force-fields, namely GROMOS, CHARMM General Force Field (CGenFF), CHARMM Drude, General Amber Force-Field (GAFF), Martini, and polarized Martini. Briefly, we tested the force-fields by simulating: (i) spontaneous self-assembly (none form a fiber within 500 ns), (ii) stability of the fiber (observed for CHARMM Drude, GAFF, MartiniP), (iii) dimerization (observed for GROMOS, GAFF, and MartiniP), and (iv) oligomerization (observed for CHARMM Drude and MartiniP). This system shows that knowledge of the force-field behavior regarding interactions in oligomer and larger self-assembled structures is crucial for designing efficient simulation protocols for self-assembling systems.
Collapse
Affiliation(s)
- Tomasz
K. Piskorz
- Department
of Chemical Engineering, Delft University
of Technology, van der Maasweg 9, Delft, 2629 HZ, The Netherlands
| | - A. H. de Vries
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Jan H. van Esch
- Department
of Chemical Engineering, Delft University
of Technology, van der Maasweg 9, Delft, 2629 HZ, The Netherlands
| |
Collapse
|
58
|
Nguyen PH, Tufféry P, Derreumaux P. Dynamics of Amyloid Formation from Simplified Representation to Atomistic Simulations. Methods Mol Biol 2022; 2405:95-113. [PMID: 35298810 DOI: 10.1007/978-1-0716-1855-4_5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Amyloid fibril formation is an intrinsic property of short peptides, non-disease proteins, and proteins associated with neurodegenerative diseases. Aggregates of the Aβ and tau proteins, the α-synuclein protein, and the prion protein are observed in the brain of Alzheimer's, Parkinson's, and prion disease patients, respectively. Due to the transient short-range and long-range interactions of all species and their high aggregation propensities, the conformational ensemble of these devastating proteins, the exception being for the monomeric prion protein, remains elusive by standard structural biology methods in bulk solution and in lipid membranes. To overcome these limitations, an increasing number of simulations using different sampling methods and protein models have been performed. In this chapter, we first review our main contributions to the field of amyloid protein simulations aimed at understanding the early aggregation steps of short linear amyloid peptides, the conformational ensemble of the Aβ40/42 dimers in bulk solution, and the stability of Aβ aggregates in lipid membrane models. Then we focus on our studies on the interactions of amyloid peptides/inhibitors to prevent aggregation, and long amyloid sequences, including new results on a monomeric tau construct.
Collapse
Affiliation(s)
- Phuong Hoang Nguyen
- Laboratoire de Biochimie Théorique, CNRS, Université de Paris, UPR 9080, Paris, France
- Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | - Pierre Tufféry
- Université de Paris, BFA, UMR 8251, CNRS, ERL U1133, Inserm, RPBS, Paris, France
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, CNRS, Université de Paris, UPR 9080, Paris, France.
- Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, Paris, France.
| |
Collapse
|
59
|
Fatafta H, Kav B, Bundschuh BF, Loschwitz J, Strodel B. Disorder-to-order transition of the amyloid-β peptide upon lipid binding. Biophys Chem 2021; 280:106700. [PMID: 34784548 DOI: 10.1016/j.bpc.2021.106700] [Citation(s) in RCA: 55] [Impact Index Per Article: 18.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2021] [Accepted: 10/15/2021] [Indexed: 12/18/2022]
Abstract
There is mounting evidence that Alzheimer's disease progression and severity are linked to neuronal membrane damage caused by aggregates of the amyloid-β (Aβ) peptide. However, the detailed mechanism behind the membrane damage is not well understood yet. Recently, the lipid-chaperone hypothesis has been put forward, based on which the formation of complexes between Aβ and free lipids enables an easy insertion of Aβ into membranes. In order to test this hypothesis, we performed numerous all-atom molecular dynamics simulations. We studied the complex formation between individual lipids, considering both POPC and DPPC, and Aβ and examined whether the resulting complexes would be able to insert into lipid membranes. Complex formation at a one-to-one ratio was readily observed, yet with minimal effects on Aβ's characteristics. Most importantly, the peptide remains largely disordered in 1:1 complexes, and the complex does not insert into the membrane; instead, it is adsorbed to the membrane surface. The results change considerably once Aβ forms a complex with a POPC cluster composed of three lipid molecules. The hydrophobic interactions between Aβ and the lipid tails cause the peptide to fold into either a helical or a β-sheet structure. These observations provide atomic insight into the disorder-to-order transition that is needed for membrane insertion or amyloid aggregation to proceed.
Collapse
Affiliation(s)
- Hebah Fatafta
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, 52425 Jülich, Germany
| | - Batuhan Kav
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, 52425 Jülich, Germany
| | - Bastian F Bundschuh
- Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Universitütstrasse 1, 40225 Düsseldorf, Germany
| | - Jennifer Loschwitz
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, 52425 Jülich, Germany; Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Universitütstrasse 1, 40225 Düsseldorf, Germany
| | - Birgit Strodel
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, 52425 Jülich, Germany; Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Universitütstrasse 1, 40225 Düsseldorf, Germany.
| |
Collapse
|
60
|
From molecular dynamics to quantum mechanics of misfolded proteins and amyloid-like macroaggregates applied to neurodegenerative diseases. J Mol Graph Model 2021; 110:108046. [PMID: 34736057 DOI: 10.1016/j.jmgm.2021.108046] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2021] [Revised: 09/14/2021] [Accepted: 10/05/2021] [Indexed: 11/24/2022]
Abstract
A misfolded protein compared with its native state lacks its biological function resulting in cell dysregulations and often death. Outdated hypotheses on protein folding must be revised: More realistic molecular models, focusing not only on classical molecular dynamics (MD) but also on ab initio quantum mechanics (QM) at the molecular orbitals (MOs) scale, which is not experimentally achievable, are presented to improve our understanding of the thermodynamics of the protein-protein interactions leading to misfolding and neurodegenerative diseases for future drug design. Protein misfolding is characterized by the formation of highly reactive beta-sheets oligomers leading to fibrillar macroscopic aggregates, which are studied with the models given herein that can be useful for the development of new immunotherapies against the Alzheimer's disease and prion, e.g. The example of the prion - an intrinsically disordered protein - is studied, but the models can be generalized to other misfolding diseases. The binding free energy and interactions in a complex of a misfolded prion with a native prion are first analyzed by MD and compared to a complex of two native conformers. A conversion of residues to toxic beta-sheets is observed in the optimized misfolded complex. Then, QM is used to compute, with a much better accuracy than that of MD, the binding free energy of the hydrophobic binding site, responsible of the aggregation, between the bound misfolded and native conformers in the misfolded complex. The latter quantity is significantly negative, so that aggregation is strong and fast. The frontier MOs from QM are used for docking to determine how the first repetitive beta-sheets building blocks of the nanofibrils can be assembled from initial cleaved complexes of the native and misfolded proteins. Successive aggregation of multiple monomers leads to an amyloid-like nanofibril that grows along a principal elongation direction, as also observed experimentally.
Collapse
|
61
|
Man VH, He X, Gao J, Wang J. Effects of All-Atom Molecular Mechanics Force Fields on Amyloid Peptide Assembly: The Case of PHF6 Peptide of Tau Protein. J Chem Theory Comput 2021; 17:6458-6471. [PMID: 34491058 DOI: 10.1021/acs.jctc.1c00028] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Molecular dynamics (MD) simulations play a vital role in revealing the mechanism of amyloid aggregation that is crucial to the therapeutic agent development for Alzheimer's Disease. However, the accuracy of MD simulation results strongly depends on the force field employed. In our previous benchmark for 17 all-atom force fields on modeling of amyloid aggregation using the Aβ16-22 dimer, we showed that AMBER14SB and CHARMM36m are suitable force fields for amyloid aggregation simulation, while GROMOS54a7 and OPLSAA are not good for the task. In this work, we continue assessing the applicability of atomistic force fields on amyloid aggregation using the VQIVYK (PHF6) peptide which is essential for tau-protein aggregation. Although, both Aβ16-22 and PHF6 peptides formed fibrils in vitro, the PHF6 fibrils are parallel β-sheets, while the Aβ16-22 fibrils are antiparallel β-sheets. We performed an all-atom large-scale MD simulation in explicit water on the PHF6 dimer and octa-peptides systems using five mainstream force fields, including AMBER99SB-disp, AMBER14SB, CHARMM36m, GROMOS54a7, and OPLSAA. The accumulated simulation time is 0.2 ms. Our result showed that the β-sheet structures of PHF6 peptides sampled by AMBER99SB-disp, AMBER14SB, GROMOS54a7, and OPLSAA are in favor of the antiparallel β-sheets, while the dominant type of β-sheet structures is parallel β-sheet by using CHARMM36m. Among the five force fields, CHARMM36m provides the strongest CH-π interaction that was observed in an NMR study. The comparison between our results and experimental observation indicates that CHARMM36m achieved the best performance on modeling the aggregation of PHF6 peptides. In summary, CHARMM36m is currently the most suitable force field for studying the aggregation of both amyloid-β and Tau through MD simulations.
Collapse
Affiliation(s)
- Viet Hoang Man
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Jie Gao
- Department of Neuroscience, The Ohio State University Wexner Medical Center, Columbus, Ohio 43210, United States
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| |
Collapse
|
62
|
Boopathi S, Poma AB, Garduño-Juárez R. An Overview of Several Inhibitors for Alzheimer's Disease: Characterization and Failure. Int J Mol Sci 2021; 22:10798. [PMID: 34639140 PMCID: PMC8509255 DOI: 10.3390/ijms221910798] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2021] [Revised: 10/01/2021] [Accepted: 10/03/2021] [Indexed: 01/04/2023] Open
Abstract
Amyloid beta (Aβ) oligomers are the most neurotoxic aggregates causing neuronal death and cognitive damage. A detailed elucidation of the aggregation pathways from oligomers to fibril formation is crucial to develop therapeutic strategies for Alzheimer's disease (AD). Although experimental techniques rely on the measure of time- and space-average properties, they face severe difficulties in the investigation of Aβ peptide aggregation due to their intrinsically disorder character. Computer simulation is a tool that allows tracing the molecular motion of molecules; hence it complements Aβ experiments, as it allows to explore the binding mechanism between metal ions and Aβ oligomers close to the cellular membrane at the atomic resolution. In this context, integrated studies of experiments and computer simulations can assist in mapping the complete pathways of aggregation and toxicity of Aβ peptides. Aβ oligomers are disordered proteins, and due to a rapid exploration of their intrinsic conformational space in real-time, they are challenging therapeutic targets. Therefore, no good drug candidate could have been identified for clinical use. Our previous investigations identified two small molecules, M30 (2-Octahydroisoquinolin-2(1H)-ylethanamine) and Gabapentin, capable of Aβ binding and inhibiting molecular aggregation, synaptotoxicity, intracellular calcium signaling, cellular toxicity and memory losses induced by Aβ. Thus, we recommend these molecules as novel candidates to assist anti-AD drug discovery in the near future. This review discusses the most recent research investigations about the Aβ dynamics in water, close contact with cell membranes, and several therapeutic strategies to remove plaque formation.
Collapse
Affiliation(s)
- Subramanian Boopathi
- Instituto de Ciencias Físicas, Universidad Nacional Autónoma de México, Cuernavaca 62210, Mexico;
| | - Adolfo B. Poma
- Department of Biosystems and Soft Matter, Institute of Fundamental Technological Research Polish Academy of Science, Pawińskiego 5B, 02-106 Warsaw, Poland
- International Center for Research on Innovative Biobased Materials (ICRI-BioM)—International Research Agenda, Lodz University of Technology, Zeromskiego 116, 90-924 Lodz, Poland;
| | - Ramón Garduño-Juárez
- Instituto de Ciencias Físicas, Universidad Nacional Autónoma de México, Cuernavaca 62210, Mexico;
| |
Collapse
|
63
|
Li Y, Peng X. Comparison of the force fields on monomeric and fibrillar PHF6 of tau protein. Biophys Chem 2021; 277:106631. [PMID: 34116358 DOI: 10.1016/j.bpc.2021.106631] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2021] [Revised: 05/18/2021] [Accepted: 06/01/2021] [Indexed: 11/19/2022]
Abstract
The hexapeptide 306VQIVYK311 (PHF6) plays an important role in the aggregation of Tau protein, which is a hallmark of the Alzheimer's disease (AD). In this article, we systematically compare the effects of eight popular all-atom force fields on the monomeric and fibrillar PHF6 in the molecular dynamics (MD) simulations, which could be helpful in the computer-aided drug design against PHF6. We show that the fibrillar PHF6 prefers β-strand-like structures in all the force fields while the monomer has different structural preferences depending on the force fields. The interactions for stabilizing the fibril are further investigated. In the end, according to the interactions revealed by NMR and the stability of the fibril in the literature, we benchmark the force fields.
Collapse
Affiliation(s)
- Yanchun Li
- Center for Quantum Technology Research, Key Laboratory of Advanced Optoelectronic Quantum Architecture and Measurements (MOE), School of Physics, Beijing Institute of Technology, Beijing 100081, China
| | - Xubiao Peng
- Center for Quantum Technology Research, Key Laboratory of Advanced Optoelectronic Quantum Architecture and Measurements (MOE), School of Physics, Beijing Institute of Technology, Beijing 100081, China; Beijing Academy of Quantum Information Sciences, Beijing 100193, China.
| |
Collapse
|
64
|
Fatafta H, Khaled M, Owen MC, Sayyed-Ahmad A, Strodel B. Amyloid-β peptide dimers undergo a random coil to β-sheet transition in the aqueous phase but not at the neuronal membrane. Proc Natl Acad Sci U S A 2021; 118:e2106210118. [PMID: 34544868 PMCID: PMC8488611 DOI: 10.1073/pnas.2106210118] [Citation(s) in RCA: 56] [Impact Index Per Article: 18.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/16/2021] [Indexed: 11/21/2022] Open
Abstract
Mounting evidence suggests that the neuronal cell membrane is the main site of oligomer-mediated neuronal toxicity of amyloid-β peptides in Alzheimer's disease. To gain a detailed understanding of the mutual interference of amyloid-β oligomers and the neuronal membrane, we carried out microseconds of all-atom molecular dynamics (MD) simulations on the dimerization of amyloid-β (Aβ)42 in the aqueous phase and in the presence of a lipid bilayer mimicking the in vivo composition of neuronal membranes. The dimerization in solution is characterized by a random coil to β-sheet transition that seems on pathway to amyloid aggregation, while the interactions with the neuronal membrane decrease the order of the Aβ42 dimer by attenuating its propensity to form a β-sheet structure. The main lipid interaction partners of Aβ42 are the surface-exposed sugar groups of the gangliosides GM1. As the neurotoxic activity of amyloid oligomers increases with oligomer order, these results suggest that GM1 is neuroprotective against Aβ-mediated toxicity.
Collapse
Affiliation(s)
- Hebah Fatafta
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, 52425 Jülich, Germany
| | - Mohammed Khaled
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, 52425 Jülich, Germany
| | - Michael C Owen
- Central European Institute of Technology, Masaryk University, Brno 625 00, Czech Republic
- Institute of Chemistry, University of Miskolc, 3515 Miskolc-Egyetemváros, Hungary
| | | | - Birgit Strodel
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, 52425 Jülich, Germany;
- Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| |
Collapse
|
65
|
Molecular dynamics study of conformation transition from helix to sheet of Aβ42 peptide. J Mol Graph Model 2021; 109:108027. [PMID: 34534891 DOI: 10.1016/j.jmgm.2021.108027] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2021] [Revised: 09/04/2021] [Accepted: 09/05/2021] [Indexed: 11/23/2022]
Abstract
Aβ42 peptides can form helix and sheet structure under different conditions. The conformational conversion is closely associated with Aβ peptides aggregation and their neurotoxicity. But the transition from helix to sheet is not be clearly understood. In this study we performed microsecond timescale MD simulations of Aβ42 peptide to investigate the conformation transition from α-helix to β-sheet. Markov state model (MSM) was built to facilitate identification of crucial intermediate states and possible transition pathway. Based on the analysis, we found that the region Y10-A21 in the middle of Aβ42 peptide plays an initial role in this transition. MSM model revealed that the collapse of helical structure in this region might trigger the formation of sheet structure. Moreover, we further simulated the aggregation of Aβ42 peptides with different conformations. We found that the Aβ42 peptides forming sheet structure have higher aggregation potential compared with peptides with helix structure. These results demonstrate that we can prevent the aggregation of Aβ42 peptides by stabilizing the helix structure in the region of Y10-A21. In addition, this study provides new insight into better understanding the conformational transition and aggregation of Aβ42 peptides.
Collapse
|
66
|
Ferreira ESC, Voroshylova IV, Figueiredo NM, Cordeiro MNDS. Molecular dynamic study of alcohol-based deep eutectic solvents. J Chem Phys 2021; 155:064506. [PMID: 34391364 DOI: 10.1063/5.0058561] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
The applicability of deep eutectic solvents is determined by their physicochemical properties. In turn, the properties of eutectic mixtures are the result of the components' molar ratio and chemical composition. Owing to the relatively low viscosities displayed by alcohol-based deep eutectic solvents (DESs), their application in industry is more appealing. Modeling the composition-property relationships established in polyalcohol-based mixtures is crucial for both understanding and predicting their behavior. In this work, a physicochemical property-structure comparison study is made between four choline chloride polyalcohol-based DESs, namely, ethaline, propeline, propaneline, and glyceline. Physicochemical properties obtained from molecular dynamic simulations are compared to experimental data, whenever possible. The simulations cover the temperature range from 298.15 to 348.15 K. The simulated and literature experimental data are generally in good agreement for all the studied DESs. Structural properties, such as radial and spatial distribution functions, coordination numbers, hydrogen bond donor (HBD)-HBD aggregate formation, and hydrogen bonding are analyzed in detail. The higher prevalence of HBD:HBD and HBD:anion hydrogen bonds is likely to be the major reason for the relatively high density and viscosity of glyceline as well as for lower DES self-diffusions.
Collapse
Affiliation(s)
- Elisabete S C Ferreira
- LAQV@REQUIMTE, Faculdade de Ciências, Universidade do Porto, Departamento de Química e Bioquímica, Rua do Campo Alegre, 4169-007 Porto, Portugal
| | - Iuliia V Voroshylova
- LAQV@REQUIMTE, Faculdade de Ciências, Universidade do Porto, Departamento de Química e Bioquímica, Rua do Campo Alegre, 4169-007 Porto, Portugal
| | - Nádia M Figueiredo
- LAQV@REQUIMTE, Faculdade de Ciências, Universidade do Porto, Departamento de Química e Bioquímica, Rua do Campo Alegre, 4169-007 Porto, Portugal
| | - M Natália D S Cordeiro
- LAQV@REQUIMTE, Faculdade de Ciências, Universidade do Porto, Departamento de Química e Bioquímica, Rua do Campo Alegre, 4169-007 Porto, Portugal
| |
Collapse
|
67
|
Samantray S, Strodel B. The Effects of Different Glycosaminoglycans on the Structure and Aggregation of the Amyloid-β (16-22) Peptide. J Phys Chem B 2021; 125:5511-5525. [PMID: 34027669 DOI: 10.1021/acs.jpcb.1c00868] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Aggregates of the amyloid-β (Aβ) peptide are implicated as a causative substance in Alzheimer's disease. Molecular dynamics simulations provide valuable contributions for elucidating the conformational transitions of monomeric and aggregated forms of Aβ be it in solution or in the presence of other molecules. Here, we study the effects of four different glycosaminoglycans (GAGs), three sulfated ones and a nonsulfated one, on the aggregation of Aβ16-22. From experiments, it has been suggested that GAGs, which belong to the main components of the brain's extracellular space, favor amyloid fibril formation. Our simulation results reveal that the binding of Aβ16-22 to the GAGs is driven by electrostatic attraction between the negative GAG charges and the positively charged K16 of Aβ16-22. While these interactions have only minor effects on the GAG and Aβ16-22 conformations at the 1 Aβ16-22/1 GAG ratio, at the 2:2 stoichiometry the aggregation of Aβ16-22 is considerably changed. In solution, the aggregation of Aβ16-22 is strongly influenced by K16-E22 attraction, leading to antiparallel β-sheets. In the presence of GAGs, on the other hand, the interaction of K16 with the GAGs increases the importance of the hydrophobic interactions during Aβ16-22 aggregation, which in turn yields parallel alignments. A templating and ordering effect of the GAGs on the Aβ16-22 aggregates is observed. In summary, this study provides new insight at the atomic level on GAG-amyloid interactions, strengthening the view that sulfation of the GAGs plays a major role in this context.
Collapse
Affiliation(s)
- Suman Samantray
- Institute of Biological Information Processing: Structural Biochemistry (IBI-7), Forschungszentrum Jülich, 52428 Jülich, Germany.,AICES Graduate School, RWTH Aachen University, Schinkelstraße 2, 52062 Aachen, Germany
| | - Birgit Strodel
- Institute of Biological Information Processing: Structural Biochemistry (IBI-7), Forschungszentrum Jülich, 52428 Jülich, Germany.,Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| |
Collapse
|
68
|
Assessment of transferable forcefields for protein simulations attests improved description of disordered states and secondary structure propensities, and hints at multi-protein systems as the next challenge for optimization. Comput Struct Biotechnol J 2021; 19:2626-2636. [PMID: 34025949 PMCID: PMC8120800 DOI: 10.1016/j.csbj.2021.04.050] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2020] [Revised: 04/20/2021] [Accepted: 04/22/2021] [Indexed: 01/02/2023] Open
Abstract
Continuous assessment of transferable forcefields for molecular simulations is essential to identify their weaknesses and direct improvement efforts. The latest efforts focused on better describing disordered proteins while retaining proper description of folded domains, important because forcefields of the previous generations produce overly compact disordered states. Such improvements should additionally alleviate the related problem of over-stabilized protein–protein interactions, which has been largely overlooked. Here we evaluated three state-of-the-art forcefields, current flagships of their respective developers, optimized for ordered and disordered proteins: CHARMM36m with its recommended corrected TIP3P* water, ff19SB with the recommended OPC water, and the 2019 a99SBdisp forcefield by D. E. Shaw Research with its modified TIP4P water; plus ff14SB with TIP3P as an example of the former generation of forcefields. Our evaluation entailed simulations of (i) multiple copies of a protein that is highly soluble yet undergoes weak dimerization, (ii) a disordered peptide with low, well-characterized alpha helical propensity, and (iii) a peptide known to form insoluble β-aggregates. Our results recapitulate ff14SB-TIP3P over-stabilizing aggregates and secondary structures and place a99SBdisp-TIP4PD at the other end i.e. predicting overly weak intermolecular interactions despite reasonably predicting secondary structure propensities. In-between, CHARMM36m-TIP3P* still over-stabilizes aggregates but predicts residue-wise alpha helical propensities in solution slightly better than ff19SB-OPC, while ff19SB-OPC poses the best prediction of weak dimerization of the soluble protein still predicting aggregation of the β-peptides. This independent assessment shows that the claimed forcefield improvements are real, but also that a right balance between noncovalent attraction and repulsion has not yet been reached. We thus propose developers to consider systems like those tested here in their forcefield tuning protocols. Last, the good performance of CHARMM36m-TIP3P* further shows that tuning 3-point water models might still be an alternative to the more costly 4-point models like OPC and TIP4PD.
Collapse
|
69
|
Roy R, Paul S. Potential of ATP toward Prevention of hIAPP Oligomerization and Destabilization of hIAPP Protofibrils: An In Silico Perspective. J Phys Chem B 2021; 125:3510-3526. [PMID: 33792323 DOI: 10.1021/acs.jpcb.1c00313] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
The aggregation of an intrinsically disordered protein, human islet amyloid polypeptide (hIAPP), leads to one of the most prevalent endocrine disorders, type II diabetes mellitus (T2DM). Hence inhibition of hIAPP aggregation provides a possible therapeutic approach for the treatment of T2DM. In this regard, a new aspect of adenosine triphosphate (ATP), which is widely known as the energy source for biological reactions, has recently been discovered, where it can inhibit the formation of protein aggregates and simultaneously dissolve preformed aggregates at a millimolar concentration scale. In this work, we investigate the effect of ATP on the aggregation of an amyloidogenic segment of hIAPP, hIAPP22-28, and also of the full length sequence. Using all-atom classical molecular dynamics simulations, we observe that the tendency of hIAPP to oligomerize into β-sheet conformers is inhibited by ATP, due to which the peptides remain distant, loosely packed random monomers. Moreover, it can also disassemble preformed hIAPP protofibrils. ATP preferentially interacts with the hydrophobic residues of hIAPP22-28 fragment and the terminal and turn residues of the full length peptide. The hydrogen bonding, hydrophobic, π-π, and N-H-π stacking interactions are the driving forces for the ATP induced inhibition of hIAPP aggregation. Interestingly, the hydrophobic adenosine of ATP is found to be more in contact with the peptide residues than the hydrophilic triphosphate moiety. The insight into the inhibitory mechanism of ATP on hIAPP aggregation can prove to be beneficial for the design of novel amyloid inhibitors in the future.
Collapse
Affiliation(s)
- Rituparna Roy
- Department of Chemistry, Indian Institute of Technology, Guwahati, Assam, India 781039
| | - Sandip Paul
- Department of Chemistry, Indian Institute of Technology, Guwahati, Assam, India 781039
| |
Collapse
|
70
|
Ngo ST, Nguyen PH, Derreumaux P. Cholesterol Molecules Alter the Energy Landscape of Small Aβ1-42 Oligomers. J Phys Chem B 2021; 125:2299-2307. [PMID: 33646777 DOI: 10.1021/acs.jpcb.1c00036] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Small amyloid-β (Aβ) oligomers are believed to be key pathogenic species in Alzheimer's disease (AD). One suggested toxicity mechanism is the detergent model where oligomers remove lipid molecules from the bilayer. Senile plaques of AD patients also accumulate a 1:1 ratio of cholesterol/Aβ. What are the dominant structures of small Aβ42 oligomers with cholesterol molecules in aqueous solution? Here, we answer this question by performing atomistic replica exchange molecular dynamics simulations of Aβ42 dimers and trimers. Our simulations demonstrate that the interactions with cholesterol molecules change completely the energy landscape of small Aβ42 oligomers. This result shows that simulations in the bulk solution cannot recapitulate aggregation in the brain extracellular space.
Collapse
Affiliation(s)
- Son Tung Ngo
- Laboratory of Theoretical and Computational Biophysics & Faculty of Applied Sciences, Ton Duc Thang University, 33000 Ho Chi Minh City, Vietnam
| | - Phuong H Nguyen
- Laboratoire de Biochimie Théorique, CNRS, Université de Paris, UPR 9080, 13 rue Pierre et Marie Curie, 75005 Paris, France.,Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, 75000 Paris, France
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, CNRS, Université de Paris, UPR 9080, 13 rue Pierre et Marie Curie, 75005 Paris, France.,Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, 75000 Paris, France.,Laboratory of Theoretical Chemistry, Ton Duc Thang University, 33000 Ho Chi Minh City, Vietnam.,Faculty of Pharmacy, Ton Duc Thang University, 33000 Ho Chi Minh City, Vietnam
| |
Collapse
|
71
|
Gomes GN, Levine ZA. Defining the Neuropathological Aggresome across in Silico, in Vitro, and ex Vivo Experiments. J Phys Chem B 2021; 125:1974-1996. [PMID: 33464098 PMCID: PMC8362740 DOI: 10.1021/acs.jpcb.0c09193] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The loss of proteostasis over the life course is associated with a wide range of debilitating degenerative diseases and is a central hallmark of human aging. When left unchecked, proteins that are intrinsically disordered can pathologically aggregate into highly ordered fibrils, plaques, and tangles (termed amyloids), which are associated with countless disorders such as Alzheimer's disease, Parkinson's disease, type II diabetes, cancer, and even certain viral infections. However, despite significant advances in protein folding and solution biophysics techniques, determining the molecular cause of these conditions in humans has remained elusive. This has been due, in part, to recent discoveries showing that soluble protein oligomers, not insoluble fibrils or plaques, drive the majority of pathological processes. This has subsequently led researchers to focus instead on heterogeneous and often promiscuous protein oligomers. Unfortunately, significant gaps remain in how to prepare, model, experimentally corroborate, and extract amyloid oligomers relevant to human disease in a systematic manner. This Review will report on each of these techniques and their successes and shortcomings in an attempt to standardize comparisons between protein oligomers across disciplines, especially in the context of neurodegeneration. By standardizing multiple techniques and identifying their common overlap, a clearer picture of the soluble neuropathological aggresome can be constructed and used as a baseline for studying human disease and aging.
Collapse
Affiliation(s)
- Gregory-Neal Gomes
- Department of Pathology, Yale School of Medicine, New Haven, CT, 06520, USA
- Department of Molecular Biophysics & Biochemistry, Yale University, New Haven, CT 06511, USA
| | - Zachary A. Levine
- Department of Pathology, Yale School of Medicine, New Haven, CT, 06520, USA
- Department of Molecular Biophysics & Biochemistry, Yale University, New Haven, CT 06511, USA
| |
Collapse
|
72
|
Mu J, Liu H, Zhang J, Luo R, Chen HF. Recent Force Field Strategies for Intrinsically Disordered Proteins. J Chem Inf Model 2021; 61:1037-1047. [PMID: 33591749 PMCID: PMC8256680 DOI: 10.1021/acs.jcim.0c01175] [Citation(s) in RCA: 42] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Intrinsically disordered proteins (IDPs) are widely distributed across eukaryotic cells, playing important roles in molecular recognition, molecular assembly, post-translational modification, and other biological processes. IDPs are also associated with many diseases such as cancers, cardiovascular diseases, and neurodegenerative diseases. Due to their structural flexibility, conventional experimental methods cannot reliably capture their heterogeneous structures. Molecular dynamics simulation becomes an important complementary tool to quantify IDP structures. This review covers recent force field strategies proposed for more accurate molecular dynamics simulations of IDPs. The strategies include adjusting dihedral parameters, adding grid-based energy correction map (CMAP) parameters, refining protein-water interactions, and others. Different force fields were found to perform well on specific observables of specific IDPs but also are limited in reproducing all available experimental observables consistently for all tested IDPs. We conclude the review with perspective areas for improvements for future force fields for IDPs.
Collapse
Affiliation(s)
- Junxi Mu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Hao Liu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Jian Zhang
- Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, School of Medicine, Shanghai Jiao Tong University, Shanghai 20025, China
| | - Ray Luo
- Departments of Molecular Biology and Biochemistry, Chemical and Molecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine, California 92697-3900, United States
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| |
Collapse
|
73
|
Fatafta H, Samantray S, Sayyed-Ahmad A, Coskuner-Weber O, Strodel B. Molecular simulations of IDPs: From ensemble generation to IDP interactions leading to disorder-to-order transitions. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2021; 183:135-185. [PMID: 34656328 DOI: 10.1016/bs.pmbts.2021.06.003] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/26/2023]
Abstract
Intrinsically disordered proteins (IDPs) lack a well-defined three-dimensional structure but do exhibit some dynamical and structural ordering. The structural plasticity of IDPs indicates that entropy-driven motions are crucial for their function. Many IDPs undergo function-related disorder-to-order transitions upon by their interaction with specific binding partners. Approaches that are based on both experimental and theoretical tools enable the biophysical characterization of IDPs. Molecular simulations provide insights into IDP structural ensembles and disorder-to-order transition mechanisms. However, such studies depend strongly on the chosen force field parameters and simulation techniques. In this chapter, we provide an overview of IDP characteristics, review all-atom force fields recently developed for IDPs, and present molecular dynamics-based simulation methods that allow IDP ensemble generation as well as the characterization of disorder-to-order transitions. In particular, we introduce metadynamics, replica exchange molecular dynamics simulations, and also kinetic models resulting from Markov State modeling, and provide various examples for the successful application of these simulation methods to IDPs.
Collapse
Affiliation(s)
- Hebah Fatafta
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany
| | - Suman Samantray
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany; AICES Graduate School, RWTH Aachen University, Aachen, Germany
| | | | - Orkid Coskuner-Weber
- Molecular Biotechnology, Turkish-German University, Sahinkaya Caddesi, Istanbul, Turkey
| | - Birgit Strodel
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany; Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany.
| |
Collapse
|
74
|
Strodel B. Amyloid aggregation simulations: challenges, advances and perspectives. Curr Opin Struct Biol 2020; 67:145-152. [PMID: 33279865 DOI: 10.1016/j.sbi.2020.10.019] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2020] [Accepted: 10/18/2020] [Indexed: 10/22/2022]
Abstract
In amyloid aggregation diseases soluble proteins coalesce into a wide array of undesirable structures, ranging through oligomers and prefibrillar assemblies to highly ordered amyloid fibrils and plaques. Explicit-solvent all-atom molecular dynamics (MD) simulations of amyloid aggregation have been performed for almost 20 years, revealing valuable information about this phenomenon. However, these simulations are challenged by three main problems. Firstly, current force fields modeling amyloid aggregation are insufficiently accurate. Secondly, the protein concentrations in MD simulations are usually orders of magnitude higher than those used in vitro or found in vivo, which has direct consequences on the aggregates that form. Finally, the third problem is the well-known time-scale limit of MD simulations. In this review I highlight recent approaches to overcome these three limitations.
Collapse
Affiliation(s)
- Birgit Strodel
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, 52425 Jülich, Germany; Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Universitätstrasse 1, 40225 Düsseldorf, Germany.
| |
Collapse
|