Active site water molecules revealed in the 2.1 A resolution structure of a site-directed mutant of isocitrate dehydrogenase

Capturing the Dynamic Conformational Changes of Human Isocitrate Dehydrogenase 1 (IDH1) upon Ligand and Metal Binding Using Hydrogen-Deuterium Exchange Mass Spectrometry

Human isocitrate dehydrogenase 1 (IDH1) is a highly conserved metabolic enzyme that catalyzes the interconversion of isocitrate and α-ketoglutarate. Kinetic and structural studies with IDH1 have revealed evidence of striking conformational changes that occur upon binding of its substrates, isocitrate and NADP⁺, and its catalytic metal cation. Here, we used hydrogen-deuterium exchange mass spectrometry (HDX-MS) to build a comprehensive map of the dynamic conformational changes experienced by IDH1 upon ligand binding. IDH1 proved well-suited for HDX-MS analysis, allowing us to capture profound changes in solvent accessibility at substrate binding sites and at a known regulatory region, as well as at more distant local subdomains that appear to support closure of this protein into its active conformation. HDX-MS analysis suggested that IDH1 is primarily purified with NADP(H) bound in the absence of its metal cation. Subsequent metal cation binding, even in the absence of isocitrate, was critical for driving large conformational changes. WT IDH1 folded into its fully closed conformation only when the full complement of substrates and metal was present. Finally, we show evidence supporting a previously hypothesized partially open conformation that forms prior to the catalytically active state, and we propose this conformation is driven by isocitrate binding in the absence of metal.

Tricarboxylic acid cycle dehydrogenases inhibition by naringenin: experimental and molecular modelling evidence

The study of polyphenols’ effects on health has been gaining attention lately. In addition to reacting with important enzymes, altering the cell metabolism, these substances can present either positive or negative metabolic alterations depending on their consumption levels. Naringenin, a citrus flavonoid, already presents diverse metabolic effects. The objective of this work was to evaluate the effect of maternal naringenin supplementation during pregnancy on the tricarboxylic acid cycle activity in offspring’s cerebellum. Adult female Wistar rats were divided into two groups: (1) vehicle (1 ml/kg by oral administration (p.o.)) or (2) naringenin (50 mg/kg p.o.). The offspring were euthanised at 7th day of life, and the cerebellum was dissected to analyse citrate synthase, isocitrate dehydrogenase (IDH), α-ketoglutarate dehydrogenase (α-KGDH) and malate dehydrogenase (MDH) activities. Molecular docking used SwissDock web server and FORECASTER Suite, and the proposed binding pose image was created on UCSF Chimera. Data were analysed by Student’s t test. Naringenin supplementation during pregnancy significantly inhibited IDH, α-KGDH and MDH activities in offspring’s cerebellum. A similar reduction was observed in vitro, using purified α-KGDH and MDH, subjected to pre-incubation with naringenin. Docking simulations demonstrated that naringenin possibly interacts with dehydrogenases in the substrate and cofactor binding sites, inhibiting their function. Naringenin administration during pregnancy may affect cerebellar development and must be evaluated with caution by pregnant women and their physicians.

The role of hydration effects in 5-fluorouridine binding to SOD1: insight from a new 3D-RISM-KH based protocol for including structural water in docking simulations

Misfolded Cu/Zn superoxide dismutase enzyme (SOD1) shows prion-like propagation in neuronal cells leading to neurotoxic aggregates that are implicated in amyotrophic lateral sclerosis (ALS). Tryptophan-32 (W32) in SOD1 is part of a potential site for templated conversion of wild type SOD1. This W32 binding site is located on a convex, solvent exposed surface of the SOD1 suggesting that hydration effects can play an important role in ligand recognition and binding. A recent X-ray crystal structure has revealed that 5-Fluorouridine (5-FUrd) binds at the W32 binding site and can act as a pharmacophore scaffold for the development of anti-ALS drugs. In this study, a new protocol is developed to account for structural (non-displaceable) water molecules in docking simulations and successfully applied to predict the correct docked conformation binding modes of 5-FUrd at the W32 binding site. The docked configuration is within 0.58 Å (RMSD) of the observed configuration. The docking protocol involved calculating a hydration structure around SOD1 using molecular theory of solvation (3D-RISM-KH, 3D-Reference Interaction Site Model-Kovalenko-Hirata) whereby, non-displaceable water molecules are identified for docking simulations. This protocol was also used to analyze the hydrated structure of the W32 binding site and to explain the role of solvation in ligand recognition and binding to SOD1. Structural water molecules mediate hydrogen bonds between 5-FUrd and the receptor, and create an environment favoring optimal placement of 5-FUrd in the W32 binding site.

Crystal Structures of the Putative Isocitrate Dehydrogenase from Sulfolobus tokodaii Strain 7 in the Apo and NADP+-Bound Forms

Isocitrate dehydrogenase is a catabolic enzyme that acts during the third step of the tricarboxylic acid cycle. The hypothetical protein ST2166 from the archaeon Sulfolobus tokodaii was isolated and crystallized. It shares high primary structure homology with prokaryotic NADP⁺-dependent IDHs, suggesting that these enzymes share a common enzymatic mechanism. The crystal structure of ST2166 was determined at 2.0 Å resolution in the apo form, and then the structure of the crystal soaked with NADP⁺ was also determined at 2.4 Å resolution, which contained NADP⁺ bound at the putative active site. Comparisons between the structures of apo and NADP⁺-bound forms and NADP-IDHs from other prokaryotes suggest that prokaryotic NADP-IDHs recognize their cofactors using conserved Lys335, Tyr336, and Arg386 in ST2166 at the opening cleft before the domain closure.

Understanding the regulatory mechanisms of milk production using integrative transcriptomic and proteomic analyses: improving inefficient utilization of crop by-products as forage in dairy industry

BACKGROUND

Bovine milk is an important nutrient source for humans. Forage plays a vital role in dairy husbandry via affecting milk quality and quantity. However, the differences in mammary metabolism of dairy cows fed different forages remain elucidated. In this study, we utilized transcriptomic RNA-seq and iTRAQ proteomic techniques to investigate and integrate the differences of molecular pathways and biological processes in the mammary tissues collected from 12 lactating cows fed corn stover (CS, low-quality, n = 6) and alfalfa hay (AH, high-quality, n = 6).

RESULTS

A total of 1631 differentially expressed genes (DEGs; 1046 up-regulated and 585 down-regulated) and 346 differentially expressed proteins (DEPs; 138 increased and 208 decreased) were detected in the mammary glands between the CS- and AH-fed animals. Expression patterns of 33 DEPs (18 increased and 15 decreased) were consistent with the expression of their mRNAs. Compared with the mammary gland of AH-fed cows, the marked expression changes found in the mammary gland of CS group were for genes involved in reduced mammary growth/development (COL4A2, MAPK3, IKBKB, LGALS3), less oxidative phosphorylation (ATPsynGL, ATP6VOA1, ATP5H, ATP6VOD1, NDUFC1), enhanced lipid uptake/metabolism (SLC27A6, FABP4, SOD2, ACADM, ACAT1, IDH1, SCP2, ECHDC1), more active fatty acid beta-oxidation (HMGCS1), less amino acid/protein transport (SLC38A2, SLC7A8, RAB5a, VPS18), reduced protein translation (RPS6, RPS12, RPS16, RPS19, RPS20, RPS27), more proteasome- (PSMC2, PSMC6, PSMD14, PSMA2, PSMA3) and ubiquitin-mediated protein degradation (UBE2B, UBE2H, KLHL9, HSPH1, DNAJA1 and CACYBP), and more protein disassembly-related enzymes (SEC63, DNAJC3, DNAJB1, DNAJB11 and DNAJC12).

CONCLUSION

Our results indicate that the lower milk production in the CS-fed dairy cows compared with the AH-fed cows was associated with a network of mammary gene expression changes, importantly, the prime factors include decreased energy metabolism, attenuated protein synthesis, enhanced protein degradation, and the lower mammary cell growth. The present study provides insights into the effects of the varying quality of forages on mammary metabolisms, which can help the improvement of strategies in feeding dairy cows with CS-based diet.

Methods for Studying Iron Regulatory Protein 1: An Important Protein in Human Iron Metabolism

Iron regulatory proteins 1 and 2 (IRP1 and IRP2) are two cytosolic proteins that maintain cellular iron homeostasis by regulating the expression of genes involved in iron metabolism. IRPs respond to cellular iron deficiency by binding to iron-responsive elements (IREs) found in the mRNAs of iron metabolism transcripts, enhancing iron import, and reducing iron storage, utilization, and export. IRP1, a bifunctional protein, exists in equilibrium between a [Fe₄S₄] cluster containing cytosolic aconitase, and an apoprotein that binds to IREs. At high cellular iron levels, this equilibrium is shifted more toward iron-sulfur cluster containing aconitase, whereas IRP2 undergoes proteasomal degradation by an E3 ubiquitin ligase complex that contains an F-box protein, FBXL5. Irp1^-/- mice develop polycythemia and pulmonary hypertension, whereas Irp2^-/- mice develop microcytic anemia and progressive neurodegeneration, indicating that Irp1 has important functions in the erythropoietic and pulmonary systems, and Irp2 has essential roles in supporting erythropoiesis and nervous system functions. Mice lacking both Irp1 and Irp2 die during embryogenesis, suggesting that functions of Irp1 and Irp2 are redundant. In this review, we will focus on the methods for studying IRP1 activities and function in cells and animals.

Activity Based Protein Profiling Leads to Identification of Novel Protein Targets for Nerve Agent VX

Organophosphorus (OP) nerve agents continue to be a threat at home and abroad during the war against terrorism. Human exposure to nerve agents such as VX results in a cascade of toxic effects relative to the exposure level including ocular miosis, excessive secretions, convulsions, seizures, and death. The primary mechanism behind these overt symptoms is the disruption of cholinergic pathways. While much is known about the primary toxicity mechanisms of nerve agents, there remains a paucity of information regarding impacts on other pathways and systemic effects. These are important for establishing a comprehensive understanding of the toxic mechanisms of OP nerve agents. To identify novel proteins that interact with VX, and that may give insight into these other mechanisms, we used activity-based protein profiling (ABPP) employing a novel VX-probe on lysates from rat heart, liver, kidney, diaphragm, and brain tissue. By making use of a biotin linked VX-probe, proteins covalently bound by the probe were isolated and enriched using streptavidin beads. The proteins were then digested, labeled with isobarically distinct tandem mass tag (TMT) labels, and analyzed by liquid chromatography tandem mass spectrometry (LC-MS/MS). Quantitative analysis identified 132 bound proteins, with many proteins found in multiple tissues. As with previously published ABPP OP work, monoacylglycerol lipase associated proteins and fatty acid amide hydrolase (FAAH) were shown to be targets of VX. In addition to these two and other predicted neurotransmitter-related proteins, a number of proteins involved with energy metabolism were identified. Four of these enzymes, mitochondrial isocitrate dehydrogenase 2 (IDH2), isocitrate dehydrogenase 3 (IDH3), malate dehydrogenase (MDH), and succinyl CoA (SCS) ligase, were assayed for VX inhibition. Only IDH2 NADP+ activity was shown to be inhibited directly. This result is consistent with other work reporting animals exposed to OP compounds exhibit reduced IDH activity. Though clearly a secondary mechanism for toxicity, this is the first time VX has been shown to directly interfere with energy metabolism. Taken together, the ABPP work described here suggests the discovery of novel protein-agent interactions, which could be useful for the development of novel diagnostics or potential adjuvant therapeutics.

Effect of livingstonepotato (Plectranthus esculenthus N.E.Br) on hyperglycemia, antioxidant activity and lipid metabolism of streptozotocin induced diabetic rats

The effect of livingstone potato (Plectranthus esculenthus N.E.Br) on serum glucose, glycated hemoglobin (HbA1C), serum triglyceride, total cholesterol, high density lipoprotein (HDL), low density lipoprotein (LDL), very low density lipoprotein (VLDL), hepatic malic enzyme (ME), isocitrate dehydrogenase (IDH) and catalase activities of Streptozotocin induced diabetic rats were investigated using standard techniques. The atherogenic index (AI) and coronary risk index (CRI) of the rats were calculated as the ratios of LDL to HDL and total cholesterol to HDL, respectively. The serum glucose of the non-diabetic, diabetic control and diabetic rats given livingstone potato incorporated feeds (test feed) were 92.58 ± 3.97, 352.30 ± 4.88 and 165.50 ± 7.88 mg/dl, respectively. Intake of the test feed by the diabetic rats of group 3, resulted in significant (P < 0.05) decrease of their serum glucose, HbA1c, triglyceride, cholesterol, LDL, VLDL, AI and CRI but significant increase (P < 0.05) of hepatic levels of ME, IDH, catalase and serum HDL compared with the diabetic control rats that had significant alteration of these parameters (P < 0.05) compared with the non-diabetic rats. The feed intakes of the non-diabetic, diabetic control and diabetic rats given the test feed were 133.34 ± 1.32, 137.84 ± 5.77 and 146.38 ± 4.33 g/rat/week by the last week of experimentation. The diabetic control rats recorded significant loss of weight (P < 0.05) compared with the non diabetic rats despite increased feed intake. Chemical analysis of the standard and test feeds showed that the standard rat feed contained 15.00 ± 0.78% protein, 7.24 ± 1.20% fat, 31.55 ± 2.62% carbohydrates, energy value of 290.65 ± 4.77 kcal/100 g, 10% crude fiber and 0.12 ± 0.04 mg Gallic Acid Equivalent while the test feed contained 40.10 ± 0.16% carbohydrates, 17.22 ± 0.40% protein, 22.16 ± 0.34% fat, energy value of 428.70 ± 2.12 kcal/100 g, 8.51 ± 0.16% crude fiber, 1.3 ± 0.2 mg Gallic Acid Equivalent/g of sample and strong antioxidant activity comparable to standard quercetin. The study shows the potentials of livingstone potato in the management of diabetes and hyperlipidemia.

Differential lysine acetylation profiles of Erwinia amylovora strains revealed by proteomics

Protein lysine acetylation (LysAc) has recently been demonstrated to be widespread in E. coli and Salmonella, and to broadly regulate bacterial physiology and metabolism. However, LysAc in plant pathogenic bacteria is largely unknown. Here we first report the lysine acetylome of Erwinia amylovora, an enterobacterium causing serious fire blight disease of apples and pears. Immunoblots using generic anti-lysine acetylation antibodies demonstrated that growth conditions strongly affected the LysAc profiles in E. amylovora. Differential LysAc profiles were also observed for two E. amylovora strains, known to have differential virulence in plants, indicating translational modification of proteins may be important in determining virulence of bacterial strains. Proteomic analysis of LysAc in two E. amylovora strains identified 141 LysAc sites in 96 proteins that function in a wide range of biological pathways. Consistent with previous reports, 44% of the proteins are involved in metabolic processes, including central metabolism, lipopolysaccharide, nucleotide and amino acid metabolism. Interestingly, for the first time, several proteins involved in E. amylovora virulence, including exopolysaccharide amylovoran biosynthesis- and type III secretion-associated proteins, were found to be lysine acetylated, suggesting that LysAc may play a major role in bacterial virulence. Comparative analysis of LysAc sites in E. amylovora and E. coli further revealed the sequence and structural commonality for LysAc in the two organisms. Collectively, these results reinforce the notion that LysAc of proteins is widespread in bacterial metabolism and virulence.

Induced fit and the catalytic mechanism of isocitrate dehydrogenase

NADP(+) dependent isocitrate dehydrogenase (IDH; EC 1.1.1.42) belongs to a large family of α-hydroxyacid oxidative β-decarboxylases that catalyze similar three-step reactions, with dehydrogenation to an oxaloacid intermediate preceding β-decarboxylation to an enol intermediate followed by tautomerization to the final α-ketone product. A comprehensive view of the induced fit needed for catalysis is revealed on comparing the first "fully closed" crystal structures of a pseudo-Michaelis complex of wild-type Escherichia coli IDH (EcoIDH) and the "fully closed" reaction product complex of the K100M mutant with previously obtained "quasi-closed" and "open" conformations. Conserved catalytic residues, binding the nicotinamide ring of NADP(+) and the metal-bound substrate, move as rigid bodies during domain closure by a hinge motion that spans the central β-sheet in each monomer. Interactions established between Thr105 and Ser113, which flank the "phosphorylation loop", and the nicotinamide mononucleotide moiety of NADP(+) establish productive coenzyme binding. Electrostatic interactions of a Lys100-Leu103-Asn115-Glu336 tetrad play a pivotal role in assembling a catalytically competent active site. As predicted, Lys230* is positioned to deprotonate/reprotonate the α-hydroxyl in both reaction steps and Tyr160 moves into position to protonate C3 following β-decarboxylation. A proton relay from the catalytic triad Tyr160-Asp307-Lys230* connects the α-hydroxyl of isocitrate to the bulk solvent to complete the picture of the catalytic mechanism.

Role of bridging water molecules in GSK3β-inhibitor complexes: insights from QM/MM, MD, and molecular docking studies

The role of water molecules is increasingly gaining interest in drug design, and several studies have highlighted their paramount contributions to the specificity and the affinity of ligand binding. In this study, we employ the two-layer ONIOM-based quantum mechanics/molecular mechanics (QM/MM) calculations, molecular dynamics (MD) simulations, and molecular docking studies to investigate the effect of bridging water molecules at the GSK3β-inhibitors interfaces. The results obtained from the ONIOM geometry optimization and AIM analysis corroborated the presence of bridging water molecules that form hydrogen bonds with protein side chain of Thr138 and/or backbone of Gln185, and mediate interactions with inhibitors in the 10 selected GSK3β-inhibitor complexes. Subsequently, MD simulations carried out on a representative system of 1R0E demonstrated that the bridging water molecule is stable at the GSK3β-inhibitor interface and appears to contribute to the stability of the protein-inhibitor interactions. Furthermore, molecular docking studies of GSK3β-inhibitor complexes indicated that the inhibitors can increase binding affinities and the better docked conformation of inhibitors can be obtained by inclusion of the bridging water molecules, especially for the flexible inhibitors, in docking experiments into individual protein conformations. Our results elucidate the importance of bridging water molecules at the GSK3β-inhibitor interfaces and suggest that they might prove useful in rational drug design.

Structural characterization of tartrate dehydrogenase: a versatile enzyme catalyzing multiple reactions

The first structure of an NAD-dependent tartrate dehydrogenase (TDH) has been solved to 2 A resolution by single anomalous diffraction (SAD) phasing as a complex with the intermediate analog oxalate, Mg(2+) and NADH. This TDH structure from Pseudomonas putida has a similar overall fold and domain organization to other structurally characterized members of the hydroxy-acid dehydrogenase family. However, there are considerable differences between TDH and these functionally related enzymes in the regions connecting the core secondary structure and in the relative positioning of important loops and helices. The active site in these complexes is highly ordered, allowing the identification of the substrate-binding and cofactor-binding groups and the ligands to the metal ions. Residues from the adjacent subunit are involved in both the substrate and divalent metal ion binding sites, establishing a dimer as the functional unit and providing structural support for an alternating-site reaction mechanism. The divalent metal ion plays a prominent role in substrate binding and orientation, together with several active-site arginines. Functional groups from both subunits form the cofactor-binding site and the ammonium ion aids in the orientation of the nicotinamide ring of the cofactor. A lysyl amino group (Lys192) is the base responsible for the water-mediated proton abstraction from the C2 hydroxyl group of the substrate that begins the catalytic reaction, followed by hydride transfer to NAD. A tyrosyl hydroxyl group (Tyr141) functions as a general acid to protonate the enolate intermediate. Each substrate undergoes the initial hydride transfer, but differences in substrate orientation are proposed to account for the different reactions catalyzed by TDH.

Each conserved active site tyr in the three subunits of human isocitrate dehydrogenase has a different function

The human NAD-dependent isocitrate dehydrogenase (IDH) is a heterotetrameric mitochondrial enzyme with 2alpha:1beta:1gamma subunit ratio. The three subunits share 40-52% identity in amino acid sequence and each includes a tyrosine in a comparable position: alphaY126, betaY137, and gammaY135. To study the role of the corresponding tyrosines of each of the subunits of human NAD-IDH, the tyrosines were mutated (one subunit at a time) to Ser, Phe, or Glu. Enzymes were expressed with one mutant and two wild-type subunits. The results of characterization of the mutant enzymes suggest that betaY137 is involved in NAD binding and allosteric activation by ADP. The alphaY126 is required for catalytic activity and likely acts as a general acid in the reaction. The gammaY135 is also required for catalytic activity and may be involved in proper folding of the enzyme. The corresponding tyrosines in the three dissimilar subunits of NAD-IDH thus have distinctive functions.

Ab initio protein structure prediction with force field parameters derived from water-phase quantum chemical calculation

Molecular dynamics (MD) simulations are extensively used in the study of the structures and functions of proteins. Ab initio protein structure prediction is one of the most important subjects in computational biology, and many trials have been performed using MD simulation so far. Since the results of MD simulations largely depend on the force field, reliable force field parameters are indispensable for the success of MD simulation. In this work, we have modified atom charges in a standard force field on the basis of water-phase quantum chemical calculations. The modified force field turned out appropriate for ab initio protein structure prediction by the MD simulation with the generalized Born method. Detailed analysis was performed in terms of the conformational stability of amino acid residues, the stability of secondary structure of proteins, and the accuracy for prediction of protein tertiary structure, comparing the modified force field with a standard one. The energy balance between alpha-helix and beta-sheet structures was significantly improved by the modification of charge parameters.

Ligand-protein docking with water molecules

The presence of water molecules plays an important role in the accuracy of ligand-protein docking predictions. Comprehensive docking simulations have been performed on a large set of ligand-protein complexes whose crystal structures contain water molecules in their binding sites. Only those water molecules found in the immediate vicinity of both the ligand and the protein were considered. We have investigated whether prior optimization of the orientation of water molecules in either the presence or absence of the bound ligand has any effect on the accuracy of docking predictions. We have observed a statistically significant overall increase in accuracy when water molecules are included during docking simulations and have found this to be independent of the method of optimization of the orientation of water molecules. These results confirm the importance of including water molecules whenever possible in a ligand-protein docking simulation. Our findings also reveal that prior optimization of the orientation of water molecules, in the absence of any bound ligand, does not have a detrimental effect on the improved accuracy of ligand-protein docking. This is important, given the use of docking simulations to predict the binding modes of new ligands or drug molecules.

Significance of water molecules in the inhibition of cylin-dependent kinase 2 and 5 complexes

Interest in CDK2 and CDK5 has stemmed mainly from their association with cancer and neuronal migration or differentiation related diseases and the need to design selective inhibitors for these kinases. In the present paper, eight Molecular Dynamics (MD) simulations are carried out to examine the importance of structure and dynamics of water in the active site of both CDK2 and CDK5 complexes with roscovitine and indirubin analogues. Together with previous results, the current work shows a highly conserved water-involved hydrogen bonding (HB) network in both CDK2- and CDK5-indirubin combinations to complete information from the X-ray crystallography. The simulations suggest the importance of such a network for combining the inhibitor to the host protein as well as the significance of using an activated CDK as a template when designing new inhibitors. Different binding patterns of roscovitine in CDK2 and CDK5 are detected during the simulations because of the different binding conformations of the group on the C2 side chain, which might offer a clue toward finding highly selective inhibitors with regards to CDK2 and CDK5.

Complete reversal of coenzyme specificity of isocitrate dehydrogenase from Haloferax volcanii

Haloferax volcanii Ds-threo-isocitrate dehydrogenase (ICDH) was highly expressed in bacteria as inclusion bodies. The recombinant enzyme was refolded, purified and characterized, and was found to be NADP-dependent like the wild-type protein. Sequence alignment of several isocitrate dehydrogenases from evolutionarily divergent organisms including H. volcanii revealed that the amino acid residues involved in coenzyme specificity are highly conserved. Our objective was to switch the coenzyme specificity of halophilic ICDH by altering these conserved amino acids. We were able to switch coenzyme specificity from NADP+ to NAD+ by changing five amino acids by site-directed mutagenesis (Arg291, Lys343, Tyr344, Val350 and Tyr390). The five mutants of ICDH were overexpressed in Escherichia coli as inclusion bodies and each recombinant ICDH protein was refolded and purified, and its kinetic parameters were determined. Coenzyme specificity did not switch until all five amino acids were substituted.

The effect of a tightly bound water molecule on scaffold diversity in the computer-aided de novo ligand design of CDK2 inhibitors

We have determined the effects that tightly bound water molecules have on the de novo design of cyclin-dependent kinase-2 (CDK2) ligands. In particular, we have analyzed the impact of a specific structural water molecule on the chemical diversity and binding mode of ligands generated through a de novo structure-based ligand generation method in the binding site of CDK2. The tightly bound water molecule modifies the size and shape of the binding site and we have found that it also imposed constraints on the observed binding modes of the generated ligands. This in turn had the indirect effect of reducing the chemical diversity of the underlying molecular scaffolds that were able to bind to the enzyme satisfactorily. [Figure: see text].

Including Tightly-Bound Water Molecules in de Novo Drug Design. Exemplification through the in Silico Generation of Poly(ADP-ribose)polymerase Ligands

Different strategies for the in silico generation of ligand molecules in the binding site of poly(ADP-ribose)polymerase (PARP) were studied in order to observe the effect of the targeting and displacement of tightly bound water molecules. Several molecular scaffolds were identified as having better interactions in the binding site when targeting one or two tightly bound water molecules in the NAD binding site. Energy calculations were conducted in order to assess the ligand-protein and ligand-water-protein interactions of different functional groups of the generated ligands. These calculations were used to evaluate the energetic consequences of the presence of tightly bound water molecules and to identify those that contribute favorably to the binding of ligands.

Critical role of Lys212 and Tyr140 in porcine NADP-dependent isocitrate dehydrogenase

Lys212 and Tyr140 are close to the enzyme-bound isocitrate in the recently determined crystal structure of porcine NADP-specific isocitrate dehydrogenase (Ceccarelli, C., Grodsky, N. B., Ariyaratne, N., Colman, R. F., and Bahnson, B. J. (2002) J. Biol. Chem. 277, 43454-43462). We have constructed mutant enzymes in which Lys212 is replaced by Gln, Tyr, and Arg, and Tyr140 is replaced by Phe, Thr, Glu, and Lys. Wild type and mutant enzymes were each expressed in Escherichia coli and purified to homogeneity. At pH 7.4, the specific activity is decreased in K212Q, K212Y, and K212R, respectively, to 0.01-9% of wild type. The most striking change is in the pH-V(max) curves. Wild type depends on the deprotonated form of a group of pKaes 5.7, whereas this pKaes is increased to 7.4 in neutral K212Q and to 8.3 in K212Y. In contrast, the positive K212R has a pKaes of 5.9. These results indicate that (by electrostatic repulsion) a positively charged residue at position 212 lowers the pK of the nearby ionizable group in the enzyme-substrate complex. Lys212 may also stabilize the carbanion formed initially on substrate decarboxylation. The Tyr140 mutants have specific activities at pH 7.4 that are reduced to 0.2-0.5% of those of wild type, whereas their Km values for isocitrate and NADP are not increased. Most notable are the altered pH-V(max) profiles. V(max) is constant from pH 5.3 to 8 for Y140F and Y140T and increases as pH is decreased for Y140E and Y140K. These results suggest that in wild type enzyme, Tyr140 is the general acid that protonates the substrate after decarboxylation and that the carboxyl and ammonium forms of Y140E and Y140K provide partial substitutes. Relative to wild type, the Y140T enzyme is specifically activated 106-fold by exogenous addition of acetic acid and 88-fold by added phenol; and the K212Q enzyme is activated 4-fold by added ethylamine. These chemical rescue experiments support the conclusion that Tyr140 and Lys212 are required for the catalytic activity of porcine NADP-dependent isocitrate dehydrogenase.

Crystal structure of the monomeric isocitrate dehydrogenase in the presence of NADP+: insight into the cofactor recognition, catalysis, and evolution

NADP+-dependent monomeric isocitrate dehydrogenase (IDH) from the nitrogen-fixing bacterium Azotobacter vinelandii (AvIDH) is one of members of the beta-decarboxylating dehydrogenase family and catalyzes the dehydration and decarboxylation of isocitrate to yield 2-oxoglutrate and CO2 in the Krebs cycle. We solved the crystal structure of the AvIDH in complex with cofactor NADP+ (AvIDH-NADP+ complex). The final refined model shows the closed form that has never been detected in any previously solved structures of beta-decarboxylating dehydrogenases. The structure also reveals all of the residues that interact with NADP+. The structure-based sequence alignment reveals that these residues were not conserved in any other dimeric NADP+-dependent IDHs. Therefore the NADP+ specificity of the monomeric and dimeric IDHs was independently acquired through the evolutional process. The AvIDH was known to show an exceptionally high turnover rate. The structure of the AvIDH-NADP+ complex indicates that one loop, which is not present in the Escherichia coli IDHs, reliably stabilizes the conformation of the nicotinamide mononucleotide of the bound NADP+ by forming a few hydrogen bonds, and such interactions are considered to be important for the monomeric enzyme to initiate the hydride transfer reaction immediately. Finally, the structure of the AvIDH is compared with that of other dimeric NADP-IDHs. Several structural features demonstrate that the monomeric IDHs are structurally more related to the eukaryotic dimeric IDHs than to the bacterial dimeric IDHs.

WaterScore: a novel method for distinguishing between bound and displaceable water molecules in the crystal structure of the binding site of protein-ligand complexes

We have performed a multivariate logistic regression analysis to establish a statistical correlation between the structural properties of water molecules in the binding site of a free protein crystal structure, with the probability of observing the water molecules in the same location in the crystal structure of the ligand-complexed form. The temperature B-factor, the solvent-contact surface area, the total hydrogen bond energy and the number of protein-water contacts were found to discriminate between bound and displaceable water molecules in the best regression functions obtained. These functions may be used to identify those bound water molecules that should be included in structure-based drug design and ligand docking algorithms. FIGURE The binding site ( thin sticks) of penicillopepsin (3app) with its crystallographically determined water molecules ( spheres) and superimposed ligand (in thick sticks, from complexed structure 1ppk). Water molecules sterically displaced by the ligand upon complexation are shown in cyan. Bound water molecules are shown in blue. Displaced water molecules are shown in yellow. Water molecules removed from the analysis due to a lack of hydrogen bonds to the protein are shown in white. WaterScore correctly predicted waters in blue as Probability=1 to remain bound and waters in yellow as Probability<1x10(-20) to remain bound.

Crystal structure of porcine mitochondrial NADP+-dependent isocitrate dehydrogenase complexed with Mn2+ and isocitrate. Insights into the enzyme mechanism

The crystal structure of porcine heart mitochondrial NADP+-dependent isocitrate dehydrogenase (IDH) complexed with Mn2+ and isocitrate was solved to a resolution of 1.85 A. The enzyme was expressed in Escherichia coli, purified as a fusion protein with maltose binding protein, and cleaved with thrombin to yield homogeneous enzyme. The structure was determined by multiwavelength anomalous diffraction phasing using selenium substitution in the form of selenomethionine as the anomalous scatterer. The porcine NADP+-IDH enzyme is structurally compared with the previously solved structures of IDH from E. coli and Bacillus subtilis that share 16 and 17% identity, respectively, with the mammalian enzyme. The porcine enzyme has a protein fold similar to the bacterial IDH structures with each monomer folding into two domains. However, considerable differences exist between the bacterial and mammalian forms of IDH in regions connecting core secondary structure. Based on the alignment of sequence and structure among the porcine, E. coli, and B. subtilis IDH, a putative phosphorylation site has been identified for the mammalian enzyme. The active site, including the bound Mn2+-isocitrate complex, is highly ordered and, therefore, mechanistically informative. The consensus IDH mechanism predicts that the Mn2+-bound hydroxyl of isocitrate is deprotonated prior to its NADP+-dependent oxidation. The present crystal structure has an active site water that is well positioned to accept the proton and ultimately transfer the proton to solvent through an additional bound water.

NADP-dependent isocitrate dehydrogenase from the halophilic archaeon Haloferax volcanii: cloning, sequence determination and overexpression in Escherichia coli

A gene encoding NADP-dependent Ds-threo-isocitrate dehydrogenase was isolated from Haloferax volcanii genomic DNA by using a combination of polymerase chain reaction and screening of a lambda EMBL3 library. Analysis of the nucleotide sequence revealed an open reading frame of 1260 bp encoding a protein of 419 amino acids with 45837 Da molecular mass. This sequence is highly similar to previously sequenced isocitrate dehydrogenases. In the alignment of the amino acid sequences with those from several archaeal and mesophilic NADP-dependent isocitrate dehydrogenases, the residues involved in dinucleotide binding and isocitrate binding are well conserved. We have developed methods for the expression in Escherichia coli and purification of the enzyme from H. volcanii. This expression was carried out in E. coli as inclusion bodies using the cytoplasmic expression vector pET3a. The enzyme was refolded by solubilisation in 8 M urea followed by dilution into a buffer containing EDTA, MgCl(2) and 3 M NaCl. Maximal activity was obtained after several hours incubation at room temperature.

Bacillus subtilis isocitrate dehydrogenase. A substrate analogue for Escherichia coli isocitrate dehydrogenase kinase/phosphatase

In Escherichia coli, the homodimeric Krebs cycle enzyme isocitrate dehydrogenase (EcIDH) is regulated by reversible phosphorylation of a sequestered active site serine. The phosphorylation cycle is catalyzed by a bifunctional protein, IDH kinase/phosphatase (IDH-K/P). To better understand the nature of the interaction between EcIDH and IDH-K/P, we have examined the ability of an IDH homologue from Bacillus subtilis (BsIDH) to serve as a substrate for the kinase and phosphatase activities. BsIDH exhibits extensive sequence and structural similarities with EcIDH, particularly around the phosphorylated serine. Our previous crystallographic analysis revealed that the active site architecture of these two proteins is almost completely conserved. We now expand the comparison to include a number of biochemical properties. Both IDHs display nearly equivalent steady-state kinetic parameters for the dehydrogenase reaction. Both proteins are also phosphorylated by IDH-K/P in the same ratio (1 mole of phosphate per mole of monomer), and this stoichiometric phosphorylation correlates with an equivalent inhibition of IDH activity. Furthermore, tandem electrospray mass spectrometry demonstrates that BsIDH, like EcIDH, is phosphorylated on the corresponding active site serine residue (Ser-104). Despite the high degree of sequence, functional, and structural congruence between these two proteins, BsIDH is surprisingly a much poorer substrate of IDH-K/P than is EcIDH, with Michaelis constants for the kinase and phosphatase activities elevated by 60- and 3,450-fold, respectively. These drastically disparate values might result from restricted access to the active site cavity and/or from the lack of a potential docking site for IDH-K/P.

On the role of conserved histidine 106 in 10-formyltetrahydrofolate dehydrogenase catalysis: connection between hydrolase and dehydrogenase mechanisms

The enzyme, 10-formyltetrahydrofolate dehydrogenase (FDH), converts 10-formyltetrahydrofolate (10-formyl-THF) to tetrahydrofolate in an NADP(+)-dependent dehydrogenase reaction or an NADP(+)-independent hydrolase reaction. The hydrolase reaction occurs in a 310-amino acid long amino-terminal domain of FDH (N(t)-FDH), whereas the dehydrogenase reaction requires the full-length enzyme. The amino-terminal domain of FDH shares some sequence identity with several other enzymes utilizing 10-formyl-THF as a substrate. These enzymes have two strictly conserved residues, aspartate and histidine, in the putative catalytic center. We have shown recently that the conserved aspartate is involved in FDH catalysis. In the present work we studied the role of the conserved histidine, His(106), in FDH function. Site-directed mutagenesis experiments showed that replacement of the histidine with alanine, asparagine, aspartate, glutamate, glutamine, or arginine in N(t)-FDH resulted in expression of insoluble proteins. Replacement of the histidine with another positively charged residue, lysine, produced a soluble mutant with no hydrolase activity. The insoluble mutants refolded from inclusion bodies adopted a conformation inherent to the wild-type N(t)-FDH, but they did not exhibit any hydrolase activity. Substitution of alanine for three non-conserved histidines located close to the conserved one did not reveal any significant changes in the hydrolase activity of N(t)-FDH. Expressed full-length FDH with the substitution of lysine for the His(106) completely lost both the hydrolase and dehydrogenase activities. Thus, our study showed that His(106), besides being an important structural residue, is also directly involved in both the hydrolase and dehydrogenase mechanisms of FDH. Modeling of the putative hydrolase catalytic center/folate-binding site suggested that the catalytic residues, aspartate and histidine, are unlikely to be adjacent to the catalytic cysteine in the aldehyde dehydrogenase catalytic center. We hypothesize that 10-formyl-THF dehydrogenase reaction is not an independent reaction but is a combination of hydrolase and aldehyde dehydrogenase reactions.