Structural and hydration properties of the partially unfolded states of the prion protein

Flexibility and binding affinity in protein-ligand, protein-protein and multi-component protein interactions: limitations of current computational approaches

The recognition process between a protein and a partner represents a significant theoretical challenge. In silico structure-based drug design carried out with nothing more than the three-dimensional structure of the protein has led to the introduction of many compounds into clinical trials and numerous drug approvals. Central to guiding the discovery process is to recognize active among non-active compounds. While large-scale computer simulations of compounds taken from a library (virtual screening) or designed de novo are highly desirable in the post-genomic area, many technical problems remain to be adequately addressed. This article presents an overview and discusses the limits of current computational methods for predicting the correct binding pose and accurate binding affinity. It also presents the performances of the most popular algorithms for exploring binary and multi-body protein interactions.

Intrinsic determinants of Aβ(12-24) pH-dependent self-assembly revealed by combined computational and experimental studies

The propensity of amyloid- (A) peptide to self-assemble into highly ordered amyloid structures lies at the core of their accumulation in the brain during Alzheimer's disease. By using all-atom explicit solvent replica exchange molecular dynamics simulations, we elucidated at the atomic level the intrinsic determinants of the pH-dependent dimerization of the central hydrophobic segment A and related these with the propensity to form amyloid fibrils measured by experimental tools such as atomic force microscopy and fluorescence. The process of A dimerization was evaluated in terms of free energy landscape, side-chain two-dimensional contact probability maps, -sheet registries, potential mean force as a function of inter-chain distances, secondary structure development and radial solvation distributions. We showed that dimerization is a key event in A amyloid formation; it is highly prompted in the order of pH 5.02.98.4 and determines further amyloid growth. The dimerization is governed by a dynamic interplay of hydrophobic, electrostatic and solvation interactions permitting some variability of -sheets at each pH. These results provide atomistic insight into the complex process of molecular recognition detrimental for amyloid growth and pave the way for better understanding of the molecular basis of amyloid diseases.

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.

PrP assemblies: spotting the responsible regions in prion propagation

The "protein only" hypothesis states that the key phenomenon in prion pathogenesis is the conversion of the host protein (PrPC) into a b-sheet enriched polymeric and pathogenic conformer (PrPSc). However the region of PrP bearing the information for structural transfer is still controversial. In a recent report, we highlighted the role of the C terminal part i.e. the helixes H2 and H3, using mutation approaches on recombinant PrP. The H2H3 was shown to be the minimal region necessary to reproduce the oligomerisation pattern of the full-length protein. The oligomers produced from isolated H2H3 domain presented the same structural characteristics as the oligomers formed from the full-length PrP. Combining other groups' results, this paper further discusses the relative, direct or indirect role of different PrP regions in assembly. The H2H3 region represents the core of PrP oligomers and fibrils, whereas the N terminus could explain divergences among different aggregates. Finally this review evocates the possibility to separate the domain involved in prion information transference (i.e. prion replication) from the domain bearing the cytotoxicity properties.

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.

Highly polar environments catalyze the unfolding of PrP C helix 1

The first alpha-helix (H1) likely plays an important role in the conversion of the cellular prion protein (PrP(C)) into its pathogenic isoform (PrP(Sc)). In this conversion, H1 may either have to unfold or may represent a site of intermolecular contact. A recent molecular dynamics simulation suggested that H1 can unfold if it is detached from the protein core (Hirschberger et al. in Biophys J 90:3908, 2006). It has been hypothesized that the high dielectric constant epsilon (S) of the bulk water environment facilitates the unfolding of H1. To check this hypothesis, we performed a number of replica exchange molecular dynamics simulations of an H1 peptide in solvents of different epsilon (S). We found that the equilibrium helix fraction in water is less than 40%, in agreement with previous experimental findings, and that the helix unfolds much faster in water than in less polar solvents. The kinetically stabilizing effect of the organic solvents is largely unspecific and correlates well with their dielectric constant epsilon (S).

Association thermodynamics and conformational stability of beta-sheet amyloid beta(17-42) oligomers: effects of E22Q (Dutch) mutation and charge neutralization

Amyloid fibrils are associated with many neurodegenerative diseases. It was found that amyloidogenic oligomers, not mature fibrils, are neurotoxic agents related to these diseases. Molecular mechanisms of infectivity, pathways of aggregation, and molecular structure of these oligomers remain elusive. Here, we use all-atom molecular dynamics, molecular mechanics combined with solvation analysis by statistical-mechanical, three-dimensional molecular theory of solvation (also known as 3D-RISM-KH) in a new MM-3D-RISM-KH method to study conformational stability, and association thermodynamics of small wild-type Abeta(17-42) oligomers with different protonation states of Glu(22), as well the E22Q (Dutch) mutants. The association free energy of small beta-sheet oligomers shows near-linear trend with the dimers being thermodynamically more stable relative to the larger constructs. The linear (within statistical uncertainty) dependence of the association free energy on complex size is a consequence of the unilateral stacking of monomers in the beta-sheet oligomers. The charge reduction of the wild-type Abeta(17-42) oligomers upon protonation of the solvent-exposed Glu(22) at acidic conditions results in lowering the association free energy compared to the wild-type oligomers at neutral pH and the E22Q mutants. The neutralization of the peptides because of the E22Q mutation only marginally affects the association free energy, with the reduction of the direct electrostatic interactions mostly compensated by the unfavorable electrostatic solvation effects. For the wild-type oligomers at acidic conditions such compensation is not complete, and the electrostatic interactions, along with the gas-phase nonpolar energetic and the overall entropic effects, contribute to the lowering of the association free energy. The differences in the association thermodynamics between the wild-type Abeta(17-42) oligomers at neutral pH and the Dutch mutants, on the one hand, and the Abeta(17-42) oligomers with protonated Glu(22), on the other, may be explained by destabilization of the inter- and intrapeptide salt bridges between Asp(23) and Lys(28). Peculiarities in the conformational stability and the association thermodynamics for the different models of the Abeta(17-42) oligomers are rationalized based on the analysis of the local physical interactions and the microscopic solvation structure.

Characterizing the denatured state of human prion 121-230

Misfolding and aggregation of the prion protein (PrP) are responsible for the development of fatal transmissible neurodegenerative diseases. PrP undergoes structural conversion from a natively folded state into a misfolded state, resulting in insoluble amyloid fibrils. Partial unfolding has been recognized as an essential step in fibrillation. The strong correlation of unfolding and fibrillation emphasizes the importance of denatured states. To gain insight into possible aggregation-prone denatured states, we characterized the denatured state of human prion (huPrP) 121-230 near extended conformation by self-guided Langevin dynamics simulations. Our results revealed that denatured huPrP is partially folded with alpha-helical structure.

The oligomerization properties of prion protein are restricted to the H2H3 domain

The propensity of the prion protein (PrP) to adopt different structures is a clue to its pathological behavior. The determination of the region involved in the PrP(C) to PrP(Sc) conversion is fundamental for the understanding of the mechanisms underlying this process at the molecular level. In this paper, the polymerization of the helical H2H3 domain of ovine PrP (OvPrP) was compared to the full-length construct (using chromatography and light scattering). We show that the oligomerization patterns are identical, although the H2H3 domain has a higher polymerization rate. Furthermore, the depolymerization kinetics of purified H2H3 oligomers compared to those purified from the full-length PrP reveal that regions outside H2H3 do not significantly contribute to the oligomerization process. By combining rational mutagenesis and molecular dynamics to investigate the early stages of H2H3 oligomerization, we observe a conformationally stable beta-sheet structure that we propose as a possible nucleus for oligomerization; we also show that single point mutations in H2 and H3 present structural polymorphisms and oligomerization properties that could constitute the basis of species or strain variability.

Insights into prion protein function from atomistic simulations

Computer simulations are a powerful tool for studies of biological systems. They have often been used to study prion protein (PrP), a protein responsible for neurodegenerative diseases, which include "mad cow disease" in cattle and Creutzfeldt-Jacob disease in humans. An important aspect of the prion protein is its interaction with copper ion, which is thought to be relevant for PrP's yet undetermined function and also potentially play a role in prion diseases. for studies of copper attachment to the prion protein, computer simulations have often been used to complement experimental data and to obtain binding structures of Cu-PrP complexes. This paper summarizes the results of recent ab initio calculations of copper-prion protein interactions focusing on the recently discovered concentration-dependent binding modes in the octarepeat region of this protein. In addition to determining the binding structures, computer simulations were also used to make predictions about PrP's function and the role of copper in prion diseases. The results demonstrate the predictive power and applicability of ab initio simulations for studies of metal-biomolecular complexes.

The unfolding of the prion protein sheds light on the mechanisms of prion susceptibility and species barrier

Prion diseases are a group of fatal neurodegenerative disorders that manifest as infectious, sporadic, or familial and are all associated with the misfolding of the prion protein (PrP). Disease-modulating polymorphisms in the PrP amino acid sequence can make an individual more or less susceptible to infection. One example is the presence of arginine in place of glutamine at position 171 in sheep, which confers resistance to scrapie. To investigate whether the physical folding properties of PrP are influenced by the presence of arginine at codon 171, we have introduced the mutation at the equivalent position (codon 167) in recombinant mouse PrP. We have then compared the unfolding properties of wild-type PrP and the Q167R mutant by monitoring the fluorescence and circular dichroism of folding-sensitive tryptophan mutants. For both wild-type PrP and the Q167R mutant the formation of secondary structure and tertiary structure is concurrent, which indicates that unfolding proceeds without the accumulation of an equilibrium intermediate. The major effect of the mutation is the destabilization of the protein as shown by the shift of the unfolding transition, which can be rationalized from high-resolution structures of PrP. Comparison of the unfolding pathways of mouse and hamster PrP highlights dramatic differences in the mechanisms of folding, which may contribute to the species barrier effect that is observed in the transmission of prion disease.

Hydration profiles of amyloidogenic molecular structures

Hydration shells of normal proteins display regions of highly structured water as well as patches of less structured bulk-like water. Recent studies suggest that isomers with larger surface densities of patches of bulk-like water have an increased propensity to aggregate. These aggregates are toxic to the cellular environment. Hence, the early detection of these toxic deposits is of paramount medical importance. We show that various morphological states of association of such isomers can be differentiated from the normal protein background based on the characteristic partition between bulk, caged, and surface hydration water and the magnetic resonance (MR) signals of this water. We derive simple mathematical equations relating the compartmentalization of water to the local hydration fraction and the packing density of the newly formed molecular assemblies. Then, we employ these equations to predict the MR response of water constrained by protein aggregation. Our results indicate that single units and compact aggregates that contain no water between constituents induce a shift of the MR signal from normal protein background to values in the hyperintensity domain (bright spots), corresponding to bulk water. In contrast, large plaques that cage significant amounts of water between constituents are likely to generate MR responses in the hypointensity domain (dark spots), typical for strongly correlated water. The implication of these results is that amyloids can display both dark and bright spots when compared to the normal gray background tissue on MR images. In addition, our findings predict that the bright spots are more likely to correspond to amyloids in their early stage of development. The results help explain the MR contrast patterns of amyloids and suggest a new approach for identifying unusual protein aggregation related to disease.

Stability of single sheet GNNQQNY aggregates analyzed by replica exchange molecular dynamics: antiparallel versus parallel association

Protein and peptide aggregation into amyloid plaques is associated with a large variety of neurodegenerative diseases. The definition of the molecular bases of these pathologies is hampered by the transient nature of pre-fibrillar small-oligomers that are considered the toxic species. The ability of the peptide GNNQQNY to form amyloid-like structures makes it a good model to investigate the complex processes involved into amyloid fiber formation. By employing full atomistic replica exchange molecular dynamics simulations, we constructed the free energy surface of small assemblies of GNNQQNY to gain novel insights into the fiber formation process. The calculations suggest that the peptide exhibits a remarkable tendency to form both parallel and antiparallel beta-sheets. The data show that GNNQQNY preference for parallel or antiparallel beta-sheets is governed by a subtle balance of factors including assemblies' size, sidechain-sidechain interactions and pH. The samplings analysis provides a rationale to the observed trends.

Hydration effects on the HET-s prion and amyloid-beta fibrillous aggregates, studied with three-dimensional molecular theory of solvation

We study the thermodynamic properties of the experimental fragments of the amyloid fibril made of the HET-s prion proteins (the infectious element of the filamentous fungus Podospora anserina) and of amyloid-beta proteins (the major component of Alzheimer's disease-associated plaques) by using the three-dimensional molecular theory of solvation. The full quantitative picture of hydration effects, including the hydration thermodynamics and hydration structure around the fragments, is presented. For both the complexes, the hydration entropic effects dominate, which results in the entropic part offsetting the unfavorable energetic part of the free energy change upon the association. This is in accord with the fact that the hydrophobic cooperativity plays an essential role in the formation of amyloid fibrils. By calculating the partial molar volume of the proteins, we found that the volume change upon the association in both the systems is large and positive, with the implication that high pressure causes destabilization of the fibril. This observation is in good agreement with the recent experimental results. We also found that both the HET-s and amyloid-beta pentamers have loose intermolecular packing with voids. The three-dimensional molecular theory of solvation predicts that water molecules can be locked in the interior cavities along the fibril axis for both the HET-s and amyloid-beta proteins. We provide a detailed molecular picture of the structural water localized in the interior of the fibrils. Our results suggest that the interior hydration plays an important role in the structural stability of fibrils.

Expansion of the octarepeat domain alters the misfolding pathway but not the folding pathway of the prion protein

A misfolded conformation of the prion protein (PrP), PrP (Sc), is the essential component of prions, the infectious agents that cause transmissible neurodegenerative diseases. Insertional mutations that lead to an increase in the number of octarepeats (ORs) in PrP are linked to familial human prion disease. In this study, we investigated how expansion of the OR domain causes PrP to favor a prion-like conformation. Therefore, we compared the conformational and aggregation modulating properties of wild-type versus expanded OR domains, either as a fusion construct with the protein G B1 domain (GB1-OR) or as an integral part of full-length mouse PrP (MoPrP). Using circular dichroism spectroscopy, we first demonstrated that ORs are not unfolded but exist as an ensemble of three distinct conformers: polyproline helix-like, beta-turn, and "Trp-related". Domain expansion had little effect on the conformation of GB1-OR fusion proteins. When part of MoPrP however, OR domain expansion changed PrP's folding landscape, not by hampering the production of native alpha-helical monomers but by greatly reducing the propensity to form amyloid and by altering the assembly of misfolded, beta-rich aggregates. These features may relate to subtle pH-dependent conformational differences between wild-type and mutant monomers. In conclusion, we propose that PrP insertional mutations are pathogenic because they enhance specific misfolding pathways of PrP rather than by undermining native folding. This idea was supported by a trial bioassay in transgenic mice overexpressing wild-type MoPrP, where intracerebral injection of recombinant MoPrP with an expanded OR domain but not wild-type MoPrP caused prion disease.

Insights into stability and toxicity of amyloid-like oligomers by replica exchange molecular dynamics analyses

Deposition of insoluble amyloid plaques is frequently associated with a large variety of neurodegenerative diseases. However, data collected in the last decade have suggested that the neurotoxic action is exerted by prefibrillar, soluble assemblies of amyloid-forming proteins and peptides. The scarcity of structural data available for both amyloid-like fibrils and soluble oligomers is a major limitation for the definition of the molecular mechanisms linked to the onset of these diseases. Recently, the structural characterization of GNNQQNY and other peptides has shown a general feature of amyloid-like fibers, the so-called steric zipper motif. However, very little is known still about the prefibrillar oligomeric forms. By using replica exchange molecular dynamics we carried out extensive analyses of the properties of several small and medium GNNQQNY aggregates arranged through the steric zipper motif. Our data show that the assembly formed by two sheets, each made of two strands, arranged as in the crystalline states are highly unstable. Conformational free energy surfaces indicate that the instability of the model can be ascribed to the high reactivity of edge backbone hydrogen bonding donors/acceptors. On the other hand, data on larger models show that steric zipper interactions may keep small oligomeric forms in a stable state. These models simultaneously display two peculiar structural motifs: a tightly packed steric zipper interface and a large number of potentially reactive exposed strands. The presence of highly reactive groups on these assemblies likely generates two distinct evolutions. On one side the reactive groups quickly lead, through self-association, to the formation of ordered fibrils, on the other they may interfere with several cellular components thereby generating toxic effects. In this scenario, fiber formation propensity and toxicity of oligomeric states are two different manifestations of the same property: the hyper-reactivity of the exposed strands.

The beta-strand-loop-beta-strand conformation is marginally populated in beta2-microglobulin (20-41) peptide in solution as revealed by replica exchange molecular dynamics simulations

Solid-state NMR study shows that the 22-residue K3 peptide (Ser(20)-Lys(41)) from beta(2)-microglobulin (beta(2)m) adopts a beta-strand-loop-beta-strand conformation in its fibril state. Residue Pro(32) has a trans conformation in the fibril state of the peptide, while it adopts a cis conformation in the native state of full-length beta(2)m. To get insights into the structural properties of the K3 peptide, and determine whether the strand-loop-strand conformation is encoded at the monomeric level, we run all-atom explicit solvent replica exchange molecular dynamics on both the cis and trans variants. Our simulations show that the conformational space of the trans- and cis-K3 peptides is very different, with 1% of the sampled conformations in common at room temperature. In addition, both variants display only 0.3-0.5% of the conformations with beta-strand-loop-beta-strand character. This finding, compared to results on the Alzheimer's Abeta peptide, suggests that the biases toward aggregation leading to the beta-strand-loop-beta-strand conformation in fibrils are peptide-dependent.