Computational simulation of the lifetime of the methoxymethyl cation in water. A simple model for a glycosyl cation: when is an intermediate an intermediate?

A two-dimensional free-energy surface is constructed for transfer of the methoxymethyl cation between two water molecules. These atoms are treated quantum mechanically within a box of >1000 classical solvent water molecules, and the molecular dynamics of the whole system is considered at 300 K. This provides a simple model for glycosyl transfer in water. The best surface obtained (MPWB1K/6-31+G(d,p) corrected AMI/TIP3P) contains a shallow free-energy well corresponding to an oxacarbenium ion intermediate in a stepwise mechanism. Molecular dynamics analysis at three temperatures leads to a classical estimate of the lifetime of the methoxymethyl cation in water; when quantum corrections for vibrational zero-point energy are included, the lifetime is estimated to be 1 ps. This result is in complete agreement with the best experimental estimate and suggests that computational simulation is a reliable tool for elucidation of glycosyl-transfer mechanisms in enzymes and whether these involve glycosyl cations as intermediates

Does glycosyl transfer involve an oxacarbenium intermediate? Computational simulation of the lifetime of the methoxymethyl cation in water

2D free-energy surfaces for transfer of the methoxymethyl cation between two water molecules are constructed from molecular dynamics (MD) simulations in which these atoms are treated quantum-mechanically within a box of 1030 classical solvent water molecules at 300 K. This provides a simple model for glycosyl transfer in water. The AM1/TIP3P surfaces with 2D-spline corrections at either MPWB1K/6-31+G(d,p) or MP2/6-31+G(d,p) contain a shallow free-energy well corresponding to an oxacarbenium ion intermediate in a DN*AN mechanism. MD analysis at three temperatures leads to a classical estimate of the lifetime of the methoxymethyl cation in water; when quantum corrections for vibrational zero-point energy are included, the lifetime is estimated to be about 1 ps, in agreement with the best experimental estimate. This suggests that computational simulation, with appropriate high-level correction, is a reliable tool to obtain detailed and reliable mechanistic descriptions for glycosidases. In view of the importance of developing improved anti-influenza drugs, simulations of sialidases that considered both sialyl oxacarbenium ion and covalent sialyl-enzyme as possible intermediates could provide particular insight

Unraveling the role of protein dynamics in dihydrofolate reductase catalysis

Protein dynamics have controversially been proposed to be at the heart of enzyme catalysis, but identification and analysis of dynamical effects in enzyme-catalyzed reactions have proved very challenging. Here, we tackle this question by comparing an enzyme with its heavy (15N, 13C, 2H substituted) counterpart, providing a subtle probe of dynamics. The crucial hydride transfer step of the reaction (the chemical step) occurs more slowly in the heavy enzyme. A combination of experimental results, quantum mechanics/molecular mechanics simulations, and theoretical analyses identify the origins of the observed differences in reactivity. The generally slightly slower reaction in the heavy enzyme reflects differences in environmental coupling to the hydride transfer step. Importantly, the barrier and contribution of quantum tunneling are not affected, indicating no significant role for “promoting motions” in driving tunneling or modulating the barrier. The chemical step is slower in the heavy enzyme because protein motions coupled to the reaction coordinate are slower. The fact that the heavy enzyme is only slightly less active than its light counterpart shows that protein dynamics have a small, but measurable, effect on the chemical reaction rate

Chemical Ligation and Isotope Labeling to Locate Dynamic Effects during Catalysis by Dihydrofolate Reductase

Chemical ligation has been used to alter motions in specific regions of dihydrofolate reductase from E. coli and to investigate the effects of localized motional changes on enzyme catalysis. Two isotopic hybrids were prepared; one with the mobile N-terminal segment containing heavy isotopes (2H, 13C, 15N) and the remainder of the protein with natural isotopic abundance, and the other one with only the C-terminal segment isotopically labeled. Kinetic investigations indicated that isotopic substitution of the N-terminal segment affected only a physical step of catalysis, whereas the enzyme chemistry was affected by protein motions from the C-terminal segment. QM/MM studies support the idea that dynamic effects on catalysis mostly originate from the C-terminal segment. The use of isotope hybrids provides insights into the microscopic mechanism of dynamic coupling, which is difficult to obtain with other studies, and helps define the dynamic networks of intramolecular interactions central to enzyme catalysis

Minimization of dynamic effects in the evolution of dihydrofolate reductase

Protein isotope labeling is a powerful technique to probe functionally important motions in enzyme catalysis and can be applied to investigate the conformational dynamics of proteins. Previous investigations have indicated that dynamic coupling is detrimental to catalysis by dihydrofolate reductase (DHFR) from the mesophile Escherichia coli (EcDHFR). Comparison of DHFRs from organisms adapted to survive at a wide range of temperatures suggests that dynamic coupling in DHFR catalysis has been minimized during evolution; it arises from reorganizational motions needed to facilitate charge transfer events. Contrary to the behaviour observed for the DHFR from the moderate thermophile Geobacillus stearothermophilus (BsDHFR), the chemical transformation catalyzed by the cold-adapted bacterium Moritella profunda (MpDHFR) is only weakly affected by protein isotope substitutions at low temperatures, but the isotopically substituted enzyme is a substantially inferior catalyst at higher, non-physiological temperatures. QM/MM studies revealed that this behaviour is caused by the enzyme’s structural sensitivity to temperature changes, which enhances unfavorable dynamic coupling at higher temperatures by promoting additional recrossing trajectories on the transition state dividing surface. We postulate that these motions are minimized by fine-tuning DHFR flexibility through optimization of the free energy surface of the reaction, such that a nearly static reaction-ready configuration with optimal electrostatic properties is maintained under physiological conditions

Loop Interactions during Catalysis by Dihydrofolate Reductase fromMoritella profunda

Dihydrofolate reductase (DHFR) is often used as a model system to study the relation between protein dynamics and catalysis. We have studied a number of variants of the cold-adapted DHFR from Moritella profunda (MpDHFR), in which the catalytically important M20 and FG loops have been altered, and present a comparison with the corresponding variants of the wellstudied DHFR from Escherichia coli (EcDHFR). Mutations in the M20 loop do not affect the actual chemical step of transfer of hydride from reduced nicotinamide adenine dinucleotide phosphate to the substrate 7,8-dihydrofolate in the catalytic cycle in either enzyme; they affect the steady state turnover rate in EcDHFR but not in MpDHFR. Mutations in the FG loop also have different effects on catalysis by the two DHFRs. Despite the two enzymes most likely sharing a common catalytic cycle at pH 7, motions of these loops, known to be important for progression through the catalytic cycle in EcDHFR, appear not to play a significant role in MpDHFR

Protein motions and dynamic effects in enzyme catalysis

The role of protein motions in promoting the chemical step of enzyme catalysed reactions remains a subject of considerable debate. Here, a unified view of the role of protein dynamics in dihydrofolate reductase catalysis is described. Recently the role of such motions has been investigated by characterising the biophysical properties of isotopically substituted enzymes through a combination of experimental and computational analyses. Together with previous work, these results suggest that dynamic coupling to the chemical coordinate is detrimental to catalysis and may have been selected against during DHFR evolution. The full catalytic power of Nature's catalysts appears to depend on finely tuning protein motions in each step of the catalytic cycle

Search for top squarks in the four-body decay mode with single lepton final states in proton-proton collisions at s \sqrt{s} = 13 TeV

A search for the pair production of the lightest supersymmetric partner of the top quark, the top squark (t∼1), is presented. The search targets the four-body decay of the t∼1, which is preferred when the mass difference between the top squark and the lightest supersymmetric particle is smaller than the mass of the W boson. This decay mode consists of a bottom quark, two other fermions, and the lightest neutralino (χ∼01), which is assumed to be the lightest supersymmetric particle. The data correspond to an integrated luminosity of 138 fb−1 of proton-proton collisions at a center-of-mass energy of 13 TeV collected by the CMS experiment at the CERN LHC. Events are selected using the presence of a high-momentum jet, an electron or muon with low transverse momentum, and a significant missing transverse momentum. The signal is selected based on a multivariate approach that is optimized for the difference between m(t∼1) and m(χ∼01). The contribution from leading background processes is estimated from data. No significant excess is observed above the expectation from standard model processes. The results of this search exclude top squarks at 95% confidence level for masses up to 480 and 700 GeV for m(t∼1) − m(χ∼01) = 10 and 80 GeV, respectively

Search for Higgs Boson Decay to a Charm Quark-Antiquark Pair in Proton-Proton Collisions at √s = 13 TeV

A search for the standard model Higgs boson decaying to a charm quark-antiquark pair, H→c¯c, produced in association with a leptonically decaying V (W or Z) boson is presented. The search is performed with proton-proton collisions at √s=13  TeV collected by the CMS experiment, corresponding to an integrated luminosity of 138  fb−1. Novel charm jet identification and analysis methods using machine learning techniques are employed. The analysis is validated by searching for Z→c¯c in VZ events, leading to its first observation at a hadron collider with a significance of 5.7 standard deviations. The observed (expected) upper limit on σ(VH)B(H→c¯c) is 0.94 (0.50+0.22−0.15)pb at 95% confidence level (C.L.), corresponding to 14 (7.6+3.4−2.3) times the standard model prediction. For the Higgs-charm Yukawa coupling modifier, κc, the observed (expected) 95% C.L. interval is 1.1<|κc|<5.5 (|κc|<3.4), the most stringent constraint to date