Molecular dynamics study of the conformational stability of esterase 2 from Alicyclobacillus acidocaldarius

Circular dichroism and differential scanning calorimetry measurements showed that esterase 2 from the thermophilic microorganism Alicyclobacillus acidocaldarius, EST2, and its variant in which the first 35 residues have been deleted, EST2-36del, unfold reversibly on increasing temperature, and possess two cooperative and coupled domains [12]. Structural features of the α/β hydrolase fold of EST2, with nine α-helices packed against the central twisted β-sheet, do not allow a straightforward identification of these two cooperative and coupled domains. Molecular dynamics simulations, each one 20 ns long, have been performed at 300, 400 and 500 K, on both proteins in explicit water. Suitable analysis of MD trajectories has allowed a reliable identification of the two cooperative domains (i.e., the less stable one corresponds to external α-helices, whereas the more stable one corresponds to the central twisted β-sheet) and the attribution of the key coupling role to the last and long α-helix of EST2.     

Modulation of the mobility of a key region in human galactokinase: Impacts on catalysis and stability

Galactokinase catalyses the phosphorylation of α-d-galactose and some structurally related monosaccharides. The enzyme is of interest due to its potential as a biocatalyst for the production of sugar 1-phosphates and due to its involvement in the inherited metabolic disease type II galactosemia. It has been previously shown that a region (residues 231–245) in human galactokinase often has altered mobility when active site residues are varied. We hypothesised that the reverse may be true and that designing changes to this region might affect the functioning of the active site of the enzyme. Focussing on four residues (Leu-231, Gln-242, Glu-244 and Glu-245) we conducted molecular dynamics simulations to explore the effects of changing these residues to glycine or serine. In most cases the variations resulted in local changes to the 231–245 region and global changes to the root mean squared fluctuation (RMSF) of the protein. The four serine variants were expressed as recombinant proteins. All had altered steady state enzyme kinetic parameters with α-d-galactose as a substrate. However, these changes were generally less than ten-fold in magnitude. Changes were also observed with 2-deoxy-α-d-galactose, α-d-galactosamine and α-d-talose as substrates, including (in some cases) loss of detectable activity, suggesting that these variations can tune the specificity of the enzyme. This study demonstrates that activity and specificity of human galactokinase can be modulated by variations designed to affect active site flexibility. It is likely that this principle can be generalised to other enzymes.     

Effect of external stresses on protein conformation: a computer modelling study

The increasing use of digital technologies such as mobile phones has led to major health concerns about the effects of non-ionizing pulsed radiation exposure. We believe that the health implications of exposure to radiation cannot be fully understood without establishing the molecular mechanisms of biological effects of pulsed microwaves. We aim to establish methods for studying the molecular mechanisms of protein structural and energetic changes occurring due to external stresses related to non-ionizing radiation by using a combination of experimental and theoretical approaches. In this paper, we present the results from our fully atomistic simulation study of chemical and thermal stress response of a prototype protein, insulin. We performed a series of molecular dynamics simulations of insulin in solution under equilibrium conditions, under chemical stress (imitated by reducing the disulfide bonds in the protein molecule), and under short-lived thermal stress (imitated by increasing simulation temperature for up to 2 ns). The resultant protein conformational behaviour was analysed for various properties with the aim of establishing analysis routines for classification of protein unfolding pathways and associated molecular mechanisms.Submitted as a record of the 2002 Australian Biophysical Society meeting     

Active-site residues move independently from the rest of the protein in a 200 ns molecular dynamics simulation of cytochrome P450 CYP119

The conformational dynamics of cytochrome P450 enzymes are critical to their catalytic activity. In this study, the correlated motion between residues in a 200 ns molecular dynamics trajectory of the thermophilic CYP119 was analyzed to parse out conformational relationships. Residues that are structurally related, for example residues within a helix, generally have highly correlated motion. In addition, clusters of non-adjacent residues that show correlated motion (“hot spots”) are seen in various regions, including at the base of the F and G helices that make up the most dynamic region of the enzyme. A modified k-means algorithm that clusters residues based on their correlated motion indicates that functionally related residues are in the same cluster (e.g., the catalytic threonines and the heme). Tightly coupled clusters form a solvent-exposed “shell” around the enzyme, whereas less coupling between clusters is seen in regions that are critical to ligand interactions, redox partner interactions, and catalysis. Most notably, we find that residues in the active site move independently from the rest of the enzyme, effectively insulating the catalytic machinery from other regions of the protein.     

Dynamics fingerprint and inherent asymmetric flexibility of a cold-adapted homodimeric enzyme. A case study of the Vibrio alkaline phosphatase

Background

Protein dynamics influence protein function and stability and modulate conformational changes. Such motions depend on the underlying networks of intramolecular interactions and communicating residues within the protein structure. Here, we provide the first characterization of the dynamic fingerprint of the dimeric alkaline phosphatase (AP) from the cold-adapted Vibrio strain G15-21 (VAP), which is among the APs with the highest known k_cat at low temperatures.

Methods

Multiple all-atom explicit solvent molecular dynamics simulations were employed in conjunction with different metrics to analyze the dynamical patterns and the paths of intra- and intermolecular communication.

Results

Interactions and coupled motions at the interface between the two VAP subunits have been characterized, along with the networks of intramolecular interactions. It turns out a low number of intermolecular interactions and coupled motions, which result differently distributed in the two monomers. The paths of long-range communication mediated from the catalytic residues to distal sites were also characterized, pointing out a different information flow in the two subunits.

Conclusions

A pattern of asymmetric flexibility is evident in the two identical subunits of the VAP dimer that is intimately linked to a different distribution of intra- and intermolecular interactions. The asymmetry was also evident in pairs of cross-correlated residues during the dynamics.

General significance

The results here discussed provide a structural rationale to the half-of-site mechanism previously proposed for VAP and other APs, as well as a general framework to characterize asymmetric dynamics in homomeric enzymes.     

Conformational stability of the full‐atom hexameric model of the ClpB chaperone from Escherichia coli

The Escherichia coli heat shock protein ClpB, a member of the Hsp100 family, plays a crucial role in cellular thermotolerance. In co‐operation with the Hsp70 chaperone system, it is able to solubilize proteins aggregated by heat shock conditions and refold them into the native state in an ATP‐dependent way. It was established that the mechanism of ClpB action depends on the formation of a ring‐shaped hexameric structure and the translocation of a protein substrate through an axial channel. The structural aspects of this process are not fully known. By means of homology modeling and protein–protein docking, we obtained a model of the hexameric arrangement of the full‐length ClpB protein complexed with ATP. A molecular dynamics simulation of this model was performed to assess its flexibility and conformational stability. The high mobility of the “linker” M‐domain, essential for the renaturing activity of ClpB, was demonstrated, and the size and shape of central channel were analyzed. In this model, we propose the coordinates for a loop between b4 and B6 structural elements, not defined in previous structural research, which faces the inside of the channel and may therefore play a role in substrate translocation. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 47–60, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Impact of the tail and mutations G131V and M129V on prion protein flexibility

Within the "protein-only" hypothesis, a detailed mechanism for the conversion of a alpha-helix to beta-sheet structure is unclear. We have investigated the effects of the tail 90-123 and the point mutations G131V and M129V on prion protein conformational plasticity at neutral pH. Molecular dynamics simulations show that the dynamics of the core 124-226 is essentially independent of the tail and that the point mutation G131V does not affect PrP thermodynamic stability. Both mutations, however, enhance the flexibility of residues that participate in the two-step process for prion propagation. They also extend the short beta-sheet in the normal protein into a larger sheet at neutral pH. This finding suggests a critical role of the tail for triggering the topological change.     

Tyr74 is essential for the formation, stability and function of Plasmodium falciparum triosephosphate isomerase dimer

Plasmodium falciparum triosephosphate isomerase (PfTIM) is known to be functional only as a homodimer. Although many studies have shown that the interface Cys13 plays a major role in the stability of the dimer, a few reports have demonstrated that structurally conserved Tyr74 may be essential for the stability of PfTIM dimer. To understand the role of Tyr74, we have performed molecular dynamics (MD) simulations of monomeric and dimeric PfTIM mutated to glycine and cysteine at position 74. Simulations of the monomer revealed that mutant Tyr74Gly does not produce changes in folding and stability of the monomer. Interestingly, comparison of the flexibility of Tyr74 in the monomer and dimer revealed that this residue possesses an intrinsic restricted mobility, indicating that Tyr74 is an anchor residue required for homodimerization. Tyr74 also appears to play an important role in binding by facilitating the disorder-to-order transitions of loops 1 and 3, which allows Cys13 to form favorable interactions with loop 3 and Lys12 to be locked in a favorable position for catalysis. High-temperature MD simulations of the wild-type and Tyr74Gly PfTIM dimers showed that the aromatic moiety of Tyr74 is necessary to preserve the geometry and native contacts between loops 1 and 3 at the interface of the dimer. Disulfide cross-linking between mutant Tyr74Cys and Cys13 further revealed that Tyr74 stabilizes the geometry of loop 1 (which contains the catalytic residue Lys12) and the interactions between loops 1 and 3 via aromatic-aromatic interactions with residues Phe69, Tyr101, and Phe102. Principal component analysis showed that Tyr74 is also necessary to preserve the collective motions in the dimer that contribute to the catalytic efficiency of PfTIM dimer. We conclude that Tyr74 not only plays a role in the stability of the dimer, but also participates in the dimerization process and collective motions via coupled disorder-to-order transitions of intrinsically disordered regions, necessary for efficiency in the catalytic function of PfTIM.     

Molecular dynamics simulations of trehalose as a 'dynamic reducer' for solvent water molecules in the hydration shell

Systematic computational work for a series of 13 disaccharides was performed to provide an atomic-level insight of unique biochemical role of the alpha,alpha-(1-->1)-linked glucopyranoside dimer over the other glycosidically linked sugars. Superior osmotic and cryoprotective abilities of trehalose were explained on the basis of conformational and hydration characteristics of the trehalose molecule. Analyses of the hydration number and radial distribution function of solvent water molecules showed that there was very little hydration adjacent to the glycosidic oxygen of trehalose and that the dynamic conformation of trehalose was less flexible than any of the other sugars due to this anisotropic hydration. The remarkable conformational rigidity that allowed trehalose to act as a sugar template was required for stable interactions with hydrogen-bonded water molecules. Trehalose made an average of 2.8 long-lived hydrogen bonds per each MD step, which was much larger than the average of 2.1 for the other sugars. The stable hydrogen-bond network is derived from the formation of long-lived water bridges at the expense of decreasing the dynamics of the water molecules. Evidence for this dynamic reduction of water by trehalose was also established based on each of the lowest translational diffusion coefficients and the lowest intermolecular coulombic energy of the water molecules around trehalose. Overall results indicate that trehalose functions as a 'dynamic reducer' for solvent water molecules based on its anisotropic hydration and conformational rigidity, suggesting that macroscopic solvent properties could be modulated by changes in the type of glycosidic linkages in sugar molecules.     

Peroxisomal proteostasis involves a lon family protein that functions as protease and chaperone

Proteins are subject to continuous quality control for optimal proteostasis. The knowledge of peroxisome quality control systems is still in its infancy. Here we show that peroxisomes contain a member of the Lon family of proteases (Pln). We show that Pln is a heptameric protein and acts as an ATP-fueled protease and chaperone. Hence, Pln is the first chaperone identified in fungal peroxisomes. In cells of a PLN deletion strain peroxisomes contain protein aggregates, a major component of which is catalase-peroxidase. We show that this enzyme is sensitive to oxidative damage. The oxidatively damaged, but not the native protein, is a substrate of the Pln protease. Cells of the pln strain contain enhanced levels of catalase-peroxidase protein but reduced catalase-peroxidase enzyme activities. Together with the observation that Pln has chaperone activity in vitro, our data suggest that catalase-peroxidase aggregates accumulate in peroxisomes of pln cells due to the combined absence of Pln protease and chaperone activities.     

Rabbit indolethylamine N-methyltransferase three-dimensional structure prediction: a model approach to bridge sequence to function in pharmacogenomic studies

Pharmacogenomics is the study of the genetic basis for individual variation in response to drugs and other xenobiotics. Successful prediction of effects of genetic variations that change encoded amino acid sequences on protein function and their consequent biomedical implications depends on three-dimensional (3D) structures of the encoded amino acid sequences. To bridge sequence to function, thus facilitating an in-depth pharmacogenomic study, we tested the feasibility of the use of a semi-computational approach to predict 3D structures of rabbit and human indolethylamine N-methyltransferases (INMTs) from their amino acid sequences, which share less than 26% sequence identity with known protein 3D structures. Herein, we report 3D models of INMTs predicted by using the crystal structure of rat catechol O-methyltransferase as a template, testing of the models both computationally and experimentally, and successful use of the models in retrospective prediction of the effects of genetic polymorphisms and in identification of residues that contribute to observed species-specific differences in substrate affinity. The results encourage the use of the semi-computational approach to predict 3D protein structures for use in pharmacogenomic studies when de novo prediction of protein 3D structures from their amino acid sequences is still not feasible and X-ray crystallography and/or solution nuclear magnetic resonance spectroscopy can only determine 3D structures for a small number of known amino acid sequences.Electronic Supplementary Material available.     

Domain stability in the AAA+ ATPase ClpB from Escherichia coli

ClpB is a heat-shock protein that reactivates aggregated proteins in cooperation with the DnaK chaperone system. ClpB belongs to the family of AAA+ ATPases and forms ring-shaped oligomers: heptamers in the absence of nucleotides and hexamers in the presence of nucleotides. We investigated the thermodynamic stability of ClpB in its monomeric and oligomeric forms. ClpB contains six distinct structural domains: the N-terminal domain involved in substrate binding, two AAA+ ATP-binding modules, each consisting of two domains, and a coiled-coil domain inserted between the AAA+ modules. We produced seven variants of ClpB, each containing a single Trp located in each of the ClpB domains and measured the changes in Trp fluorescence during the equilibrium urea-induced unfolding of ClpB. We found that two structural domains: the small domain of the C-terminal AAA+ module and the coiled-coil domain were destabilized in the oligomeric form of ClpB, which indicates that only those domains change their conformation and/or interactions during formation of the ClpB rings.     

Structural stability of myoglobin and glycomyoglobin: a comparative molecular dynamics simulation study

Glycoproteins are formed as the result of enzymatic glycosylation or chemical glycation in the body, and produced in vitro in industrial processes. The covalently attached carbohydrate molecule(s) confer new properties to the protein, including modified stability. In the present study, the structural stability of a glycoprotein form of myoglobin, bearing a glucose unit in the N-terminus, has been compared with its native form by the use of molecular dynamics simulation. Both structures were subjected to temperatures of 300 and 500 K in an aqueous environment for 10 ns. Changes in secondary structures and RMSD were then assessed. An overall higher stability was detected for glycomyoglobin, for which the most stable segments/residues were highlighted and compared with the native form. The simple addition of a covalently bound glucose is suggested to exert its stabilizing effect via increased contacts with surrounding water molecules, as well as a different pattern of interactions with neighbor residues.

Electronic supplementary material

The online version of this article (doi:10.1007/s10867-015-9383-2) contains supplementary material, which is available to authorized users.     

Disulfide bonds convert small heat shock protein Hsp16.3 from a chaperone to a non-chaperone: implications for the evolution of cysteine in molecular chaperones

Molecular chaperones mainly function in assisting newly synthesized polypeptide folding and protect non-native proteins from aggregation, with known structural features such as the ability of spontaneous folding/refolding and high conformational flexibility. In this report, we verified the assumption that the lack of disulfide bonds in molecular chaperones is a prerequisite for such unique structural features. Using small heat shock protein (one sub-class of chaperones) Hsp16.3 as a model system, our results show the following: (1) Cysteine-free Hsp16.3 wild type protein can efficiently exhibit chaperone activity and spontaneously refold/reassemble with high conformational flexibility. (2) Whereas Hsp16.3 G89C mutant with inter-subunit disulfide bonds formed seems to lose the nature of chaperone proteins, i.e., under stress conditions, it neither acts as molecular chaperone nor spontaneously refolds/reassembles. Structural analysis indicated that the mutant exists as an unstable molten globule-like state, which incorrectly exposes hydrophobic surfaces and irreversibly tends to form aggregates that can be suppressed by the other molecular chaperone (alpha-crystallin). By contrast, reduction of disulfide bond in the Hsp16.3 G89C mutant can significantly recover its character as a molecular chaperone. In light of these results, we propose that disulfide bonds could severely disturb the structure/function of molecular chaperones like Hsp16.3. Our results might not only provide insights into understanding the structural basis of chaperone upon binding substrates, but also explain the observation that the occurrence of cysteine in molecular chaperones is much lower than that in other protein families, subsequently being helpful to understand the evolution of protein family.     

Structure-activity relationship study of tetrapeptide inhibitors of the Vascular Endothelial Growth Factor A binding to Neuropilin-1

Neuropilin-1 is considered as one of the key receptors responsible for signaling pathways involved in pathological angiogenesis necessary for tumor progression, therefore targeting of VEGF₁₆₅ binding to NRP-1 could be a relevant strategy for antiangiogenic treatment. It was shown before that the VEGF₁₆₅/NRP-1 interaction can be inhibited by short tetrapeptides with K/RXXR sequence.Here, we present a structure–activity relationship study of the systematic optimization of amino acid residues in positions 1–3 in the above tetrapeptides. All the 13 synthesized analogs possessed C-terminal arginine that is a necessary element for interaction with NRP-1. The obtained results of the inhibitory activity and modeling by molecular dynamics indicate that simultaneous interactions of the basic amino acid residues in position 1 and 4 (Arg) with Neuropilin-1 are crucial and their cooperation strongly affects the inhibitory activity. In addition, the binding strength is modulated by the flexibility of the peptide backbone (in the central part of the peptide), and the nature of the side chain of the amino acids at the second or third position. A dramatic decrease in the activity to the receptor is observed in flexible derivatives that are missing proline residues. The results described in this paper should prove useful for future studies aimed at establishing the best pharmacophore for inhibitors of VEGF₁₆₅ binding to NRP-1.     

Stability and fluctuations of amide hydrogen bonds in a bacterial cytochrome c: a molecular dynamics study

Molecular dynamics (MD) simulations on a bacterial cytochrome c were performed to investigate the lifetime and fluctuations of backbone hydrogen bonds and to correlate these data with protection factors for hydrogen exchange measured by NMR spectroscopy (Bartalesi et al. in Biochemistry, 42:10923–10930, 2003). The MD simulations provide a consistent pattern in that long lifetimes of hydrogen bonds go along with small amplitude fluctuations. In agreement with experiments, differences in stability were found with a rather flexible N-terminal segment as compared with a more rigid C-terminal part. Protection factors of backbone hydrogen exchange correlate strongly with the number of contacts but also with hydrogen-bond occupancy, hydrogen-bond survival times, as well as the inverse of fluctuations of backbone atoms and hydrogen-bond lengths derived from MD simulation data. We observed a conformational transition in the C-terminal loop, and significant motion in the N-terminal loop, which can be interpreted as being the structural units involved in the onset of the protein unfolding process in agreement with experimental evidence on mitochondrial cytochrome c. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users. Gernot Kieseritzky and Giulia Morra both contributed equally to this work.     

Kinetics of cyanide binding as a probe of local stability/flexibility of cytochrome c

Effect of anions of the Hofmeister series (thiocyanate, perchlorate, iodide, bromide, nitrate, chloride, sulfate, and phosphate) on local and global stability and flexibility of horse heart ferricytochrome c (cyt c) has been studied. Global stability of cyt c was determined by iso/thermal denaturations monitored by change in ellipticity in the far-UV region and its local stability was determined from absorbance changes in the Soret region. Particularly, relative stability/flexibility of the Met80–heme iron bond has been assessed by analysis of binding of cyanide into the heme iron. Both global and local stabilities of cyt c exhibited monotonous increase induced by a change of anion from chaotropic to kosmotropic species. However, this monotonous dependence was not observed for the rate constants of cyanide association with cyt c. As expected more chaotropic ions induced lower stability of protein and faster binding of cyanide but this correlation was reversed for kosmotropic anions. We propose that the unusual bell-shaped dependence of the rate constant of cyanide association is a result of modulation of Met80–heme iron bond strength and/or flexibility of heme region by Hofmeister anions independently on global stability of cyt c. Further, our results demonstrate sensitivity of cyanide binding to local change in stability/flexibility in the heme region of cyt c.     

A structural model of the Sgt2 protein and its interactions with chaperones and the Get4/Get5 complex

The insertion of tail-anchored transmembrane (TA) proteins into the appropriate membrane is a post-translational event that requires stabilization of the transmembrane domain and targeting to the proper destination. Sgt2 is a heat-shock protein cognate (HSC) co-chaperone that preferentially binds endoplasmic reticulum-destined TA proteins and directs them to the GET pathway via Get4 and Get5. Here, we present the crystal structure from a fungal Sgt2 homolog of the tetratrico-repeat (TPR) domain and part of the linker that connects to the C-terminal domain. The linker extends into the two-carboxylate clamp of the TPR domain from a symmetry-related molecule mimicking the binding to HSCs. Based on this structure, we provide biochemical evidence that the Sgt2 TPR domain has the ability to directly bind multiple HSC family members. The structure allows us to propose features involved in this lower specificity relative to other TPR containing co-chaperones. We further show that a dimer of Sgt2 binds a single Get5 and use small angle x-ray scattering to characterize the domain arrangement of Sgt2 in solution. These results allow us to present a structural model of the Sgt2-Get4/Get5-HSC complex.     

Molecular dynamics simulations to investigate the stability and aggregation behaviour of the amyloid-forming peptide VQIVYK from tau protein

The formation of paired helical filaments arising from the short hexapeptide in the third repeat of tau protein, ³⁰⁶VQIVYK³¹¹, is critical for tau polymerisation. The atomic structure of the VQIVYK oligomer has revealed a dry, tightly self-complementing structure between the neighbouring β-sheet layers, termed as ‘steric zipper’. In this study, several molecular dynamics simulations with all-atom explicit water were conducted to investigate the structural stability and aggregation behaviour of the VQIVYK peptide with various sizes and its single alanine replacement mutations. Our results indicate that the van der Waals interaction between side chains of Q2, the π–π stacking interaction between aromatic rings of Y5, and the electrostatic interaction between K6 and the C-terminus play an important role in stabilising the VQIVYK oligomers within the same β-sheet layer, while hydrophobic steric zipper involving V1, I3 and Y5 is responsible for holding the neighbouring β-sheet layers together. The twisted angles of the VQIVYK oligomers were also analysed and shown to be size dependent. The present results not only provide atomic insights into amyloid formation, but are also helpful for designing new or modified capping peptides and inhibitors to prevent fibril formation of the VQIVYK peptide from tau protein.