Juan Aranda

The Catalytic Mechanism of Carboxylesterases: A Computational Study

The catalytic mechanism of carboxylesterases (CEs, EC 3.1.1.1) is explored by computational means. CEs hydrolyze ester, amide, and carbamate bonds found in xenobiotics and endobiotics. They can also perform transesterification, a reaction important, for instance, in cholesterol homeostasis. The catalytic mechanisms with three different substrates (ester, thioester, and amide) have been established at the M06-2X/6-311++G**//B3LYP/6-31G* level of theory. It was found that the reactions proceed through a mechanism involving four steps instead of two as is generally proposed: (i) nucleophilic attack of serine to the substrate, forming the first tetrahedral intermediate, (ii) formation of the acyl-enzyme complex concomitant with the release of the alcohol product, (iii) nucleophilic attack of a water or alcohol molecule forming the second tetrahedral intermediate, and (iv) the release of the second product of the reaction. The results agree very well with the available experimental data and show that the hydrolytic and the transesterification reactions are competitive processes when the substrate is an ester. In all the other studied substrates (thioester or amide), the hydrolytic and transesterification process are less favorable and some of them might not even take place under in vivo conditions.

Dynamics and Reactivity in Thermus aquaticus N6-Adenine Methyltransferase

M.TaqI is a DNA methyltransferase from Thermus aquaticus that catalyzes the transfer of a methyl group from S-adenosyl-L-methionine to the N6 position of an adenine, a process described only in prokaryotes. We have used full atomistic classical molecular dynamics simulations to explore the protein-SAM-DNA ternary complex where the target adenine is flipped out into the active site. Key protein-DNA interactions established by the target adenine in the active site are described in detail. The relaxed structure was used for a combined quantum mechanics/molecular mechanics exploration of the reaction mechanism using the string method. According to our free energy calculations the reaction takes place through a stepwise mechanism where the methyl transfer precedes the abstraction of the proton from the exocyclic amino group. The methyl transfer is the rate-determining step, and the obtained free energy barrier is in good agreement with the value derived from the experimental rate constant. Two possible candidates to extract the leftover proton have been explored: a water molecule found in the active site and Asn105, a residue activated by the hydrogen bonds formed through the amide hydrogens. The barrier for the proton abstraction is smaller when Asn105 acts as a base. The reaction mechanisms can be different in other N6-DNA-methyltransferases, as determined from the exploration of the reaction mechanism in the Asn105Asp M.TaqI mutant.

Singlet Oxygen Attack on Guanine: Reactivity and Structural Signature within the B‐DNA Helix

Oxidatively generated DNA lesions are numerous and versatile, and have been the subject of intensive research since the discovery of 8-oxoguanine in 1984. Even for this prototypical lesion, the precise mechanism of formation remains elusive due to the inherent difficulties in characterizing high-energy intermediates. We have probed the stability of the guanine endoperoxide in B-DNA as a key intermediate and determined a unique activation free energy of around 6 kcal mol(-1) for the formation of the first C-O covalent bond upon the attack of singlet molecular oxygen ((1) O2 ) on the central guanine of a solvated 13 base-pair poly(dG-dC), described by means of quantum mechanics/molecular mechanics (QM/MM) simulations. The B-helix remains stable upon oxidation in spite of the bulky character of the guanine endoperoxide. Our modeling study has revealed the nature of the versatile (1) O2 attack in terms of free energy and shows a sensitivity to electrostatics and solvation as it involves a charge-separated intermediate.

Substrate promiscuity in DNA methyltransferase M.PvuII. A mechanistic insight

M.PvuII is a DNA methyltransferase from the bacterium Proteus vulgaris that catalyzes methylation of cytosine at the N4 position. This enzyme also displays promiscuous activity catalyzing methylation of adenine at the N6 position. In this work we use QM/MM methods to investigate the reaction mechanism of this promiscuous activity. We found that N6 methylation in M.PvuII takes place by means of a stepwise mechanism in which deprotonation of the exocyclic amino group is followed by the methyl transfer. Deprotonation involves two residues of the active site, Ser53 and Asp96, while methylation takes place directly from the AdoMet cofactor to the target nitrogen atom. The same reaction mechanism was described for cytosine methylation in the same enzyme, while the reversal timing, that is methylation followed by deprotonation, has been described in M.TaqI, an enzyme that catalyzes the N6-adenine DNA methylation from Thermus aquaticus. These mechanistic findings can be useful to understand the evolutionary paths followed by N-methyltransferases.

Theoretical Study of the Catalytic Mechanism of DNA-(N4-Cytosine)-Methyltransferase from the Bacterium Proteus vulgaris

In this paper the reaction mechanism for methylation of cytosine at the exocyclic N4 position catalyzed by M.PvuII has been explored by means of hybrid quantum mechanics/molecular mechanics (QM/MM) methods. A reaction model was prepared by placing a single cytosine base in the active site of the enzyme. In this model the exocyclic amino group of the base establishes hydrogen bond interactions with the hydroxyl oxygen atom of Ser53 and the carbonyl oxygen atom of Pro54. The reaction mechanism involves a direct methyl transfer from AdoMet to the N4 atom and a proton transfer from this atom to Ser53, which in turn transfers a proton to Asp96. Different timings for the proton transfers and methylation steps have been explored at the AM1/MM and B3LYP/MM levels including localization and characterization of stationary structures. At our best estimate the reaction proceeds by means of a simultaneous but asynchronous proton transfer from Ser53 to Asp96 and from N4 of cytosine to Ser53 followed by a direct methyl transfer from AdoMet to the exocyclic N4 of cytosine.

Modeling methods for studying post-translational and transcriptional modifying enzymes

Biological catalysis is a complex chemical process that involves not only electronic reorganization in the substrate but also the reorganization of the catalyst. This complexity is even larger in the case of post-transcriptional and post-translational modifications which may involve the interaction between two biomacromolecules. However, the development over the past decades of new computational methods and strategies is offering a detailed molecular picture of the catalytic event and a deep understanding of the mechanisms of chemical reactions in biological environments. Here we review the efforts made in the last years to model catalysis in post-transcriptional and post-translational processes. We stress on the advantages and problems of the different computational strategies and their applicability in different cases.

Regioselectivity of the OH Radical Addition to Uracil in Nucleic Acids. A Theoretical Approach Based on QM/MM Simulations

Oxidation of nucleic acids is ubiquitous in living beings under metabolic impairments and/or exposed to external agents such as radiation, pollutants, or drugs, playing a central role in the development of many diseases mediated by DNA/RNA degeneration. Great efforts have been devoted to unveil the molecular mechanisms behind the OH radical additions to the double bonds of nucleobases; however, the specific role of the biological environment remains relatively unexplored. The present contribution tackles the study of the OH radical addition to uracil from the gas phase to a full RNA macromolecule by means of quantum-chemistry methods combined with molecular dynamics simulations. It is shown that, in addition to the intrinsic reactivity of each position driven by the electronic effects, the presence of bridge water molecules intercalated into the RNA structure favors the addition to the C5 position of uracil in biological conditions. The results also suggest that diffusion of the OH radical does not play a relevant role in the regioselectivity of the reaction, which is mainly controlled at the chemical stage of the addition process.