Mimicking the action of folding chaperones in molecular dynamics simulations: Application to the refinement of homology-based protein structures

Small Ubiquitin-Like Modifier Protein 3 Enhances the Solubilization of Human Bone Morphogenetic Protein 2 in E. coli

Small ubiquitin-like modifier (SUMO) fusion technology is widely used in the production of heterologous proteins from prokaryotic system to aid in protein solubilization and refolding. Due to an extensive clinical application of human bone morphogenetic protein 2 (hBMP2) in bone augmentation, total RNA was isolated from human gingival tissue and mature gene was amplified through RT-PCR, cloned (pET21a), sequence analyzed, and submitted to GenBank (Accession no. KF250425). To obtain soluble expression, SUMO3 was tagged at the N-terminus of hBMP2 gene (pET21a/SUMO3-hBMP2), transferred in BL21 codon+, and ~ 40% soluble expression was obtained on induction with IPTG. The dimerized hBMP2 was confirmed with Western blot, native PAGE analysis, and purified by fast protein liquid chromatography with 0.5 M NaCl elution. The cleavage of SUMO3 tag from hBMP2 converted it to an insoluble form. Computational 3D structural analysis of the SUMO3-hBMP2 was performed and optimized by molecular dynamic simulation. Protein-protein interaction of SUMO3-hBMP2 with BMP2 receptor was carried out using HADDOCK and inferred stable interaction. The alkaline phosphatase assay of SUMO3-hBMP2 on C2C12 cells showed maximum 200-ng/ml dose-dependent activity. We conclude that SUMO3-tagged hBMP2 is more suited for generation of soluble form of the protein and addition of SUMO3 tag does not affect the functional activity of hBMP2.

Identification of a 3rd Na+ Binding Site of the Glycine Transporter, GlyT2

The Na+/Cl- dependent glycine transporters GlyT1 and GlyT2 regulate synaptic glycine concentrations. Glycine transport by GlyT2 is coupled to the co-transport of three Na+ ions, whereas transport by GlyT1 is coupled to the co-transport of only two Na+ ions. These differences in ion-flux coupling determine their respective concentrating capacities and have a direct bearing on their functional roles in synaptic transmission. The crystal structures of the closely related bacterial Na+-dependent leucine transporter, LeuTAa, and the Drosophila dopamine transporter, dDAT, have allowed prediction of two Na+ binding sites in GlyT2, but the physical location of the third Na+ site in GlyT2 is unknown. A bacterial betaine transporter, BetP, has also been crystallized and shows structural similarity to LeuTAa. Although betaine transport by BetP is coupled to the co-transport of two Na+ ions, the first Na+ site is not conserved between BetP and LeuTAa, the so called Na1' site. We hypothesized that the third Na+ binding site (Na3 site) of GlyT2 corresponds to the BetP Na1' binding site. To identify the Na3 binding site of GlyT2, we performed molecular dynamics (MD) simulations. Surprisingly, a Na+ placed at the location consistent with the Na1' site of BetP spontaneously dissociated from its initial location and bound instead to a novel Na3 site. Using a combination of MD simulations of a comparative model of GlyT2 together with an analysis of the functional properties of wild type and mutant GlyTs we have identified an electrostatically favorable novel third Na+ binding site in GlyT2 formed by Trp263 and Met276 in TM3, Ala481 in TM6 and Glu648 in TM10.

Mimicking the action of folding chaperones by Hamiltonian replica-exchange molecular dynamics simulations: application in the refinement of de novo models

The efficiency of using a variant of Hamiltonian replica-exchange molecular dynamics (Chaperone H-replica-exchange molecular dynamics [CH-REMD]) for the refinement of protein structural models generated de novo is investigated. In CH-REMD, the interaction between the protein and its environment, specifically, the electrostatic interaction between the protein and the solvating water, is varied leading to cycles of partial unfolding and refolding mimicking some aspects of folding chaperones. In 10 of the 15 cases examined, the CH-REMD approach sampled structures in which the root-mean-square deviation (RMSD) of secondary structure elements (SSE-RMSD) with respect to the experimental structure was more than 1.0 Å lower than the initial de novo model. In 14 of the 15 cases, the improvement was more than 0.5 Å. The ability of three different statistical potentials to identify near-native conformations was also examined. Little correlation between the SSE-RMSD of the sampled structures with respect to the experimental structure and any of the scoring functions tested was found. The most effective scoring function tested was the DFIRE potential. Using the DFIRE potential, the SSE-RMSD of the best scoring structures was on average 0.3 Å lower than the initial model. Overall the work demonstrates that targeted enhanced-sampling techniques such as CH-REMD can lead to the systematic refinement of protein structural models generated de novo but that improved potentials for the identification of near-native structures are still needed.

Mechanisms of constitutive activation of Janus kinase 2-V617F revealed at the atomic level through molecular dynamics simulations

BACKGROUND

The tyrosine kinase Janus kinase 2 (JAK2) is important in triggering nuclear translocation and regulation of target genes expression through signal transducer and activator of transcription pathways. The valine-to-phenylalanine mutation at amino acid 617 (V617F), which results in the deregulation of JAK2, has been implicated in the oncogenesis of chronic myeloproliferative disease. However, both the mechanism of JAK2 autoinhibition and the mechanism of V617F constitutive activation remain unclear.

METHOD

In this work, the authors used molecular dynamics simulation techniques to establish plausible mechanisms of JAK2 autoinhibition and V617F constitutive activation at the atomic level.

RESULTS

In wild-type JAK2, the activation loop of JAK2-homology domain 1 (JH1) is pulled toward the JH1/JH2 interface through interactions with key residues of JH2, especially S591, F595, and V617, and stabilizes the inactivated form of JH1. In the case of V617F, through the aromatic ring-ring stacking interaction, F617 blocks the interaction of JH1 the activation loop, S591, and F595, thus causing the JH1 activation loop to move back to its activated form.

CONCLUSIONS

The current results indicated that this simulation-derived mechanism of JAK2 autoregulation is consistent with current available experimental evidence and may lead to a deeper understanding of JAK2 and other kinase systems that are regulated by pseudokinases.

Refining homology models by combining replica-exchange molecular dynamics and statistical potentials

A protocol is presented for the global refinement of homology models of proteins. It combines the advantages of temperature-based replica-exchange molecular dynamics (REMD) for conformational sampling and the use of statistical potentials for model selection. The protocol was tested using 21 models. Of these 14 were models of 10 small proteins for which high-resolution crystal structures were available, the remainder were targets of the recent CASPR exercise. It was found that REMD in combination with currently available force fields could sample near-native conformational states starting from high-quality homology models. Conformations in which the backbone RMSD of secondary structure elements (SSE-RMSD) was lower than the starting value by 0.5-1.0 A were found for 15 out of the 21 cases (average 0.82 A). Furthermore, when a simple scoring function consisting of two statistical potentials was used to rank the structures, one or more structures with SSE-RMSD of at least 0.2 A lower than the starting value was found among the five best ranked structures in 11 out of the 21 cases. The average improvement in SSE-RMSD for the best models was 0.42 A. However, none of the scoring functions tested identified the structures with the lowest SSE-RMSD as the best models although all identified the native conformation as the one with lowest energy. This suggests that while the proposed protocol proved effective for the refinement of high-quality models of small proteins scoring functions remain one of the major limiting factors in structure refinement. This and other aspects by which the methodology could be further improved are discussed.

Structural studies of the Enterococcus faecalis SufU [Fe-S] cluster protein

Background

Iron-sulfur clusters are ubiquitous and evolutionarily ancient inorganic prosthetic groups, the biosynthesis of which depends on complex protein machineries. Three distinct assembly systems involved in the maturation of cellular Fe-S proteins have been determined, designated the NIF, ISC and SUF systems. Although well described in several organisms, these machineries are poorly understood in Gram-positive bacteria. Within the Firmicutes phylum, the Enterococcus spp. genus have recently assumed importance in clinical microbiology being considered as emerging pathogens for humans, wherein Enterococcus faecalis represents the major species associated with nosocomial infections. The aim of this study was to carry out a phylogenetic analysis in Enterococcus faecalis V583 and a structural and conformational characterisation of it SufU protein.

Results

BLAST searches of the Enterococcus genome revealed a series of genes with sequence similarity to the Escherichia coli SUF machinery of [Fe-S] cluster biosynthesis, namely sufB, sufC, sufD and SufS. In addition, the E. coli IscU ortholog SufU was found to be the scaffold protein of Enterococcus spp., containing all features considered essential for its biological activity, including conserved amino acid residues involved in substrate and/or co-factor binding (Cys^50,76,138 and Asp⁵²) and, phylogenetic analyses showed a close relationship with orthologues from other Gram-positive bacteria. Molecular dynamics for structural determinations and molecular modeling using E. faecalis SufU primary sequence protein over the PDB:1su0 crystallographic model from Streptococcus pyogenes were carried out with a subsequent 50 ns molecular dynamic trajectory. This presented a stable model, showing secondary structure modifications near the active site and conserved cysteine residues. Molecular modeling using Haemophilus influenzae IscU primary sequence over the PDB:1su0 crystal followed by a MD trajectory was performed to analyse differences in the C-terminus region of Gram-positive SufU and Gram-negative orthologous proteins, in which several modifications in secondary structure were observed.

Conclusion

The data describe the identification of the SUF machinery for [Fe-S] cluster biosynthesis present in the Firmicutes genome, showing conserved sufB, sufC, sufD and sufS genes and the presence of the sufU gene coding for scaffold protein, instead of sufA; neither sufE nor sufR are present. Primary sequences and structural analysis of the SufU protein demonstrated its structural-like pattern to the scaffold protein IscU nearby on the ISC machinery. E. faecalis SufU molecular modeling showed high flexibility over the active site regions, and demonstrated the existence of a specific region in Firmicutes denoting the Gram positive region (GPR), suggested as a possible candidate for interaction with other factors and/or regulators.

In silico chaperonin-like cycle helps folding of proteins for structure prediction

Currently, one of the most serious problems in protein-folding simulations for de novo structure prediction is conformational sampling of medium-to-large proteins. In vivo, folding of these proteins is mediated by molecular chaperones. Inspired by the functions of chaperonins, we designed a simple chaperonin-like simulation protocol within the framework of the standard fragment assembly method: in our protocol, the strength of the hydrophobic interaction is periodically modulated to help the protein escape from misfolded structures. We tested this protocol for 38 proteins and found that, using a certain defined criterion of success, our method could successfully predict the native structures of 14 targets, whereas only those of 10 targets were successfully predicted using the standard protocol. In particular, for non-alpha-helical proteins, our method yielded significantly better predictions than the standard approach. This chaperonin-inspired protocol that enhanced de novo structure prediction using folding simulations may, in turn, provide new insights into the working principles underlying the chaperonin system.

The Proteomic Code: a molecular recognition code for proteins

Background

The Proteomic Code is a set of rules by which information in genetic material is transferred into the physico-chemical properties of amino acids. It determines how individual amino acids interact with each other during folding and in specific protein-protein interactions. The Proteomic Code is part of the redundant Genetic Code.

Review

The 25-year-old history of this concept is reviewed from the first independent suggestions by Biro and Mekler, through the works of Blalock, Root-Bernstein, Siemion, Miller and others, followed by the discovery of a Common Periodic Table of Codons and Nucleic Acids in 2003 and culminating in the recent conceptualization of partial complementary coding of interacting amino acids as well as the theory of the nucleic acid-assisted protein folding.

Methods and conclusions

A novel cloning method for the design and production of specific, high-affinity-reacting proteins (SHARP) is presented. This method is based on the concept of proteomic codes and is suitable for large-scale, industrial production of specifically interacting peptides.

Hydrophobic aided replica exchange: an efficient algorithm for protein folding in explicit solvent

A hydrophobic aided replica exchange method (HAREM) is introduced to accelerate the simulation of all-atom protein folding in explicit solvent. This method is based on exaggerating the hydrophobic effect of various protein amino acids in water by attenuating the protein-water attractive interactions (mimicking the Chaperon effect) while leaving other interactions among protein atoms and water molecules unchanged. The method is applied to a small representative protein, the alpha-helix 3K(I), and it is found that the HAREM method successfully folds the protein within 4 ns, while the regular replica exchange method does not fold the same protein within 5 ns, even with many more replicas.

Characterization of the structure of RAMP1 by mutagenesis and molecular modeling

Receptor activity modifying proteins (RAMPs) are a family of single-pass transmembrane proteins that dimerize with G-protein-coupled receptors. They may alter the ligand recognition properties of the receptors (particularly for the calcitonin receptor-like receptor, CLR). Very little structural information is available about RAMPs. Here, an ab initio model has been generated for the extracellular domain of RAMP1. The disulfide bond arrangement (Cys27-Cys82, Cys40-Cys72, and Cys57-Cys104) was determined by site-directed mutagenesis. The secondary structure (alpha-helices from residues 29-51, 60-80, and 87-100) was established from a consensus of predictive routines. Using these constraints, an assemblage of 25,000 structures was constructed and these were ranked using an all-atom statistical potential. The best 1000 conformations were energy minimized. The lowest scoring model was refined by molecular dynamics simulation. To validate our strategy, the same methods were applied to three proteins of known structure; PDB:1HP8, PDB:1V54 chain H (residues 21-85), and PDB:1T0P. When compared to the crystal structures, the models had root mean-square deviations of 3.8 A, 4.1 A, and 4.0 A, respectively. The model of RAMP1 suggested that Phe93, Tyr100, and Phe101 form a binding interface for CLR, whereas Trp74 and Phe92 may interact with ligands that bind to the CLR/RAMP1 heterodimer.

Amino acid size, charge, hydropathy indices and matrices for protein structure analysis

BACKGROUND

Prediction of protein folding and specific interactions from only the sequence (ab initio) is a major challenge in bioinformatics. It is believed that such prediction will prove possible if Anfinsen's thermodynamic principle is correct for all kinds of proteins, and all the information necessary to form a concrete 3D structure is indeed present in the sequence.

RESULTS

We indexed the 200 possible amino acid pairs for their compatibility regarding the three major physicochemical properties--size, charge and hydrophobicity--and constructed Size, Charge and Hydropathy Compatibility Indices and Matrices (SCI & SCM, CCI & CCM, and HCI & HCM). Each index characterized the expected strength of interaction (compatibility) of two amino acids by numbers from 1 (not compatible) to 20 (highly compatible). We found statistically significant positive correlations between these indices and the propensity for amino acid co-locations in real protein structures (a sample containing total 34630 co-locations in 80 different protein structures): for HCI: p < 0.01, n = 400 in 10 subgroups; for SCI p < 1.3E-08, n = 400 in 10 subgroups; for CCI: p < 0.01, n = 175). Size compatibility between residues (well known to exist in nucleic acids) is a novel observation for proteins. Regression analyzes indicated at least 7 well distinguished clusters regarding size compatibility and 5 clusters of charge compatibility. We tried to predict or reconstruct simple 2D representations of 3D structures from the sequence using these matrices by applying a dot plot-like method. The location and pattern of the most compatible subsequences was very similar or identical when the three fundamentally different matrices were used, which indicates the consistency of physicochemical compatibility. However, it was not sufficient to choose one preferred configuration between the many possible predicted options.

CONCLUSION

Indexing of amino acids for major physico-chemical properties is a powerful approach to understanding and assisting protein design. However, it is probably insufficient itself for complete ab initio structure prediction.

Mimicking the action of GroEL in molecular dynamics simulations: application to the refinement of protein structures

Bacterial chaperonin, GroEL, together with its co-chaperonin, GroES, facilitates the folding of a variety of polypeptides. Experiments suggest that GroEL stimulates protein folding by multiple cycles of binding and release. Misfolded proteins first bind to an exposed hydrophobic surface on GroEL. GroES then encapsulates the substrate and triggers its release into the central cavity of the GroEL/ES complex for folding. In this work, we investigate the possibility to facilitate protein folding in molecular dynamics simulations by mimicking the effects of GroEL/ES namely, repeated binding and release, together with spatial confinement. During the binding stage, the (metastable) partially folded proteins are allowed to attach spontaneously to a hydrophobic surface within the simulation box. This destabilizes the structures, which are then transferred into a spatially confined cavity for folding. The approach has been tested by attempting to refine protein structural models generated using the ROSETTA procedure for ab initio structure prediction. Dramatic improvements in regard to the deviation of protein models from the corresponding experimental structures were observed. The results suggest that the primary effects of the GroEL/ES system can be mimicked in a simple coarse-grained manner and be used to facilitate protein folding in molecular dynamics simulations. Furthermore, the results support the assumption that the spatial confinement in GroEL/ES assists the folding of encapsulated proteins.