1
|
Torielli L, Serapian SA, Mussolin L, Moroni E, Colombo G. Integrating Protein Interaction Surface Prediction with a Fragment-Based Drug Design: Automatic Design of New Leads with Fragments on Energy Surfaces. J Chem Inf Model 2023; 63:343-353. [PMID: 36574607 PMCID: PMC9832486 DOI: 10.1021/acs.jcim.2c01408] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Protein-protein interactions (PPIs) have emerged in the past years as significant pharmacological targets in the development of new therapeutics due to their key roles in determining pathological pathways. Herein, we present fragments on energy surfaces, a simple and general design strategy that integrates the analysis of the dynamic and energetic signatures of proteins to unveil the substructures involved in PPIs, with docking, selection, and combination of drug-like fragments to generate new PPI inhibitor candidates. Specifically, structural representatives of the target protein are used as inputs for the blind physics-based prediction of potential protein interaction surfaces using the matrix of low coupling energy decomposition method. The predicted interaction surfaces are subdivided into overlapping windows that are used as templates to direct the docking and combination of fragments representative of moieties typically found in active drugs. This protocol is then applied and validated using structurally diverse, important PPI targets as test systems. We demonstrate that our approach facilitates the exploration of the molecular diversity space of potential ligands, with no requirement of prior information on the location and properties of interaction surfaces or on the structures of potential lead compounds. Importantly, the hit molecules that emerge from our ab initio design share high chemical similarity with experimentally tested active PPI inhibitors. We propose that the protocol we describe here represents a valuable means of generating initial leads against difficult targets for further development and refinement.
Collapse
Affiliation(s)
- Luca Torielli
- Department
of Chemistry, University of Pavia, Via Taramelli 12, Pavia27100, Italy
| | - Stefano A. Serapian
- Department
of Chemistry, University of Pavia, Via Taramelli 12, Pavia27100, Italy
| | - Lara Mussolin
- Department
of Woman’s and Child’s Health, Pediatric Hematology,
Oncology and Stem Cell Transplant Center, University of Padua, Via Giustiniani, 3, Padua35128, Italy,Istituto
di Ricerca Pediatrica Città della Speranza, Corso Stati Uniti, 4 F, Padova35127, Italy
| | | | - Giorgio Colombo
- Department
of Chemistry, University of Pavia, Via Taramelli 12, Pavia27100, Italy,
| |
Collapse
|
2
|
Marchetti F, Capelli R, Rizzato F, Laio A, Colombo G. The Subtle Trade-Off between Evolutionary and Energetic Constraints in Protein-Protein Interactions. J Phys Chem Lett 2019; 10:1489-1497. [PMID: 30855965 DOI: 10.1021/acs.jpclett.9b00191] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
Life machinery, although overwhelmingly complex, is rooted on a rather limited number of molecular processes. One of the most important is protein-protein interaction. Metabolic regulation, protein folding control, and cellular motility are examples of processes based on the fine-tuned interaction of several protein partners. The region on the protein surface devoted to the recognition of a specific partner is essential for the function of the protein and is, therefore, likely to be conserved during evolution. On the other hand, the physical chemistry of amino acids underlies the mechanism of interactions. Both evolutionary and energetic constraints can then be used to build scoring functions capable of recognizing interaction sites. Our working hypothesis is that residues within the interaction interface tend at the same time to be evolutionarily conserved (to preserve their function) and to provide little contribution to the internal stabilization of the structure of their cognate protein, to facilitate conformational adaptation to the partner. Here, we show that for some classes of protein partners (for example, those involved in signal transduction and in enzymes) evolutionary constraints play the key role in defining the interaction surface. In contrast, energetic constraints emerge as more important in protein partners involved in immune response, in inhibitor proteins, and in structural proteins. Our results indicate that a general-purpose scoring function for protein-protein interaction should not be agnostic of the biological function of the partners.
Collapse
Affiliation(s)
- Filippo Marchetti
- Istituto di Chimica del Riconoscimento Molecolare , CNR Via Mario Bianco 9 , 20131 Milano , Italy
- Dipartimento di Chimica , Università degli Studi di Milano , Via Venezian 21 , I-20133 Milano , Italy
| | - Riccardo Capelli
- INM-9/IAS-5 Computational Biomedicine , Forschungszentrum Jülich , Wilhelm-Johnen-Straße , D-54245 Jülich , Germany
| | - Francesca Rizzato
- SISSA, Scuola Internazionale Superiore Studi Avanzati , Via Bonomea 265 , I-34136 Trieste , Italy
| | - Alessandro Laio
- SISSA, Scuola Internazionale Superiore Studi Avanzati , Via Bonomea 265 , I-34136 Trieste , Italy
- ICTP, International Centre for Theoretical Physics , Strada Costiera 11 , I-34100 Trieste , Italy
| | - Giorgio Colombo
- Istituto di Chimica del Riconoscimento Molecolare , CNR Via Mario Bianco 9 , 20131 Milano , Italy
- Dipartimento di Chimica , Università di Pavia , V.le Taramelli 12 , 27100 Pavia , Italy
| |
Collapse
|
3
|
Application of the Fuzzy Oil Drop Model Describes Amyloid as a Ribbonlike Micelle. ENTROPY 2017. [DOI: 10.3390/e19040167] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
|
4
|
Verma S, Singh A, Mishra A. The effect of fulvic acid on pre‐ and postaggregation state of Aβ17–42: Molecular dynamics simulation studies. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2013; 1834:24-33. [DOI: 10.1016/j.bbapap.2012.08.016] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/27/2012] [Revised: 07/25/2012] [Accepted: 08/20/2012] [Indexed: 11/25/2022]
|
5
|
Bieler NS, Knowles TPJ, Frenkel D, Vácha R. Connecting macroscopic observables and microscopic assembly events in amyloid formation using coarse grained simulations. PLoS Comput Biol 2012; 8:e1002692. [PMID: 23071427 PMCID: PMC3469425 DOI: 10.1371/journal.pcbi.1002692] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2011] [Accepted: 07/26/2012] [Indexed: 01/06/2023] Open
Abstract
The pre-fibrillar stages of amyloid formation have been implicated in cellular toxicity, but have proved to be challenging to study directly in experiments and simulations. Rational strategies to suppress the formation of toxic amyloid oligomers require a better understanding of the mechanisms by which they are generated. We report Dynamical Monte Carlo simulations that allow us to study the early stages of amyloid formation. We use a generic, coarse-grained model of an amyloidogenic peptide that has two internal states: the first one representing the soluble random coil structure and the second one the -sheet conformation. We find that this system exhibits a propensity towards fibrillar self-assembly following the formation of a critical nucleus. Our calculations establish connections between the early nucleation events and the kinetic information available in the later stages of the aggregation process that are commonly probed in experiments. We analyze the kinetic behaviour in our simulations within the framework of the theory of classical nucleated polymerisation, and are able to connect the structural events at the early stages in amyloid growth with the resulting macroscopic observables such as the effective nucleus size. Furthermore, the free-energy landscapes that emerge from these simulations allow us to identify pertinent properties of the monomeric state that could be targeted to suppress oligomer formation. A number of normally soluble proteins can form amyloid structures in a process associated with neurodegenerative diseases such as Alzheimer's and Parkinson's diseases. Mature amyloid structures consist of large fibrils containing thousands of individual proteins aggregated into linear nanostructures; there is increasing evidence, however, that the toxic species responsible for neurodegeneration are not the mature fibrils themselves but rather lower molecular weight precursors commonly known as amyloid oligomers. Unfortunately, these early oligomers are commonly thermodynamically unstable and of nanometer scale dimensions, factors which make them highly challenging to probe in detail in experiments. We have used computer simulations of a model inspired by Alzheimer's Abeta peptide to investigate the early stages of protein aggregation. The results that we obtain were shown to fit Oosawa's polymerization theory, a finding which allows us to provide a connection between the microscopic molecular parameters and macroscopic growth. One crucial parameter is size of the nucleus, i.e. the basic oligomer existing at origin of the formation of each fiber. We have revealed a path for the formation of this nucleus and validate its size by several methods. Our results provide fundamental information for influencing the early stages of amyloid formation in a rational manner.
Collapse
Affiliation(s)
- Noah S. Bieler
- Laboratory for Physical Chemistry, ETH Zürich, Zurich, Switzerland
| | | | - Daan Frenkel
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Robert Vácha
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
- National Centre for Biomolecular Research, Faculty of Science and CEITEC - Central European Institute of Technology, Masaryk University, Brno-Bohunice, Czech Republic
- * E-mail:
| |
Collapse
|
6
|
Schreck JS, Yuan JM. A statistical mechanical approach to protein aggregation. J Chem Phys 2011; 135:235102. [DOI: 10.1063/1.3666837] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
|
7
|
Schmit JD, Ghosh K, Dill K. What drives amyloid molecules to assemble into oligomers and fibrils? Biophys J 2011; 100:450-8. [PMID: 21244841 PMCID: PMC3021675 DOI: 10.1016/j.bpj.2010.11.041] [Citation(s) in RCA: 99] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2010] [Revised: 11/06/2010] [Accepted: 11/09/2010] [Indexed: 11/22/2022] Open
Abstract
We develop a theory for three states of equilibrium of amyloid peptides: the monomer, oligomer, and fibril. We assume that the oligomeric state is a disordered micellelike collection of a few peptide chains held together loosely by hydrophobic interactions into a spherical hydrophobic core. We assume that fibrillar amyloid chains are aligned and further stabilized by steric zipper interactions-hydrogen bonding, steric packing, and specific hydrophobic side-chain contacts. The model makes a broad set of predictions that are consistent with experimental results: 1), Similar to surfactant micellization, amyloid oligomerization should increase with peptide concentration in solution. 2), The onset of fibrillization limits the concentration of oligomers in the solution. 3), The extent of Aβ fibrillization increases with peptide concentration. 4), The predicted average fibril length versus monomer concentration agrees with data on α-synuclein. 5), Full fibril length distributions agree with data on α-synuclein. 6), Denaturants should melt out fibrils. And finally, 7), added salt should stabilize fibrils by reducing repulsions between amyloid peptide chains. It is of interest that small changes in solvent conditions can tip the equilibrium balance between oligomer and fibril and cause large changes in rates through effects on the transition-state barrier. This model may provide useful insights into the physical processes underlying amyloid diseases.
Collapse
Affiliation(s)
- Jeremy D Schmit
- Department of Pharmaceutical Chemistry, University of California, San Francisco, USA.
| | | | | |
Collapse
|
8
|
Lee CF. Self-assembly of protein amyloids: a competition between amorphous and ordered aggregation. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 80:031922. [PMID: 19905161 DOI: 10.1103/physreve.80.031922] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/09/2008] [Revised: 08/12/2009] [Indexed: 05/28/2023]
Abstract
Protein aggregation in the form of amyloid fibrils has important biological and technological implications. Although the self-assembly process is highly efficient, aggregates not in the fibrillar form would also occur and it is important to include these disordered species when discussing the thermodynamic equilibrium behavior of the system. Here, we initiate such a task by considering a mixture of monomeric proteins and the corresponding aggregates in the disordered form (micelles) and in the fibrillar form (amyloid fibrils). Starting with a model on the respective binding free energies for these species, we calculate their concentrations at thermal equilibrium. We then discuss how the incorporation of the disordered structure furthers our understanding on the various amyloid promoting factors observed empirically, and on the kinetics of fibrilization.
Collapse
Affiliation(s)
- Chiu Fan Lee
- Physics Department, Clarendon Laboratory, Oxford University, Parks Road, Oxford OX1 3PU, United Kingdom.
| |
Collapse
|
9
|
Nayak A, Sorci M, Krueger S, Belfort G. A universal pathway for amyloid nucleus and precursor formation for insulin. Proteins 2009; 74:556-65. [PMID: 18655073 DOI: 10.1002/prot.22169] [Citation(s) in RCA: 58] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
To help identify the etiological agents for amyloid-related diseases, attention is focused here on the fibrillar precursors, also called oligomers and protofibrils, and on modeling the reaction kinetics of the formation of the amyloid nucleus. Insulin is a favored model for amyloid formation, not only because amyloidosis can be a problem in diabetes, but also because aggregation and fibrillation causes problems during production, storage, and delivery. Small angle neutron scattering (SANS) is used to measure the temporal formation of insulin oligomers in H(2)O- and D(2)O-based solvents and obtain consistent evidence of the composition of the insulin nucleus that comprised three dimers or six monomers similar to that recently proposed in the literature. A simple molecular structural model that describes the growth of oligomers under a wide range of environmental conditions is proposed. The model first involves lengthening or end-on-end association of dimers to form three-dimer nuclei, and then exhibits broadening or side-on-side association of nuclei. Using different additives to demonstrate their influence on the kinetics of oligomer formation, we showed that, although the time required to form the nucleus was dependent on a specific system, they all followed a universal pathway for nucleus and precursor formation. The methods and analyses presented here provide the first experimental molecular size description of the details of amyloid nucleus formation and subsequent propagation to fibril precursors independent of kinetics.
Collapse
Affiliation(s)
- Arpan Nayak
- Howard P. Isermann Department of Chemical and Biological Engineering, Rensselaer Polytechnic Institute, Troy, New York 12180, USA
| | | | | | | |
Collapse
|
10
|
Pawar A, Favrin G. Finite size effects in simulations of protein aggregation. PLoS One 2008; 3:e2641. [PMID: 18612385 PMCID: PMC2441439 DOI: 10.1371/journal.pone.0002641] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2008] [Accepted: 05/26/2008] [Indexed: 11/30/2022] Open
Abstract
It is becoming increasingly clear that the soluble protofibrillar species that proceed amyloid fibril formation are associated with a range of neurodegenerative disorders such as Alzheimer's and Parkinson diseases. Computer simulations of the processes that lead to the formation of these oligomeric species are starting to make significant contributions to our understanding of the determinants of protein aggregation. We simulate different systems at constant concentration but with a different number of peptides and we study the how the finite number of proteins affects the underlying free energy of the system and therefore the relative stability of the species involved in the process. If not taken into account, this finite size effect can undermine the validity of theoretical predictions regarding the relative stability of the species involved and the rates of conversion from one to the other. We discuss the reasons that give rise to this finite size effect form both a probabilistic and energy fluctuations point of view and also how this problem can be dealt by a finite size scaling analysis.
Collapse
Affiliation(s)
- Amol Pawar
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Giorgio Favrin
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
- * E-mail:
| |
Collapse
|
11
|
Carare RO, Bernardes-Silva M, Newman TA, Page AM, Nicoll JAR, Perry VH, Weller RO. Solutes, but not cells, drain from the brain parenchyma along basement membranes of capillaries and arteries: significance for cerebral amyloid angiopathy and neuroimmunology. Neuropathol Appl Neurobiol 2008; 34:131-44. [PMID: 18208483 DOI: 10.1111/j.1365-2990.2007.00926.x] [Citation(s) in RCA: 435] [Impact Index Per Article: 27.2] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
UNLABELLED Elimination of interstitial fluid and solutes plays a role in homeostasis in the brain, but the pathways are unclear. Previous work suggests that interstitial fluid drains along the walls of arteries. AIMS to define the pathways within the walls of capillaries and arteries for drainage of fluid and solutes out of the brain. METHODS Fluorescent soluble tracers, dextran (3 kDa) and ovalbumin (40 kDa), and particulate fluospheres (0.02 microm and 1.0 microm in diameter) were injected into the corpus striatum of mice. Brains were examined from 5 min to 7 days by immunocytochemistry and confocal microscopy. RESULTS soluble tracers initially spread diffusely through brain parenchyma and then drain out of the brain along basement membranes of capillaries and arteries. Some tracer is takenf up by vascular smooth muscle cells and by perivascular macrophages. No perivascular drainage was observed when dextran was injected into mouse brains following cardiac arrest. Fluospheres expand perivascular spaces between vessel walls and surrounding brain, are ingested by perivascular macrophages but do not appear to leave the brain even following an inflammatory challenge with lipopolysaccharide or kainate. CONCLUSIONS capillary and artery basement membranes act as 'lymphatics of the brain' for drainage of fluid and solutes; such drainage appears to require continued cardiac output as it ceases following cardiac arrest. This drainage pathway does not permit migration of cells from brain parenchyma to the periphery. Amyloid-beta is deposited in basement membrane drainage pathways in cerebral amyloid angiopathy, and may impede elimination of amyloid-beta and interstitial fluid from the brain in Alzheimer's disease. Soluble antigens, but not cells, drain from the brain by perivascular pathways. This atypical pattern of drainage may contribute to partial immune privilege of the brain and play a role in neuroimmunological diseases such as multiple sclerosis.
Collapse
Affiliation(s)
- R O Carare
- Clinical Neurosciences, University of Southampton, UK
| | | | | | | | | | | | | |
Collapse
|
12
|
Abstract
Amyloid fibrils are structures consisting of many proteins with a well-defined conformation. The formation of these fibrils has been the subject of intense research, largely due to their connection to several diseases. We focus here on the computational studies and discuss these from a free-energy point of view. The fibrillogenic properties of many proteins can be predicted and understood by taking the relevant free energies into account in an appropriate way. This is because both the equilibrium and the kinetic properties of the protein system depend on its free-energy landscape. Advanced simulation techniques can be used to understand the relationship between the free-energy landscape of a protein and its three-dimensional structure and propensity to form amyloid fibrils. We give an overview of existing simulation techniques that operate at a molecular level of detail and that are capable of generating relevant free-energy values. The free energies obtained with these methods can be inserted into a statistical-mechanical or kinetic framework to predict mean fibril properties on length scales and time scales that are inaccessible by molecular-scale simulation methods.
Collapse
Affiliation(s)
- Maarten G Wolf
- DelftChemTech, Delft University of Technology, Delft, The Netherlands
| | | | | |
Collapse
|
13
|
Milojevic J, Esposito V, Das R, Melacini G. Understanding the Molecular Basis for the Inhibition of the Alzheimer's Aβ-Peptide Oligomerization by Human Serum Albumin Using Saturation Transfer Difference and Off-Resonance Relaxation NMR Spectroscopy. J Am Chem Soc 2007; 129:4282-90. [PMID: 17367135 DOI: 10.1021/ja067367+] [Citation(s) in RCA: 97] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Human serum albumin (HSA) inhibits the formation of amyloid beta-peptide (Abeta) fibrils in human plasma. However, currently it is not known how HSA affects the formation of the highly toxic soluble diffusible oligomers that occur in the initial stages of Abeta fibrillization. We have therefore investigated by solution NMR the interaction of HSA with the Abeta(12-28) peptide, which has been previously shown to provide a reliable and stable model for the early prefibrillar oligomers as well as to contain key determinants for the recognition by albumin. For this purpose we propose a novel NMR approach based on the comparative analysis of Abeta in its inhibited and filtrated states monitored through both saturation transfer difference and recently developed nonselective off-resonance relaxation experiments. This combined NMR strategy reveals a mechanism for the oligomerization inhibitory function of HSA, according to which HSA targets preferentially the soluble oligomers of Abeta(12-28) rather than its monomeric state. Specifically, HSA caps the exposed hydrophobic patches located at the growing and/or transiently exposed sites of the Abeta oligomers, thereby blocking the addition of further monomers and the growth of the prefibrillar assemblies. The proposed model has implications not only for the pharmacological treatment of Alzheimer's disease specifically but also for the inhibition of oligomerization in amyloid-related diseases in general. In addition, the proposed NMR approach is expected to be useful for the investigation of the mechanism of action of other oligomerization inhibitors as well as of other amyloidogenic systems.
Collapse
Affiliation(s)
- Julijana Milojevic
- Department of Chemistry, McMaster University, 1280 Main Street West, Hamilton, Ontario, L8S 4M1, Canada
| | | | | | | |
Collapse
|
14
|
Hills RD, Brooks CL. Hydrophobic cooperativity as a mechanism for amyloid nucleation. J Mol Biol 2007; 368:894-901. [PMID: 17368485 PMCID: PMC1997311 DOI: 10.1016/j.jmb.2007.02.043] [Citation(s) in RCA: 88] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2006] [Revised: 02/12/2007] [Accepted: 02/14/2007] [Indexed: 11/18/2022]
Abstract
The kinetics of amyloid fibril formation are in most cases explained by classical nucleation theory, yet the mechanisms behind nucleation are not well understood. We show using molecular dynamics simulations that the hydrophobic cooperativity in the self-association of the model amyloidogenic peptide STVIYE is sufficient to allow for nucleation-dependent polymerization with a pentamer critical nucleus. The role of electrostatics was also investigated. Novel considerations of the electrostatic solvation energy using the Born-Onsager equation are put forth to rationalize the aggregation of charged peptides and provide new insight into the energetic differences between parallel and antiparallel beta-sheets. Together these results help explain the influence of molecular charge in the class of fibril-forming hexapeptides recently designed by Serrano and collaborators.
Collapse
|
15
|
Fogolari F, Corazza A, Viglino P, Zuccato P, Pieri L, Faccioli P, Bellotti V, Esposito G. Molecular dynamics simulation suggests possible interaction patterns at early steps of beta2-microglobulin aggregation. Biophys J 2006; 92:1673-81. [PMID: 17158575 PMCID: PMC1796822 DOI: 10.1529/biophysj.106.098483] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Early events in aggregation of proteins are not easily accessible by experiments. In this work, we perform a 5-ns molecular dynamics simulation of an ensemble of 27 copies of beta(2)-microglobulin in explicit solvent. During the simulation, the formation of intermolecular contacts is observed. The simulation highlights the importance of apical residues and, in particular, of those at the N-terminus end of the molecule. The most frequently found pattern of interaction involves a head-to-head contact arrangement of molecules. Hydrophobic contacts appear to be important for the establishment of long-lived (on the simulation timescale) contacts. Although early events on the pathway to aggregation and fibril formation are not directly related to the end-state of the process, which is reached on a much longer timescale, simulation results are consistent with experimental data and in general with a parallel arrangement of intermolecular beta-strand pairs.
Collapse
Affiliation(s)
- Federico Fogolari
- Dipartimento di Scienze e Tecnologie Biomediche, Università di Udine, Udine, Italy.
| | | | | | | | | | | | | | | |
Collapse
|
16
|
Baumketner A, Shea JE. Folding Landscapes of the Alzheimer Amyloid-β(12-28) Peptide. J Mol Biol 2006; 362:567-79. [PMID: 16930617 DOI: 10.1016/j.jmb.2006.07.032] [Citation(s) in RCA: 70] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2006] [Revised: 06/12/2006] [Accepted: 07/17/2006] [Indexed: 11/17/2022]
Abstract
The energy landscape for folding of the 12-28 fragment of the Alzheimer amyloid beta (Abeta) peptide is characterized using replica-exchange molecular dynamics simulations with an all-atom peptide model and explicit solvent. At physiological temperatures, the peptide exists mostly as a collapsed random coil, populating a small fraction (less than 10%) of hairpins with a beta-turn at position V18F19, with another 10% of hairpin-like conformations possessing a bend rather than a turn in the central VFFA positions. A small fraction of the populated states, approximately 14%, adopt polyproline II (PPII) conformations. Folding of the structured hairpin states proceeds through the assembly of two locally stable segments, VFFAE and EDVGS. The interactions stabilizing these locally folded structural motifs are in conflict with those stabilizing the global fold of A12-28, a signature of underlying residual frustration in this peptide. At increased temperature, the population of both beta-strand and PPII conformations diminishes in favor of beta-turn and random-coil states. On the basis of the conformational preferences of Abeta 12-28 monomers, two models for the molecular structure of amyloid fibrils formed by this peptide are proposed.
Collapse
Affiliation(s)
- Andrij Baumketner
- Department of Chemistry and Biochemistry, University of California Santa Barbara, CA 93106, USA
| | | |
Collapse
|
17
|
Chen YR, Huang HB, Chyan CL, Shiao MS, Lin TH, Chen YC. The Effect of Aβ Conformation on the Metal Affinity and Aggregation Mechanism Studied by Circular Dichroism Spectroscopy. ACTA ACUST UNITED AC 2006; 139:733-40. [PMID: 16672274 DOI: 10.1093/jb/mvj083] [Citation(s) in RCA: 49] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Abstract
The conformational change and associated aggregation of beta amyloid (Abeta) with or without metals is the main cause of Alzheimer's disease (AD). In order to further understand the effects of Abeta and its associated metals on the aggregation mechanism, the influence of Abeta conformation on the metal affinity and aggregation was investigated using circular dichroism (CD) spectroscopy. The Abeta conformation is dependent on pH and trifluoroethanol (TFE). The binding of metals to Abeta was found to be dependent on the Abeta conformation. The aggregation induced by Abeta itself or its associated metals is completely diminished for Abeta in 40% TFE. Only in 5% and 25% TFE can Abeta undergo an alpha-helix to beta-sheet aggregation, which involve a three-state mechanism for the metal-free state, and a two-state transition for the metal-bound state, respectively. The aggregation-inducing activity of metals is in the order, Cu2+ > Fe3+ > or = Al3+ > Zn2+.
Collapse
Affiliation(s)
- Y R Chen
- Institute of Medical Science, Tzu Chi University, Hualien, Taiwan
| | | | | | | | | | | |
Collapse
|
18
|
Abstract
Soluble oligomers of Alzheimer's amyloid beta protein (Abeta) may act as effectors of neurotoxicity in early stages of Alzheimer's disease. Detailed information about the structure of Abeta in atomistic level and the dynamics of assembly of monomeric Abeta into oligomeric structures is rather elusive. We have performed replica exchange molecular dynamics (REMD) simulations on the formation of the dimer and trimer of Abeta10-35 peptide. We have observed spontaneous formation of several basic structural units that may act as a template or an intermediate for further aggregation of Alzheimer's Abeta protein. Various conformers, including interlocking structures of experimentally known bend double beta strands, are identified.
Collapse
Affiliation(s)
- Soonmin Jang
- School of Chemistry, Seoul National University, Seoul 151-747, Korea
| | | |
Collapse
|
19
|
López de la Paz M, de Mori GMS, Serrano L, Colombo G. Sequence Dependence of Amyloid Fibril Formation: Insights from Molecular Dynamics Simulations. J Mol Biol 2005; 349:583-96. [PMID: 15882870 DOI: 10.1016/j.jmb.2005.03.081] [Citation(s) in RCA: 72] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2004] [Revised: 03/23/2005] [Accepted: 03/29/2005] [Indexed: 11/19/2022]
Abstract
The clarification of the physico-chemical determinants underlying amyloid deposition is critical for our understanding of misfolding diseases. With this purpose we have performed a systematic all-atom molecular dynamics (MD) study of a series of single point mutants of the de novo designed amyloidogenic peptide STVIIE. Sixteen different 50ns long simulations using explicit solvent have been carried out starting from four different conformations of a polymeric six-stranded beta-sheet. The simulations have provided evidence for the influence of a small number of site-specific hydrophobic interactions on the packing and stabilization of nascent aggregates, as well as the interplay between side-chain interactions and the net charge of the molecule on the strand arrangement of polymeric beta-sheets. This MD analysis has also shed light into the origin of the position dependence on mutation of beta-sheet polymerization that was found experimentally for this model system. Our results suggest that MD can be applied to detect critical positions for beta-sheet aggregation within a given amyloidogenic stretch. Studies similar to the one presented here can guide site-directed mutations or the design of drugs that specifically disrupt the key stabilizing interactions of beta-sheet aggregates.
Collapse
|
20
|
Simona F, Tiana G, Broglia RA, Colombo G. Modeling the α-helix to β-hairpin transition mechanism and the formation of oligomeric aggregates of the fibrillogenic peptide A(12–28): insights from all-atom molecular dynamics simulations. J Mol Graph Model 2004; 23:263-73. [PMID: 15530822 DOI: 10.1016/j.jmgm.2004.07.004] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2004] [Revised: 07/26/2004] [Accepted: 07/26/2004] [Indexed: 11/25/2022]
Abstract
In this paper, all-atom molecular dynamics simulations in explicit solvent are used to investigate the structural and dynamical determinants of the alpha-helical to beta-hairpin conformational transition of the 12-28 fragment from the full length Abeta Alzheimer's peptide. The transition from alpha-helical to beta-structure requires the peptide to populate intermediate beta-bend geometries in which several mainly hydrophobic interactions are partially formed. This is followed by the sudden collapse to ordered beta-hairpin structures and the simultaneous disruption of the hydrophobic side-chain interactions with a consequent increase in solvent exposure. The solvent exposure of hydrophobic side-chains belonging to a sequence of five consecutive residues in the beta-hairpin defines a possible starting point for the onset of the aggregation mechanisms. Several different conformations of model oligomeric (dimeric and tetrameric) aggregates are then investigated. These simulations show that while hydrophobic contacts are important to bring together different monomers with a beta-hairpin like conformation, more specific interactions such as hydrogen-bonding and coulombic interactions, should be considered necessary to provide further stabilization and ordering to the nascent fibrillar aggregates.
Collapse
Affiliation(s)
- Fabio Simona
- Istituto di Chimica del Riconoscimento Molecolare, CNR, Via Mario Bianco 9, 20131 Milano, Italy
| | | | | | | |
Collapse
|
21
|
Cecchini M, Rao F, Seeber M, Caflisch A. Replica exchange molecular dynamics simulations of amyloid peptide aggregation. J Chem Phys 2004; 121:10748-56. [PMID: 15549960 DOI: 10.1063/1.1809588] [Citation(s) in RCA: 177] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The replica exchange molecular dynamics (REMD) approach is applied to four oligomeric peptide systems. At physiologically relevant temperature values REMD samples conformation space and aggregation transitions more efficiently than constant temperature molecular dynamics (CTMD). During the aggregation process the energetic and structural properties are essentially the same in REMD and CTMD. A condensation stage toward disordered aggregates precedes the beta-sheet formation. Two order parameters, borrowed from anisotropic fluid analysis, are used to monitor the aggregation process. The order parameters do not depend on the peptide sequence and length and therefore allow to compare the amyloidogenic propensity of different peptides
Collapse
Affiliation(s)
- M Cecchini
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | | | | | | |
Collapse
|
22
|
Santini S, Mousseau N, Derreumaux P. In silico assembly of Alzheimer's Abeta16-22 peptide into beta-sheets. J Am Chem Soc 2004; 126:11509-16. [PMID: 15366896 DOI: 10.1021/ja047286i] [Citation(s) in RCA: 114] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Recent studies suggest that soluble oligomers of amyloid-forming peptides have toxic effects in cell cultures. In this study, the folding of three Alzheimer's A beta(16-22) peptides have been simulated with the activation-relaxation technique and a generic energy model. Starting from randomly chosen states, the predicted lowest energy structure superposes within 1 A rms deviation from its conformation within the fibrils. This antiparallel structure is found to be in equilibrium with several out-of-register antiparallel beta-sheets and mixed parallel-antiparallel beta-sheets, indicating that full structural order in the fibrils requires larger aggregates. Folding involves the formation of dimers followed by the addition of a monomer and proceeds through a generalized mechanism between disordered and native alignments of beta-strands.
Collapse
Affiliation(s)
- Sébastien Santini
- Contribution from the Information Génomique et Structurale, UPR 2589 CNRS, 31 Chemin Joseph Aiguier, 13402 Marseille Cedex 20, France
| | | | | |
Collapse
|