One gene, two diseases and three conformations: Molecular dynamics simulations of mutants of human prion protein at room temperature and elevated temperatures

Effects of the pathological E200K mutation on human prion protein: A computational screening and molecular dynamics approach

The human prion protein gene (PRNP) is mapped to the short arm of chromosome 20 (20pter-12). Prion disease is associated with mutations in the prion protein-encoding gene sequence. Earlier studies found that the mutation G127V in the PRNP increases protein stability. In contrast, the mutation E200K, which has the highest mutation rate in the prion protein, causes Creutzfeldt-Jakob disease (CJD) in humans and induces protein aggregation. We aimed to identify the structural mechanisms of E200k and G127V mutations causing CJD. We used a variety of bioinformatic algorithms, including SIFT, PolyPhen, I-Mutant, PhD-SNP, and SNP& GO, to predict the association of the E200K mutation with prion disease. MD simulation is performed, and graphs for root mean square deviation, root mean square fluctuation, radius of gyration, DSSP, principal component analysis, porcupine, and free energy landscape are generated to confirm and prove the stability of the wild-type and mutant protein structures. The protein is analyzed for aggregation, and the results indicate more fluctuations in the protein structure during the simulation owing to the E200K mutation; however, the G127V mutation makes the protein structure stable against aggregation during the simulation.

Human prion protein: exploring the thermodynamic stability and structural dynamics of its pathogenic mutants

Human familial prion diseases are known to be associated with different single-point mutants of the gene coding for prion protein with a primary focus at several locations of the globular domain. We have identified 12 different single-point pathogenic mutants of human prion protein (HuPrP) with the help of extensive perturbations/mutation technique at multiple locations of HuPrP sequence related to potentiality towards conformational disorders. Among these, some of the mutants include pathogenic variants that corroborate well with the literature reported proteins while majority include some unique single-point mutants that are either not explicitly studied early or studied for variants with different residues at the specific position. Primarily, our study sheds light on the unfolding mechanism of the above mentioned mutants in depth. Besides, we could identify some mutants under investigation that demonstrates not only unfolding of the helical structures but also extension and generation of the β-sheet structures and or simultaneously have highly exposed hydrophobic surface which is assumed to be linked with the production of aggregate/fibril structures of the prion protein. Among the identified mutants, Q212E needs special attention due to its maximum exposure of hydrophobic core towards solvent and E200Q is found to be important due to its maximum extent of β-content. We are also able to identify different respective structural conformations of the proteins according to their degree of structural unfolding and those conformations can be extracted and further studied in detail. Communicated by Ramaswamy H. Sarma.

Advances in the understanding of protein misfolding and aggregation through molecular dynamics simulation

Aberrant protein folding known as protein misfolding is counted as one of the striking factors of neurodegenerative diseases. The extensive range of pathologies caused by protein misfolding, aggregation and subsequent accumulation are mainly classified into either gain of function diseases or loss of function diseases. In order to seek for novel strategies for treatment and diagnosis of neurodegenerative diseases, insights into the mechanism of misfolding and aggregation is essential. A comprehensive knowledge on the factors influencing misfolding and aggregation is required as well. An extensive experimental study on protein aggregation is somewhat challenging due to the insoluble and noncrystalline nature of amyloid fibrils. Thus there has been a growing use of computational approaches including Monte Carlo simulation, docking simulation, molecular dynamics simulation in the study of protein misfolding and aggregation. The review presents a discussion on molecular dynamics simulation alone as to how it has emerged as a promising tool in the understanding of protein misfolding and aggregation in general, detailing upon three different aspects considering four misfold prone proteins in particular. It is noticeable that all four proteins considered in this review i.e prion, superoxide dismutase1, huntingtin and amyloid β are linked to chronic neurodegenerative diseases with debilitating effects. Initially the review elaborates on the factors influencing the misfolding and aggregation. Next, it addresses our current understanding of the amyloid structures and the associated aggregation mechanisms, finally, summarizing the contribution of this computational tool in the search for therapeutic strategies against the respective protein-deposition diseases.

Human Aquaporin-4 and Molecular Modeling: Historical Perspective and View to the Future

Among the different aquaporins (AQPs), human aquaporin-4 (hAQP4) has attracted the greatest interest in recent years as a new promising therapeutic target. Such a membrane protein is, in fact, involved in a multiple sclerosis-like immunopathology called Neuromyelitis Optica (NMO) and in several disorders resulting from imbalanced water homeostasis such as deafness and cerebral edema. The gap of knowledge in its functioning and dynamics at the atomistic level of detail has hindered the development of rational strategies for designing hAQP4 modulators. The application, lately, of molecular modeling has proved able to fill this gap providing a breeding ground to rationally address compounds targeting hAQP4. In this review, we give an overview of the important advances obtained in this field through the application of Molecular Dynamics (MD) and other complementary modeling techniques. The case studies presented herein are discussed with the aim of providing important clues for computational chemists and biophysicists interested in this field and looking for new challenges.

Molecular Dynamics Simulations of the Temperature Induced Unfolding of Crambin Follow the Arrhenius Equation

Molecular dynamics simulations have been used extensively to model the folding and unfolding of proteins. The rates of folding and unfolding should follow the Arrhenius equation over a limited range of temperatures. This study shows that molecular dynamic simulations of the unfolding of crambin between 500K and 560K do follow the Arrhenius equation. They also show that while there is a large amount of variation between the simulations the average values for the rate show a very high degree of correlation.

Requirements for mutant and wild-type prion protein misfolding in vitro

Misfolding of the prion protein (PrP) plays a central role in the pathogenesis of infectious, sporadic, and inherited prion diseases. Here we use a chemically defined prion propagation system to study misfolding of the pathogenic PrP mutant D177N in vitro. This mutation causes PrP to misfold spontaneously in the absence of cofactor molecules in a process dependent on time, temperature, pH, and intermittent sonication. Spontaneously misfolded mutant PrP is able to template its unique conformation onto wild-type PrP substrate in a process that requires a phospholipid activity distinct from that required for the propagation of infectious prions. Similar results were obtained with a second pathogenic PrP mutant, E199K, but not with the polymorphic substitution M128V. Moreover, wild-type PrP inhibits mutant PrP misfolding in a dose-dependent manner, and cofactor molecules can antagonize this effect. These studies suggest that interactions between mutant PrP, wild-type PrP, and other cellular factors may control the rate of PrP misfolding in inherited prion diseases.

Structural and dynamic properties of the human prion protein

Prion diseases involve the conformational conversion of the cellular prion protein (PrP(C)) to its misfolded pathogenic form (PrP(Sc)). To better understand the structural mechanism of this conversion, we performed extensive all-atom, explicit-solvent molecular-dynamics simulations for three structures of the wild-type human PrP (huPrP) at different pH values and temperatures. Residue 129 is polymorphic, being either Met or Val. Two of the three structures have Met in position 129 and the other has Val. Lowering the pH or raising the temperature induced large conformational changes of the C-terminal globular domain and increased exposure of its hydrophobic core. In some simulations, HA and its preceding S1-HA loop underwent large displacements. The C-terminus of HB was unstable and sometimes partially unfolded. Two hydrophobic residues, Phe-198 and Met-134, frequently became exposed to solvent. These conformational changes became more dramatic at lower pH or higher temperature. Furthermore, Tyr-169 and the S2-HB loop, or the X-loop, were different in the starting structures but converged to common conformations in the simulations for the Met-129, but not the Val-129, protein. α-Strands and β-strands formed in the initially unstructured N-terminus. α-Strand propensity in the N-terminus was different between the Met-129 and Val129 proteins, but β-strand propensity was similar. This study reveals detailed structural and dynamic properties of huPrP, providing insight into the mechanism of the conversion of PrP(C) to PrP(Sc).

Reversibility of prion misfolding: insights from constant-pH molecular dynamics simulations

The prion protein (PrP) is the cause of a group of diseases known as transmissible spongiform encephalopathies (TSEs). Creutzfeldt-Jakob disease and bovine spongiform encephalopathy are examples of TSEs. Although the normal form of PrP (PrP(C)) is monomeric and soluble, it can misfold into a pathogenic form (PrP(Sc)) that has a high content of β-structure and can aggregate forming amyloid fibrils. The mechanism of conversion of PrP(C) into PrP(Sc) is not known but different triggers have been proposed. It can be catalyzed by a PrP(Sc) sample, or it can be induced by an external factor, such as low pH. The pH effect on the structure of PrP was recently studied by computational methods [Campos et al. J. Phys. Chem. B 2010, 114, 12692-12700], and an evident trend of loss of helical structure was observed with pH decrease, together with a gain of β-structures. In particular, one simulation at pH 2 showed an evident misfolding transition. The main goal of the present work was to study the effects of a change in pH to 7 in several transient conformations of this simulation, in order to draw some conclusions about the reversibility of PrP misfolding. Although the most significant effect caused by the change of pH to 7 was a global stabilization of the protein structure, we could also observe that some conformational transitions induced by pH 2 were reversible in many of our simulations, namely those started from the early moments of the misfolding transition. This observation is in good agreement with experiments showing that, even at pH as low as 1.7, it is possible to revert the misfolding process [Bjorndahl et al. Biochemistry 2011, 50, 1162-1173].

NMR structure and CD titration with metal cations of human prion alpha2-helix-related peptides

The 173–195 segment corresponding to the helix 2 of the C-globular prion protein domain could be one of several “spots” of intrinsic conformational flexibility. In fact, it possesses chameleon conformational behaviour and gathers several disease-associated point mutations. We have performed spectroscopic studies on the wild-type fragment 173–195 and on its D178N mutant dissolved in trifluoroethanol to mimic the in vivo system, both in the presence and in the absence of metal cations. NMR data showed that the structure of the D178N mutant is characterized by two short helices separated by a kink, whereas the wild-type peptide is fully helical. Both peptides retained these structural organizations, as monitored by CD, in the presence of metal cations. NMR spectra were however not in favour of the formation of definite ion-peptide complexes. This agrees with previous evidence that other regions of the prion protein are likely the natural target of metal cation binding.

Molecular dynamics as an approach to study prion protein misfolding and the effect of pathogenic mutations

Computer simulation of protein dynamics offers unique high-resolution information that complements experiment. Using experimentally derived structures of the natively folded prion protein (PrP), physically realistic dynamics and conformational changes can be simulated, including the initial steps of misfolding. By introducing mutations in silico, the effect of pathogenic mutations on PrP conformation and dynamics can be assessed. Here, we briefly introduce molecular dynamics methods and review the application of molecular dynamics simulations to obtain insight into various aspects of the PrP, including the mechanism of misfolding, the response to changes in the environment, and the influence of disease-related mutations.

Similarities between forms of sheep scrapie and Creutzfeldt-Jakob disease are encoded by distinct prion types

Transmissible spongiform encephalopathies such as scrapie in sheep, Creutzfeldt-Jakob disease (CJD) in humans, and bovine sporadic encephalopathy in cattle are characterized by the accumulation of a misfolded protein: the pathological prion protein. Ever since bovine sporadic encephalopathy was discovered as the likely cause of the new variant of CJD in humans, parallels between human and animal transmissible spongiform encephalopathies must be viewed under the aspect of a disease risk for humans. In our study we have compared prion characteristics of different forms of sheep scrapie with those of different phenotypes of sporadic CJD. The disease characteristics of sporadic CJD depend considerably on the prion type 1 or 2. Our results show that there are obvious parallels between sporadic CJD type 1 and the so-called atypical/Nor98 scrapie. These parelleles apply to the deposition form of pathological prion protein in the brain, detected by the paraffin-embedded-tissue blot and the prion aggregate stability with regard to denaturation by the chaotropic salt guanidine hydrochloride. The same applies to sporadic CJD type 2 and classical scrapie. The observed parallels between types of sporadic CJD and types of sheep scrapie demonstrate that distinct groups of prion disease exist in different species. This should be taken into consideration when discussing interspecies transmission.

Prion protein expression and processing in human mononuclear cells: the impact of the codon 129 prion gene polymorphism

Background

So far, all clinical cases of new variant Creutzfeldt-Jakob disease (vCJD), thought to result from the Bovine Spongiform Encephalopathy (BSE) prion agent, have shown Methionine–Methionine (M/M) homozygosity at the M129V polymorphism of the PRNP gene. Although established, this relationship is still not understood. In both vCJD and experimental BSE models prion agents do reach the bloodstream, raising concerns regarding disease transmission through blood transfusion.

Methodology/Principal Findings

We investigated the impact of the M129V polymorphism on the expression and processing of the prion protein in human peripheral blood mononuclear cells (PBMCs) from three blood donor populations with Methionine-Methionine (M/M), Valine-Valine (V/V) and M/V genotypes. Using real-time PCR, ELISA and immunoblot assays we were unable to find differences in prion protein expression and processing relating to the M129V polymorphism.

Conclusions/Significance

These results suggest that in PBMCs, the M129V PrP polymorphism has no significant impact on PrP expression, processing and the apparent glycoform distribution. Prion propagation should be investigated further in other cell types or tissues.

Structural characterization of a neurotoxic threonine-rich peptide corresponding to the human prion protein alpha 2-helical 180-195 segment, and comparison with full-length alpha 2-helix-derived peptides

The 173-195 segment corresponding to the helix 2 of the globular PrP domain is a good candidate to be one of the several 'spots' of intrinsic structural flexibility, which might induce local destabilization and concur to protein transformation, leading to aggregation-prone conformations. Here, we report CD and NMR studies on the alpha2-helix-derived peptide of maximal length (hPrP[180-195]) that is able to exhibit a regular structure different from the prevalently random arrangement of other alpha2-helix-derived peptides. This peptide, which has previously been shown to be affected by buffer composition via the ion charge density dependence typical of Hofmeister effects, corresponds to the C-terminal sequence of the PrP(C) full-length alpha2-helix and includes the highly conserved threonine-rich 188-195 segment. At neutral pH, its conformation is dominated by beta-type contributions, which only very strong environmental modifications are able to modify. On TFE addition, an increase of alpha-helical content can be observed, but a fully helical conformation is only obtained in neat TFE. However, linking of the 173-179 segment, as occurring in wild-type and mutant peptides corresponding to the full-length alpha2-helix, perturbs these intrinsic structural propensities in a manner that depends on whether the environment is water or TFE. Overall, these results confirm that the 180-195 parental region in hPrP(C) makes a strong contribution to the chameleon conformational behavior of the segment corresponding to the full-length alpha2-helix, and could play a role in determining structural rearrangements of the entire globular domain.

Dynamic alpha-helices: conformations that do not conform

Structural transitions are important for the stability and function of proteins, but these phenomena are poorly understood. An extensive analysis of Protein Data Bank entries reveals 103 regions in proteins with a tendency to transform from helical to nonhelical conformation and vice versa. We find that these dynamic helices, unlike other helices, are depleted in hydrophobic residues. Furthermore, the dynamic helices have higher surface accessibility and conformational mobility (P-value = 3.35e-07) than the rigid helices. Contact analyses show that these transitions result from protein-ligand, protein-nucleic acid, and crystal-contacts. The immediate structural environment differs quantitatively (P-value = 0.003) as well as qualitatively in the two alternate conformations. Often, dynamic helix experiences more contacts in its helical conformation than in the nonhelical counterpart (P-value = 0.001). There is differential preference for the type of short contacts observed in two conformational states. We also demonstrate that the regions in protein that can undergo such large conformational transitions can be predicted with a reasonable accuracy using logistic regression model of supervised learning. Our findings have implications in understanding the molecular basis of structural transitions that are coupled with binding and are important for the function and stability of the protein. Based on our observations, we propose that several functionally relevant regions on the protein surface can switch over their conformation from coil to helix and vice-versa, to regulate the recognition and binding of their partner and hence these may work as "molecular switches" in the proteins to regulate certain biological process. Our results supports the idea that protein structure-function paradigm should transform from static to a highly dynamic one.

Physiological role of the cellular prion protein

The prion protein (PrP) plays a key role in the pathogenesis of prion diseases. However, the normal function of the protein remains unclear. The cellular isoform (PrP(C)) is expressed most abundantly in the brain, but has also been detected in other non-neuronal tissues as diverse as lymphoid cells, lung, heart, kidney, gastrointestinal tract, muscle, and mammary glands. Cell biological studies of PrP contribute to our understanding of PrP(C) function. Like other membrane proteins, PrP(C) is post-translationally processed in the endoplasmic reticulum and Golgi on its way to the cell surface after synthesis. Cell surface PrP(C) constitutively cycles between the plasma membrane and early endosomes via a clathrin-dependent mechanism, a pathway consistent with a suggested role for PrP(C) in cellular trafficking of copper ions. Although PrP(-/-) mice have been reported to have only minor alterations in immune function, PrP(C) is up-regulated in T cell activation and may be expressed at higher levels by specialized classes of lymphocytes. Furthermore, antibody cross-linking of surface PrP(C) modulates T cell activation and leads to rearrangements of lipid raft constituents and increased phosphorylation of signaling proteins. These findings appear to indicate an important but, as yet, ill-defined role in T cell function. Recent work has suggested that PrP(C) is required for self-renewal of haematopoietic stem cells. PrP(C) is highly expressed in the central nervous system, and since this is the major site of prion pathology, most interest has focused on defining the role of PrP(C) in neurones. Although PrP(-/-) mice have a grossly normal neurological phenotype, even when neuronal PrP(C) is knocked out postnatally, they do have subtle abnormalities in synaptic transmission, hippocampal morphology, circadian rhythms, and cognition and seizure threshold. Other postulated neuronal roles for PrP(C) include copper-binding, as an anti- and conversely, pro-apoptotic protein, as a signaling molecule, and in supporting neuronal morphology and adhesion. The prion protein may also function as a metal binding protein such as copper, yielding cellular antioxidant capacity suggesting a role in the oxidative stress homeostasis. Finally, recent observations on the role of PrP(C) in long-term memory open a challenging field.

Beta-sheet containment by flanking prolines: molecular dynamic simulations of the inhibition of beta-sheet elongation by proline residues in human prion protein

Previous molecular dynamic simulations have reported elongation of the existing beta-sheet in prion proteins. Detailed examination has shown that these elongations do not extend beyond the proline residues flanking these beta-sheets. In addition, proline has also been suggested to possess a possible structural role in preserving protein interaction sites by preventing invasion of neighboring secondary structures. In this work, we have studied the possible structural role of the flanking proline residues by simulating mutant structures with alternate substitution of the proline residues with valine. Simulations showed a directional inhibition of elongation, with the elongation progressing in the direction of valine including evident inhibition of elongation by existing proline residues. This suggests that the flanking proline residues in prion proteins may have a containment role and would confine the beta-sheet within a specific length.

Does tetracycline bind helix 2 of prion? An integrated spectroscopical and computational study of the interaction between the antibiotic and alpha helix 2 human prion protein fragments

We demonstrate here that tetracycline (TC) can strongly interact (KD' = 189 +/- 7 nM) with model peptides derived from the C-terminal globular domain of the prion protein, hPrP [173-195], and that interaction concerns residues within the C-terminal half of the helix 2, a short region previously indicated as endowed with ambivalent conformational behavior and implicated in PrP conversion to the beta-sheet-rich, infective scrapie variant. Data have been confirmed by binding studies with the N-terminal truncated 180-195 variant that displays a dissociation constant of 483 +/- 30 nM. Remarkably, TC does not influence the structure of the N-terminally fluoresceinated peptides that both show alpha-helical conformations. Docking calculations and molecular dynamics simulations suggest a direct, strong interaction of the antibiotic with exposed side chain functional groups of threonines 190-193 on the solvent-exposed surface of helix 2.

Short-term memory formation and long-term memory consolidation are enhanced by cellular prion association to stress-inducible protein 1

Cellular prion protein (PrP(C)) is a cell surface glycoprotein that interacts with several ligands such as laminin, NCAM (Neural-Cell Adhesion Molecule) and the stress-inducible protein 1 (STI1). PrP(C) association with these proteins in neurons mediates adhesion, differentiation and protection against programmed cell death. Herein, we used an aversively motivated learning paradigm in rats to investigate whether STI1 interaction with PrP(C) affects short-term memory (STM) formation and long-term memory (LTM) consolidation. Blockage of PrP(C)-STI1 interaction with intra-hippocampal infusion of antibodies against PrP(C) or STI1 immediately after training impaired both STM and LTM. Furthermore, infusion of PrP(C) peptide 106-126, which competes for PrP(C)-STI1 interaction, also inhibited both forms of memory. Remarkably, STI1 peptide 230-245, which includes the PrP(C) binding site, had a potent enhancing effect on memory performance, which could be blocked by co-treatment with the competitive PrP(C) peptide 106-126. Taken together, these results demonstrate that PrP(C)-STI1 interaction modulates both STM and LTM and suggests a potential use of ST11 peptide 230-245 as a pharmacological agent.

Mutation analysis in the FKRP gene provides an explanation for a rare cause of intrafamilial clinical variability in LGMD2I

We report a limb-girdle muscular dystrophy 2I family with three affected sisters and a highly variable clinical course. FKRP gene sequencing showed that all three sisters carried a nonsense paternal mutation (W225X). The two oldest sisters with a severe phenotype carried two maternal mutations V79M and P89A. However, the youngest sister with a milder course carried the paternal and only the V79M maternal mutation, due to an intragenic recombination.

Effects of temperature and salt concentration on the structural stability of human lymphotactin: insights from molecular simulations

Extensive molecular dynamics (MD) simulations ( approximately 70 ns total) with explicit solvent molecules and salt ions are carried out to probe the effects of temperature and salt concentration on the structural stability of the human Lymphotactin (hLtn). The distribution of ions near the protein surface and the stability of various structural motifs are observed to exhibit interesting dependence on the local sequence and structure. Whereas chloride association to the protein is overall enhanced as the temperature increases, the sodium distribution in the C-terminal helical region and, to a smaller degree, the chloride distribution in the same region are found higher at the lower temperature. The similar trend is also observed in nonlinear Poisson-Boltzmann calculations with a temperature-dependent water dielectric constant, once conformational averaging over a series of MD snapshots is done. The unexpected temperature dependence in the ion distribution is explained on the basis of the cancellation of association entropy for ion-side chain pairs of opposite-charge and like-charge characters, which have positive and negative contributions, respectively. The C-terminal helix is observed to partially melt whereas a short beta strand forms at the higher temperature with little salt dependence. The N-terminal region, by contrast, develops partial helical structure at a higher salt concentration. These observed behaviors are consistent with solvent and salt screening playing an important role in stabilizing the canonical chemokine fold of hLtn.

Modeling of tumor necrosis factor receptor superfamily 1A mutants associated with tumor necrosis factor receptor–associated periodic syndrome indicates misfolding consistent with abnormal function

OBJECTIVE

To investigate the effect of mutations in the tumor necrosis factor receptor superfamily 1A (TNFRSF1A) gene on the conformation and behavior of the TNFRSF1A protein. Mutations in TNFRSF1A cause the autosomal-dominant, autoinflammatory TNFR-associated periodic syndrome (TRAPS).

METHODS

The expression of recombinant TNFRSF1A was compared in SK-HEp-1 endothelial cells and HEK 293 epithelial cells stably transfected with full-length R347A or Deltasig constructs of wild-type or TRAPS-associated mutant TNFRSF1A. TNF binding was assessed in HEK 293 cell lines expressing R347A wild-type or mutant TNFRSF1A. Homology modeling of the 3-dimensional structure of the ectodomains of wild-type and mutant TNFRSF1A was performed.

RESULTS

TRAPS-associated mutant and wild-type TNFRSF1A behaved differently and had different localization properties within the cell, as a direct result of mutations in the ectodomains of TNFRSF1A. From a structural perspective, mutants with a predicted structure similar to that of the wild-type protein (e.g., R92Q) behaved similarly to wild-type TNFRSF1A, whereas forms of TNFRSF1A with mutations predicted to drastically destabilize the protein structure (e.g., cysteine mutations) showed defects in cell surface expression and TNF binding.

CONCLUSION

The results obtained from the in vitro experiments, in combination with the modeled structures, indicate that the phenotype and clinical differences between different TRAPS-associated mutants of TNFRSF1A result from different conformations of the TNFRSF1A ectodomains.