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.

Nucleases catalyze the cleavage of phosphodiester bonds in nucleic acids using a range of metal cofactors. Although it is well accepted that many nucleases rely on two metal ions, the...

Glutaminyl-tRNA synthetase (GlnRS) catalyzes the aminoacylation of glutamine to the corresponding tRNA. However, most bacteria and all archaea lack GlnRS and thus an indirect noncanonical aminoacylation is required. With the assistance of a non-discriminating version of Glutamyl-tRNA synthetases (ND-GluRS) the tRNA is misaminoacylated by glutamate. In this study, we have computationally investigated the aminoacylation mechanism in GlnRS and ND-GluRS employing Molecular Dynamics (MD) simulations, Quantum Mechanics (QM) cluster and Quantum Mechanics/Molecular Mechanics (QM/MM) calculations. Our investigations demonstrated the feasibility of a water-mediated, substrate-assisted catalysis pathway with rate limiting steps occurring at energy barriers of 25.0 and 25.4 kcal mol for GlnRS and ND-GluRS, respectively. A conserved lysine residue participates in a second proton transfer to facilitate the departure of the adenosine monophosphate (AMP) group. Thermodynamically stable (-29.9 and -9.3 kcal mol for GlnRS and ND-GluRS) product complexes are obtained only when the AMP group is neutral.

Polyethylene terephthalate (PET), the most extensively used plastic, is one of the significant contributors to global plastic pollution. Enzymatic biodegradation of PET using different hydrolases has been previously reported as a promising biodegradation strategy for closed-loop recycling. Among the different hydrolases known to depolymerize PET to its soluble building blocks, the PETase and cutinase family of enzymes have notable PET biodegradation activities. In fact, they exhibit different thermostabilities and efficiencies in hydrolyzing PET polyesters despite sharing high structural similarities. Herein, we employed quantum mechanics/molecular mechanics calculations to identify the key factors necessary for efficient PET hydrolysis. Our results show that in both PETase and cutinase (Tfcut2 as a model system), the PET hydrolysis reaction pathway proceeds through a multi-step process with rate-limiting steps having energy barriers of ∼18.0 and ∼20 kcal/mol for PETase and TfCut2, respectively, which agrees well with the experimental data. A deeper inspection of the structural complexes revealed that the bent conformation adopted by PET and the tighter H-bond interaction between the catalytic triad residues, mediated by the unique disulfide bridge, contribute to the lower barrier (i.e., better catalytic performance) of PETase. The intrinsic molecular features identified in this work will also be useful for rational engineering of more efficient cutinases for PET hydrolysis.

New research and development efforts using computational chemistry in studying an assessment of the validity of different quantum chemical methods to describe the molecular and electronic structures of some corrosion inhibitors were introduced. The standard and the highly accurate CCSD method with 6-311++G(d,p), ab initio calculations using the HF/6-31G++(d,p) and MP2 with 6-311G(d,p), 6-31++G(d,p), and 6-311++G(2df,p) methods as well as DFT method at the B3LYP, BP86, B3LYP*, M06L, and M062x/6-31G++(d,p) basis set level were performed on some triazole derivatives and sulfur containing compounds used as corrosion inhibitors. Quantum chemical parameters, such as the energy of the highest occupied molecular orbital energy (E(HOMO)), the energy of the lowest unoccupied molecular orbital energy (E(LUMO)), energy gap (ΔE), dipole moment (μ), sum of total negative charges (TNC), chemical potential (Pi), electronegativity (χ), hardness (η), softness (σ), local softness (s), Fukui functions (f (+),f (-)), electrophilicity (ω), the total energy change (∆E(T)) and the solvation energy (S.E), were calculated. Furthermore, the accuracy and the applicability of these methods were estimated relative to the highest accuracy and standard CCSD with 6-311++G(d,p) method. Good correlations between the quantum chemical parameters and the corresponding inhibition efficiency (IE%) were found.