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.

α-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.

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.

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.

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.

(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.

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.

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.

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.

The structure of H5N1 avian influenza neuraminidase suggests new opportunities for drug design.

The worldwide spread of H5N1 avian influenza has raised concerns that this virus might acquire the ability to pass readily among humans and cause a pandemic. Two anti-influenza drugs currently being used to treat infected patients are oseltamivir (Tamiflu) and zanamivir (Relenza), both of which target the neuraminidase enzyme of the virus. Reports of the emergence of drug resistance make the development of new anti-influenza molecules a priority. Neuraminidases from influenza type A viruses form two genetically distinct groups: group-1 contains the N1 neuraminidase of the H5N1 avian virus and group-2 contains the N2 and N9 enzymes used for the structure-based design of current drugs. Here we show by X-ray crystallography that these two groups are structurally distinct. Group-1 neuraminidases contain a cavity adjacent to their active sites that closes on ligand binding. Our analysis suggests that it may be possible to exploit the size and location of the group-1 cavity to develop new anti-influenza drugs.

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.

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.

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.

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)(-).

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.

Anionic charge is prioritized over geometry in aluminum and magnesium fluoride transition state analogs of phosphoryl transfer enzymes.

Phosphoryl transfer reactions are ubiquitous in biology and metal fluoride complexes have played a central role in structural approaches to understanding how they are catalyzed. In particular, numerous structures of AlFx-containing complexes have been reported to be transition state analogs (TSAs). A survey of nucleotide kinases has proposed a correlation between the pH of the crystallization solution and the number of coordinated fluorides in the resulting aluminum fluoride TSA complexes formed. Enzyme ligands crystallized above pH 7.0 were attributed to AlF3, whereas those crystallized at or below pH 7.0 were assigned as AlF4-. We use 19F NMR to show that for beta-phosphoglucomutase from Lactococcus lactis, the pH-switch in fluoride coordination does not derive from an AlF4- moiety converting into AlF3. Instead, AlF4- is progressively replaced by MgF3- as the pH increases. Hence, the enzyme prioritizes anionic charge at the expense of preferred native trigonal geometry over a very broad range of pH. We demonstrate similar behavior for two phosphate transfer enzymes that represent typical biological phosphate transfer catalysts: an amino acid phosphatase, phosphoserine phosphatase from Methanococcus jannaschii and a nucleotide kinase, phosphoglycerate kinase from Geobacillus stearothermophilus. Finally, we establish that at near-physiological ratios of aluminum to magnesium, aluminum can dominate over magnesium in the enzyme-metal fluoride inhibitory TSA complexes, and hence is the more likely origin of some of the physiological effects of fluoride.

Fhit proteins can also recognize substrates other than dinucleoside polyphosphates.

We show here that Fhit proteins, in addition to their function as dinucleoside triphosphate hydrolases, act similarly to adenylylsulfatases and nucleoside phosphoramidases, liberating nucleoside 5'-monophosphates from such natural metabolites as adenosine 5'-phosphosulfate and adenosine 5'-phosphoramidate. Moreover, Fhits recognize synthetic nucleotides, such as adenosine 5'-O-phosphorofluoridate and adenosine 5'-O-(gamma-fluorotriphosphate), and release AMP from them. With respect to the former, Fhits behave like a phosphodiesterase I concomitant with cleavage of the P-F bond. Some kinetic parameters and implications of the novel reactions catalyzed by the human and plant (Arabidopsis thaliana) Fhit proteins are presented.

Carbonic anhydrase inhibitors. Interaction of 2-N,N-dimethylamino-1,3,4-thiadiazole-5-methanesulfonamide with 12 mammalian isoforms: kinetic and X-ray crystallographic studies.

2-N,N-Dimethylamino-1,3,4-thiadiazole-5-methanesulfonamide was tested for its interaction with the 12 catalytically active mammalian carbonic anhydrase (CA, EC 4.2.1.1) isozymes, CA I-XIV. The compound is a potent inhibitor of CA IV, VII, IX, XII, and XIII (K(I)s of 0.61-39 nM), a medium potency inhibitor of CA II and VA (K(I)s of 121-438 nM), and a weak inhibitor against the other isoforms (CA III, VB, VI, and XIV), making it a very interesting candidate for situations in which a strong/selective inhibition of certain isozymes is needed. The crystal structure of the hCA II adduct of this sulfonamide revealed interesting interactions between the inhibitor and the enzyme which are quite different from those observed in the adducts of CA II with the structurally related aliphatic derivatives zonisamide, 2-amino-1,3,4-thiadiazolyl-5-difluoromethanesulfonamide, and 2-dimethylamino-5-[sulfonamido-(aminomethyl)]-1,3,4-thiadiazole reported earlier.

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.

Turnover-based in vitro selection and evolution of biocatalysts from a fully synthetic antibody library.

This report describes the selection of highly efficient antibody catalysts by combining chemical selection from a synthetic library with directed in vitro protein evolution. Evolution started from a naive antibody library displayed on phage made from fully synthetic, antibody-encoding genes (the Human Combinatorial Antibody Library; HuCAL-scFv). HuCAL-scFv was screened by direct selection for catalytic antibodies exhibiting phosphatase turnover. The substrate used was an aryl phosphate, which is spontaneously transformed into an electrophilic trapping reagent after cleavage. Chemical selection identified an efficient biocatalyst that then served as a template for error-prone PCR (epPCR) to generate randomized repertoires that were subjected to further selection cycles. The resulting superior catalysts displayed cumulative mutations throughout the protein sequence; the ten-fold improvement of their catalytic proficiencies (>10(10) M(-1)) resulted from increased kcat values, thus demonstrating direct selection for turnover. The strategy described here makes the search for new catalysts independent of the immune system and the antibody framework.