1
|
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
|
2
|
Unpolarized laser method for infrared spectrum calculation of amide I CO bonds in proteins using molecular dynamics simulation. Comput Biol Med 2023; 159:106902. [PMID: 37086661 PMCID: PMC10186340 DOI: 10.1016/j.compbiomed.2023.106902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2022] [Revised: 03/07/2023] [Accepted: 04/09/2023] [Indexed: 04/24/2023]
Abstract
The investigation of the strong infrared (IR)-active amide I modes of peptides and proteins has received considerable attention because a wealth of detailed information on hydrogen bonding, dipole-dipole interactions, and the conformations of the peptide backbone can be derived from the amide I bands. The interpretation of experimental spectra typically requires substantial theoretical support, such as direct ab-initio molecular dynamics simulation or mixed quantum-classical description. However, considering the difficulties associated with these theoretical methods and their applications are limited in small peptides, it is highly desirable to develop a simple yet efficient approach for simulating the amide I modes of any large proteins in solution. In this work, we proposed a comprehensive computational method that extends the well-established molecular dynamics (MD) simulation method to include an unpolarized IR laser for exciting the CO bonds of proteins. We showed the amide I frequency corresponding to the frequency of the laser pulse which resonated with the CO bond vibration. At this frequency, the protein energy and the CO bond length fluctuation were maximized. Overall, the amide I bands of various single proteins and amyloids agreed well with experimental data. The method has been implemented into the AMBER simulation package, making it widely available to the scientific community. Additionally, the application of the method to simulate the transient amide I bands of amyloid fibrils during the IR laser-induced disassembly process was discussed in details.
Collapse
|
3
|
Disassembly of Amyloid Fibril with Infrared Free Electron Laser. Int J Mol Sci 2023; 24:ijms24043686. [PMID: 36835098 PMCID: PMC9967569 DOI: 10.3390/ijms24043686] [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: 12/19/2022] [Revised: 02/03/2023] [Accepted: 02/09/2023] [Indexed: 02/15/2023] Open
Abstract
Amyloid fibril causes serious amyloidosis such as neurodegenerative diseases. The structure is composed of rigid β-sheet stacking conformation which makes it hard to disassemble the fibril state without denaturants. Infrared free electron laser (IR-FEL) is an intense picosecond pulsed laser that is oscillated through a linear accelerator, and the oscillation wavelengths are tunable from 3 μm to 100 μm. Many biological and organic compounds can be structurally altered by the mode-selective vibrational excitations due to the wavelength variability and the high-power oscillation energy (10-50 mJ/cm2). We have found that several different kinds of amyloid fibrils in amino acid sequences were commonly disassembled by the irradiation tuned to amide I (6.1-6.2 μm) where the abundance of β-sheet decreased while that of α-helix increased by the vibrational excitation of amide bonds. In this review, we would like to introduce the IR-FEL oscillation system briefly and describe combination studies of experiments and molecular dynamics simulations on disassembling amyloid fibrils of a short peptide (GNNQQNY) from yeast prion and 11-residue peptide (NFLNCYVSGFH) from β2-microglobulin as representative models. Finally, possible applications of IR-FEL for amyloid research can be proposed as a future outlook.
Collapse
|
4
|
Probing protein misfolding and dissociation with an infrared free-electron laser. Methods Enzymol 2022; 679:65-96. [PMID: 36682873 DOI: 10.1016/bs.mie.2022.08.047] [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: 12/24/2022]
Abstract
Misfolding is observed in the mutant proteins that are causative for neurodegenerative disorders such as polyglutamine diseases. These proteins are prone to aggregate in the cytoplasm and nucleus of cells. To reproduce cells with the aggregated proteins, gene expression system is usually applied, in which the expression construct having the mutated DNA sequence of the interest is transfected into cells. The transfected DNA is finally converted into the mutant protein, which is gradually aggregated in the cells. In addition, a simple method to prepare the cells having aggregates inside has been recently applied. Peptides were first aggregated by incubating them in water. The aggregates are spontaneously taken up by cells because aggregated proteins generally transfer between cells. Peptides with different degrees of aggregation can be made by changing the incubation times and temperatures, which enables to examine contribution of aggregation to the toxicity to the recipient cells. Moreover, such cells can be used for therapeutic researches of diseases in which aggregates are involved. In this chapter, we show methods to induce aggregation of peptides. The functional analyses of the cells with aggregates are also described. Then, experimental dissociation of the aggregates produced using this method by mid infrared free electron laser irradiation and its theoretical support by molecular dynamics simulation are introduced as the therapeutic research for neurodegenerative disorders.
Collapse
|
5
|
Application study of infrared free-electron lasers towards the development of amyloidosis therapy. JOURNAL OF SYNCHROTRON RADIATION 2022; 29:1133-1140. [PMID: 36073871 PMCID: PMC9455209 DOI: 10.1107/s1600577522007330] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/24/2022] [Accepted: 07/15/2022] [Indexed: 05/27/2023]
Abstract
Amyloidosis is known to be caused by the deposition of amyloid fibrils into various biological tissues; effective treatments for the disease are little established today. An infrared free-electron laser (IR-FEL) is an accelerator-based picosecond-pulse laser having tunable infrared wavelengths. In the current study, the irradiation effect of an IR-FEL was tested on an 11-residue peptide (NFLNCYVSGFH) fibril from β2-microglobulin (β2M) with the aim of applying IR-FELs to amyloidosis therapy. Infrared microspectroscopy (IRM) and scanning electron microscopy showed that a fibril of β2M peptide was clearly dissociated by IR-FEL at 6.1 µm (amide I) accompanied by a decrease of the β-sheet and an increase of the α-helix. No dissociative process was recognized at 6.5 µm (amide II) as well as at 5.0 µm (non-specific wavelength). Equilibrium molecular dynamics simulations indicated that the α-helix can exist stably and the probability of forming interchain hydrogen bonds associated with the internal asparagine residue (N4) is notably reduced compared with other amino acids after the β-sheet is dissociated by amide I specific irradiation. This result implies that N4 plays a key role for recombination of hydrogen bonds in the dissociation of the β2M fibril. In addition, the β-sheet was disrupted at temperatures higher than 340 K while the α-helix did not appear even though the fibril was heated up to 363 K as revealed by IRM. The current study gives solid evidence for the laser-mediated conversion from β-sheet to α-helix in amyloid fibrils at the molecular level.
Collapse
|
6
|
Stable Cavitation Interferes with Aβ 16-22 Oligomerization. J Chem Inf Model 2022; 62:3885-3895. [PMID: 35920625 DOI: 10.1021/acs.jcim.2c00764] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Ultrasound and microbubbles are used for many medical applications nowadays. Scanning ultrasound can remove amyloid-β (Aβ) aggregates in the mouse brain and restores memory in an Alzheimer's disease mouse model. In vitro studies showed that amyloid fibrils are fragmented due to the ultrasound-induced bubble inertial cavitation, and ultrasonic pulses accelerate the depolymerization of Aβ fibrils into monomers at 1 μM of concentration. Under applied ultrasound, microbubbles can be in a stable oscillating state or unstable inertial cavitation state. The latter occurs when ultrasound causes a dramatic change of bubble sizes above a certain acoustic pressure. We have developed and implemented a nonequilibrium molecular dynamics simulation algorithm to the AMBER package, to facilitate the investigation of the molecular mechanism of Aβ oligomerization under stable cavitation. Our results indicated that stable cavitation not only inhibited oligomeric formation, but also prevented the formation of β-rich oligomers. The network analysis of state transitions revealed that stable cavitation altered the oligomerization pathways of Aβ16-22 peptides. Our simulation tool may be applied to optimize the experimental conditions to achieve the best therapeutical effect.
Collapse
|
7
|
Abstract
Protein aggregation can lead to well-defined structures that are functional, but is also the cause of the death of neuron cells in many neurodegenerative diseases. The complexity of the molecular events involved in the aggregation kinetics of amyloid proteins and the transient and heterogeneous characters of all oligomers prevent high-resolution structural experiments. As a result, computer simulations have been used to determine the atomic structures of amyloid proteins at different association stages as well as to understand fibril dissociation. In this chapter, we first review the current computer simulation methods used for aggregation with some atomistic and coarse-grained results aimed at better characterizing the early formed oligomers and amyloid fibril formation. Then we present the applications of non-equilibrium molecular dynamics simulations to comprehend the dissociation of protein assemblies.
Collapse
|
8
|
Dynamics of Amyloid Formation from Simplified Representation to Atomistic Simulations. Methods Mol Biol 2022; 2405:95-113. [PMID: 35298810 DOI: 10.1007/978-1-0716-1855-4_5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Amyloid fibril formation is an intrinsic property of short peptides, non-disease proteins, and proteins associated with neurodegenerative diseases. Aggregates of the Aβ and tau proteins, the α-synuclein protein, and the prion protein are observed in the brain of Alzheimer's, Parkinson's, and prion disease patients, respectively. Due to the transient short-range and long-range interactions of all species and their high aggregation propensities, the conformational ensemble of these devastating proteins, the exception being for the monomeric prion protein, remains elusive by standard structural biology methods in bulk solution and in lipid membranes. To overcome these limitations, an increasing number of simulations using different sampling methods and protein models have been performed. In this chapter, we first review our main contributions to the field of amyloid protein simulations aimed at understanding the early aggregation steps of short linear amyloid peptides, the conformational ensemble of the Aβ40/42 dimers in bulk solution, and the stability of Aβ aggregates in lipid membrane models. Then we focus on our studies on the interactions of amyloid peptides/inhibitors to prevent aggregation, and long amyloid sequences, including new results on a monomeric tau construct.
Collapse
|
9
|
Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design. Pharmaceuticals (Basel) 2021; 14:1062. [PMID: 34681286 PMCID: PMC8537670 DOI: 10.3390/ph14101062] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2021] [Revised: 10/14/2021] [Accepted: 10/15/2021] [Indexed: 11/17/2022] Open
Abstract
We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.
Collapse
|
10
|
Kinetic and thermodynamic stability comparison for the fibrillar form of small amyloid-β(1-42) oligomers using scaled molecular dynamics. Phys Chem Chem Phys 2021; 23:16897-16908. [PMID: 34328153 DOI: 10.1039/d1cp01866c] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Amyloid-β (Aβ) oligomers act as intermediates for several neurodegenerative disease-relevant fibril formations. However, gaining insight into the oligomer to fibril conversion process remains a challenge due to the transient nature of small Aβ. In this study, we probe the kinetic and thermodynamic stabilities of small Aβ(1-42) oligomers in fibrillar conformations to understand from what size these aggregates start forming stable fibrils. With no definite structures available for small Aβ42 aggregates, we have started with oligomers extracted from mature fibrils having four, five, six and nine chains stacked together, and have performed order-to-disorder transition on these systems. Using scaled molecular dynamics (sMD) simulation, the timescale for breaking the native contacts of fibrils has been compared. The results indicate that the kinetic stability of oligomers increases with size, especially at the C-terminus end beyond five-chain oligomers. The free energy of breaking the contacts at the β-sheet regions in the structures has been obtained on an unscaled potential from a free energy extrapolation (FEE) approach. The values show that although stable minima are obtained for larger oligomers due to the enhanced stability of the C-terminus ends, fully stable fibril formation may require aggregates larger than the ones considered in our study. Additionally, dissimilar kinetics for the unbinding of terminal chains across all the oligomers has been observed. The interaction energy values calculated from unscaled MD simulations reveal the crucial role of water in our observations. Our work provides the application of an easy-to-deploy method that sheds light on interactions which could be significant in the early stages of Aβ42 fibril formation.
Collapse
|
11
|
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: 12] [Impact Index Per Article: 4.0] [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
|
12
|
Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis. Chem Rev 2021; 121:2545-2647. [PMID: 33543942 PMCID: PMC8836097 DOI: 10.1021/acs.chemrev.0c01122] [Citation(s) in RCA: 355] [Impact Index Per Article: 118.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.
Collapse
|
13
|
Application of mid-infrared free-electron laser for structural analysis of biological materials. JOURNAL OF SYNCHROTRON RADIATION 2021; 28:28-35. [PMID: 33399549 DOI: 10.1107/s160057752001406x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/05/2020] [Accepted: 10/21/2020] [Indexed: 06/12/2023]
Abstract
A mid-infrared free-electron laser (MIR-FEL) is a synchrotron-radiation-based femto- to pico-second pulse laser. It has unique characteristics such as variable wavelengths in the infrared region and an intense pulse energy. So far, MIR-FELs have been utilized to perform multi-photon absorption reactions against various gas molecules and protein aggregates in physical chemistry and biomedical fields. However, the applicability of MIR-FELs for the structural analysis of solid materials is not well recognized in the analytical field. In the current study, an MIR-FEL is applied for the first time to analyse the internal structure of biological materials by using fossilized inks from cephalopods as the model sample. Two kinds of fossilized inks that were collected from different strata were irradiated at the dry state by tuning the oscillation wavelengths of the MIR-FEL to the phosphoryl stretching mode of hydroxyapatite (9.6 µm) and to the carbonyl stretching mode of melanin (5.8 µm), and the subsequent structural changes in those materials were observed by using infrared microscopy and far-infrared spectroscopy. The structural variation of these biological fossils is discussed based on the infrared-absorption spectral changes that were enhanced by the MIR-FEL irradiation, and the potential use of MIR-FELs for the structural evaluation of biomaterials is suggested.
Collapse
|
14
|
Nonequilibrium molecular dynamics simulations of infrared laser-induced dissociation of a tetrameric Aβ42 β-barrel in a neuronal membrane model. Chem Phys Lipids 2020; 234:105030. [PMID: 33347835 DOI: 10.1016/j.chemphyslip.2020.105030] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2020] [Revised: 12/04/2020] [Accepted: 12/10/2020] [Indexed: 11/29/2022]
Abstract
Experimental studies have reported that the amyloid-β proteins can form pores in cell membranes, and this could be one possible source of toxicity in Alzheimer's disease. Dissociation of these pores could therefore be a potential therapeutic approach. It is known that high photon density free-electron laser experiments and laser-induced nonequilibrium molecular dynamics simulations (NEMD) can dissociate amyloid fibrils at specific frequencies in vitro. Our question is whether NEMD simulations can dissociate amyloid pores in a bilayer mimicking a neuronal membrane, and as an example, we select a tetrameric Aβ42 β-barrel. Our simulations shows that the resonance between the laser field and the amide I vibrational mode of the barrel destabilises all intramolecular and intermolecular hydrogen bonds of Aβ42 and converts the β-barrel to a random/coil disordered oligomer. Starting from this disordered oligomer, extensive standard MD simulations shows sampling of disordered Aβ42 states without any increase of β-sheet and reports that the orientational order of lipids is minimally disturbed. Interestingly, the frequency to be employed to dissociate this beta-barrel is specific to the amino acid sequence. Taken together with our previous simulation results, this study indicates that infrared laser irradiation can dissociate amyloid fibrils and oligomers in bulk solution and in a membrane environment without affecting the surrounding molecules, offering therefore a promising way to retard the progression of AD.
Collapse
|
15
|
Irradiation effect of a submillimeter wave from 420 GHz gyrotron on amyloid peptides in vitro. BIOMEDICAL OPTICS EXPRESS 2020; 11:5341-5351. [PMID: 33014618 PMCID: PMC7510884 DOI: 10.1364/boe.395218] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/22/2020] [Revised: 08/07/2020] [Accepted: 08/07/2020] [Indexed: 06/11/2023]
Abstract
On using the far-infrared radiation system, whether the irradiation effect is thermal or non-thermal is controversial. We irradiated amyloid peptides that are causal factors for amyloidosis by using a submillimeter wave from 420 GHz gyrotron. Fluorescence reagent assay, optical and electron microscopies, and synchrotron-radiation infrared microscopy showed that the irradiation increased the fibrous conformation of peptides at room temperature for 30 min. The temperature increase on the sample was only below 5 K, and a simple heating up to 318 K hardly induced the fibril formation. Therefore, the amyloid aggregation was driven by the far-infrared radiation with little thermal effect.
Collapse
|