1
|
Akter S, Islam MJ, Ali MA, Zakaria Tashrif M, Uddin MJ, Ullah MO, Halim MA. Structure and dynamics of whole-sequence homology model of ORF3a protein of SARS-CoV-2: An insight from microsecond molecular dynamics simulations. J Biomol Struct Dyn 2024; 42:6726-6739. [PMID: 37528650 DOI: 10.1080/07391102.2023.2236715] [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] [Received: 09/22/2022] [Accepted: 07/08/2023] [Indexed: 08/03/2023]
Abstract
The ORF3a is a large accessory protein in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which plays an important role in virulence and viral replication; especially in inflammasome activation and apoptosis. However,, the existing cryo-EM structure of SARS-CoV-2 ORF3a is incomplete, . making it challenging to understand its structural and functional features. The aim of this study is to investigate the dynamic behaviors of the full-sequence homology model of ORF3a and compare it with the cryo-EM structure using microsecond molecular dynamics simulations. The previous studies indicated that the unresolved residues of the cryo-EM structure are not only involved in the pathogenesis of the SARS-CoV-2 but also exhibit a significant antigenicity. The dynamics scenario of homology model revealed higher RMSD, Rg, and SASA values with stable pattern when compared to the cryo-EM structure. Moreover, the RMSF analysis demonstrated higher fluctuations at specific positions (1-43, 97-110, 172-180, 219-243) in the model structure, whereas the cryo-EM structure displayed lower overall drift (except 1-43) in comparison to the model structure.Secondary structural features indicated that a significant unfolding in the transmembrane domains and β-strand at positions 166 to 172, affecting the stability and compactness of the cryo-EM structure , whereas the model exhibited noticeable unfolding in transmembrane domains and small-coiled regions in the N-terminal. , The results from molecular docking and steered molecular dynamics investigations showed the model structure had a greater number of non-bonding interactions, leading to enhanced stability when compared to the cryo-EM structure. Consequently, higher forces were necessary for unbinding of the baricitinib and ruxolitinib inhibitors from the model structure.. Our findings can help better understanding of the significance of unresolved residues at the molecular level. Additionally, this information can guide researchers for experimental endeavors aimed at completing the full-sequence structure of the ORF3a.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Shaila Akter
- Division of Infectious Diseases and Division of Computer-Aided Drug Design, The Red-Green Research Centre, BICCB, Tejgaon, Dhaka, Bangladesh
| | - Md Jahirul Islam
- Division of Infectious Diseases and Division of Computer-Aided Drug Design, The Red-Green Research Centre, BICCB, Tejgaon, Dhaka, Bangladesh
| | - Md Ackas Ali
- Department of Chemistry and Biochemistry, Kennesaw State University, Kennesaw, GA, USA
| | - Md Zakaria Tashrif
- Division of Infectious Diseases and Division of Computer-Aided Drug Design, The Red-Green Research Centre, BICCB, Tejgaon, Dhaka, Bangladesh
| | - Md Jaish Uddin
- Division of Infectious Diseases and Division of Computer-Aided Drug Design, The Red-Green Research Centre, BICCB, Tejgaon, Dhaka, Bangladesh
| | - M Obayed Ullah
- Division of Infectious Diseases and Division of Computer-Aided Drug Design, The Red-Green Research Centre, BICCB, Tejgaon, Dhaka, Bangladesh
| | - Mohammad A Halim
- Department of Chemistry and Biochemistry, Kennesaw State University, Kennesaw, GA, USA
| |
Collapse
|
2
|
Sunda AP, Sharma AK. Molecular Insights into Cu/Zn Metal Response to the Amyloid β-Peptide (1-42). ACS PHYSICAL CHEMISTRY AU 2024; 4:57-66. [PMID: 38283784 PMCID: PMC10811771 DOI: 10.1021/acsphyschemau.3c00041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Revised: 10/03/2023] [Accepted: 10/04/2023] [Indexed: 01/30/2024]
Abstract
Aβ1-40 peptide and Aβ1-42 peptide are the building units of beta-amyloid plaques present in Alzheimer's disease (AD)-affected brain. The binding affinity of various divalent metal ions such as Cu and Zn present in AD-affected brain with different amino acids available in Aβ-peptide became the focus to explore their role in soluble neurotoxic oligomer formation. Cu2+ metal ions are known to enhance the neurotoxicity of the Aβ1-42 peptide by catalyzing the formation of soluble neurotoxic oligomers. The competitive preference of both Cu2+ and Zn2+ simultaneously to interact with the Aβ-peptide is unknown. The divalent Cu and Zn ions were inserted in explicit aqueous Aβ1-42 peptide configurations to get insights into the binding competence of these metal ions with peptides using classical molecular dynamics (MD) simulations. The metal-ion interactions reveal that competitive binding preferences of various peptide sites become metal-ion-specific and differ significantly. For Cu2+, interactions are found to be more significant with respect to those of Asp-7, His-6, Glu-11, and His-14. Asp-1, Glu-3, Asp-7, His-6, Glu-11, and His-13 amino acid residues show higher affinity toward Zn2+ ions. MD simulations show notable variation in the solvent-accessible surface area in the hydrophobic region of the peptide. Infinitesimal mobility was obtained for Zn2+ compared to Cu2+ in an aqueous solution and Cu2+ diffusivity deviated significantly at different time scales, proving its labile features in aqueous Aβ1-42 peptides.
Collapse
Affiliation(s)
- Anurag Prakash Sunda
- Department
of Chemistry, J. C. Bose University of Science
and Technology, YMCA, Faridabad 121006, India
| | - Anuj Kumar Sharma
- Department
of Chemistry, School of Chemical Sciences and Pharmacy, Central University of Rajasthan, Ajmer 305817, India
| |
Collapse
|
3
|
Bhattacharya P, Mandal A. Identification of amentoflavone as a potent SARS-CoV-2 M pro inhibitor: a combination of computational studies and in vitro biological evaluation. J Biomol Struct Dyn 2024:1-19. [PMID: 38263736 DOI: 10.1080/07391102.2024.2304676] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2023] [Accepted: 01/07/2024] [Indexed: 01/25/2024]
Abstract
Small-molecule inhibitors of SARS-CoV-2 Mpro that block the active site pocket of the viral main protease have been considered potential therapeutics for the development of drugs against SARS-CoV-2. Here, we report the identification of amentoflavone (a biflavonoid) through docking-based virtual screening of a library comprised of 231 compounds consisting of flavonoids and isoflavonoids. The docking results were further substantiated through extensive analysis of the data obtained from all-atom 150 ns MD simulation. End-state effective free energy calculations using MM-PBSA calculations further suggested that (Ra)-amentoflavone (C3'-C8''-atropisomer) may show a greater binding affinity towards the Mpro than (Sa)-amentoflavone. In vitro cytotoxicity assay established that amentoflavone showed a high CC50 value indicating much lower toxicity. Further, potent inhibition of the Mpro by amentoflavone was established by studying the effect on HEK293T cells treated with SARS-CoV-2 Mpro expressing plasmid.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
| | - Anirban Mandal
- Department of Microbiology, Mrinalini Datta Mahavidyapith, Kolkata, India
| |
Collapse
|
4
|
Fukuhara D, Yamauchi M, Itoh SG, Okumura H. Ingenuity in performing replica permutation: How to order the state labels for improving sampling efficiency. J Comput Chem 2023; 44:534-545. [PMID: 36346137 PMCID: PMC10099539 DOI: 10.1002/jcc.27020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2022] [Revised: 09/19/2022] [Accepted: 09/21/2022] [Indexed: 11/11/2022]
Abstract
In the replica-permutation method, an advanced version of the replica-exchange method, all combinations of replicas and parameters are considered for parameter permutation, and a list of all the combinations is prepared. Here, we report that the temperature transition probability depends on how the list is created, especially in replica permutation with solute tempering (RPST). We found that the transition probabilities decrease at large replica indices when the combinations are sequentially assigned to the state labels as in the originally proposed list. To solve this problem, we propose to modify the list by randomly assigning the combinations to the state labels. We performed molecular dynamics simulations of amyloid-β(16-22) peptides using RPST with the "randomly assigned" list (RPST-RA) and RPST with the "sequentially assigned" list (RPST-SA). The results show the decreases in the transition probabilities in RPST-SA are eliminated, and the sampling efficiency is improved in RPST-RA.
Collapse
Affiliation(s)
- Daiki Fukuhara
- Department of Structural Molecular Science, SOKENDAI (The Graduate University for Advanced Studies), Okazaki, Japan.,Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Japan
| | - Masataka Yamauchi
- Department of Structural Molecular Science, SOKENDAI (The Graduate University for Advanced Studies), Okazaki, Japan.,Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Japan
| | - Satoru G Itoh
- Department of Structural Molecular Science, SOKENDAI (The Graduate University for Advanced Studies), Okazaki, Japan.,Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Japan.,Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki, Japan
| | - Hisashi Okumura
- Department of Structural Molecular Science, SOKENDAI (The Graduate University for Advanced Studies), Okazaki, Japan.,Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Japan.,Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki, Japan
| |
Collapse
|
5
|
Lui LH, Egbu R, Graver T, Williams GR, Brocchini S, Velayudhan A. Computational and Experimental Evaluation of the Stability of a GLP-1-like Peptide in Ethanol–Water Mixtures. Pharmaceutics 2022; 14:pharmaceutics14071462. [PMID: 35890357 PMCID: PMC9321252 DOI: 10.3390/pharmaceutics14071462] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2022] [Revised: 07/05/2022] [Accepted: 07/12/2022] [Indexed: 12/04/2022] Open
Abstract
Aggregation resulting from the self-association of peptide molecules remains a major challenge during preformulation. Whereas certain organic solvents are known to promote aggregation, ethanol (EtOH) is capable of disrupting interactions between peptide molecules. It is unclear whether it is beneficial or counterproductive to include EtOH in formulations of short peptides. Here, we employed molecular dynamics simulations using the DAFT protocol and MARTINI force field to predict the formation of self-associated dimers and to estimate the stability of a GLP-1-like peptide (G48) in 0–80% aqueous EtOH solutions. Both simulation and experimental data reveal that EtOH leads to a remarkable increase in the conformational stability of the peptide when stored over 15 days at 27 °C. In the absence of EtOH, dimerisation and subsequent loss in conformational stability (α-helix → random coil) were observed. EtOH improved conformational stability by reducing peptide–peptide interactions. The data suggest that a more nuanced approach may be applied in formulation decision making and, if the native state of the peptide is an α-helix organic solvent, such as EtOH, may enhance stability and improve prospects of long-term storage.
Collapse
Affiliation(s)
- Lok Hin Lui
- UCL School of Pharmacy, University College London, London WC1N 1AX, UK; (L.H.L.); (R.E.); (T.G.); (G.R.W.); (S.B.)
| | - Raphael Egbu
- UCL School of Pharmacy, University College London, London WC1N 1AX, UK; (L.H.L.); (R.E.); (T.G.); (G.R.W.); (S.B.)
| | - Thomas Graver
- UCL School of Pharmacy, University College London, London WC1N 1AX, UK; (L.H.L.); (R.E.); (T.G.); (G.R.W.); (S.B.)
| | - Gareth R. Williams
- UCL School of Pharmacy, University College London, London WC1N 1AX, UK; (L.H.L.); (R.E.); (T.G.); (G.R.W.); (S.B.)
| | - Steve Brocchini
- UCL School of Pharmacy, University College London, London WC1N 1AX, UK; (L.H.L.); (R.E.); (T.G.); (G.R.W.); (S.B.)
| | - Ajoy Velayudhan
- Department of Biochemical Engineering, University College London, London WC1E 6BT, UK
- Correspondence:
| |
Collapse
|
6
|
Molecular Dynamics Simulation Studies on the Aggregation of Amyloid-β Peptides and Their Disaggregation by Ultrasonic Wave and Infrared Laser Irradiation. Molecules 2022; 27:molecules27082483. [PMID: 35458686 PMCID: PMC9030874 DOI: 10.3390/molecules27082483] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2022] [Revised: 03/29/2022] [Accepted: 04/07/2022] [Indexed: 01/02/2023] Open
Abstract
Alzheimer’s disease is understood to be caused by amyloid fibrils and oligomers formed by aggregated amyloid-β (Aβ) peptides. This review article presents molecular dynamics (MD) simulation studies of Aβ peptides and Aβ fragments on their aggregation, aggregation inhibition, amyloid fibril conformations in equilibrium, and disruption of the amyloid fibril by ultrasonic wave and infrared laser irradiation. In the aggregation of Aβ, a β-hairpin structure promotes the formation of intermolecular β-sheet structures. Aβ peptides tend to exist at hydrophilic/hydrophobic interfaces and form more β-hairpin structures than in bulk water. These facts are the reasons why the aggregation is accelerated at the interface. We also explain how polyphenols, which are attracting attention as aggregation inhibitors of Aβ peptides, interact with Aβ. An MD simulation study of the Aβ amyloid fibrils in equilibrium is also presented: the Aβ amyloid fibril has a different structure at one end from that at the other end. The amyloid fibrils can be destroyed by ultrasonic wave and infrared laser irradiation. The molecular mechanisms of these amyloid fibril disruptions are also explained, particularly focusing on the function of water molecules. Finally, we discuss the prospects for developing treatments for Alzheimer’s disease using MD simulations.
Collapse
|
7
|
Fukuhara D, Itoh SG, Okumura H. Replica permutation with solute tempering for molecular dynamics simulation and its application to the dimerization of amyloid-β fragments. J Chem Phys 2022; 156:084109. [DOI: 10.1063/5.0081686] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
We propose the replica permutation with solute tempering (RPST) by combining the replica-permutation method (RPM) and the replica exchange with solute tempering (REST). Temperature permutations are performed among more than two replicas in RPM, whereas temperature exchanges are performed between two replicas in the replica-exchange method (REM). The temperature transition in RPM occurs more efficiently than in REM. In REST, only the temperatures of the solute region, the solute temperatures, are exchanged to reduce the number of replicas compared to REM. Therefore, RPST is expected to be an improved method taking advantage of these methods. For comparison, we applied RPST, REST, RPM, and REM to two amyloid-β(16–22) peptides in explicit water. We calculated the transition ratio and the number of tunneling events in the temperature space and the number of dimerization events of amyloid-β(16–22) peptides. The results indicate that, in RPST, the number of replicas necessary for frequent random walks in the temperature and conformational spaces is reduced compared to the other three methods. In addition, we focused on the dimerization process of amyloid-β(16–22) peptides. The RPST simulation with a relatively small number of replicas shows that the two amyloid-β(16–22) peptides form the intermolecular antiparallel β-bridges due to the hydrophilic side-chain contact between Lys and Glu and hydrophobic side-chain contact between Leu, Val, and Phe, which stabilizes the dimer of the peptides.
Collapse
Affiliation(s)
- Daiki Fukuhara
- Department of Structural Molecular Science, SOKENDAI (The Graduate University for Advanced Studies), Okazaki, Aichi 444-8787, Japan
- Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan
| | - Satoru G. Itoh
- Department of Structural Molecular Science, SOKENDAI (The Graduate University for Advanced Studies), Okazaki, Aichi 444-8787, Japan
- Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan
- Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan
| | - Hisashi Okumura
- Department of Structural Molecular Science, SOKENDAI (The Graduate University for Advanced Studies), Okazaki, Aichi 444-8787, Japan
- Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan
- Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan
| |
Collapse
|
8
|
Tachi Y, Itoh SG, Okumura H. Molecular dynamics simulations of amyloid-β peptides in heterogeneous environments. Biophys Physicobiol 2022; 19:1-18. [PMID: 35666692 PMCID: PMC9135617 DOI: 10.2142/biophysico.bppb-v19.0010] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2021] [Accepted: 03/31/2022] [Indexed: 12/01/2022] Open
Affiliation(s)
- Yuhei Tachi
- Department of Physics, Graduate school of Science, Nagoya University
| | - Satoru G. Itoh
- Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences
| | - Hisashi Okumura
- Institute for Molecular Science, National Institutes of Natural Sciences
| |
Collapse
|
9
|
Bhattacharya S, Xu L, Thompson D. Characterization of Amyloidogenic Peptide Aggregability in Helical Subspace. Methods Mol Biol 2022; 2340:401-448. [PMID: 35167084 DOI: 10.1007/978-1-0716-1546-1_18] [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
Prototypical amyloidogenic peptides amyloid-β (Aβ) and α-synuclein (αS) can undergo helix-helix associations via partially folded helical conformers, which may influence pathological progression to Alzheimer's (AD) and Parkinson's disease (PD), respectively. At the other extreme, stable folded helical conformers have been reported to resist self-assembly and amyloid formation. Experimental characterisation of such disparities in aggregation profiles due to subtle differences in peptide stabilities is precluded by the conformational heterogeneity of helical subspace. The diverse physical models used in molecular simulations allow sampling distinct regions of the phase space and are extensive in capturing the ensemble of rich helical subspace. Robust and powerful computational predictive methods utilizing network theory and free energy mapping can model the origin of helical population shifts in amyloidogenic peptides, which highlight their inherent aggregability. In this chapter, we discuss computational models, methods, design rules, and strategies to identify the driving force behind helical self-assembly and the molecular origin of aggregation resistance in helical intermediates of Aβ42 and αS. By extensive multiscale mapping of intrapeptide interactions, we show that the computational models can capture features that are otherwise imperceptible to experiments. Our models predict that targeting terminal residues may allow modulation and control of initial pathogenic aggregability of amyloidogenic peptides.
Collapse
Affiliation(s)
- Shayon Bhattacharya
- Department of Physics, Bernal Institute, University of Limerick, Limerick, Ireland
| | - Liang Xu
- Department of Physics, Bernal Institute, University of Limerick, Limerick, Ireland
| | - Damien Thompson
- Department of Physics, Bernal Institute, University of Limerick, Limerick, Ireland.
| |
Collapse
|
10
|
Williams AE, Hammer NI, Fortenberry RC, Reinemann DN. Tracking the Amide I and αCOO- Terminal ν(C=O) Raman Bands in a Family of l-Glutamic Acid-Containing Peptide Fragments: A Raman and DFT Study. Molecules 2021; 26:4790. [PMID: 34443382 PMCID: PMC8399447 DOI: 10.3390/molecules26164790] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2021] [Revised: 08/02/2021] [Accepted: 08/04/2021] [Indexed: 11/16/2022] Open
Abstract
The E-hook of β-tubulin plays instrumental roles in cytoskeletal regulation and function. The last six C-terminal residues of the βII isotype, a peptide of amino acid sequence EGEDEA, extend from the microtubule surface and have eluded characterization with classic X-ray crystallographic techniques. The band position of the characteristic amide I vibration of small peptide fragments is heavily dependent on the length of the peptide chain, the extent of intramolecular hydrogen bonding, and the overall polarity of the fragment. The dependence of the E residue's amide I ν(C=O) and the αCOO- terminal ν(C=O) bands on the neighboring side chain, the length of the peptide fragment, and the extent of intramolecular hydrogen bonding in the structure are investigated here via the EGEDEA peptide. The hexapeptide is broken down into fragments increasing in size from dipeptides to hexapeptides, including EG, ED, EA, EGE, EDE, DEA, EGED, EDEA, EGEDE, GEDEA, and, finally, EGEDEA, which are investigated with experimental Raman spectroscopy and density functional theory (DFT) computations to model the zwitterionic crystalline solids (in vacuo). The molecular geometries and Boltzmann sum of the simulated Raman spectra for a set of energetic minima corresponding to each peptide fragment are computed with full geometry optimizations and corresponding harmonic vibrational frequency computations at the B3LYP/6-311++G(2df,2pd) level of theory. In absence of the crystal structure, geometry sampling is performed to approximate solid phase behavior. Natural bond order (NBO) analyses are performed on each energetic minimum to quantify the magnitude of the intramolecular hydrogen bonds. The extent of the intramolecular charge transfer is dependent on the overall polarity of the fragment considered, with larger and more polar fragments exhibiting the greatest extent of intramolecular charge transfer. A steady blue shift arises when considering the amide I band position moving linearly from ED to EDE to EDEA to GEDEA and, finally, to EGEDEA. However, little variation is observed in the αCOO- ν(C=O) band position in this family of fragments.
Collapse
Affiliation(s)
- Ashley E. Williams
- Department of Chemistry and Biochemistry, University of Mississippi, University, MS 38677, USA; (A.E.W.); (N.I.H.)
| | - Nathan I. Hammer
- Department of Chemistry and Biochemistry, University of Mississippi, University, MS 38677, USA; (A.E.W.); (N.I.H.)
| | - Ryan C. Fortenberry
- Department of Chemistry and Biochemistry, University of Mississippi, University, MS 38677, USA; (A.E.W.); (N.I.H.)
| | - Dana N. Reinemann
- Department of Biomedical Engineering, University of Mississippi, University, MS 38677, USA
- Department of Chemical Engineering, University of Mississippi, University, MS 38677, USA
| |
Collapse
|
11
|
La Sala G, Gunnarsson A, Edman K, Tyrchan C, Hogner A, Frolov AI. Unraveling the Allosteric Cross-Talk between the Coactivator Peptide and the Ligand-Binding Site in the Glucocorticoid Receptor. J Chem Inf Model 2021; 61:3667-3680. [PMID: 34156843 DOI: 10.1021/acs.jcim.1c00323] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The glucocorticoid receptor (GR) is a nuclear receptor that controls critical biological processes by regulating the transcription of specific genes. There is a known allosteric cross-talk between the ligand and coregulator binding sites within the GR ligand-binding domain that is crucial for the control of the functional response. However, the molecular mechanisms underlying such an allosteric control remain elusive. Here, molecular dynamics (MD) simulations, bioinformatic analysis, and biophysical measurements are integrated to capture the structural and dynamic features of the allosteric cross-talk within the GR. We identified a network of evolutionarily conserved residues that enables the allosteric signal transduction, in agreement with experimental data. MD simulations clarify how such a network is dynamically interconnected and offer a mechanistic explanation of how different peptides affect the intensity of the allosteric signal. This study provides useful insights to elucidate the GR allosteric regulation, ultimately providing a foundation for designing novel drugs.
Collapse
Affiliation(s)
- Giuseppina La Sala
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
| | - Anders Gunnarsson
- Discovery Science, BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
| | - Karl Edman
- Discovery Science, BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
| | - Christian Tyrchan
- Medicinal Chemistry, Research and Early Development, Respiratory & Immunology, BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
| | - Anders Hogner
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
| | - Andrey I Frolov
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
| |
Collapse
|
12
|
Okumura H, Itoh SG. Molecular dynamics simulations of amyloid-β(16-22) peptide aggregation at air-water interfaces. J Chem Phys 2021; 152:095101. [PMID: 33480728 DOI: 10.1063/1.5131848] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
Oligomers of amyloid-β (Aβ) peptides are known to be related to Alzheimer's disease, and their formation is accelerated at hydrophilic-hydrophobic interfaces, such as the cell membrane surface and air-water interface. Here, we report molecular dynamics simulations of aggregation of Aβ(16-22) peptides at air-water interfaces. First, 100 randomly distributed Aβ(16-22) peptides moved to the interface. The high concentration of peptides then accelerated their aggregation and formation of antiparallel β-sheets. Two layers of oligomers were observed near the interface. In the first layer from the interface, the oligomer with less β-bridges exposed the hydrophobic residues to the air. The second layer consisted of oligomers with more β-bridges that protruded into water. They are more soluble in water because the hydrophobic residues are covered by N- and C-terminal hydrophilic residues that are aligned well along the oligomer edge. These results indicate that amyloid protofibril formation mainly occurs in the second layer.
Collapse
Affiliation(s)
- Hisashi Okumura
- Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan
| | - Satoru G Itoh
- Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan
| |
Collapse
|
13
|
Promotion and Inhibition of Amyloid-β Peptide Aggregation: Molecular Dynamics Studies. Int J Mol Sci 2021; 22:ijms22041859. [PMID: 33668406 PMCID: PMC7918115 DOI: 10.3390/ijms22041859] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2021] [Revised: 02/10/2021] [Accepted: 02/11/2021] [Indexed: 01/06/2023] Open
Abstract
Aggregates of amyloid-β (Aβ) peptides are known to be related to Alzheimer’s disease. Their aggregation is enhanced at hydrophilic–hydrophobic interfaces, such as a cell membrane surface and air-water interface, and is inhibited by polyphenols, such as myricetin and rosmarinic acid. We review molecular dynamics (MD) simulation approaches of a full-length Aβ peptide, Aβ40, and Aβ(16–22) fragments in these environments. Since these peptides have both hydrophilic and hydrophobic amino acid residues, they tend to exist at the interfaces. The high concentration of the peptides accelerates the aggregation there. In addition, Aβ40 forms a β-hairpin structure, and this structure accelerates the aggregation. We also describe the inhibition mechanism of the Aβ(16–22) aggregation by polyphenols. The aggregation of Aβ(16–22) fragments is caused mainly by the electrostatic attraction between charged amino acid residues known as Lys16 and Glu22. Since polyphenols form hydrogen bonds between their hydroxy and carboxyl groups and these charged amino acid residues, they inhibit the aggregation.
Collapse
|
14
|
Illig AM, Strodel B. Performance of Markov State Models and Transition Networks on Characterizing Amyloid Aggregation Pathways from MD Data. J Chem Theory Comput 2020; 16:7825-7839. [DOI: 10.1021/acs.jctc.0c00727] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Alexander-Maurice Illig
- Institute of Biological Information Processing: Structural Biochemistry (IBI-7), Forschungszentrum Jülich, 52428 Jülich, 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
|
15
|
Samantray S, Yin F, Kav B, Strodel B. Different Force Fields Give Rise to Different Amyloid Aggregation Pathways in Molecular Dynamics Simulations. J Chem Inf Model 2020; 60:6462-6475. [DOI: 10.1021/acs.jcim.0c01063] [Citation(s) in RCA: 51] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Suman Samantray
- Institute of Biological Information Processing: Structural Biochemistry (IBI-7), Forschungszentrum Jülch, 52428 Jülich, Germany
- AICES Graduate School, RWTH Aachen University, Schinkelstraße 2, 52062 Aachen, Germany
| | - Feng Yin
- Institute of Biological Information Processing: Structural Biochemistry (IBI-7), Forschungszentrum Jülch, 52428 Jülich, Germany
| | - Batuhan Kav
- Institute of Biological Information Processing: Structural Biochemistry (IBI-7), Forschungszentrum Jülch, 52428 Jülich, Germany
| | - Birgit Strodel
- Institute of Biological Information Processing: Structural Biochemistry (IBI-7), Forschungszentrum Jülch, 52428 Jülich, Germany
- Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| |
Collapse
|
16
|
Katyal N, Deep S. A computational approach to get insights into multiple faces of additives in modulation of protein aggregation pathways. Phys Chem Chem Phys 2019; 21:24269-24285. [PMID: 31670327 DOI: 10.1039/c9cp03763b] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
An enormous population worldwide is presently confronted with debilitating neurodegenerative diseases. The etiology of the disease is connected to protein aggregation and the events involved therein. Thus, a complete understanding of an inhibitor at different stages in the process is imperative for the formulation of a drug molecule. This review presents a detailed summary of the current status of different cosolvents. It further develops how the complex aggregation pathway can be simplified into three steps common to all proteins and the way computer simulations can be exploited to gain insights into the ways by which known inhibitors can affect all these stages. Computation of theoretical parameters in this regard and their correlation with experimental techniques is accentuated. In addition to providing an outline of the scope of different additives, this review showcases the way by which the problem of analyzing an effect of an additive can be addressed effectively via MD simulations.
Collapse
Affiliation(s)
- Nidhi Katyal
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi-110016, Delhi, India.
| | | |
Collapse
|
17
|
Nguyen HL, Krupa P, Hai NM, Linh HQ, Li MS. Structure and Physicochemical Properties of the Aβ42 Tetramer: Multiscale Molecular Dynamics Simulations. J Phys Chem B 2019; 123:7253-7269. [DOI: 10.1021/acs.jpcb.9b04208] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Hoang Linh Nguyen
- Institute for Computational Science and Technology, SBI Building, Quang Trung Software
City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City 700000, Vietnam
- Biomedical Engineering Department, Ho Chi Minh City University of Technology-VNU HCM, 268 Ly Thuong Kiet Street, Distr. 10, Ho Chi Minh City 700000, Vietnam
| | - Pawel Krupa
- Institute of Physics Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| | - Nguyen Minh Hai
- Faculty of Physics and Engineering Physics, University of Science-VNU HCM, Ho Chi Minh City 700000, Vietnam
| | - Huynh Quang Linh
- Biomedical Engineering Department, Ho Chi Minh City University of Technology-VNU HCM, 268 Ly Thuong Kiet Street, Distr. 10, Ho Chi Minh City 700000, Vietnam
| | - Mai Suan Li
- Institute of Physics Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland
| |
Collapse
|
18
|
Itoh SG, Yagi-Utsumi M, Kato K, Okumura H. Effects of a Hydrophilic/Hydrophobic Interface on Amyloid-β Peptides Studied by Molecular Dynamics Simulations and NMR Experiments. J Phys Chem B 2019; 123:160-169. [PMID: 30543290 DOI: 10.1021/acs.jpcb.8b11609] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Oligomer formation of amyloid-β peptides (Aβ) is accelerated at a hydrophilic/hydrophobic interface. However, details of the acceleration mechanism have not been elucidated. To understand the effects of the interface on oligomerization at the atomic level, we performed molecular dynamics simulations for an Aβ40 monomer in the presence and absence of the hydrophilic/hydrophobic interface. Nuclear magnetic resonance experiments of Aβ40 peptides with gangliosidic micelles were also carried out. In the simulations and experiments, the hydrophobic residues of Aβ40 bound to the interface stably. Moreover, we found that Aβ40 formed a hairpin structure at the interface more readily than in bulk water. From these results, we discussed the acceleration mechanism of the oligomer formation at the interface.
Collapse
Affiliation(s)
- Satoru G Itoh
- Institute for Molecular Science (IMS) , National Institutes of Natural Sciences , Okazaki , Aichi 444-8585 , Japan.,Exploratory Research Center on Life and Living Systems (ExCELLS) , National Institutes of Natural Sciences , Okazaki , Aichi 444-8787 , Japan.,Department of Structural Molecular Science , SOKENDAI (The Graduate University for Advanced Studies) , Okazaki , Aichi 444-8585 , Japan
| | - Maho Yagi-Utsumi
- Institute for Molecular Science (IMS) , National Institutes of Natural Sciences , Okazaki , Aichi 444-8585 , Japan.,Exploratory Research Center on Life and Living Systems (ExCELLS) , National Institutes of Natural Sciences , Okazaki , Aichi 444-8787 , Japan.,Department of Functional Molecular Science , SOKENDAI (The Graduate University for Advanced Studies) , Okazaki , Aichi 444-8787 , Japan.,Graduate School of Pharmaceutical Sciences , Nagoya City University , Nagoya , Aichi 465-8603 , Japan
| | - Koichi Kato
- Institute for Molecular Science (IMS) , National Institutes of Natural Sciences , Okazaki , Aichi 444-8585 , Japan.,Exploratory Research Center on Life and Living Systems (ExCELLS) , National Institutes of Natural Sciences , Okazaki , Aichi 444-8787 , Japan.,Department of Functional Molecular Science , SOKENDAI (The Graduate University for Advanced Studies) , Okazaki , Aichi 444-8787 , Japan.,Graduate School of Pharmaceutical Sciences , Nagoya City University , Nagoya , Aichi 465-8603 , Japan
| | - Hisashi Okumura
- Institute for Molecular Science (IMS) , National Institutes of Natural Sciences , Okazaki , Aichi 444-8585 , Japan.,Exploratory Research Center on Life and Living Systems (ExCELLS) , National Institutes of Natural Sciences , Okazaki , Aichi 444-8787 , Japan.,Department of Structural Molecular Science , SOKENDAI (The Graduate University for Advanced Studies) , Okazaki , Aichi 444-8585 , Japan
| |
Collapse
|
19
|
Carballo-Pacheco M, Ismail AE, Strodel B. On the Applicability of Force Fields To Study the Aggregation of Amyloidogenic Peptides Using Molecular Dynamics Simulations. J Chem Theory Comput 2018; 14:6063-6075. [PMID: 30336669 DOI: 10.1021/acs.jctc.8b00579] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Molecular dynamics simulations play an essential role in understanding biomolecular processes such as protein aggregation at temporal and spatial resolutions which are not attainable by experimental methods. For a correct modeling of protein aggregation, force fields must accurately represent molecular interactions. Here, we study the effect of five different force fields on the oligomer formation of Alzheimer's Aβ16-22 peptide and two of its mutants: Aβ16-22(F19V,F20V), which does not form fibrils, and Aβ16-22(F19L) which forms fibrils faster than the wild type. We observe that while oligomer formation kinetics depends strongly on the force field, structural properties, such as the most relevant protein-protein contacts, are similar between them. The oligomer formation kinetics obtained with different force fields differ more from each other than the kinetics between aggregating and nonaggregating peptides simulated with a single force field. We discuss the difficulties in comparing atomistic simulations of amyloid oligomer formation with experimental observables.
Collapse
Affiliation(s)
- Martín Carballo-Pacheco
- Institute of Complex Systems: Structural Biochemistry (ICS-6) , Forschungszentrum Jülich GmbH , 52425 Jülich , Germany.,AICES Graduate School , RWTH Aachen University , Schinkelstraße 2 , 52062 Aachen , Germany
| | - Ahmed E Ismail
- AICES Graduate School , RWTH Aachen University , Schinkelstraße 2 , 52062 Aachen , Germany.,Aachener Verfahrenstechnik, Faculty of Mechanical Engineering , RWTH Aachen University , Schinkelstraße 2 , 52062 Aachen , Germany
| | - Birgit Strodel
- Institute of Complex Systems: Structural Biochemistry (ICS-6) , Forschungszentrum Jülich GmbH , 52425 Jülich , Germany.,Institute of Theoretical and Computational Chemistry , Heinrich Heine University Düsseldorf , Universitätstrasse 1 , 40225 Düsseldorf , Germany
| |
Collapse
|
20
|
Bhattacharya S, Xu L, Thompson D. Revisiting the earliest signatures of amyloidogenesis: Roadmaps emerging from computational modeling and experiment. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1359] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Affiliation(s)
- Shayon Bhattacharya
- Department of Physics, Bernal InstituteUniversity of LimerickLimerickIreland
| | - Liang Xu
- Department of Physics, Bernal InstituteUniversity of LimerickLimerickIreland
| | - Damien Thompson
- Department of Physics, Bernal InstituteUniversity of LimerickLimerickIreland
| |
Collapse
|
21
|
Katyal N, Deep S. Inhibition of GNNQQNY prion peptide aggregation by trehalose: a mechanistic view. Phys Chem Chem Phys 2018; 19:19120-19138. [PMID: 28702592 DOI: 10.1039/c7cp02912h] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Deposition of amyloid fibrils is the seminal event in the pathogenesis of numerous neurodegenerative diseases. The formation of this amyloid assembly is the manifestation of a cascade of structural transitions including toxic oligomer formation in the early stages of aggregation. Thus a viable therapeutic strategy involves the use of small molecular ligands to interfere with this assembly. In this perspective, we have explored the kinetics of aggregate formation of the fibril forming GNNQQNY peptide fragment from the yeast prion protein SUP35 using multiple all atom MD simulations with explicit solvent and provided mechanistic insights into the way trehalose, an experimentally known aggregation inhibitor, modulates the aggregation pathway. The results suggest that the assimilation process is impeded by different barriers at smaller and larger oligomeric sizes: the initial one being easily surpassed at higher temperatures and peptide concentrations. The kinetic profile demonstrates that trehalose delays the aggregation process by increasing both these activation barriers, specifically the latter one. It increases the sampling of small-sized aggregates that lack the beta sheet conformation. Analysis reveals that the barrier in the growth of larger stable oligomers causes the formation of multiple stable small oligomers which then fuse together bimolecularly. The PCA of 26 properties was carried out to deconvolute the events within the temporary lag phases, which suggested dynamism in lags involving an increase in interchain contacts and burial of SASA. The predominant growth route is monomer addition, which changes to condensation on account of a large number of depolymerisation events in the presence of trehalose. The favourable interaction of trehalose specifically with the sidechain of the peptide promotes crowding of trehalose molecules in its vicinity - the combination of both these factors imparts the observed behaviour. Furthermore, increasing trehalose concentration leads to faster expulsion of water molecules than interpeptide interactions. These expelled water molecules have larger translational movement, suggesting an entropy factor to favor the assembly process. Different conformations observed under this condition suggest the role of water molecules in guiding the morphology of the aggregates as well. A similar scenario exists on increasing peptide concentration.
Collapse
Affiliation(s)
- Nidhi Katyal
- Department of Chemistry, Indian Institute of Technology, Delhi, Hauzkhas, New Delhi, India.
| | - Shashank Deep
- Department of Chemistry, Indian Institute of Technology, Delhi, Hauzkhas, New Delhi, India.
| |
Collapse
|
22
|
Barz B, Liao Q, Strodel B. Pathways of Amyloid-β Aggregation Depend on Oligomer Shape. J Am Chem Soc 2017; 140:319-327. [PMID: 29235346 DOI: 10.1021/jacs.7b10343] [Citation(s) in RCA: 107] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
One of the main research topics related to Alzheimer's disease is the aggregation of the amyloid-β peptide, which was shown to follow different pathways for the two major alloforms of the peptide, Aβ40 and the more toxic Aβ42. Experimental studies emphasized that oligomers of specific sizes appear in the early aggregation process in different quantities and might be the key toxic agents for each of the two alloforms. We use transition networks derived from all-atom molecular dynamics simulations to show that the oligomers leading to the type of oligomer distributions observed in experiments originate from compact conformations. Extended oligomers, on the other hand, contribute more to the production of larger aggregates thus driving the aggregation process. We further demonstrate that differences in the aggregation pathways of the two Aβ alloforms occur as early as during the dimer stage. The higher solvent-exposure of hydrophobic residues in Aβ42 oligomers contributes to the different aggregation pathways of both alloforms and also to the increased cytotoxicity of Aβ42.
Collapse
Affiliation(s)
- Bogdan Barz
- Institute of Complex Systems: Structural Biochemistry (ICS-6), Forschungszentrum Jülich GmbH , 52425 Jülich, Germany.,Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf , 40225 Düsseldorf, Germany
| | - Qinghua Liao
- Institute of Complex Systems: Structural Biochemistry (ICS-6), Forschungszentrum Jülich GmbH , 52425 Jülich, Germany.,Department of Cell and Molecular Biology, Uppsala University , S-75124 Uppsala, Sweden
| | - Birgit Strodel
- Institute of Complex Systems: Structural Biochemistry (ICS-6), Forschungszentrum Jülich GmbH , 52425 Jülich, Germany.,Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf , 40225 Düsseldorf, Germany
| |
Collapse
|
23
|
Matthes D, Gapsys V, Griesinger C, de Groot BL. Resolving the Atomistic Modes of Anle138b Inhibitory Action on Peptide Oligomer Formation. ACS Chem Neurosci 2017; 8:2791-2808. [PMID: 28906103 DOI: 10.1021/acschemneuro.7b00325] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
The diphenyl-pyrazole compound anle138b is a known inhibitor of oligomeric aggregate formation in vitro and in vivo. Therefore, anle138b is considered a promising drug candidate to beneficially interfere with neurodegenerative processes causing devastating pathologies in humans. The atomistic details of the aggregation inhibition mechanism, however, are to date unknown since the ensemble of small nonfibrillar aggregates is structurally heterogeneous and inaccessible to direct structural characterization. Here, we set out to elucidate anle138b's mode of action using all-atom molecular dynamics simulations on the multi-microsecond time scale. By comparing simulations of dimeric to tetrameric aggregates from fragments of four amyloidogenic proteins (Aβ, hTau40, hIAPP, and Sup35N) in the presence and absence of anle138b, we show that the compound reduces the overall number of intermolecular hydrogen bonds, disfavors the sampling of the aggregated state, and remodels the conformational distributions within the small oligomeric peptide aggregates. Most notably, anle138b preferentially interacts with the disordered structure ensemble via its pyrazole moiety, thereby effectively blocking interpeptide main chain interactions and impeding the spontaneous formation of ordered β-sheet structures, in particular those with out-of-register antiparallel β-strands. The structurally very similar compound anle234b was previously identified as inactive by in vitro experiments. Here, we show that anle234b has no significant effect on the aggregation process in terms of reducing the β-structure content. Moreover, we demonstrate that the hydrogen bonding capabilities are autoinhibited due to steric effects imposed by the molecular geometry of anle234b and thereby indirectly confirm the proposed inhibitory mechanism of anle138b. We anticipate that the prominent binding of anle138b to partially disordered and dynamical aggregate structures is a generic basis for anle138b's ability to suppress toxic oligomer formation in a wide range of amyloidogenic peptides and proteins.
Collapse
Affiliation(s)
- Dirk Matthes
- Computational
Biomolecular Dynamics Group, Department of Theoretical and Computational
Biophysics, Max Planck Institute for Biophysical Chemistry, Am Fassberg
11, 37077 Göttingen, Germany
| | - Vytautas Gapsys
- Computational
Biomolecular Dynamics Group, Department of Theoretical and Computational
Biophysics, Max Planck Institute for Biophysical Chemistry, Am Fassberg
11, 37077 Göttingen, Germany
| | - Christian Griesinger
- Department
of Structural Biology, Max Planck Institute for Biophysical Chemistry, Am Fassberg 11, 37077 Göttingen, Germany
| | - Bert L. de Groot
- Computational
Biomolecular Dynamics Group, Department of Theoretical and Computational
Biophysics, Max Planck Institute for Biophysical Chemistry, Am Fassberg
11, 37077 Göttingen, Germany
| |
Collapse
|
24
|
Cysteine to Serine Conversion at 111th Position Renders the Disaggregation and Retains the Stabilization of Detrimental SOD1 A4V Mutant Against Amyotrophic Lateral Sclerosis in Human-A Discrete Molecular Dynamics Study. Cell Biochem Biophys 2017; 76:231-241. [PMID: 28952073 DOI: 10.1007/s12013-017-0830-5] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2017] [Accepted: 09/18/2017] [Indexed: 12/11/2022]
Abstract
Protein aggregation is a hallmark of various neurodegenerative disorders, such as amyotrophic lateral sclerosis (ALS) in humans. Mutations in Cu/Zn superoxide dismutase (SOD1) protein were found to be a prominent cause behind the majority of the familial ALS cases with abnormal protein aggregates. Herein, we report the biophysical characterization of the beneficial mutation C111S that stabilizes the SOD1 harboring A4V mutation, one of the most lethal diseases causing mutant that leads to protein destabilization and aggregation. In this study, we utilized discrete molecular dynamics (DMD) simulations, which stipulated an outlook over the systematic action of C111S mutation in the A4V mutant that stabilizes the protein and impedes the formation of protein aggregation. Herewith, the findings from our study manifested that the mutation of C111S in SOD1 could aid in regaining the protein structural conformations that protect against the formation of toxic aggregates, thereby hindering the disease pathogenicity subtly. Hence, our study provides a feasible pharmaceutical strategy in developing the treatment for incurable ALS affecting the mankind.
Collapse
|
25
|
Jong K, Grisanti L, Hassanali A. Hydrogen Bond Networks and Hydrophobic Effects in the Amyloid β30–35 Chain in Water: A Molecular Dynamics Study. J Chem Inf Model 2017; 57:1548-1562. [DOI: 10.1021/acs.jcim.7b00085] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- KwangHyok Jong
- Condensed
Matter and Statistical Physics, International Centre for Theoretical Physics, Strada Costiera 11, Trieste 34151, Italy
- SISSA-Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, Trieste 34136, Italy
- Department
of Physics, Kim II Sung University, RyongNam Dong, TaeSong District, Pyongyang, D.P.R., Korea
| | - Luca Grisanti
- Condensed
Matter and Statistical Physics, International Centre for Theoretical Physics, Strada Costiera 11, Trieste 34151, Italy
- SISSA-Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, Trieste 34136, Italy
| | - Ali Hassanali
- Condensed
Matter and Statistical Physics, International Centre for Theoretical Physics, Strada Costiera 11, Trieste 34151, Italy
| |
Collapse
|
26
|
Liao Q, Owen MC, Olubiyi OO, Barz B, Strodel B. Conformational Transitions of the Amyloid-β Peptide Upon Copper(II) Binding and pH Changes. Isr J Chem 2017. [DOI: 10.1002/ijch.201600108] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Qinghua Liao
- Institute of Complex Systems: Structural Biochemistry (ICS-6); Forschungszentrum Jülich GmbH; 52425 Jülich Germany
| | - Michael C. Owen
- Institute of Complex Systems: Structural Biochemistry (ICS-6); Forschungszentrum Jülich GmbH; 52425 Jülich Germany
| | - Olujide O. Olubiyi
- Department of Pharmacology and Therapeutics; College of Medicine and Health Sciences; Afe Babalola University; Nigeria
| | - Bogdan Barz
- Institute of Complex Systems: Structural Biochemistry (ICS-6); Forschungszentrum Jülich GmbH; 52425 Jülich Germany
- Institute of Theoretical and Computational Chemistry; Heinrich Heine University Düsseldorf; 40225 Düsseldorf Germany
| | - Birgit Strodel
- Institute of Complex Systems: Structural Biochemistry (ICS-6); Forschungszentrum Jülich GmbH; 52425 Jülich Germany
- Institute of Theoretical and Computational Chemistry; Heinrich Heine University Düsseldorf; 40225 Düsseldorf Germany
| |
Collapse
|
27
|
Riccardi L, Arencibia JM, Bono L, Armirotti A, Girotto S, De Vivo M. Lid domain plasticity and lipid flexibility modulate enzyme specificity in human monoacylglycerol lipase. Biochim Biophys Acta Mol Cell Biol Lipids 2017; 1862:441-451. [PMID: 28088576 DOI: 10.1016/j.bbalip.2017.01.002] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2016] [Revised: 12/20/2016] [Accepted: 01/09/2017] [Indexed: 12/17/2022]
Abstract
Human monoacylglycerol lipase (MAGL) is a membrane-interacting enzyme that generates pro-inflammatory signaling molecules. For this reason, MAGL inhibition is a promising strategy to treat pain, cancer, and neuroinflammatory diseases. MAGL can hydrolyze monoacylglycerols bearing an acyl chain of different lengths and degrees of unsaturation, cleaving primarily the endocannabinoid 2-arachidonoylglycerol. Importantly, the enzymatic binding site of MAGL is confined by a 75-amino-acid-long, flexible cap domain, named 'lid domain', which is structurally similar to that found in several other lipases. However, it is unclear how lid domain plasticity affects catalysis in MAGL. By integrating extensive molecular dynamics simulations and free-energy calculations with mutagenesis and kinetic experiments, we here define a lid-domain-mediated mechanism for substrate selection and binding in MAGL catalysis. In particular, we clarify the key role of Phe159 and Ile179, two conserved residues within the lid domain, in regulating substrate specificity in MAGL. We conclude by proposing that other structurally related lipases may share this lid-domain-mediated mechanism for substrate specificity.
Collapse
Affiliation(s)
- Laura Riccardi
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy.
| | - Jose M Arencibia
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy.
| | - Luca Bono
- D3-PharmaChemistry, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy.
| | - Andrea Armirotti
- D3-PharmaChemistry, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy.
| | - Stefania Girotto
- CompuNet, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy.
| | - Marco De Vivo
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy; IAS-5/INM-9 Computational Biomedicine Forschungszentrum Jülich, Wilhelm-Johnen-Straße, 52428 Jülich, Germany.
| |
Collapse
|
28
|
Itoh SG, Okumura H. Oligomer Formation of Amyloid-β(29-42) from Its Monomers Using the Hamiltonian Replica-Permutation Molecular Dynamics Simulation. J Phys Chem B 2016; 120:6555-61. [PMID: 27281682 DOI: 10.1021/acs.jpcb.6b03828] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Oligomers of amyloid-β peptides (Aβ) are formed during the early stage of the amyloidogenesis process and exhibit neurotoxicity. The oligomer formation process of Aβ and even that of Aβ fragments are still poorly understood, though understanding of these processes is essential for remedying Alzheimer's disease. In order to better understand the oligomerization process of the C-terminal Aβ fragment Aβ(29-42) at the atomic level, we performed the Hamiltonian replica-permutation molecular dynamics simulation with Aβ(29-42) molecules using the explicit water solvent model. We observed that oligomers increased in size through the sequential addition of monomers to the oligomer, rather than through the assembly of small oligomers. Moreover, solvent effects played an important role in this oligomerization process.
Collapse
Affiliation(s)
- Satoru G Itoh
- Department of Theoretical and Computational Molecular Science, Institute for Molecular Science , Okazaki, Aichi 444-8585, Japan.,Department of Structural Molecular Science, The Graduate University for Advanced Studies , Okazaki, Aichi 444-8585, Japan
| | - Hisashi Okumura
- Department of Theoretical and Computational Molecular Science, Institute for Molecular Science , Okazaki, Aichi 444-8585, Japan.,Department of Structural Molecular Science, The Graduate University for Advanced Studies , Okazaki, Aichi 444-8585, Japan
| |
Collapse
|
29
|
Carballo-Pacheco M, Strodel B. Advances in the Simulation of Protein Aggregation at the Atomistic Scale. J Phys Chem B 2016; 120:2991-9. [PMID: 26965454 DOI: 10.1021/acs.jpcb.6b00059] [Citation(s) in RCA: 91] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Protein aggregation into highly structured amyloid fibrils is associated with various diseases including Alzheimer's disease, Parkinson's disease, and type II diabetes. Amyloids can also have normal biological functions and, in the future, could be used as the basis for novel nanoscale materials. However, a full understanding of the physicochemical forces that drive protein aggregation is still lacking. Such understanding is crucial for the development of drugs that can effectively inhibit aberrant amyloid aggregation and for the directed design of functional amyloids. Atomistic simulations can help understand protein aggregation. In particular, atomistic simulations can be used to study the initial formation of toxic oligomers which are hard to characterize experimentally and to understand the difference in aggregation behavior between different amyloidogenic peptides. Here, we review the latest atomistic simulations of protein aggregation, concentrating on amyloidogenic protein fragments, and provide an outlook for the future in this field.
Collapse
Affiliation(s)
- Martín Carballo-Pacheco
- Institute of Complex Systems: Structural Biochemistry , Forschungszentrum Jülich, 52425 Jülich, Germany.,AICES Graduate School, RWTH Aachen University , Schinkelstraße 2, 52062 Aachen, Germany
| | - Birgit Strodel
- Institute of Complex Systems: Structural Biochemistry , Forschungszentrum Jülich, 52425 Jülich, Germany.,Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf , Universitätsstrasse 1, 40225 Düsseldorf, Germany
| |
Collapse
|
30
|
Carballo-Pacheco M, Ismail AE, Strodel B. Oligomer Formation of Toxic and Functional Amyloid Peptides Studied with Atomistic Simulations. J Phys Chem B 2015; 119:9696-705. [PMID: 26130191 DOI: 10.1021/acs.jpcb.5b04822] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Amyloids are associated with diseases, including Alzheimer's, as well as functional roles such as storage of peptide hormones. It is still unclear what differences exist between aberrant and functional amyloids. However, it is known that soluble oligomers formed during amyloid aggregation are more toxic than the final fibrils. Here, we perform molecular dynamics simulations to study the aggregation of the amyloid-β peptide Aβ25-35, associated with Alzheimer's disease, and two functional amyloid-forming tachykinin peptides: kassinin and neuromedin K. Although the three peptides have similar primary sequences, tachykinin peptides, in contrast to Aβ25-35, form nontoxic amyloids. Our simulations reveal that the charge of the C-terminus is essential to controlling the aggregation process. In particular, when the kassinin C-terminus is not amidated, the aggregation kinetics decreases considerably. In addition, we observe that the monomeric peptides in extended conformations aggregate faster than those in collapsed hairpin-like conformations.
Collapse
Affiliation(s)
- Martín Carballo-Pacheco
- †AICES Graduate School and Aachener Verfahrenstechnik: Molecular Simulations and Transformations, RWTH Aachen University, Schinkelstraße 2, 52062 Aachen, Germany.,‡Institute of Complex Systems: Structural Biochemistry (ICS-6), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Ahmed E Ismail
- †AICES Graduate School and Aachener Verfahrenstechnik: Molecular Simulations and Transformations, RWTH Aachen University, Schinkelstraße 2, 52062 Aachen, Germany
| | - Birgit Strodel
- ‡Institute of Complex Systems: Structural Biochemistry (ICS-6), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany.,¶Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Universitätstraße 1, 40225 Düsseldorf, Germany
| |
Collapse
|
31
|
Yang YI, Gao YQ. Computer Simulation Studies of Aβ37–42 Aggregation Thermodynamics and Kinetics in Water and Salt Solution. J Phys Chem B 2014; 119:662-70. [DOI: 10.1021/jp502169b] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Y. Isaac Yang
- Institute of Theoretical
and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Yi Qin Gao
- Institute of Theoretical
and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
32
|
Smeeton LC, Oakley MT, Johnston RL. Visualizing energy landscapes with metric disconnectivity graphs. J Comput Chem 2014; 35:1481-90. [PMID: 24866379 PMCID: PMC4285870 DOI: 10.1002/jcc.23643] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2013] [Revised: 03/12/2014] [Accepted: 04/14/2014] [Indexed: 11/24/2022]
Abstract
The visualization of multidimensional energy landscapes is important, providing insight into the kinetics and thermodynamics of a system, as well the range of structures a system can adopt. It is, however, highly nontrivial, with the number of dimensions required for a faithful reproduction of the landscape far higher than can be represented in two or three dimensions. Metric disconnectivity graphs provide a possible solution, incorporating the landscape connectivity information present in disconnectivity graphs with structural information in the form of a metric. In this study, we present a new software package, PyConnect, which is capable of producing both disconnectivity graphs and metric disconnectivity graphs in two or three dimensions. We present as a test case the analysis of the 69-bead BLN coarse-grained model protein and show that, by choosing appropriate order parameters, metric disconnectivity graphs can resolve correlations between structural features on the energy landscape with the landscapes energetic and kinetic properties.
Collapse
Affiliation(s)
- Lewis C Smeeton
- School of Chemistry, University of Birmingham, Edgbaston, Birmingham, B15 2TT, United Kingdom
| | | | | |
Collapse
|
33
|
Agbo JK, Gnanasekaran R, Leitner DM. Communication Maps: Exploring Energy Transport through Proteins and Water. Isr J Chem 2014. [DOI: 10.1002/ijch.201300139] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
|
34
|
Nguyen PH, Li MS, Derreumaux P. Amyloid oligomer structure characterization from simulations: A general method. J Chem Phys 2014; 140:094105. [DOI: 10.1063/1.4866902] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
|
35
|
Barz B, Wales DJ, Strodel B. A kinetic approach to the sequence-aggregation relationship in disease-related protein assembly. J Phys Chem B 2014; 118:1003-11. [PMID: 24401100 PMCID: PMC3908877 DOI: 10.1021/jp412648u] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
It is generally accepted that oligomers of aggregating proteins play an important role in the onset of neurodegenerative diseases. While in silico aggregation studies of full length amyloidogenic proteins are computationally expensive, the assembly of short protein fragments derived from these proteins with similar aggregating properties has been extensively studied. In the present work, molecular dynamics simulations are performed to follow peptide aggregation on the microsecond time scale. By defining aggregation states, we identify transition networks, disconnectivity graphs, and first passage time distributions to describe the kinetics of the assembly process. This approach unravels differences in the aggregation into hexamers of two peptides with different primary structures. The first is GNNQQNY, a hydrophilic fragment from the prion protein Sup35, and the second is KLVFFAE, a fragment from amyloid-β protein, with a hydrophobic core delimited by two charged amino acids. The assembly of GNNQQNY suggests a mechanism of monomer addition, with a bias toward parallel peptide pairs and a gradual increase in the amount of β-strand content. For KLVFFAE, a mechanism involving dimers rather than monomers is revealed, involving a generally higher β-strand content and a transition toward a larger number of antiparallel peptide pairs during the rearrangement of the hexamer. The differences observed for the aggregation of the two peptides suggests the existence of a sequence-aggregation relationship.
Collapse
Affiliation(s)
- Bogdan Barz
- Forschungszentrum Jülich GmbH Institute of Complex Systems: Structural Biochemistry (ICS-6), 52425 Jülich, Germany
| | - David J. Wales
- University Chemical Laboratories, Lensfield Road, Cambridge CB2 1EW, UK
| | - Birgit Strodel
- Forschungszentrum Jülich GmbH Institute of Complex Systems: Structural Biochemistry (ICS-6), 52425 Jülich, Germany
- Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| |
Collapse
|
36
|
Jain A, Stock G. Hierarchical Folding Free Energy Landscape of HP35 Revealed by Most Probable Path Clustering. J Phys Chem B 2014; 118:7750-60. [DOI: 10.1021/jp410398a] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Abhinav Jain
- Biomolecular
Dynamics, Institute
of Physics and Freiburg Institute for Advanced Studies (FRIAS), Albert Ludwigs University, 79104 Freiburg, Germany
| | - Gerhard Stock
- Biomolecular
Dynamics, Institute
of Physics and Freiburg Institute for Advanced Studies (FRIAS), Albert Ludwigs University, 79104 Freiburg, Germany
| |
Collapse
|
37
|
Baftizadeh F, Pietrucci F, Biarnés X, Laio A. Nucleation process of a fibril precursor in the C-terminal segment of amyloid-β. PHYSICAL REVIEW LETTERS 2013; 110:168103. [PMID: 23679641 DOI: 10.1103/physrevlett.110.168103] [Citation(s) in RCA: 45] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/04/2012] [Indexed: 06/02/2023]
Abstract
By extended atomistic simulations in explicit solvent and bias-exchange metadynamics, we study the aggregation process of 18 chains of the C-terminal segment of amyloid-β, an intrinsically disordered protein involved in Alzheimer's disease and prone to form fibrils. Starting from a disordered aggregate, we are able to observe the formation of an ordered nucleus rich in beta sheets. The rate limiting step in the nucleation pathway involves crossing a barrier of approximately 40 kcal/mol and is associated with the formation of a very specific interdigitation of the side chains belonging to different sheets. This structural pattern is different from the one observed experimentally in a microcrystal of the same system, indicating that the structure of a "nascent" fibril may differ from the one of an "extended" fibril.
Collapse
|