Search

The domain-based local pair natural orbital coupled-cluster with single, double, and perturbative triples excitation (DLPNO-CCSD(T)) method was employed to portray the activation and reaction energies of four ubiquitous enzymatic reactions, and its performance was confronted to CCSD(T)/complete basis set (CBS) to assess its accuracy and robustness in this specific field. The DLPNO-CCSD(T) results were also confronted to those of a set of density functionals (DFs) to understand the benefit of implementing this technique in enzymatic quantum mechanics/molecular mechanics (QM/MM) calculations as a second QM component, which is often treated with DF theory (DFT). On average, the DLPNO-CCSD(T)/aug-cc-pVTZ results were 0.51 kcal·mol -1 apart from the canonic CCSD(T)/CBS, without noticeable biases toward any of the reactions under study. All DFs fell short to the DLPNO-CCSD(T), both in terms of accuracy and robustness, which suggests that this method is advantageous to characterize enzymatic reactions and that its use in QM/MM calculations, either alone or in conjugation with DFT, in a two-region QM layer (DLPNO-CCSD(T):DFT), should enhance the quality and faithfulness of the results., (© 2020 Wiley Periodicals LLC.)

In this study, a set of 50 transition-metal complexes of Cu(I) and Cu(II), were used in the evaluation of 18 density functionals in geometry determination. In addition, 14 different basis sets were considered, including four commonly used Pople's all-electron basis sets; four basis sets including popular types of effective-core potentials: Los Alamos, Steven-Basch-Krauss, and Stuttgart-Dresden; and six triple-ζ basis sets. The results illustrate the performance of different methodological alternatives for the treatment of geometrical properties in relevant copper complexes, pointing out Double-Hybrid (DH) and Long-range Correction (LC) Generalized Gradient Approximation (GGA) methods as better descriptors of the geometry of the evaluated systems. These however, are associated with a computational cost several times higher than some of the other methods employed, such as the M06 functional, which has also demonstrated a comparable performance. Regarding the basis sets, 6-31+G(d) and 6-31+G(d,p) were the best performing approaches. In addition, the results show that the use of effective-core potentials has a limited impact, in terms of the accuracy in the determination of metal-ligand bond-lengths and angles in our dataset of copper complexes. Hence, these could become a good alternative for the geometrical description of these systems, particularly CEP-121G and SDD basis sets, if one is considering larger copper complexes where the computational cost could be an issue., (Copyright © 2013 Wiley Periodicals, Inc.)

Protein-protein binding is one of the critical events in biology, and knowledge of proteic complexes three-dimensional structures is of fundamental importance for the biochemical study of pharmacologic compounds. In the past two decades there was an emergence of a large variety of algorithms designed to predict the structures of protein-protein complexes--a procedure named docking. Computational methods, if accurate and reliable, could play an important role, both to infer functional properties and to guide new experiments. Despite the outstanding progress of the methodologies developed in this area, a few problems still prevent protein-protein docking to be a widespread practice in the structural study of proteins. In this review we focus our attention on the principles that govern docking, namely the algorithms used for searching and scoring, which are usually referred as the docking problem. We also focus our attention on the use of a flexible description of the proteins under study and the use of biological information as the localization of the hot spots, the important residues for protein-protein binding. The most common docking softwares are described too., (Copyright 2009 Wiley Periodicals, Inc.)

A set of 44 Zinc-ligand bond-lengths and of 60 ligand-metal-ligand bond angles from 10 diverse transition-metal complexes, representative of the coordination spheres of typical biological Zn systems, were used to evaluate the performance of a total of 18 commonly available density functionals in geometry determination. Five different basis sets were considered for each density functional, namely two all-electron basis sets (a double-zeta and triple-zeta formulation) and three basis sets including popular types of effective-core potentials: Los Alamos, Steven-Basch-Krauss, and Stuttgart-Dresden. The results show that there are presently several better alternatives to the popular B3LYP density functional for the determination of Zn-ligand bond-lengths and angles. BB1K, MPWB1K, MPW1K, B97-2 and TPSS are suggested as the strongest alternatives for this effect presently available in most computational chemistry software packages. In addition, the results show that the use of effective-core potentials (in particular Stuttgart-Dresden) has a very limited impact, in terms of accuracy, in the determination of metal-ligand bond-lengths and angles in Zinc-complexes, and is a good and safe alternative to the use of an all-electron basis set such as 6-31G(d) or 6-311G(d,p)., ((c) 2009 Wiley Periodicals, Inc.)

The enzymes of the thioredoxin family fulfill a wide range of physiological functions. Although they possess a similar CXYC active site motif, with identical environment and stereochemical properties, the redox potential and pK(a) of the cysteine pair varies widely across the family. As a consequence, each family member promotes oxidation or reduction reactions, or even isomerization reactions. The analysis of the three-dimensional structures gives no clues to identify the molecular source for the different active site properties. Therefore, we carried out a set of quantum mechanical calculations in active site models to gain more understanding on the elusive molecular-level origin of the differentiation of the properties across the family. The obtained results, together with earlier quantum mechanical calculations performed in our laboratories, gave rise to a consistent line of evidence, which points to the fact that both active site cysteines play an important role in the differentiation. In contrary to what was assumed, differentiation is not achieved through a different stabilization of the solvent exposed cysteine but, instead, through a fine tuning of the nucleophilicity of both active site cysteines. Reductant enzymes have both cysteine thiolates poorly stabilized, oxidant proteins have both cysteine thiolates highly stabilized, and isomerases have one thiolate (solvent exposed) poorly stabilized and the other (buried) thiolate highly stabilized. The feasibility of shifting the chemical equilibrium toward oxidation, reduction, or isomerization only through subtle electrostatic effects is quite unusual, and it relies on the inherent thermoneutrality of the catalytic steps carried out by a set of chemically equivalent entities all of which are cysteine thiolates. Such pattern of stabilization/destabilization, detected in our calculations is fully consistent with the observed physiological roles of this family of enzymes., (2008 Wiley Periodicals, Inc.)

Glycoside hydrolases catalyze the breaking of the glycosidic bond. This type of bond fashioned between two monosaccharides is very stable, and the polymers created are involved in multiple cellular processes, being crucial to life. In this article, computational methods were used to study the first step of the mechanism of reaction of retaining glycoside hydrolases in atomic detail. The systems modeled included a simplified reaction center and a small substrate/inhibitor. Using DFT calculations we were able to corroborate and provide molecular-level detail to the dissociative mechanism proposed in the literature. The role of the hydrogen bridge between the nucleophile and the C(2)--OH group of the ring was also investigated. Therefore, we concluded that this bridge is responsible for lowering the activation barrier by 5.1 kcal mol(-1) with functional BB1K/6-311+G(2d,2p), and the absence of the bridge explains, at least in part, the inhibitory effect of fluoro-substituted glycosides in the -2 position. The hydrogen bridge could also be involved in favoring the ring distortion verified in the transition state, and the dissociative character of the reaction mechanism. Using the NBO method, point atomic charges were calculated. In the transition state, the positive charge generated in the sugar ring is distributed nearly equally between the anomeric carbon and the ring oxygen, through a partial double bond involving the two atoms., (Copyright 2008 Wiley Periodicals, Inc.)

Farnesyltransferase (FTase), an interesting zinc metaloenzyme, has been the subject of great attention in anticancer research over the last decade. However, despite the major accomplishments in the field, some very pungent questions on the farnesylation mechanism still persist. In this study, the authors have analyzed a mechanistic paradox that arises from the existence of several contradicting and inconclusive experimental evidence regarding the existence of direct coordination between the active-site zinc cation and the thioether from the farnesylated peptide product, which include UV-vis spectroscopy data on a Co(2+)-substituted FTase, two X-ray crystallographic structures of the FTase-product complex, and extended X-ray absorption fine structure results. Using high-level theoretical calculations on two models of different sizes, and QM/MM calculations on the full enzyme, the authors have shown that the farnesylated product is Zn coordinated, and that a subsequent step where this Zn bond is broken is coherent with the available kinetic results. Furthermore, an explanation for the contradicting experimental evidence is suggested.

Alanine scanning mutagenesis of protein-protein interfacial residues can be applied to a wide variety of protein complexes to understand the structural and energetic characteristics of the hot-spots. Binding free energies have been estimated with reasonable accuracy with empirical methods, such as Molecular Mechanics/Poisson-Boltzmann surface area (MM-PBSA), and with more rigorous computational approaches like Free Energy Perturbation (FEP) and Thermodynamic Integration (TI). The main objective of this work is the development of an improved methodological approach, with less computational cost, that predicts accurately differences in binding free energies between the wild-type and alanine mutated complexes (DeltaDeltaG(binding)). The method was applied to three complexes, and a mean unsigned error of 0.80 kcal/mol was obtained in a set of 46 mutations. The computational method presented here achieved an overall success rate of 80% and an 82% success rate in residues for which alanine mutation causes an increase in the binding free energy > 2.0 kcal/mol (warm- and hot-spots). This fully atomistic computational methodological approach consists in a computational Molecular Dynamics simulation protocol performed in a continuum medium using the Generalized Born model. A set of three different internal dielectric constants, to mimic the different degree of relaxation of the interface when different types of amino acids are mutated for alanine, have to be used for the proteins, depending on the type of amino acid that is mutated. This method permits a systematic scanning mutagenesis of protein-protein interfaces and it is capable of anticipating the experimental results of mutagenesis, thus guiding new experimental investigations.

Thioredoxin superfamily members share a considerable degree of structural similarity, with a conserved CX(i)X(j)C motif at the active site, where C stand for two cysteines that alternate between a reduced thiol and oxidized disulfide states, and X(i)and X(j) are two amino acids different in each family member. Despite these similarities, they display very different redox potentials and pKas for the active site dithiol, and fulfill different physiological roles. Thioredoxin, for example, promotes the reduction of disulfide bonds, while DsbA promotes their oxidation in prokaryotic cells. The factors that promote these differences are still not fully understood. However, it is generally accepted that the different stabilities of the redox active disulfide bond depends on the degree of stabilization, in the reduced state, of the thiolate of one of the active site cysteines (nucleophilic cysteine). In this work we have used QM/MM methods to compare and characterize the active site dithiols of both enzymes, and to shed some light on the structural features responsible for the large differences in pKa and redox potential between two homologous enzymes, thioredoxin and DsbA. We have also analyzed the main factors pointed out in the literature as responsible for their different properties. We obtained the value of 4.5 for pKa difference (DeltapKa) between the nucleophilic cysteines of both enzymes, which is in excellent agreement with most of the experimental values. Additionally, we found that the principal differentiating factor responsible for this observed DeltapKa are the alpha2-alpha helices, which greatly contribute to the mentioned value, by stabilizing the DsbA thiolate in a much greater extend than the thioredoxin thiolate. A double mutation of the conserved residues Asp26 and Lys57, in thioredoxin, and Glu24 Lys58, in DsbA, by alanines did not change the DeltapKa value; this supports the hypothesis that these residues are not involved in the differentiation of the properties of the active centre dithiol. However, we found out that these residues are important for the stabilization of the nucleophilic thiolate. The X(i) and X(j) residues also do not seem to promote the stabilization of the thiolates. In fact, the corresponding double alanine mutants are more stable than the wild-type enzymes. However, these residues are involved in the differentiation between thioredoxin and DsbA, stabilizing the DsbA thiolate by a larger extent than the thioredoxin thiolate.

This article focuses on the first step of the catalytic mechanism for the reduction of ribonucleotides catalyzed by the enzyme Ribonucleotide Reductase (RNR). This corresponds to the activation of the substrate. In this work a large model of the active site region involving 130 atoms was used instead of the minimal gas phase models used in previous works. The ONIOM method was employed to deal with such a large system. The results gave additional information, which previous small models could not provide, allowing a much clearer evaluation of the role of the enzyme in this step. Enzyme-substrate interaction energies, specific transition state stabilization, and substrate steric strain energies were obtained. It was concluded that the transition state is stabilized in 4.0 kcal/mol by specific enzyme-substrate interactions. However, this stabilization is cancelled by the cost in conformational energy for the enzyme to adopt the transition state geometry; the overall result is that the enzyme machinery does not lead to a rate enhancement in this step. It was also found that the substrate binds to the active site with almost no steric strain, emphasizing the complementarity and specificity of the RNR active site for nucleotide binding. The main role of the enzyme at the very beginning of the catalytic cycle was concluded to be to impose stereospecifity upon substrate activation and to protect the enzyme radical from the solvent, rather than to be an reaction rate enhancement., ((c) 2004 Wiley Periodicals, Inc)

Gemcitabine (2',2'-difluoro-2'-deoxycytidine, dFdC) is a very promising anticancer drug, already approved for clinical use in three therapeutic indications. It is metabolized intracellularly to 5'-diphosphate (dFdCDP), which is known to be a potent inhibitor of ribonucleotide reductase (RNR). Although several nucleotide analogs show in vitro capacity of RNR inactivation, none has shown the in vivo efficacy of dFdCDP. Accordingly, the experimental data suggests that its mechanism of inhibition is different from the other known RNR suicide inhibitors. Enzyme inhibition in the absence of reductive species leads to complete loss of the essential radical in subunit R2, and formation of a new nucleotide-based radical. Interestingly, however, the presence of the reductants does not prevent inhibition--the radical is not lost but the targeted subunit of RNR becomes R1, which is inactivated possibly by alkylation. We have conducted a theoretical study, which led us to the first proposal of a possible mechanism for RNR inhibition by dFdCDP in the absence of reductants. This mechanism turned out to be very similar to the natural substrate reduction pathway and only deviates from the natural course after the formation of the well-known disulphide bridge. This deviation is caused precisely by the F atom in the beta-face, only present in this inhibitor. The essential radical in R2 is lost, and so is the enzyme catalytic activity. The nucleotide-based radical that constitutes the end product of our mechanism has been suggested in the literature as a possible candidate for the one detected experimentally. In fact, all experimental data available has been reproduced by the theoretical calculations performed here., (Copyright 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1286-1294, 2004)

2'-azido-2'-deoxyribonucleoside 5'-diphosphates are mechanism-based inhibitors of Ribonucleotide Reductase. Considerable effort has been made to elucidate their mechanism of inhibition, which is still controversial and not fully understood. Previous studies have detected the formation of a radical intermediate when the inhibitors interact with the enzyme, and several authors have proposed possible structures for this radical. We have conducted a theoretical study of the possible reactions involved, which allowed us to identify the structure of the new radical among the several proposals. A new reactional path is also proposed that is the most kinetically favored to yield this radical and ultimately inactivate the enzyme. The energetic involved in this mechanism, both for radical formation and radical decay, as well as the calculated Hyperfine Coupling Constants for the radical intermediate, are in agreement with the correspondent experimental values. This mechanistic alternative is fully coherent with remaining experimental data., (Copyright 2003 Wiley Periodicals, Inc.)