DNA cleavage by endonuclease I-DmoI: a QM/MM study and comparison with experimental data provide indications on the environmental effects

Unveiling the Cleavage Mechanism of an RNA Model Compound on the whole pH Scale: Computations Meet Experiments in the Determination of Reaction Rates

The knowledge of the mechanism of reactions occurring in solution is a primary research line both in the context of theoretical-computational chemistry and in the field of organic and bio-organic chemistry. Given the importance of the hydrolysis of nucleic acids in life-related phenomena, here we present a combined experimental and computational study on the cleavage of an RNA model compound. This phosphodiester features a cleavage rate strictly dependent on the pH with three different dependence domains. Such experimental evidence, highlighted by an in-depth kinetic investigation, unequivocally suggests a change in the reaction mechanism along the pH scale. In order to interpret the data and to explain the experimental behavior, we have applied a theoretical-computational procedure, involving a hybrid quantum/classical approach, able to model chemical reactions in complex environments, i. e. in solution. This study turns out to quantitatively reproduce the experimental data with accuracy and, in addition, provides useful mechanistic insight into the transesterification process of the investigated compound. The study indicates that the cleavage can occur through anA N D N ${A_N D_N }$ , anA N + D N ${A_N + D_N }$ , and aD N A N ${D_N A_N }$ mechanism depending on the pH values.

A Theoretical-Computational Study of Phosphodiester Bond Cleavage Kinetics as a Function of the Temperature

The hydrolysis of the phosphodiester bond is an important chemical reaction involved in several biological processes. Here, we study the cleavage of this bond by means of a theoretical-computational method in a model system, the dineopentyl phosphate. By such an approach, we reconstructed the kinetics and related thermodynamics of this chemical reaction along an isochore. In particular, we evaluated the kinetic constants of all the reaction steps within a wide range of temperatures, mostly corresponding to conditions where no experimental measures are available due to the extremely slow kinetics. Our results, in good agreement with the experimental data, show the robustness of our theoretical-computational methodology which can be easily extended to more complex systems.

Molecular Modeling of the Deamidation Reaction in Solution: A Theoretical-Computational Study

In this work, a theoretical-computational method is applied to study the deamidation reaction, a critical post-translational modification in proteins, using a simple model molecule in solution. The method allows one to comprehensively address the environmental effect, thereby enabling one to accurately derive the kinetic rate constants for the three main steps of the deamidation process. The results presented, in rather good agreement with the available experimental data, underline the necessity for a rigorous treatment of environmental factors and a precise kinetic model to correctly assess the overall kinetics of the deamidation reaction.

Modelling Complex Bimolecular Reactions in a Condensed Phase: The Case of Phosphodiester Hydrolysis

(1) Background: the theoretical modelling of reactions occurring in liquid phase is a research line of primary importance both in theoretical-computational chemistry and in the context of organic and biological chemistry. Here we present the modelling of the kinetics of the hydroxide-promoted hydrolysis of phosphoric diesters. (2) Method: the theoretical-computational procedure involves a hybrid quantum/classical approach based on the perturbed matrix method (PMM) in conjunction with molecular mechanics. (3) Results: the presented study reproduces the experimental data both in the rate constants and in the mechanistic aspects (C-O bond vs. O-P bond reactivity). The study suggests that the basic hydrolysis of phosphodiesters occurs through a concerted ANDN mechanism, with no formation of penta-coordinated species as reaction intermediates. (4) Conclusions: the presented approach, despite the approximations, is potentially applicable to a large number of bimolecular transformations in solution and therefore leads the way to a fast and general method to predict the rate constants and reactivities/selectivities in complex environments.

Theoretical Evaluation of Sulfur-Based Reactions as a Model for Biological Antioxidant Defense

Sulfur-containing amino acids, Methionine (Met) and Cysteine (Cys), are very susceptible to Reactive Oxygen Species (ROS). Therefore, sulfur-based reactions regulate many biological processes, playing a key role in maintaining cellular redox homeostasis and modulating intracellular signaling cascades. In oxidative conditions, Met acts as a ROS scavenger, through Met sulfoxide formation, while thiol/disulfide interchange reactions take place between Cys residues as a response to many environmental stimuli. In this work, we apply a QM/MM theoretical-computational approach, which combines quantum-mechanical calculations with classical molecular dynamics simulations to estimate the free energy profile for the above-mentioned reactions in solution. The results obtained, in good agreement with experimental data, show the validity of our approach in modeling sulfur-based reactions, enabling us to study these mechanisms in more complex biological systems.

Effects of Environmental and Electric Perturbations on the pKa of Thioredoxin Cysteine 35: A Computational Study

Here we present a theoretical-computational study dealing with the evaluation of the pKa of the Cysteine residues in Thioredoxin (TRX) and in its complex with the Thioredoxin-interacting protein (TXNIP). The free energy differences between the anionic and neutral form of the Cysteine 32 and 35 have been evaluated by means of the Perturbed Matrix Method with classical perturbations due to both the environment and an exogenous electric field as provided by Molecular Dynamics (MD) simulations. The evaluation of the free energies allowed us to show that the effect of the perturbing terms is to lower the pKa of Cysteine 32 and Cysteine 35 with respect to the free amino-acid. On the other hand, in the complex TRX-TXNIP, our data show an enhanced stabilization of the neutral reduced form of Cys 35. These results suggest that external electric stimuli higher than 0.02 V/nm can modulate the Cysteine pKa, which can be connected to the tight regulation of the TRX acting as an antioxidant agent.