Characterization of cell-surface prion protein relative to its recombinant analogue: insights from molecular dynamics simulations of diglycosylated, membrane-bound human prion protein

The prion protein (PrP) is responsible for several fatal neurodegenerative diseases via conversion from its normal to disease-related isoform. The recombinant form of the protein is typically studied to investigate the conversion process. This constructs lacks the co- and post-translational modifications present in vivo , there the protein has two N-linked glycans and is bound to the outer leaflet of the plasma membrane via a glycosylphosphatidylinositol (GPI) anchor. The inherent flexibility and heterogeneity of the glycans, the plasticity of the GPI anchor, and the localization of the protein in a membrane make experimental structural characterization of biological constructs of cellular prion protein (PrP^C) challenging. Yet this characterization is central in determining not only the suitability of recombinant (rec)-PrP^C as a model for biological forms of the protein but also the potential role of co- and post-translational modifications on the disease process. Here, we present molecular dynamics simulations of three human prion protein constructs: (i) a protein-only construct modeling the recombinant form, (ii) a diglycosylated and soluble construct, and (iii) a diglycosylated and GPI-anchored construct bound to a lipid bilayer. We found that glycosylation and membrane anchoring do not significantly alter the structure or dynamics of PrP^C, but they do appreciably modify the accessibility of the polypeptide surface PrP^C. In addition, the simulations of membrane-bound PrP^C revealed likely recognition domains for the disease-initiating PrP^C:PrP^Sc (infectious and/or misfolded form of the prion protein) binding event and a potential mechanism for the observed inefficiency of conversion associated with differentially glycosylated PrP species.     

The pH dependence of flavivirus envelope protein structure: insights from molecular dynamics simulations

The flavivirus membrane fusion is triggered by the acid pH of the endosomes after virus endocytosis. The proposed mechanism involves changes in the protonation state of conserved histidine residues of the E protein present in the viral surface that undergoes a series of structural rearrangements that result in the fusion between the endosome and viral bilayers. We studied the pH dependence of E protein rearrangements of dengue virus type 2, used as a model, in the pH range experimented by the virus along the fusion process. We employed a low computational cost scheme to explore the behavior of the E protein by molecular dynamics (MD) simulations of complete systems that include the protein, the solvent, and ions. The procedure alternates cyclically the update of the ionization states of the protein residues with common MD steps applied to the new ionization configuration. Important pH-dependent protein structure rearrangements consistent with the changes of the protonation states of conserved histidine residues were observed. The involvement of other conserved residues in the flavivirus in the rearrangements was also identified. The results show interesting correlations with a proposed model for the fusion mechanism, as well as the experimentally identified key residues, contributing to a better understanding of the structural changes in protein E that lead to the fusion process.     

Sequence dependence of amyloid fibril formation: insights from molecular dynamics simulations

The clarification of the physico-chemical determinants underlying amyloid deposition is critical for our understanding of misfolding diseases. With this purpose we have performed a systematic all-atom molecular dynamics (MD) study of a series of single point mutants of the de novo designed amyloidogenic peptide STVIIE. Sixteen different 50ns long simulations using explicit solvent have been carried out starting from four different conformations of a polymeric six-stranded beta-sheet. The simulations have provided evidence for the influence of a small number of site-specific hydrophobic interactions on the packing and stabilization of nascent aggregates, as well as the interplay between side-chain interactions and the net charge of the molecule on the strand arrangement of polymeric beta-sheets. This MD analysis has also shed light into the origin of the position dependence on mutation of beta-sheet polymerization that was found experimentally for this model system. Our results suggest that MD can be applied to detect critical positions for beta-sheet aggregation within a given amyloidogenic stretch. Studies similar to the one presented here can guide site-directed mutations or the design of drugs that specifically disrupt the key stabilizing interactions of beta-sheet aggregates.     

The determinants of stability in the human prion protein: insights into folding and misfolding from the analysis of the change in the stabilization energy distribution in different conditions

The dynamic evolution of the PrP(C) from its NMR-derived conformation to a beta-sheet-rich, aggregation-prone conformation is studied through all-atom, explicit solvent molecular dynamics in different temperature and pH conditions. The trajectories are analyzed by means of a recently introduced energy decomposition approach aimed at identifying the key residues for the stabilization and folding of the protein. It is shown that under native conditions the stabilization energy is concentrated in regions of the helices H1 and H3, whereas under misfolding conditions (low pH, high temperature, or mutations in selected sites) it is spread out over helix H2. Misfolding appears to be a rearrangement of the chain that disrupts most of the native secondary structure of the protein, producing some beta-rich conformations with an energy distribution similar to that of the native state.     

CGDB: A database of membrane protein/lipid interactions by coarse-grained molecular dynamics simulations

Membrane protein function and stability has been shown to be dependent on the lipid environment. Recently, we developed a high-throughput computational approach for the prediction of membrane protein/lipid interactions. In the current study, we enhanced this approach with the addition of a new measure of the distortion caused by membrane proteins on a lipid bilayer. This is illustrated by considering the effect of lipid tail length and headgroup charge on the distortion caused by the integral membrane proteins MscS and FLAP, and by the voltage sensing domain from the channel KvAP. Changing the chain length of lipids alters the extent but not the pattern of distortion caused by MscS and FLAP; lipid headgroups distort in order to interact with very similar but not identical regions in these proteins for all bilayer widths investigated. Introducing anionic lipids into a DPPC bilayer containing the KvAP voltage sensor does not affect the extent of bilayer distortion.     

Molecular dynamics study on the stability of wild-type and the R220K mutant of human prion protein

Prion diseases are invariably fatal and highly infectious neurodegenerative diseases related to the structure transition of α-helix into β-sheet. In order to gain more direct insight into the molecular basis of the disease, the stability of the wild-type human prion protein (hPrPc) and the R220K mutant (m-hPrPc) was studied by molecular dynamics (MD) and flow MD simulation. Both the thermodynamic stability and the mechanical properties of hPrPc were investigated in this work. It was found that β-sheet was more readily to be unfolded in m-hPrPc. In the case of hPrPc, less content of helix was preserved after water turbulence. The H-bond network formed by the mutation-related residue 220 was found to play a key role in the stability of hPrPc.     

Folding and stability of the three-stranded beta-sheet peptide Betanova: insights from molecular dynamics simulations

The dynamics of the three-stranded beta-sheet peptide Betanova has been studied at four different temperatures (280, 300, 350, and 450 K by molecular dynamics simulation techniques, in explicit water. Two 20-ns simulations at 280 K indicate that the peptide remains very flexible under "folding" conditions sampling a range of conformations that together satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints. Two simulations at 300 K (above the experimental folding temperature) of 20 ns each show partial formation of "native"-like structure, which also satisfies most of the NOE constraints at 280 K. At higher temperature, the presence of compact states, in which a series of hydrophobic contacts remain present, are observed. This is consistent with experimental observations regarding the role of hydrophobic contacts in determining the peptide's stability and in initiating the formation of turns and loops. A set of different structures is shown to satisfy NMR-derived distance restraints and a possible mechanism for the folding of the peptide into the NMR-determined structure is proposed.     

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 3₁₀ 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.     

The dominant-negative effect of the Q218K variant of the prion protein does not require protein X

Previous studies identified several single-point mutants of the prion protein that displayed dominant-negative effects on prion replication. The dominant-negative effect was assumed to be mediated by protein X, an as-yet-unknown cellular cofactor that is believed to be essential for prion replication. To gain insight into the mechanism that underlies the dominant-negative phenomena, we evaluated the effect of the Q218K variant of full-length recombinant prion protein (Q218K rPrP), one of the dominant-negative mutants, on cell-free polymerization of wild-type rPrP into amyloid fibrils. We found that both Q218K and wild-type (WT) rPrPs were incorporated into fibrils when incubated as a mixture; however, the yield of polymerization was substantially decreased in the presence of Q218K rPrP. Furthermore, in contrast to fibrils produced from WT rPrP, the fibrils generated in the mixture of WT and Q218K rPrPs did not acquire the proteinase K-resistant core of 16 kDa that was shown previously to encompass residues 97-230 and was similar to that of PrP(Sc). Our studies demonstrate that the Q218K variant exhibits the dominant-negative effect in cell-free conversion in the absence of protein X, and that this effect is, presumably, mediated by physical interaction between Q218K and WT rPrP during the polymerization process.     

Quaternary structure effects on the hexacoordination equilibrium in rice hemoglobin rHb1: Insights from molecular dynamics simulations

Nonsymbiotic hemoglobins (nsHbs) form a widely distributed class of plant proteins, which function remains unknown. Despite the fact that class 1 plant nonsymbiotic hemoglobins are hexacoordinate (6c) heme proteins (hxHbs), their hexacoordination equilibrium constants are much lower than in hxHbs from animals or bacteria. In addition, they are characterized by having very high oxygen affinities and low oxygen dissociation rate constants. Rice hemoglobin 1 (rHb1) is a class 1 nonsymbiotic hemoglobin. It crystallizes as a fully associated homodimer with both subunits in 6c state, but showing slightly different conformations, thus leading to an asymmetric crystallographic homodimer. The residues that constitute the dimeric interface are conserved among all nsHbs, suggesting that the quaternary structure could be relevant to explain the chemical behavior and biological function of this family of proteins. In this work, we analyze the molecular basis that determine the hexacoordination equilibrium in rHb1. Our results indicate that dynamical features of the quaternary structure significantly affect the hexacoordination process. Specifically, we observe that the pentacoordinate state is stabilized in the dimer with respect to the isolated monomers. Moreover, the dimer behaves asymmetrically, in a negative cooperative scheme. The results presented in this work are fully consistent with our previous hypothesis about the key role played by the nature of the CD region in determining the coordination state of globins. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

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.     

Studies on the structural stability of rabbit prion probed by molecular dynamics simulations of its wild-type and mutants

Prion diseases are invariably fatal and highly infectious neurodegenerative diseases that affect humans and animals. Rabbits are the only mammalian species reported to be resistant to infection from prion diseases isolated from other species (Vorberg et al., 2003). Fortunately, the NMR structure of rabbit prion (124-228) (PDB entry 2FJ3), the NMR structure of rabbit prion protein mutation S173N (PDB entry 2JOH) and the NMR structure of rabbit prion protein mutation I214V (PDB entry 2JOM) were released recently. This paper studies these NMR structures by molecular dynamics simulations. Simulation results confirm the structural stability of wild-type rabbit prion, and show that the salt bridge between D177 and R163 greatly contributes to the structural stability of rabbit prion protein.     

Comparison studies of the structural stability of rabbit prion protein with human and mouse prion proteins

Background: Prion diseases are fatal and infectious neurodegenerative diseases affecting humans and animals. Rabbits are one of the few mammalian species reported to be resistant to infection from prion diseases isolated from other species (I. Vorberg et al., Journal of Virology 77 (3) (2003) 2003-2009). Thus the study of rabbit prion protein structure to obtain insight into the immunity of rabbits to prion diseases is very important.Findings: The paper is a straight forward molecular dynamics simulation study of wild-type rabbit prion protein (monomer cellular form) which apparently resists the formation of the scrapie form. The comparison analyses with human and mouse prion proteins done so far show that the rabbit prion protein has a stable structure. The main point is that the enhanced stability of the C-terminal ordered region especially helix 2 through the D177-R163 salt-bridge formation renders the rabbit prion protein stable. The salt bridge D201-R155 linking helixes 3 and 1 also contributes to the structural stability of rabbit prion protein. The hydrogen bond H186-R155 partially contributes to the structural stability of rabbit prion protein.Conclusions: Rabbit prion protein was found to own the structural stability, the salt bridges D177-R163, D201-R155 greatly contribute and the hydrogen bond H186-R155 partially contributes to this structural stability. The comparison of the structural stability of prion proteins from the three species rabbit, human and mouse showed that the human and mouse prion protein structures were not affected by the removing these two salt bridges. Dima et al. (Biophysical Journal 83 (2002) 1268-1280 and Proceedings of the National Academy of Sciences of the United States of America 101 (2004) 15335-15340) also confirmed this point and pointed out that “correlated mutations that reduce the frustration in the second half of helix 2 in mammalian prion proteins could inhibit the formation of PrP^Sc”.     

Reversible aggregation of mouse prion protein derivatives with PrPSC-like structural properties

Three carbamylated derivatives of reduced mouse prion protein (mPrP) were isolated during the aborted oxidative folding in the presence of urea. These three prion protein derivatives (mPrP-a, mPrP-b, and mPrP-c) exist as monomer in the acidic solution (pH < 2.0) and exhibit prevalent random coil structure. However, they undergo rapid aggregation and transformation to a predominant -sheet structure upon exposure to ionic buffer with pH greater than 3.0. The stability of aggregates of mPrP conformers is in part dependent upon the time that they were allowed to develop. The nascent aggregates comprise a significant fraction of loosely packed mPrP monomers that can be dissociated by treatment with strong acidic solution. Matured aggregates acquired through prolonged sample incubation contain more tightly packed mPrP monomers that cannot be dissociated by strong acid but can be disaggregated by denaturant. The properties of reversible aggregation of mPrP-a, mPrP-b, and mPrP-c bear a striking resemblance to that observed with aggregates of hamster PrP^SC.     

Partition of protein solvation into group contributions from molecular dynamics simulations

Linear response theory coupled to molecular dynamics simulations with an explicit solvent representation is used to derive fractional contributions of amino acid residues to the solvation of proteins. The new fractional methods developed here are compared with standard approaches based on empirical 1D and 3D statistical potentials, as well as with estimates obtained from the analysis of classical molecular interaction potentials. The new fractional methods, which have a clear physical basis and explicitly account for the effects due to protein structure and flexibility, provide an accurate picture of the contribution to solvation of different regions of the protein.     

Effects of ligand binding on the dynamics of rice nonspecific lipid transfer protein 1: a model from molecular simulations

Plant nonspecific lipid transfer proteins (nsLTPs) are small, basic proteins constituted mainly of alpha-helices and stabilized by four conserved disulfide bridges. They are characterized by the presence of a tunnel-like hydrophobic cavity, capable of transferring various lipid molecules between lipid bilayers in vitro. In this study, molecular dynamics (MD) simulations were performed at room temperature to investigate the effects of lipid binding on the dynamic properties of rice nsLTP1. Rice nsLTP1, either in the free form or complexed with one or two lipids was subjected to MD simulations. The C-terminal loop was very flexible both before and after lipid binding, as revealed by calculating the root-mean-square fluctuation. After lipid binding, the flexibility of some residues that were not in direct contact with lipid molecules increased significantly, indicating an increase of entropy in the region distal from the binding site. Essential dynamics analysis revealed clear differences in motion between unliganded and liganded rice nsLTP1s. In the free form of rice nsLTP1, loop1 exhibited the largest directional motion. This specific essential motion mode diminished after binding one or two lipid molecules. To verify the origin of the essential motion observed in the free form of rice nsLTP1, we performed multiple sequence alignments to probe the intrinsic motion encoded in the primary sequence. We found that the amino acid sequence of loop1 is highly conserved among plant nsLTP1s, thus revealing its functional importance during evolution. Furthermore, the sequence of loop1 is composed mainly of amino acids with short side chains. In this study, we show that MD simulations, together with essential dynamics analysis, can be used to determine structural and dynamic differences of rice nsLTP1 upon lipid binding.