Conserved water mediated recognition and the dynamics of active site Cys 331 and Tyr 411 in hydrated structure of human IMPDH‐II

Inosine monophosphate dehydrogenase (IMPDH) of human is involved in GMP biosynthesis pathway, increased level of IMPDH‐II (an isoform of enzyme) activity have found in leukemic and sarcoma cells. Modeling and extensive molecular dynamics simulation (15 ns) studies of IMPDH‐II (1B3O PDB structure) have indicated the intricate involvement of four conserved water molecules (W 1, W 2, W 3, and W 4) in the conformational transition or the mobilities of “flap” (residues 400–450) and “loop” (residues 325–342) regions in enzyme. The stabilization of active site residues Asn 303, Gly 324, Ser 329, Cys 331, Asp 364, and Tyr 411 through variable H‐bonding coordination from the conserved water molecular center seems interesting in the uninhibited hydrated form of human IMPDH‐II structures. This conformational transition or the flexibility of mobile regions, water molecular recognition to active site residues Cys 331 and Tyr 411, and the presence of a hydrophilic cavity ～540 ų (enclaved by the loop and flap region) near the C‐terminal surface of this enzyme may explore a rational hope toward the water mimic inhibitor or anticancer agent design for human. Copyright © 2010 John Wiley & Sons, Ltd.     

Conserved water-mediated recognition and dynamics of NAD+ (carboxamide group) to hIMPDH enzyme: water mimic approach toward the design of isoform-selective inhibitor

Inosine monophosphate dehydrogenase (IMPDH) enzyme involves in GMP biosynthesis pathway. Type I hIMPDH is expressed at lower levels in all cells, whereas type II is especially observed in acute myelogenous leukemia, chronic myelogenous leukemia cancer cells, and 10?ns simulation of the IMP–NAD⁺ complex structures (PDB ID. 1B3O and 1JCN) have revealed the presence of a few conserved hydrophilic centers near carboxamide group of NAD⁺. Three conserved water molecules (W1, W, and W1′) in di-nucleotide binding pocket of enzyme have played a significant role in the recognition of carboxamide group (of NAD⁺) to D274 and H93 residues. Based on H-bonding interaction of conserved hydrophilic (water molecular) centers within IMP–NAD⁺-enzyme complexes and their recognition to NAD⁺, some covalent modification at carboxamide group of di-nucleotide (NAD⁺) has been made by substituting the –CONH₂group by –CONHNH₂ (carboxyl hydrazide group) using water mimic inhibitor design protocol. The modeled structure of modified ligand may, though, be useful for the development of antileukemic agent or it could be act as better inhibitor for hIMPDH-II.     

Reconciling Mediating and Slaving Roles of Water in Protein Conformational Dynamics

Proteins accomplish their physiological functions with remarkably organized dynamic transitions among a hierarchical network of conformational substates. Despite the essential contribution of water molecules in shaping functionally important protein dynamics, their exact role is still controversial. Water molecules were reported either as mediators that facilitate or as masters that slave protein dynamics. Since dynamic behaviour of a given protein is ultimately determined by the underlying energy landscape, we systematically analysed protein self energies and protein-water interaction energies obtained from extensive molecular dynamics simulation trajectories of barstar. We found that protein-water interaction energy plays the dominant role when compared with protein self energy, and these two energy terms on average have negative correlation that increases with increasingly longer time scales ranging from 10 femtoseconds to 100 nanoseconds. Water molecules effectively roughen potential energy surface of proteins in the majority part of observed conformational space and smooth in the remaining part. These findings support a scenario wherein water on average slave protein conformational dynamics but facilitate a fraction of transitions among different conformational substates, and reconcile the controversy on the facilitating and slaving roles of water molecules in protein conformational dynamics.     

172 Role of conserved water molecular triad in the recognition of IMP,NAD+ with Asp 274, Asn 303, Arg 322, and Asp 364 in both the isoform of hIMPDH

Inosine monophosphate dehydrogenase (IMPDH) plays an important role in the Guanosine monophosphate (GMP) biosynthesis pathway. As hIMPDH-II is involved in CML-Cancer, it is thought to be an active target for leukemic drug design. The importance of conserved water molecules in the salt-bridge-mediated interdomain recognition and loop-flap recognition of hIMPDH has already been indicated in some simulation studies (Bairagya et al., 2009, 2011a, 2011b, 2012; Mishra et al., 2012). In this work, the role of conserved water molecules in the recognition of Inosine monophosphate (IMP) and NAD⁺ (co-factor) to active site residues of both the isoforms has been investigated by all atoms MD-Simulation studies. During 25-ns dynamics of the solvated hIMPDH-II and I (1B3O and 1JCN PDB structures), the involvement of conserved water molecular triad (W _M, W _L and W _C) in the recognition of active site residues (Asp 274, Asn 303, Arg 322, and Asp 364), IMP and NAD⁺ has been observed (Figure 1). The H-bonding co-ordination of all three conserved water molecular centers is within 4–7 and their occupation frequency is 1.0. The H-bonding geometry and the electronic consequences of the water molecular interaction at the different residues (and also IMP and NAD⁺) may put forward some rational clues on antileukemic agent design.     

Conserved water-mediated H-bonding dynamics of catalytic Asn 175 in plant thiol protease

The role of invariant water molecules in the activity of plant cysteine protease is ubiquitous in nature. On analysing the 11 different Protein DataBank (PDB) structures of plant thiol proteases, the two invariant water molecules W1 and W2 (W220 and W222 in the template 1PPN structure) were observed to form H-bonds with the O_b atom of Asn 175. Extensive energy minimization and molecular dynamics simulation studies up to 2 ns on all the PDB and solvated structures clearly revealed the involvement of the H-bonding association of the two water molecules in fixing the orientation of the asparagine residue of the catalytic triad. From this study, it is suggested that H-bonding of the water molecule at the W1 invariant site better stabilizes the Asn residue at the active site of the catalytic triad.     

Conserved water molecular dynamics of the different X-ray structures of rusticyanin: an unique aquation potentiality of the ligand bonded Cu++ center

The invariant water molecular interaction involving in the Rusticyanin of Thiobacillus ferrooxidans is thought to be important for its molecular complexation with other proteins at differential acidophilic situation. The comparative analysis of the different x-ray, energy minimized, and auto solvated structures of Rusticyanin revealed the presence of five specific invariant bound water molecules (among the approximately 150 water molecules per monomer) in the crystals. The five W 205, W 206, W 112, W 214, and W 221 water molecules (in Rusticyanin PDB code: 1RCY) were seem to be invariant in all the seven structures (PDB codes: 1RCY, 1A3Z, 1A8Z, 1E3O, 1GY1, 1GY2, 2CAL). Among the five conserved water molecules the W 221 (of 1 RCY or the equivalent water molecules in the other oxidized form of Rusticyanin structures) had endowed an interesting coordination potentiality to Cu(+2) ion during the energy minimization. The W 221 was observed to approach toward the tetrahedrally bonded Cu(+2) ion through the opposite (or trans) route of metal-bonded Met 148. This direct water molecular coordination affected the tetrahedral geometry of Cu(+2) to trigonal bipyramidal. Presumably this structural dynamics at the Cu(+2) center could involve in the electron transport process during protein-protein complexation.     

Molecular dynamics study of a complex between the human histocompatibility antigen HLA-A2 and the IMP58-66 nonapeptide from influenza virus matrix protein.

The structure of the influenza-virus-matrix-protein (IMP) 58-66 nonapeptide, bound to the major-histocompatibility-complex-encoded human leukocyte antigen (HLA) A2 protein was studied by molecular dynamics simulation. Starting from the extra electron density map of peptides co-crystallized with HLA-A2, the nonapeptide IMP58-66 was docked residue by residue in the protein binding cleft. The complex was simulated for 100 ps in a shell of 1372 water molecules. The averaged simulated HLA-A2 conformation was found to be similar to the crystal structure (0.182 nm RMS deviation, for the backbone atoms of the alpha 1-alpha 2 domain). Nine out of the 14 hydrogen bonds observed in the antigen-binding site were reproduced in the simulation. The IMP58-66 peptide exhibits an extended conformation with kinks at positions 3 and 5. The side chains of residues 2, 3 and 9 develop van der Waals' interactions with hydrophobic pockets of HLA-A2, corresponding to polymorphic residues of the major-histocompatibility-complex-encoded proteins. Both the N-terminus and C-terminus of the nonapeptide were anchored in the antigen-binding groove by hydrogen bonds with conserved amino acids. The N-terminus was more flexible and contacts four HLA-A2 conserved tyrosines (Tyr7, Tyr59, Tyr159 and Tyr171) and Glu63 by direct or water-relayed hydrogen bonds. Water intercalation occurred only around the N-terminus of the peptide, the C-terminal carboxylate forming strong hydrogen bonds with polar residues (Tyr84 and Thr143) and a salt bridge with Lys146 all over the molecular dynamics simulation. This model is fully compatible with the recently published crystal structure of the HLA-B27 protein, complexed by a mixture of self nonapeptides.     

Autophosphorylation dependent destabilization of the insulin receptor kinase domain: tryptophan-1175 reports changes in the catalytic cleft

Protein kinases are regulated by conformational or chemical changes which facilitate access of substrates to the active site and promote correct orientations of catalytically essential residues and water molecules. The switch between basal and activated states of the insulin receptor's kinase domain (IRKD) results from autophosphorylation. We investigated the effects of IRKD autophosphorylation on the conformational stability by guanidine hydrochloride (GdnHCl) dependent denaturation and by iodide quenching of intrinsic fluorescence. Tryptophan residues of the recombinant soluble IRKD (residues R953-S1355) were excited at a lambdaex of 295 nm, and emission spectra were analyzed for centroid (a characteristic of average polarity of the indole rings' environments) and integrated fluorescence intensity over the lambdaem range of 310-420 nm. Denaturation profiles of both apo- and phospho-IRKD forms are complex with at least three distinct unfolding transitions. The first and last transitions were reversible and cooperative and had midpoints at 0.4 or 0.7 M GdnHCl and 2.4 or 2.7 M GdnHCl, respectively; transitions of phospho-IRKD occurred at lower GdnHCl concentrations. Calculations of free energy of unfolding suggested a loss of approximately 2.3 kcal/mol of stabilization for the first transition and approximately 1.5 kcal/mol for the third transition. Circular dichroism showed subtle changes in secondary structure over the first transition and global unfolding over the last transition. The first transition reports changes primarily in the local environment of W1175, which is near the catalytic loop and is conserved among protein tyrosine kinases. W1175 is also the dominant fluorophore of the native emission spectrum. Iodide quenching of W1175 was virtually undetectable in the apo-IRKD but significant in the phospho-IRKD, suggesting that W1175 exposure to small solutes is strongly dependent on the conformation of the activation loop. These studies indicate that autophosphorylation, while exposing the catalytic center, also produces a conformer less stable than the apoenzyme.     

Conformational transitions in eosinophil cationic protein: a molecular dynamics study in aqueous environment

Extensive molecular dynamics simulations have been performed on eosinophil cationic protein (ECP). The two structures found in the crystallographic dimer (ECPA and ECPB) have been independently simulated. A small difference in the pattern of the sidechain hydrogen bonds in the starting structure has resulted in interesting differences in the conformations accessed during the simulations. In one simulation (ECPB), a stable equilibrium conformation was obtained and in the other (ECPA), conformational transitions at the level of sidechain interactions were observed. The conformational transitions exhibit the involvement of the solvent (water) molecules with a pore-like construct in the equilibrium conformation and an opening for a large number of water molecules during the transition phase. The details of these transitions are examined in terms of intra-protein hydrogen bonds, protein-water networks and the residence times of water molecules on the polar atoms of the protein. These properties show some significant differences in the region between the N-terminal helix and the loop before the C-terminal strand as a function of different conformations accessed during the simulations. However, the stable hydrogen bonds, the protein-water networks, and the hydration patterns in most part of the protein including the active site are very much similar in both the simulations, indicating the fact that these are intrinsic properties of proteins.     

An insight to conserved water molecular dynamics of catalytic and structural Zn(+2) ions in matrix metalloproteinase 13 of human

Matrix Metalloproteinase (MMP)--13 or Collagenase--3 plays a significant role in the formation and remodeling of bone, tumor invasion and causes osteoarthritis. Water molecular dynamic studies of the five (1XUC, 1XUD, 1XUR, 456C, 830C) PDB and solvated structures of MMP-13 in human have been carried out upto 5 ns on assigning the differential charges (+2, +1, +0.5 e) to both the Zinc ions. The MM and MD-studies have revealed the coordination of three water molecules (W(H), W(I) and W(S)) to Zn(c) and one water to Zn(s). The transition of geometry around the Znc from tetrahedral to octahedral via trigonal bipyramidal, and for Zn(s) from tetrahedral to trigonal bipyramidal are seem interesting. Recognition of two zinc ions through water molecular bridging (Zn(c) - W(H) (W(1))...W(2)....W(3)....H(187) Zn(s)) and the stabilization of variable coordination geometries around metal ions may indicate the possible involvement of Zn(c) ...Zn(s) coupled mechanism in the catalytic process. So the hydrophilic topology and stereochemistry of water mediated coupling between Zn-ions may provide some plausible hope towards the design of some bidentate/polydentate bridging ligands or inhibitors for MMP-13.     

Role of invariant water molecules in retaining the active site geometry of β-lactamase: a molecular dynamics simulation study

Water molecules play a critical role in stabilising the three-dimensional architecture, dynamics and function of biological macromolecules. Comparative analysis of structurally similar proteins has shown that there are water molecules conserved in the same relative positions and make similar hydrogen bonds with proteins in all crystal structures. These invariant water molecules are essential for the maintenance of the native structure of proteins. The present study explores the role of invariant water molecules to maintain the active site geometry of β-lactamase enzyme. Thirteen crystal structures of class-A β-lactamase from Staphylococcus aureus have been used in this study. Molecular dynamics simulations of the protein structures were performed in hydrated as well as in dehydrated conditions. The analysis showed that significant changes occur in the active site geometry due to dehydration. These changes can be attributed to the removal of water molecules at the active site.     

Molecular dynamics simulations of heme reorientational motions in myoglobin.

Molecular dynamics simulations of 2-ns duration were performed on carbonmonoxymyoglobin and deoxymyoglobin in vacuo to study the reorientational dynamics of the heme group. The heme in both simulations undergoes reorientations of approximately 5 degrees amplitude on a subpicosecond time scale, which produce a rapid initial decay in the reorientational correlation function to about 0.99. The heme also experiences infrequent changes in average orientation of approximately 10 degrees amplitude, which lead to a larger slow decay of the reorientational correlation function over a period of hundreds of picoseconds. The simulations have not converged with respect to these infrequent transitions. However, an estimate of the order parameter for rapid internal motions of the heme from those orientations which are sampled by the simulations suggests that the subnanosecond orientational dynamics of the heme accounts for at least 30% of the unresolved initial anisotropy decay observed in the nanosecond time-resolved optical absorption experiments on myoglobin reported by Ansari et al. in a companion paper (Ansari, A., C.M. Jones, E.R. Henry, J. Hofrichter, and W.A. Eaton. 1992. Biophys. J. 64:852-868.). A more complete sampling of the accessible heme orientations would most likely increase this fraction further. The simulation of the liganded molecule also suggests that the conformational dynamics of the CO ligand may contribute significantly to discrepancies between the ligand conformation as probed by x-ray diffraction and by infrared-optical photoselection experiments. The protein back-bone explores multiple conformations during the simulations, with the largest structural changes appearing in the E and F helices, which are in contact with the heme. The variations in the heme orientation correlate with the conformational dynamics of the protein on a time scale of hundreds of picoseconds, suggesting that the heme orientation may provide a useful probe of dynamical processes in the protein.     

Molecular dynamics study of water penetration in staphylococcal nuclease

The ionization properties of Lys and Glu residues buried in the hydrophobic core of staphylococcal nuclease (SN) suggest that the interior of this protein behaves as a highly polarizable medium with an apparent dielectric constant near 10. This has been rationalized previously in terms of localized conformational relaxation concomitant with the ionization of the internal residue, and with contributions by internal water molecules. Paradoxically, the crystal structure of the SN V66E variant shows internal water molecules and the structure of the V66K variant does not. To assess the structural and dynamical character of interior water molecules in SN, a series of 10-ns-long molecular dynamics (MD) simulations was performed with wild-type SN, and with the V66E and V66K variants with Glu66 and Lys66 in the neutral form. Internal water molecules were identified based on their coordination state and characterized in terms of their residence times, average location, dipole moment fluctuations, hydrogen bonding interactions, and interaction energies. The locations of the water molecules that have residence times of several nanoseconds and display small mean-square displacements agree well with the locations of crystallographically observed water molecules. Additional, relatively disordered water molecules that are not observed crystallographically were found in internal hydrophobic locations. All of the interior water molecules that were analyzed in detail displayed a distribution of interaction energies with higher mean value and narrower width than a bulk water molecule. This underscores the importance of protein dynamics for hydration of the protein interior. Further analysis of the MD trajectories revealed that the fluctuations in the protein structure (especially the loop elements) can strongly influence protein hydration by changing the patterns or strengths of hydrogen bonding interactions between water molecules and the protein. To investigate the dynamical response of the protein to burial of charged groups in the protein interior, MD simulations were performed with Glu66 and Lys66 in the charged state. Overall, the MD simulations suggest that a conformational change rather than internal water molecules is the dominant determinant of the high apparent polarizability of the protein interior.     

Nucleotide sequence of the gene encoding the GMP reductase of Escherichia coli K12.

(1) The nucleotide sequence of a 1991 bp segment of DNA that expresses the GMP reductase (guaC) gene of Escherichia coli K12 was determined. (2) This gene comprises 1038 bp, 346 codons (including the initiation codon but excluding the termination codon), and it encodes a polypeptide of Mr 37,437 which is in good agreement with previous maxicell studies. (3) The sequence contains a putative promoter 102 bp upstream of the translational start codon, and this is immediately followed by a (G + C)-rich discriminator sequence suggesting that guaC expression may be under stringent control (4) The GMP reductase exhibits a high degree of sequence identity (34%) with IMP dehydrogenase (the guaB gene product) indicative of a close evolutionary relationship between the salvage pathway and the biosynthetic enzymes, GMP reductase and IMP dehydrogenase, respectively. (5) A single conserved cysteine residue, possibly involved in IMP binding to IMP dehydrogenase, was located within a region that possesses some of the features of a nucleotide binding site. (6) The IMP dehydrogenase polypeptide contains an internal segment of 123 amino acid residues that has no counterpart in GMP reductase and may represent an independent folding domain flanked by (alanine + glycine)-rich interdomain linkers.     

Conserved Water Molecular Dynamics of the Different X-ray Structures of Rusticyanin: An Unique Aquation Potentiality of the Ligand Bonded Cu++ Center

Abstract

The invariant water molecular interaction involving in the Rusticyanin of Thiobacillus ferrooxidans is thought to be important for its molecular complexation with other proteins at differential acidophilic situation. The comparative analysis of the different x-ray, energy minimized, and auto solvated structures of Rusticyanin revealed the presence of five specific invariant bound water molecules (among the ~ 150 water molecules per monomer) in the crystals. The five W 205, W 206, W 112, W 214, and W 221 water molecules (in Rusticyanin PDB code: 1RCY) were seem to be invariant in all the seven structures (PDB codes: 1RCY, 1A3Z, 1A8Z, 1E3O, 1GY1, 1GY2, 2CAL). Among the five conserved water molecules the W 221 (of 1 RCY or the equivalent water molecules in the other oxidized form of Rusticyanin structures) had endowed an interesting coordination potentiality to Cu⁺² ion during the energy minimization. The W 221 was observed to approach toward the tetrahedrally bonded Cu⁺² ion through the opposite (or trans) route of metal-bonded Met 148. This direct water molecular coordination affected the tetrahedral geometry of Cu⁺² to trigonal bipyramidal. Presumably this structural dynamics at the Cu⁺² center could involve in the electron transport process during protein-protein complexation.     

Recognition dynamics of trinuclear copper cluster and associated histidine residues through conserved or semi-conserved water molecules in human Ceruloplasmin: The involvement of aspartic and glutamic acid gates

Human Ceruloplasmin belongs to the family of multi-copper oxidases and it is involved in different physiological processes, copper ion transport, iron metabolism, iron homeostasis, and biogenic amine metabolism. MD-simulation studies have indicated the higher hydrophilic susceptibility of the trinuclear copper cluster in native CP compared to its oxygen bound form. The copper (T2/T3) atom Cu3047 of the cluster, which is close to T1 copper center Cu3052 (~13 Å) has a higher affinity for water molecules compared to other copper centers. The water molecules of W3, W4, W5, W9, and W12 conserved water sites are coordinated to Cu3047, where W3, W9, and W12 centers are found to play some crucial role in the stabilization of native trinuclear copper cluster. The hydrogen bonding interaction of Asp169, Glu112, Asp995, and Glu1032 residues with the copper-bound conserved water molecules (W3, W4, W5, W10, and W11) in native CP is observed to be unique. The conformational flexibility of Asp169 and Glu112 and their association with the copper-bound water molecules, but the absence of such interaction in O₂-bound simulated structure of the enzyme is indicating some plausible rational on the role of these acidic residues in the gating of O₂ molecule in the native trinuclear Cu cluster of CP. The simulation results may shade some new light on the biochemistry/chemistry of CP, specially on the hydration dynamics of the trinuclear copper cluster.     

An insight to the conserved water mediated dynamics of catalytic His88 and its recognition to thyroxin and RBP binding residues in human transthyretin

Human transthyretin (hTTR) is a multifunctional protein involved in several amyloidogenic diseases. Besides transportation of thyroxin and vitamin-A, its role towards the catalysis of apolipoprotein-A1 and Aβ-peptide are also drawing interest. The role of water molecules in the catalytic mechanism is still unknown. Extensive analyses of 14 high-resolution X-ray structures of human transthyretin and MD simulation studies have revealed the presence of eight conserved hydrophilic centres near its catalytic zone which may be indispensable for the function, dynamics and stability of the protein. Three water molecules (W1, W2 and W3) form a cluster and play an important role in the recognition of the catalytic and RBP-binding residues. They also induce the reorganisation of the His88 for coupling with other catalytic residues (His90, Glu92). Another water molecule (W5) participate in inter-monomer recognition between the catalytic and thyroxin binding sites. The rest four water molecules (W6, W*, W^# and W^?) form a distorted tetrahedral cluster and impart stability to the catalytic core of hTTR. The conserved water mediated recognition dynamics of the different functional sites may provide some rational clues towards the understanding of the activity and mechanism of hTTR.