Despite remarkable stability, the phosphodiester bond of nucleic acids is hydrolytically cleaved in critical biological processes. Although this reaction is commonly accepted to take place via a two-metal-assisted mechanism, recent experimental evidence suggests that several enzymes use a single-metal ion, but the precise catalytic mechanism is unknown. In the present work, we employ a multiscale computational approach to decipher the phosphodiester cleavage mechanism for this unique pathway by focusing on the human APE1 repair enzyme, which catalyzes the incision of phosphodiester bonds adjacent to DNA lesions. To resolve ambiguity in the literature regarding the role of the single-metal (Mg­(II)) center, several catalytic mechanisms were carefully examined. Our predicted preferred hydrolysis pathway proceeds in two steps via a pentacovalent phosphorane intermediate in the absence of substrate ligation to Mg­(II), with a rate-limiting barrier (19.3 kcal/mol) in close agreement with experiment (18.3 kcal/mol). In this mechanism, D210 promotes catalysis by activating water for nucleophilic attack at the 5′-phosphate group with respect to the damaged site. Subsequently, a Mg­(II)-bound water triggers leaving group departure by neutralizing the 3′-hydroxyl of the neighboring nucleotide. Consistent with experimental kinetic and mutational data, several other active site residues (N212, Y171, and H309) play multiple roles throughout the reaction to facilitate this challenging chemistry. In addition to revealing previously unknown mechanistic features of the APE1 catalyzed reaction, our work sets the stage for exploring the phosphodiester bond cleavage catalyzed by other single-metal-dependent enzymes, as well as different pharmaceutical and biotechnological applications.

Promoting selective interactions between a nucleophile and electrophilic dye in complex environments is a central goal in nucleophilic chemosensor development. Commonly employed dyes are hemicyanines containing either the N-methylbenzothiazolium (Btz) or the N-methyl-3,3-dimethylindolium (Ind) acceptors. The dyes are related to α,β-unsaturated carbonyls and contain two sites of reactivity (C2 vs C4) with the C2-site directly attached to the quaternary nitrogen possessing greater electrophilicity. We demonstrate the regioselectivity between reactions of sodium thiomethoxide (NaSMe) with two electrophilic hemicyanine dyes bearing Btz (1) or Ind (2) in dipolar aprotic solvent–water mixtures. Adduct complexation was followed by NMR spectroscopy, and structures were optimized in the gas phase to estimate relative adduct stability. The key results include finding a preference for thiolate attachment at the C4-site to generate an enamine adduct with no evidence for attachment at the more electrophilic C2-position. Equilibration between NaSMe and water also affords NaOH that displays a thermodynamic preference for C2-attachment. Dye 1 containing the Btz moiety exhibits greater selectivity for the thiolate addition, with dye 2 being more reactive toward adventitious water to generate OH-adducts. Our data affords diagnostic 1H/13C NMR adduct signals, regioselectivity for various dye/nucleophile combinations, and suggests use of the Btz acceptor for direct thiolate detection.

Aminoacyl-tRNA synthetases (aaRSs) are central to a number of physiological processes, including protein biosynthesis. In particular, they activate and then transfer their corresponding amino acid to the cognate tRNA. This is achieved with a generally remarkably high fidelity by editing against incorrect standard and nonstandard amino acids. Using docking, molecular dynamics (MD), and hybrid quantum mechanical/molecular mechanics methods, we have investigated mechanisms by which methionyl-tRNA synthetase (MetRS) may edit against the highly toxic, noncognate, amino acids homocysteine (Hcy) and its oxygen analogue, homoserine (Hse). Substrate-assisted editing of Hcy-AMP in which its own phosphate acts as the mechanistic base occurs with a rate-limiting barrier of 98.2 kJ mol(-1). This step corresponds to nucleophilic attack of the Hcy side-chain sulfur at its own carbonyl carbon (CCarb). In contrast, a new possible editing mechanism is identified in which an active site aspartate (Asp259) acts as the base. The rate-limiting step is now rotation about the substrate's aminoacyl Cβ-Cγ bond with a barrier of 27.5 kJ mol(-1), while for Hse-AMP, the rate-limiting step is cleavage of the CCarb-OP bond with a barrier of 30.9 kJ mol(-1). A similarly positioned aspartate or glutamate also occurs in the homologous enzymes LeuRS, IleRS, and ValRS, which also discriminate against Hcy. Docking and MD studies suggest that at least in the case of LeuRS and ValRS, a similar editing mechanism may be possible.

Phosphodiester bond hydrolysis in nucleic acids is a ubiquitous reaction that can be facilitated by enzymes called nucleases, which often use metal ions to achieve catalytic function. While a two-metal-mediated pathway has been well established for many enzymes, there is growing support that some enzymes require only one metal for the catalytic step. Using human apurinic/ apyrimidinic endonuclease (APE1) as a prototypical example and cluster models, this study clarifies the impact of DFT functional, cluster model size, and implicit solvation on single-metal-mediated phosphodiester bond cleavage and provides insight into how to efficiently model this chemistry. Initially, a model containing 69 atoms built from a high-resolution X-ray crystal structure is used to explore the reaction pathway mapped by a range of DFT functionals and basis sets, which provides support for the use of standard functionals (M06-2X and B3LYP-D3) to study this reaction. Subsequently, systematically increasing the model size to 185 atoms by including additional amino acids and altering residue truncation points highlights that small models containing only a few amino acids or β carbon truncation points introduce model strains and lead to incorrect metal coordination. Indeed, a model that contains all key residues (general base and acid, residues that stabilize the substrate, and amino acids that maintain the metal coordination) is required for an accurate structural depiction of the one-metal-mediated phosphodiester bond hydrolysis by APE1, which results in 185 atoms. The additional inclusion of the broader enzyme environment through continuum solvation models has negligible effects. The insights gained in the present work can be used to direct future computational studies of other one-metal-dependent nucleases to provide a greater understanding of how nature achieves this difficult chemistry.