1
|
Nandi S, Sarkar N. Interactions between Lipid Vesicle Membranes and Single Amino Acid Fibrils: Probable Origin of Specific Neurological Disorders. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024; 40:1971-1987. [PMID: 38240221 DOI: 10.1021/acs.langmuir.3c02429] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/31/2024]
Abstract
Amyloid fibrils are known to be responsible for several neurological disorders, like Alzheimer's disease (AD), Parkinson's disease (PD), etc. For decades, mostly proteins and peptide-based amyloid fibrils have been focused on, and the topic has acknowledged the rise, development, understanding of, and controversy, as well. However, the single amino acid based amyloid fibrils, responsible for several disorders, such as phenylketonuria, tyrosenimia type II, hypermethioninemia, etc., have gotten scientific attention lately. To understand the molecular level pathogenesis of such disorders originated due to the accumulation of single amino acid-based amyloid fibrils, interaction of these fibrils with phospholipid vesicle membranes is found to be an excellent cell-free in vitro setup. Based on such an in vitro setup, these fibrils show a generic mechanism of membrane insertion driven by electrostatic and hydrophobic effects inside the membrane that reduces the integral rigidity of the membrane. Alteration of such fundamental properties of the membrane, therefore, might be referred to as one of the prime pathological factors for the development of these neurological disorders. Hence, such interactions must be investigated in cellular and intracellular compartments to design suitable therapeutic modulators against fibrils.
Collapse
Affiliation(s)
- Sourav Nandi
- Yale School of Medicine, Yale University, New Haven, Connecticut 06510, United States
| | - Nilmoni Sarkar
- Department of Chemistry, Indian Institute of Technology, Kharagpur, 721302, West Bengal, India
| |
Collapse
|
2
|
Okumura H. Perspective for Molecular Dynamics Simulation Studies of Amyloid-β Aggregates. J Phys Chem B 2023; 127:10931-10940. [PMID: 38109338 DOI: 10.1021/acs.jpcb.3c06051] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2023]
Abstract
The cause of Alzheimer's disease is related to aggregates such as oligomers and amyloid fibrils consisting of amyloid-β (Aβ) peptides. Molecular dynamics (MD) simulation studies have been conducted to understand the molecular mechanism of the formation and disruption of Aβ aggregates. In this Perspective, the MD simulation studies are classified into four categories, focusing on the target systems: aggregation of Aβ peptides in bulk solution, Aβ aggregation at the interface, aggregation inhibitor against Aβ peptides, and nonequilibrium MD simulation of Aβ aggregates. MD simulation studies in these categories are first reviewed. Future perspectives in each category are then presented. Finally, the overall perspective is presented on how MD simulations of Aβ aggregates can be utilized for developing Alzheimer's disease treatment.
Collapse
Affiliation(s)
- Hisashi Okumura
- Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan
- Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan
- Graduate Institute for Advanced Studies, SOKENDAI, Okazaki, Aichi 444-8787, Japan
| |
Collapse
|
3
|
Okumura H, Itoh SG, Zen H, Nakamura K. Dissociation process of polyalanine aggregates by free electron laser irradiation. PLoS One 2023; 18:e0291093. [PMID: 37683014 PMCID: PMC10491298 DOI: 10.1371/journal.pone.0291093] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2023] [Accepted: 08/22/2023] [Indexed: 09/10/2023] Open
Abstract
Polyalanine (polyA) disease-causative proteins with an expansion of alanine repeats can be aggregated. Although curative treatments for polyA diseases have not been explored, the dissociation of polyA aggregates likely reduces the cytotoxicity of polyA. Mid-infrared free electron laser (FEL) successfully dissociated multiple aggregates. However, whether the FEL dissociates polyA aggregates like other aggregates has not been tested. Here, we show that FEL at 6.1 μm experimentally weakened the extent of aggregation of a peptide with 13 alanine repeats (13A), and the irradiated 13A exerted lesser cytotoxicity to neuron-like cells than non-irradiated 13A. Then, we applied molecular dynamics (MD) simulation to follow the dissociation process by FEL. We successfully observed how the intermolecular β-sheet of polyA aggregates was dissociated and separated into monomers with helix structures upon FEL irradiation. After the dissociation by FEL, water molecules inhibited the reformation of polyA aggregates. We recently verified the same dissociation process using FEL-treated amyloid-β aggregates. Thus, a common mechanism underlies the dissociation of different protein aggregates that cause different diseases, polyA disease and Alzheimer's disease. However, MD simulation indicated that polyA aggregates are less easily dissociated than amyloid-β aggregates and require longer laser irradiation due to hydrophobic alanine repeats.
Collapse
Affiliation(s)
- Hisashi Okumura
- Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki, Aichi, Japan
- Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Aichi, Japan
- Graduate Institute for Advanced Studies, SOKENDAI, Okazaki, Aichi, Japan
| | - Satoru G Itoh
- Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki, Aichi, Japan
- Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Aichi, Japan
- Graduate Institute for Advanced Studies, SOKENDAI, Okazaki, Aichi, Japan
| | - Heishun Zen
- Institute of Advanced Energy, Kyoto University, Gokasho, Uji, Kyoto, Japan
| | - Kazuhiro Nakamura
- Department of Laboratory Sciences, Gunma University Graduate School of Health Sciences, Maebashi, Gunma, Japan
| |
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
|
Gomari MM, Tarighi P, Choupani E, Abkhiz S, Mohamadzadeh M, Rostami N, Sadroddiny E, Baammi S, Uversky VN, Dokholyan NV. Structural evolution of Delta lineage of SARS-CoV-2. Int J Biol Macromol 2023; 226:1116-1140. [PMID: 36435470 PMCID: PMC9683856 DOI: 10.1016/j.ijbiomac.2022.11.227] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2022] [Revised: 11/19/2022] [Accepted: 11/21/2022] [Indexed: 11/27/2022]
Abstract
One of the main obstacles in prevention and treatment of COVID-19 is the rapid evolution of the SARS-CoV-2 Spike protein. Given that Spike is the main target of common treatments of COVID-19, mutations occurring at this virulent factor can affect the effectiveness of treatments. The B.1.617.2 lineage of SARS-CoV-2, being characterized by many Spike mutations inside and outside of its receptor-binding domain (RBD), shows high infectivity and relative resistance to existing cures. Here, utilizing a wide range of computational biology approaches, such as immunoinformatics, molecular dynamics (MD), analysis of intrinsically disordered regions (IDRs), protein-protein interaction analyses, residue scanning, and free energy calculations, we examine the structural and biological attributes of the B.1.617.2 Spike protein. Furthermore, the antibody design protocol of Rosetta was implemented for evaluation the stability and affinity improvement of the Bamlanivimab (LY-CoV55) antibody, which is not capable of interactions with the B.1.617.2 Spike. We observed that the detected mutations in the Spike of the B1.617.2 variant of concern can cause extensive structural changes compatible with the described variation in immunogenicity, secondary and tertiary structure, oligomerization potency, Furin cleavability, and drug targetability. Compared to the Spike of Wuhan lineage, the B.1.617.2 Spike is more stable and binds to the Angiotensin-converting enzyme 2 (ACE2) with higher affinity.
Collapse
Affiliation(s)
- Mohammad Mahmoudi Gomari
- Student Research Committee, Iran University of Medical Sciences, Tehran 1449614535, Iran,Department of Medical Biotechnology, Faculty of Allied Medicine, Iran University of Medical Sciences, Tehran 1449614535, Iran
| | - Parastoo Tarighi
- Department of Medical Biotechnology, Faculty of Allied Medicine, Iran University of Medical Sciences, Tehran 1449614535, Iran
| | - Edris Choupani
- Department of Medical Biotechnology, Faculty of Allied Medicine, Iran University of Medical Sciences, Tehran 1449614535, Iran
| | - Shadi Abkhiz
- Department of Medical Biotechnology, Faculty of Allied Medicine, Iran University of Medical Sciences, Tehran 1449614535, Iran
| | - Masoud Mohamadzadeh
- Department of Chemistry, Faculty of Sciences, University of Hormozgan, Bandar Abbas 7916193145, Iran
| | - Neda Rostami
- Department of Chemical Engineering, Faculty of Engineering, Arak University, Arak 3848177584, Iran
| | - Esmaeil Sadroddiny
- Medical Biotechnology Department, School of Advanced Technologies in Medicine, Tehran University of Medical Sciences, Tehran 1417613151, Iran
| | - Soukayna Baammi
- African Genome Centre (AGC), Mohammed VI Polytechnic University, Benguerir 43150, Morocco
| | - Vladimir N. Uversky
- Department of Molecular Medicine and USF Health Byrd Alzheimer's Research Institute, Morsani College of Medicine, University of South Florida, Tampa, FL 33620, USA,Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, 141700 Dolgoprudny, Russia,Correspondence to: V.N. Uversky, Department of Molecular Medicine and USF Health Byrd Alzheimer's Research Institute, Morsani College of Medicine, University of South Florida, Tampa, FL 33620, USA
| | - Nikolay V. Dokholyan
- Department of Pharmacology, Department of Biochemistry & Molecular Biology, Pennsylvania State University College of Medicine, Hershey, PA 16802, USA,Corresponding author
| |
Collapse
|
6
|
Shao L, Ma J, Prelesnik JL, Zhou Y, Nguyen M, Zhao M, Jenekhe SA, Kalinin SV, Ferguson AL, Pfaendtner J, Mundy CJ, De Yoreo JJ, Baneyx F, Chen CL. Hierarchical Materials from High Information Content Macromolecular Building Blocks: Construction, Dynamic Interventions, and Prediction. Chem Rev 2022; 122:17397-17478. [PMID: 36260695 DOI: 10.1021/acs.chemrev.2c00220] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Hierarchical materials that exhibit order over multiple length scales are ubiquitous in nature. Because hierarchy gives rise to unique properties and functions, many have sought inspiration from nature when designing and fabricating hierarchical matter. More and more, however, nature's own high-information content building blocks, proteins, peptides, and peptidomimetics, are being coopted to build hierarchy because the information that determines structure, function, and interfacial interactions can be readily encoded in these versatile macromolecules. Here, we take stock of recent progress in the rational design and characterization of hierarchical materials produced from high-information content blocks with a focus on stimuli-responsive and "smart" architectures. We also review advances in the use of computational simulations and data-driven predictions to shed light on how the side chain chemistry and conformational flexibility of macromolecular blocks drive the emergence of order and the acquisition of hierarchy and also on how ionic, solvent, and surface effects influence the outcomes of assembly. Continued progress in the above areas will ultimately usher in an era where an understanding of designed interactions, surface effects, and solution conditions can be harnessed to achieve predictive materials synthesis across scale and drive emergent phenomena in the self-assembly and reconfiguration of high-information content building blocks.
Collapse
Affiliation(s)
- Li Shao
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Jinrong Ma
- Molecular Engineering and Sciences Institute, University of Washington, Seattle, Washington 98195, United States
| | - Jesse L Prelesnik
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Yicheng Zhou
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Mary Nguyen
- Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, United States.,Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Mingfei Zhao
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Samson A Jenekhe
- Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, United States.,Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Sergei V Kalinin
- Department of Materials Science and Engineering, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Andrew L Ferguson
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Jim Pfaendtner
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States.,Materials Science and Engineering, University of Washington, Seattle, Washington 98195, United States
| | - Christopher J Mundy
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States.,Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, United States
| | - James J De Yoreo
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States.,Materials Science and Engineering, University of Washington, Seattle, Washington 98195, United States
| | - François Baneyx
- Molecular Engineering and Sciences Institute, University of Washington, Seattle, Washington 98195, United States.,Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, United States
| | - Chun-Long Chen
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States.,Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
7
|
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
|
8
|
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
|
9
|
All-Atom Molecular Dynamics Simulation Methods for the Aggregation of Protein and Peptides: Replica Exchange/Permutation and Nonequilibrium Simulations. METHODS IN MOLECULAR BIOLOGY (CLIFTON, N.J.) 2022; 2340:197-220. [PMID: 35167076 DOI: 10.1007/978-1-0716-1546-1_10] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Protein aggregates are associated with more than 40 serious human diseases. To understand the formation mechanism of protein aggregates at atomic level, all-atom molecular dynamics (MD) simulation is a powerful computational tool. In this chapter, we review the all-atom MD simulation methods that are useful for study on the protein aggregation. We first explain conventional MD simulation methods in physical statistical ensembles, such as the canonical and isothermal-isobaric ensembles. We then describe the generalized-ensemble algorithms such as replica-exchange and replica-permutation MD methods. These methods can overcome a difficulty, in which simulations tend to get trapped in local-minimum free-energy states. Finally we explain the nonequilibrium MD method. Some simulation results based on these methods are also presented.
Collapse
|
10
|
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
|
11
|
Okumura H, Itoh SG, Nakamura K, Kawasaki T. Role of Water Molecules and Helix Structure Stabilization in the Laser-Induced Disruption of Amyloid Fibrils Observed by Nonequilibrium Molecular Dynamics Simulations. J Phys Chem B 2021; 125:4964-4976. [PMID: 33961416 DOI: 10.1021/acs.jpcb.0c11491] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Water plays a crucial role in the formation and destruction of biomolecular structures. The mechanism for destroying biomolecular structures was thought to be an active breaking of hydrogen bonds by water molecules. However, using nonequilibrium molecular dynamics simulations, in which an amyloid-β amyloid fibril was destroyed via infrared free-electron laser (IR-FEL) irradiation, we discovered a new mechanism, in which water molecules disrupt protein aggregates. The intermolecular hydrogen bonds formed by C═O and N-H in the fibril are broken at each pulse of laser irradiation. These bonds spontaneously re-form after the irradiation in many cases. However, when a water molecule happens to enter the gap between C═O and N-H, it inhibits the re-formation of the hydrogen bonds. Such sites become defects in the regularly aligned hydrogen bonds, from which all hydrogen bonds in the intermolecular β-sheet are broken as the fraying spreads. This role of water molecules is entirely different from other known mechanisms. This new mechanism can explain the recent experiments showing that the amyloid fibrils are not destroyed by laser irradiation under dry conditions. Additionally, we found that helix structures form more after the amyloid disruption; this is because the resonance frequency is different in a helix structure. Our findings provide a theoretical basis for the application of IR-FEL to the future treatment of amyloidosis.
Collapse
Affiliation(s)
- Hisashi Okumura
- Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan.,Institute for Molecular Science, 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-8787, Japan
| | - Satoru G Itoh
- Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan.,Institute for Molecular Science, 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-8787, Japan
| | - Kazuhiro Nakamura
- Department of Laboratory Sciences, Graduate School of Health Sciences, Gunma University, Maebashi, Gunma 371-8514, Japan
| | - Takayasu Kawasaki
- IR Free Electron Laser Research Center, Research Institute for Science and Technology, Organization for Research Advancement, Tokyo University of Science, Noda, Chiba 278-8510, Japan
| |
Collapse
|
12
|
Yamauchi M, Okumura H. Dimerization of α-Synuclein Fragments Studied by Isothermal-Isobaric Replica-Permutation Molecular Dynamics Simulation. J Chem Inf Model 2021; 61:1307-1321. [PMID: 33625841 DOI: 10.1021/acs.jcim.0c01056] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Aggregates and fibrils of intrinsically disordered α-synuclein are associated with Parkinson's disease. Within a non-amyloid β component (NAC) spanning from the 61st to the 95th residue of α-synuclein, an 11-residue segment called NACore (68GAVVTGVTAVA78) is an essential region for both fibril formation and cytotoxicity. Although NACore peptides alone are known to form aggregates and amyloid fibrils, the mechanisms of aggregation and fibrillation remain unknown. This study investigated the dimerization process of NACore peptides as the initial stage of the aggregation and fibrillation processes. We performed an isothermal-isobaric replica-permutation molecular dynamics simulation, which is one of the efficient sampling methods, for the two NACore peptides in explicit water over 96 μs. The simulation succeeded in sampling a variety of dimer structures. An analysis of secondary structure revealed that most of the NACore dimers form intermolecular β-bridges. In particular, more antiparallel β-bridges were observed than parallel β-bridges. We also found that intramolecular secondary structures such as α-helix and antiparallel β-bridge are stabilized in the pre-dimer state. However, we identified that the intermolecular β-bridges tend to form directly between residues with no specific structure rather than via the intramolecular β-bridges. This is because the NACore peptides still have a low propensity to form the intramolecular secondary structures even though they are stabilized in the pre-dimer state.
Collapse
Affiliation(s)
- Masataka Yamauchi
- Department of Structural Molecular Science, The Graduate University for Advanced Studies(SOKENDAI), Okazaki, Aichi 444-8787, Japan.,Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan.,Institute for Molecular Science (IMS), National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan
| | - Hisashi Okumura
- Department of Structural Molecular Science, The Graduate University for Advanced Studies(SOKENDAI), Okazaki, Aichi 444-8787, Japan.,Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan.,Institute for Molecular Science (IMS), 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
|
Ding F, Peng W, Peng YK, Liu BQ. Elucidating the potential neurotoxicity of chiral phenthoate: Molecular insight from experimental and computational studies. CHEMOSPHERE 2020; 255:127007. [PMID: 32416396 DOI: 10.1016/j.chemosphere.2020.127007] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/08/2020] [Revised: 05/05/2020] [Accepted: 05/06/2020] [Indexed: 06/11/2023]
Abstract
Chiral organophosphorus pollutants are existed ubiquitously in the ecological environment, but the enantioselective toxicities of these nerve agents to humans and their molecular bases have not been fully elucidated. Using experimental and computational approaches, this story was to explore the neurotoxic response process of the target acetylcholinesterase (AChE) to chiral phenthoate and further decipher the microscopic mechanism of such toxicological effect at the enantiomeric level. The results showed that the toxic reaction of AChE with chiral phenthoate exhibited significant enantioselectivity, and (R)-phenthoate (K=1.486 × 105 M-1) has a bioaffinity for the nerve enzyme nearly three times that of (S)-phenthoate (K=4.503 × 104 M-1). Dynamic research outcomes interpreted the wet experiments, and the inherent conformational flexibility of the target enzyme has a great influence on the enantioselective neurotoxicological action processes, especially reflected in the conformational changes of the three key loop regions (i.e. residues His-447, Gly-448, and Tyr-449; residues Gly-122, Phe-123, and Tyr-124; and residues Thr-75, Leu-76, and Tyr-77) around the reaction patch. This was supported by the quantitative results of conformational studies derived from circular dichroism spectroscopy (α-helix: 34.7%→30.2%/31.6%; β-sheet: 23.6%→19.5%/20.7%; turn: 19.2%→22.4%/21.9%; and random coil: 22.5%→27.9%/25.8%). Meanwhile, via analyzing the modes of toxic action and free energies, we can find that (R)-phenthoate has a strong inhibitory effect on the enzymatic activity of AChE, as compared with (S)-phenthoate, and electrostatic energy (-23.79/-17.77 kJ mol-1) played a critical role in toxicological reactions. These points were the underlying causes of chiral phenthoate displaying different degrees of enantioselective neurotoxicity.
Collapse
Affiliation(s)
- Fei Ding
- Department of Environmental Science and Engineering, School of Water and Environment, Chang'an University, Xi'an, 710054, China; Key Laboratory of Subsurface Hydrology and Ecological Effect in Arid Region of Ministry of Education, Chang'an University, No. 126 Yanta Road, Yanta District, Xi'an, 710054, China
| | - Wei Peng
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, 361005, China.
| | - Yu-Kui Peng
- Center for Food Quality Supervision, Inspection & Testing, Ministry of Agriculture and Rural Affairs, Northwest A&F University, Yangling, 712100, China
| | - Bing-Qi Liu
- Department of Agricultural Chemistry, Qingdao Agricultural University, Qingdao, 266109, China
| |
Collapse
|
15
|
2019-A year in Biophysical Reviews. Biophys Rev 2019; 11:833-839. [PMID: 31741173 DOI: 10.1007/s12551-019-00607-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2019] [Accepted: 11/07/2019] [Indexed: 02/07/2023] Open
|
16
|
Yamauchi M, Okumura H. Replica sub-permutation method for molecular dynamics and monte carlo simulations. J Comput Chem 2019; 40:2694-2711. [PMID: 31365132 DOI: 10.1002/jcc.26030] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2019] [Revised: 06/03/2019] [Accepted: 06/21/2019] [Indexed: 11/12/2022]
Abstract
We propose an improvement of the replica-exchange and replica-permutation methods, which we call the replica sub-permutation method (RSPM). Instead of considering all permutations, this method uses a new algorithm referred to as sub-permutation to perform parameter transition. The RSPM succeeds in reducing the number of combinations between replicas and parameters without the loss of sampling efficiency. For comparison, we applied the replica sub-permutation, replica-permutation, and replica-exchange methods to a β-hairpin mini protein, chignolin, in explicit water. We calculated the transition ratio and number of tunneling events in the parameter space, the number of folding-unfolding events, the autocorrelation function, and the autocorrelation time as measures of sampling efficiency. The results indicate that among the three methods, the proposed RSPM is the most efficient in both parameter and conformational spaces. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Masataka Yamauchi
- Department of Structural Molecular Science, SOKENDAI (The Graduate University for Advanced Studies), Okazaki, Aichi, 444-8585, Japan.,Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki, Aichi, 444-8787, Japan.,Institute for Molecular Science (IMS), National Institutes of Natural Sciences, Okazaki, Aichi, 444-8585, Japan
| | - Hisashi Okumura
- Department of Structural Molecular Science, SOKENDAI (The Graduate University for Advanced Studies), Okazaki, Aichi, 444-8585, Japan.,Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki, Aichi, 444-8787, Japan.,Institute for Molecular Science (IMS), National Institutes of Natural Sciences, Okazaki, Aichi, 444-8585, Japan
| |
Collapse
|
17
|
Battle AR, Norton RS, Böcking T, Noji H, Kim KK, Nagayama K. Editorial: Special issue of Biophysical Reviews dedicated to the joint 10th Asian Biophysics Association Symposium and 42nd Australian Society for Biophysics Meeting, Melbourne, Australia, December 2-6, 2018. Biophys Rev 2019; 11:245-247. [PMID: 31115863 PMCID: PMC6557946 DOI: 10.1007/s12551-019-00553-y] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2019] [Accepted: 05/07/2019] [Indexed: 12/15/2022] Open
Affiliation(s)
- Andrew R Battle
- School of Biomedical Sciences, Queensland University of Technology, Brisbane, 4000, Australia.
- Translational Research Institute and Institute for Biomedical Innovation (QUT), 37 Kent Street, Woolloongabba, 4102, Australia.
- The University of Queensland Diamantina Institute, Faculty of Medicine, The University of Queensland, Brisbane, 4102, Australia.
| | - Raymond S Norton
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University (Parkville Campus), 381 Royal Parade, Parkville, VIC, 3052, Australia
| | - Till Böcking
- EMBL Australia Node in Single Molecule Science, School of Medical Sciences, UNSW Sydney, Sydney, NSW, Australia
| | - Hiroyuki Noji
- Department of Applied Chemistry, School of Engineering, The University of Tokyo, Tokyo, 113-8656, Japan
| | - Kyeong Kyu Kim
- Department of Molecular Cell Biology, School of Medicine, Sungkyunkwan University, Suwon, 16419, South Korea
| | - Kuniaki Nagayama
- National Institute for Physiological Sciences, Myodaiji-cho, Okazaki, 444-8585, Japan
| |
Collapse
|