Arginine Kinase Activates Arginine for Phosphorylation by Pyramidalization and Polarization.

Arginine phosphorylation plays numerous roles throughout biology. Arginine kinase (AK) catalyzes the delivery of an anionic phosphoryl group (PO3-) from ATP to a planar, trigonal nitrogen in a guanidinium cation. Density functional theory (DFT) calculations have yielded a model of the transition state (TS) for the AK-catalyzed reaction. They reveal a network of over 50 hydrogen bonds that delivers unprecedented pyramidalization and out-of-plane polarization of the arginine guanidinium nitrogen (Nη2) and aligns the electron density on Nη2 with the scissile P-O bond, leading to in-line phosphoryl transfer via an associative mechanism. In the reverse reaction, the hydrogen-bonding network enforces the conformational distortion of a bound phosphoarginine substrate to increase the basicity of Nη2. This enables Nη2 protonation, which triggers PO3- migration to generate ATP. This polarization-pyramidalization of nitrogen in the arginine side chain is likely a general phenomenon that is exploited by many classes of enzymes mediating the post-translational modification of arginine.

Octahedral Trifluoromagnesate, an Anomalous Metal Fluoride Species, Stabilizes the Transition State in a Biological Motor.

Isoelectronic metal fluoride transition state analogue (TSA) complexes, MgF3 - and AlF4 -, have proven to be immensely useful in understanding mechanisms of biological motors utilizing phosphoryl transfer. Here we report a previously unobserved octahedral TSA complex, MgF3(H2O)-, in a 1.5 Å resolution Zika virus NS3 helicase crystal structure. 19F NMR provided independent validation and also the direct observation of conformational tightening resulting from ssRNA binding in solution. The TSA stabilizes the two conformations of motif V of the helicase that link ATP hydrolysis with mechanical work. DFT analysis further validated the MgF3(H2O)- species, indicating the significance of this TSA for studies of biological motors.

Multiscale computation delivers organophosphorus reactivity and stereoselectivity to immunoglobulin scavengers.

Quantum mechanics/molecular mechanics (QM/MM) maturation of an immunoglobulin (Ig) powered by supercomputation delivers novel functionality to this catalytic template and facilitates artificial evolution of biocatalysts. We here employ density functional theory-based (DFT-b) tight binding and funnel metadynamics to advance our earlier QM/MM maturation of A17 Ig-paraoxonase (WTIgP) as a reactibody for organophosphorus toxins. It enables regulation of biocatalytic activity for tyrosine nucleophilic attack on phosphorus. The single amino acid substitution l-Leu47Lys results in 340-fold enhanced reactivity for paraoxon. The computed ground-state complex shows substrate-induced ionization of the nucleophilic l-Tyr37, now H-bonded to l-Lys47, resulting from repositioning of l-Lys47. Multiple antibody structural homologs, selected by phenylphosphonate covalent capture, show contrasting enantioselectivities for a P-chiral phenylphosphonate toxin. That is defined by crystallographic analysis of phenylphosphonylated reaction products for antibodies A5 and WTIgP. DFT-b analysis using QM regions based on these structures identifies transition states for the favored and disfavored reactions with surprising results. This stereoselection analysis is extended by funnel metadynamics to a range of WTIgP variants whose predicted stereoselectivity is endorsed by experimental analysis. The algorithms used here offer prospects for tailored design of highly evolved, genetically encoded organophosphorus scavengers and for broader functionalities of members of the Ig superfamily, including cell surface-exposed receptors.

Thermokinetic profile of NDM-1 and its inhibition by small carboxylic acids.

The New Delhi metallo-ß-lactamase (NDM-1) is an important clinical target for antimicrobial research, but there are insufficient clinically useful inhibitors and the details of NDM-1 enzyme catalysis remain unclear. The aim of this work is to provide a thermodynamic profile of NDM-1 catalysed hydrolysis of ß-lactams using an isothermal titration calorimetry (ITC) approach and to apply this new method to the identification of new low-molecular-weight dicarboxylic acid inhibitors. The results reveal that hydrolysis of penicillin G and imipenem by NDM-1 share the same thermodynamic features with a significant intrinsic enthalpy change and the release of one proton into solution, while NDM-1 hydrolysis of cefazolin exhibits a different mechanism with a smaller enthalpy change and the release of two protons. The inhibitory constants of four carboxylic acids are found to be in the micromolar range. The compounds pyridine-2,6-dicarboxylic acid and thiazolidine-2,4-dicarboxylic acid show the best inhibitory potency and are confirmed to inhibit NDM-1 using a clinical strain of Escherichia coli The pyridine compound is further shown to restore the susceptibility of this E. coli strain to imipenem, at an inhibitor concentration of 400 µM, while the thiazoline compound also shows a synergistic effect with imipenem. These results provide valuable information to enrich current understanding on the catalytic mechanism of NDM-1 and to aid the future optimisation of ß-lactamase inhibitors based on these scaffolds to tackle the problem of antibiotic resistance.

Assessing the Influence of Mutation on GTPase Transition States by Using X-ray Crystallography, 19 F†NMR, and DFT Approaches.

We report X-ray crystallographic and 19 F NMR studies of the G-protein RhoA complexed with MgF3- , GDP, and RhoGAP, which has the mutation Arg85'Ala. When combined with DFT calculations, these data permit the identification of changes in transition state (TS) properties. The X-ray data show how Tyr34 maintains solvent exclusion and the core H-bond network in the active site by relocating to replace the missing Arg85' sidechain. The 19 F NMR data show deshielding effects that indicate the main function of Arg85' is electronic polarization of the transferring phosphoryl group, primarily mediated by H-bonding to O3G and thence to PG . DFT calculations identify electron-density redistribution and pinpoint why the TS for guanosine 5'-triphosphate (GTP) hydrolysis is higher in energy when RhoA is complexed with RhoGAPArg85'Ala relative to wild-type (WT) RhoGAP. This study demonstrates that 19 F NMR measurements, in combination with X-ray crystallography and DFT calculations, can reliably dissect the response of small GTPases to site-specific modifications.

Metal Fluorides: Tools for Structural and Computational Analysis of Phosphoryl Transfer Enzymes.

The phosphoryl group, PO3-, is the dynamic structural unit in the biological chemistry of phosphorus. Its transfer from a donor to an acceptor atom, with oxygen much more prevalent than nitrogen, carbon, or sulfur, is at the core of a great majority of enzyme-catalyzed reactions involving phosphate esters, anhydrides, amidates, and phosphorothioates. The serendipitous discovery that the phosphoryl group could be labeled by "nuclear mutation," by substitution of PO3- by MgF3- or AlF4-, has underpinned the application of metal fluoride (MF x ) complexes to mimic transition states for enzymatic phosphoryl transfer reactions, with sufficient stability for experimental analysis. Protein crystallography in the solid state and 19F NMR in solution have enabled direct observation of ternary and quaternary protein complexes embracing MF x transition state models with precision. These studies have underpinned a radically new mechanistic approach to enzyme catalysis for a huge range of phosphoryl transfer processes, as varied as kinases, phosphatases, phosphomutases, and phosphohydrolases. The results, without exception, have endorsed trigonal bipyramidal geometry (tbp) for concerted, "in-line" stereochemistry of phosphoryl transfer. QM computations have established the validity of tbp MF x complexes as reliable models for true transition states, delivering similar bond lengths, coordination to essential metal ions, and virtually identical hydrogen bond networks. The emergence of protein control of reactant orbital overlap between bond-forming species within enzyme transition states is a new challenging theme for wider exploration.

Metal Fluorides as Analogues for Studies on Phosphoryl Transfer Enzymes.

The 1994 structure of a transition-state analogue with AlF4- and GDP complexed to G1α, a small G protein, heralded a new field of research into the structure and mechanism of enzymes that manipulate the transfer of phosphoryl (PO3- ) groups. The number of enzyme structures in the PDB containing metal fluorides (MFx ) as ligands that imitate either a phosphoryl or a phosphate group was 357 at the end of 2016. They fall into three distinct geometrical classes: 1) Tetrahedral complexes based on BeF3- that mimic ground-state phosphates; 2) octahedral complexes, primarily based on AlF4- , which mimic "in-line" anionic transition states for phosphoryl transfer; and 3) trigonal bipyramidal complexes, represented by MgF3- and putative AlF30 moieties, which mimic the geometry of the transition state. The interpretation of these structures provides a deeper mechanistic understanding into the behavior and manipulation of phosphate monoesters in molecular biology. This Review provides a comprehensive overview of these structures, their uses, and their computational development.

Robotic QM/MM-driven maturation of antibody combining sites.

In vitro selection of antibodies from large repertoires of immunoglobulin (Ig) combining sites using combinatorial libraries is a powerful tool, with great potential for generating in vivo scavengers for toxins. However, addition of a maturation function is necessary to enable these selected antibodies to more closely mimic the full mammalian immune response. We approached this goal using quantum mechanics/molecular mechanics (QM/MM) calculations to achieve maturation in silico. We preselected A17, an Ig template, from a naïve library for its ability to disarm a toxic pesticide related to organophosphorus nerve agents. Virtual screening of 167,538 robotically generated mutants identified an optimum single point mutation, which experimentally boosted wild-type Ig scavenger performance by 170-fold. We validated the QM/MM predictions via kinetic analysis and crystal structures of mutant apo-A17 and covalently modified Ig, thereby identifying the displacement of one water molecule by an arginine as delivering this catalysis.

(19)F†NMR and DFT Analysis Reveal Structural and Electronic Transition State Features for RhoA-Catalyzed GTP Hydrolysis.

Molecular details for RhoA/GAP catalysis of the hydrolysis of GTP to GDP are poorly understood. We use (19)F NMR chemical shifts in the MgF3(-) transition state analogue (TSA) complex as a spectroscopic reporter to indicate electron distribution for the γ-PO3(-) oxygens in the corresponding TS, implying that oxygen coordinated to Mg has the greatest electron density. This was validated by QM calculations giving a picture of the electronic properties of the transition state (TS) for nucleophilic attack of water on the γ-PO3(-) group based on the structure of a RhoA/GAP-GDP-MgF3(-) TSA complex. The TS model displays a network of 20 hydrogen bonds, including the GAP Arg85' side chain, but neither phosphate torsional strain nor general base catalysis is evident. The nucleophilic water occupies a reactive location different from that in multiple ground state complexes, arising from reorientation of the Gln-63 carboxamide by Arg85' to preclude direct hydrogen bonding from water to the target γ-PO3(-) group.

A novel expression cassette delivers efficient production of exclusively tetrameric human butyrylcholinesterase with improved pharmacokinetics for protection against organophosphate poisoning.

Butyrylcholinesterase is a stoichiometric bioscavenger against poisoning by organophosphorus pesticides and nerve agents. The low level of expression and extremely rapid clearance of monomeric recombinant human butyrylcholinesterase (rhBChE) from bloodstream (t½≈2 min) limits its pharmaceutical application. Recently (Ilyushin at al., PNAS, 2013) we described a long-acting polysialylated recombinant butyrylcholinesterase (rhBChE-CAO), stable in the bloodstream, that protects mice against 4.2 LD50 of VR. Here we report a set of modifications of the initial rhBChE expression vector to improve stability of the enzyme in the bloodstream and increase its production in CHO cells by introducing in the expression cassette: (i) the sequence of the natural human PRAD-peptide in frame with rhBChE gene via "self-processing" viral F2A peptide under control of an hEF/HTLV promoter, and (ii) previously predicted in silico MAR 1-68 and MAR X-29 sequences. This provides fully tetrameric rhBChE (4rhBChE) at 70 mg/l, that displays improved pharmacokinetics (t½ = 32 ± 1.2 h, MRT = 43 ± 2 h). 3D Fluorescent visualization and distribution of (125)I-labeled enzyme reveals similar low level 4rhBChE and rhBChE-CAO accumulation in muscle, fat, and brain. Administered 4rhBChE was mainly catabolized in the liver and breakdown products were excreted in kidney. Injection of 1.2 LD50 and 1.1 LD50 of paraoxon to BALB/c and knockout BChE-/- mice pre-treated with 4rhBChE (50 mg/kg) resulted in 100% and 78% survival, respectively, without perturbation of long-term behavior. In contrast, 100% mortality of non-pre-treated mice was observed. The high expression level of 4rhBChE in CHO cells permits consideration of this new expression system for manufacturing BChE as a biopharmaceutical.

α-Fluorophosphonates reveal how a phosphomutase conserves transition state conformation over hexose recognition in its two-step reaction.

ß-Phosphoglucomutase (ßPGM) catalyzes isomerization of ß-D-glucose 1-phosphate (ßG1P) into D-glucose 6-phosphate (G6P) via sequential phosphoryl transfer steps using a ß-D-glucose 1,6-bisphosphate (ßG16BP) intermediate. Synthetic fluoromethylenephosphonate and methylenephosphonate analogs of ßG1P deliver novel step 1 transition state analog (TSA) complexes for ßPGM, incorporating trifluoromagnesate and tetrafluoroaluminate surrogates of the phosphoryl group. Within an invariant protein conformation, the ß-D-glucopyranose ring in the ßG1P TSA complexes (step 1) is flipped over and shifted relative to the G6P TSA complexes (step 2). Its equatorial hydroxyl groups are hydrogen-bonded directly to the enzyme rather than indirectly via water molecules as in step 2. The (C)O-P bond orientation for binding the phosphate in the inert phosphate site differs by ∼ 30° between steps 1 and 2. By contrast, the orientations for the axial O-Mg-O alignment for the TSA of the phosphoryl group in the catalytic site differ by only ∼ 5°, and the atoms representing the five phosphorus-bonded oxygens in the two transition states (TSs) are virtually superimposable. The conformation of ßG16BP in step 1 does not fit into the same invariant active site for step 2 by simple positional interchange of the phosphates: the TS alignment is achieved by conformational change of the hexose rather than the protein.

Chemical polysialylation of human recombinant butyrylcholinesterase delivers a long-acting bioscavenger for nerve agents in vivo.

The creation of effective bioscavengers as a pretreatment for exposure to nerve agents is a challenging medical objective. We report a recombinant method using chemical polysialylation to generate bioscavengers stable in the bloodstream. Development of a CHO-based expression system using genes encoding human butyrylcholinesterase and a proline-rich peptide under elongation factor promoter control resulted in self-assembling, active enzyme multimers. Polysialylation gives bioscavengers with enhanced pharmacokinetics which protect mice against 4.2 LD(50) of S-(2-(diethylamino)ethyl) O-isobutyl methanephosphonothioate without perturbation of long-term behavior.

Near attack conformers dominate ß-phosphoglucomutase complexes where geometry and charge distribution reflect those of substrate.

Experimental observations of fluoromagnesate and fluoroaluminate complexes of ß-phosphoglucomutase (ß-PGM) have demonstrated the importance of charge balance in transition-state stabilization for phosphoryl transfer enzymes. Here, direct observations of ground-state analog complexes of ß-PGM involving trifluoroberyllate establish that when the geometry and charge distribution closely match those of the substrate, the distribution of conformers in solution and in the crystal predominantly places the reacting centers in van der Waals proximity. Importantly, two variants are found, both of which satisfy the criteria for near attack conformers. In one variant, the aspartate general base for the reaction is remote from the nucleophile. The nucleophile remains protonated and forms a nonproductive hydrogen bond to the phosphate surrogate. In the other variant, the general base forms a hydrogen bond to the nucleophile that is now correctly orientated for the chemical transfer step. By contrast, in the absence of substrate, the solvent surrounding the phosphate surrogate is arranged to disfavor nucleophilic attack by water. Taken together, the trifluoroberyllate complexes of ß-PGM provide a picture of how the enzyme is able to organize itself for the chemical step in catalysis through the population of intermediates that respond to increasing proximity of the nucleophile. These experimental observations show how the enzyme is capable of stabilizing the reaction pathway toward the transition state and also of minimizing unproductive catalysis of aspartyl phosphate hydrolysis.

Chiral phosphoproline-catalyzed asymmetric Michael addition of ketones to nitroolefins: an experimental and theoretical study.

A novel pyrrolidine-based chiral phosphoproline is an effective bifunctional organocatalyst for the asymmetric Michael addition of ketones to nitroolefins giving high levels of diastereo- and enantio-selectivities (up to > 99 : 1 dr and 96% ee). anti-SR Transition state has the lowest barrier which controls the stereoselectivity, in agreement with experimental results.

Prioritization of charge over geometry in transition state analogues of a dual specificity protein kinase.

The direct observation of a transition state analogue (TSA) complex for tyrosine phosphorylation by a signaling kinase has been achieved using (19)F NMR analysis of MEK6 in complex with tetrafluoroaluminate (AlF(4)(-)), ADP, and p38α MAP kinase (acceptor residue: Tyr182). Solvent-induced isotope shifts and chemical shifts for the AlF(4)(-) moiety indicate that two fluorine atoms are coordinated by the two catalytic magnesium ions of the kinase active site, while the two remaining fluorides are liganded by protein residues only. An equivalent, yet distinct, AlF(4)(-) complex involving the alternative acceptor residue in p38α (Thr180) is only observed when the Tyr182 is mutated to phenylalanine. The formation of octahedral AlF(4)(-) species for both acceptor residues, rather than the trigonal bipyramidal AlF(3)(0) previously identified in the only other metal fluoride complex with a protein kinase, shows the requirement of MEK6 for a TSA that is isoelectronic with the migrating phosphoryl group. This requirement has hitherto only been demonstrated for proteins having a single catalytic magnesium ion.

Transition state analogue structures of human phosphoglycerate kinase establish the importance of charge balance in catalysis.

Transition state analogue (TSA) complexes formed by phosphoglycerate kinase (PGK) have been used to test the hypothesis that balancing of charge within the transition state dominates enzyme-catalyzed phosphoryl transfer. High-resolution structures of trifluoromagnesate (MgF(3)(-)) and tetrafluoroaluminate (AlF(4)(-)) complexes of PGK have been determined using X-ray crystallography and (19)F-based NMR methods, revealing the nature of the catalytically relevant state of this archetypal metabolic kinase. Importantly, the side chain of K219, which coordinates the alpha-phosphate group in previous ground state structures, is sequestered into coordinating the metal fluoride, thereby creating a charge environment complementary to the transferring phosphoryl group. In line with the dominance of charge balance in transition state organization, the substitution K219A induces a corresponding reduction in charge in the bound aluminum fluoride species, which changes to a trifluoroaluminate (AlF(3)(0)) complex. The AlF(3)(0) moiety retains the octahedral geometry observed within AlF(4)(-) TSA complexes, which endorses the proposal that some of the widely reported trigonal AlF(3)(0) complexes of phosphoryl transfer enzymes may have been misassigned and in reality contain MgF(3)(-).

Atomic details of near-transition state conformers for enzyme phosphoryl transfer revealed by MgF-3 rather than by phosphoranes.

Prior evidence supporting the direct observation of phosphorane intermediates in enzymatic phosphoryl transfer reactions was based on the interpretation of electron density corresponding to trigonal species bridging the donor and acceptor atoms. Close examination of the crystalline state of beta-phosphoglucomutase, the archetypal phosphorane intermediate-containing enzyme, reveals that the trigonal species is not PO-3 , but is MgF-3 (trifluoromagnesate). Although MgF-3 complexes are transition state analogues rather than phosphoryl group transfer reaction intermediates, the presence of fluorine nuclei in near-transition state conformations offers new opportunities to explore the nature of the interactions, in particular the independent measures of local electrostatic and hydrogen-bonding distributions using 19F NMR. Measurements on three beta-PGM-MgF-3 -sugar phosphate complexes show a remarkable relationship between NMR chemical shifts, primary isotope shifts, NOEs, cross hydrogen bond F...H-N scalar couplings, and the atomic positions determined from the high-resolution crystal structure of the beta-PGM-MgF--3 -G6P complex. The measurements provide independent validation of the structural and isoelectronic MgF--3 model of near-transition state conformations.

Structural tightening and interdomain communication in the catalytic cycle of phosphoglycerate kinase.

Changes in amide-NH chemical shift and hydrogen exchange rates as phosphoglycerate kinase progresses through its catalytic cycle have been measured to assess whether they correlate with changes in hydrogen bonding within the protein. Four representative states were compared: the free enzyme, a product complex containing 3-phosphoglyceric acid (3PG), a substrate complex containing ADP and a transition-state analogue (TSA) complex containing a 3PG-AlF(4)(-)-ADP moiety. There are an overall increases in amide protection from hydrogen exchange when the protein binds the substrate and product ligands and an additional increase when the TSA complex is formed. This is consistent with stabilisation of the protein structure by ligand binding. However, there is no correlation between the chemical shift changes and the protection factor changes, indicating that the protection factor changes are not associated with an overall shortening of hydrogen bonds in the protected ground state, but rather can be ascribed to the properties of the high-energy, exchange-competent state. Therefore, an overall structural tightening mechanism is not supported by the data. Instead, we observed that some cooperativity is exhibited in the N-domain, such that within this domain the changes induced upon forming the TSA complex are an intensification of those induced by binding 3PG. Furthermore, chemical shift changes induced by 3PG binding extend through the interdomain region to the C-domain beta-sheet, highlighting a network of hydrogen bonds between the domains that suggests interdomain communication. Interdomain communication is also indicated by amide protection in one domain being significantly altered by binding of substrate to the other, even where no associated change in the structure of the substrate-free domain is indicated by chemical shifts. Hence, the communication between domains is also manifested in the accessibility of higher-energy, exchange-competent states. Overall, the data that are consistent with structural tightening relate to defined regions and are close to the 3PG binding site and in the hinge regions of 3-phosphoglycerate kinase.

MgF(3)(-) and alpha-galactose 1-phosphate in the active site of beta-phosphoglucomutase form a transition state analogue of phosphoryl transfer.

(19)F-based NMR analysis and hydrogen/deuterium primary isotope shifts establish the formation of a highly populated solution-state trigonal bipyramidal complex involving beta-phosphoglucomutase (beta-PGM), alpha-galactose 1-phosphate (alphaGal1P), and trifluoromagnesate (MgF(3)(-)), PGM-MgF(3)-alphaGal1P, that is a transition state analogue for phosphoryl transfer. Full backbone resonance assignment of the protein shows that its structure is in the closed conformation required for catalytic activity and is closely related to the corresponding complex with glucose 6-phosphate, which we have recently identified using NMR analysis in solution and X-ray crystallography in the solid state. The previous identification of three structural waters in a PGM-alphaGal1P binary substrate complex had indicated that, in the presence of alphaGal1P, magnesium ions, and fluoride, beta-PGM should indeed form a PGM-MgF(3)-alphaGal1P-TSA complex whereas, in the solid-state, apparently it did not. This cast doubt on the validity of the interpretation of MgF(3)(-) complexes. The present work establishes that, in solution, the expectation that a PGM-MgF(3)-alphaGal1P-TSA complex should readily form is fulfilled. These results thus refute the final evidence used to claim that the trigonal bipyramidal species observed in some solid-state structures of complexes involving beta-PGM are pentaoxyphosphorane intermediates.