Transition Metals as Probes of Metal Cofactors in Nucleic Acid Biochemistry

Assessment of Density Functional Theory Methods for the Structural Prediction of Transition and Post-Transition Metal-Nucleic Acid Complexes

Understanding the structure of metal-nucleic acid systems is important for many applications such as the design of new pharmaceuticals, metal detection platforms, and nanomaterials. Herein, we explore the ability of 20 density functional theory (DFT) functionals to reproduce the crystal structure geometry of transition and post-transition metal-nucleic acid complexes identified in the Protein Data Bank and Cambridge Structural Database. The environmental extremes of the gas phase and implicit water were considered, and analysis focused on the global and inner coordination geometry, including the coordination distances. Although gas-phase calculations were unable to describe the structure of 12 out of the 53 complexes in our test set regardless of the DFT functional considered, accounting for the broader environment through implicit solvation or constraining the model truncation points to crystallographic coordinates generally afforded agreement with the experimental structure, suggesting that functional performance for these systems is likely due to the models rather than the methods. For the remaining 41 complexes, our results show that the reliability of functionals depends on the metal identity, with the magnitude of error varying across the periodic table. Furthermore, minimal changes in the geometries of these metal-nucleic acid complexes occur upon use of the Stuttgart-Dresden effective core potential and/or inclusion of an implicit water environment. The overall top three performing functionals are ωB97X-V, ωB97X-D3(BJ), and MN15, which reliably describe the structure of a broad range of metal-nucleic acid systems. Other suitable functionals include MN15-L, which is a cheaper alternative to MN15, and PBEh-3c, which is commonly used in QM/MM calculations of biomolecules. In fact, these five methods were the only functionals tested to reproduce the coordination sphere of Cu²⁺-containing complexes. For metal-nucleic acid systems that do not contain Cu²⁺, ωB97X and ωB97X-D are also suitable choices. These top-performing methods can be utilized in future investigations of diverse metal-nucleic acid complexes of relevance to biology and material science.

Effect of substitution inert metal complexes on calcineurin

As a substitute for M(H2O)2+6, Co(NH3)3+6 was found to activate calcineurin with para-nitrophenyl phosphate as substrate. Kinetics for calcineurin catalyzed hydrolysis of para-nitrophenyl phosphate at pH 7.0 with Mn2+, Mg2+, Co2+, and Co(NH3)3+6 were compared. Although kcat and Km were different with the metals, values of kcat/Km were nearly identical for Mn2+ and Mg2+, but lower for Co2+ and Co(NH3)3+6. The concentration of each metal providing half-maximal activation, designated Kact, was evaluated as 15.9 mM for Co(NH3)3+6, compared to Kact = 0.17 mM for Mn2+ and Co2+ and 6.3 mM for Mg2+, respectively. Comparing kcat/Kcat showed that Co(NH3)3+6 was a 170-fold poorer activator of calcineurin than was Mn2+, but only 1.5-fold poorer than Mg2+. Activation by Co(NH3)3+6 indicated that activation of calcineurin by exogenous metal ions can proceed via an outer coordination sphere reaction mechanism with no requirement for the direct coordination of substrate by metal. Because Co(NH3)3+6 was found to support calcineurin activity, the related compound [Co-(ethylenediamine)3]3+ (or Co(en)3+3) was tested as a possible activator. Co(en)3+3 did not support calcineurin activity but did inhibit calcineurin. Co(en)3+3 showed competitive inhibition kinetics with either Mn2+ or pNPP as the varied ligand and the other at a fixed, subsaturating concentration. Inorganic phosphate was used as a known competitive inhibitor to pNPP (B. L. Martin and D. J. Graves, J. Biol. Chem. 261, 14545-14550, 1986) and showed uncompetitive inhibition with Mn2+ as the varied ligand. These patterns are consistent with the mechanism of ligand binding to calcineurin being ordered with metal preceding substrate. Prior formation of a metal-substrate complex was not required for association with calcineurin.

Metallobiochemistry of the magnesium ion. Characterization of the essential metal-binding site in Escherichia coli ribonuclease H

Ribonuclease H (Escherichia coli) contains one strong magnesium-binding site, as determined by metal-titration experiments monitored by high field 1H-NMR and also by direct titration calorimetry. Kinetic and thermodynamic parameters were evaluated by 25Mg-NMR and were as follows: dissociation constant Kd, approximately 60 +/- 10 microM; activation free energy delta G*, approximately 49.8 +/- 0.9 kJ; on/off-rate for magnesium binding Kon, approximately 1.8 x 10(8) M-1 s-1, koff, approximately 1.1 x 10(4) s-1; quadrupole coupling constant chi B, 1.2 +/- 0.2 MHz. The dissociation constant was independently determined by standard analysis of 1H chemical shifts in magnesium-titration experiments and by microcalorimetry (Kd approximately 200 +/- 20 microM). Cobalt hexaamine, which also activates RNase H [Jou, R. & Cowan, J. A. (1991) J. Am. Chem. Soc. 113, 6685-6686], appears to bind at the same location as Mg2+(aqueous). Assignments of C2H and C4H protons to specific histidine residues have been made by two-dimensional correlated spectroscopy experiments. Direct 25Mg-NMR pH titrations show that an ionizable residue (pKa approximately 5.8), most likely one of the carboxylates in the active site, influences magnesium binding. On the basis of the magnesium coordination chemistry elucidated herein, recent proposals on active-site chemistry are critically assessed and general physicochemical aspects of magnesium-binding sites on proteins and enzymes are discussed.

Sequence selective coordination of Mg2+(aq) to DNA

Thermodynamic parameters for magnesium binding to a series of DNA molecules of defined sequence have been evaluated by 25Mg NMR spectroscopy. These results demonstrate that G/C-DNA binds Mg2+(aq) up to 40 to 100-fold more strongly than A/T-DNA, i.e., coordination of Mg2+(aq) to G/C-DNA is ca. 2.1-2.7 kcal (mole Mg2+)-1 more stable relative to A/T-DNA. Activation free energies [delta G* approximately (12.7-13.3) x 10(3) kcal] and exchange rates [k(ex) approximately (0.5-3.0) x 10(3)s-1] were estimated by variable temperature experiments. The low value of the quadrupole coupling constants (chi B = 0.2-0.6 MHz) is indicative of outer-sphere coordination by Mg(H2O)6(2+).