Structure-based dissection of the active site chemistry of leukotriene A4 hydrolase: implications for M1 aminopeptidases and inhibitor design

M1 aminopeptidases comprise a large family of biologically important zinc enzymes. We show that peptide turnover by the M1 prototype, leukotriene A4 hydrolase/aminopeptidase, involves a shift in substrate position associated with exchange of zinc coordinating groups, while maintaining the overall coordination geometry. The transition state is stabilized by residues conserved among M1 members and in the final reaction step, Glu-296 of the canonical zinc binding HEXXH motif shuffles a proton from the hydrolytic water to the leaving group. Tripeptide substrates bind along the conserved GXMEN motif, precisely occupying the distance between Glu-271 and Arg-563, whereas the Arg specificity is governed by a narrow S1 pocket capped with Asp-375. Our data provide detailed insights to the active site chemistry of M1 aminopeptidases and will aid in the development of novel enzyme inhibitors.     

Theoretical studies on pyridoxal 5'-phosphate-dependent transamination of alpha-amino acids

Density functional methods have been applied to investigate the irreversible transamination between glyoxylic acid and pyridoxamine analog and the catalytic mechanism for the critical [1,3] proton transfer step in aspartate aminotransferase (AATase). The results indicate that the catalytic effect of pyridoxal 5'-phosphate (PLP) may be attributed to its ability to stabilize related transition states through structural resonance. Additionally, the PLP hydroxyl group and the carboxylic group of the amino acid can shuttle proton, thereby lowering the barrier. The rate-limiting step is the tautomeric conversion of the aldimine to ketimine by [1,3] proton transfer, with a barrier of 36.3 kcal/mol in water solvent. A quantum chemical model consisting 142 atoms was constructed based on the crystal structure of the native AATase complex with the product L-glutamate. The electron-withdrawing stabilization by various residues, involving Arg386, Tyr225, Asp222, Asn194, and peptide backbone, enhances the carbon acidity of 4'-C of PLP and Calpha of amino acid. The calculations support the proposed proton transfer mechanism in which Lys258 acts as a base to shuttle a proton from the 4'-C of PLP to Calpha of amino acid. The first step (proton transfer from 4'-C to lysine) is shown to be the rate-limiting step. Furthermore, we provided an explanation for the reversibility and specificity of the transamination in AATase.     

The Catalytic Mechanism of Fluoroacetate Dehalogenase: A Computational Exploration of Biological Dehalogenation

The biological dehalogenation of fluoroacetate carried out by fluoroacetate dehalogenase is discussed by using quantum mechanical/molecular mechanical (QM/MM) calculations for a whole‐enzyme model of 10 800 atoms. Substrate fluoroacetate is anchored by a hydrogen‐bonding network with water molecules and the surrounding amino acid residues of Arg105, Arg108, His149, Trp150, and Tyr212 in the active site in a similar way to haloalkane dehalogenase. Asp104 is likely to act as a nucleophile to attack the α‐carbon of fluoroacetate, resulting in the formation of an ester intermediate, which is subsequently hydrolyzed by the nucleophilic attack of a water molecule to the carbonyl carbon atom. The cleavage of the strong C? F bond is greatly facilitated by the hydrogen‐bonding interactions between the leaving fluorine atom and the three amino acid residues of His149, Trp150, and Tyr212. The hydrolysis of the ester intermediate is initiated by a proton transfer from the water molecule to His271 and by the simultaneous nucleophilic attack of the water molecule. The transition state and produced tetrahedral intermediate are stabilized by Asp128 and the oxyanion hole composed of Phe34 and Arg105.     

Exploring nucleoside hydrolase catalysis in silico: molecular dynamics study of enzyme-bound substrate and transition state

The mechanism of action of inosine-uridine nucleoside hydrolase has been investigated by long-term molecular dynamics (MD) simulation in TIP3P water using stochastic boundary conditions. Five MD studies have been performed with enzyme substrate complex (E.S), enzyme substrate complex with protonated His241 (EH.S), enzyme transition state complex (E.TS), enzyme transition state complex with protonated His241 (EH.TS), and His241Ala transition state complex E(H241A).TS. Special attention has been given to the role of His241, which has been considered as the general acid catalyst to assist departure of the leaving nucleobase on the basis of its location in the active site in the X-ray crystal structure (). Yet on the basis of the location in the active site, Tyr229 is closer to the aniline ring of pAPIR as compared to His241. On initiation of MD simulations, His241 does not approach the nucleobase in the structures of EH.S, E.S, EH.TS, and E.TS. In the solvated enzyme, Tyr229, which is a member of the hydrogen bonding network inosine O2'.Asp14.His241.Tyr229.inosine N7, serves as a proton source to the leaving nucleobase. The loss of significant activity of His241Ala mutant is shown to be related to the disruption of the above hydrogen bonded network and the distancing of Tyr229 from inosine N7. The structures of the enzyme complexes with substrate or TS are not visibly altered on protonation of His241, a most unusual outcome. The bell-shaped pH dependence upon pK(app)'s of 7.1 and 9.1 may be attributed to the necessity of the dissociation of Asp10 or Asp15 and the acid form of Tyr229, respectively. In TS, the residue Ile81 migrated closer, whereas Arg233 moved away from the nucleobase. The probability of ribooxocarbenium ion stabilization by Asn168 and Asp14 is discussed. The Asp14-CO(2)(-) is hydrogen bonded to the ribose 2'-OH for 96% of the MD simulation time. Nucleophilic addition of water138 to ribooxocarbenium ion is suggested to be assisted by the proton shuttle from water138 --> Asp10 --> Asp15 --> water pool. An anticorrelation motion between Tyr229-OH and Asn168-OD1 in EH.S and E.S is observed. The relationship of this anticorrelated motion to mechanism, if any, deserves further exploration, perhaps the formation of a near attack conformation.     

The role of the putative catalytic base in the phosphoryl transfer reaction in a protein kinase: first-principles calculations

Protein kinases are important enzymes controlling the majority of cellular signaling events via a transfer of the gamma-phosphate of ATP to a target protein. Even after many years of study, the mechanism of this reaction is still poorly understood. Among many factors that may be responsible for the 1011-fold rate enhancement due to this enzyme, the role of the conserved aspartate (Asp166) has been given special consideration. While the essential presence of Asp166 has been established by mutational studies, its function is still debated. The general base catalyst role assigned to Asp166 on the basis of its position in the active site has been brought into question by the pH dependence of the reaction rate, isotope measurements, and pre-steady-state kinetics. Recent semiempirical calculations have added to the controversy surrounding the role of Asp166 in the catalytic mechanism. No major role for Asp166 has been found in these calculations, which have predicted the reaction process consisting of an early transfer of a substrate proton onto the phosphate group. These conclusions were inconsistent with experimental observations. To address these differences between experimental results and theory with a more reliable computational approach and to provide a theoretical platform for understanding catalysis in this important enzyme family, we have carried out first-principles structural and dynamical calculations of the reaction process in cAPK kinase. To preserve the essential features of the reaction, representations of all of the key conserved residues (82 atoms) were included in the calculation. The structural calculations were performed using the local basis density functional (DFT) approach with both hybrid B3LYP and PBE96 generalized gradient approximations. This kind of calculation has been shown to yield highly accurate structural information for a large number of systems. The optimized reactant state structure is in good agreement with X-ray data. In contrast to semiempirical methods, the lowest energy product state places the substrate proton on Asp166. First-principles molecular dynamics simulations provide additional support for the stability of this product state. The latter also demonstrate that the proton transfer to Asp166 occurs at a point in the reaction where bond cleavage at the PO bridging position is already advanced. This mechanism is further supported by the calculated structure of the transition state in which the substrate hydroxyl group is largely intact. A metaphoshate-like structure is present in the transition state, which is consistent with the X-ray structures of transition state mimics. On the basis of the calculated structure of the transition state, it is estimated to be 85% dissociative. Our analysis also indicates an increase in the hydrogen bond strength between Asp166 and substrate hydroxyl and a small decrease in the bond strength of the latter in the transition state. In summary, our calculations demonstrate the importance of Asp166 in the enzymatic mechanism as a proton acceptor. However, the proton abstraction from the substrate occurs late in the reaction process. Thus, in the catalytic mechanism of cAPK protein kinase, Asp166 plays a role of a "proton trap" that locks the transferred phosphoryl group to the substrate. These results resolve prior inconsistencies between theory and experiment and bring new understanding of the role of Asp166 in the protein kinase catalytic mechanism.     

Molecular dynamics simulations of the mononuclear zinc-beta-lactamase from Bacillus cereus complexed with benzylpenicillin and a quantum chemical study of the reaction mechanism

Herein, we present results from MD simulations of the Michaelis complex formed between the B. cereus zinc-beta-lactamase enzyme and benzylpenicillin. The structural and dynamical effects induced by substrate-binding, the specific role of the conserved residues, and the near attack conformers of the Michaelis complex are discussed. Quantum chemical methods (HF/6-31G* and B3LYP/6-31G*) are also applied to study the hydrolysis reaction of N-methylazetidinone catalyzed by a monozinc system consisting of the side chains of the histidine residues (His86, His88, and His149) complexed with Zn-OH and the side chains of Asp90 and His210. From this model system, we built molecular-mechanics representations of the prereactive complex and transition state configurations docked into the active site. Linear-scaling semiempirical calculations coupled with a continuum solvent model were then performed on these static models. We propose that the experimental rate data for the B. cereus enzyme is compatible with a one-step mechanism for the hydrolysis of beta-lactam substrates in which His210 acts as a proton donor.     

Multiple Proton Confinement in the M2 Channel from the Influenza A Virus

The tetrameric M2 protein bundle of the influenza A virus is the proton channel responsible for the acidification of the viral interior, a key step in the infection cycle. Selective proton transport is achieved by successive protonation of the conserved histidine amino acids at position 37. A recent X-ray structure of the tetrameric transmembrane (TM) domain of the protein (residues 22-46) resolved several water clusters in the channel lumen, which suggest possible proton pathways to the His37 residues. To explore this hypothesis, we have carried out molecular dynamics (MD) simulations of a proton traveling towards the His37 side chains using MD with classical and quantum force fields. Diffusion through the first half of the channel to the "entry" water cluster near His37 may be hampered by significant kinetic barriers due to electrostatic repulsion. However, once in the entry cluster, a proton can move to one of the acceptor His37 in a nearly barrierless fashion, as evidenced both by MD simulations and a scan of the potential energy surface (PES). Water molecules of the entry cluster, although confined in the M2 pore and restricted in their motions, can conduct protons with a rate very similar to that of bulk water.     

Theoretical investigation on the dissociation of (R)‐benzoin catalyzed by benzaldehyde lyase

Benzaldehyde lyase (BAL) is a versatile thiamin diphosphate (THDP)‐dependent enzyme with widespread synthetic applications in industry. Besides lyase activity, BAL also performs the functions as carboligase and decarboxylase. Unlike many other THDP‐dependent enzymes, the active center of BAL is devoid of any acid‐base amino acid residues except Glu50 and His29, and therefore, the catalytic mechanism of BAL is unusual. In this article, the dissociation mechanism of (R)‐benzoin to benzaldehyde catalyzed by BAL has been studied by using density functional theory method. The calculation results indicate that the whole reaction consists of four elementary steps, and at least two steps contribute to rate‐limiting. A big difference with other THDP‐dependent enzymes is that, in the first stage of the reaction, the ligation of substrate and THDP ylide is not companied by proton transfer, and in the subsequent transition states and intermediates, the carbonyl oxygen always exists in the form of anion. Gln113, His29, and 4′‐amino group of THDP are found to have the function to stabilize the transition states and intermediates. His29 acts as the proton acceptor in step 2 and proton donor in step 3 using one water molecule as mediator. © 2013 Wiley Periodicals, Inc.     

Spectroscopic and computational studies of α-keto acid binding to Dke1: understanding the role of the facial triad and the reactivity of β-diketones

The O(2) activating mononuclear nonheme iron enzymes generally have a common facial triad (two histidine and one carboxylate (Asp or Glu) residue) ligating Fe(II) at the active site. Exceptions to this motif have recently been identified in nonheme enzymes, including a 3His triad in the diketone cleaving dioxygenase Dke1. This enzyme is used to explore the role of the facial triad in directing reactivity. A combination of spectroscopic studies (UV-vis absorption, MCD, and resonance Raman) and DFT calculations is used to define the nature of the binding of the α-keto acid, 4-hydroxyphenlpyruvate (HPP), to the active site in Dke1 and the origin of the atypical cleavage (C2-C3 instead of C1-C2) pattern exhibited by this enzyme in the reaction of α-keto acids with dioxygen. The reduced charge of the 3His triad induces α-keto acid binding as the enolate dianion, rather than the keto monoanion, found for α-keto acid binding to the 2His/1 carboxylate facial triad enzymes. The mechanistic insight from the reactivity of Dke1 with the α-keto acid substrate is then extended to understand the reaction mechanism of this enzyme with its native substrate, acac. This study defines a key role for the 2His/1 carboxylate facial triad in α-keto acid-dependent mononuclear nonheme iron enzymes in stabilizing the bound α-keto acid as a monoanion for its decarboxylation to provide the two additional electrons required for O(2) activation.     

Role of Asp102 in the catalytic relay system of serine proteases: a theoretical study

The role of Asp102 in the catalytic relay system of serine proteases is studied theoretically by calculating the free energy profiles of the single proton-transfer reaction by the Asn102 mutant trypsin and the concerted double proton-transfer reaction (so-called the charge-relay mechanism) of the wild-type trypsin. For each reaction, the reaction free energy profile of the rate-determining step (the tetrahedral intermediate formation step) is calculated by using ab initio QM/MM electronic structure calculations combined with molecular dynamics-free energy perturbation method. In the mutant reaction, the free energy monotonically increases along the reaction path. The rate-determining step of the mutant reaction is the formation of tetrahedral intermediate complex, not the base (His57) abstraction of the proton from Ser195. In contrast to the single proton-transfer reaction of the wild-type, MD simulations of the enzyme-substrate complex show that the catalytically favorable alignment of the relay system (the hydrogen bonding network between the mutant triad, His57, Asn102, and Ser195) is rarely observed even in the presence of a substrate at the active site. In the double proton-transfer reaction, the energy barrier is observed at the proton abstraction step, which corresponds to the rate-determining step of the single proton-transfer reaction of the wild-type. Although both reaction profiles show an increase of the activation barrier by several kcals/mol, these increases have different energetic origins: a large energetic loss of the electrostatic stabilization between His57 and Asn102 in the mutant reaction, while the lack of stabilization by the protein environment in the double proton-transfer reaction. Comparing the present results with the single proton transfer of the wild-type, Asp102 is proven to play two important roles in the catalytic process. One is to stabilize the protonated His57, or ionic intermediate, formed during the acylation, and the other is to fix the configuration around the active site, which is favorable to promote the catalytic process. These two factors are closely related to each other and are indispensable for the efficient catalysis. Also the present calculations suggest the importance of the remote site interaction between His57 and Val213-Ser214 at the catalytic transition state.     

His75-Asp97 cluster in green proteorhodopsin

The proteorhodopsin (PR) family found in bacteria near the ocean's surface consists of hundreds of PR variants color-tuned to their environment. PR contains a highly conserved single histidine at position 75, which is not found in most other retinal proteins. Using (13)C and (15)N MAS NMR, we were able to prove for green PR that His75 forms a pH-dependent H-bond with the primary proton acceptor Asp97, which explains its unusually high pK(a). The functional role of His75 has been studied using site-directed mutagenesis and time-resolved optical spectroscopy: Ultrafast vis-pump/vis-probe experiments on PR(H75N) showed that the primary reaction dynamics is retained, while flash photolysis experiments revealed an accelerated photocycle. Our data show the formation of a pH-dependent His-Asp cluster which might be typical for eubacterial retinal proteins. Despite its stabilizing function, His75 was found to slow the photocycle in wild-type PR. This means that PR was not optimized by evolution for fast proton transfer, which raises questions about its true function in vivo.     

Conjugation and deconjugation of ubiquitin regulating the destiny of proteins

The homeostasis for a number of cellular proteins is regulated by not only phosphorylation and dephosphorylation, but also ubiquitination and deubiquitination. A number of proteins involved in the degradation of polypeptides have been isolated in various eukaryotic organisms from Saccharomyces cerevisiae to human. Recently, several deubiquitinating enzymes, classified into either the Ub C-terminal hydrolase (UCH) or the Ub-specific processing protease (UBP), have been reported. It has been shown that they contain conserved domains including Cys, His, and Asp residues throughout the enzyme. These proteins have been demonstrated that Cys and His domains are critical for deubiquitinating enzymatic activity. Recently, we have shown that the Asp domain localized between Cys and His domains is also essential for cleaving the ubiquitin from protein substrates. Mouse deubiquitinating enzymes including DUB-1, DUB-2, and DUB-2A have been isolated and they showed the expression specificity. Of these, DUB- 1 and DUB-2 are expressed in lymphocytes depending on the presence of cytokines (interleukin-3 in B-lymphocytes and interleukin-2 in T- lymphocytes, respectively), indicating that they are involved in cytokine signaling pathways. Isolation of all putative DUBs will help to identify their substrates and to regulate the homeostasis of cellular proteins, especially in proliferative cells.     

Ammonium transporters achieve charge transfer by fragmenting their substrate

Proteins of the Amt/MEP family facilitate ammonium transport across the membranes of plants, fungi, and bacteria and are essential for growth in nitrogen-poor environments. Some are known to facilitate the diffusion of the neutral NH(3), while others, notably in plants, transport the positively charged NH(4)(+). On the basis of the structural data for AmtB from Escherichia coli , we illustrate the mechanism by which proteins from the Amt family can sustain electrogenic transport. Free energy calculations show that NH(4)(+) is stable in the AmtB pore, reaching a binding site from which it can spontaneously transfer a proton to a pore-lining histidine residue (His168). The substrate diffuses down the pore in the form of NH(3), while the excess proton is cotransported through a highly conserved hydrogen-bonded His168-His318 pair. This constitutes a novel permeation mechanism that confers to the histidine dyad an essential mechanistic role that was so far unknown.     

A role for internal water molecules in proton affinity changes in the Schiff base and Asp85 for one-way proton transfer in bacteriorhodopsin

Light-induced proton pumping in bacteriorhodospin is carried out through five proton transfer steps. We propose that the proton transfer to Asp85 from the Schiff base in the L-to-M transition is accompanied by the relocation of a water cluster on the cytoplasmic side of the Schiff base from a site close to the Schiff base in L to the Phe219-Thr46 region in M. The water cluster present in L, formed at 170 K, is more rigid than that at room temperature. This may be responsible for blocking the conversion of L to M at 170 K. In the photocycle at room temperature, this water cluster returns to the site close to the Schiff base in N, with a rigid structure similar to that of L at 170 K. The increase in the proton affinity of Asp85, which is a prerequisite for the one-way proton transfer in the M-to-N transition, is suggested to be facilitated by a structural change which disrupts interactions between Asp212 and the Schiff base, and between Asp212 and Arg82. We propose that this liberation of Asp212 is accompanied by a rearrangement of the structure of water molecules between Asp85 and Asp212, stabilizing the protonated Asp85 in M.     

On the role of the conserved aspartate in the hydrolysis of the phosphocysteine intermediate of the low molecular weight tyrosine phosphatase

The usual rate-determining step in the catalytic mechanism of the low molecular weight tyrosine phosphatases involves the hydrolysis of a phosphocysteine intermediate. To explain this hydrolysis, general base-catalyzed attack of water by the anion of a conserved aspartic acid has sometimes been invoked. However, experimental measurements of solvent deuterium kinetic isotope effects for this enzyme do not reveal a rate-limiting proton transfer accompanying dephosphorylation. Moreover, base activation of water is difficult to reconcile with the known gas-phase proton affinities and solution phase pK(a)'s of aspartic acid and water. Alternatively, hydrolysis could proceed by a direct nucleophilic attack by a water molecule. To understand the hydrolysis mechanism, we have used high-level density functional methods of quantum chemistry combined with continuum electrostatics models of the protein and the solvent. Our calculations do not support a catalytic activation of water by the aspartate. Instead, they indicate that the water oxygen directly attacks the phosphorus, with the aspartate residue acting as a H-bond acceptor. In the transition state, the water protons are still bound to the oxygen. Beyond the transition state, the barrier to proton transfer to the base is greatly diminished; the aspartate can abstract a proton only after the transition state, a result consistent with experimental solvent isotope effects for this enzyme and with established precedents for phosphomonoester hydrolysis.     

Density-functional study of mechanisms for the cofactor-free decarboxylation performed by uroporphyrinogen III decarboxylase

Uroporphyrinogen III decarboxylase catalyzes the fifth step in heme biosynthesis: the elimination of carboxyl groups from the four acetate side chains of uroporphyrinogen III to yield coproporphyrinogen III. The enzyme acts by successively protonating each of the four pyrrole rings present in the substrate, thereby allowing decarboxylation of their side chains, but the identity of the proton donors has not been established yet. Tyr164 has been suggested as a proton donor, and Asp86 has been proposed to act either as a proton donor or as an intermediate-stabilizing residue. We have performed density-functional calculations to study this reaction mechanism, and found that the rate-limiting step is substrate protonation, rather than decarboxylation. Surprisingly, whereas Tyr164 is unable to protonate the substrate, this protonation can be effected by a nearby arginine residue (Arg37), with a free energy barrier of 21.4 kcal.mol(-1), in remarkable agreement with the experimental value of 19.5 kcal.mol(-1). The central positioning of this residue in close proximity to all four pyrrole rings in the substrate may play a key role in the sequential activation of each of these moieties.     

Fluorescence quantum yield and photochemistry of bacteriophytochrome constructs

Bacteriophytochromes (Bphs) are red-light photoreceptor proteins with a photosensory core that consists of three distinct domains, PAS, GAF and PHY, and covalently binds biliverdin (BV) to a conserved cysteine in the PAS domain. In a recent development, PAS-GAF variants were engineered for use as a near-infrared fluorescent marker in mammalian tissues (Tsien and co-workers, Science, 2009, 324, 804-807). Here, we report the fluorescence quantum yield and photochemistry of two highly-related Bphs from Rps. palustris, RpBphP2 (P2) and RpBphP3 (P3) with distinct photoconversion and fluorescence properties. We applied ultrafast spectroscopy to wild type P3 and P2 PAS-GAF proteins and their P3 D216A, Y272F and P2 D202A PAS-GAF-PHY mutant proteins. In these mutants hydrogen-bond interactions between a conserved aspartate (Asp) which connects the BV chromophore with the PHY domains are disrupted. The excited-state lifetime of the truncated P3 and P2 PAS-GAF proteins was significantly longer than in their PAS-GAF-PHY counterparts that constitute the full photosensory core. Mutation of the conserved Asp to Ala in the PAS-GAF-PHY protein had a similar but larger effect. The fluorescence quantum yields of the P3 D216A and Y272F mutants were 0.066, higher than that of wild type P3 (0.043) and similar to the engineered Bph of Tsien and co-workers. We conclude that elimination of a key hydrogen-bond interaction between Asp and a conserved Arg in the PHY domain is responsible for the excited-state lifetime increase in all Bph variants studied here. H/D exchange resulted in a 1.4-1.7 fold increase of excited-state lifetime. The results support a reaction model in which deactivation of the BV chromophore proceeds via excited-state proton transfer from the BV pyrrole nitrogens to the backbone of the conserved Asp or to a bound water. This work may aid in rational structure- and mechanism-based conversion of constructs based on P3 and other BPhs into efficient near-IR, deep tissue, fluorescent markers.     

A novel genetic selection system for PLP-dependent threonine aldolases

Threonine aldolases are versatile pyridoxal-5′-phosphate (PLP)-dependent enzymes key to glycine, serine and threonine metabolism. Because they catalyze the reversible addition of glycine to an aldehyde to give β-hydroxy-α-amino acids, they are also attractive as biotechnological catalysts for the diastereoselective synthesis of many pharmaceutically useful compounds. To study and evolve such enzymes, we have developed a simple selection system based on the simultaneous inactivation of four genes involved in glycine biosynthesis in Escherichia coli. Glycine prototrophy in the deletion strain is restored by expression of a gene encoding an aldolase that converts β-hydroxy-α-amino acids, provided in the medium, to glycine and the corresponding aldehyde. Combinatorial mutagenesis and selection experiments with a previously uncharacterized l-threonine aldolase from Caulobacter crescentus CB15 (Cc-LTA) illustrate the power of this system. The codons for four active site residues, His91, Asp95, Glu96, and Asp176, were simultaneously randomized and active variants selected. The results show that only His91, which π-stacks against the PLP cofactor and probably serves as the catalytic base in the carbon-carbon bond cleavage step, is absolutely required for aldolase activity. In contrast, Asp176, one of the most conserved residues in this enzyme superfamily, can be replaced conservatively by glutamate, albeit with a >5000-fold decrease in efficiency. Though neither Asp95 nor Glu96 is catalytically essential, they appear to modulate substrate binding and His91 activity, respectively. The broad dynamic range of this novel selection system should make it useful for mechanistic investigations and directed evolution of many natural and artificial aldolases.     

A water-mediated and substrate-assisted catalytic mechanism for Sulfolobus solfataricus DNA polymerase IV

DNA polymerases are enzymes responsible for the synthesis of DNA from nucleotides. Understanding their molecular fundamentals is a prerequisite for elucidating their aberrant activities in diseases such as cancer. Here we have carried out ab initio quantum mechanical/molecular mechanical (QM/MM) studies on the nucleotidyl-transfer reaction catalyzed by the lesion-bypass DNA polymerase IV (Dpo4) from Sulfolobus solfataricus, with template guanine and Watson-Crick paired dCTP as the nascent base pair. The results suggested a novel water-mediated and substrate-assisted (WMSA) mechanism: the initial proton transfer to the alpha-phosphate of the substrate via a bridging crystal water molecule is the rate-limiting step, the nucleotidyl-transfer step is associative with a metastable pentacovalent phosphorane intermediate, and the pyrophosphate leaving is facilitated by a highly coordinated proton relay mechanism through mediation of water which neutralizes the evolving negative charge. The conserved carboxylates, which retain their liganding to the two Mg2+ ions during the reaction process, are found to be essential in stabilizing transition states. This WMSA mechanism takes specific advantage of the unique structural features of this low-fidelity lesion-bypass Y-family polymerase, which has a more spacious and solvent-exposed active site than replicative and repair polymerases.