Spermine: an "invisible" component in the crystals of B-DNA. A grand canonical Monte Carlo and molecular dynamics simulation study

The association of spermine(4+) (Spm(4+)), Mg(2+) and monovalent (M(+)) ions with DNA in crystal form, have been studied using grand canonical Monte Carlo (GCMC) and molecular dynamics (MD) computer simulations. GCMC calculations were used to calculate the distribution of Spm(4+), Mg(2+), and M(+) between the equilibrating solvent and the DNA crystal under conditions mimicking the crystal-growing protocols reported in a number of recent X-ray diffraction studies of DNA oligomers. The GCMC simulations show that the composition of ions neutralizing the negative charge of DNA can vary in a broad range. The GCMC simulations were used to provide appropriate conditions for subsequent 6 ns constant pressure and temperature MD simulations of DNA in a typical crystalline environment consisting of three DNA double helix decamers in a periodic hexagonal cell, containing 1200 water molecules, eight Spm(4+), 32 Na(+) and four Cl(-) ions. Based on the simulation results, it seems possible to give an explanation why spermine molecules are usually not detected in X-ray studies in spite of their high concentration in the preparatory samples used as the crystallizing agent. It appears that this flexible polyamine molecule has several binding modes, interacting in fairly irregular manner with different sites on DNA and showing no regular ordering in the DNA crystals. Ions of Na(+) and Spm(4+) compete with each other and with water molecules in binding to bases in the minor groove and they influence the structure of the DNA hydration shell in different ways.     

Solution probing of metal ion binding by helix 27 from Escherichia coli 16S rRNA

Helix (H)27 from Escherichia coli 16S ribosomal (r)RNA is centrally located within the small (30S) ribosomal subunit, immediately adjacent to the decoding center. Bacterial 30S subunit crystal structures depicting Mg(2+) binding sites resolve two magnesium ions within the vicinity of H27: one in the major groove of the G886-U911 wobble pair, and one within the GCAA tetraloop. Binding of such metal cations is generally thought to be crucial for RNA folding and function. To ask how metal ion-RNA interactions in crystals compare with those in solution, we have characterized, using solution NMR spectroscopy, Tb(3+) footprinting and time-resolved fluorescence resonance energy transfer (tr-FRET), location, and modes of metal ion binding in an isolated H27. NMR and Tb(3+) footprinting data indicate that solution secondary structure and Mg(2+) binding are generally consistent with the ribosomal crystal structures. However, our analyses also suggest that H27 is dynamic in solution and that metal ions localize within the narrow major groove formed by the juxtaposition of the loop E motif with the tandem G894-U905 and G895-U904 wobble pairs. In addition, tr-FRET studies provide evidence that Mg(2+) uptake by the H27 construct results in a global lengthening of the helix. We propose that only a subset of H27-metal ion interactions has been captured in the crystal structures of the 30S ribosomal subunit, and that small-scale structural dynamics afforded by solution conditions may contribute to these differences. Our studies thus highlight an example for differences between RNA-metal ion interactions observed in solution and in crystals.     

Interactions of cations with the cytoplasmic pores of inward rectifier K(+) channels in the closed state

Ion channels gate at membrane-embedded domains by changing their conformation along the ion conduction pathway. Inward rectifier K(+) (Kir) channels possess a unique extramembrane cytoplasmic domain that extends this pathway. However, the relevance and contribution of this domain to ion permeation remain unclear. By qualitative x-ray crystallographic analysis, we found that the pore in the cytoplasmic domain of Kir3.2 binds cations in a valency-dependent manner and does not allow the displacement of Mg(2+) by monovalent cations or spermine. Electrophysiological analyses revealed that the cytoplasmic pore of Kir3.2 selectively binds positively charged molecules and has a higher affinity for Mg(2+) when it has a low probability of being open. The selective blocking of chemical modification of the side chain of pore-facing residues by Mg(2+) indicates that the mode of binding of Mg(2+) is likely to be similar to that observed in the crystal structure. These results indicate that the Kir3.2 crystal structure has a closed conformation with a negative electrostatic field potential at the cytoplasmic pore, the potential of which may be controlled by conformational changes in the cytoplasmic domain to regulate ion diffusion along the pore.     

X-ray structural analysis of magnesium-containing Serratia marcescens endonuclease]

The three-dimensional crystal structure of the DNA/RNA nonspecific endonuclease from Serratia marcescens was refined at the resolution of 1.07 A to R factor of 12.4% and Rfree factor of 15.3% using the anisotropic approximation. The structure includes 3924 non-hydrogen atoms, 715 protein-bound water molecules, and a Mg2+ ion in each binding site of each subunit of the nuclease homodimeric globular molecule. The 3D topological model of the enzyme was revealed, the inner symmetry of the monomers in its N- and C-termini was found, and the local environment of the magnesium cofactor in the nuclease active site was defined. Mg2+ ion was found to be bound to the Asn119 residue and surrounded by five associated water molecules that form an octahedral configuration. The coordination distances for the water molecules and the O delta 1 atom of Asn119 were shown to be within a range of 2.01-2.11 A. The thermal factors for the magnesium ion in subunits are 7.08 and 4.60 A2, and the average thermal factors for the surrounding water molecules are 11.14 and 10.30 A2, respectively. The region of the nuclease subunit interactions was localized, and the alternative side chain conformations were defined for 51 amino acid residues of the nuclease dimer.     

The crystal structure of yeast phenylalanine tRNA at 2.0 A resolution: cleavage by Mg(2+) in 15-year old crystals

We have re-determined the crystal structure of yeast tRNA(Phe) to 2. 0 A resolution using 15 year old crystals. The accuracy of the new structure, due both to higher resolution data and formerly unavailable refinement methods, consolidates the previous structural information, but also reveals novel details. In particular, the water structure around the tightly bound Mg(2+) is now clearly resolved, and hence provides more accurate information on the geometry of the magnesium-binding sites and the role of water molecules in coordinating the metal ions to the tRNA. We have assigned a total of ten magnesium ions and identified a partly conserved geometry for high-affinity Mg(2+ )binding. In the electron density map there is also clear density for a spermine molecule binding in the major groove of the TPsiC arm and also contacting a symmetry-related tRNA molecule. Interestingly, we have also found that two specific regions of the tRNA in the crystals are partially cleaved. The sites of hydrolysis are within the D and anticodon loops in the vicinity of Mg(2+).     

An orthorhombic crystal form of cyclohexaicosaose, CA26.32.59 H(2)O: comparison with the triclinic form.

Cycloamylose containing 26 glucose residues (cyclohexaicosaose, CA26) crystallized from water and 30% (v/v) polyethyleneglycol 400 in the orthorhombic space group P2(1)2(1)2(1) in the highly hydrated form CA26.32.59 H(2)O. X-ray analysis of the crystals at 0.85 A resolution shows that the macrocycle of CA26 is folded into two short left-handed V-amylose helices in antiparallel arrangement and related by a twofold rotational pseudosymmetry as reported recently for the (CA26)(2).76.75 H(2)O triclinic crystal form [Gessler, K. et al. Proc. Natl. Acad. Sci. USA 1999, 96, 4246-4251]. In the orthorhombic crystal form, CA26 molecules are packed in motifs reminiscent of V-amylose in hydrated and anhydrous forms. The intramolecular interface between the V-helices in CA26 is dictated by formation of an extended network of interhelical C-H...O hydrogen bonds; a comparable molecular arrangement is also evident for the intermolecular packing, suggesting that it is a characteristic feature of V-amylose interaction. The hydrophobic channels of CA26 are filled with disordered water molecules arranged in chains and held in position by multiple C-H...O hydrogen bonds. In the orthorhombic and triclinic crystal forms, the structures of CA26 molecules are equivalent but the positions of the individual water molecules are different, suggesting that the patterns of water chains are perturbed even by small structural changes associated with differences in packing arrangements in the two crystal lattices rather than with differences in the CA26 geometry.     

Biochemical and structural studies on the high affinity of Hsp70 for ADP

The molecular chaperone 70-kDa heat shock protein (Hsp70) is driven by ATP hydrolysis and ADP-ATP exchange. ADP dissociation from Hsp70 is reportedly slow in the presence of inorganic phosphate (P(i) ). In this study, we investigated the interaction of Hsp70 and its nucleotide-binding domain (NBD) with ADP in detail, by isothermal titration calorimetry measurements and found that Mg(2+) ion dramatically elevates the affinity of Hsp70 for ADP. On the other hand, P(i) increased the affinity in the presence of Mg(2+) ion, but not in its absence. Thus, P(i) enhances the effect of the Mg(2+) ion on the ADP binding. Next, we determined the crystal structures of the ADP-bound NBD with and without Mg(2+) ion. As compared with the Mg(2+) ion-free structure, the ADP- and Mg(2+) ion-bound NBD contains one Mg(2+) ion, which is coordinated with the β-phosphate group of ADP and associates with Asp10, Glu175, and Asp199, through four water molecules. The Mg(2+) ion is also coordinated with one P(i) molecule, which interacts with Lys71, Glu175, and Thr204. In fact, the mutations of Asp10 and Asp199 reduced the affinity of the NBD for ADP, in both the presence and the absence of P(i) . Therefore, the Mg(2+) ion-mediated network, including the P(i) and water molecules, increases the affinity of Hsp70 for ADP, and thus the dissociation of ADP is slow. In ADP-ATP exchange, the slow ADP dissociation might be rate-limiting. However, the nucleotide-exchange factors actually enhance ADP release by disrupting the Mg(2+) ion-mediated network.     

A minimized human integrin alpha(5)beta(1) that retains ligand recognition

Two isolated recombinant fragments from human integrin alpha(5)beta(1) encompassing the FG-GAP repeats III to VII of alpha(5) and the insertion-type domain from beta(1), respectively, are structurally well defined in solution, based on CD evidence. Divalent cation binding induces a conformational adaptation that is achieved by Ca(2+) or Mg(2+) (or Mn(2+)) with alpha(5) and only by Mg(2+) (or Mn(2+)) with beta(1). Mn(2+) bound to beta(1) is highly hydrated ( approximately 3 water molecules), based on water NMR relaxation, in agreement with a metal ion-dependent adhesion site-type metal coordination. Each fragment saturated with Mg(2+) (or Mn(2+)) binds a recombinant fibronectin ligand in an RGD-dependent manner. A conformational rearrangement is induced on the fibronectin ligand upon binding to the alpha(5), but not to the beta(1) fragment, based on CD. Ligand binding results in metal ion displacement from beta(1). Both alpha(5) and beta(1) fragments form a stable heterodimer (alpha(5)beta(1) mini-integrin) that retains ligand recognition to form a 1:1:1 ternary complex, in the presence of Mg(2+), and induces a specific conformational adaptation of the fibronectin ligand. A two-site model for RGD binding to both alpha and beta integrin components is inferred from our data using low molecular weight RGD mimetics.     

The crystal structure of R-specific alcohol dehydrogenase from Lactobacillus brevis suggests the structural basis of its metal dependency

The crystal structure of the apo-form of an R-specific alcohol dehydrogenase from Lactobacillus brevis (LB-RADH) was solved and refined to 1.8A resolution. LB-RADH is a member of the short-chain dehydrogenase/reductase (SDR) enyzme superfamily. It is a homotetramer with 251 amino acid residues per subunit and uses NADP(H) as co-enzyme. NADPH and the substrate acetophenone were modelled into the active site. The enantiospecificity of the enzyme can be explained on the basis of the resulting hypothetical ternary complex. In contrast to most other SDR enzymes, the catalytic activity of LB-RADH depends strongly on the binding of Mg(2+). Mg(2+) removal by EDTA inactivates the enzyme completely. In the crystal structure, the Mg(2+)-binding site is well defined. The ion has a perfect octahedral coordination sphere and occupies a special position concerning crystallographic and molecular point symmetry, meaning that each RADH tetramer contains two magnesium ions. The magnesium ion is no direct catalytic cofactor. However, it is structurally coupled to the putative C-terminal hinge of the substrate-binding loop and, via an extended hydrogen bonding network, to some side-chains forming the substrate binding region. Therefore, the presented structure of apo-RADH provides plausible explanations for the metal dependence of the enzyme.     

The structure of sitting-atop complexes of metalloporphyrins studied by theoretical methods

The metallation of tetrapyrroles is believed to proceed via a sitting-atop (SAT) complex, in which some of the pyrrole nitrogen atoms are still protonated and the metal ion resides above the ring plane. No crystal structure of such a complex has been presented, but NMR and extended X-ray absorption fine structure (EXAFS) data has been reported for Cu(2+) in acetonitrile. We have used density functional calculations to obtain reasonable models for SAT complexes of porphyrins with Mg(2+), Fe(2+), and Cu(2+). The results show that there are many possible SAT complexes with 1-5 solvent molecules, one or two metal ions, and cis or trans protonation of the porphyrin ring. Many of these have similar energies and their relative stabilities vary with the metal ion. A complex with two cis pyrrolenine nitrogens atoms and 2-4 solvent molecules coordinated to Cu(2+) fits the NMR and EXAFS data best. However, we cannot fully exclude the possibility that what is observed is rather a mixture of a doubly protonated porphyrin and the copper porphyrin. Mg(2+) has a lower affinity for porphyrin and stronger affinity for water, so a complex with five water molecules and only one bond to porphyrin seems to be most stable. For Fe(2+), a cis structure with two first-sphere water molecules and four interactions to the porphyrin seems to be most likely.     

Magnesium ion-dependent activation of the RecA protein involves the C terminus

Optimal conditions for RecA protein-mediated DNA strand exchange include 6-8 mm Mg(2+) in excess of that required to form complexes with ATP. We provide evidence that the free magnesium ion is required to mediate a conformational change in the RecA protein C terminus that activates RecA-mediated DNA strand exchange. In particular, a "closed" (low Mg(2+)) conformation of a RecA nucleoprotein filament restricts DNA pairing by incoming duplex DNA, although single-stranded overhangs at the ends of a duplex allow limited DNA pairing to occur. The addition of excess Mg(2+) results in an "open" conformation, which can promote efficient DNA pairing and strand exchange regardless of DNA end structure. The removal of 17 amino acid residues at the Escherichia coli RecA C terminus eliminates a measurable requirement for excess Mg(2+) and permits efficient DNA pairing and exchange similar to that seen with the wild-type protein at high Mg(2+) levels. Thus, the RecA C terminus imposes the need for the high magnesium ion concentrations requisite in RecA reactions in vitro. We propose that the C terminus acts as a regulatory switch, modulating the access of double-stranded DNA to the presynaptic filament and thereby inhibiting homologous DNA pairing and strand exchange at low magnesium ion concentrations.     

Mg2+ cation initiates the conversion of nucleoside triphosphate to nucleoside monophospate through a radical mechanism. Quantum molecular dynamics simulation

DFT:B3LYP (6-311G** basis set) quantum molecular dynamics simulation was used to study the conversion of guanosine triphosphate (GTP) to guanosine monophosphate upon the action of Mg2+. The computations were carried out at 310 K in a basin of 178 water molecules surrounding a Mg(2+)-GTP complex and imitating the behavior of the solvent. The cleavage of the Mg(2+)-GTP complex occurs over the 5-ps period and gives rise to two inorganic phosphates (Pi), a hydrated Mg2+ complex, atomic oxygen, and highly reactive GMP radical. The appearance of this radical is a result of action of the Mg2+ cation, which initiates the radical mechanism of GTP cleavage. At the very early step of interaction with GTP, Mg2+ is reduced to Mg+, thus producing an ion radical pair (+)Mg-GTP(3-). In the absence of Mg2+, a non-reactive form of GMP is formed rather than GMP; the process corresponds to hydrolytic cleavage of GTP through the ionic mechanism. The formation of GMP and its analogues with adenosine, cytidine, thymidine, and uridine is, seemingly, a key point in DNA and RNA synthesis.     

Atomic-resolution crystal structures of B-DNA reveal specific influences of divalent metal ions on conformation and packing.

Crystal structures of B-form DNA have provided insights into the global and local conformational properties of the double helix, the solvent environment, drug binding and DNA packing. For example, structures of the duplex with sequence CGCGAATTCGCG, the Dickerson-Drew dodecamer (DDD), established a unique geometry of the central A-tract and a hydration spine in the minor groove. However, our knowledge of the various interaction modes between metal ions and DNA is very limited and almost no information exists concerning the origins of the different effects on DNA conformation and packing exerted by individual metal ions.Crystallization of the DDD duplex in the presence of Mg(2+)and Ca(2+)yields different crystal forms. The structures of the new Ca(2+)-form and isomorphous structures of oligonucleotides with sequences GGCGAATTCGCG and GCGAATTCGCG were determined at a maximum resolution of 1.3 A. These and the 1.1 A structure of the DDD Mg(2+)-form have revealed the most detailed picture yet of the ionic environment of B-DNA. In the Mg(2+)and Ca(2+)-forms, duplexes in the crystal lattice are surrounded by 13 magnesium and 11 calcium ions, respectively.Mg(2+)and Ca(2+)generate different DNA crystal lattices and stabilize different end-to-end overlaps and lateral contacts between duplexes, thus using different strategies for reducing the effective repeat length of the helix to ten base-pairs. Mg(2+)crystals allow the two outermost base-pairs at either end to interact laterally via minor groove H-bonds, turning the 12-mer into an effective 10-mer. Ca(2+)crystals, in contrast, unpair the outermost base-pair at each end, converting the helix into a 10-mer that can stack along its axis. This reduction of a 12-mer into a functional 10-mer is followed no matter what the detailed nature of the 5'-end of the chain: C-G-C-G-A-ellipsis, G-G-C-G-A-ellipsis, or a truncated G-C-G-A-ellipsis Rather than merely mediating close contacts between phosphate groups, ions are at the origin of many well-known features of the DDD duplex structure. A Mg(2+)coordinates in the major groove, contributing to kinking of the duplex at one end. While Ca(2+)resides in the minor groove, coordinating to bases via its hydration shell, two magnesium ions are located at the periphery of the minor groove, bridging phosphate groups from opposite strands and contracting the groove at one border of the A-tract.     

Magnesium-sugar interaction. Synthesis, spectroscopic and structural characterization of Mg-sugar complexes containing beta-D-fructose

The interaction between beta-D-fructose and hydrated magnesium salts has been studied and complexes of the type Mg(beta-D-fructose)Cl2.4H2O and Mg(beta-D-fructose)Br2.4H2O have been isolated and characterized. On the basis of comparisons of the spectroscopic and other chemical properties of several structurally known calcium-fructose compounds with those of the corresponding magnesium complexes, it is concluded that Mg2+ binds to two sugar moieties via O(2), O(3) of the first and O(4), O(5) of the second and to two water molecules, resulting in a six-coordinate geometry around the Mg2+. The strong sugar hydrogen-bonding network is rearranged upon sugar metallation and the sugar moiety shows the beta-anomer conformation in these magnesium-sugar complexes.     

Mutagenesis of E477 or K505 in the B' domain of human topoisomerase II beta increases the requirement for magnesium ions during strand passage

A type II topoisomerase is essential for decatenating DNA replication products, and it accomplishes this task by passing one DNA duplex through a transient break in a second duplex. The B' domain of topoisomerase II contains three highly conserved motifs, EGDSA, PL(R/K)GK(I/L/M)LNVR, and IMTD(Q/A)DXD. We have investigated these motifs in topoisomerase II beta by mutagenesis, and report that they play a critical role in establishing the DNA cleavage-religation equilibrium. In addition, the mutations E477Q (EGDSA) and K505E (PLRGKILNVR) increase the optimal magnesium ion concentration for strand passage, without affecting the Mg(2+) dependence of ATP hydrolysis. It is likely that the binding affinity of the magnesium ion(s) specifically required for DNA cleavage has been reduced by these mutations. The crystal structure of yeast topo II indicates that residues E477 and K505 may help to position the three aspartate residues of the IMTD(Q/A)DXD motif for magnesium ion coordination, and we propose two possible locations for the magnesium ion binding site(s). These observations are consistent with a previous model in which the B' domain is positioned such that these acidic residues lie next to the active site tyrosine residue. A magnesium ion bound by these aspartate residues could therefore mediate the DNA cleavage-religation reaction.     

Metal ions and phosphate binding in the H-N-H motif: crystal structures of the nuclease domain of ColE7/Im7 in complex with a phosphate ion and different divalent metal ions

H-N-H is a motif found in the nuclease domain of a subfamily of bacteria toxins, including colicin E7, that are capable of cleaving DNA nonspecifically. This H-N-H motif has also been identified in a subfamily of homing endonucleases, which cleave DNA site specifically. To better understand the role of metal ions in the H-N-H motif during DNA hydrolysis, we crystallized the nuclease domain of colicin E7 (nuclease-ColE7) in complex with its inhibitor Im7 in two different crystal forms, and we resolved the structures of EDTA-treated, Zn(2+)-bound and Mn(2+)-bound complexes in the presence of phosphate ions at resolutions of 2.6 A to 2.0 A. This study offers the first determination of the structure of a metal-free and substrate-free enzyme in the H-N-H family. The H-N-H motif contains two antiparallel beta-strands linked to a C-terminal alpha-helix, with a divalent metal ion located in the center. Here we show that the metal-binding sites in the center of the H-N-H motif, for the EDTA-treated and Mg(2+)-soaked complex crystals, were occupied by water molecules, indicating that an alkaline earth metal ion does not reside in the same position as a transition metal ion in the H-N-H motif. However, a Zn(2+) or Mn(2+) ions were observed in the center of the H-N-H motif in cases of Zn(2+) or Mn(2+)-soaked crystals, as confirmed in anomalous difference maps. A phosphate ion was found to bridge between the divalent transition metal ion and His545. Based on these structures and structural comparisons with other nucleases, we suggest a functional role for the divalent transition metal ion in the H-N-H motif in stabilizing the phosphoanion in the transition state during hydrolysis.     

Distinct reaction pathway promoted by non-divalent-metal cations in a tertiary stabilized hammerhead ribozyme

Divalent ion sensitivity of hammerhead ribozymes is significantly reduced when the RNA structure includes appropriate tertiary stabilization. Therefore, we investigated the activity of the tertiary stabilized "RzB" hammerhead ribozyme in several nondivalent ions. Ribozyme RzB is active in spermidine and Na(+) alone, although the cleavage rates are reduced by more than 1,000-fold relative to the rates observed in Mg(2+) and in transition metal ions. The trivalent cobalt hexammine (CoHex) ion is often used as an exchange-inert analog of hydrated magnesium ion. Trans-cleavage rates exceeded 8 min(-1) in 20 mM CoHex, which promoted cleavage through outersphere interactions. The stimulation of catalysis afforded by the tertiary structural interactions within RzB does not require Mg(2+), unlike other extended hammerhead ribozymes. Site-specific interaction with at least one Mg(2+) ion is suggested by CoHex competition experiments. In the presence of a constant, low concentration of Mg(2+), low concentrations of CoHex decreased the rate by two to three orders of magnitude relative to the rate in Mg(2+) alone. Cleavage rates increased as CoHex concentrations were raised further, but the final fraction cleaved was lower than what was observed in CoHex or Mg(2+) alone. These observations suggest that Mg(2+) and CoHex compete for binding and that they cause misfolded structures when they are together. The results of this study support the existence of an alternate catalytic mechanism used by nondivalent ions (especially CoHex) that is distinct from the one promoted by divalent metal ions, and they imply that divalent metals influence catalysis through a specific nonstructural role.     

X-ray Structures of Magnesium and Manganese Complexes with the N-Terminal Domain of Calmodulin: Insights into the Mechanism and Specificity of Metal Ion Binding to an EF-Hand

Calmodulin (CaM), a member of the EF-hand superfamily, regulates many aspects of cell function by responding specifically to micromolar concentrations of Ca(2+) in the presence of an ～1000-fold higher concentration of cellular Mg(2+). To explain the structural basis of metal ion binding specificity, we have determined the X-ray structures of the N-terminal domain of calmodulin (N-CaM) in complexes with Mg(2+), Mn(2+), and Zn(2+). In contrast to Ca(2+), which induces domain opening in CaM, octahedrally coordinated Mg(2+) and Mn(2+) stabilize the closed-domain, apo-like conformation, while tetrahedrally coordinated Zn(2+) ions bind at the protein surface and do not compete with Ca(2+). The relative positions of bound Mg(2+) and Mn(2+) within the EF-hand loops are similar to those of Ca(2+); however, the Glu side chain at position 12 of the loop, whose bidentate interaction with Ca(2+) is critical for domain opening, does not bind directly to either Mn(2+) or Mg(2+), and the vacant ligand position is occupied by a water molecule. We conclude that this critical interaction is prevented by specific stereochemical constraints imposed on the ligands by the EF-hand β-scaffold. The structures suggest that Mg(2+) contributes to the switching off of calmodulin activity and possibly other EF-hand proteins at the resting levels of Ca(2+). The Mg(2+)-bound N-CaM structure also provides a unique view of a transiently bound hydrated metal ion and suggests a role for the hydration water in the metal-induced conformational change.     

Salt or ion bridges in biological systems: a study employing quantum and molecular mechanics

Equilibrium geometries and binding energies of model "salt" or "ion" bridge systems have been computed by ab initio quantum chemistry techniques (GAUSSIAN82) and by empirical force field techniques (AMBER2.0). Formate and dimethyl phosphate served as anions in the model compounds while interacting with several organic cations, including methyl ammonium, methyl guanidinium, and divalent metal ion (either Mg2+ or Ca2+) without and with an additional chloride; and a divalent metal ion (either Mg2+ or Ca2+), chloride, and four water molecules of hydration about the metal ion. The majority of the quantum chemical computations were performed using a split-valence basis set. For the model compounds studied we find that the ab initio optimized geometries are in remarkably good agreement with the molecular mechanics geometries. Several calculations were also performed using diffuse fractions. The formate anion binds these model cations more strongly than does dimethyl phosphate, while the organic cation methyl ammonium binds model anions more strongly than does methyl guanidinium. Finally, in model compounds including organic anions, Mg2+ or Ca2+ and four molecules of water, and a chloride anion, we find that the equilibrium structure of the magnesium complex involves a solvent separated ion pair (the magnesium ion is six coordinate), whereas the calcium ion complex remains seven coordinate. Molecular mechanics overestimates binding energies, but the estimates may be close enough to actual binding energies to give useful insight into the details of salt bridges in biological systems.