1
|
Wu M, Wille H, Stepanova M. Essential collective dynamics analysis reveals nonlocal interactions of alpha-synuclein 38-95 monomers with fibrillar seeds. J Chem Phys 2022; 157:235101. [PMID: 36550041 DOI: 10.1063/5.0116032] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
Abstract
Alpha-synuclein (α-syn) is a cytoplasmic protein commonly found in the nervous system. In solution, α-syn adopts disordered unfolded conformations, although it can also form α-helices upon binding to lipid membranes. Under conditions that are not yet fully understood, α-syn can misfold and aggregate, giving rise to β-sheet rich amyloid fibrils, which then tend to accumulate in degenerating neurons. This leads to Parkinson's disease (PD) and several other conditions collectively termed synucleinopathies. Development of disease-modifying treatments requires detailed understanding of structure and dynamics of α-syn's misfolded aggregates. We have employed 1000 ns long, all-atom molecular dynamics simulations to investigate the interaction of monomeric α-syn38-95 fragments, which contain the most important amyloidogenic regions, with preformed fibrillar seeds composed of staggered, β-sheet rich α-syn chains of matching length. The simulations indicate that α-syn38-95 monomers tend to form aggregates with the fibrillar seeds, although we have not observed alignment of the monomeric chains with β-strands of the fibril. To analyze the stability of these aggregates, we have employed the essential collective dynamics method, which allows making accurate assessment of dynamical coupling across individual atoms in macromolecules and supramolecular complexes. The analysis revealed extensive dynamical coupling across initially monomeric α-syn chains and the fibrillar seeds including distal regions thereof that did not contact the monomer directly. We have discussed structural origins of these long-range interactions, their impacts for the stability of α-syn aggregates, and potential implications for the development of anti-PD treatments.
Collapse
Affiliation(s)
- Min Wu
- Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Alberta T6G 2E1, Canada
| | - Holger Wille
- Department of Biochemistry, University of Alberta, Edmonton, Alberta T6G 2H7, Canada
| | - Maria Stepanova
- Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Alberta T6G 2E1, Canada
| |
Collapse
|
2
|
Wu M, Dorosh L, Schmitt-Ulms G, Wille H, Stepanova M. Aggregation of Aβ40/42 chains in the presence of cyclic neuropeptides investigated by molecular dynamics simulations. PLoS Comput Biol 2021; 17:e1008771. [PMID: 33711010 PMCID: PMC7990313 DOI: 10.1371/journal.pcbi.1008771] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2020] [Revised: 03/24/2021] [Accepted: 02/04/2021] [Indexed: 11/18/2022] Open
Abstract
Alzheimer’s disease is associated with the formation of toxic aggregates of amyloid beta (Aβ) peptides. Despite tremendous efforts, our understanding of the molecular mechanisms of aggregation, as well as cofactors that might influence it, remains incomplete. The small cyclic neuropeptide somatostatin-14 (SST14) was recently found to be the most selectively enriched protein in human frontal lobe extracts that binds Aβ42 aggregates. Furthermore, SST14’s presence was also found to promote the formation of toxic Aβ42 oligomers in vitro. In order to elucidate how SST14 influences the onset of Aβ oligomerization, we performed all-atom molecular dynamics simulations of model mixtures of Aβ42 or Aβ40 peptides with SST14 molecules and analyzed the structure and dynamics of early-stage aggregates. For comparison we also analyzed the aggregation of Aβ42 in the presence of arginine vasopressin (AVP), a different cyclic neuropeptide. We observed the formation of self-assembled aggregates containing the Aβ chains and small cyclic peptides in all mixtures of Aβ42–SST14, Aβ42–AVP, and Aβ40–SST14. The Aβ42–SST14 mixtures were found to develop compact, dynamically stable, but small aggregates with the highest exposure of hydrophobic residues to the solvent. Differences in the morphology and dynamics of aggregates that comprise SST14 or AVP appear to reflect distinct (1) regions of the Aβ chains they interact with; (2) propensities to engage in hydrogen bonds with Aβ peptides; and (3) solvent exposures of hydrophilic and hydrophobic groups. The presence of SST14 was found to impede aggregation in the Aβ42–SST14 system despite a high hydrophobicity, producing a stronger “sticky surface” effect in the aggregates at the onset of Aβ42–SST14 oligomerization. Improper folding of proteins causes disorders known as protein misfolding diseases. Under normal conditions most proteins adopt particular folds, which allow them functioning properly. However, for reasons that are not yet fully understood, proteins may misfold and aggregate, forming deposits known as amyloid fibrils, which accumulate in the brain or other tissues. This process affects functioning of the nervous system, gradually causing loss of cognitive abilities. Alzheimer’s disease is one of the most common diseases from this group. A better understanding of the aggregation of peptides implicated in Alzheimer’s disease, known as amyloid beta (Aβ) peptides, may facilitate the development of treatments that ameliorate or prevent the disease. We use detailed molecular dynamics simulations to investigate the influence of somatostatin-14 (SST14), a cyclic neuropeptide that might be involved in the amyloidogenic aggregation of Aβ, on molecular processes occurring at the onset of Aβ aggregation. Results of these simulations explain how the presence of SST14 might alter pathways of aggregation of Aβ, shedding light upon the possible role of extrinsic factors in the aggregation at a molecular level.
Collapse
Affiliation(s)
- Min Wu
- Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Canada
| | - Lyudmyla Dorosh
- Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Canada
| | - Gerold Schmitt-Ulms
- Laboratory Medicine and Pathobiology, University of Toronto, Toronto, Canada
| | - Holger Wille
- Department of Biochemistry, University of Alberta, Edmonton, Canada
- Centre for Prions and Protein Folding Diseases, Edmonton, Canada
- Neuroscience and Mental Health Institute, University of Alberta, Edmonton, Canada
| | - Maria Stepanova
- Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Canada
- * E-mail:
| |
Collapse
|
3
|
Wille H, Dorosh L, Amidian S, Schmitt-Ulms G, Stepanova M. Combining molecular dynamics simulations and experimental analyses in protein misfolding. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2020; 118:33-110. [PMID: 31928730 DOI: 10.1016/bs.apcsb.2019.10.001] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The fold of a protein determines its function and its misfolding can result in loss-of-function defects. In addition, for certain proteins their misfolding can lead to gain-of-function toxicities resulting in protein misfolding diseases such as Alzheimer's, Parkinson's, or the prion diseases. In all of these diseases one or more proteins misfold and aggregate into disease-specific assemblies, often in the form of fibrillar amyloid deposits. Most, if not all, protein misfolding diseases share a fundamental molecular mechanism that governs the misfolding and subsequent aggregation. A wide variety of experimental methods have contributed to our knowledge about misfolded protein aggregates, some of which are briefly described in this review. The misfolding mechanism itself is difficult to investigate, as the necessary timescale and resolution of the misfolding events often lie outside of the observable parameter space. Molecular dynamics simulations fill this gap by virtue of their intrinsic, molecular perspective and the step-by-step iterative process that forms the basis of the simulations. This review focuses on molecular dynamics simulations and how they combine with experimental analyses to provide detailed insights into protein misfolding and the ensuing diseases.
Collapse
Affiliation(s)
- Holger Wille
- Department of Biochemistry, University of Alberta, Edmonton, Canada; Centre for Prions and Protein Folding Diseases, University of Alberta, Edmonton, Canada; Neuroscience and Mental Health Institute, University of Alberta, Edmonton, Canada
| | - Lyudmyla Dorosh
- Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Canada
| | - Sara Amidian
- Department of Biochemistry, University of Alberta, Edmonton, Canada; Centre for Prions and Protein Folding Diseases, University of Alberta, Edmonton, Canada
| | - Gerold Schmitt-Ulms
- Laboratory Medicine and Pathobiology, University of Toronto, Toronto, Canada
| | - Maria Stepanova
- Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Canada
| |
Collapse
|
4
|
Abstract
AbstractThe dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.
Collapse
|
5
|
Jani V, Sonavane U, Joshi R. Detecting early stage structural changes in wild type, pathogenic and non-pathogenic prion variants using Markov state model. RSC Adv 2019; 9:14567-14579. [PMID: 35519320 PMCID: PMC9064127 DOI: 10.1039/c9ra01507h] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2019] [Accepted: 04/28/2019] [Indexed: 12/19/2022] Open
Abstract
The conversion of prion protein from normal to scrapie followed by the aggregation and deposition of this scrapie form leads to various neurodegenerative diseases. A few studies carried out by researchers suggest that E219K prion mutant (glutamate to lysine mutation at residue position 219) is more stable than wild type protein. However a similar point mutation E200K (glutamate to lysine mutation at residue position 200) is pathogenic. In this study we have carried out detailed atomistic simulation of the wild type, pathogenic mutant E200K and E219K mutant which provides more stability. The aim of the study was to detect the early structural changes present in all the three variants which might be responsible for the stability or for their conversion from PrPC to PrPSc. MSM based analyses have been carried out to find out the differences between WT, E200K and E219K systems. Markov state model (MSM) analysis was able to predict the intermediate states which helped to understand the effect of same mutation at two different locations. The MSM analysis was able to show that the extra stability of E219K mutant may be a result of the increase in number of native contacts, strong salt bridges and less random motions. While pathogenicity of E200K mutant can be attributed to loss of some crucial salt-bridge interactions, increased random motions between helix 2 and helix 3. Markov state model to find out the differences between WT, E200K and E219K systems.![]()
Collapse
Affiliation(s)
- Vinod Jani
- High Performance Computing-Medical & Bioinformatics Applications Group
- Centre for Development of Advanced Computing (C-DAC)
- Savitribai Phule Pune University Campus
- Pune 411007
- India
| | - Uddhavesh Sonavane
- High Performance Computing-Medical & Bioinformatics Applications Group
- Centre for Development of Advanced Computing (C-DAC)
- Savitribai Phule Pune University Campus
- Pune 411007
- India
| | - Rajendra Joshi
- High Performance Computing-Medical & Bioinformatics Applications Group
- Centre for Development of Advanced Computing (C-DAC)
- Savitribai Phule Pune University Campus
- Pune 411007
- India
| |
Collapse
|
6
|
Mercer RCC, Daude N, Dorosh L, Fu ZL, Mays CE, Gapeshina H, Wohlgemuth SL, Acevedo-Morantes CY, Yang J, Cashman NR, Coulthart MB, Pearson DM, Joseph JT, Wille H, Safar JG, Jansen GH, Stepanova M, Sykes BD, Westaway D. A novel Gerstmann-Sträussler-Scheinker disease mutation defines a precursor for amyloidogenic 8 kDa PrP fragments and reveals N-terminal structural changes shared by other GSS alleles. PLoS Pathog 2018; 14:e1006826. [PMID: 29338055 PMCID: PMC5786331 DOI: 10.1371/journal.ppat.1006826] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2017] [Revised: 01/26/2018] [Accepted: 12/18/2017] [Indexed: 11/29/2022] Open
Abstract
To explore pathogenesis in a young Gerstmann-Sträussler-Scheinker Disease (GSS) patient, the corresponding mutation, an eight-residue duplication in the hydrophobic region (HR), was inserted into the wild type mouse PrP gene. Transgenic (Tg) mouse lines expressing this mutation (Tg.HRdup) developed spontaneous neurologic syndromes and brain extracts hastened disease in low-expressor Tg.HRdup mice, suggesting de novo formation of prions. While Tg.HRdup mice exhibited spongiform change, PrP aggregates and the anticipated GSS hallmark of a proteinase K (PK)-resistant 8 kDa fragment deriving from the center of PrP, the LGGLGGYV insertion also imparted alterations in PrP's unstructured N-terminus, resulting in a 16 kDa species following thermolysin exposure. This species comprises a plausible precursor to the 8 kDa PK-resistant fragment and its detection in adolescent Tg.HRdup mice suggests that an early start to accumulation could account for early disease of the index case. A 16 kDa thermolysin-resistant signature was also found in GSS patients with P102L, A117V, H187R and F198S alleles and has coordinates similar to GSS stop codon mutations. Our data suggest a novel shared pathway of GSS pathogenesis that is fundamentally distinct from that producing structural alterations in the C-terminus of PrP, as observed in other prion diseases such as Creutzfeldt-Jakob Disease and scrapie. Prion diseases can be sporadic, infectious or genetic. The central event of all prion diseases is the structural conversion of the cellular prion protein (PrPC) to its disease associated conformer, PrPSc. Gerstmann-Sträussler-Scheinker Disease (GSS) is a genetic prion disease presenting as a multi-systemic neurological syndrome. A novel mutation, an eight amino acid insertion, was discovered in a young GSS patient. We created transgenic mice expressing this mutation and found that they recapitulate key features of the disease; namely PrP deposition in the brain and a low molecular weight proteinase K (PK) resistant internal PrP fragment. While structural investigations did not reveal a gross alteration in the conformation of this mutant PrP, the insertion lying at the boundary of the globular domain causes alterations in the unstructured amino terminal portion of the protein such that it becomes resistant to digestion by the enzyme thermolysin. We demonstrate by kinetic analysis and sequential digestion that this novel thermolysin resistant species is a precursor to the pathognomonic PK resistant fragment. Analysis of samples from other GSS patients revealed this same signature, suggesting a common molecular pathway.
Collapse
Affiliation(s)
- Robert C. C. Mercer
- Centre for Prions and Protein Folding Diseases, University of Alberta, Edmonton, Alberta, Canada
- Department of Medicine (Neurology), University of Alberta, Edmonton, Alberta, Canada
| | - Nathalie Daude
- Centre for Prions and Protein Folding Diseases, University of Alberta, Edmonton, Alberta, Canada
| | - Lyudmyla Dorosh
- National Research Council of Canada, Edmonton, Alberta, Canada
- Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Alberta, Canada
| | - Ze-Lin Fu
- Centre for Prions and Protein Folding Diseases, University of Alberta, Edmonton, Alberta, Canada
- Department of Biochemistry, University of Alberta, Edmonton, Alberta, Canada
| | - Charles E. Mays
- Centre for Prions and Protein Folding Diseases, University of Alberta, Edmonton, Alberta, Canada
| | - Hristina Gapeshina
- Centre for Prions and Protein Folding Diseases, University of Alberta, Edmonton, Alberta, Canada
| | - Serene L. Wohlgemuth
- Centre for Prions and Protein Folding Diseases, University of Alberta, Edmonton, Alberta, Canada
| | | | - Jing Yang
- Centre for Prions and Protein Folding Diseases, University of Alberta, Edmonton, Alberta, Canada
| | - Neil R. Cashman
- Brain Research Centre, University of British Columbia, Vancouver, British Columbia, Canada
| | - Michael B. Coulthart
- Canadian Creutzfeldt-Jakob Disease Surveillance System, Centre for Foodborne, Environmental and Zoonotic Infectious Diseases, Public Health Agency of Canada, Ottawa, Ontario, Canada
| | - Dawn M. Pearson
- Department of Clinical Neurosciences, University of Calgary, Calgary, Alberta, Canada
| | - Jeffrey T. Joseph
- Hotchkiss Brain Institute and Calgary Laboratory Services, University of Calgary, Calgary, Alberta, Canada
| | - Holger Wille
- Centre for Prions and Protein Folding Diseases, University of Alberta, Edmonton, Alberta, Canada
- Department of Biochemistry, University of Alberta, Edmonton, Alberta, Canada
| | - Jiri G. Safar
- Departments of Pathology and Neurology, School of Medicine Case Western Reserve University, Cleveland, Ohio, United States of America
| | - Gerard H. Jansen
- Canadian Creutzfeldt-Jakob Disease Surveillance System, Centre for Foodborne, Environmental and Zoonotic Infectious Diseases, Public Health Agency of Canada, Ottawa, Ontario, Canada
- Division of Anatomical Pathology, University of Ottawa, Ottawa, Ontario, Canada
| | - Maria Stepanova
- National Research Council of Canada, Edmonton, Alberta, Canada
- Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Alberta, Canada
| | - Brian D. Sykes
- Centre for Prions and Protein Folding Diseases, University of Alberta, Edmonton, Alberta, Canada
- Department of Biochemistry, University of Alberta, Edmonton, Alberta, Canada
| | - David Westaway
- Centre for Prions and Protein Folding Diseases, University of Alberta, Edmonton, Alberta, Canada
- Department of Medicine (Neurology), University of Alberta, Edmonton, Alberta, Canada
- Department of Biochemistry, University of Alberta, Edmonton, Alberta, Canada
- * E-mail:
| |
Collapse
|
7
|
Dorosh L, Stepanova M. Probing oligomerization of amyloid beta peptide in silico. MOLECULAR BIOSYSTEMS 2017; 13:165-182. [PMID: 27844078 DOI: 10.1039/c6mb00441e] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Aggregation of amyloid β (Aβ) peptide is implicated in fatal Alzheimer's disease, for which no cure is available. Understanding the mechanisms responsible for this aggregation is required in order for therapies to be developed. In an effort to better understand the molecular mechanisms involved in spontaneous aggregation of Aβ peptide, extensive molecular dynamics simulations are reported, and the results are analyzed through a combination of structural biology tools and a novel essential collective dynamics method. Several model systems composed of ten or twelve Aβ17-42 chains in water are investigated, and the influence of metal ions is probed. The results suggest that Aβ monomers tend to aggregate into stable globular-like oligomers with 13-23% of β-sheet content. Two stages of oligomer formation have been identified: quick collapse within the first 40 ns of the simulation, characterized by a decrease in inter-chain separation and build-up of β-sheets, and the subsequent slow relaxation of the oligomer structure. The resulting oligomers comprise a stable, coherently moving sub-aggregate of 6-9 strongly inter-correlated chains. Cu2+ and Fe2+ ions have been found to develop coordination bonds with carboxylate groups of E22, D23 and A42, which remain stable during 200 ns simulations. The presence of Fe2+, and particularly Cu2+ ions, in negatively charged cavities has been found to cause significant changes in the structure and dynamics of the oligomers. The results indicate, in particular, that formation of non-fibrillar oligomers might be involved in early template-free aggregation of Aβ17-42 monomers, with charged species such as Cu2+ or Fe2+ ions playing an important role.
Collapse
Affiliation(s)
- L Dorosh
- Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Alberta, Canada. and National Research Council of Canada, Edmonton, Alberta, Canada
| | - M Stepanova
- Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Alberta, Canada. and National Research Council of Canada, Edmonton, Alberta, Canada and Department of Physics, Astronomy, and Materials Science, Missouri State University, Springfield, MO, USA
| |
Collapse
|
8
|
Mane JY, Stepanova M. Understanding the dynamics of monomeric, dimeric, and tetrameric α-synuclein structures in water. FEBS Open Bio 2016; 6:666-86. [PMID: 27398307 PMCID: PMC4932447 DOI: 10.1002/2211-5463.12069] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2016] [Revised: 04/05/2016] [Accepted: 04/06/2016] [Indexed: 01/20/2023] Open
Abstract
Human α‐synuclein (αS) is an intrinsically disordered protein associated with Parkinson's disease. Molecular mechanisms of corruptive misfolding and aggregation of αS resulting in the disease, as well as the structure and other properties of the corresponding oligomers are not entirely understood yet, preventing the development of efficient therapies. In this study, we investigate the folding dynamics of initially unfolded hypothetical αS constructs in water using all‐atom molecular dynamics simulations. We also employ the novel essential collective dynamics method to analyze the results obtained from the simulations. Our comparative analysis of monomeric, dimeric, and tetrameric αS models reveals pronounced differences in their structure and stability, emphasizing the importance of small oligomers, particularly dimers, in the process of misfolding.
Collapse
Affiliation(s)
- Jonathan Y Mane
- Department of Electrical and Computer Engineering University of Alberta Edmonton Canada; National Institute for Nanotechnology National Research Council Canada Edmonton Canada
| | - Maria Stepanova
- Department of Electrical and Computer Engineering University of Alberta Edmonton Canada; National Institute for Nanotechnology National Research Council Canada Edmonton Canada; Department of Physics, Astronomy, and Materials Science Missouri State University Springfield MO USA
| |
Collapse
|
9
|
Stueker O, Ortega VA, Goss GG, Stepanova M. Understanding interactions of functionalized nanoparticles with proteins: a case study on lactate dehydrogenase. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2014; 10:2006-2021. [PMID: 24591162 DOI: 10.1002/smll.201303639] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/23/2013] [Revised: 01/10/2014] [Indexed: 06/03/2023]
Abstract
Nanomaterials in biological solutions are known to interact with proteins and have been documented to affect protein function, such as enzyme activity. Understanding the interactions of nanoparticles with biological components at the molecular level will allow for rational designs of nanomaterials for use in medical technologies. Here we present the first detailed molecular mechanics model of functionalized gold nanoparticle (NP) interacting with an enzyme (L-lactate dehydrogenase (LDH) enzyme). Molecular dynamics (MD) simulations of the response of LDH to the NP binding demonstrate that although atomic motions (dynamics) of the main chain exhibit only a minor response to the binding, the dynamics of side chains are significantly constrained in all four active sites that predict alteration in kinetic properties of the enzyme. It is also demonstrated that the 5 nm gold NPs cause a decrease in the maximal velocity of the enzyme reaction (V(max)) and a trend towards a reduced affinity (increased K(m)) for the β-NAD binding site, while pyruvate enzyme kinetics (K(m) and V(max)) are not significantly altered in the presence of the gold NPs. These results demonstrate that modeling of NP:protein interactions can be used to understand alterations in protein function.
Collapse
Affiliation(s)
- Oliver Stueker
- National Institute for Nanotechnology, National Research Council of Canada, 11421 Saskatchewan Drive, EdmontonAlberta, T6G 2M9, Canada; Department of Electrical and Computer Engineering, University of Alberta, 9107-116 Street, Edmonton, Alberta, T6G 2V4, Canada
| | | | | | | |
Collapse
|
10
|
Evolutionary implications of metal binding features in different species' prion protein: an inorganic point of view. Biomolecules 2014; 4:546-65. [PMID: 24970230 PMCID: PMC4101497 DOI: 10.3390/biom4020546] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2014] [Revised: 04/29/2014] [Accepted: 05/06/2014] [Indexed: 12/21/2022] Open
Abstract
Prion disorders are a group of fatal neurodegenerative conditions of mammals. The key molecular event in the pathogenesis of such diseases is the conformational conversion of prion protein, PrPC, into a misfolded form rich in β-sheet structure, PrPSc, but the detailed mechanistic aspects of prion protein conversion remain enigmatic. There is uncertainty on the precise physiological function of PrPC in healthy individuals. Several evidences support the notion of its role in copper homeostasis. PrPC binds Cu2+ mainly through a domain composed by four to five repeats of eight amino acids. In addition to mammals, PrP homologues have also been identified in birds, reptiles, amphibians and fish. The globular domain of protein is retained in the different species, suggesting that the protein carries out an essential common function. However, the comparison of amino acid sequences indicates that prion protein has evolved differently in each vertebrate class. The primary sequences are strongly conserved in each group, but these exhibit a low similarity with those of mammals. The N-terminal domain of different prions shows tandem amino acid repeats with an increasing amount of histidine residues going from amphibians to mammals. The difference in the sequence affects the number of copper binding sites, the affinity and the coordination environment of metal ions, suggesting that the involvement of prion in metal homeostasis may be a specific characteristic of mammalian prion protein. In this review, we describe the similarities and the differences in the metal binding of different species' prion protein, as revealed by studies carried out on the entire protein and related peptide fragments.
Collapse
|
11
|
Dorosh L, Rajagopalan N, Loewen MC, Stepanova M. Molecular mechanisms in the selective basal activation of pyrabactin receptor 1: Comparative analysis of mutants. FEBS Open Bio 2014; 4:496-509. [PMID: 24944884 PMCID: PMC4060014 DOI: 10.1016/j.fob.2014.05.001] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2014] [Revised: 05/10/2014] [Accepted: 05/13/2014] [Indexed: 01/29/2023] Open
Abstract
Pyrabactin receptors (PYR) play a central role in abscisic acid (ABA) signal transduction; they are ABA receptors that inhibit type 2C protein phosphatases (PP2C). Molecular aspects contributing to increased basal activity of PYR against PP2C are studied by molecular dynamics (MD) simulations. An extensive series of MD simulations of the apo-form of mutagenized PYR1 as a homodimer and in complex with homology to ABA-insensitive 1 (HAB1) phosphatase are reported. In order to investigate the detailed molecular mechanisms mediating PYR1 activity, the MD data was analyzed by essential collective dynamics (ECD), a novel approach that allows the identification, with atomic resolution, of persistent dynamic correlations based on relatively short MD trajectories. Employing the ECD method, the effects of select mutations on the structure and dynamics of the PYR1 complexes were investigated and considered in the context of experimentally determined constitutive activities against HAB1. Approaches to rationally design constitutively active PYR1 constructs to increase PP2C inhibition are discussed.
Collapse
Key Words
- ABA, A8S abscisic acid
- Abscisic acid signaling
- CA, constitutively active
- Constitutively active mutations
- ECD, essential collective dynamics
- Essential collective dynamics analysis
- HAB1, homology to ABA insensitive 1
- MD, molecular dynamics
- Molecular dynamics simulations
- P2M, N-(pyridin-2-ylmethyl) naphthalene-1-sulfonamide or C16H14N2O2S
- PCA, principal component analysis
- PDB, Protein Data Bank
- PP2C, phosphatase type 2C
- PYL, PYR1-like
- PYR1, pyrabactin resistance 1
- PYV, pyrabactin or C16H13BrN2O2S
- Pyrabactin resistance
- RCAR, regulatory component of ABA response
- WT, wild type
Collapse
Affiliation(s)
- Lyudmyla Dorosh
- National Research Council of Canada, Edmonton, Alberta, Canada ; Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Alberta, Canada
| | | | - Michele C Loewen
- National Research Council of Canada, Saskatoon, Saskatchewan, Canada ; Department of Biochemistry, University of Saskatchewan, Saskatoon, Saskatchewan, Canada
| | - Maria Stepanova
- National Research Council of Canada, Edmonton, Alberta, Canada ; Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Alberta, Canada
| |
Collapse
|
12
|
Molecular dynamics simulation of temperature induced unfolding of animal prion protein. J Mol Model 2013; 19:4433-41. [DOI: 10.1007/s00894-013-1955-0] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2013] [Accepted: 07/22/2013] [Indexed: 01/22/2023]
|
13
|
NMR spectroscopy on domain dynamics in biomacromolecules. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2013; 112:58-117. [DOI: 10.1016/j.pbiomolbio.2013.05.001] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/03/2013] [Revised: 05/06/2013] [Accepted: 05/07/2013] [Indexed: 12/22/2022]
|
14
|
Dorosh L, Kharenko OA, Rajagopalan N, Loewen MC, Stepanova M. Molecular mechanisms in the activation of abscisic acid receptor PYR1. PLoS Comput Biol 2013; 9:e1003114. [PMID: 23825939 PMCID: PMC3694813 DOI: 10.1371/journal.pcbi.1003114] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2013] [Accepted: 05/07/2013] [Indexed: 11/26/2022] Open
Abstract
The pyrabactin resistance 1 (PYR1)/PYR1-like (PYL)/regulatory component of abscisic acid (ABA) response (RCAR) proteins comprise a well characterized family of ABA receptors. Recent investigations have revealed two subsets of these receptors that, in the absence of ABA, either form inactive homodimers (PYR1 and PYLs 1–3) or mediate basal inhibition of downstream target type 2C protein phosphatases (PP2Cs; PYLs 4–10) respectively in vitro. Addition of ABA has been shown to release the apo-homodimers yielding ABA-bound monomeric holo-receptors that can interact with PP2Cs; highlighting a competitive-interaction process. Interaction selectivity has been shown to be mediated by subtle structural variations of primary sequence and ligand binding effects. Now, the dynamical contributions of ligand binding on interaction selectivity are investigated through extensive molecular dynamics (MD) simulations of apo and holo-PYR1 in monomeric and dimeric form as well as in complex with a PP2C, homology to ABA insensitive 1 (HAB1). Robust comparative interpretations were enabled by a novel essential collective dynamics approach. In agreement with recent experimental findings, our analysis indicates that ABA-bound PYR1 should efficiently bind to HAB1. However, both ABA-bound and ABA-extracted PYR1-HAB1 constructs have demonstrated notable similarities in their dynamics, suggesting that apo-PYR1 should also be able to make a substantial interaction with PP2Cs, albeit likely with slower complex formation kinetics. Further analysis indicates that both ABA-bound and ABA-free PYR1 in complex with HAB1 exhibit a higher intra-molecular structural stability and stronger inter-molecular dynamic correlations, in comparison with either holo- or apo-PYR1 dimers, supporting a model that includes apo-PYR1 in complex with HAB1. This possibility of a conditional functional apo-PYR1-PP2C complex was validated in vitro. These findings are generally consistent with the competitive-interaction model for PYR1 but highlight dynamical contributions of the PYR1 structure in mediating interaction selectivity suggesting added degrees of complexity in the regulation of the competitive-inhibition. Protein pyrabactin resistance 1 (PYR1) belongs to a group of PYR1-like (PYL) proteins that regulate plant development and responses to conditions of drought and salinity. Recent studies have reported characterization of their molecular structures as well as elucidation of important aspects of their function; highlighting their roles as receptors for the stress responsive phytohormone, abscisic acid (ABA). However details of the molecular mechanisms regulating their receptor signalling remain enigmatic. In this work, we use molecular dynamics simulations complemented by a sophisticated statistical-mechanical analysis to investigate structural and dynamical properties of PYR1 protein and how its interaction with ABA modifies receptor-protein complex formation. Our results provide detailed insight into how the PYR1-mediated inactivation of its downstream phosphatase target is regulated by homodimer formation and yield new hypotheses, supported by in vitro experiments, for further investigation. Ultimately, this knowledge provides insight into how plants respond to stress, with potential applications in the development of crops with improved growth characteristics and higher stress tolerance.
Collapse
Affiliation(s)
- Lyudmyla Dorosh
- National Research Council of Canada, Edmonton, Alberta, Canada
- Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Alberta, Canada
| | | | | | - Michele C. Loewen
- National Research Council of Canada, Saskatoon, Saskatchewan, Canada
- Department of Biochemistry, University of Saskatchewan, Saskatoon, Saskatchewan, Canada
| | - Maria Stepanova
- National Research Council of Canada, Edmonton, Alberta, Canada
- Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Alberta, Canada
- * E-mail:
| |
Collapse
|
15
|
Chen J, Thirumalai D. Helices 2 and 3 are the initiation sites in the PrP(C) → PrP(SC) transition. Biochemistry 2012; 52:310-9. [PMID: 23256626 DOI: 10.1021/bi3005472] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
It is established that prion protein is the sole causative agent in a number of diseases in humans and animals. However, the nature of conformational changes that the normal cellular form, PrP(C), undergoes in its conversion to a self-replicating state is still not fully understood. The ordered C-terminus of PrP(C) proteins has three helices (H1-H3). Here, we use statistical coupling analysis (SCA) to infer covariations at various locations using a family of evolutionarily related sequences and the response of mouse and human PrP(C)s to mechanical force to decipher the initiation sites for the transition from PrP(C) to an aggregation-prone PrP* state. Sequence-based SCA predicts that the clustered residues in nonmammals are localized in the stable core (near H1) of PrP(C), whereas in mammalian PrP(C), they are localized in frustrated helices H2 and H3 where most of the pathogenic mutations are found. Force-extension curves and free energy profiles as a function of extension of mouse and human PrP(C) in the absence of a disulfide (SS) bond between residues Cys179 and Cys214, generated by applying mechanical force to the ends of the molecule, show a sequence of unfolding events starting first with rupture of H2 and H3. This is followed by disruption of structure in two strands. Helix H1, stabilized by three salt bridges, resists substantial force before unfolding. Force extension profiles and the dynamics of rupture of tertiary contacts also show that even in the presence of an SS bond the instabilities in most of H3 and parts of H2 still determine the propensity to form the PrP* state. In mouse PrP(C) with an SS bond, there are ∼10 residues that retain their order even at high forces. Both SCA and single-molecule force simulations show that in the conversion from PrP(C) to PrP(SC) major conformational changes occur (at least initially) in H2 and H3, which because of their sequence compositions are frustrated in the helical state. Implications of our findings for the structural model for the scrapie form of PrP(C) are discussed.
Collapse
Affiliation(s)
- Jie Chen
- Biophysics Program, Institute for Physical Science and Technology, and ‡Department of Chemistry and Biochemistry, University of Maryland , College Park, Maryland 20742, USA
| | | |
Collapse
|
16
|
Gendoo DMA, Harrison PM. The landscape of the prion protein's structural response to mutation revealed by principal component analysis of multiple NMR ensembles. PLoS Comput Biol 2012; 8:e1002646. [PMID: 22912570 PMCID: PMC3415401 DOI: 10.1371/journal.pcbi.1002646] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2011] [Accepted: 07/04/2012] [Indexed: 11/18/2022] Open
Abstract
Prion Proteins (PrP) are among a small number of proteins for which large numbers of NMR ensembles have been resolved for sequence mutants and diverse species. Here, we perform a comprehensive principle components analysis (PCA) on the tertiary structures of PrP globular proteins to discern PrP subdomains that exhibit conformational change in response to point mutations and clade-specific evolutionary sequence mutation trends. This is to our knowledge the first such large-scale analysis of multiple NMR ensembles of protein structures, and the first study of its kind for PrPs. We conducted PCA on human (n = 11), mouse (n = 14), and wildtype (n = 21) sets of PrP globular structures, from which we identified five conformationally variable subdomains within PrP. PCA shows that different non-local patterns and rankings of variable subdomains arise for different pathogenic mutants. These subdomains may thus be key areas for initiating PrP conversion during disease. Furthermore, we have observed the conformational clustering of divergent TSE-non-susceptible species pairs; these non-phylogenetic clusterings indicate structural solutions towards TSE resistance that do not necessarily coincide with evolutionary divergence. We discuss the novelty of our approach and the importance of PrP subdomains in structural conversion during disease. Prion Proteins (PrP) cause a variety of incurable TSE diseases, and are among a small number of proteins for which large numbers of NMR ensembles have been resolved for sequence mutants and diverse species. Here, we perform a comprehensive PCA study to assess conformational variation and discern the landscape of the PrP structural response to sequence mutation. This is to our knowledge the first large-scale analysis of multiple NMR ensembles for a specific protein, and the first study to perform a multivariate PCA on the native globular structures of PrP. We conducted exhaustive PCA on three PrP subsets: human and mouse subsets that include structures of sequence mutants, and the set of wild-type PrP (16 PrP species). PCA shows that different non-local patterns of variable subdomains arise for different pathogenic mutants. These subdomains may thus be key areas for initiating PrP conversion during disease. Furthermore, we observed that some evolutionarily divergent species that are non-susceptible to TSEs have surprising structural similarities in their PrPs. We discuss the novelty of our approach with respect to prions, and the advantage of this analysis as a fast, reliable starting point to identify interesting domains that may warrant further experimental and computational analysis.
Collapse
Affiliation(s)
- Deena M. A. Gendoo
- Department of Biology, McGill University, Montreal, Quebec, Canada
- McGill Center for Bioinformatics, McGill University, Montreal, Quebec, Canada
| | - Paul M. Harrison
- Department of Biology, McGill University, Montreal, Quebec, Canada
- McGill Center for Bioinformatics, McGill University, Montreal, Quebec, Canada
- * E-mail:
| |
Collapse
|
17
|
Issack BB, Berjanskii M, Wishart DS, Stepanova M. Exploring the essential collective dynamics of interacting proteins: application to prion protein dimers. Proteins 2012; 80:1847-65. [PMID: 22488640 DOI: 10.1002/prot.24082] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2011] [Revised: 03/07/2012] [Accepted: 03/18/2012] [Indexed: 11/11/2022]
Abstract
Essential collective dynamics is a promising and robust approach for analysing the slow motions of macromolecules from short molecular dynamics trajectories. In this study, an extension of the method to treat a collection of interacting protein molecules is presented. The extension is applied to investigate the effects of dimerization on the collective dynamics of ovine prion protein molecules in two different arrangements. Examination of the structural plasticity shows that aggregation has a restricting effect on the local mobility of the prion protein molecules in the interfacial regions. Domain motions of the two dimeric ovine prion protein conformations are distinctly different and can be related to interatomic correlations at the interface. Correlated motions are among the slow collective modes extensively analysed by considering both main-chain and side-chain atoms. Correlation maps reveal the existence of a vast network of dynamically correlated side groups, which extends beyond individual subunits via interfacial interconnections. The network is formed by a core of hydrophobic side chains surrounded by hydrophilic groups at the periphery. The relevance of these findings are discussed in the context of mutations associated with prion diseases. The binding free energy of the dimeric conformations is evaluated to probe their thermodynamic stability. The descriptions afforded by the analysis of the essential collective dynamics of the prion dimers are consistent with their binding free energies. The agreement validates the extension of the methodology and provides a means of interpreting the collective dynamics in terms of the thermodynamic stability of ovine prion proteins.
Collapse
Affiliation(s)
- Bilkiss B Issack
- National Institute for Nanotechnology, National research Council, Edmonton, AB, Canada
| | | | | | | |
Collapse
|
18
|
Barakat K, Issack BB, Stepanova M, Tuszynski J. Effects of temperature on the p53-DNA binding interactions and their dynamical behavior: comparing the wild type to the R248Q mutant. PLoS One 2011; 6:e27651. [PMID: 22110706 PMCID: PMC3218007 DOI: 10.1371/journal.pone.0027651] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2011] [Accepted: 10/21/2011] [Indexed: 12/20/2022] Open
Abstract
Background The protein p53 plays an active role in the regulation of cell cycle. In about half of human cancers, the protein is inactivated by mutations located primarily in its DNA-binding domain. Interestingly, a number of these mutations possess temperature-induced DNA-binding characteristics. A striking example is the mutation of Arg248 into glutamine or tryptophan. These mutants are defective for binding to DNA at 310 K although they have been shown to bind specifically to several p53 response elements at sub-physiological temperatures (298–306 K). Methodology/Principal Findings This important experimental finding motivated us to examine the effects of temperature on the structure and configuration of R248Q mutant and compare it to the wild type protein. Our aim is to determine how and where structural changes of mutant variants take place due to temperature changes. To answer these questions, we compared the mutant to the wild-type proteins from two different aspects. First, we investigated the systems at the atomistic level through their DNA-binding affinity, hydrogen bond networks and spatial distribution of water molecules. Next, we assessed changes in their long-lived conformational motions at the coarse-grained level through the collective dynamics of their side-chain and backbone atoms separately. Conclusions The experimentally observed effect of temperature on the DNA-binding properties of p53 is reproduced. Analysis of atomistic and coarse-grained data reveal that changes in binding are determined by a few key residues and provide a rationale for the mutant-loss of binding at physiological temperatures. The findings can potentially enable a rescue strategy for the mutant structure.
Collapse
Affiliation(s)
- Khaled Barakat
- Department of Physics, University of Alberta, Edmonton, Alberta, Canada
- Department of Engineering Mathematics and Physics, Fayoum University, Fayoum, Egypt
| | - Bilkiss B. Issack
- Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Alberta, Canada
- National Institute for Nanotechnology, National Research Council, Edmonton, Alberta, Canada
| | - Maria Stepanova
- Department of Electrical and Computer Engineering, University of Alberta, Edmonton, Alberta, Canada
- National Institute for Nanotechnology, National Research Council, Edmonton, Alberta, Canada
| | - Jack Tuszynski
- Department of Physics, University of Alberta, Edmonton, Alberta, Canada
- Department of Oncology, University of Alberta, Edmonton, Alberta, Canada
- * E-mail:
| |
Collapse
|