Faculty Publications

We report results from a computational investigation of the oxidative deboronation of boroglycine, H2N-CH2-B(OH)(2), using H2O and H2O2 as the reactive oxygen species (ROS) to yield aminomethanol, H2N-CH2-OH; these results complement Our study on the protodeboronation of boroglycine to produce methylamine, H2N-CH3 (Larkin et al. J. Phys. Chem. A 2007, 111, 6489-6500). Second-order Moller-Plesset (MP2) perturbation theory with Dunning-Woon correlation-consistent (cc) basis sets were used for the calculations with comparisons made to results from density functional theory (DFT) at the PBE1PBE/6-311++G(d,p)(cc-pVDZ) levels. The effects of a bulk aqueous environment were also incorporated into the calculations employing PCM and CPCM methodology. Using H2O as the ROS, the reaction H2O + H2N-CH2-B(OH)(2) -> H2N-CH2-OH + H-B(OH)(2) was calculated to be endothermic; the value of Delta H-298(0) was + 12.0 kcal/mol at the MP2(FC)/cc-pVTZ computational level in vacuo and + 13.7 kcal/mol in PCM aqueous media; the corresponding value for the activation barrier, Delta H double dagger, was +94.3 kcal/mol relative to the separated reactants in vacuo and +89.9 kcal/mol in PCM aqueous media. In contrast, the reaction H2O2 + H2N-CH2-B(OH)(2) -> H2N-CH2-OH + B(OH)(3) was calculated to be highly exothermic with an Delta H-298(0) value of -100.9 kcal/mol at the MP2(FC)/cc-pVTZ computational level in vacuo and -99.6 kcal/mol in CPCM aqueous media; the highest-energy transition state for the multistep process associated with this reaction involved the rearrangement of H2N-CH2-B(OH)(OOH) to H2N-CH2-O-B(OH)(2) with a Delta H double dagger value of +23.2 kcal/mol in vacuo relative to the separated reactants. These computational results for boroglycine are in accord with the experimental observations for the deboronation of the FDA approved anticancer drug bortezomib (Velcade, PS-341), where it was found to be the principle deactivation pathway (Labutti et al. Chem. Res. Toxicol. 2006, 19, 539-546).

Methionine adenosyltransferases (MATs) are the family of enzymes that synthesize the main biological methyl donor, S-adenosylmethionine. The high sequence conservation among catalytic subunits from bacteria and eukarya preserves key residues that control activity and oligomerization, which is reflected in the protein structure. However, structural differences among complexes with substrates and products have led to proposals of several reaction mechanisms. In parallel, folding studies begin to explain how the three intertwined domains of the catalytic subunit are produced, and to highlight the importance of certain intermediates in attaining the active final conformation. This review analyzes the available structural data and proposes a consensus interpretation that facilitates an understanding of the pathological problems derived from impairment of MAT function. In addition, new research opportunities directed toward clarification of aspects that remain obscure are also identified.

Inosine 5'-monophosphate dehydrogenase (IMPDH) catalyses the rate-limiting step in guanine nucleotide biosynthesis. IMPDH has an evolutionary conserved CBS subdomain of unknown function. The subdomain can be deleted without impairing the in vitro IMPDH catalytic activity and is the site for mutations associated with human retinitis pigmentosa. A guanine-prototrophic Escherichia coli strain, MP101, was constructed with the subdomain sequence deleted from the chromosomal gene for IMPDH. The ATP content was substantially elevated in MP101 whereas the GTP content was slighty reduced. The activities of IMPDH, adenylosuccinate synthetase and GMP reductase were two to threefold lower in MP101 crude extracts compared with the BW25113 wild-type strain. Guanine induced a threefold reduction in the MP101 ATP pool and a fourfold increase in the GTP pool within 10 min of addition to growing cells; this response does not result from the reduced IMPDH activity or starvation for guanylates. In vivo kinetic analysis using 14-C tracers and 33-P pulse-chasing revealed mutation-associated changes in purine nucleotide fluxes and turnover rates. We conclude that the CBS subdomain of IMPDH may coordinate the activities of the enzymes of purine nucleotide metabolism and is essential for maintaining the normal ATP and GTP pool sizes in E. coli.

Enzymes that regulate their activity by modulating an equilibrium of alternate, nonadditive, functionally distinct oligomeric assemblies (morpheeins) constitute a recently described mode of allostery. The oligomeric equilibrium for porphobilinogen synthase (PBGS) consists of high-activity octamers, low-activity hexamers, and two dimer conformations. A phylogenetically diverse allosteric site specific to hexamers is proposed as an inhibitor binding site. Inhibitor binding is predicted to draw the oligomeric equilibrium toward the low-activity hexamer. In silico docking enriched a selection from a small-molecule library for compounds predicted to bind to this allosteric site. In vitro testing of selected compounds identified one compound whose inhibition mechanism is species-specific conversion of PBGS octamers to hexamers. We propose that this strategy for inhibitor discovery can be applied to other proteins that use the morpheein model for allosteric regulation.

In this article the geometrical structure of the simple, achiral, alpha-amino boronic acid boroglycine, H2N-CH2-B(OH)2, was investigated using density functional theory (DFT), second-order Moller-Plesset (MP2) perturbation theory, and coupled cluster methodology with single- and double-excitations (CCSD); the effects of an aqueous environment were incorporated into the results by using a few explicit water molecules and/or self-consistent reaction field (SCRF) calculations with the IEF polarizable continuum model (PCM). Neutral reaction mechanisms were investigated for the direct protodeboronation (hydrolysis) of boroglycine (H2O+H2N-CH2-B(OH)2-->B(OH)3+H2N-CH3), for which DeltaH degrees 298 was -21.9 kcal/mol at the MP2(FC)/aug-cc-pVDZ level, and for the 1,2-carbon-to-nitrogen shift of the -B(OH)2 moiety (H2N-CH2-B(OH)2-->H3C-NH-B(OH)2), for which the corresponding value of DeltaH degrees 298 was -18.2 kcal/mol. A boron-oxygen double-bonded intermediate was found to play an important role in the 1,2-rearrangement mechanism.

Despite the widespread use of boronic acids in materials science and as pharmaceutical agents, many aspects of their structure and reactivity are not well understood. In this research the boronic acid dimer, [HB(OH)(2)](2), was studied by second-order Moller-Plesset (MP2) perturbation theory and coupled cluster methodology with single and double excitations (CCSD). Pople split-valence 6-31+G, 6-311G, and 6-311++G and Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets were employed for the calculations. A doubly hydrogen-bonded conformer (1) of the dimer was consistently found to be lowest in energy; the structure of 1 was planar (C(2)(h)()) at most computational levels employed but was significantly nonplanar (C(2)) at the MP2/6-311++G and CCSD/6-311++G levels, the result of an intrinsic problem with Pople-type sp-diffuse basis functions on heavy atoms. The dimerization energy, enthalpy, and free energy for the formation of (1) from the exo-endo conformer of the monomer were -10.8, -9.2, and +1.2 kcal/mol, respectively, at the MP2/aug-cc-pVTZ level. Several other hydrogen-bonded conformers of the dimer were local minima on the potential energy surface (PES) and ranged from 2 to 5 kcal/mol higher in energy than 1. Nine doubly OH-bridged conformers, in which the boron atoms were tetracoordinated, were also local minima on the PES, but they were all greater than 13 kcal/mol higher in energy than 1; doubly H-bridged structures proved to be transition states. MP2 and CCSD results were compared to those from the BLYP, B3LYP, OLYP, O3LYP, PBE1PBE, and TPSS functionals with the 6-311++G and aug-cc-pVTZ basis sets; the PBE1PBE functional performed best relative to the MP2 and CCSD results. Self-consistent reaction field (SCRF) calculations predict that boronic acid dimerization is less favorable in solution than in vacuo.

S-Adenosylmethionine synthetase (ATP: L-methionine S-adenosyltransferase) catalyzes a two-step reaction in which tripolyphosphate (PPPi) is a tightly bound intermediate. Diimidotriphosphate (O3P-NH-PO2-NH-PO3; PNPNP), a non-hydrolyzable analog of PPPi, is the most potent known inhibitor of AdoMet synthetase with a Ki of 2 nM. The structural basis for the slow, tight-binding inhibition by PNPNP has been investigated by spectroscopic methods. UV difference spectra reveal environmental alterations of arom. protein residues upon PNPNP binding to form the enzyme.2Mg2+.PNPNP complex, and more extensive changes upon formation of the enzyme.2Mg2+.PNPNP.AdoMet complex. Stopped-flow kinetic studies of complex formation revealed that two slow isomerizations follow PNPNP binding in the presence of AdoMet, in contrast to the lower affinity, rapid-equil. binding in the absence of AdoMet. 31P NMR spectra of enzyme complexes with PNPNP revealed electronic perturbations of each phosphorus atom by distinct upfield chem. shifts for each of the three phosphoryl groups in the enzyme.2Mg2+.PNPNP complex, and further upfield shifts of at least 2 resonances in the complex with AdoMet. Comparison of the chem. shifts for the enzyme-bound PNPNP with the enzyme complexes contg. either the product analog O3P-NH-PO3 or O3P-O-PO2-NH-PO3 indicates that the shifts on binding are largest at the binding sites corresponding to those for the a and g phosphoryl groups of the nucleotide (-3.1 to -4.1 ppm), while the resonance at the b phosphoryl group position shifts by -2.1 ppm. EPR spectra of Mn2+ complexes demonstrate spin coupling between the two Mn2+ in both enzyme.2Mn2+.PNPNP and enzyme.2Mn2+.PNPNP.AdoMet, indicating that the metal ions have comparable distances in both cases. The combined results indicate that formation of the highest affinity complex is assocd. with protein side chain rearrangements and increased electron d. at the ligand phosphorus atoms, likely due to ionization of an -NH- group of the inhibitor. The energetic feasibility of ionization of a -NH- group when two Mg2+ ions are bound to O3P-NH-PO3 is supported by d. functional theor. calcns. on model chelates. This mode of interaction is uniquely available to compds. with P-NH-P linkages and may be possible with other proteins in which multiple cations coordinate a polyphosphate chain. [on SciFinder (R)]