Complexes of Pentahydrated Zn2+ with Guanine, Adenine, and the Guanine−Cytosine and Adenine−Thymine Base Pairs. Structures and Energies Characterized by Polarizable Molecular Mechanics and ab Initio Calculations

Facile Approach to Fabricate a Chemical Sensor Array Based on Nanocurcumin-Metal Ions Aggregates: Detection and Identification of DNA Nucleobases

Here, a three-channel absorbance sensor array based on the nanocurcumin-metal ion (NCur-MI) aggregates is designed for the detection and identification of deoxyribonucleic acid nucleobases (DNA NBs) for the first time. For this purpose, the binding affinities of some of MIs (i.e., Co2+, Cr3+, Cu2+, Fe2+, Fe3+, Hg2+, Mn2+, Ni2+, V3+, and Zn2+) to the NCur to induce the aggregation were evaluated under various experimental conditions. Further studies reveal that in the presence of DNA NBs, the aggregates of NCur-Co2+, NCur-Ni2+, and NCur-Zn2+ show the diverse absorbance responses to the deaggregation of NCur depending on the binding affinity of each of DNA NBs to the metal ions Co2+, Ni2+, and Zn2+. These responses are distinguishable from one another. Thus, clear differentiation among the DNA NBs is achieved by linear discriminant analysis and hierarchical clustering analysis to generate clustering maps. The discriminatory capacity of the sensor array for the identification of the DNA NBs is tested in the ranges of 2.4-16 and 5.6-10.4 μM. Furthermore, a mixed set of the DNA NBs was prepared for multivariate multicomponent analysis. Finally, the practicability of the sensor array is confirmed by the discrimination of the DNA NBs in an animal DNA sample. It should be noted that the proposed array is the first example to fabricate an NCur-based sensor array for the simultaneous detection of DNA NBs.

Binding of metal ions and water molecules to nucleic acid bases: the influence of water molecule coordination to a metal ion on water-nucleic acid base hydrogen bonds

The interactions of nucleic acid bases with non-coordinated and coordinated water molecules were studied by analyzing data in the Protein Data Bank (PDB) and by quantum chemical calculations. The analysis of the data in the crystal structures from the PDB indicates that hydrogen bonds involving oxygen or nitrogen atoms of nucleic acid bases and water molecules are shorter when water is bonded to a metal ion. These results are in agreement with the quantum chemical calculations on geometries and interaction energies of hydrogen bonds; the calculations on model systems show that hydrogen bonds of nucleic acid bases with water bonded to a metal ion are stronger than hydrogen bonds with non-coordinated water. These calculated values are similar to the strength of hydrogen bonds between nucleic acid bases. The results presented in this paper may be relevant to understand the role of water molecules and metal ions in the process of replication and stabilization of nucleic acids and also to understand the possible toxicity of metal ion interactions with nucleic acids.

Calcium and Zinc Differentially Affect the Structure of Lipid Membranes

Interactions of calcium (Ca²⁺) and zinc (Zn²⁺) cations with biomimetic membranes made of dipalmitoylphosphatidylcholine (DPPC) were studied by small angle neutron diffraction (SAND). Experiments show that the structure of these lipid bilayers is differentially affected by the two divalent cations. Initially, both Ca²⁺ and Zn²⁺ cause DPPC bilayers to thicken, while further increases in Ca²⁺ concentration result in the bilayer thinning, eventually reverting to having the same thickness as pure DPPC. The binding of Zn²⁺, on the other hand, causes the bilayers to swell to a maximum thickness, and the addition of more Zn²⁺ does not result in a further thickening of the membrane. Agreement between our results obtained using oriented planar membranes and those from vesicular samples implies that the effect of cations on bilayer thickness is the result of electrostatic interactions, rather than geometrical constraints due to bilayer curvature. This notion is further reinforced by MD simulations. Finally, the radial distribution functions reveal a strong interaction between Ca²⁺ and the phosphate oxygens, while Zn²⁺ shows a much weaker binding specificity.

Hydrogen-bonded proton transfer in the hydrated adenine–thymine anion

The proton transfer processes of microhydrated adenine–thymine anions are studied using density functional theory with the B3LYP method and DZP++ basis set. The microhydration effects on the geometrical structures, adsorption site, and the proton transfer reaction of the adenine–thymine anion are investigated. The site N10 atom of the adenine moiety has a larger proton affinity than the site O24 atom of thymine, which facilitates the proton H26 transfers from the N25 site of thymine to the N10 site of adenine. Therefore, the first single-proton transfer pathway (SPT1) is found for the all of the monohydrated adenine–thymine anions (AN4T)−·H2O, (AN13T)−·H2O, (ATO24)−·H2O, and (ATO28)−·H2O and tetrahydrated adenine–thymine anions (AT)−·4H2O. The proton H9 at the N7 site of adenine is also found to transfer to the O24 site of thymine for (AN4T)−·H2O and (AN13T)−·H2O in the gas phase. The double-proton transferred pathway is found when one water molecule interacts with the O28 atom of thymine. The reactant structures before the proton transfer are more stable than the product structures, and the structural changes mainly occur in thymine. The reaction energies of the microhydrated adenine–thymine anion in the gas phase and in the aqueous environment predict that the proton transfer process of the microhydrated adenine–thymine anion are more favorable in the gas phase than in aqueous solution.

Spectral study on the unique enhanced fluorescence of guanosine triphosphate by zinc ions

Binding effect of guanosine triphosphate (GTP) with metal ions is involved in many biologically important processes, and so its investigation has been one interesting research focus for many chemical and biochemical research groups. In this contribution, we presented the unique fluorescence recovery and enhancement of GTP induced by Zn(II) based on the spectrofluorometric method. When excited at 280 nm, GTP is hardly fluorescent at the alkaline condition. However, the presence of Zn(II) caused an obvious fluorescence emission of GTP at 346 nm, and the binding molar ratio between GTP and Zn(II) had been proved to be 1. The investigations of binding property of various nucleotides with metal ions demonstrated that this fluorescence recovery and enhancement of GTP with Zn(II) was highly specific, which could successfully discriminate GTP from other structurally similar nucleotides including GDP and GMP. Furthermore, similar fluorescence response of the bacterial alarmone ppGpp to Zn(II) had also been identified.

A ratiometric fluorescence recognition of guanosine triphosphate on the basis of Zn(II) complex of 1,4-bis(imidazol-1-ylmethyl) benzene

As one vital member among the family of phosphates, guanosine triphosphate (GTP) not only plays a very important role in many critical biological processes but also closely associates with definite pathological states. Based on the ratiometric fluorescence response of the zinc complex of 1,4-bis(imidazol-1-ylmethyl) benzene (bix) in this contribution, a highly selective recognition of GTP has been successfully developed. The fluorescence of bix-Zn(II) at 289 nm decreased in the presence of GTP with the appearance of one new emission band at 341 nm, resulting in ratiometric fluorescence changes with the concentration of GTP. With that, ratiometric fluorescence recognition for GTP could be effectively established, and so GTP could be successfully discriminated from other structurally similar anions, including ATP and PPi. Furthermore, bix-Zn(II) also has a ratiometric fluorescence response to DNA sequences containing guanine.

The reaction mechanism of type I phosphomannose isomerases: new information from inhibition and polarizable molecular mechanics studies

Type I phosphomannose isomerases (PMIs) are zinc-dependent metalloenzymes involved in the reversible isomerization of D-mannose 6-phosphate (M6P) and D-fructose 6-phosphate (F6P). 5-Phospho-D-arabinonohydroxamic acid (5PAH), an inhibitor endowed with nanomolar affinity for yeast (Type I) and Pseudomonas aeruginosa (Type II) PMIs (Roux et al., Biochemistry 2004; 43:2926-2934), strongly inhibits human (Type I) PMI (for which we report an improved expression and purification procedure), as well as Escherichia coli (Type I) PMI. Its K(i) value of 41 nM for human PMI is the lowest value ever reported for an inhibitor of PMI. 5-Phospho-D-arabinonhydrazide, a neutral analogue of the reaction intermediate 1,2-cis-enediol, is about 15 times less efficient at inhibiting both enzymes, in accord with the anionic nature of the postulated high-energy reaction intermediate. Using the polarizable molecular mechanics, sum of interactions between fragments ab initio computed (SIBFA) procedure, computed structures of the complexes between Candida albicans (Type I) PMI and the cyclic substrate β-D-mannopyranose 6-phosphate (β-M6P) and between the enzyme and the high-energy intermediate analogue inhibitor 5PAH are reported. Their analysis allows us to identify clearly the nature of each individual active site amino acid and to formulate a hypothesis for the overall mechanism of the reaction catalyzed by Type I PMIs, that is, the ring-opening and isomerization steps, respectively. Following enzyme-catalyzed ring-opening of β-M6P by zinc-coordinated water and Gln111 ligands, Lys136 is identified as the probable catalytic base involved in proton transfer between the two carbon atoms C1 and C2 of the substrate D-mannose 6-phosphate.

Quantum chemical studies of nucleic acids: can we construct a bridge to the RNA structural biology and bioinformatics communities?

In this feature article, we provide a side-by-side introduction for two research fields: quantum chemical calculations of molecular interaction in nucleic acids and RNA structural bioinformatics. Our main aim is to demonstrate that these research areas, while largely separated in contemporary literature, have substantial potential to complement each other that could significantly contribute to our understanding of the exciting world of nucleic acids. We identify research questions amenable to the combined application of modern ab initio methods and bioinformatics analysis of experimental structures while also assessing the limitations of these approaches. The ultimate aim is to attain valuable physicochemical insights regarding the nature of the fundamental molecular interactions and how they shape RNA structures, dynamics, function, and evolution.

The effect of CH3, F and NO2 substituents on the individual hydrogen bond energies in the adenine-thymine and guanine-cytosine base pairs

The substituent effects on the geometrical parameters and the individual hydrogen bond (HB) energies of base pairs such as X-adenine-thymine (X-A-T), X-thymine-adenine (X-T-A), X-guanine-cytosine (X-G-C), and X-cytosine-guanine (X-C-G) have been studied by the quantum mechanical calculations at the B3LYP and MP2 levels with the 6-311++G(d,p) basis set. The electron withdrawing (EW) substituents (F and NO(2)) increase the total binding energy (DeltaE) of X-G-C derivatives and the electron donating (ED) substituent (CH(3)) decreases it when they are introduced in the 8 and 9 positions of G. The effects of substituents are reversed when they are located in the 1, 5, and 6 positions of C, with exception of CH(3) in the 1 position and F in the 5 position, which in both cases the DeltaE value decreases negligibly small. With minor exceptions (X=8-CH(3), 8-F, and 9-NO(2)), both ED and EW substituents increase slightly the DeltaE values of X-A-T derivatives. The individual HB energies (E (HB)s) have been estimated using electron densities that calculated at the hydrogen bond critical points (HBCPs) by the atoms in molecules (AIM) method. Most of changes of individual HBs are in consistent with the ED/EW nature of substituents and the role of atoms entered H-bonding. The remarkable change is observed for NO(2) substituted derivative in each case.

Frequency and effect of the binding of Mg2+, Mn2+, and Co2+ ions on the guanine base in Watson-Crick and reverse Watson-Crick base pairs

We performed MP2 calculations to elucidate the structure and energetics of the Mg(2+), Mn(2+), and Co(2+) hexahydrated aquaions, and the effect of the metal binding to the N7 atom of (i) a single guanine, (ii) a guanine involved in a Watson-Crick pair, and (iii) a guanine involved in a reverse Watson-Crick base pair. Our comparative analysis of the three aquaions indicates a clear inverse correlation between the radius of the cation and the binding energy, that indeed increases in the order Mn(2+) < Co(2+) < Mg(2+). The trend in the binding energies of the pentahydrated cations to the N7 atom of the guanine is instead Mg(2+) < Mn(2+) < Co(2+), suggesting a rather different bonding scheme that, for the two transition metals, involves back-donation from the aromatic ring of the guanine to their empty d orbitals. In the gas phase, the three hydrated metals significantly stabilize both G-C base pair geometries, Watson-Crick and reverse Watson-Crick, we investigated. Inclusion of a continuous solvent model, however, remarkably reduces this additional stabilization, which becomes almost negligible in the case of the Mg(2+) cation coordinated to the guanine in the standard Watson-Crick geometry. Conversely, all three metal ions sensibly stabilize the reverse Watson-Crick geometry, also in water. Our results are supported by a screening of the structures available in the Protein Data Bank, which clearly indicates that the two transition metals we investigated have a tendency greater than Mg(2+) to coordinate to the N7 atom of guanines, and that there is no clear correlation between the number of guanines in experimental structures with a metal bound to N7 atom and their involvement in Watson-Crick base pairs.

Infrared absorption detection of metal ion-deoxyguanosine monophosphate binding: experimental and theoretical study

Metal ion interactions with nucleic acids attract attention because of the environmental and biological consequences. The formation of the complex is often monitored by the vibrational spectroscopy. To identify characteristic binding patterns and marker bands on a model DNA component, infrared absorption spectra of the deoxyguanosine monophosphate complexes with Na(+), Mg(2+), Ca(2+), Ni(2+), Cu(2+), Zn(2+), and Cd(2+) cations were recorded and interpreted on the basis of density-functional computations. The aqueous environment was simulated by continuum and combined continuum-explicit solvent models. For the binding to the N7 position of the guanine base, the computation predicted a characteristic frequency upshift and splitting of the 1578 cm(-1) band, which is in accord with available experimental data. Contrary to the expectation, the modeling suggests that the binding to the carbonyl group might not be detectable, as the metal causes smaller spectral changes if compared to the hydrogen-bound water molecules. The binding to the phosphate group causes significant spectral changes in the sugar-phosphate vibrating region ( approximately 800-1200 cm(-1)), but also notable frequency shifts of the carbonyl vibrations. The Cu(2+) and Zn(2+) ions induced the largest alterations in observed vibrational absorption, which corresponds to the calculated strong interaction energies in the N7-complexes and to previous experimental experience. Additional changes in the vibrational spectra of the copper complexes were observed under high metal concentration, corresponding to the simultaneous binding to the phosphate residue. The two-step Cu(2+) binding process was also confirmed by the microcalorimetry titration curve. The computations and combination of more techniques thus help us to assign and localize spectral changes caused by the metal ion binding to nucleic acids.

Molecular dynamic simulations of the metallo-beta-lactamase from Bacteroides fragilis in the presence and absence of a tight-binding inhibitor

The beta-lactam-based antibiotics are among the most prescribed and effective antibacterial agents. Widespread use of these antibiotics, however, has created tremendous pressure for the emergence of resistance mechanisms in bacteria. The most common cause of antibiotic resistance is bacterial production of actamases that efficiently degrade antibiotics. The metallo-beta-lactamases are of particular clinical concern due to their transference between bacterial strains. We used molecular dynamics (MD) simulations to further study the conformational changes that occur due to binding of an inhibitor to the dicanzinc metallo-beta-lactamase from Bacteroides fragilis. Our studies confirm previous findings that the major flap is a major source of plasticity within the active site, therefore its dynamic response should be considered in drug development. However, our results also suggest the need for care in using MD simulations in evaluating loop mobility, both due to relaxation times and to the need to accurately model the zinc active site. Our study also reveals two new robust responses to ligand binding. First, there are specific localized changes in the zinc active site--a local loop flip--due to ligand intercalation that may be critical to the function of this enzyme. Second, inhibitor binding perturbs the dynamics throughout the protein, without otherwise perturbing the enzyme structure. These dynamic perturbations radiate outward from the active site and their existence suggests that long-range communication and dynamics may be important in the activity of this enzyme.

Cu2+/+ cation coordination to adenine--thymine base pair. Effects on intermolecular proton-transfer processes

Intermolecular proton-transfer processes in the Watson & Crick adenine-thymine Cu+ and Cu2+ cationized base pairs have been studied using the density functional theory (DFT) methods. Cationized systems subject to study are those resulting from cation coordination to the main basic sites of the base pair, N7 and N3 of adenine and O2 of thymine. For Cu+ coordinated to N7 or N3 of adenine, only the double proton-transferred product is found to be stable, similarly to the neutral system. However, when Cu+ interacts with thymine, through the O2 carbonyl atom, the single proton transfer from thymine to adenine becomes thermodynamically spontaneous, and thus rare forms of the DNA bases may spontaneously appear. For Cu2+ cation, important effects on proton-transfer processes appear due to oxidation of the base pair, which stabilizes the different single proton-transfer products. Results for hydrated systems show that the presence of the water molecules interacting with the metal cation (and their mode of coordination) can strongly influence the ability of Cu2+ to induce oxidation on the base pair.

Anisotropic, Polarizable Molecular Mechanics Studies of Inter- and Intramolecular Interactions and Ligand-Macromolecule Complexes. A Bottom-Up Strategy

We present an overview of the SIBFA polarizable molecular mechanics procedure, which is formulated and calibrated on the basis of quantum chemistry (QC). It embodies nonclassical effects such as electrostatic penetration, exchange-polarization, and charge transfer. We address the issues of anisotropy, nonadditivity, and transferability by performing parallel QC computations on multimolecular complexes. These encompass multiply H-bonded complexes and polycoordinated complexes of divalent cations. Recent applications to the docking of inhibitors to Zn-metalloproteins are presented next, namely metallo-beta-lactamase, phosphomannoisomerase, and the nucleocapsid of the HIV-1 retrovirus. Finally, toward third-generation intermolecular potentials based on density fitting, we present the development of a novel methodology, the Gaussian electrostatic model (GEM), which relies on ab initio-derived fragment electron densities to compute the components of the total interaction energy. As GEM offers the possibility of a continuous electrostatic model going from distributed multipoles to densities, it allows an inclusion of short-range quantum effects in the molecular mechanics energies. The perspectives of an integrated SIBFA/GEM/QM procedure are discussed.

Binding of 5-phospho-D-arabinonohydroxamate and 5-phospho-D-arabinonate inhibitors to zinc phosphomannose isomerase from Candida albicans studied by polarizable molecular mechanics and quantum mechanics

Type I phosphomannose isomerase (PMI) is a Zn-dependent metalloenzyme involved in the isomerization of D-fructose 6-phosphate to D-mannose 6-phosphate. One of our laboratories has recently designed and synthesized 5-phospho-D-arabinonohydroxamate (5PAH), an inhibitor endowed with a nanomolar affinity for PMI (Roux et al., Biochemistry 2004, 43, 2926). By contrast, the 5-phospho-D-arabinonate (5PAA), in which the hydroxamate moiety is replaced by a carboxylate one, is devoid of inhibitory potency. Subsequent biochemical studies showed that in its PMI complex, 5PAH binds Zn(II) through its hydroxamate moiety rather than through its phosphate. These results have stimulated the present theoretical investigation in which we resort to the SIBFA polarizable molecular mechanics procedure to unravel the structural and energetical aspects of 5PAH and 5PAA binding to a 164-residue model of PMI. Consistent with the experimental results, our theoretical studies indicate that the complexation of PMI by 5PAH is much more favorable than by 5PAA, and that in the 5PAH complex, Zn(II) ligation by hydroxamate is much more favorable than by phosphate. Validations by parallel quantum-chemical computations on model of the recognition site extracted from the PMI-inhibitor complexes, and totaling up to 140 atoms, showed the values of the SIBFA intermolecular interaction energies in such models to be able to reproduce the quantum-chemistry ones with relative errors < 3%. On the basis of the PMI-5PAH SIBFA energy-minimized structure, we report the first hypothesis of a detailed view of the active site of the zinc PMI complexed to the high-energy intermediate analogue inhibitor, which allows us to identify active site residues likely involved in the proton transfer between the two adjacent carbons of the substrates.

An ab initio investigation of zinc chloro complexes

A series of geometry, frequency, and energy calculations of chloroaquazinc(II) complexes were carried out at up to the MP2/6-31+G* level. A thorough examination of all species up to and including hexacoordinate species, and with up to six chlorides, was carried out. The structures of the complexes are compared with experimental data where available. The solution chemistry of zinc(II) in the presence of chloride is discussed, and Raman spectra of zinc perchlorate with increasing amount of chloride are presented.

Representation of Zn(II) complexes in polarizable molecular mechanics. Further refinements of the electrostatic and short-range contributions. Comparisons with parallel ab initio computations

We present refinements of the SIBFA molecular mechanics procedure to represent the intermolecular interaction energies of Zn(II). The two first-order contributions, electrostatic (E(MTP)), and short-range repulsion (E(rep)), are refined following the recent developments due to Piquemal et al. (Piquemal et al. J Phys Chem A 2003, 107, 9800; and Piquemal et al., submitted). Thus, E(MTP) is augmented with a penetration component, E(pen), which accounts for the effects of reduction in electronic density of a given molecular fragment sensed by another interacting fragment upon mutual overlap. E(pen) is fit in a limited number of selected Zn(II)-mono-ligated complexes so that the sum of E(MTP) and E(pen) reproduces the Coulomb contribution E(c) from an ab initio Hartree-Fock energy decomposition procedure. Denoting by S, the overlap matrix between localized orbitals on the interacting monomers, and by R, the distance between their centroids, E(rep) is expressed by a S(2)/R term now augmented with an S(2)/R(2) one. It is calibrated in selected monoligated Zn(II) complexes to fit the corresponding exchange repulsion E(exch) from ab initio energy decomposition, and no longer as previously the difference between (E(c) + E(exch)) and E(MTP). Along with the reformulation of the first-order contributions, a limited recalibration of the second-order contributions was carried out. As in our original formulation (Gresh, J Comput Chem 1995, 16, 856), the Zn(II) parameters for each energy contribution were calibrated to reproduce the radial behavior of its ab initio HF counterpart in monoligated complexes with N, O, and S ligands. The SIBFA procedure was subsequently validated by comparisons with parallel ab initio computations on several Zn(II) polyligated complexes, including binuclear Zn(II) complexes as in models for the Gal4 and beta-lactamase metalloproteins. The largest relative error with respect to the RVS computations is 3%, and the ordering in relative energies of competing structures reproduced even though the absolute numerical values of the ab initio interaction energies can be as large as 1220 kcal/mol. A term-to-term identification of the SIBFA contributions to their ab initio counterparts remained possible even for the largest sized complexes.

Marked Variations of Dissociation Energy and H-Bond Character of the Guanine-Cytosine Base Pair Induced by One-Electron Oxidation and Li+ Cation Coupling

The variation of dissociation energy and H-bond character of the G-C cation and the Li-GC cation have been investigated by employing density functional theory (B3LYP) with the 6-31+G* basis set. The one-electron oxidation and the coupling of Li(+) to the guanine-cytosine base pair can strengthen the interaction between guanine and cytosine. The interaction of the cation Li(+) with guanine is attractive and is attributed to the polarization of the H-bonds between G-C that enhances G-C interaction. The cooperativity of the three H-bonds in the GC and Li-GC cations is different from that in the neutral GC base pair. The proton-transfer process between N(1) of the guanine and N(3) of the cytosine can occur in the GC cation and the Li-GC cation. The geometries of the transition state are out of plane, especially for the transition state of the Li-GC cation. The analysis of the activation energy for the proton-transfer process shows that the GC(+) before and after proton transfer can exist simultaneously in the gas phase, but for the Li-GC(+) system, the Li-GC(+) without proton transfer is the dominating species in the gas phase.

Accurate interaction energies of hydrogen-bonded nucleic acid base pairs

Hydrogen-bonded nucleic acids base pairs substantially contribute to the structure and stability of nucleic acids. The study presents reference ab initio structures and interaction energies of selected base pairs with binding energies ranging from -5 to -47 kcal/mol. The molecular structures are obtained using the RI-MP2 (resolution of identity MP2) method with extended cc-pVTZ basis set of atomic orbitals. The RI-MP2 method provides results essentially identical with the standard MP2 method. The interaction energies are calculated using the Complete Basis Set (CBS) extrapolation at the RI-MP2 level. For some base pairs, Coupled-Cluster corrections with inclusion of noniterative triple contributions (CCSD(T)) are given. The calculations are compared with selected medium quality methods. The PW91 DFT functional with the 6-31G basis set matches well the RI-MP2/CBS absolute interaction energies and reproduces the relative values of base pairing energies with a maximum relative error of 2.6 kcal/mol when applied with Becke3LYP-optimized geometries. The Becke3LYP DFT functional underestimates the interaction energies by few kcal/mol with relative error of 2.2 kcal/mol. Very good performance of nonpolarizable Cornell et al. force field is confirmed and this indirectly supports the view that H-bonded base pairs are primarily stabilized by electrostatic interactions.

Conformation-dependent intermolecular interaction energies of the triphosphate anion with divalent metal cations. Application to the ATP-binding site of a binuclear bacterial enzyme. A parallel quantum chemical and polarizable molecular mechanics investigation

We have explored the conformation-dependent interaction energy of the triphosphate moiety, a key constituent of ATP and GTP, with a closed-shell divalent cation, Zn2+, used as a probe. This was done using the SIBFA polarizable molecular mechanics procedure. We have resorted to a previously developed approach in which triphosphate is built out from its elementary constitutive fragments, and the intramolecular, interfragment, interaction energies are computed simultaneously with their intermolecular interactions with the divalent cation. This approach has enabled reproduction of the values of the intermolecular interaction energies from ab initio quantum-chemistry with relative errors <3%. It was extended to the complex of a nonhydrolyzable analog of ATP with the active site of a bacterial enzyme having two Mg2+ cations as cofactors. We obtained following energy-minimization a very close overlap of the ATP analog over its position from X-ray crystallography. For models of the ATP analog-enzyme complex encompassing up to 169 atoms, the values of the SIBFA interaction energies were found to match their DFT counterparts with relative errors of <2%.

Theoretical Calculation of the NMR Spin−Spin Coupling Constants and the NMR Shifts Allow Distinguishability between the Specific Direct and the Water-Mediated Binding of a Divalent Metal Cation to Guanine

The calculated intermolecular and intramolecular indirect NMR spin-spin coupling constants and NMR shifts were used for the discrimination between the inner-shell and the outer-shell binding motif of hydrated divalent cations Mg(2+) or Zn(2+) with a guanine base. The intermolecular coupling constants (1)J(X,O6) and (1)J(X,N7) (X = Mg(2+), Zn(2+)) can be unambiguously assigned to the specific inner-shell binding motif of the hydrated cation either with oxygen O6 or with nitrogen N7 of guanine. The calculated coupling constants (1)J(Mg,O6) and (1)J(Zn,O6) were 6.2 and -17.5 Hz, respectively, for the inner-shell complex of cation directly interacting with oxygen O6 of guanine. For the inner-shell coordination of the cation at nitrogen N7, the calculated coupling constants (1)J(Mg,N7) and (1)J(Zn,N7) were 5.6 and -36.5 Hz, respectively. When the binding of the cation is water-mediated, the coupling constant is zero. To obtain reliable shifts in NMR parameters, hydrated guanine was utilized as the reference state. The calculated change of NMR spin-spin coupling constants due to the hydration and coordination of the cation with guanine is caused mainly by the variation of Fermi-contact coupling contribution while the variation of diamagnetic spin-orbit, paramagnetic spin-orbit, and spin-dipolar coupling contributions is small. The change of s-character of guanine sigma bonding, sigma antibonding, and lone pair orbitals upon the hydration and cation coordination (calculated using the Natural Bond Orbital analysis) correlates with the variation of the Fermi-contact term. The calculated NMR shifts delta(N7) of -15.3 and -12.2 ppm upon the coordination of Mg(2+) and Zn(2+) ion are similar to the NMR shift of 19.6 ppm toward the high field measured by Tanaka for N7 of guanine upon the coordination of the Cd(2+) cation (Tanaka, Y.; Kojima, C.; Morita, E. H.; Kasai. Y.; Yamasaki, K.; Ono, A.; Kainosho, M.; Taira, K. J. Am. Chem. Soc. 2002, 124, 4595-4601). The present data indicate that measurements of NMR intermolecular coupling constants may be used to discriminate between the specific inner- and outer-shell binding of divalent cations to nucleobases in DNA and RNA.

Substrate binding to mononuclear metallo-β-lactamase from Bacillus cereus

Structure and dynamics of substrate binding (cefotaxime) to the catalytic pocket of the mononuclear zinc-beta-lactamase from Bacillus cereus are investigated by molecular dynamics simulations. The calculations, which are based on the hydrogen-bond pattern recently proposed by Dal Peraro et al. (J Biol Inorg Chem 2002; 7:704-712), are carried out for both the free and the complexed enzyme. In the resting state, active site pattern and temperature B-factors are in agreement with crystallographic data. In the complexed form, cefotaxime is accommodated into a stable orientation in the catalytic pocket within the nanosecond timescale, interacting with the enzyme zinc-bound hydroxide and the surrounding loops. The beta-lactam ring remains stable and very close to the hydroxide nucleophile agent, giving a stable representation of the productive enzyme-substrate complex.

The effect of metal cations on the phase behavior and hydration characteristics of phospholipid membranes

To characterize the specificity of ion binding to phospholipids in terms of headgroup structure, hydration and lyotropic phase behavior we studied 1-palmitoyl-2-oleoyl-phosphatidylcholine as a function of relative humidity (RH) at 25 degrees C in the presence and absence of Li+, Na+, K+, Be2+, Mg2+, Ca2+, Sr2+, Ba2+, Zn2+ and Cu2+ ions by means of infrared (IR) spectroscopy. All divalent cations and Li+ shift the gel-to-liquid crystalline phase transition towards bigger RH values indicating stabilization of the gel state. The observed shift correlates in a linearly fashion with the electrostatic solvation free energy for most of the ions in water that in turn, is inversely related to the ionic radius. This interesting result was interpreted in terms of the excess chemical potential of mixing of hydrated ions and lipids. Calcium, zinc and partially lithium, cause a positive deviation from the linear relationship. IR spectral analysis shows that the carbonyl groups become more accessible to the water in the presence of Mg2+, Ca2+, Sr2+ and Ba2+ probably because of their involvement into the hydration shell of the ions. In contrast, Be2+, Zn2+ and Cu2+ dehydrate the carbonyl groups at small and medium RH. The ability of the lipid to take up water is distinctly reduced in the presence of Zn2+ and, partially, of Cu2+ meaning that the headgroups have become less hydrophilic. The binding mode of Be2+ to lipid headgroups involves hydrolyzed water. Polarized IR spectra show that complex formation of the phosphate groups with divalent ions gives rise to conformational changes and immobilization of the headgroups. The results are discussed in terms of the lyotropic Hofmeister series and of fusogenic activity of the ionic species.

Electronic properties, hydrogen bonding, stacking, and cation binding of DNA and RNA bases

This review summarizes results concerning molecular interactions of nucleic acid bases as revealed by advanced ab initio quantum chemical (QM) calculations published in last few years. We first explain advantages and limitations of modern QM calculations of nucleobases and provide a brief history of this still rather new field. Then we provide an overview of key electronic properties of standard and selected modified nucleobases, such as their charge distributions, dipole moments, polarizabilities, proton affinities, tautomeric equilibria, and amino group hybridization. Then we continue with hydrogen bonding of nucleobases, by analyzing energetics of standard base pairs, mismatched base pairs, thio-base pairs, and others. After this, the nature of aromatic stacking interactions is explained. Also, nonclassical interactions in nucleic acids such as interstrand bifurcated hydrogen bonds, interstrand close amino group contacts, C [bond] H...O interbase contacts, sugar-base stacking, intrinsically nonplanar base pairs, out-of-plane hydrogen bonds, and amino-acceptor interactions are commented on. Finally, we overview recent calculations on interactions between nucleic acid bases and metal cations. These studies deal with effects of cation binding on the strength of base pairs, analysis of specific differences among cations, such as the difference between zinc and magnesium, the influence of metalation on protonation and tautomeric equlibria of bases, and cation-pi interactions involving nucleobases. In this review, we do not provide methodological details, as these can be found in our preceding reviews. The interrelation between advanced QM approaches and classical molecular dynamics simulations is briefly discussed.

Influence of sodium ions on the dynamics and structure of single-stranded DNA oligomers: a molecular dynamics study

The effects of sodium counterion presence and chain length on the structure and dynamics of single DNA strands of polythymidylate were studied by means of molecular dynamics simulations. The importance of the base-base stacking phenomenon increases with the chain length and partially reduces the flexibility of the strand. Sodium ions directly interact with the phosphate groups and keto oxygens of the thymine bases, complexes showing lifetimes below 400 ps. Simultaneous phosphate and keto complexes were observed for one of the sodium ions with lifetimes around 1 ns. The implications of such complexes in the folding process experienced by the strand are considered. Structurally, cation inner- and outer-sphere complexes were observed in the coordination of phosphate groups. For the inner-sphere complexes, the structural information retrieved from the simulations is in very good agreement with experimental data. The diffusion properties of the sodium ions also reflect both types of coordination modes.