Comparative assessment of theoretical methods for the determination of geometrical properties in biological zinc complexes

The Role of the CuCl Active Complex in the Stereoselectivity of the Salt-Induced Peptide Formation Reaction: Insights from Density Functional Theory Calculations

The salt-induced peptide formation (SIPF) reaction is a prebiotically plausible mechanism for the spontaneous polymerization of amino acids into peptides on early Earth. Experimental investigations of the SIPF reaction have found that in certain conditions, the l enantiomer is more reactive than the d enantiomer, indicating its potential role in the rise of biohomochirality. Previous work hypothesized that the distortion of the CuCl active complex toward a tetrahedral-like structure increases the central chirality on the Cu ion, which amplifies the inherent parity-violating energy differences between l- and d-amino acid enantiomers, leading to stereoselectivity. Computational evaluations of this theory have been limited to the protonated-neutral l + l forms of the CuCl active complex. Here, density functional theory methods were used to compare the energies and geometries of the homochiral (l + l and d + d) and heterochiral (l + d) CuCl-amino acid complexes for both the positive-neutral and neutral-neutral forms for alanine, valine, and proline. Significant energy differences were not observed between different chiral active complexes (i.e., d + d, l + l vs. l + d), and the distortions of active complexes between stereoselective systems and non-selective systems were not consistent, indicating that the geometry of the active complex is not the primary driver of the observed stereoselectivity of the SIPF reaction.

Roles of the Active Site Zn(II) and Residues in Substrate Discrimination by Threonyl-tRNA Synthetase: An MD and QM/MM Investigation

Threonyl-tRNA synthetase (ThrRS) is a Zn(II) containing enzyme that catalyzes the activation of threonine and its subsequent transfer to the cognate tRNA. This process is accomplished with remarkable fidelity, with ThrRS being able to discriminate its cognate substrate from similar analogues such as serine and valine. Molecular dynamics (MD) simulations and hybrid quantum mechanics/molecular mechanics (QM/MM) methods have been used to elucidate the role of Zn(II) in the aminoacylation mechanism of ThrRS. More specifically, the role of Zn(II) and active site residues in ThrRS's ability to discriminate between its cognate substrate l-threonine and the noncognate l-serine, l-valine, and d-threonine has been examined. The present results suggest that a role of the Zn(II) ion, with its Lewis acidity, is to facilitate deprotonation of the side chain hydroxyl groups of the aminoacyl moieties of cognate Thr-AMP and noncognate Ser-AMP substrates. In their deprotonated forms, these substrates are able to adopt a conformation preferable for aminoacyl transfer from aa-AMP onto the Ado-3'OH of the tRNA^Thr cosubstrate. Relative to the neutral substrates, when the substrates are first deprotonated with the assistance of the Zn(II) ion, the barrier for the rate-limiting step is decreased significantly by 42.0 and 39.2 kJ mol^-1 for l-Thr-AMP and l-Ser-AMP, respectively. An active site arginyl also plays a key role in stabilizing the buildup of negative charge on the substrate's bridging phosphate oxygen during the mechanism. For the enantiomeric substrate analogue d-Thr-AMP, product formation is highly disfavored, and as a result, the reverse reaction has a very low barrier of 16.0 kJ mol^-1.

The influence of zinc(ii) on thioredoxin/glutathione disulfide exchange: QM/MM studies to explore how zinc(ii) accelerates exchange in higher dielectric environments

QM/MM calculations on thiolate disulfide exchange reveal that a polar, 4-centered, cyclic transition state is formed when Zn(ii) is present, helping to explain faster exchange rates in higher dielectric solvents for metal-assisted exchange.

Theoretical study on the chemical mechanism of enoyl-CoA hydratase and the form of inhibitor binding

Enoyl-CoA hydratase (ECH) catalyzes the second step in the vital β-oxidation pathway of fatty acid metabolism. This enzyme catalyzes the syn-addition of a water molecule across the double bond of 4-(N,N-dimethylamino) cinnamoyl-CoA (DAC-CoA). In this work, the reaction mechanisms of ECH were investigated using the density functional theory (DFT) methods. The different protonation states in which the important residues Glu164 and Glu144 are either neutral or ionized were considered. Four models of the active site were designed based on the X-ray crystal structure of the enzyme. The calculations gave strong support to the proposed mechanism and confirmed that both Glu164 and Glu144 are in a deprotonated state in the reaction mechanism of ECH. In addition, we constructed a model of the active site with the inhibitor acetoacetyl-CoA based on the crystal structure. Caomparison of the calculated energy barriers showed that binding of the keto-enol form of the inhibitor is more reasonable than that of the di-keto form in the inhibition process. Moreover, acetoacetyl-CoA was found to exhibit a keto-enol tautomerism when it acts as an inhibitor in the reaction. The present theoretical results indicated that both residues Glu164 and Glu144 are unprotonated in ECH with the substrate bound, while only Glu164 is unprotonated when the inhibitor binds ECH.

Comparative analysis of the performance of commonly available density functionals in the determination of geometrical parameters for copper complexes

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.

High- and low-spin Fe(III) complexes of various AGE inhibitors

Density functional theory calculations [CPCM/UM06/6-31+G(d,p)] were used to elucidate the structures and relative stability of Fe(III) complexes with various ligands that inhibit the formation of advanced glycation end products (AGEs) or iron overloaded disease (viz. aminoguanidine, pyridoxamine, LR-74, Amadori compounds, and ascorbic acid). EDTA was used as the free energy reference ligand. The distorted neutral octahedral complex containing one iron atom and three molecules of pyridoxamine [Fe(PM)(3)] was found to be the most stable. The stability of the complexes decreases in the following chelate sequence: pyridoxamine, Amadori complex, aminoguanidine, LR inhibitor, and ascorbic acid.

The Accuracy of Density Functional Theory in the Description of Cation-π and π-Hydrogen Bond Interactions

Cation-π and π-hydrogen bond interactions are ubiquitous in protein folding, molecular recognition, and ligand-receptor associations. As such systems are routinely studied at the DFT level, it becomes essential to understand the underlying accuracy of the plethora of density functionals currently available for the description of these interactions. For that purpose, we carried out theoretical calculations on two small model systems (benzene-Na(+) and benzene-H2O) that represent a paradigm for those intermolecular interactions and systematically tested 46 density functionals against the results of high-level post-HF methods, ranging from MP2 to extrapolated CCSD(T)/CBS. A total of 13 basis sets were also tested to examine the convergence of the interaction energy with basis set size. The convergence was surprisingly fast, with deviations below 0.2 kcal/mol for double-ζ polarized basis sets with diffuse functions. Concerning functional benchmarking, the Truhlar group functionals were particularly well suited for the description of the π-hydrogen bond interactions. In the case of cation-π interactions, there was not a clear correlation between accuracy and functional sophistication. Despite the large number of functionals predicting interaction energies within chemical accuracy (five for π-hydrogen bond and 20 for cation-π interactions), not a single functional has shown chemical accuracy in both cases. Moreover, if we calculate the average error for these two interactions, only two density functionals resulted in an average error below 1.0 kcal/mol (M06 and HCTH, with average errors of 0.6 and 0.8 kcal/mol). The obtained results serve as a guide for future computer simulations on this kind of system.

Effect of metal Ions (Ni²⁺, Cu²⁺ and Zn²⁺) and water coordination on the structure of L-phenylalanine, L-tyrosine, L-tryptophan and their zwitterionic forms

Methods of quantum chemistry have been applied to double-charged complexes involving the transition metals Ni²⁺, Cu²⁺ and Zn²⁺ with the aromatic amino acids (AAA) phenylalanine, tyrosine and tryptophan. The effect of hydration on the relative stability and geometry of the individual species studied has been evaluated within the supermolecule approach. The interaction enthalpies, entropies and Gibbs energies of nine complexes Phe•M, Tyr•M, Trp•M, (M = Ni²⁺, Cu²⁺ and Zn²⁺) were determined at the Becke3LYP density functional level of theory. Of the transition metals studied the bivalent copper cation forms the strongest complexes with AAAs. For Ni²⁺and Cu²⁺ the most stable species are the NO coordinated cations in the AAA metal complexes, Zn²⁺cation prefers a binding to the aromatic part of the AAA (complex II). Some complexes energetically unfavored in the gas-phase are stabilized upon microsolvation.

Phosphate mono- and diesterase activities of the trinuclear zinc enzyme nuclease P1--insights from quantum chemical calculations

Nuclease P1 is a trinuclear zinc enzyme that catalyzes the hydrolysis of single-stranded DNA and RNA. Density functional calculations are used to elucidate the reaction mechanism of this enzyme with a model of the active site designed on the basis of the X-ray crystal structure. 2-Tetrahydrofuranyl phosphate and methyl 2-tetrahydrofuranyl phosphate substrates are used to explore the phosphomonoesterase and phosphodiesterase activities of this enzyme, respectively. The calculations reveal that for both activities, a bridging hydroxide performs an in-line attack on the phosphorus center, resulting in inversion of the configuration. Simultaneously, the P-O bond is cleaved, and Zn2 stabilizes the negative charge of the leaving alkoxide anion and assists its departure. All three zinc ions, together with Arg48, provide electrostatic stabilization to the penta-coordinated transition state, thereby lowering the reaction barrier.

Effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+ and Zn2+) and water coordination on the structure and properties of L-histidine and zwitterionic L-histidine

Interactions between metal ions and amino acids are common both in solution and in the gas phase. The effect of metal ions and water on the structure of L-histidine is examined. The effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+ and Zn2+) and water on structures of His·M(H2O)m, m=0.1 complexes have been determined theoretically employing density functional theories using extended basis sets. Of the five stable complexes investigated the relative stability of the gas-phase complexes computed with DFT methods (with one exception of K+ systems) suggest metallic complexes of the neutral L-histidine to be the most stable species. The calculations of monohydrated systems show that even one water molecule has a profound effect on the relative stability of individual complexes. Proton dissociation enthalpies and Gibbs energies of L-histidine in the presence of the metal cations Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+ and Zn2+ were also computed. Its gas-phase acidity considerably increases upon chelation. Of the Lewis acids investigated, the strongest affinity to L-histidine is exhibited by the Cu2+ cation. The computed Gibbs energies ΔG are negative, span a rather broad energy interval (from -130 to -1,300 kJ/mol), and upon hydration are appreciably lowered.

Comparative analysis of the performance of commonly available density functionals in the determination of geometrical parameters for zinc complexes

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).

Effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water coordination on the structure and properties of L-arginine and zwitterionic L-arginine

Interactions between metal ions and amino acids are common both in solution and in the gas phase. The effect of metal ions and water on the structure of L-arginine is examined. The effects of metal ions (Li(+), Na(+), K(+), Mg(2+), Ca(2+), Ni(2+), Cu(2+), and Zn(2+)) and water on structures of Arg x M(H2O)m , m = 0, 1 complexes have been determined theoretically by employing the density functional theories (DFT) and using extended basis sets. Of the three stable complexes investigated, the relative stability of the gas-phase complexes computed with DFT methods (with the exception of K(+) systems) suggests metallic complexes of the neutral L-arginine to be the most stable species. The calculations of monohydrated systems show that even one water molecule has a profound effect on the relative stability of individual complexes. Proton dissociation enthalpies and Gibbs energies of arginine in the presence of the metal cations Li(+), Na(+), K(+), Mg(2+), Ca(2+), Ni(2+), Cu(2+), and Zn(2+) were also computed. Its gas-phase acidity considerably increases upon chelation. Of the Lewis acids investigated, the strongest affinity to arginine is exhibited by the Cu(2+) cation. The computed Gibbs energies DeltaG(o) are negative, span a rather broad energy interval (from -150 to -1500 kJ/mol), and are appreciably lowered upon hydration.

Computational enzymatic catalysis

Computational methodologies are playing increasingly important roles in elucidating and presenting the complete and detailed mechanisms of enzymatic reactions because of their capacity to determine and characterize intermediates and transition states from both structural and energetics points of view, independent of their reduced lifetimes and without interfering with the natural reactional flux. These features are turning the field into an active and interesting area of research, involving a diverse range of studies, mostly directed at understanding the ways in which enzymes function under certain circumstances and predicting how they will behave under others. The accuracy of the computational data obtained for a given mechanistic hypothesis depends essentially on three mutually exclusive factors: the accuracy of the Hamiltonian of the reaction mechanism, consideration of the modulating aspect of the enzyme's structure in the energetics of the active center, and consideration of the enzyme's conformational fluctuations and dynamics. Although, unfortunately, it is impossible at present to optimize these crucial factors simultaneously, the success of any enzymatic mechanistic study depends on the level of equilibrium achieved among them. Different authors adopt different solutions, and this Account summarizes the most favored, with emphasis placed on our own preferences. Another crucial aspect in computational enzymatic catalysis is the model used in the calculations. Our aim is to build the simplest model that captures the essence of the catalytic power of an enzyme, allowing us to apply the highest possible theoretical level and minimize accidental errors. The choice is, however, far from obvious, ranging from simple models containing tens of atoms up to models of full enzymes plus solvent. Many factors underlie the choice of an appropriate model; here, examples are presented of very different modeling strategies that have been employed to obtain meaningful results. One particular case study, that of enzyme ribonucleotide reductase (RNR), a radical enzyme that catalyzes the reduction of ribonucleotides into deoxyribonucleotides, is one of the examples illustrating how the successive increase of the system's size does not dramatically change the thermodynamics and kinetics of the reaction. The values obtained and presented speak for themselves in that the only ones that are distinctly different are those calculated using an exceedingly small model, which omitted the amino acids that establish hydrogen bonds with the reactive unit of the substrate. This Account also describes our computational analysis of the mechanism of farnesyltransferase, a heterodimeric zinc metalloenzyme that is currently one of the most fascinating targets in cancer research. We focus on the present methodologies that we have been using, our models and understanding of the problem, and the accuracy of results and associated problems within this area of research.