Quantum mechanics/molecular mechanics investigation of the chemical reaction in Dpo4 reveals water-dependent pathways and requirements for active site reorganization

The nucleotidyl-transfer reaction coupled with the conformational transitions in DNA polymerases is critical for maintaining the fidelity and efficiency of DNA synthesis. We examine here the possible reaction pathways of a Y-family DNA polymerase, Sulfolobus solfataricus DNA polymerase IV (Dpo4), for the correct insertion of dCTP opposite 8-oxoguanine using the quantum mechanics/molecular mechanics (QM/MM) approach, both from a chemistry-competent state and a crystal closed state. The latter examination is important for understanding pre-chemistry barriers to interpret the entire enzyme mechanism, since the crystal closed state is not an ideal state for initiating the chemical reaction. The most favorable reaction path involves initial deprotonation of O3'H via two bridging water molecules to O1A, overcoming an overall potential energy barrier of approximately 20.0 kcal/mol. The proton on O1A-P(alpha) then migrates to the gamma-phosphate oxygen of the incoming nucleotide as O3' attacks P(alpha), and the P(alpha)-O3A bond breaks. The other possible pathway in which the O3'H proton is transferred directly to O1A on P(alpha) has an overall energy barrier of 25.0 kcal/mol. In both reaction paths, the rate-limiting step is the initial deprotonation, and the trigonal-bipyramidal configuration for P(alpha) occurs during the concerted bond formation (O3'-P(alpha)) and breaking (P(alpha)-O3A), indicating the associative nature of the chemical reaction. In contrast, the Dpo4/DNA complex with an imperfect active-site geometry corresponding to the crystal state must overcome a much higher activation energy barrier (29.0 kcal/mol) to achieve a tightly organized site due to hindered O3'H deprotonation stemming from larger distances and distorted conformation of the proton acceptors. This significant difference demonstrates that the pre-chemistry reorganization in Dpo4 costs approximately 4.0 to 9.0 kcal/mol depending on the primer terminus environment. Compared to the higher fidelity DNA polymerase beta from the X-family, Dpo4 has a higher chemical reaction barrier (20.0 vs 15.0 kcal/mol) due to the more solvent-exposed active site.     

Time-resolved events on the reaction pathway of transcript initiation by a single-subunit RNA polymerase: Raman crystallographic evidence

The nucleotidyl transfer reaction leading to formation of the first phosphodiester bond has been followed in real time by Raman microscopy, as it proceeds in single crystals of the N4 phage virion RNA polymerase (RNAP). The reaction is initiated by soaking nucleoside triphosphate (NTP) substrates and divalent cations into the RNAP and promoter DNA complex crystal, where the phosphodiester bond formation is completed in about 40 min. This slow reaction allowed us to monitor the changes of the RNAP and DNA conformations as well as bindings of substrate and metal through Raman spectra taken every 5 min. Recently published snapshot X-ray crystal structures along the same reaction pathway assisted the spectroscopic assignments of changes in the enzyme and DNA, while isotopically labeled NTP substrates allowed differentiation of the Raman spectra of bases in substrates and DNA. We observed that substrates are bound at 2-7 min after soaking is commenced, the O-helix completes its conformational change, and binding of both divalent metals required for catalysis in the active site changes the conformation of the ribose triphosphate at position +1. These are followed by a slower decrease of NTP triphosphate groups due to phosphodiester bond formation that reaches completion at about 15 min and even slower complete release of the divalent metals at about 40 min. We have also shown that the O-helix movement can be driven by substrate binding only. The kinetics of the in crystallo nucleotidyl transfer reaction revealed in this study suggest that soaking the substrate and metal into the RNAP-DNA complex crystal for a few minutes generates novel and uncharacterized intermediates for future X-ray and spectroscopic analysis.     

Single-molecule Förster resonance energy transfer reveals an innate fidelity checkpoint in DNA polymerase I

Enzymatic reactions typically involve complex dynamics during substrate binding, conformational rearrangement, chemistry, and product release. The noncovalent steps provide kinetic checkpoints that contribute to the overall specificity of enzymatic reactions. DNA polymerases perform DNA replication with outstanding fidelity by actively rejecting noncognate nucleotide substrates early in the reaction pathway. Substrates are delivered to the active site by a flexible fingers subdomain of the enzyme, as it converts from an open to a closed conformation. The conformational dynamics of the fingers subdomain might also play a role in nucleotide selection, although the precise role is currently unknown. Using single-molecule F?rster resonance energy transfer, we observed individual Escherichia coli DNA polymerase I (Klenow fragment) molecules performing substrate selection. We discovered that the fingers subdomain actually samples through three distinct conformations--open, closed, and a previously unrecognized intermediate conformation. We measured the overall dissociation rate of the polymerase-DNA complex and the distribution among the various conformational states in the absence and presence of nucleotide substrates, which were either correct or incorrect. Correct substrates promote rapid progression of the polymerase to the catalytically competent closed conformation, whereas incorrect nucleotides block the enzyme in the intermediate conformation and induce rapid dissociation from DNA. Remarkably, incorrect nucleotide substrates also promote partitioning of DNA to the spatially separated 3'-5' exonuclease domain, providing an additional mechanism to prevent misincorporation at the polymerase active site. These results reveal the existence of an early innate fidelity checkpoint, rejecting incorrect nucleotide substrates before the enzyme encloses the nascent base pair.     

Opposed steric constraints in human DNA polymerase beta and E. coli DNA polymerase I

DNA polymerase selectivity is crucial for the survival of any living species, yet varies significantly among different DNA polymerases. Errors within DNA polymerase-catalyzed DNA synthesis result from the insertion of noncanonical nucleotides and extension of misaligned DNA substrates. The substrate binding characteristics among DNA polymerases are believed to vary in properties such as shape and tightness of the binding pocket, which might account for the observed differences in fidelity. Here, we employed 4'-alkylated nucleotides and primer strands bearing 4'-alkylated nucleotides at the 3'-terminal position as steric probes to investigate differential active site properties of human DNA polymerase beta (Pol beta) and the 3'-->5'-exonuclease-deficient Klenow fragment of E. coli DNA polymerase I (KF(exo-)). Transient kinetic measurements indicate that both enzymes vary significantly in active site tightness at both positions. While small 4'-methyl and -ethyl modifications of the nucleoside triphosphate perturb Pol beta catalysis, extension of modified primer strands is only marginally affected. Just the opposite was observed for KF(exo-). Here, incorporation of the modified nucleotides is only slightly reduced, whereas size augmentation of the 3'-terminal nucleotide in the primer reduces the catalytic efficiency by more than 7000- and 260,000-fold, respectively. NMR studies support the notion that the observed effects derive from enzyme substrate interactions rather than inherent properties of the modified substrates. These findings are consistent with the observed differential capability of the investigated DNA polymerases in fidelity such as processing misaligned DNA substrates. The results presented provide direct evidence for the involvement of varied steric effects among different DNA polymerases on their fidelity.     

Prechemistry barriers and checkpoints do not contribute to fidelity and catalysis as long as they are not rate limiting

In the preceding article, “Perspective: Pre-chemistry conformational changes in DNA polymerase mechanisms” contributed by Schlick and coworkers as well as previous studies of these workers (Schlick et al. in Theor Chem Acc 131:1287, 2012; Radhakrishnan and Schlick in J Am Chem Soc 127:13245–13252, 2005; Radhakrishnan and Schlick in Biochem Biophys Res Commun 350:521–529, 2006; Radhakrishnan et al. in Biochemistry 45:15142–15156, 2006; Radhakrishnan and Schlick in Proc Natl Acad Sci USA 101:5970–5975, 2004) have argued that the conformational changes preceding the chemical step contribute to DNA synthesis and to the fidelity of DNA polymerases. In one of our previous investigations (Ram Prasad and Warshel in Proteins 79:2900–2919, 2011), we argued and showed that as long as the free energy barriers associated with any of the prechemistry steps are not rate limiting, they could not contribute to the catalysis and then to the fidelity. Though all our arguments are based on exact and well-defined scientific logics, Schlick and coworkers seem to overlook some of the clear conditions in these arguments and in particular the requirement that the chemical step is rate limiting in their arguments that the prechemistry barriers contribute to the catalysis. In fact, as long as the prechemistry steps are not rate limiting, we have shown that the enzymes cannot carry the memory of the previous steps. We also address other potential misunderstandings about several key issues; First, we clarify that it is misleading to relate the prechemistry proposal to the clear fact that the substrate-induced conformational changes determine the final preorganization (the issue is the height of the barrier of the enzyme substrate system and not the trivial fact that the enzyme has to change its structure when the substrate binds). Second, we address the presumed role of dynamical effects in enzyme catalysis and the assumption that any observable should be explored in studies of biological function even if they are not relevant to the given effect. Third, we clarify that the fidelity cannot be explained or quantified by invoking the induced fit or conformational selection effects but by evaluating the free energy contributions to the rate-limiting steps from the structures of the corresponding systems (that of course can reflect the induce fit structural changes). Overall, we put a major emphasis on clarifying what is the prechemistry proposal and thus on trying to force the reader to focus on the only real controversy. We of course dismiss any implication that our studies cannot explore mutational effects as we actually pioneered such computational studies and we clarify that in studies of chemical rates, the focus must be placed on evaluating the chemical barriers, rather than on irrelevant factors, but that the calculations of the chemical barriers must consider all the factors that determine this barrier (including metal ions) and also examine if needed different problematic proposals such as dynamical effects, tunneling, and prechemistry.     

Conformational transition pathway of polymerase beta/DNA upon binding correct incoming substrate

The closing conformational transition of wild-type polymerase beta bound to DNA template/primer before the chemical step (nucleotidyl transfer reaction) is simulated using the stochastic difference equation (in length version, "SDEL") algorithm that approximates long-time dynamics. The order of the events and the intermediate states during pol beta's closing pathway are identified and compared to a separate study of pol beta using transition path sampling (TPS) (Radhakrishnan, R.; Schlick, T. Proc. Natl. Acad. Sci. USA 2004, 101, 5970-5975). Results highlight the cooperative and subtle conformational changes in the pol beta active site upon binding the correct substrate that may help explain DNA replication and repair fidelity. These changes involve key residues that differentiate the open from the closed conformation (Asp192, Arg258, Phe272), as well as residues contacting the DNA template/primer strand near the active site (Tyr271, Arg283, Thr292, Tyr296) and residues contacting the beta and gamma phosphates of the incoming nucleotide (Ser180, Arg183, Gly189). This study compliments experimental observations by providing detailed atomistic views of the intermediates along the polymerase closing pathway and by suggesting additional key residues that regulate events prior to or during the chemical reaction. We also show general agreement between two sampling methods (the stochastic difference equation and transition path sampling) and identify methodological challenges involved in the former method relevant to large-scale biomolecular applications. Specifically, SDEL is very quick relative to TPS for obtaining an approximate path of medium resolution and providing qualitative information on the sequence of events; however, associated free energies are likely very costly to obtain because this will require both successful further refinement of the path segments close to the bottlenecks and large computational time.     

Insights from the sea: Structural biology of marine polyketide synthases

Covering: up to the end of 2011The world's oceans are a rich source of natural products with extremely interesting chemistry. Biosynthetic pathways have been worked out for a few, and the story is being enriched with crystal structures of interesting pathway enzymes. By far, the greatest number of structural insights from marine biosynthetic pathways has originated with studies of curacin A, a poster child for interesting marine chemistry with its cyclopropane and thiazoline rings, internal cis double bond, and terminal alkene. Using the curacin A pathway as a model, structural details are now available for a novel loading enzyme with remarkable dual decarboxylase and acetyltransferase activities, an Fe(2+)/α-ketoglutarate-dependent halogenase that dictates substrate binding order through conformational changes, a decarboxylase that establishes regiochemistry for cyclopropane formation, and a thioesterase with specificity for β-sulfated substrates that lead to terminal alkene offloading. The four curacin A pathway dehydratases reveal an intrinsic flexibility that may accommodate bulky or stiff polyketide intermediates. In the salinosporamide A pathway, active site volume determines the halide specificity of a halogenase that catalyzes for the synthesis of a halogenated building block. Structures of a number of putative polyketide cyclases may help in understanding reaction mechanisms and substrate specificities although their substrates are presently unknown.     

DNA polymerase template interactions probed by degenerate isosteric nucleobase analogs

The development of novel artificial nucleobases and detailed X-ray crystal structures for primer/template/DNA polymerase complexes provide opportunities to assess DNA-protein interactions that dictate specificity. Recent results have shown that base pair shape recognition in the context of DNA polymerase must be considered a significant component. The isosteric azole carboxamide nucleobases (compounds 1-5; ) differ only in the number and placement of nitrogen atoms within a common shape and therefore present unique electronic distributions that are shown to dictate the selectivity of template-directed nucleotide incorporation by DNA polymerases. The results demonstrate how nucleoside triphosphate substrate selection by DNA polymerase is a complex phenomenon involving electrostatic interactions in addition to hydrogen bonding and shape recognition. These azole nucleobase analogs offer unique molecular tools for probing nonbonded interactions dictating substrate selection and fidelity of DNA polymerases.     

A quantum mechanical investigation of possible mechanisms for the nucleotidyl transfer reaction catalyzed by DNA polymerase beta

Several quantum mechanical (QM) and hybrid quantum/molecular mechanical (QM/MM) studies have been employed recently to analyze the nucleotidyl transfer reaction in DNA polymerase beta (pol beta). Our examination reveals strong dependence of the reported mechanism on the initial molecular model. Thus, we explore here several model systems by QM methods to investigate pol beta's possible pathway variations. Although our most favorable pathway involves a direct proton transfer from O3'(primer) to O2alpha(Palpha), we also discuss other initial proton-transfer steps--to an adjacent water, to triphosphate, or to aspartic units--and the stabilizing effect of crystallographic water molecules in the active site. Our favored reaction route has an energetically undemanding initial step of less than 1.0 kcal/mol (at the B3LYP/6-31G(d,p) level), and involves a slight rearrangement in the geometry of the active site. This is followed by two major steps: (1) direct proton transfer from O3'(primer) to O2alpha(Palpha) leading to the formation of a pentavalent, trigonal bipyramidal Palpha center, via an associative mechanism, at a cost of about 28 kcal/mol, and (2) breakage of the triphosphate unit (exothermic process, approximately 22 kcal/mol) that results in the full transfer of the nucleotide to the DNA and the formation of pyrophosphate. These energy values are expected to be lower in the physical system when full protein effects are incorporated. We also discuss variations from this dominant pathway, and their impact on the overall repair process. Our calculated barrier for the chemical reaction clearly indicates that chemistry is rate-limiting overall for correct nucleotide insertion in pol beta, in accord with other studies. Protonation studies on relevant intermediates suggest that, although protonation at a single aspartic residue may occur, the addition of a second proton to the system significantly disturbs the active site. We conclude that the active site rearrangement step necessary to attain a reaction-competent geometry is essential and closely related to the "pre-chemistry" avenue described recently as a key step in the overall kinetic cycle of DNA polymerases. Thus, our work emphasizes the many possible ways for DNA polymerase beta's chemical reaction to occur, determined by the active site environment and initial models.     

EPR Spectroscopic Studies of Lipoxygenases

Polyunsaturated fatty acids are sources of diverse natural, and chemically designed products. The enzyme lipoxygenase selectively oxidizes fatty acid acyl chains using controlled free radical chemistry; the products are regio‐ and stereo‐chemically unique hydroperoxides. A conserved structural fold of ≈600 amino acids harbors a long and narrow substrate channel and a well‐shielded catalytic iron. Oxygen, a co‐substrate, is blocked from the active site until a hydrogen atom is abstracted from substrate bis‐allylic carbon, in a non‐heme iron redox cycle. EPR spectroscopy of ferric intermediates in lipoxygenase catalysis reveals changes in the metal coordination and leads to a proposal on the nature of the reactive intermediate. Remarkably, free radicals are so well controlled in lipoxygenase chemistry that spin label technology can be applied as well. The current level of understanding of steps in lipoxygenase catalysis, from the EPR perspective, will be reviewed.     

Fidelity discrimination in DNA polymerase beta: differing closing profiles for a mismatched (G:A) versus matched (G:C) base pair

Understanding fidelity-the faithful replication or repair of DNA by polymerases-requires tracking of the structural and energetic changes involved, including the elusive transient intermediates, for nucleotide incorporation at the template/primer DNA junction. We report, using path sampling simulations and a reaction network model, strikingly different transition states in DNA polymerase beta's conformational closing for correct dCTP versus incorrect dATP incoming nucleotide opposite a template G. The cascade of transition states leads to differing active-site assembly processes toward the "two-metal-ion catalysis" geometry. We demonstrate that these context-specific pathways imply different selection processes: while active-site assembly occurs more rapidly with the correct nucleotide and leads to primer extension, the enzyme remains open longer, has a more transient closed state, and forms product more slowly when an incorrect nucleotide is present. Our results also suggest that the rate-limiting step in pol beta's conformational closing is not identical to that for overall nucleotide insertion and that the rate-limiting step in the overall nucleotide incorporation process for matched as well as mismatched systems occurs after the closing conformational change.     

Hydrophobic amino acid and single-atom substitutions increase DNA polymerase selectivity

DNA polymerase fidelity is of immense biological importance due to the fundamental requirement for accurate DNA synthesis in both replicative and repair processes. Subtle hydrogen-bonding networks between DNA polymerases and their primer/template substrates are believed to have impact on DNA polymerase selectivity. We show that deleting defined interactions of that kind by rationally designed hydrophobic substitution mutations can result in a more selective enzyme. Furthermore, a single-atom replacement within the DNA substrate through chemical modification, which leads to an altered acceptor potential and steric demand of the DNA substrate, further increased the selectivity of the developed systems. Accordingly, this study about the impact of hydrophobic alterations on DNA polymerase selectivity--enzyme and substrate wise--further highlights the relevance of shape complementary and polar interactions on DNA polymerase selectivity.     

Deconstructing Covalent Organocatalysis

Modern organocatalysis has rapidly evolved into an essential component of contemporary organic synthesis. One of the most distinctive aspects of organocatalytic processes is the biomimetic nature in which the catalyst engages the substrate, often forming covalently bound intermediates in a manner reminiscent of enzyme catalysis. Indeed, the process of intramolecularization is often accompanied by a conformational change of the catalyst scaffold, further accentuating this analogy with biological systems. The isolation and study of these catalytic intermediates facilitate the rapid generation of conformation and reactivity profiles to assist in organocatalytic reaction development and/or clarify reaction outcomes. Emulating the formative advances that have derived from studying reaction intermediates in mechanistic organometallic and enzymatic catalysis, the deconstruction of covalently bound organocatalysis intermediates is gaining momentum as a design strategy.     

Probing the substrate specificity of Golgi alpha-mannosidase II by use of synthetic oligosaccharides and a catalytic nucleophile mutant

Inhibition of Golgi alpha-mannosidase II (GMII), which acts late in the N-glycan processing pathway, provides a route to blocking cancer-induced changes in cell surface oligosaccharide structures. To probe the substrate requirements of GMII, oligosaccharides were synthesized that contained an alpha(1,3)- or alpha(1,6)-linked 1-thiomannoside. Surprisingly, these oligosaccharides were not observed in X-ray crystal structures of native Drosophila GMII (dGMII). However, a mutant enzyme in which the catalytic nucleophilic aspartate was changed to alanine (D204A) allowed visualization of soaked oligosaccharides and led to the identification of the binding site for the alpha(1,3)-linked mannoside of the natural substrate. These studies also indicate that the conformational change of the bound mannoside to a high-energy B 2,5 conformation is facilitated by steric hindrance from, and the formation of strong hydrogen bonds to, Asp204. The observation that 1-thio-linked mannosides are not well tolerated by the catalytic site of dGMII led to the synthesis of a pentasaccharide containing the alpha(1,6)-linked Man of the natural substrate and the beta(1,2)-linked GlcNAc moiety proposed to be accommodated by the extended binding site of the enzyme. A cocrystal structure of this compound with the D204A enzyme revealed the molecular interactions with the beta(1,2)-linked GlcNAc. The structure is consistent with the approximately 80-fold preference of dGMII for the cleavage of substrates containing a nonreducing beta(1,2)-linked GlcNAc. By contrast, the lysosomal mannosidase lacks an equivalent GlcNAc binding site and kinetic analysis indicates oligomannoside substrates without non-reducing-terminal GlcNAc modifications are preferred, suggesting that selective inhibitors for GMII could exploit the additional binding specificity of the GlcNAc binding site.     

Protein structural dynamics of photoactive yellow protein in solution revealed by pump-probe X-ray solution scattering

Photoreceptor proteins play crucial roles in receiving light stimuli that give rise to the responses required for biological function. However, structural characterization of conformational transition of the photoreceptors has been elusive in their native aqueous environment, even for a prototype photoreceptor, photoactive yellow protein (PYP). We employ pump-probe X-ray solution scattering to probe the structural changes that occur during the photocycle of PYP in a wide time range from 3.16 μs to 300 ms. By the analysis of both kinetics and structures of the intermediates, the structural progression of the protein in the solution phase is vividly visualized. We identify four structurally distinct intermediates and their associated five time constants and reconstructed the molecular shapes of the four intermediates from time-independent, species-associated difference scattering curves. The reconstructed structures of the intermediates show the large conformational changes such as the protrusion of N-terminus, which is restricted in the crystalline phase due to the crystal contact and thus could not be clearly observed by X-ray crystallography. The protrusion of the N-terminus and the protein volume gradually increase with the progress of the photocycle and becomes maximal in the final intermediate, which is proposed to be the signaling state. The data not only reveal that a common kinetic mechanism is applicable to both the crystalline and the solution phases, but also provide direct evidence for how the sample environment influences structural dynamics and the reaction rates of the PYP photocycle.     

Probing the interaction between prostacyclin synthase and prostaglandin H2 analogues or inhibitors via a combination of resonance Raman spectroscopy and molecular dynamics simulation approaches

In an aim to probe the structure-function relationship of prostacyclin synthase (PGIS), resonance Raman (RR) spectroscopy and molecular dynamic (MD) simulation approaches have been exploited to characterize the heme conformation and heme-protein matrix interactions for human PGIS (hPGIS) and zebrafish PGIS (zPGIS) in the presence and absence of ligands. The high-frequency RR (1300-1700 cm(-1)) indicates that the heme group is in the ferric, six-coordinate, low-spin state for both resting and ligand-bound hPGIS/zPGIS. The low-frequency RR (300-500 cm(-1)) and MD simulation reveal a salient difference in propionate-protein matrix interactions between hPGIS and zPGIS, as evident by a predominant propionate bending vibration at 386 cm(-1) in resting hPGIS, but two vibrations near 370 and 387 cm(-1) in resting zPGIS. Upon binding of a substrate analogue (U46619, U51605, or U44069), both hPGIS and zPGIS induce a distinctive perturbation of the propionate-protein matrix interactions, resulting in similar Raman shifts to ~381 cm(-1). On the contrary, the bending vibration remains unchanged upon binding of inhibitor/ligand (minoxidil, clotrimazole, or miconazole), indicating that these inhibitors/ligands do not interfere with the propionate-protein matrix interactions. These results, together with subtle changes in vinyl bending modes, demonstrate drastically different RR shifts with heme conformational changes in both hPGIS and zPGIS upon different ligand bindings, suggesting that PGIS exhibits a ligand-specific heme conformational change to accommodate the substrate binding. This substrate-induced modulation of the heme conformation may confer high product fidelity upon PGIS catalysis.     

Theoretical investigation of the origins of catalysis of a retro-Diels-Alder reaction by antibody 10F11

The antibody 10F11 catalyzes a retro-Diels-Alder reaction that forms HNO. Deductions about the mechanism of catalysis were made by Reymond, Baumann et al. from X-ray crystal structures and from kinetic measurements for mutated antibodies. We report a study of these reactions with quantum mechanical methods and a study of the substrate and transition state binding to the active site of the antibody 10F11 using density functional theory and empirical docking methods. We have quantitated the likely contributions to catalysis of three residues identified as possible causes of catalysis: Trp H104, Phe H101, and Ser H100. Trp H104 can make a significant contribution to catalysis through dispersive interactions (pi-stacking aromatic-aromatic stabilization). On its own, Phe H101 makes only a small contribution to catalysis. When both aromatic residues are present, they act cooperatively and can make greater contributions to catalysis than expected for each residue alone. Ser H100 and the backbone NH of Phe H101 are expected to act through hydrogen bonding to speed up the reaction, but our calculations suggest that they make only a small contribution to catalysis. Reymond's studies suggest that the hydrogen-bonding network may be mediated through a water molecule in the binding site.     

Mobility of the active site bound paraoxon and sarin in zinc-phosphotriesterase by molecular dynamics simulation and quantum chemical calculation.

The kinetic data published on phosphotriesterase (PTE), with various complexed metals, clearly indicates that the P=O and P=S bonds of phosphotriester and thiophosphotriester substrates, respectively, are strongly polarized by one or both of the active site complexed metal ions. However, this observation is not consistent with the three-dimensional X-ray crystal structure of zinc-substituted PTE with active site bound substrate analogue diethyl 4-methylbenzylphosphonate. In this structure, the distance between the phosphoryl oxygen and the nearest zinc is 3.4 A, a distance too large to afford strong polarization. In the present paper, the geometry and mobility of various PTE active site-substrate complexes are examined by performing both molecular dynamics (MD) simulations and quantum mechanical calculations. Two known substrates are considered, paraoxon and sarin, although their turnover rates vary about 100-fold. The results indicate that PTE forms a complex with either substrate in which the phosphoryl oxygen becomes strongly coordinated with the less buried zinc atom. It is shown that the geometry of the active site is changed when the protein is immersed in a water bath and relaxed by MD. The most substantial conformational change is the opening of the gateway in a pocket where the location of the leaving group is expected. The opening is observed for the pure enzyme as well as for the enzyme/substrate complexes and it ranges from 11 to 18 A. It is also shown that the pockets, in which the substrate substituents are localized, exhibit different flexibility and interact with the substrate with coordinated conformational adjustments.     

Probing Substrate Diffusion in Interstitial MOF Chemistry with Kinetic Isotope Effects

Metal–organic frameworks (MOFs) have garnered substantial interest as platforms for site‐isolated catalysis. Efficient diffusion of small‐molecule substrates to interstitial lattice‐confined catalyst sites is critical to leveraging unique opportunities of these materials as catalysts. Understanding the rates of substrate diffusion in MOFs is challenging, and few in situ chemical tools are available to evaluate substrate diffusion during interstitial MOF chemistry. Herein, we demonstrate nitrogen atom transfer (NAT) from a lattice‐confined Ru₂ nitride to toluene to generate benzylamine. We use the comparison of the intramolecular deuterium kinetic isotope effect (KIE), determined for amination of a partially deuterated substrate, with the intermolecular KIE, determined by competitive amination of a mixture of perdeuterated and undeuterated substrates, to establish the relative rates of substrate diffusion and interstitial chemistry. We anticipate that the developed KIE‐based experiments will contribute to the development of porous materials for group‐transfer catalysis.