Mapping the early steps in the pH-induced conformational conversion of the prion protein

Intrinsic determinants of prion protein neurotoxicity in Drosophila: from sequence to (dys)function

Prion diseases are fatal brain disorders characterized by deposition of insoluble isoforms of the prion protein (PrP). The normal and pathogenic structures of PrP are relatively well known after decades of studies. Yet our current understanding of the intrinsic determinants regulating PrP misfolding are largely missing. A 3D subdomain of PrP comprising the β2-α2 loop and helix 3 contains high sequence and structural variability among animals and has been proposed as a key domain regulating PrP misfolding. We combined in vivo work in Drosophila with molecular dynamics (MD) simulations, which provide additional insight to assess the impact of candidate substitutions in PrP from conformational dynamics. MD simulations revealed that in human PrP WT the β2-α2 loop explores multiple β-turn conformations, whereas the Y225A (rabbit PrP-like) substitution strongly favors a 310-turn conformation, a short right-handed helix. This shift in conformational diversity correlates with lower neurotoxicity in flies. We have identified additional conformational features and candidate amino acids regulating the high toxicity of human PrP and propose a new strategy for testing candidate modifiers first in MD simulations followed by functional experiments in flies. In this review we expand on these new results to provide additional insight into the structural and functional biology of PrP through the prism of the conformational dynamics of a 3D domain in the C-terminus. We propose that the conformational dynamics of this domain is a sensitive measure of the propensity of PrP to misfold and cause toxicity. This provides renewed opportunities to identify the intrinsic determinants of PrP misfolding through the contribution of key amino acids to different conformational states by MD simulations followed by experimental validation in transgenic flies.

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

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

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.

Molecular simulations of proteins: From simplified physical interactions to complex biological phenomena

Molecular dynamics simulation is the most popular computational technique for investigating the structural and dynamical behaviour of proteins, in search of the molecular basis of their function. Far from being a completely settled field of research, simulations are still evolving to best capture the essential features of the atomic interactions that govern a protein's inner motions. Modern force fields are becoming increasingly accurate in providing a physical description adequate to this purpose, and allow us to model complex biological systems under fairly realistic conditions. Furthermore, the use of accelerated sampling techniques is improving our access to the observation of progressively larger molecular structures, longer time scales, and more hidden functional events. In this review, the basic principles of molecular dynamics simulations and a number of key applications in the area of protein science are summarized, and some of the most important results are discussed. Examples include the study of the structure, dynamics and binding properties of 'difficult' targets, such as intrinsically disordered proteins and membrane receptors, and the investigation of challenging phenomena like hydration-driven processes and protein aggregation. The findings described provide an overall picture of the current state of this research field, and indicate new perspectives on the road ahead to the upcoming future of molecular simulations.

Identification of Two Early Folding Stage Prion Non-Local Contacts Suggested to Serve as Key Steps in Directing the Final Fold to Be Either Native or Pathogenic

The initial steps of the folding pathway of the C-terminal domain of the murine prion protein mPrP(90–231) are predicted based on the sequential collapse model (SCM). A non-local dominant contact is found to form between the connecting region between helix 1 and β-sheet 1 and the C-terminal region of helix 3. This non-local contact nucleates the most populated molten globule-like intermediate along the folding pathway. A less stable early non-local contact between segments 120–124 and 179–183, located in the middle of helix 2, promotes the formation of a less populated molten globule-like intermediate. The formation of the dominant non-local contact constitutes an example of the postulated Nature’s Shortcut to the prion protein collapse into the native structure. The possible role of the less populated molten globule-like intermediate is explored as the potential initiation point for the folding for three pathogenic mutants (T182A, I214V, and Q211P in mouse prion numbering) of the prion protein.

Mutation-Dependent Refolding of Prion Protein Unveils Amyloidogenic-Related Structural Ramifications: Insights from Molecular Dynamics Simulations

The main focus of prion structural biology studies is to understand the molecular basis of prion diseases caused by misfolding, and aggregation of the cellular prion protein PrP^C remains elusive. Several genetic mutations are linked with human prion diseases and driven by the conformational conversion of PrP^C to the toxic PrP^Sc. The main goal of this study is to gain a better insight into the molecular effect of disease-associated V210I mutation on this process by molecular dynamics simulations. This inherited mutation elicited copious structural changes in the β1-α1-β2 subdomain, including an unfolding of a helix α1 and the elongation of the β-sheet. These unusual structural changes likely appeared to detach the β1-α1-β2 subdomain from the α2-α3 core, an early misfolding event necessary for the conformational conversion of PrP^C to PrP^Sc. Ultimately, the unfolded α1 and its prior β1-α1 loop further engaged with unrestrained conformational dynamics and were widely considered as amyloidogenic-inducing traits. Furthermore, the resulting folding intermediate possesses a highly unstable β1-α1-β2 subdomain, thereby enhancing the aggregation of misfolded PrP^C through intermolecular interactions between frequently refolding regions. Briefly, these remarkable changes as seen in the mutant β1-α1-β2 subdomain are consistent with previous experimental results and thus provide a molecular basis of PrP^C misfolding associated with the conformational conversion of PrP^C to PrP^Sc.

A study on the secondary structure of the metalloregulatory protein CueR: effect of pH, metal ions and DNA

The response of CueR towards environmental changes in solution was investigated. CueR is a bacterial metal ion selective transcriptional metalloregulator protein, which controls the concentration of copper ions in the cell. Although several articles have been devoted to the discussion of the structural and functional features of this protein, CueR has not previously been extensively characterized in solution. Here, we studied the effect of change in pH, temperature, and the presence of specific or non-specific binding partners on the secondary structure of CueR with circular dichroism (CD) spectroscopy. A rather peculiar reversible pH-dependent secondary structure transformation was observed, elucidated and supplemented with pK_a estimation by PROPKA and CpHMD simulations suggesting an important role of His(76) and His(94) in this process. CD experiments revealed that the presence of DNA prevents this structural switch, suggesting that DNA locks CueR in the α-helical-rich form. In contrast to the non-cognate metal ions Hg^II, Cd^II and Zn^II, the presence of the cognate Ag^I ion affects the secondary structure of CueR, most probably by stabilizing the metal ion and DNA-binding domains of the protein.

Identification of critical amino acid residues in the regulatory N-terminal domain of PMEL

The pigment cell-specific protein PMEL forms a functional amyloid matrix in melanosomes onto which the pigment melanin is deposited. The amyloid core consists of a short proteolytic fragment, which we have termed the core-amyloid fragment (CAF) and perhaps additional parts of the protein, such as the PKD domain. A highly O-glycosylated repeat (RPT) domain also derived from PMEL proteolysis associates with the amyloid and is necessary to establish the sheet-like morphology of the assemblies. Excluded from the aggregate is the regulatory N-terminus, which nevertheless must be linked in cis to the CAF in order to drive amyloid formation. The domain is then likely cleaved away immediately before, during, or immediately after the incorporation of a new CAF subunit into the nascent amyloid. We had previously identified a 21 amino acid long region, which mediates the regulatory activity of the N-terminus towards the CAF. However, many mutations in the respective segment caused misfolding and/or blocked PMEL export from the endoplasmic reticulum, leaving their phenotype hard to interpret. Here, we employ a saturating mutagenesis approach targeting the motif at single amino acid resolution. Our results confirm the critical nature of the PMEL N-terminal region and identify several residues essential for PMEL amyloidogenesis.

Edge Strand Dissociation and Conformational Changes in Transthyretin under Amyloidogenic Conditions

During amyloidogenesis, proteins undergo conformational changes that allow them to aggregate and assemble into insoluble, fibrillar structures. Soluble oligomers that form during this process typically contain 2-24 monomeric subunits and are cytotoxic. Before the formation of these soluble oligomers, monomeric species first adopt aggregation-competent conformations. Knowledge of the structures of these intermediate states is invaluable to the development of molecular strategies to arrest pathological amyloid aggregation. However, the highly dynamic and interconverting nature of amyloidogenic species limits biophysical characterization of their structures during amyloidogenesis. Here, we use molecular dynamics simulations to probe conformations sampled by monomeric transthyretin under amyloidogenic conditions. We show that certain β-strands in transthyretin tend to unfold and sample nonnative conformations and that the edge strands in one β-sheet (the DAGH sheet) are particularly susceptible to conformational changes in the monomeric state. We also find that changes in the tertiary structure of transthyretin can be associated with disruptions to the secondary structure. We evaluated the conformations produced by molecular dynamics by calculating how well molecular-dynamics-derived structures reproduced NMR-derived interatomic distances. Finally, we leverage our computational results to produce experimentally testable hypotheses that may aid experimental explorations of pathological conformations of transthyretin.

Insight From Animals Resistant to Prion Diseases: Deciphering the Genotype - Morphotype - Phenotype Code for the Prion Protein

Prion diseases are a group of neurodegenerative diseases endemic in humans and several ruminants caused by the misfolding of native prion protein (PrP) into pathological conformations. Experimental work and the mad-cow epidemic of the 1980s exposed a wide spectrum of animal susceptibility to prion diseases, including a few highly resistant animals: horses, rabbits, pigs, and dogs/canids. The variable susceptibility to disease offers a unique opportunity to uncover the mechanisms governing PrP misfolding, neurotoxicity, and transmission. Previous work indicates that PrP-intrinsic differences (sequence) are the main contributors to disease susceptibility. Several residues have been cited as critical for encoding PrP conformational stability in prion-resistant animals, including D/E159 in dog, S167 in horse, and S174 in rabbit and pig PrP (all according to human numbering). These amino acids alter PrP properties in a variety of assays, but we still do not clearly understand the structural correlates of PrP toxicity. Additional insight can be extracted from comparative structural studies, followed by molecular dynamics simulations of selected mutations, and testing in manipulable animal models. Our working hypothesis is that protective amino acids generate more compact and stable structures in a C-terminal subdomain of the PrP globular domain. We will explore this idea in this review and identify subdomains within the globular domain that may hold the key to unravel how conformational stability and disease susceptibility are encoded in PrP.

Theoretical and computational advances in protein misfolding

Misfolded proteins escape the cellular quality control mechanism and fail to fold properly or remain correctly folded leading to a loss in their functional specificity. Thus misfolding of proteins cause a large number of very different diseases ranging from errors in metabolism to various types of complex neurodegenerative diseases. A theoretical and computational perspective of protein misfolding is presented with a special emphasis on its salient features, mechanism and consequences. These insights quantitatively analyze different determinants of misfolding, that may be applied to design disease specific molecular targets.

Drivers of α-Sheet Formation in Transthyretin under Amyloidogenic Conditions

Amyloid diseases make up a set of fatal disorders in which proteins aggregate to form fibrils that deposit in tissues throughout the body. Amyloid-associated diseases are challenging to study because amyloid formation occurs on time scales that span several orders of magnitude and involve heterogeneous, interconverting protein conformations. The development of more effective technologies to diagnose and treat amyloid disease requires both a map of the conformations sampled during amyloidogenesis and an understanding of the molecular mechanisms that drive this process. In prior molecular dynamics simulations of amyloid proteins, we observed the formation of a nonstandard type of secondary structure, called α-sheet, that we proposed is associated with the pathogenic conformers in amyloid disease, the soluble oligomers. However, the detailed molecular interactions that drive the conversion to α-sheet remain elusive. Here we use molecular dynamics simulations to interrogate a critical event in transthyretin aggregation, the formation of aggregation-competent, monomeric species. We show that conformational changes in one of the two β-sheets in transthyretin enable solvent molecules and polar side chains to form electrostatic interactions with main-chain peptide groups to facilitate and modulate conversion to α-sheet secondary structure. Our results shed light on the early conformational changes that drive transthyretin toward the α-sheet structure associated with toxicity. Delineation of the molecular events that lead to aggregation at atomic resolution can aid strategies to target the early, critical toxic soluble oligomers.

Strengths and Weaknesses of Molecular Simulations of Electrosprayed Droplets

The origin and the magnitude of the charge in a macroion are critical questions in mass spectrometry analysis coupled to electrospray and other ionization techniques that transfer analytes from the bulk solution into the gaseous phase via droplets. In many circumstances, it is the later stages of the existence of a macroion in the containing solvent drop before the detection that determines the final charge state. Experimental characterization of small (with linear dimensions of several nanometers) and short-lived droplets is quite challenging. Molecular simulations in principle may provide insight exactly in this challenging for experiments regime. We discuss the strengths and weaknesses of the molecular modeling of electrosprayed droplets using molecular dynamics. We illustrate the limitations of the molecular modeling in the analysis of large macroions and specifically proteins away from their native states. Graphical Abstract ᅟ.

An Evaluation of the Potential of NMR Spectroscopy and Computational Modelling Methods to Inform Biopharmaceutical Formulations

Protein-based therapeutics are considered to be one of the most important classes of pharmaceuticals on the market. The growing need to prolong stability of high protein concentrations in liquid form has proven to be challenging. Therefore, significant effort is being made to design formulations which can enable the storage of these highly concentrated protein therapies for up to 2 years. Currently, the excipient selection approach involves empirical high-throughput screening, but does not reveal details on aggregation mechanisms or the molecular-level effects of the formulations under storage conditions. Computational modelling approaches have the potential to elucidate such mechanisms, and rapidly screen in silico prior to experimental testing. Nuclear Magnetic Resonance (NMR) spectroscopy can also provide complementary insights into excipient–protein interactions. This review will highlight the underpinning principles of molecular modelling and NMR spectroscopy. It will also discuss the advancements in the applications of computational and NMR approaches in investigating excipient–protein interactions.

Effects of pH and aggregation in the human prion conversion into scrapie form: a study using molecular dynamics with excited normal modes

There are two different prion conformations: (1) the cellular natural (PrP^C) and (2) the scrapie (PrP^Sc), an infectious form that tends to aggregate under specific conditions. PrP^C and PrP^Sc are widely different regarding secondary and tertiary structures. PrP^Sc contains more and longer β-strands compared to PrP^C. The lack of solved PrP^Sc structures precludes a proper understanding of the mechanisms related to the transition between cellular and scrapie forms, as well as the aggregation process. In order to investigate the conformational transition between PrP^C and PrP^Sc, we applied MDeNM (molecular dynamics with excited normal modes), an enhanced sampling simulation technique that has been recently developed to probe large structural changes. These simulations yielded new structural rearrangements of the cellular prion that would have been difficult to obtain with standard MD simulations. We observed an increase in β-sheet formation under low pH (≤ 4) and upon oligomerization, whose relevance was discussed on the basis of the energy landscape theory for protein folding. The characterization of intermediate structures corresponding to transition states allowed us to propose a conversion model from the cellular to the scrapie prion, which possibly ignites the fibril formation. This model can assist the design of new drugs to prevent neurological disorders related to the prion aggregation mechanism.

Fatal familial insomnia and sporadic fatal insomnia

Fatal familial insomnia (FFI) and sporadic fatal insomnia (sFI), or thalamic form of sporadic Creutzfeldt-Jakob disease MM2 (sCJDMM2T), are prion diseases originally named and characterized in 1992 and 1999, respectively. FFI is genetically determined and linked to a D178N mutation coupled with the M129 genotype in the prion protein gene (PRNP) at chromosome 20. sFI is a phenocopy of FFI and likely its sporadic form. Both diseases are primarily characterized by progressive sleep impairment, disturbances of autonomic nervous system, and motor signs associated with severe loss of nerve cells in medial thalamic nuclei. Both diseases harbor an abnormal disease-associated prion protein isoform, resistant to proteases with relative mass of 19 kDa identified as resPrP^TSE type 2. To date at least 70 kindreds affected by FFI with 198 members and 18 unrelated carriers along with 25 typical cases of sFI have been published. The D178N-129M mutation is thought to cause FFI by destabilizing the mutated prion protein and facilitating its conversion to PrP^TSE. The thalamus is the brain region first affected. A similar mechanism triggered spontaneously may underlie sFI.

Sequence-dependent dynamical instability of the human prion protein: a comparative simulation study

The present study aimed to explore the most probable regions of the human prion protein backbone for which the initial steps of conformational transitions as a result of intrinsic and extrinsic perturbing factors on the protein structure can be assigned. A total of 0.3-μs molecular dynamics simulations on several analog structures of the protein have been performed. To mimic the impact of the extrinsic and intrinsic destructive parameters on the dynamical characteristics of the protein, mild acidic conditions and R208H mutation have been simulated. The findings indicated that distribution of conformational flexibilities along the protein chain was almost independent of the induced perturbing factors, and was mostly centralized on certain distinct parts of the structure comprising residues 132-145 and 187-203. Analyses also revealed that the segment comprising residues 187-203 may be considered as a peptide sequence, possessing high potential to start the initial steps of conformational rearrangements due to the induced physicochemical alterations. Sequence alignment and molecular dynamics data also revealed that segment 178-203 prefers to accommodate in extended structures rather than α-helices. Region 178-203 may be considered as a peptide switch capable of initiating the conformational transitions due to the introduced modifications and perturbing parameters.

Extensive Assessment of Various Computational Methods for Aspartate's pKa Shift

A series of computational methods for pK_a shift prediction are extensively tested on a set of benchmark protein systems, aiming at identifying pitfalls and evaluating their performance on high variants. Including 19 ASP residues in 10 protein systems, the benchmark set consists of both residues with highly shifted pK_a values as well as those varying little from the reference value, with an experimental RMS free energy differences of 2.49 kcal/mol with respect to blocked amino acid, namely the RMS pK_a shift being 1.82 pK_a units. The constant pH molecular dynamics (MD), alchemical methods, PROPKA3.1, and multiconformation continuum electrostatics give RMSDs of 1.52, 2.58, 1.37, and 3.52 pK_a units, respectively, on the benchmark set. The empirical scoring method is the most accurate one with extremely low computational cost, and the pH-dependent model is also able to provide accurate results, while the accuracy of MD sampling incorporating alchemical free energy simulation is prohibited by convergence achievement and the performance of conformational search incorporating multiconformation continuum electrostatics is bad. Former research works did not define statistical uncertainty with care and yielded the questionable conclusion that alchemical methods perform well in most benchmarks. In this work the traditional alchemical methods are thoroughly tested for high variants. We also performed the first application of nonequilibrium alchemical methods to the pK_a cases.

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

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

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

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

Melanosomal formation of PMEL core amyloid is driven by aromatic residues

PMEL is a pigment cell protein that forms physiological amyloid in melanosomes. Many amyloids and/or their oligomeric precursors are toxic, causing or contributing to severe, incurable diseases including Alzheimer’s and prion diseases. Striking similarities in intracellular formation pathways between PMEL and various pathological amyloids including Aβ and PrP^Sc suggest PMEL is an excellent model system to study endocytic amyloid. Learning how PMEL fibrils assemble without apparent toxicity may help developing novel therapies for amyloid diseases. Here we identify the critical PMEL domain that forms the melanosomal amyloid core (CAF). An unbiased alanine-scanning screen covering the entire region combined with quantitative electron microscopy analysis of the full set of mutants uncovers numerous essential residues. Many of these rely on aromaticity for function suggesting a role for π-stacking in melanosomal amyloid assembly. Various mutants are defective in amyloid nucleation. This extensive data set informs the first structural model of the CAF and provides insights into how the melanosomal amyloid core forms.

Insights from molecular dynamics simulations for computational protein design

A grand challenge in the field of structural biology is to design and engineer proteins that exhibit targeted functions. Although much success on this front has been achieved, design success rates remain low, an ever-present reminder of our limited understanding of the relationship between amino acid sequences and the structures they adopt. In addition to experimental techniques and rational design strategies, computational methods have been employed to aid in the design and engineering of proteins. Molecular dynamics (MD) is one such method that simulates the motions of proteins according to classical dynamics. Here, we review how insights into protein dynamics derived from MD simulations have influenced the design of proteins. One of the greatest strengths of MD is its capacity to reveal information beyond what is available in the static structures deposited in the Protein Data Bank. In this regard simulations can be used to directly guide protein design by providing atomistic details of the dynamic molecular interactions contributing to protein stability and function. MD simulations can also be used as a virtual screening tool to rank, select, identify, and assess potential designs. MD is uniquely poised to inform protein design efforts where the application requires realistic models of protein dynamics and atomic level descriptions of the relationship between dynamics and function. Here, we review cases where MD simulations was used to modulate protein stability and protein function by providing information regarding the conformation(s), conformational transitions, interactions, and dynamics that govern stability and function. In addition, we discuss cases where conformations from protein folding/unfolding simulations have been exploited for protein design, yielding novel outcomes that could not be obtained from static structures.

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.

Structural and functional diversity among amyloid proteins: Agents of disease, building blocks of biology, and implications for molecular engineering

Amyloids have long been associated with protein dysfunction and neurodegenerative diseases, but recent research has demonstrated that some organisms utilize the unique properties of the amyloid fold to create functional structures with important roles in biological processes. Additionally, new engineering approaches have taken advantage of amyloid structures for implementation in a wide variety of materials and devices. In this review, the role of amyloid in human disease is discussed and compared to the functional amyloids, which serve a largely structural purpose. We then consider the use of amyloid constructs in engineering applications, including their utility as building blocks for synthetic biology and molecular engineering. Biotechnol. Bioeng. 2017;114: 7-20. © 2016 Wiley Periodicals, Inc.

Free energy and hidden barriers of the β-sheet structure of prion protein

On-the-fly free-energy parametrization is a new collective variable biasing approach akin to metadynamics with one important distinction: rather than acquiring an accelerated distribution via a history-dependent bias potential, sampling on this distribution is achieved from the beginning of the simulation using temperature-accelerated molecular dynamics. In the present work, we compare the performance of both approaches to compute the free-energy profile along a scalar collective variable measuring the H-bond registry of the β-sheet structure of the mouse Prion protein. Both methods agree on the location of the free-energy minimum, but free-energy profiles from well-tempered metadynamics are subject to a much higher degree of statistical noise due to hidden barriers. The sensitivity of metadynamics to hidden barriers is shown to be a consequence of the history dependence of the bias potential, and we detail the nature of these barriers for the prion β-sheet. In contrast, on-the-fly parametrization is much less sensitive to these barriers and thus displays improved convergence behavior relative to that of metadynamics. While hidden barriers are a frequent and central issue in free-energy methods, on-the-fly free-energy parametrization appears to be a robust and preferable method to confront this issue.

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.

The molecular mechanism of conformational changes of the triplet prion fibrils for pH

The HET-s prion fibril, which is found in the filamentous fungus Podospora anserina, exhibits conformational changes due to variations in pH.

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

Nucleic acid induced unfolding of recombinant prion protein globular fragment is pH dependent

Nucleic acid can catalyze the conversion of α-helical cellular prion protein to β-sheet rich Proteinase K resistant prion protein oligomers and amyloid polymers in vitro and in solution. Because unfolding of a protein molecule from its ordered α-helical structure is considered to be a necessary step for the structural conversion to its β-sheet rich isoform, we have studied the unfolding of the α-helical globular 121-231 fragment of mouse recombinant prion protein in the presence of different nucleic acids at neutral and acid pH. Nucleic acids, either single or double stranded, do not have any significant effect on the secondary structure of the protein fragment at neutral pH; however the protein secondary structure is modified by the nucleic acids at pH 5. Nucleic acids do not show any significant effect on the temperature induced unfolding of the globular prion protein domain at neutral pH which, however, undergoes a gross conformational change at pH 5 as evidenced from the lowering of the midpoint of thermal denaturation temperatures, Tm, of the protein. The extent of Tm decrease shows a dependence on the nature of nucleic acid. The interaction of nucleic acid with the nonpolar groups exposed from the protein interior at pH 5 probably contributes substantially to the unfolding process of the protein.

Polymorphisms at amino acid residues 141 and 154 influence conformational variation in ovine PrP

Polymorphisms in ovine PrP at amino acid residues 141 and 154 are associated with susceptibility to ovine prion disease: Leu141Arg154 with classical scrapie and Phe141Arg154 and Leu141His154 with atypical scrapie. Classical scrapie is naturally transmissible between sheep, whereas this may not be the case with atypical scrapie. Critical amino acid residues will determine the range or stability of structural changes within the ovine prion protein or its functional interaction with potential cofactors, during conversion of PrPC to PrPSc in these different forms of scrapie disease. Here we computationally identified that regions of ovine PrP, including those near amino acid residues 141 and 154, displayed more conservation than expected based on local structural environment. Molecular dynamics simulations showed these conserved regions of ovine PrP displayed genotypic differences in conformational repertoire and amino acid side-chain interactions. Significantly, Leu141Arg154 PrP adopted an extended beta sheet arrangement in the N-terminal palindromic region more frequently than the Phe141Arg154 and Leu141His154 variants. We supported these computational observations experimentally using circular dichroism spectroscopy and immunobiochemical studies on ovine recombinant PrP. Collectively, our observations show amino acid residues 141 and 154 influence secondary structure and conformational change in ovine PrP that may correlate with different forms of scrapie.

Structural diversity and initial oligomerization of PrP106-126 studied by replica-exchange and conventional molecular dynamics simulations

Prion diseases are marked by cerebral accumulation of the abnormal isoform of the prion protein. A fragment of prion protein composed of residues 106–126 (PrP106–126) exhibits similar properties to full length prion and plays a key role in the conformational conversion from cellular prion to its pathogenic pattern. Soluble oligomers of PrP106–126 have been proposed to be responsible for neurotoxicity. However, the monomeric conformational space and initial oligomerization of PrP106–126 are still obscure, which are very important for understanding the conformational conversion of PrP106–126. In this study, replica exchange molecular dynamics simulations were performed to investigate monomeric and dimeric states of PrP106–126 in implicit solvent. The structural diversity of PrP106–126 was observed and this peptide did not acquire stable structure. The dimeric PrP106–126 also displayed structural diversity and hydrophobic interaction drove the dimerization. To further study initial oligomerization of PrP106–126, 1 µs conventional molecular dynamics simulations of trimer and tetramer formation were carried out in implicit solvent. We have observed the spontaneous formation of several basic oligomers and stable oligomers with high β-sheet contents were sampled in the simulations of trimer and tetramer formation. The β-hairpin formed in hydrophobic tail of PrP106–126 with residues 118–120 in turn may stabilize these oligomers and seed the formation oligomers. This study can provide insight into the detailed information about the structure of PrP106–126 and the dynamics of aggregation of monomeric PrP106–126 into oligomers in atomic level.

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.

Impact of methionine oxidation as an initial event on the pathway of human prion protein conversion

Prion diseases comprise a group of fatal neurodegenerative disorders characterized by the autocatalytic conversion of the cellular prion protein PrP(C) into the infectious misfolded isoform PrP(Sc). Increasing evidence supports a specific role of oxidative stress in the onset of pathogenesis. Although the associated molecular mechanisms remain to be elucidated in detail, several studies currently suggest that methionine oxidation already detected in misfolded PrP(Sc) destabilizes the native PrP fold as an early event in the conversion pathway. To obtain more insights about the specific impact of surface-exposed methionine residues on the oxidative-induced conversion of human PrP we designed, produced, and comparatively investigated two new pseudosulfoxidation mutants of human PrP 121-231 that comprises the well-folded C-terminal domain. Applying circular dichroism spectroscopy and dynamic light scattering techniques we showed that pseudosulfoxidation of all surface exposed Met residues formed a monomeric molten globule-like species with striking similarities to misfolding intermediates recently reported by other groups. However, individual pseudosulfoxidation at the polymorphic M129 site did not significantly contribute to the structural destabilization. Further metal-induced oxidation of the partly unfolded pseudosulfoxidation mutant resulted in the formation of an oligomeric state that shares a comparable size and stability with PrP oligomers detected after the application of different other triggers for structural conversion, indicating a generic misfolding pathway of PrP. The obtained results highlight the specific importance of methionine oxidation at surface exposed residues for PrP misfolding, strongly supporting the hypothesis that increased oxidative stress could be one causative event for sporadic prion diseases and other neurodegenerative disorders.

Probing early misfolding events in prion protein mutants by NMR spectroscopy

The post-translational conversion of the ubiquitously expressed cellular form of the prion protein, PrPC, into its misfolded and pathogenic isoform, known as prion or PrPSc, plays a key role in prion diseases. These maladies are denoted transmissible spongiform encephalopathies (TSEs) and affect both humans and animals. A prerequisite for understanding TSEs is unraveling the molecular mechanism leading to the conversion process whereby most α-helical motifs are replaced by β-sheet secondary structures. Importantly, most point mutations linked to inherited prion diseases are clustered in the C-terminal domain region of PrPC and cause spontaneous conversion to PrPSc. Structural studies with PrP variants promise new clues regarding the proposed conversion mechanism and may help identify "hot spots" in PrPC involved in the pathogenic conversion. These investigations may also shed light on the early structural rearrangements occurring in some PrPC epitopes thought to be involved in modulating prion susceptibility. Here we present a detailed overview of our solution-state NMR studies on human prion protein carrying different pathological point mutations and the implications that such findings may have for the future of prion research.