Checking the pH-induced conformational transition of prion protein by molecular dynamics simulations: effect of protonation of histidine residues

RedOx regulation of LRRK2 kinase activity by active site cysteines

Mutations of the human leucine-rich repeat kinase 2 (LRRK2) have been associated with both, idiopathic and familial Parkinson's disease (PD). Most of these pathogenic mutations are located in the kinase domain (KD) or GTPase domain of LRRK2. In this study we describe a mechanism in which protein kinase activity can be modulated by reversible oxidation or reduction, involving a unique pair of adjacent cysteines, the "CC" motif. Among all human protein kinases, only LRRK2 contains this "CC" motif (C2024 and C2025) in the Activation Segment (AS) of the kinase domain. In an approach combining site-directed mutagenesis, biochemical analyses, cell-based assays, and Gaussian accelerated Molecular Dynamics (GaMD) simulations we could attribute a role for each of those cysteines. We employed reducing and oxidizing agents with potential clinical relevance to investigate effects on kinase activity and microtubule docking. We find that each cysteine gives a distinct contribution: the first cysteine, C2024, is essential for LRRK2 protein kinase activity, while the adjacent cysteine, C2025, contributes significantly to redox sensitivity. Implementing thiolates (R-S-) in GaMD simulations allowed us to analyse how each of the cysteines in the "CC" motif interacts with its surrounding residues depending on its oxidation state. From our studies we conclude that oxidizing agents can downregulate kinase activity of hyperactive LRRK2 PD mutations and may provide promising tools for therapeutic strategies.

Probing the role of the protonation state of a minor groove-linker histidine in Exd-Hox-DNA binding

DNA recognition and targeting by transcription factors (TFs) through specific binding are fundamental in biological processes. Furthermore, the histidine protonation state at the TF-DNA binding interface can significantly influence the binding mechanism of TF-DNA complexes. Nevertheless, the role of histidine in TF-DNA complexes remains underexplored. Here, we employed all-atom molecular dynamics simulations using AlphaFold2-modeled complexes based on previously solved co-crystal structures to probe the role of the His-12 residue in the Extradenticle (Exd)-Sex combs reduced (Scr)-DNA complex when binding to Scr and Ultrabithorax (Ubx) target sites. Our results demonstrate that the protonation state of histidine notably affected the DNA minor-groove width profile and binding free energy. Examining flanking sequences of various binding affinities derived from SELEX-seq experiments, we analyzed the relationship between binding affinity and specificity. We uncovered how histidine protonation leads to increased binding affinity but can lower specificity. Our findings provide new mechanistic insights into the role of histidine in modulating TF-DNA binding.

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.

Destabilization of polar interactions in the prion protein triggers misfolding and oligomerization

The prion protein (PrP) misfolds and oligomerizes at pH 4 in the presence of physiological salt concentrations. Low pH and salt cause structural perturbations in the monomeric prion protein that lead to misfolding and oligomerization. However, the changes in stability within different regions of the PrP prior to oligomerization are poorly understood. In this study, we have characterized the local stability in PrP at high resolution using amide temperature coefficients (T_C ) measured by nuclear magnetic resonance (NMR) spectroscopy. The local stability of PrP was investigated under native as well as oligomerizing conditions. We have also studied the rapidly oligomerizing PrP variant (Q216R) and the protective PrP variant (A6). We report that at low pH, salt destabilizes PrP at several polar residues, and the hydrogen bonds in helices α2 and α3 are weakened. In addition, salt changes the curvature of the α3 helix, which likely disrupts α2-α3 contacts and leads to oligomerization. These results are corroborated by the T_C values of rapidly oligomerizing Q216R-PrP. The poly-alanine substitution in A6-PrP stabilizes α2, which prevents oligomerization. Altogether, these results highlight the importance of native polar interactions in determining the stability of PrP and reveal the structural disruptions in PrP that lead to misfolding and oligomerization.

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.

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.

Unraveling the Molecular Mechanism of Prion H2 C-Terminus Misfolding by Metadynamics Simulations

Conformational transition from the normal cellular form of prion protein (PrP^C) to the pathogenic "scrapie" form (PrP^Sc) is considered to be a key event in the occurrence of prion disease. Additionally, the H2 C-terminus is widely considered to be a vital site for PrP conformational transition, which can be used as an important region to explore the potential mechanism of PrP misfolding. Therefore, to study the misfolding mechanism of PrP, 500 ns well-tempered metadynamics simulations were performed by focusing on the H2 C-terminus of PrP. For comparison, three systems were designed in total, including PrP in neutral and acidic conditions, as well as H187R mutant. The resulting free energy surfaces (FESs) obtained from metadynamics simulations reveal that acidic conditions and H187R mutation can facilitate PrP misfolding by decreasing free energy barriers for conformational transition and forming energy stable conformational states. Further analyses aimed at H2 C-terminus show that due to the increase of positive charge on residue 187 in both acidic and H187R systems, the electrostatic repulsion of residue 187 and R136/R156 increases greatly, which disrupts the electrostatic interaction network around H2 C-terminus and exposes the hydrophobic core to the solvent. Taken together, acidic conditions and H187R mutation can accelerate PrP misfolding mainly by forming more energetically stable metastable conformations with lower free energy barriers, and electrostatic network disruption involving residue 187 drives the initial misfolding of H2 C-terminus. This study provides quantitative insight into the related function of the H2 C-terminus in the PrP misfolding process, which may guide H2 C-terminus mediated drug design in the future.

Conformation and Domain Movement Analysis of Human Matrix Metalloproteinase-2: Role of Associated Zn2+ and Ca2+ Ions

Matrix metaloproteinase-2 (MMP-2) is an extracellular Zn²⁺ protease specific to type I and IV collagens. Its expression is associated with several inflammatory, degenerative, and malignant diseases. Conformational properties, domain movements, and interactions between MMP-2 and its associated metal ions were characterized using a 1.0 µs molecular dynamics simulation. Dihedral principle component analysis revealed ten families of conformations with the greatest degree of variability occurring in the link region connecting the catalytic and hemopexin domains. Dynamic cross-correlation analysis indicated domain movements corresponding to the opening and closing of the hemopexin domain in relation to the fibronectin and catalytic domains facilitated by the link region. Interaction energies were calculated using the molecular mechanics Poisson Boltzman surface area-interaction entropy (MMPBSA-IE) analysis method and revealed strong binding energies for the catalytic Zn²⁺ ion 1, Ca²⁺ ion 1, and Ca²⁺ ion 3 with significant conformational stability at the binding sites of Zn²⁺ ion 1 and Ca²⁺ ion 1. Ca²⁺ ion 2 diffuses freely away from its crystallographically defined binding site. Zn²⁺ ion 2 plays a minor role in conformational stability of the catalytic domain while Ca²⁺ ion 3 is strongly attracted to the highly electronegative sidechains of the Asp residues around the central β-sheet core of the hemopexin domain; however, the interacting residue sidechain carboxyl groups are outside of Ca²⁺ ion 3′s coordination sphere.

The protonation state of an evolutionarily conserved histidine modulates domainswapping stability of FoxP1

Forkhead box P (FoxP) proteins are members of the versatile Fox transcription factors, which control the timing and expression of multiple genes for eukaryotic cell homeostasis. Compared to other Fox proteins, they can form domain-swapped dimers through their DNA-binding –forkhead– domains, enabling spatial reorganization of distant chromosome elements by tethering two DNA molecules together. Yet, domain swapping stability and DNA binding affinity varies between different FoxP proteins. Experimental evidence suggests that the protonation state of a histidine residue conserved in all Fox proteins is responsible for pH-dependent modulation of these interactions. Here, we explore the consequences of the protonation state of another histidine (H59), only conserved within FoxM/O/P subfamilies, on folding and dimerization of the forkhead domain of human FoxP1. Dimer dissociation kinetics and equilibrium unfolding experiments demonstrate that protonation of H59 leads to destabilization of the domain-swapped dimer due to an increase in free energy difference between the monomeric and transition states. This pH–dependence is abolished when H59 is mutated to alanine. Furthermore, anisotropy measurements and molecular dynamics evidence that H59 has a direct impact in the local stability of helix H3. Altogether, our results highlight the relevance of H59 in domain swapping and folding stability of FoxP1.

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.

Interplay of buried histidine protonation and protein stability in prion misfolding

Misofolding of mammalian prion proteins (PrP) is believed to be the cause of a group of rare and fatal neurodegenerative diseases. Despite intense scrutiny however, the mechanism of the misfolding reaction remains unclear. We perform nuclear Magnetic Resonance and thermodynamic stability measurements on the C-terminal domains (residues 90–231) of two PrP variants exhibiting different pH-induced susceptibilities to aggregation: the susceptible hamster prion (GHaPrP) and its less susceptible rabbit homolog (RaPrP). The pKa of histidines in these domains are determined from titration experiments, and proton-exchange rates are measured at pH 5 and pH 7. A single buried highly conserved histidine, H187/H186 in GHaPrP/RaPrP, exhibited a markedly down shifted pKa ~5 for both proteins. However, noticeably larger pH-induced shifts in exchange rates occur for GHaPrP versus RaPrP. Analysis of the data indicates that protonation of the buried histidine destabilizes both PrP variants, but produces a more drastic effect in the less stable GHaPrP. This interpretation is supported by urea denaturation experiments performed on both PrP variants at neutral and low pH, and correlates with the difference in disease susceptibility of the two species, as expected from the documented linkage between destabilization of the folded state and formation of misfolded and aggregated species.

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.

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.

Partially Unfolded Forms of the Prion Protein Populated under Misfolding-promoting Conditions: CHARACTERIZATION BY HYDROGEN EXCHANGE MASS SPECTROMETRY AND NMR

The susceptibility of the cellular prion protein (PrP(C)) to convert to an alternative misfolded conformation (PrP(Sc)), which is the key event in the pathogenesis of prion diseases, is indicative of a conformationally flexible native (N) state. In the present study, hydrogen-deuterium exchange (HDX) in conjunction with mass spectrometry and nuclear magnetic resonance spectroscopy were used for the structural and energetic characterization of the N state of the full-length mouse prion protein, moPrP(23-231), under conditions that favor misfolding. The kinetics of HDX of 34 backbone amide hydrogens in the N state were determined at pH 4. In contrast to the results of previous HDX studies on the human and Syrian hamster prion proteins at a higher pH, various segments of moPrP were found to undergo different extents of subglobal unfolding events at pH 4, a pH at which the protein is known to be primed to misfold to a β-rich conformation. No residual structure around the disulfide bond was observed for the unfolded state at pH 4. The N state of the prion protein was observed to be at equilibrium with at least two partially unfolded forms (PUFs). These PUFs, which are accessed by stochastic fluctuations of the N state, have altered surface area exposure relative to the N state. One of these PUFs resembles a conformation previously implicated to be an initial intermediate in the conversion of monomeric protein into misfolded oligomer at pH 4.

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.

Short-time dynamics of pH-dependent conformation and substrate binding in the active site of beta-glucosidases: A computational study

The complete degradation of cellulose to glucose is essential to carbon turnover in terrestrial ecosystems and to engineered biofuel production. A rate-limiting step in this pathway is catalyzed by beta-glucosidase (BG) enzymes, which convert cellulobiose into two glucose molecules. The activity of these enzymes has been shown to vary with solution pH. However, it is not well understood how pH influences the enzyme conformation required for catalytic action on the substrate. A structural understanding of this pH effect is important for predicting shifts in BG activity in bioreactors and environmental matrices, in addition to informing targeted protein engineering. Here we applied molecular dynamics simulations to explore conformational and substrate binding dynamics in two well-characterized BGs of bacterial (Clostridium cellulovorans) and fungal (Trichoderma reesei) origins as a function of pH. The enzymes were simulated in an explicit solvated environment, with NaCl as electrolytes, at their prominent ionization states obtained at pH 5, 6, 7, and 7.5. Our findings indicated that pH-dependent changes in the ionization states of non-catalytic residues localized outside of the immediate active site led to pH-dependent disruption of the active site conformation. This disruption interferes with favorable H-bonding interactions with catalytic residues required to initiate catalysis on the substrate. We also identified specific non-catalytic residues that are involved in stabilizing the substrate at the optimal pH for enzyme activity. The simulations further revealed the dynamics of water-bridging interactions both outside and inside the substrate binding cleft during structural changes in the enzyme-substrate complex. These findings provide new structural insights into the pH-dependent substrate binding specificity in BGs.

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

In silico studies and fluorescence binding assays of potential anti-prion compounds reveal an important binding site for prion inhibition from PrP(C) to PrP(Sc)

To understand the pharmacophore properties of 2-aminothiazoles and design novel inhibitors against the prion protein, a highly predictive 3D quantitative structure-activity relationship (QSAR) has been developed by performing comparative molecular field analysis (CoMFA) and comparative similarity analysis (CoMSIA). Both CoMFA and CoMSIA maps reveal the presence of the oxymethyl groups in meta and para positions on the phenyl ring of compound 17 (N-[4-(3,4-dimethoxyphenyl)-1,3-thiazol-2-yl]quinolin-2-amine), is necessary for activity while electro-negative nitrogen of quinoline is highly favorable to enhance activity. The blind docking results for these compounds show that the compound with quinoline binds with higher affinity than isoquinoline and naphthalene groups. Out of 150 novel compounds retrieved using finger print analysis by pharmacophoric model predicted based on five test sets of compounds, five compounds with diverse scaffolds were selected for biological evaluation as possible PrP inhibitors. Molecular docking combined with fluorescence quenching studies show that these compounds bind to pocket-D of SHaPrP near Trp145. The new antiprion compounds 3 and 6, which bind with the interaction energies of -12.1 and -13.2 kcal/mol, respectively, show fluorescence quenching with binding constant (Kd) values of 15.5 and 44.14 μM, respectively. Further fluorescence binding assays with compound 5, which is similar to 2-aminothiazole as a positive control, also show that the molecule binds to the pocket-D with the binding constant (Kd) value of 84.7 μM. Finally, both molecular docking and a fluorescence binding assay of noscapine as a negative control reveals the same binding site on the surface of pocket-A near a rigid loop between β2 and α2 interacting with Arg164. This high level of correlation between molecular docking and fluorescence quenching studies confirm that these five compounds are likely to act as inhibitors for prion propagation while noscapine might act as a prion accelerator from PrP(C) to PrP(Sc).

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.

Toward quantitative estimates of binding affinities for protein-ligand systems involving large inhibitor compounds: a steered molecular dynamics simulation route

Understanding binding mechanisms between enzymes and potential inhibitors and quantifying protein-ligand affinities in terms of binding free energy is of primary importance in drug design studies. In this respect, several approaches based on molecular dynamics simulations, often combined with docking techniques, have been exploited to investigate the physicochemical properties of complexes of pharmaceutical interest. Even if the geometric properties of a modeled protein-ligand complex can be well predicted by computational methods, it is still challenging to rank with chemical accuracy a series of ligand analogues in a consistent way. In this article, we face this issue calculating relative binding free energies of a focal adhesion kinase, an important target for the development of anticancer drugs, with pyrrolopyrimidine-based ligands having different inhibitory power. To this aim, we employ steered molecular dynamics simulations combined with nonequilibrium work theorems for free energy calculations. This technique proves very powerful when a series of ligand analogues is considered, allowing one to tackle estimation of protein-ligand relative binding free energies in a reasonable time. In our cases, the calculated binding affinities are comparable with those recovered from experiments by exploiting the Michaelis-Menten mechanism with a competitive inhibitor.

Decrypting Prion Protein Conversion into a β-Rich Conformer by Molecular Dynamics

Prion diseases are fatal neurodegenerative diseases characterized by the formation of β-rich oligomers and the accumulation of amyloid fibrillar deposits in the central nervous system. Understanding the conversion of the cellular prion protein into its β-rich polymeric conformers is fundamental to tackling the early stages of the development of prion diseases. In this paper, we have identified unfolding and refolding steps critical to the conversion into a β-rich conformer for different constructs of the ovine prion protein by molecular dynamics simulations. By combining our results with in vitro experiments, we show that the folded C-terminus of the ovine prion protein is able to recurrently undergo a drastic conformational change by displacement of the H1 helix, uncovering of the H2H3 domain, and formation of persistent β-sheets between H2 and H3 residues. The observed β-sheets refold toward the C-terminus exposing what we call a "bending region" comprising residues 204-214. This is strikingly coincident with the region harboring mutations determining the fate of the prion oligomerization process. The β-rich intermediate is used here for the construction of a putative model for the assembly into an oligomeric aggregate. The results presented here confirm the importance of the H2H3 domain for prion oligomer formation and therefore its potential use as molecular target in the design of novel prion inhibitors.

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

Role of individual histidines in the pH-dependent global stability of human chloride intracellular channel 1

Chloride intracellular channel proteins exist in both a soluble cytosolic form and a membrane-bound form. The mechanism of conversion between the two forms is not properly understood, although one of the contributing factors is believed to be the variation in pH between the cytosol (~7.4) and the membrane (~5.5). We systematically mutated each of the three histidine residues in CLIC1 to an alanine at position 74 and a phenylalanine at positions 185 and 207. We examined the effect of the histidine-mediated pH dependence on the structure and global stability of CLIC1. None of the mutations were found to alter the global structure of the protein. However, the stability of H74A-CLIC1 and H185F-CLIC1, as calculated from the equilibrium unfolding data, is no longer dependent on pH because similar trends are observed at pH 7.0 and 5.5. The crystal structures show that the mutations result in changes in the local hydrogen bond coordination. Because the mutant total free energy change upon unfolding is not different from that of the wild type at pH 7.0, despite the presence of intermediates that are not seen in the wild type, we propose that it may be the stability of the intermediate state rather than the native state that is dependent on pH. On the basis of the lower stability of the intermediate in the H74A and H185F mutants compared to that of the wild type, we conclude that both His74 and His185 are involved in triggering the pH changes to the conformational stability of wild-type CLIC1 via their protonation, which stabilizes the intermediate state.

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

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

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.

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.

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.

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.