Jessica Nasica-Labouze

Amyloid β Protein and Alzheimer’s Disease: When Computer Simulations Complement Experimental Studies

Simulations Complement Experimental Studies Jessica Nasica-Labouze,† Phuong H. Nguyen,† Fabio Sterpone,† Olivia Berthoumieu,‡ Nicolae-Viorel Buchete, Sebastien Cote, Alfonso De Simone, Andrew J. Doig, Peter Faller,‡ Angel Garcia, Alessandro Laio, Mai Suan Li, Simone Melchionna, Normand Mousseau, Yuguang Mu, Anant Paravastu, Samuela Pasquali,† David J. Rosenman, Birgit Strodel, Bogdan Tarus,† John H. Viles, Tong Zhang,†,▲ Chunyu Wang, and Philippe Derreumaux*,†,□ †Laboratoire de Biochimie Theorique, Institut de Biologie Physico-Chimique (IBPC), UPR9080 CNRS, Universite Paris Diderot, Sorbonne Paris Cite, 13 rue Pierre et Marie Curie, 75005 Paris, France ‡LCC (Laboratoire de Chimie de Coordination), CNRS, Universite de Toulouse, Universite Paul Sabatier (UPS), Institut National Polytechnique de Toulouse (INPT), 205 route de Narbonne, BP 44099, Toulouse F-31077 Cedex 4, France School of Physics & Complex and Adaptive Systems Laboratory, University College Dublin, Belfield, Dublin 4, Ireland Deṕartement de Physique and Groupe de recherche sur les proteines membranaires (GEPROM), Universite de Montreal, C.P. 6128, succursale Centre-ville, Montreal, Quebec H3C 3T5, Canada Department of Life Sciences, Imperial College London, London SW7 2AZ, United Kingdom Manchester Institute of Biotechnology, University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom Department of Physics, Applied Physics, & Astronomy, and Department of Biology, Rensselaer Polytechnic Institute, Troy, New York 12180, United States The International School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-668 Warsaw, Poland Institute for Computational Science and Technology, SBI Building, Quang Trung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam Instituto Processi Chimico-Fisici, CNR-IPCF, Consiglio Nazionale delle Ricerche, 00185 Roma, Italy School of Biological Sciences, Nanyang Technological University, 60 Nanyang Drive, 637551 Singapore Department of Chemical and Biomedical Engineering, Florida A&M University-Florida State University (FAMU-FSU) College of Engineering, 2525 Pottsdamer Street, Tallahassee, Florida 32310, United States National High Magnetic Field Laboratory, 1800 East Paul Dirac Drive, Tallahassee, Florida 32310, United States Institute of Complex Systems: Structural Biochemistry (ICS-6), Forschungszentrum Julich GmbH, 52425 Julich, Germany School of Biological and Chemical Sciences, Queen Mary University of London, London E1 4NS, United Kingdom Institut Universitaire de France, 75005 Paris, France

A Multiscale Approach to Characterize the Early Aggregation Steps of the Amyloid-Forming Peptide GNNQQNY from the Yeast Prion Sup-35

The self-organization of peptides into amyloidogenic oligomers is one of the key events for a wide range of molecular and degenerative diseases. Atomic-resolution characterization of the mechanisms responsible for the aggregation process and the resulting structures is thus a necessary step to improve our understanding of the determinants of these pathologies. To address this issue, we combine the accelerated sampling properties of replica exchange molecular dynamics simulations based on the OPEP coarse-grained potential with the atomic resolution description of interactions provided by all-atom MD simulations, and investigate the oligomerization process of the GNNQQNY for three system sizes: 3-mers, 12-mers and 20-mers. Results for our integrated simulations show a rich variety of structural arrangements for aggregates of all sizes. Elongated fibril-like structures can form transiently in the 20-mer case, but they are not stable and easily interconvert in more globular and disordered forms. Our extensive characterization of the intermediate structures and their physico-chemical determinants points to a high degree of polymorphism for the GNNQQNY sequence that can be reflected at the macroscopic scale. Detailed mechanisms and structures that underlie amyloid aggregation are also provided.

Structures of the Alzheimer's Wild-Type Aβ1-40 Dimer from Atomistic Simulations.

We have studied the dimer of amyloid beta peptide Aβ of 40 residues by means of all-atom replica exchange molecular dynamics. The Aβ-dimers have been found to be the smallest toxic species in Alzheimers disease, but their inherent flexibilities have precluded structural characterization by experimental methods. Though the 24-μs-scale simulation reveals a mean secondary structure of 18% β-strand and 10% α helix, we find transient configurations with an unstructured N-terminus and multiple β-hairpins spanning residues 17-21 and 30-36, but the antiparallel and perpendicular peptide orientations are preferred over the parallel organization. Short-lived conformational states also consist of all α topologies, and one compact peptide with β-sheet structure stabilized by a rather extended peptide with α-helical content. Overall, this first all-atom study provides insights into the equilibrium structure of the Aβ1-40 dimer in aqueous solution, opening a new avenue for a comprehensive understanding of the impact of pathogenic and protective mutations in early-stage Alzheimers disease on a molecular level.

Kinetics of Amyloid Aggregation: A Study of the GNNQQNY Prion Sequence

The small amyloid-forming GNNQQNY fragment of the prion sequence has been the subject of extensive experimental and numerical studies over the last few years. Using unbiased molecular dynamics with the OPEP coarse-grained potential, we focus here on the onset of aggregation in a 20-mer system. With a total of 16.9 of simulations at 280 K and 300 K, we show that the GNNQQNY aggregation follows the classical nucleation theory (CNT) in that the number of monomers in the aggregate is a very reliable descriptor of aggregation. We find that the critical nucleus size in this finite-size system is between 4 and 5 monomers at 280 K and 5 and 6 at 300 K, in overall agreement with experiment. The kinetics of growth cannot be fully accounted for by the CNT, however. For example, we observe considerable rearrangements after the nucleus is formed, as the system attempts to optimize its organization. We also clearly identify two large families of structures that are selected at the onset of aggregation demonstrating the presence of well-defined polymorphism, a signature of amyloid growth, already in the 20-mer aggregate.

Folding Atomistic Proteins in Explicit Solvent Using Simulated Tempering.

Following a previous report on a coarse-grained protein model in implicit solvent, we applied simulated tempering (ST) with on-the-fly Helmholtz free energy (weight factors) determination to the folding or aggregation of seven proteins with the CHARMM, OPLS, and AMBER protein, and the SPC and TIP3P water force fields. For efficiency and reliability, we also performed replica exchange molecular dynamics (REMD) simulations on the alanine di- and deca-peptide, and the dimer of the Aβ16-22 Alzheimers fragment, and used experimental data and previous simulation results on the chignolin, beta3s, Trp-cage, and WW domain peptides of 10-37 amino acids. The sampling with ST is found to be more efficient than with REMD for a much lower CPU cost. Starting from unfolded or extended conformations, the WW domain and the Trp-cage peptide fold to their NMR structures with a backbone RMSD of 2.0 and 1 Å. Remarkably, the ST simulation explores transient non-native topologies for Trp-cage that have been rarely discussed by other simulations. Our ST simulations also show that the CHARMM22* force field has limitations in describing accurately the beta3s peptide. Taken together, these results open the door to the study of the configurations of single proteins, protein aggregates, and any molecular systems at atomic details in explicit solvent using a single normal CPU. They also demonstrate that our ST scheme can be used with any force field ranging from quantum mechanics to coarse-grain and atomistic.

Impact of the A2V Mutation on the Heterozygous and Homozygous Aβ1–40 Dimer Structures from Atomistic Simulations

The A2V mutation was reported to protect from Alzheimers disease in its heterozygous form and cause an early Alzheimers disease type dementia in its homozygous form. Experiments showed that the aggregation rate follows the order A2V > WT (wild-type) > A2V-WT. To understand the impact of this mutation, we carried out replica exchange molecular dynamics simulations of Aβ1-40 WT-A2V and A2V-A2V dimers and compared to the WT dimer. Our atomistic simulations reveal that the mean secondary structure remains constant, but there are substantial differences in the intramolecular and intermolecular conformations upon single and double A2V mutation. Upon single mutation, the intrinsic disorder is reduced, the intermolecular potential energies are reduced, the population of intramolecular three-stranded β-sheets is increased, and the number of all α dimer topologies is decreased. Taken together, these results offer an explanation for the reduced aggregation rate of the Aβ1-40 A2V-WT peptides and the protective effect of A2V in heterozygotes.

Why Is Research on Amyloid-β Failing to Give New Drugs for Alzheimer’s Disease?

The two hallmarks of Alzheimers disease (AD) are the presence of neurofibrillary tangles (NFT) made of aggregates of the hyperphosphorylated tau protein and of amyloid plaques composed of amyloid-β (Aβ) peptides, primarily Aβ1-40 and Aβ1-42. Targeting the production, aggregation, and toxicity of Aβ with small molecule drugs or antibodies is an active area of AD research due to the general acceptance of the amyloid cascade hypothesis, but thus far all drugs targeting Aβ have failed. From a review of the recent literature and our own experience based on in vitro, in silico, and in vivo studies, we present some reasons to explain this repetitive failure.