Molecular dynamics simulations capture the misfolding of the bovine prion protein at acidic pH

The Role of Glial Cells in Neurobiology and Prion Neuropathology

Prion diseases are rare and neurodegenerative diseases that are characterized by the misfolding and infectious spread of the prion protein in the brain, causing progressive and irreversible neuronal loss and associated clinical and behavioral manifestations in humans and animals, ultimately leading to death. The brain has a complex network of neurons and glial cells whose crosstalk is critical for function and homeostasis. Although it is established that prion infection of neurons is necessary for clinical disease to occur, debate remains in the field as to the role played by glial cells, namely astrocytes and microglia, and whether these cells are beneficial to the host or further accelerate disease. Here, we review the current literature assessing the complex morphologies of astrocytes and microglia, and the crosstalk between these two cell types, in the prion-infected brain.

The role of α-sheet structure in amyloidogenesis: characterization and implications

Amyloid diseases are linked to protein misfolding whereby the amyloidogenic protein undergoes a conformational change, aggregates and eventually forms amyloid fibrils. While the amyloid fibrils and plaques are hallmarks of these diseases, they typically form late in the disease process and do not correlate with disease. Instead, there is growing evidence that smaller, soluble toxic oligomers form prior and appear to be early triggers of the molecular pathology underlying these diseases. Nearly 20 years ago, we proposed the α-sheet hypothesis after discovering that the early conformational changes observed during atomistic molecular dynamics simulations involve the formation of a non-standard protein structure, α-sheet. Furthermore, we proposed that toxic oligomers contain α-sheet structure and that preferentially targeting this structure could neutralize the toxicity, prevent further aggregation and serve as the basis for early detection of disease. Here, we present the origin of the α-sheet hypothesis and describe α-sheet structure and the corresponding mechanisms of conversion. We discuss experimental studies demonstrating that both mammalian and bacterial amyloid systems form α-sheet oligomers before converting to conventional β-sheet fibrils. Furthermore, we show that the process can be inhibited with de novo designed α-sheet peptides complementary to the structure in the toxic oligomers.

Epitope-specific anti-PrP antibody toxicity: a comparative in-silico study of human and mouse prion proteins

Despite having therapeutic potential, anti-PrP antibodies caused a major controversy due to their neurotoxic effects. For instance, treating mice with ICSM antibodies delayed prion disease onset, but both were found to be either toxic or innocuous to neurons by researchers following cross-linking PrP^C. In order to elucidate and understand the reasons that led to these contradictory outcomes, we conducted a comprehensive in silico study to assess the antibody-specific toxicity. Since most therapeutic anti-PrP antibodies were generated against human truncated recombinant PrP^91-231 or full-length mouse PrP^23-231, we reasoned that host specificity (human vs murine) of PrP^C might influence the nature of the specific epitopes recognized by these antibodies at the structural level possibly explaining the 'toxicity' discrepancies reported previously. Initially, molecular dynamics simulation and pro-motif analysis of full-length human (hu)PrP and mouse (mo)PrP 3D structure displayed conspicuous structural differences between huPrP and moPrP. We identified 10 huPrP and 6 moPrP linear B-cell epitopes from the prion protein 3D structure where 5 out of 10 huPrP and 3 out of 6 moPrP B-cell epitopes were predicted to be potentially toxic in immunoinformatics approaches. Herein, we demonstrate that some of the predicted potentially 'toxic' epitopes identified by the in silico analysis were similar to the epitopes recognized by the toxic antibodies such as ICSM18 (146-159), POM1 (138-147), D18 (133-157), ICSM35 (91-110), D13 (95-103) and POM3 (95-100). This in silico study reveals the role of host specificity of PrP^C in epitope-specific anti-PrP antibody toxicity.

Combining molecular dynamics simulations and experimental analyses in protein misfolding

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.

Computational insights into pH-dependence of structure and dynamics of pyrazinamidase: A comparison of wild type and mutants

The mycobacterial enzyme pyrazinamidase (PZase) is the target of key tuberculosis drug, pyrazinamide. Mutations in PZase cause drug resistance. Herein, three point mutations, W68G, L85P, and V155G, were investigated through over 8 µs of molecular dynamics simulations coupled with essential dynamics and binding pocket analysis at neutral (pH = 7) and acidic (pH = 4) ambient conditions. The 51-71 flap region exhibited drastic displacement leading to enlargement of binding cavity, especially at the lower pH. Accessibility of solvent to the active site of the mutant enzymes was also reduced. The protonation of key surface residues at low pH results in more contribution of these residues to structural stability and integrity of the enzyme and reduced interactions with solvent molecules, which acts as a cage, keeping the enzyme together. The observed results suggest a pattern of structural alterations due to point mutations in PZase, which is consistent with other experimental and theoretical investigations and, can be harnessed for drug design purposes.

Detecting early stage structural changes in wild type, pathogenic and non-pathogenic prion variants using Markov state model

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.

Structural basis for the complete resistance of the human prion protein mutant G127V to prion disease

Prion diseases are caused by the propagation of misfolded cellular prion proteins (PrPs). A completely prion disease-resistant genotype, V127M129, has been identified in Papua New Guinea and verified in transgenic mice. To disclose the structural basis of the disease-resistant effect of the G127V mutant, we determined and compared the structural and dynamic features of the G127V-mutated human PrP (residues 91-231) and the wild-type PrP in solution. HuPrP(G127V) contains α1, α2 and α3 helices and a stretch-strand (SS) pattern comprising residues Tyr128-Gly131 (SS1) and Val161-Arg164 (SS2), with extending atomic distances between the SS1 and SS2 strands, and a structural rearrangement of the Tyr128 side chain due to steric hindrance of the larger hydrophobic side chain of Val127. The extended α1 helix gets closer to the α2 and α3 helices. NMR dynamics analysis revealed that Tyr128, Gly131 and Tyr163 underwent significant conformational exchanges. Molecular dynamics simulations suggest that HuPrP(G127V) prevents the formation of stable β-sheets and dimers. Unique structural and dynamic features potentially inhibit the conformational conversion of the G127V mutant. This work is beneficial for understanding the molecular mechanisms underlying the complete resistance of the G127V mutant to prion disease and for developing new therapeutics for prion disease.

Probing the early stages of prion protein (PrP) aggregation with atomistic molecular dynamics simulations

Prions are self-replicating infectious proteinaceous agents whose conformations are capable of forming amyloid-like aggregate fibrils. Here we present molecular dynamics simulations aimed at investigating the aggregation process of the β-rich H2H3 domain of the ovine prion protein (H2H3-OvPrPSc), known to be the portion of prion protein carrying oligomerization activity.

Simulations of membrane-bound diglycosylated human prion protein reveal potential protective mechanisms against misfolding

Prion diseases are associated with the misfolding of the prion protein (PrP) from its normal cellular form (PrP^C ) to its infectious scrapie form (PrP^Sc ). Post-translational modifications in PrP in vivo can play an important role in modulating the process of misfolding. To gain more insight into the effects of post-translational modifications in PrP structure and dynamics and to test the hypothesis that such modifications can interact with the protein, we have performed molecular dynamics simulations of diglycosylated human PrP^C bound to a lipid bilayer via a glycophosphatidylinositol anchor. Multiple simulations were performed at three different pH ranges to explore pH effects on structure and dynamics. In contrast to simulations of protein-only PrP^C , no large effects were observed upon lowering the pH of the system. The protein tilted toward the membrane surface in all of the simulations and the putative PrP^Sc oligomerization sites became inaccessible, thereby offering a possible protective mechanism against PrP^Sc -induced misfolding of PrP^C .

Mechanical and vibrational characterization of amyloid-like HET-s nanosheets based on the skewed plate theory

Pathological amyloidogenic prion proteins have a toxic effect on functional cells in the human cerebrum because of poor degradability and the tendency to accumulate in an uncontrolled manner under physiological conditions. HET-s, a fungal prion protein, is known to undergo conformational variations from fibrillar to nanosheet structures during a change from low to high pH conditions. It has been said that this conformational change can lead to self-propagation by nucleating on the lateral surface of singlet fibrils. Efforts have been made toward the mechanical characterization of fibrillar amyloids, but a global understanding of amyloid-like HET-s nanosheet structures is lacking. In this study, we analyzed the mechanical and vibrational characteristics of the skewed HET-s nanosheet structures that developed under neutral pH conditions by performing various molecular dynamics simulations. By applying the skewed plate theory to HET-s nanosheets for various length scales with numerous pores inside the structures, we found that the skewed HET-s nanosheet structure has mechanical properties comparable to those of previously reported biological film materials and nanomaterials. Considering the inherent characteristics of structural stability, our observation provides valuable and detailed structural information on skewed amyloid-like HET-s nanosheets.

Structural insight into the antiprion compound inhibition mechanism of native prion folding over misfolding

Transition of a physiological folded prion (PrP^C ) into a pathogenic misfolded prion (PrP^Sc ) causes lethal neurodegenerative disorders and prion diseases. Antiprion compounds have been developed to prevent this conversion; however, their mechanism of action remains unclear. Recently, we reported two antiprion compounds, BMD29 and BMD35, identified by in silico and in vitro screening. In this study, we used extensive explicit-solvent molecular dynamics simulations to investigate ligand-binding inhibition by antiprion compounds in prion folding over misfolding behavior at acidic pH. The two antiprion compounds and the previously reported GN8 compound resulted in a remarkably stabilized intermediate by binding to the hotspot region of PrP^C , whereas free PrP^C and the inactive compound BMD01 destabilized the structure of PrP^C leading to the misfolded form. The results uncovered a secondary structural transition of free PrP^C and transition suppression by the antiprion compounds. One of the major misfolding processes in PrP^C , alternation of hydrophobic core residues, disruption of intramolecular interactions, and the increase in residue solvent exposure were significantly inhibited by both antiprion compounds. These findings provide insights into prion misfolding and inhibition by antiprion compounds.

Mechanism of Unfolding of Human Prion Protein

Misfolding and aggregation of prion proteins are associated with several neurodegenerative diseases. Therefore, understanding the mechanism of the misfolding process is of enormous interest in the scientific community. It has been speculated and widely discussed that the native cellular prion protein (PrP^C) form needs to undergo substantial unfolding to a more stable PrP^C* state, which may further oligomerize into the toxic scrapie (PrP^Sc) form. Here, we have studied the mechanism of the unfolding of the human prion protein (huPrP) using a set of extensive well-tempered metadynamics simulations. Through multiple microsecond-long metadynamics simulations, we find several possible unfolding pathways. We show that each pathway leads to an unfolded state of lower free energy than the native state. Thus, our study may point to the signature of a PrP^C* form that corresponds to a global minimum on the conformational free-energy landscape. Moreover, we find that these global minima states do not involve an increased β-sheet content, as was assumed to be a signature of PrP^Sc formation in previous simulation studies. We have further analyzed the origin of metastability of the PrP^C form through free-energy surfaces of the chopped helical segments to show that the helices, particularly H2 and H3 of the prion protein, have the tendency to form either a random coil or a β-structure. Therefore, the secondary structural elements of the prion protein are only weakly stabilized by tertiary contacts and solvation forces so that relatively weak perturbations induced by temperature, pressure, pH, and so forth can lead to substantial unfolding with characteristics of intrinsically disordered proteins.

Elucidating the pH-Dependent Structural Transition of T7 Bacteriophage Endolysin

Bacteriophages are the most abundant and diverse biological entities on earth. Bacteriophage endolysins are unique peptidoglycan hydrolases and have huge potential as effective enzybiotics in various infectious models. T7 bacteriophage endolysin (T7L), also known as N-acetylmuramoyl-l-alanine amidase or T7 lysozyme, is a 17 kDa protein that lyses a range of Gram-negative bacteria by hydrolyzing the amide bond between N-acetylmuramoyl residues and the l-alanine of the peptidoglycan layer. Although the activity profiles of several of the T7 family members have been known for many years, the molecular basis for their pH-dependent differential activity is not clear. In this study, we explored the pH-induced structural, stability, and activity characteristics of T7L by applying a variety of biophysical techniques and protein nuclear magnetic resonance (NMR) spectroscopy. Our studies established a reversible structural transition of T7L below pH 6 and the formation of a partially denatured conformation at pH 3. This low-pH conformation is thermally stable and exposed its hydrophobic pockets. Further, NMR relaxation measurements and structural analysis unraveled that T7L is highly dynamic in its native state and a network of His residues are responsible for the observed pH-dependent conformational dynamics and transitions. As bacteriophage chimeric and engineered endolysins are being developed as novel therapeutics against multiple drug resistance pathogens, we believe that our results are of great help in designing these entities as broadband antimicrobial and/or antibacterial agents.

Pathogenic mutation R959W alters recognition dynamics of dysferlin inner DysF domain

Dysferlin, a 220 kD protein, plays a major role in regulating plasma membrane repair in muscle cells. Mutations in the dysferlin inner DysF domain are known to cause different types of muscular dystrophy, including limb-girdle muscular dystrophy type 2B (LGMD2B) and Miyoshi myopathy (MM). Replacement of arginine in position 959 by tryptophan has been frequently associated with both LGMD2B and MM, but the molecular mechanisms by which this mutation alters dysferlin function remain unknown. In this study, we have used protein binding site predictions and microsecond molecular dynamics (MD) simulations to determine the effects pathogenic mutation R959W on the structural dynamics of dysferlin inner DysF domain. Analysis of 2 μs long MD trajectories revealed that mutation R959W does not induce local destabilization, unfolding or misfolding of the domain. We used a binding site predictor to discover a protein-binding site (residues T958-I966 and E1031-H1037) that resembles pincers in shape. Cartesian principal component analysis and interresidue distance distributions of the wild-type domain showed that the predicted protein-binding site undergoes a pincer motion, and populates two structural states, open and closed. We found that mutation R959W inhibits the pincer motion of the protein-binding site and completely shifts the equilibrium toward the open state. These differences in the structural dynamics of the predicted binding site suggest that mutation R959W alters recognition dynamics of the inner DysF domain. Based on these findings and on previous experimental studies, we propose a novel role for the inner DysF domain in muscle membrane repair through recruitment of dysferlin to plasma membrane. In conclusion, these findings have important implications for our understanding of the structural aspects of muscular dystrophies in atomic-level resolution.

Perturbations in inter-domain associations may trigger the onset of pathogenic transformations in PrP(C): insights from atomistic simulations

Conversion of the predominantly α-helical cellular prion protein (PrP(C)) to the misfolded β-sheet enriched Scrapie form (PrP(Sc)) is a critical event in prion pathogenesis. However, the conformational triggers that lead to the isoform conversion (PrP(C) to PrP(Sc)) remain obscure, and conjectures about the role of unusually hydrophilic, short helix H1 of the C-terminal globular domain in the transition are varied. Helix H1 is anchored to helix H3 via a few stabilizing polar interactions. We have employed fully atomistic molecular dynamics simulations to study the effects triggered by a minor perturbation in the network of these non-bonded interactions in PrP(C). The elimination of just one of the key H1-H3 hydrogen bonds led to a cascade of conformational changes that are consistent with those observed in partially unfolded intermediates of PrP(C), with pathogenic mutations and in low pH environments. Our analyses reveal that the perturbation results in the enhanced conformational flexibility of the protein. The resultant enhancement in the dynamics leads to overall increased solvent exposure of the hydrophobic core residues and concomitant disruption of the H1-H3 inter-domain salt bridge network. This study lends credence to the hypothesis that perturbing the cooperativity of the stabilizing interactions in the PrP(C) globular domain can critically affect its dynamics and may lead to structural transitions of pathological relevance.

Molecular dynamics studies on the buffalo prion protein

It was reported that buffalo is a low susceptibility species resisting to transmissible spongiform encephalopathies (TSEs) (same as rabbits, horses, and dogs). TSEs, also called prion diseases, are invariably fatal and highly infectious neurodegenerative diseases that affect a wide variety of species (except for rabbits, dogs, horses, and buffalo), manifesting as scrapie in sheep and goats; bovine spongiform encephalopathy (BSE or "mad-cow" disease) in cattle; chronic wasting disease in deer and elk; and Creutzfeldt-Jakob diseases, Gerstmann-Sträussler-Scheinker syndrome, fatal familial insomnia, and Kulu in humans etc. In molecular structures, these neurodegenerative diseases are caused by the conversion from a soluble normal cellular prion protein (PrP(C)), predominantly with α-helices, into insoluble abnormally folded infectious prions (PrP(Sc)), rich in β-sheets. In this article, we studied the molecular structure and structural dynamics of buffalo PrP(C) (BufPrP(C)), in order to understand the reason why buffalo is resistant to prion diseases. We first did molecular modeling of a homology structure constructed by one mutation at residue 143 from the NMR structure of bovine and cattle PrP(124-227); immediately we found that for BufPrP(C)(124-227), there are five hydrogen bonds (HBs) at Asn143, but at this position, bovine/cattle do not have such HBs. Same as that of rabbits, dogs, or horses, our molecular dynamics studies also revealed there is a strong salt bridge (SB) ASP178-ARG164 (O-N) keeping the β2-α2 loop linked in buffalo. We also found there is a very strong HB SER170-TYR218 linking this loop with the C-terminal end of α-helix H3. Other information, such as (i) there is a very strong SB HIS187-ARG156 (N-O) linking α-helices H2 and H1 (if mutation H187R is made at position 187, then the hydrophobic core of PrP(C) will be exposed (L.H. Zhong (2010). Exposure of hydrophobic core in human prion protein pathogenic mutant H187R. Journal of Biomolecular Structure and Dynamics 28(3), 355-361)), (ii) at D178, there is a HB Y169-D178 and a polar contact R164-D178 for BufPrP(C) instead of a polar contact Q168-D178 for bovine PrP(C) (C.J. Cheng, & V. Daggett. (2014). Molecular dynamics simulations capture the misfolding of the bovine prion protein at acidic pH. Biomolecules 4(1), 181-201), (iii) BufPrP(C) owns three 310 helices at 125-127, 152-156, and in the β2-α2 loop, respectively, and (iv) in the β2-α2 loop, there is a strong π-π stacking and a strong π-cation F175-Y169-R164.(N)NH2, has been discovered.

Molecular dynamics studies on the NMR structures of rabbit prion protein wild type and mutants: surface electrostatic charge distributions

Prion diseases are invariably fatal and highly infectious neurodegenerative diseases that affect a wide variety of mammalian species such as sheep and goats, cattle, deer and elk, and humans. But for rabbits, studies have shown that they have a low susceptibility to be infected by prion diseases. This paper does molecular dynamics (MD) studies of rabbit NMR structures (of the wild type and its two mutants of two surface residues), in order to understand the specific mechanism of rabbit prion proteins (RaPrP(C)). Protein surface electrostatic charge distributions are specially focused to analyze the MD trajectories. This paper can conclude that surface electrostatic charge distributions indeed contribute to the structural stability of wild-type RaPrP(C); this may be useful for the medicinal treatment of prion diseases.