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

Cumulative phylogenetic, sequence and structural analysis of Insulin superfamily proteins provide unique structure-function insights

The insulin superfamily proteins (ISPs), in particular, insulin, IGFs and relaxin proteins are key modulators of animal physiology. They are known to have evolved from the same ancestral gene and have diverged into proteins with varied sequences and distinct functions, but maintain a similar structural architecture stabilized by highly conserved disulphide bridges. The recent surge of sequence data and the structures of these proteins prompted a need for a comprehensive analysis, which connects the evolution of these sequences (427 sequences) in the light of available functional and structural information including representative complex structures of ISPs with their cognate receptors. This study reveals (a) unusually high sequence conservation of IGFs (>90 % conservation in 184 sequences) and provides a possible structure-based rationale for such high sequence conservation; (b) provides an updated definition of the receptor-binding signature motif of the functionally diverse relaxin family members (c) provides a probable non-canonical C-peptide cleavage site in a few insulin sequences. The high conservation of IGFs appears to represent a classic case of resistance to sequence diversity exerted by physiologically important interactions with multiple partners. We also propose a probable mechanism for C-peptide cleavage in a few distinct insulin sequences and redefine the receptor-binding signature motif of the relaxin family. Lastly, we provide a basis for minimally modified insulin mutants with potential therapeutic application, inspired by concomitant changes observed in other insulin superfamily protein members supported by molecular dynamics simulation.

Alternate conformations found in protein structures implies biological functions: A case study using cyclophilin A

Protein dynamics linked to numerous biomolecular functions, such as ligand binding, allosteric regulation, and catalysis, must be better understood at the atomic level. Reactive atoms of key residues drive a repertoire of biomolecular functions by flipping between alternate conformations or conformational substates, seldom found in protein structures. Probing such sparsely sampled alternate conformations would provide mechanistic insight into many biological functions. We are therefore interested in evaluating the instance of amino acids adopted alternate conformations, either in backbone or side-chain atoms or in both. Accordingly, over 70000 protein structures appear to contain alternate conformations only 'A' and 'B' for any atom, particularly the instance of amino acids that adopted alternate conformations are more for Arg, Cys, Met, and Ser than others. The resulting protein structure analysis depicts that amino acids with alternate conformations are mainly found in the helical and β-regions and are often seen in high-resolution X-ray crystal structures. Furthermore, a case study on human cyclophilin A (CypA) was performed to explain the pre-existing intrinsic dynamics of catalytically critical residues from the CypA and how such intrinsic dynamics perturbed upon Ser99Thr mutation using molecular dynamics simulations on the ns-μs timescale. Simulation results demonstrated that the Ser99Thr mutation had impaired the alternate conformations or the catalytically productive micro-environment of Phe113, mimicking the experimentally observed perturbation captured by X-ray crystallography. In brief, a deeper comprehension of alternate conformations adopted by the amino acids may shed light on the interplay between protein structure, dynamics, and function.

Genetic Variations in the Purinergic P2X7 Receptor Are Associated with the Immune Response to Ocular Toxoplasmosis in Colombia

Ocular toxoplasmosis (OT) is characterized by inflammation within the eye and is the most recognized clinical manifestation of toxoplasmosis. The objective of this study was to identify new single-nucleotide polymorphisms (SNPs) in the P2RX7 gene that may have significance in the immune response to OT in Colombian patients. A case-control study was conducted to investigate the associations between SNPs (rs1718119 and rs2230912) in the P2RX7 gene and OT in 64 Colombian patients with OT and 64 controls. Capillary electrophoresis was used to analyze the amplification products, and in silico algorithms were employed to predict deleterious SNPs. Stability analysis of amino acid changes indicated that both mutations could lead to decreased protein structure stability. A nonsynonymous SNP, Gln460Arg, located in the long cytoplasmic tail of the receptor, showed a significant association with OT (Bonferroni correction (BONF) = 0.029; odds ratio OR = 3.46; confidence interval CI: 1.05 to 11.39), while no significant association between rs1718119 and OT risk was observed. Based on the 3D structure analysis of the P2RX7 protein trimer, it is hypothesized that an increase in the flexibility of the cytoplasmic domain of this receptor could alter its function. This SNP could potentially serve as a biomarker for identifying Colombian patients at risk of OT.

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.

Homodimeric complexes of the 90-231 human prion: a multilayered computational study based on FMO/GRID-DRY approach

The molecular interaction properties and aggregation capabilities disclosed by PrP-E200K, a pathogenic mutant of the human prion protein, were investigated in detail using multilayered computational approaches. In a previous work, we reported that the electrostatic complementarity between region1 (negative) and region3 (positive) has been assumed to lead to a head-to tail interaction between 120 and 231 PrP-E200K units and to initiation of the aggregation process. In this work, we extended the PrP-E200K structure by including the unstructured 90-120 segment which was found to assume different conformations. Plausible models of 90-231 PrP-E200K dimers were calculated and analyzed in depth to identify the nature of the involved protein-protein interactions. The unstructured 90-120 segment was found to extend the positively charged region3 involved in the association of PrP-E200K units which resulted to be driven by hydrophobic interactions. The combination of molecular dynamics, protein-protein docking, grid-based mapping, and fragment molecular orbital approaches allowed us to provide a plausible mechanism of the early state of 90-231 PrP-E200K aggregation, considered a preliminary step of amyloid conversion.

Investigation of the conformation of human prion protein in ethanol solution using molecular dynamics simulations

When the conformation of protein is changed from its natural state to a misfolded state, some diseases will happen like prion disease. Prion diseases are a set of deadly neurodegenerative diseases caused by prion protein misfolding and aggregation. Monohydric alcohols have a strong influence on the structure of protein. However, whether monohydric alcohols inhibit amyloid fibrosis remains uncertain. Here, to elucidate the effect of ethanol on the structural stability of human prion protein, molecular dynamics simulations were employed to analyze the conformational changes and dynamics characteristics of human prion proteins at different temperatures. The results show that the extension of β-sheet occurs more easily and the α-helix is more easily disrupted at high temperatures. We found that ethanol can destroy the hydrophobic interactions and make the hydrogen bonds stable, which protects the secondary structure of the protein, especially at 500 K.Communicated by Ramaswamy H. Sarma.