Assessing the acid–base and conformational properties of histidine residues in human prion protein (125–228) by means of pK a calculations and molecular dynamics simulations

Histidine tautomerism-mediated transthyretin amyloidogenesis: A molecular insight

Characterization of the conformational alterations involved in monomer misfolding is essential for elucidating the molecular basis of the initial stage of protein accumulation. Here, we report the first structural analyses of transthyretin (TTR) (26-57) fragments with two histidine tautomeric states (δ; N_δ1H and ε; N_ε2H) using replica-exchange molecular dynamics (REMD) simulations. Explaining the organizational properties and misfolding procedure is challenging because the δ and ε configurations can occur in the free neutral state. REMD revealed that β-sheet generation is favored for the δδ (16.8%) and εδ (6.7%) tautomeric isomers, showing frequent main-chain contacts between the stable regions near the head (N-terminus) and central (middle) part compared to the εε (4.8%) and δε (2.8%) isomers. The presence of smaller and wider local energy minima may be related to the structural stability and toxicity of δδ/εδ and εε/δε. Histidines31 and 56 were the parts of regular (such as β-strand) and nonregular (such as coil) secondary structures within the highly toxic TTR isomer. For TTR amyloidosis, focusing on hazardous isomeric forms with high sheet contents may be a potent treatment strategy. Overall, our findings support the tautomerism concept and aid in our comprehension of the basic tautomeric actions of neutral histidine throughout the misfolding process.

Unraveling the Histidine Tautomerism Effect on the Initial Stages of Prion Misfolding: New Insights from a Computational Perspective

The aggregation and structural conversion of normal prion peptide (PrP^C) into the pathogenic scrapie form (PrP^Sc), which can act as a seed to enhance prion amyloid fiber formation, is believed to be a crucial event in prionopathies. Previous research suggests that the prion monomer may play an important role in oligomer generation during disease pathogenesis. In the present study, extensive replica-exchange molecular dynamics (REMD) simulations were conducted to explore the conformational characteristics of the huPrP (125-160) monomer under the histidine tautomerism effect. Investigating the structural characteristics and fibrilization process is challenging because two histidine tautomers [N_ε2-H (ε) and N_δ1-H (δ)] can occur in the open neutral state. Molecular dynamics (MD) simulation outcomes have shown that the toxic εδ and δδ isomer (containing several and broader local minima) had the highest α-helix structures, with contents of 21.11% and 21.01%, respectively, and may have a strong influence on the organizational behavior of a monomeric prion. The amino acids aspartate 20 (D20)-asparagine 29 (N29) and isoleucine 15 (I15)-histidine 16 (H16), D20-arginine 27 (R27) as well as N29 formed α-helix with the highest probabilities in the δδ and εδ isomer, accordingly. On the basis of our findings, we propose the histidine tautomerization hypothesis as a new prion accumulation mechanism, which may exist to induce the formation of prion accumulates. Overall, our tautomerism hypothesis constitutes a promising perspective for enhancing understanding of prion disease pathobiology and may help in the design of a good inhibitor.

Structural insight into conformational change in prion protein by breakage of electrostatic network around H187 due to its protonation

A conformational change from normal prion protein(PrP^C) to abnormal prion protein(PrP^SC) induces fatal neurodegenerative diseases. Acidic pH is well-known factors involved in the conformational change. Because the protonation of H187 is strongly linked to the change in PrP stability, we examined the charged residues R156, E196, and D202 around H187. Interestingly, there have been reports on pathological mutants, such as H187R, E196A, and D202N. In this study, we focused on how an acidic pH and pathological mutants disrupt this electrostatic network and how this broken network destabilizes PrP structure. To do so, we performed a temperature-based replica-exchange molecular dynamics (T-REMD) simulation using a cumulative 252 μs simulation time. We measured the distance between amino acids comprising four salt bridges (R156–E196/D202 and H187–E196/D202). Our results showed that the spatial configuration of the electrostatic network was significantly altered by an acidic pH and mutations. The structural alteration in the electrostatic network increased the RMSF value around the first helix (H1). Thus, the structural stability of H1, which is anchored to the H2–H3 bundle, was decreased. It induces separation of R156 from the electrostatic network. Analysis of the anchoring energy also shows that two salt-bridges (R156-E196/D202) are critical for PrP stability.

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.

The protonation state of histidine 111 regulates the aggregation of the evolutionary most conserved region of the human prion protein

In a group of neurodegenerative diseases, collectively termed transmissible spongiform encephalopathies, the prion protein aggregates into β-sheet rich amyloid-like deposits. Because amyloid structure has been connected to different prion strains and cellular toxicity, it is important to obtain insight into the structural properties of prion fibrils. Using a combination of solution NMR spectroscopy, thioflavin-T fluorescence and electron microscopy we here show that within amyloid fibrils of a peptide containing residues 108-143 of the human prion protein [humPrP (108-143)]-the evolutionary most conserved part of the prion protein - residue H111 and S135 are in close spatial proximity and their interaction is critical for fibrillization. We further show that residues H111 and H140 share the same microenvironment in the unfolded, monomeric state of the peptide, but not in the fibrillar form. While protonation of H140 has little influence on fibrillization of humPrP (108-143), a positive charge at position 111 blocks the conformational change, which is necessary for amyloid formation of humPrP (108-143). Our study thus highlights the importance of protonation of histidine residues for protein aggregation and suggests point mutations to probe the structure of infectious prion particles.

Unraveling the Molecular Mechanism of pH-Induced Misfolding and Oligomerization of the Prion Protein

The misfolding of the prion protein (PrP) to aggregated forms is linked to several neurodegenerative diseases. Misfolded oligomeric forms of PrP are associated with neurotoxicity and/or infectivity, but the molecular mechanism by which they form is still poorly understood. A reduction in pH is known to be a key factor that triggers misfolded oligomer formation by PrP, but the residues whose protonation is linked with misfolding remain unidentified. The structural consequences of the protonation of these residues also remain to be determined. In the current study, amino acid residues whose protonation is critical for PrP misfolding and oligomerization have been identified using site-directed mutagenesis and misfolding/oligomerization assays. It is shown that the protonation of either H186 or D201, which mimics the effects of pathogenic mutations (H186R and D201N) at both residue sites, is critically linked to the stability, misfolding and oligomerization of PrP. Hydrogen-deuterium exchange studies coupled with mass spectrometry show that the protonation of either H186 or D201 leads to the same common structural change: increased structural dynamics in helix 1 and that in the loop between helix 1 and β-strand 2. It is shown that the protonation of either of these residues is sufficient for accelerating misfolded oligomer formation, most likely because the protonation of either residue causes the same structural perturbation. Hence, the increased structural dynamics in helix 1 and that in the loop between helix 1 and β-strand 2 appear to play an early critical role in acid-induced misfolding of PrP.

Molecular Mechanism of the Misfolding and Oligomerization of the Prion Protein: Current Understanding and Its Implications

Prion diseases, also known as transmissible spongiform encephalopathies, make up a group of fatal neurodegenerative disorders linked with the misfolding and aggregation of the prion protein (PrP). Although it is not yet understood how the misfolding of PrP induces neurodegeneration, it is widely accepted that the formation of misfolded prion protein (termed PrP(Sc)) is both the triggering event in the disease and the main component of the infectious agent responsible for disease transmission. Despite the clear involvement of PrP(Sc) in prion diseases, the exact composition of PrP(Sc) is not yet well-known. Recent studies show that misfolded oligomers of PrP could, however, be responsible for neurotoxicity and/or infectivity in the prion diseases. Hence, understanding the molecular mechanism of formation of the misfolded oligomers of PrP is critical for developing an understanding about the prion diseases and for developing anti-prion therapeutics. This review discusses recent advances in understanding the molecular mechanism of misfolded oligomer formation by PrP and its implications for the development of anti-prion therapeutics.

Different misfolding mechanisms converge on common conformational changes: human prion protein pathogenic mutants Y218N and E196K

Prion diseases are caused by misfolding and aggregation of the prion protein (PrP). Pathogenic mutations such as Y218N and E196K are known to cause Gerstmann-Sträussler-Scheinker syndrome and Creutzfeldt-Jakob disease, respectively. Here we describe molecular dynamics simulations of these mutant proteins to better characterize the detailed conformational effects of these sequence substitutions. Our results indicate that the mutations disrupt the wild-type native PrP(C) structure and cause misfolding. Y218N reduced hydrophobic packing around the X-loop (residues 165-171), and E196K abolished an important wild-type salt bridge. While differences in the mutation site led PrP mutants to misfold along different pathways, we observed multiple traits of misfolding that were common to both mutants. Common traits of misfolding included: 1) detachment of the short helix (HA) from the PrP core; 2) exposure of side chain F198; and 3) formation of a nonnative strand at the N-terminus. The effect of the E196K mutation directly abolished the wild-type salt bridge E196-R156, which further destabilized the F198 hydrophobic pocket and HA. The Y218N mutation propagated its effect by increasing the HB-HC interhelical angle, which in turn disrupted the packing around F198. Furthermore, a nonnative contact formed between E221 and S132 on the S1-HA loop, which offered a direct mechanism for disrupting the hydrophobic packing between the S1-HA loop and HC. While there were common misfolding features shared between Y218N and E196K, the differences in the orientation of HB and HC and the X-loop conformation might provide a structural basis for identifying different prion strains.

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).

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

Bovine spongiform encephalopathy (BSE), or mad cow disease, is a fatal neurodegenerative disease that is transmissible to humans and that is currently incurable. BSE is caused by the prion protein (PrP), which adopts two conformers; PrPC is the native innocuous form, which is α-helix rich; and PrPSc is the β-sheet rich misfolded form, which is infectious and forms neurotoxic species. Acidic pH induces the conversion of PrPC to PrPSc. We have performed molecular dynamics simulations of bovine PrP at various pH regimes. An acidic pH environment induced conformational changes that were not observed in neutral pH simulations. Putative misfolded structures, with nonnative β-strands formed in the flexible N-terminal domain, were found in acidic pH simulations. Two distinct pathways were observed for the formation of nonnative β-strands: at low pH, hydrophobic contacts with M129 nucleated the nonnative β-strand; at mid-pH, polar contacts involving Q168 and D178 facilitated the formation of a hairpin at the flexible N-terminus. These mid- and low pH simulations capture the process of nonnative β-strand formation, thereby improving our understanding of how PrPC misfolds into the β-sheet rich PrPSc and how pH factors into the process.

Two misfolding routes for the prion protein around pH 4.5

Using molecular dynamics simulations, we show that the prion protein (PrP) exhibits a dual behavior, with two possible transition routes, upon protonation of H187 around pH 4.5, which mimics specific conditions encountered in endosomes. Our results suggest a picture in which the protonated imidazole ring of H187 experiences an electrostatic repulsion with the nearby guanidinium group of R136, to which the system responds by pushing either H187 or R136 sidechains away from their native cavities. The regions to which H187 and R136 are linked, namely the C-terminal part of H2 and the loop connecting S1 to H1, respectively, are affected in a different manner depending on which pathway is taken. Specific in vivo or in vitro conditions, such as the presence of molecular chaperones or a particular experimental setup, may favor one transition pathway over the other, which can result in very different [Formula: see text] monomers. This has some possible connections with the observation of various fibril morphologies and the outcome of prion strains. In addition, the finding that the interaction of H187 with R136 is a weak point in mammalian PrP is supported by the absence of the [Formula: see text] residue pair in non-mammalian species that are known to be resistant to prion diseases.

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].

Calculating pH-dependent free energy of proteins by using Monte Carlo protonation probabilities of ionizable residues

Protein folding, stability, and function are usually influenced by pH. And free energy plays a fundamental role in analysis of such pH-dependent properties. Electrostatics-based theoretical framework using dielectric solvent continuum model and solving Poisson-Boltzmann equation numerically has been shown to be very successful in understanding the pH-dependent properties. However, in this approach the exact computation of pH-dependent free energy becomes impractical for proteins possessing more than several tens of ionizable sites (e.g. > 30), because exact evaluation of the partition function requires a summation over a vast number of possible protonation microstates. Here we present a method which computes the free energy using the average energy and the protonation probabilities of ionizable sites obtained by the well-established Monte Carlo sampling procedure. The key feature is to calculate the entropy by using the protonation probabilities. We used this method to examine a well-studied protein (lysozyme) and produced results which agree very well with the exact calculations. Applications to the optimum pH of maximal stability of proteins and protein-DNA interactions have also resulted in good agreement with experimental data. These examples recommend our method for application to the elucidation of the pH-dependent properties of proteins.

The prion protein binds thiamine

Although highly conserved throughout evolution, the exact biological function of the prion protein is still unclear. In an effort to identify the potential biological functions of the prion protein we conducted a small-molecule screening assay using the Syrian hamster prion protein [shPrP(90-232)]. The screen was performed using a library of 149 water-soluble metabolites that are known to pass through the blood-brain barrier. Using a combination of 1D NMR, fluorescence quenching and surface plasmon resonance we identified thiamine (vitamin B1) as a specific prion ligand with a binding constant of ~60 μM. Subsequent studies showed that this interaction is evolutionarily conserved, with similar binding constants being seen for mouse, hamster and human prions. Various protein construct lengths, both with and without the unstructured N-terminal region in the presence and absence of copper, were examined. This indicates that the N-terminus has no influence on the protein's ability to interact with thiamine. In addition to thiamine, the more biologically abundant forms of vitamin B1 (thiamine monophosphate and thiamine diphosphate) were also found to bind the prion protein with similar affinity. Heteronuclear NMR experiments were used to determine thiamine's interaction site, which is located between helix 1 and the preceding loop. These data, in conjunction with computer-aided docking and molecular dynamics, were used to model the thiamine-binding pharmacophore and a comparison with other thiamine binding proteins was performed to reveal the common features of interaction.

Comparative analysis of essential collective dynamics and NMR-derived flexibility profiles in evolutionarily diverse prion proteins

Collective motions on ns-μs time scales are known to have a major impact on protein folding, stability, binding and enzymatic efficiency. It is also believed that these motions may have an important role in the early stages of prion protein misfolding and prion disease. In an effort to accurately characterize these motions and their potential influence on the misfolding and prion disease transmissibility we have conducted a combined analysis of molecular dynamic simulations and NMR-derived flexibility measurements over a diverse range of prion proteins. Using a recently developed numerical formalism, we have analyzed the essential collective dynamics (ECD) for prion proteins from 8 different species including human, cow, elk, cat, hamster, chicken, turtle and frog. We also compared the numerical results with flexibility profiles generated by the random coil index (RCI) from NMR chemical shifts. Prion protein backbone flexibility derived from experimental NMR data and from theoretical computations show strong agreement with each other, demonstrating that it is possible to predict the observed RCI profiles employing the numerical ECD formalism. Interestingly, flexibility differences in the loop between second beta strand (S2) and the second alpha helix (HB) appear to distinguish prion proteins from species that are susceptible to prion disease and those that are resistant. Our results show that the different levels of flexibility in the S2-HB loop in various species are predictable via the ECD method, indicating that ECD may be used to identify disease resistant variants of prion proteins, as well as the influence of prion proteins mutations on disease susceptibility or misfolding propensity.

Influence of pH on the human prion protein: insights into the early steps of misfolding

Transmissible spongiform encephalopathies, or prion diseases, are caused by misfolding and aggregation of the prion protein PrP. Conversion from the normal cellular form (PrP(C)) or recombinant PrP (recPrP) to a misfolded form is pH-sensitive, in that misfolding and aggregation occur more readily at lower pH. To gain more insight into the influence of pH on the dynamics of PrP and its potential to misfold, we performed extensive molecular-dynamics simulations of the recombinant PrP protein (residues 90-230) in water at three different pH regimes: neutral (or cytoplasmic) pH (∼7.4), middle (or endosomal) pH (∼5), and low pH (<4). We present five different simulations of 50 ns each for each pH regime, amounting to a total of 750 ns of simulation time. A detailed analysis and comparison with experiment validate the simulations and lead to new insights into the mechanism of pH-induced misfolding. The mobility of the globular domain increases with decreasing pH, through displacement of the first helix and instability of the hydrophobic core. At middle pH, conversion to a misfolded (PrP(Sc)-like) conformation is observed. The observed changes in conformation and stability are consistent with experimental data and thus provide a molecular basis for the initial steps in the misfolding process.

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.

Genetic Creutzfeldt-Jakob disease and fatal familial insomnia: insights into phenotypic variability and disease pathogenesis

Human prion diseases are a group of rare neurodegenerative disorders characterized by the conversion of the constitutively expressed prion protein, PrP(C), into an abnormally aggregated isoform, called PrP(Sc). While most people who develop a prion disease have no identifiable cause and a few acquire the disease through an identified source of infection, about 10-15% of patients are affected by a genetic form and carry either a point mutation or an insertion of octapeptide repeats in the prion protein gene. Prion diseases show the highest extent of phenotypic heterogeneity among neurodegenerative disorders and comprise three major disease entities with variable though overlapping phenotypic features: Creutzfeldt-Jakob disease (CJD), fatal insomnia and the Gerstmann-Sträussler-Scheinker syndrome. Both CJD and fatal insomnia are fully transmissible diseases, a feature that led to the isolation and characterization of different strains of the agent or prion showing distinctive clinical and neuropathological features after transmission to syngenic animals. Here, we review the current knowledge of the effects of the pathogenic mutations linked to genetic CJD and fatal familial insomnia on the prion protein metabolism and physicochemical properties, the disease phenotype and the strain characteristics. The data derived from studies in vitro and from those using cell and animal models are compared with those obtained from the analyses of the naturally occurring disease. The extent of phenotypic variation in genetic prion disease is analyzed in comparison to that of the sporadic disease, which has recently been the topic of a systematic and detailed characterization.

The H187R mutation of the human prion protein induces conversion of recombinant prion protein to the PrP(Sc)-like form

Prion diseases are associated with a conformational switch in the prion protein (PrP) from its normal cellular form (denoted PrP(C)) to a disease-associated "scrapie" form (PrP(Sc)). A number of PrP(Sc)-like conformations can be generated by incubating recombinant PrP(C) at low pH, indicating that protonation of key residues is likely to destabilize PrP(C), facilitating its conversion to PrP(Sc). Here, we examine the stability of human PrP(C) with pH and find that PrP(C) fold stability is significantly reduced by the protonation of two histidine residues, His187 and His155. Mutation of His187 to an arginine, which imposes a permanently positively charged residue in this region of the protein, has a dramatic effect on the folding of PrP(C), resulting in a molecule that displays a markedly increased propensity to oligomerize. The oligomeric form is characterized by an increased β-sheet content, loss of fixed side chain interactions, and partial proteinase resistance. Hence, the protonation state of H187 appears to be crucial in determining the conformation of PrP; the unprotonated form favors native PrP(C), while the protonated form favors PrP(Sc)-like conformations. These results are relevant to the pathogenic H187R mutation found in humans, which is associated with an inherited prion disease [also termed Gerstmann-Sträussler-Scheinker (GSS) syndrome] with unusual features such as childhood neuropsychiatric illness. Our data imply that the intrinsic instability of the PrP(C) conformation in this variant is caused by a positive charge at this site in the protein. This mutation is distinct from all those associated with GSS, which have much more subtle physical consequences. The degree of instability might be the cause of the unusually early onset of mental disturbance in affected individuals.

The consequences of pathogenic mutations to the human prion protein

Prion diseases, in which the conformational transition of the native prion protein (PrP) to a misfolded form causes aggregation and subsequent neurodegeneration, have fascinated the scientific community as this transmissible disease appears to be purely protein-based. Disease can arise due to genetic factors only. At least 30 single point mutations have been indicated to cause disease in humans. Somehow, these mutations must influence the stability, processing and/or cellular interactions of PrP, such that aggregation can occur and disease develops. In this review, the current evidence for such effects of single point mutations is discussed, indicating that PrP can be affected in many different ways, although questions remain about the mechanism by which mutations cause disease.

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.

Electrostatic contribution to the thermodynamic and kinetic stability of the homotrimeric coiled coil Lpp-56: A computational study

The protein moiety of the Braun's E. coli outer membrane lipoprotein (Lpp-56) is an attractive object of biophysical investigation in several aspects. It is a homotrimeric, parallel coiled coil, a class of coiled coils whose stability and folding have been studied only occasionally. Lpp-56 possesses unique structural properties and exhibits extremely low rates of folding and unfolding. It is natural to ask how the specificity of the structure determines the extraordinary physical chemical properties of this protein. Recently, a seemingly controversial data on the stability and unfolding rate of Lpp-56 have been published (Dragan et al., Biochemistry 2004;43: 14891-14900; Bjelic et al., Biochemistry 2006;45:8931-8939). The unfolding rate constant measured using GdmCl as the denaturing agent, though extremely low, was substantially higher than that obtained on the basis of thermal unfolding. If this large difference arises from the effect of screening of electrostatic interactions induced by GdmCl, electrostatic interactions would appear to be an important factor determining the unusual properties of Lpp-56. We present here a computational analysis of the electrostatic properties of Lpp-56 combining molecular dynamics simulations and continuum pK calculations. The pH-dependence of the unfolding free energy is predicted in good agreement with the experimental data: the change in DeltaG between pH 3 and pH 7 is approximately 60 kJ mol(-1). The results suggest that the difference in the stability of the protein observed using different experimental methods is mainly because of the effect of the reduction of electrostatic interactions when the salt (GdmCl) concentration increases. We also find that the occupancy of the interhelical salt bridges is unusually high. We hypothesize that electrostatic interactions, and the interhelical salt bridges in particular, are an important factor determining the low unfolding rate of Lpp-56.

The stability and aggregation of ovine prion protein associated with classical and atypical scrapie correlates with the ease of unwinding of helix-2

Susceptibility to scrapie disease in sheep, the archetypal prion disease, correlates with polymorphisms within the ovine PrP (prion-related protein) gene. The VRQ (Val136Arg154Gln171) and AL141RQ (Ala136Leu141Arg154Gln171) allelic variants are associated with classical scrapie, whereas the ARR (Ala136Arg154Arg171), AF141RQ (Ala136Phe141Arg154Gln171) and AHQ (Ala136His154Gln171) allelic variants are associated with atypical scrapie. Recent studies have suggested that there are differences in the stability of PrPSc (abnormal disease-specific conformation of PrP) associated with these different forms of scrapie. To address which structural features of ovine PrP may contribute to this difference, in the present study we have investigated the conformational stability and susceptibility to aggregation of allelic variants of ovine PrP associated with classical or atypical scrapie. We find that the melting temperature of ovine recombinant VRQ and AL141RQ PrP is higher than that of AF141RQ, AHQ and ARR. In addition, monoclonal-antibody studies show that the region around helix-1 of VRQ and AL141RQ is less accessible compared with other ovine PrP allelic variants. Furthermore, the extent of both the structural change to copper-ion-treatment and denaturant-induced aggregation was reduced in PrP associated with atypical scrapie compared with PrP associated with classical scrapie. Through the use of molecular dynamics simulations we have found that these biochemical and biophysical properties of ovine PrP correlate with the ease of unwinding of helix-2 and a concurrent conformational change of the helix-2–helix-3 loop. These results reveal significant differences in the overall stability and potential for aggregation of different allelic variants of ovine PrP and consequently have implications for the differences in stability of PrPSc associated with classical and atypical scrapie.

Computational approaches to shed light on molecular mechanisms in biological processes

Computational approaches based on Molecular Dynamics simulations, Quantum Mechanical methods and 3D Quantitative Structure-Activity Relationships were employed by computational chemistry groups at the University of Milano-Bicocca to study biological processes at the molecular level. The paper reports the methodologies adopted and the results obtained on Aryl hydrocarbon Receptor and homologous PAS proteins mechanisms, the properties of prion protein peptides, the reaction pathway of hydrogenase and peroxidase enzymes and the defibrillogenic activity of tetracyclines.

Intracellular pH sensors: design principles and functional significance

Changes in intracellular pH regulate many cell behaviors, including proliferation, migration, and transformation. However, our understanding of how physiological changes in pH affect protein conformations and macromolecular assemblies is limited. We present design principles, current modeling predictions, and examples of pH sensors or proteins that have activities or ligand-binding affinities that are regulated by changes in intracellular pH.

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.

A Novel Member of the Protein Disulfide Oxidoreductase Family from Aeropyrum pernix K1: Structure, Function and Electrostatics

The formation of disulfide bonds between cysteine residues is a rate-limiting step in protein folding. To control this oxidative process, different organisms have developed different systems. In bacteria, disulfide bond formation is assisted by the Dsb protein family; in eukarya, disulfide bond formation and rearrangement are catalyzed by PDI. In thermophilic organisms, a potential key role in disulfide bond formation has recently been ascribed to a new cytosolic Protein Disulphide Oxidoreductase family whose members have a molecular mass of about 26 kDa and are characterized by two thioredoxin folds comprising a CXXC active site motif each. Here we report on the functional and structural characterization of ApPDO, a new member of this family, which was isolated from the archaeon Aeropyrum pernix K1. Functional studies have revealed that ApPDO can catalyze the reduction, oxidation and isomerization of disulfide bridges. Structural studies have shown that this protein has two CXXC active sites with fairly similar geometrical parameters typical of a stable conformation. Finally, a theoretical calculation of the cysteine pK(a) values has suggested that the two active sites have similar functional properties and each of them can impart activity to the enzyme. Our results are evidence of functional similarity between the members of the Protein Disulphide Oxidoreductase family and the eukaryotic enzyme PDI. However, as the different three-dimensional features of these two biological systems strongly suggest significantly different mechanisms of action, further experimental studies will be needed to make clear how different three-dimensional structures can result in systems with similar functional behavior.