Host target modification as a strategy to counter pathogen hijacking of the jasmonate hormone receptor.

In the past decade, characterization of the host targets of pathogen virulence factors took a center stage in the study of pathogenesis and disease susceptibility in plants and humans. However, the impressive knowledge of host targets has not been broadly exploited to inhibit pathogen infection. Here, we show that host target modification could be a promising new approach to "protect" the disease-vulnerable components of plants. In particular, recent studies have identified the plant hormone jasmonate (JA) receptor as one of the common targets of virulence factors from highly evolved biotrophic/hemibiotrophic pathogens. Strains of the bacterial pathogen Pseudomonas syringae, for example, produce proteinaceous effectors, as well as a JA-mimicking toxin, coronatine (COR), to activate JA signaling as a mechanism to promote disease susceptibility. Guided by the crystal structure of the JA receptor and evolutionary clues, we succeeded in modifying the JA receptor to allow for sufficient endogenous JA signaling but greatly reduced sensitivity to COR. Transgenic Arabidopsis expressing this modified receptor not only are fertile and maintain a high level of insect defense, but also gain the ability to resist COR-producing pathogens Pseudomonas syringae pv. tomato and P. syringae pv. maculicola. Our results provide a proof-of-concept demonstration that host target modification can be a promising new approach to prevent the virulence action of highly evolved pathogens.

Combined QM(DFT)/MM molecular dynamics simulations of the deamination of cytosine by yeast cytosine deaminase (yCD).

Extensive combined quantum mechanical (B3LYP/6-31G*) and molecular mechanical (QM/MM) molecular dynamics simulations have been performed to elucidate the hydrolytic deamination mechanism of cytosine to uracil catalyzed by the yeast cytosine deaminase (yCD). Though cytosine has no direct binding to the zinc center, it reacts with the water molecule coordinated to zinc, and the adjacent conserved Glu64 serves as a general acid/base to shuttle protons from water to cytosine. The overall reaction consists of several proton-transfer processes and nucleophilic attacks. A tetrahedral intermediate adduct of cytosine and water binding to zinc is identified and similar to the crystal structure of yCD with the inhibitor 2-pyrimidinone. The rate-determining step with the barrier of 18.0 kcal/mol in the whole catalytic cycle occurs in the process of uracil departure where the proton transfer from water to Glu64 and nucleophilic attack of the resulting hydroxide anion to C2 of the uracil ring occurs synchronously. © 2016 Wiley Periodicals, Inc.

Structural determinants of host specificity of complement Factor H recruitment by Streptococcus pneumoniae.

Many human pathogens have strict host specificity, which affects not only their epidemiology but also the development of animal models and vaccines. Complement Factor H (FH) is recruited to pneumococcal cell surface in a human-specific manner via the N-terminal domain of the pneumococcal protein virulence factor choline-binding protein A (CbpAN). FH recruitment enables Streptococcus pneumoniae to evade surveillance by human complement system and contributes to pneumococcal host specificity. The molecular determinants of host specificity of complement evasion are unknown. In the present study, we show that a single human FH (hFH) domain is sufficient for tight binding of CbpAN, present the crystal structure of the complex and identify the critical structural determinants for host-specific FH recruitment. The results offer new approaches to the development of better animal models for pneumococcal infection and redesign of the virulence factor for pneumococcal vaccine development and reveal how FH recruitment can serve as a mechanism for both pneumococcal complement evasion and adherence.

Molecular dynamics simulations of the Escherichia coli HPPK apo-enzyme reveal a network of conformational transitions.

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the first reaction in the folate biosynthetic pathway. Comparison of its X-ray and nuclear magnetic resonance structures suggests that the enzyme undergoes significant conformational change upon binding to its substrates, especially in three catalytic loops. Experimental research has shown that even when confined by crystal contacts, loops 2 and 3 remain rather flexible when the enzyme is in its apo form, raising questions about the putative large-scale induced-fit conformational change of HPPK. To investigate the loop dynamics in a crystal-free environment, we performed conventional molecular dynamics simulations of the apo-enzyme at two different temperatures (300 and 350 K). Our simulations show that the crystallographic B-factors considerably underestimate the loop dynamics; multiple conformations of loops 2 and 3, including the open, semi-open, and closed conformations that an enzyme must adopt throughout its catalytic cycle, are all accessible to the apo-enzyme. These results revise our previous view of the functional mechanism of conformational change upon MgATP binding and offer valuable structural insights into the workings of HPPK. In this paper, conformational network analysis and principal component analysis related to the loops are discussed to support the presented conclusions.

Structure of bacterial transcription factor SpoIIID and evidence for a novel mode of DNA binding.

SpoIIID is evolutionarily conserved in endospore-forming bacteria, and it activates or represses many genes during sporulation of Bacillus subtilis. An SpoIIID monomer binds DNA with high affinity and moderate sequence specificity. In addition to a predicted helix-turn-helix motif, SpoIIID has a C-terminal basic region that contributes to DNA binding. The nuclear magnetic resonance (NMR) solution structure of SpoIIID in complex with DNA revealed that SpoIIID does indeed have a helix-turn-helix domain and that it has a novel C-terminal helical extension. Residues in both of these regions interact with DNA, based on the NMR data and on the effects on DNA binding in vitro of SpoIIID with single-alanine substitutions. These data, as well as sequence conservation in SpoIIID binding sites, were used for information-driven docking to model the SpoIIID-DNA complex. The modeling resulted in a single cluster of models in which the recognition helix of the helix-turn-helix domain interacts with the major groove of DNA, as expected. Interestingly, the C-terminal extension, which includes two helices connected by a kink, interacts with the adjacent minor groove of DNA in the models. This predicted novel mode of binding is proposed to explain how a monomer of SpoIIID achieves high-affinity DNA binding. Since SpoIIID is conserved only in endospore-forming bacteria, which include important pathogenic Bacilli and Clostridia, whose ability to sporulate contributes to their environmental persistence, the interaction of the C-terminal extension of SpoIIID with DNA is a potential target for development of sporulation inhibitors.

Feedback inhibition of deoxy-D-xylulose-5-phosphate synthase regulates the methylerythritol 4-phosphate pathway.

The 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway leads to the biosynthesis of isopentenyl diphosphate (IDP) and dimethylallyl diphosphate (DMADP), the precursors for isoprene and higher isoprenoids. Isoprene has significant effects on atmospheric chemistry, whereas other isoprenoids have diverse roles ranging from various biological processes to applications in commercial uses. Understanding the metabolic regulation of the MEP pathway is important considering the numerous applications of this pathway. The 1-deoxy-D-xylulose-5-phosphate synthase (DXS) enzyme was cloned from Populus trichocarpa, and the recombinant protein (PtDXS) was purified from Escherichia coli. The steady-state kinetic parameters were measured by a coupled enzyme assay. An LC-MS/MS-based assay involving the direct quantification of the end product of the enzymatic reaction, 1-deoxy-D-xylulose 5-phosphate (DXP), was developed. The effect of different metabolites of the MEP pathway on PtDXS activity was tested. PtDXS was inhibited by IDP and DMADP. Both of these metabolites compete with thiamine pyrophosphate for binding with the enzyme. An atomic structural model of PtDXS in complex with thiamine pyrophosphate and Mg(2+) was built by homology modeling and refined by molecular dynamics simulations. The refined structure was used to model the binding of IDP and DMADP and indicated that IDP and DMADP might bind with the enzyme in a manner very similar to the binding of thiamine pyrophosphate. The feedback inhibition of PtDXS by IDP and DMADP constitutes an important mechanism of metabolic regulation of the MEP pathway and indicates that thiamine pyrophosphate-dependent enzymes may often be affected by IDP and DMADP.

Role of glutamate 64 in the activation of the prodrug 5-fluorocytosine by yeast cytosine deaminase.

Yeast cytosine deaminase (yCD) catalyzes the hydrolytic deamination of cytosine to uracil as well as the deamination of the prodrug 5-fluorocytosine (5FC) to the anticancer drug 5-fluorouracil. In this study, the role of Glu64 in the activation of the prodrug 5FC was investigated by site-directed mutagenesis, biochemical, nuclear magnetic resonance (NMR), and computational studies. Steady-state kinetics studies showed that the mutation of Glu64 causes a dramatic decrease in k(cat) and a dramatic increase in K(m), indicating Glu64 is important for both binding and catalysis in the activation of 5FC. (19)F NMR experiments showed that binding of the inhibitor 5-fluoro-1H-pyrimidin-2-one (5FPy) to the wild-type yCD causes an upfield shift, indicating that the bound inhibitor is in the hydrated form, mimicking the transition state or the tetrahedral intermediate in the activation of 5FC. However, binding of 5FPy to the E64A mutant enzyme causes a downfield shift, indicating that the bound 5FPy remains in an unhydrated form in the complex with the mutant enzyme. (1)H and (15)N NMR analysis revealed trans-hydrogen bond D/H isotope effects on the hydrogen of the amide of Glu64, indicating that the carboxylate of Glu64 forms two hydrogen bonds with the hydrated 5FPy. ONIOM calculations showed that the wild-type yCD complex with the hydrated form of the inhibitor 1H-pyrimidin-2-one is more stable than the initial binding complex, and in contrast, with the E64A mutant enzyme, the hydrated inhibitor is no longer favored and the conversion has a higher activation energy, as well. The hydrated inhibitor is stabilized in the wild-type yCD by two hydrogen bonds between it and the carboxylate of Glu64 as revealed by (1)H and (15)N NMR analysis. To explore the functional role of Glu64 in catalysis, we investigated the deamination of cytosine catalyzed by the E64A mutant by ONIOM calculations. The results showed that without the assistance of Glu64, both proton transfers before and after the formation of the tetrahedral reaction intermediate become partially rate-limiting steps. The results of the experimental and computational studies together indicate that Glu64 plays a critical role in both the binding and the chemical transformation in the conversion of the prodrug 5FC to the anticancer drug 5-fluorouracil.

Bisubstrate analog inhibitors of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase: new lead exhibits a distinct binding mode.

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK), a key enzyme in the folate biosynthesis pathway catalyzing the pyrophosphoryl transfer from ATP to 6-hydroxymethyl-7,8-dihydropterin, is an attractive target for developing novel antimicrobial agents. Previously, we studied the mechanism of HPPK action, synthesized bisubstrate analog inhibitors by linking 6-hydroxymethylpterin to adenosine through phosphate groups, and developed a new generation of bisubstrate inhibitors by replacing the phosphate bridge with a piperidine-containing linkage. To further improve linker properties, we have synthesized a new compound, characterized its protein binding/inhibiting properties, and determined its structure in complex with HPPK. Surprisingly, this inhibitor exhibits a new binding mode in that the adenine base is flipped when compared to previously reported structures. Furthermore, the side chain of amino acid residue E77 is involved in protein-inhibitor interaction, forming hydrogen bonds with both 2' and 3' hydroxyl groups of the ribose moiety. Residue E77 is conserved among HPPK sequences, but interacts only indirectly with the bound MgATP via water molecules. Never observed before, the E77-ribose interaction is compatible only with the new inhibitor-binding mode. Therefore, this compound represents a new direction for further development.

Bisubstrate analogue inhibitors of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase: New design with improved properties.

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK), a key enzyme in the folate biosynthetic pathway, catalyzes the pyrophosphoryl transfer from ATP to 6-hydroxymethyl-7,8-dihydropterin. The enzyme is essential for microorganisms, is absent from humans, and is not the target for any existing antibiotics. Therefore, HPPK is an attractive target for developing novel antimicrobial agents. Previously, we characterized the reaction trajectory of HPPK-catalyzed pyrophosphoryl transfer and synthesized a series of bisubstrate analog inhibitors of the enzyme by linking 6-hydroxymethylpterin to adenosine through 2, 3, or 4 phosphate groups. Here, we report a new generation of bisubstrate analog inhibitors. To improve protein binding and linker properties of such inhibitors, we have replaced the pterin moiety with 7,7-dimethyl-7,8-dihydropterin and the phosphate bridge with a piperidine linked thioether. We have synthesized the new inhibitors, measured their K(d) and IC(50) values, determined their crystal structures in complex with HPPK, and established their structure-activity relationship. 6-Carboxylic acid ethyl ester-7,7-dimethyl-7,8-dihydropterin, a novel intermediate that we developed recently for easy derivatization at position 6 of 7,7-dimethyl-7,8-dihydropterin, offers a much high yield for the synthesis of bisubstrate analogs than that of previously established procedure.

Structure-Based Design of a Dual-Targeted Covalent Inhibitor Against Papain-like and Main Proteases of SARS-CoV-2.

The two proteases, PLpro and Mpro, of SARS-CoV-2 are essential for replication of the virus. Using a structure-based co-pharmacophore screening approach, we developed a novel dual-targeted inhibitor that is equally potent in inhibiting PLpro and Mpro of SARS-CoV-2. The inhibitor contains a novel warhead, which can form a covalent bond with the catalytic cysteine residue of either enzyme. The maximum rate of the covalent inactivation is comparable to that of the most potent inhibitors reported for the viral proteases and covalent inhibitor drugs currently in clinical use. The covalent inhibition appears to be very specific for the viral proteases. The inhibitor has a potent antiviral activity against SARS-CoV-2 and is also well tolerated by mice and rats in toxicity studies. These results suggest that the inhibitor is a promising lead for development of drugs for treatment of COVID-19.

Probing single-molecule enzyme active-site conformational state intermittent coherence.

The relationship between protein conformational dynamics and enzymatic reactions has been a fundamental focus in modern enzymology. Using single-molecule fluorescence resonance energy transfer (FRET) with a combined statistical data analysis approach, we have identified the intermittently appearing coherence of the enzymatic conformational state from the recorded single-molecule intensity-time trajectories of enzyme 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) in catalytic reaction. The coherent conformational state dynamics suggests that the enzymatic catalysis involves a multistep conformational motion along the coordinates of substrate-enzyme complex formation and product releasing, presenting as an extreme dynamic behavior intrinsically related to the time bunching effect that we have reported previously. The coherence frequency, identified by statistical results of the correlation function analysis from single-molecule FRET trajectories, increases with the increasing substrate concentrations. The intermittent coherence in conformational state changes at the enzymatic reaction active site is likely to be common and exist in other conformation regulated enzymatic reactions. Our results of HPPK interaction with substrate support a multiple-conformational state model, being consistent with a complementary conformation selection and induced-fit enzymatic loop-gated conformational change mechanism in substrate-enzyme active complex formation.

Dynamics in Thermotoga neapolitana adenylate kinase: 15N relaxation and hydrogen-deuterium exchange studies of a hyperthermophilic enzyme highly active at 30 degrees C.

Backbone conformational dynamics of Thermotoga neapolitana adenylate kinase in the free form (TNAK) and inhibitor-bound form (TNAK*Ap5A) were investigated at 30 degrees C using (15)N NMR relaxation measurements and NMR monitored hydrogen-deuterium exchange. With kinetic parameters identical to those of Escherichia coli AK (ECAK) at 30 degrees C, TNAK is a unique hyperthermophilic enzyme. These catalytic properties make TNAK an interesting and novel model to study the interplay between protein rigidity, stability, and activity. Comparison of fast time scale dynamics (picosecond to nanosecond) in the open and closed states of TNAK and ECAK at 30 degrees C reveals a uniformly higher rigidity across all domains of TNAK. Within this framework of a rigid TNAK structure, several residues located in the AMP-binding domain and in the core-lid hinge regions display high picosecond to nanosecond time scale flexibility. Together with the recent comparison of ECAK dynamics with those of hyperthermophilic Aquifex aeolicus AK (AAAK), our results provide strong evidence for the role of picosecond to nanosecond time scale fluctuations in both stability and activity. In the slow time scales, TNAK's increased rigidity is not uniform but localized in the AMP-binding and lid domains. The core domain amides of ECAK and TNAK in the open and closed states show comparable protection against exchange. Significantly, the hinges framing the lid domain show similar exchange data in ECAK and TNAK open and closed forms. Our NMR relaxation and hydrogen-deuterium exchange studies therefore suggest that TNAK maintains high activity at 30 degrees C by localizing flexibility to the hinge regions that are key to facilitating conformational changes.

Productive versus unproductive nucleotide binding in yeast guanylate kinase mutants: comparison of R41M with K14M by proton two dimensional transferred NOESY.

The R41M and K14M mutant enzymes of yeast guanylate kinase (GKy) were studied to investigate the effects of these site-directed mutations on bound-substrate conformations. Published X-ray crystal structures of yeast guanylate kinase indicate that K14 is part of the "P" loop involved in ATP and ADP binding, while R41 is suggested as a hydrogen bonding partner for the phosphoryl moiety of GMP. Both of these residues might be involved in transition state stabilization. Adenosine conformations of ATP and ADP and guanosine conformations of GMP bound to R41M and K14M mutant yeast guanylate kinase in the complexes GKy.MgATP, GKy.MgADP, and GKy.MgADP.[u-(13)C]GMP were determined by two-dimensional transferred nuclear Overhauser effect (TRNOESY) measurements combined with molecular dynamics simulations, and these conformations were compared with previously published conformations for the wild type. In the fully constrained, two substrate complexes, GKy.MgADP.[u-(13)C]GMP, the guanyl glycosidic torsion angle, chi, is 51 +/- 5 degrees for R41M and 47 +/- 5 degrees for K14M. Both are similar to the published 50 +/- 5 degrees published for wild type. For R41M with adenyl nucleotides, the glycosidic torsion angle, chi, was 55 +/- 5 degrees with MgATP, and 47 +/- 5 degrees with MgADP, which compares well to the 54 +/- 5 degrees published for wild type. However, for K14M with adenyl nucleotides, the glycosidic torsion angle was 30 +/- 5 degrees with MgATP and 28 +/- 5 degrees with MgADP. The results indicate that bound adenyl-nucleotides have significantly different conformations in the wild-type and K14M mutant enzymes, suggesting that K14 plays an important role in orienting the triphosphate of MgATP for catalysis.

Dynamics of the conformational transitions in the assembling of the Michaelis complex of a bisubstrate enzyme: a (15)N relaxation study of Escherichia coli 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase.

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the transfer of pyrophosphate from ATP to 6-hydroxymethyl-7,8-dihydropterin (HP), which follows an ordered bi-bi kinetic mechanism with ATP binding to the enzyme first. HPPK undergoes dramatic conformational changes during its catalytic cycle as revealed by X-ray crystallography, and the conformational changes are essential for the enzymatic catalysis as shown by site-directed mutagenesis and biochemical and crystallographic analysis of the mutants. However, the dynamic properties of the enzyme have not been measured experimentally. Here, we report a (15)N NMR relaxation study of the dynamic properties of Escherichia coli HPPK from the apo form to the binary substrate complex with MgATP (represented by MgAMPCPP, an ATP analogue) to the Michaelis complex (ternary substrate complex) with MgATP (represented by MgAMPCPP) and HP (represented by 7,7-dimethyl-6-hydroxypterin, an HP analogue). The results show that the binding of the nucleotide to HPPK does not cause major changes in the dynamic properties of the enzyme. Whereas enzymes are often more rigid when bound to the ligand or the substrate, the internal mobility of HPPK is not reduced and is even moderately increased in the binary complex, particularly in the catalytic loops. The internal mobility of the catalytic loops is significantly quenched upon the formation of the ternary complex, but some mobility remains. The enhanced motions in the catalytic loops of the binary substrate complex may be required for the assembling of the ternary complex. On the other hand, some degrees of mobility in the catalytic loops of the ternary complex may be required for the optimal stabilization of the transition state, which may need the instantaneous adjustment and alignment of the side-chain positions of catalytic residues. Such dynamic behaviors may be characteristic of bisubstrate enzymes.

Scoring ligand similarity in structure-based virtual screening.

Scoring to identify high-affinity compounds remains a challenge in virtual screening. On one hand, protein-ligand scoring focuses on weighting favorable and unfavorable interactions between the two molecules. Ligand-based scoring, on the other hand, focuses on how well the shape and chemistry of each ligand candidate overlay on a three-dimensional reference ligand. Our hypothesis is that a hybrid approach, using ligand-based scoring to rank dockings selected by protein-ligand scoring, can ensure that high-ranking molecules mimic the shape and chemistry of a known ligand while also complementing the binding site. Results from applying this approach to screen nearly 70 000 National Cancer Institute (NCI) compounds for thrombin inhibitors tend to support the hypothesis. EON ligand-based ranking of docked molecules yielded the majority (4/5) of newly discovered, low to mid-micromolar inhibitors from a panel of 27 assayed compounds, whereas ranking docked compounds by protein-ligand scoring alone resulted in one new inhibitor. Since the results depend on the choice of scoring function, an analysis of properties was performed on the top-scoring docked compounds according to five different protein-ligand scoring functions, plus EON scoring using three different reference compounds. The results indicate that the choice of scoring function, even among scoring functions measuring the same types of interactions, can have an unexpectedly large effect on which compounds are chosen from screening. Furthermore, there was almost no overlap between the top-scoring compounds from protein-ligand versus ligand-based scoring, indicating the two approaches provide complementary information. Matchprint analysis, a new addition to the SLIDE (Screening Ligands by Induced-fit Docking, Efficiently) screening toolset, facilitated comparison of docked molecules' interactions with those of known inhibitors. The majority of interactions conserved among top-scoring compounds for a given scoring function, and from the different scoring functions, proved to be conserved interactions in known inhibitors. This was particularly true in the S1 pocket, which was occupied by all the docked compounds.

Hydrogen-bond detection, configuration assignment and rotamer correction of side-chain amides in large proteins by NMR spectroscopy through protium/deuterium isotope effects.

The configuration and hydrogen-bonding network of side-chain amides in a 35 kDa protein were determined by measuring differential and trans-hydrogen-bond H/D isotope effects by using the isotopomer-selective (IS)-TROSY technique, which leads to a reliable recognition and correction of erroneous rotamers that are frequently found in protein structures. First, the differential two-bond isotope effects on carbonyl (13)C' shifts, which are defined as Delta(2)Delta(13)C'(ND) = (2)Delta(13)C'(ND(E))-(2)Delta(13)C'(ND(Z)), provide a reliable means for the configuration assignment for side-chain amides, because environmental effects (hydrogen bonds and charges, etc.) are greatly attenuated over the two bonds that separate the carbon and hydrogen atoms, and the isotope effects fall into a narrow range of positive values. Second and more importantly, the significant variations in the differential one-bond isotope effects on (15)N chemical shifts, which are defined as Delta(1)Delta(15)N(D) = (1)Delta(15)N(D(E))-(1)Delta(15)N(D(Z)) can be correlated with hydrogen-bonding interactions, particularly those involving charged acceptors. The differential one-bond isotope effects are additive, with major contributions from intrinsic differential conjugative interactions between the E and Z configurations, H-bonding interactions, and charge effects. Furthermore, the pattern of trans-H-bond H/D isotope effects can be mapped onto more complicated hydrogen-bonding networks that involve bifurcated hydrogen-bonds. Third, the correlations between Delta(1)Delta(15)N(D) and hydrogen-bonding interactions afford an effective means for the correction of erroneous rotamer assignments of side-chain amides. Rotamer correction by differential isotope effects is not only robust, but also simple and can be applied to large proteins.

Structural basis for the aldolase and epimerase activities of Staphylococcus aureus dihydroneopterin aldolase.

Dihydroneopterin aldolase (DHNA) catalyzes the conversion of 7,8-dihydroneopterin (DHNP) to 6-hydroxymethyl-7,8-dihydropterin (HP) and the epimerization of DHNP to 7,8-dihydromonopterin (DHMP). Although crystal structures of the enzyme from several microorganisms have been reported, no structural information is available about the critical interactions between DHNA and the trihydroxypropyl moiety of the substrate, which undergoes bond cleavage and formation. Here, we present the structures of Staphylococcus aureus DHNA (SaDHNA) in complex with neopterin (NP, an analog of DHNP) and with monapterin (MP, an analog of DHMP), filling the gap in the structural analysis of the enzyme. In combination with previously reported SaDHNA structures in its ligand-free form (PDB entry 1DHN) and in complex with HP (PDB entry 2DHN), four snapshots for the catalytic center assembly along the reaction pathway can be derived, advancing our knowledge about the molecular mechanism of SaDHNA-catalyzed reactions. An additional step appears to be necessary for the epimerization of DHMP to DHNP. Three active site residues (E22, K100, and Y54) function coordinately during catalysis: together, they organize the catalytic center assembly, and individually, each plays a central role at different stages of the catalytic cycle.

Mechanism of dihydroneopterin aldolase. NMR, equilibrium and transient kinetic studies of the Staphylococcus aureus and Escherichia coli enzymes.

Dihydroneopterin aldolase (DHNA) catalyzes both the cleavage of 7,8-dihydro-D-neopterin (DHNP) to form 6-hydroxymethyl-7,8-dihydropterin (HP) and glycolaldehyde and the epimerization of DHNP to form 7,8-dihydro-L-monapterin (DHMP). Whether the epimerization reaction uses the same reaction intermediate as the aldol reaction or the deprotonation and reprotonation of C2' of DHNP has been investigated by NMR analysis of the reaction products in a D2O solvent. No deuteration of C2' was observed for the newly formed DHMP. This result strongly suggests that the epimerization reaction uses the same reaction intermediate as the aldol reaction. In contrast with an earlier observation, the DHNA-catalyzed reaction is reversible, which also supports a nonstereospecific retroaldol/aldol mechanism for the epimerization reaction. The binding and catalytic properties of DHNAs from both Staphylococcus aureus (SaDHNA) and Escherichia coli (EcDHNA) were determined by equilibrium binding and transient kinetic studies. A complete set of kinetic constants for both the aldol and epimerization reactions according to a unified kinetic mechanism was determined for both SaDHNA and EcDHNA. The results show that the two enzymes have significantly different binding and catalytic properties, in accordance with the significant sequence differences between them.

Catalytic mechanism of guanine deaminase: an ONIOM and molecular dynamics study.

The catalytic mechanism of Bacillus subtilis guanine deaminase (bGD), a Zn metalloenzyme, has been investigated by a combination of quantum mechanical calculations using the multilayered ONIOM method and molecular dynamics simulations. In contrast to a previously proposed catalytic mechanism, which requires the bound guanine to assume a rare tautomeric state, the ONIOM calculations showed that the active-site residues of the enzyme do not affect the tautomeric state of guanine, and consequently the bound guanine is a tautomer that is the most abundant in aqueous solution. Two residues, Glutamate 55 and Aspartate 114, were found to play important roles in proton shuttling in the reaction. The proposed reaction path is initiated by proton transfer from a Zn-bound water to protonate Asp114. This process may be quite complex and rather dynamic in nature, as revealed by the molecular dynamics (MD) simulations, whereby another water may bridge the Zn-bound water and Asp114, which then is eliminated by positioning of guanine in the active site. The binding of guanine stabilizes protonated Asp114 by hydrogen bond formation. Asp114 can then transfer its proton to the N3 of the bound guanine, facilitating the nucleophilic attack on C2 of the guanine by the Zn-bound hydroxide to form a tetrahedral intermediate. This occurs with a rather low barrier. Glu55 then transfers a proton from the Zn-hydroxide to the amino group of the reaction intermediate and, at this point, the C2-N2 bond has lengthened by 0.2 A compared to guanine, making C2-N2 bond cleavage more facile. The C2-N2 bond breaks forming ammonia, with an energy barrier of approximately 8.8 kcal/mol. Ammonia leaves the active site, and xanthine is freed by the cleavage of the Zn-O2 bond, with a barrier approximately 8.4 kcal/mol. Along this reaction path, the highest barrier comes from C2-N2 bond cleavage, while the barrier from the cleavage of the Zn-O2 bond is slightly smaller. The Zn-O2 bond can be broken without the assistance of water during the release of xanthine.

TROSY of side-chain amides in large proteins.

By using the mixed solvent of 50% H2O/50% D2O and employing deuterium decoupling, TROSY experiments exclusively detect NMR signals from semideuterated isotopomers of carboxamide groups with high sensitivities for proteins with molecular weights up to 80 kDa. This isotopomer-selective strategy extends TROSY experiments from exclusively detecting backbone to both backbone and side-chain amides, particularly in large proteins. Because of differences in both TROSY effect and dynamics between 15N-H(E){D(Z)} and 15N-H(Z){D(E)} isotopomers of the same carboxamide, the 15N transverse magnetization of the latter relaxes significantly faster than that of the former, which provides a direct and reliable stereospecific distinction between the two configurations. The TROSY effects on the 15N-H(E){D(Z)} isotopomers of side-chain amides are as significant as on backbone amides.