Prevalent Mutations of Human Prion Protein: A Molecular Modeling and Molecular Dynamics Study

Human prion protein: exploring the thermodynamic stability and structural dynamics of its pathogenic mutants

Human familial prion diseases are known to be associated with different single-point mutants of the gene coding for prion protein with a primary focus at several locations of the globular domain. We have identified 12 different single-point pathogenic mutants of human prion protein (HuPrP) with the help of extensive perturbations/mutation technique at multiple locations of HuPrP sequence related to potentiality towards conformational disorders. Among these, some of the mutants include pathogenic variants that corroborate well with the literature reported proteins while majority include some unique single-point mutants that are either not explicitly studied early or studied for variants with different residues at the specific position. Primarily, our study sheds light on the unfolding mechanism of the above mentioned mutants in depth. Besides, we could identify some mutants under investigation that demonstrates not only unfolding of the helical structures but also extension and generation of the β-sheet structures and or simultaneously have highly exposed hydrophobic surface which is assumed to be linked with the production of aggregate/fibril structures of the prion protein. Among the identified mutants, Q212E needs special attention due to its maximum exposure of hydrophobic core towards solvent and E200Q is found to be important due to its maximum extent of β-content. We are also able to identify different respective structural conformations of the proteins according to their degree of structural unfolding and those conformations can be extracted and further studied in detail. Communicated by Ramaswamy H. Sarma.

Structural and dynamical determinants of a β-sheet-enriched intermediate involved in amyloid fibrillar assembly of human prion protein

The conformational conversion of the cellular prion protein (PrP^C) into a misfolded, aggregated and infectious scrapie isoform is associated with prion disease pathology and neurodegeneration. Despite the significant number of experimental and theoretical studies the molecular mechanism regulating this structural transition is still poorly understood. Here, via Nuclear Magnetic Resonance (NMR) methodologies we investigate at the atomic level the mechanism of the human HuPrP(90–231) thermal unfolding and characterize the conformational equilibrium between its native structure and a β-enriched intermediate state, named β-PrPI. By comparing the folding mechanisms of metal-free and Cu²⁺-bound HuPrP(23–231) and HuPrP(90–231) we show that the coupling between the N- and C-terminal domains, through transient electrostatic interactions, is the key molecular process in tuning long-range correlated μs–ms dynamics that in turn modulate the folding process. Moreover, via thioflavin T (ThT)-fluorescence fibrillization assays we show that β-PrPI is involved in the initial stages of PrP fibrillation, overall providing a clear molecular description of the initial phases of prion misfolding. Finally, we show by using Real-Time Quaking-Induced Conversion (RT-QuIC) that the β-PrPI acts as a seed for the formation of amyloid aggregates with a seeding activity comparable to that of human infectious prions.

The N-ter domain in HuPrP regulates the folding mechanism by tuning the long-range μs–ms dynamics. Removal of the N-ter domain triggers the formation of a stable β-enriched intermediate state inducing amyloid aggregates with HuPrP^Sc seeding activity.

Tuning Structure and Dynamics of Blue Copper Azurin Junctions via Single Amino-Acid Mutations

In the growing field of biomolecular electronics, blue-copper Azurin stands out as one of the most widely studied protein in single-molecule contacts. Interestingly, despite the paramount importance of the structure/dynamics of molecular contacts in their transport properties, these factors remain largely unexplored from the theoretical point of view in the context of single Azurin junctions. Here we address this issue using all-atom Molecular Dynamics (MD) of Pseudomonas Aeruginosa Azurin adsorbed to a Au(111) substrate. In particular, we focus on the structure and dynamics of the free/adsorbed protein and how these properties are altered upon single-point mutations. The results revealed that wild-type Azurin adsorbs on Au(111) along two well defined configurations: one tethered via cysteine groups and the other via the hydrophobic pocket surrounding the Cu 2 + . Surprisingly, our simulations revealed that single amino-acid mutations gave rise to a quenching of protein vibrations ultimately resulting in its overall stiffening. Given the role of amino-acid vibrations and reorientation in the dehydration process at the protein-water-substrate interface, we suggest that this might have an effect on the adsorption process of the mutant, giving rise to new adsorption configurations.

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

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

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

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

In silico analysis of prion protein mutants: a comparative study by molecular dynamics approach

Polymorphisms in the human prion proteins lead to amino acid substitutions by the conversion of PrPC to PrPSc and amyloid formation, resulting in prion diseases such as familial Creutzfeldt-Jakob disease, Gerstmann-Straussler-Scheinker disease and fatal familial insomnia. Cation-π interaction is a non-covalent binding force that plays a significant role in protein stability. Here, we employ a novel approach by combining various in silico tools along with molecular dynamics simulation to provide structural and functional insight into the effect of mutation on the stability and activity of mutant prion proteins. We have investigated impressions of prevalent mutations including 1E1S, 1E1P, 1E1U, 1E1P, 1FKC and 2K1D on the human prion proteins and compared them with wild type. Structural analyses of the models were performed with the aid of molecular dynamics simulation methods. According to our results, frequently occurred mutations were observed in conserved sequences of human prion proteins and the most fluctuation values appear in the 2K1D mutant model at around helix 4 with residues ranging from 190 to 194. Our observations in this study could help to further understand the structural stability of prion proteins.

Exploring the mechanism of a regulatory SNP of KLK3 by molecular dynamics simulation

The SNP -158G>A of KLK3 has been validated as a regulatory SNP (rSNP) by molecular biology assays, but the mechanism of how it affects the binding of an androgen receptor (AR) homodimer with DNA is unclear. In the current study, molecular dynamics simulation was adopted to explain its inner cause. Based on a recent review), three types of intermolecular forces were analyzed, and the differences among them were compared between complexes containing -158 A:T and -158 G:C. Extra hydrophobic contacts caused by the methyl group on the mutated thymine were the most crucial factor to the regulatory effect of this rSNP. Further analysis concerning the relative motion of the two recognition helixes of the AR homodimer indicated that the hydrophobic interactions between the recognition helix B and the major groove containing -158 A:T changed that helix's motion greatly from swaying in a plane at free state to vibrating slightly around an equilibrium position. A relatively full explanation on the occurrence of rSNP -158G>A is presented here.

Homology modeling and docking studies of FabH (β-ketoacyl-ACP synthase III) enzyme involved in type II fatty acid biosynthesis of Chlorella variabilis: a potential algal feedstock for biofuel production

The concept of using microalgae as an alternative renewable source of biofuel has gained much importance in recent years. However, its commercial feasibility is still an area of concern for researchers. Unraveling the fatty acid metabolic pathway and understanding structural features of various key enzymes regulating the process will provide valuable insights to target microalgae for augmented oil content. FabH (β-ketoacyl-acyl carrier protein synthase; KAS III) is a condensing enzyme catalyzing the initial elongation step of type II fatty acid biosynthetic process and acyl carrier protein (ACP) facilitates the shuttling of the fatty acyl intermediates to the active site of the respective enzymes in the pathway. In the present study, a reliable three-dimensional structure of FabH from Chlorella variabilis, an oleaginous green microalga was modeled and subsequently the key residues involved in substrate binding were determined by employing protein-protein docking and molecular dynamics (MD) simulation protocols. The FabH-ACP complex having the lowest docking energy score showed the binding of ACP to the electropositive FabH surface with strong hydrogen bond interactions. The MD simulation results indicated that the substrate-complexed FabH adopted a more stable conformation than the free enzyme. Further, the FabH structure retained its stability throughout the simulation although noticeable displacements were observed in the loop regions. Molecular simulation studies suggested the importance of crucial hydrogen bonding of the conserved Arg(91) of FabH with Glu(53) and Asp(56) of ACP for exhibiting high affinity between the enzyme and substrate. The molecular modeling results are consistent with available experimental results on the flexibility of FabH and the present study provides first in silico insights into the structural and dynamical aspect of catalytic mechanism of FabH, which could be used for further site-specific mutagenic experiments to develop engineered high oil-yielding microalgal strains for biofuel production.

The effect of fucoidan on tyrosinase: computational molecular dynamics integrating inhibition kinetics

Fucoidan is a complex sulfated polysaccharide extracted from brown seaweed and has a wide variety of biological activities. In this study, we investigated the inhibitory effect of fucoidan on tyrosinase via a combination of inhibition kinetics and computational simulations. Fucoidan reversibly inhibited tyrosinase in a mixed-type manner. Time-interval kinetics showed that the inhibition was processed as first order with biphasic processes. For further insight, we simulated dockings with various sizes of molecular models (monomer to decamer) of fucoidan and showed that the best binding energy change results were obtained from the pentamer (-1.89 kcal/mol) and the hexamer (-1.97 kcal/mol) models of AutoDock Vina. The molecular dynamics simulation confirmed the binding mechanisms between tyrosinase and fucoidan and suggested that fucoidan mostly interacts with several residues including copper ions located in the active site. Our study suggests that fucoidan might be a potential natural antipigment agent.

A molecular history of the amyloidoses

The molecular investigation of the amyloidoses began in the mid-19th century with the observation of areas in human tissues obtained at autopsy that were homogeneous and eosinophilic with conventional stains but became blue when exposed to mixtures of iodine and sulfuric acid. The foci corresponded to regions formerly identified as "waxy" or lardaceous. Subsequent identification of the characteristic staining of the same tissues with metachromatic dyes such as crystal violet or with the cotton dye Congo red (particularly under polarized light) and thioflavins allowed the pathological classification of those tissues as belonging to a set of disorders known as the amyloidoses. Not unexpectedly, progress has reflected evolving technology and parallel advances in all fields of biological science. Investigation using contemporary methods has expanded our notions of amyloid proteins from being simply agents or manifestations of systemic, largely extracellular diseases to include "protein-only infection," the concept that "normal" functional amyloids might exist in eukaryotes and prokaryotes and that aggregatability may be an intrinsic structural price to be paid for some functional protein domains. We now distinguish between the amyloidoses, that is, diseases caused by the deposition of amyloid fibrils and amyloid proteins (i.e., purified or recombinant proteins that form amyloid fibrils in vitro), which may or may not be associated with disease in vivo.