Quantitative analysis of the ion-dependent folding stability of DNA triplexes

A DNA triplex is formed through binding of a third strand to the major groove of a duplex. Due to the high charge density of a DNA triplex, metal ions are critical for its stability. We recently developed the tightly bound ion (TBI) model for ion-nucleic acids interactions. The model accounts for the potential correlation and fluctuations of the ion distribution. We now apply the TBI model to analyze the ion dependence of the thermodynamic stability for DNA triplexes. We focus on two experimentally studied systems: a 24-base DNA triplex and a pair of interacting 14-base triplexes. Our theoretical calculations for the number of bound ions indicate that the TBI model provides improved predictions for the number of bound ions than the classical Poisson-Boltzmann (PB) equation. The improvement is more significant for a triplex, which has a higher charge density than a duplex. This is possibly due to the higher ion concentration around the triplex and hence a stronger ion correlation effect for a triplex. In addition, our analysis for the free energy landscape for a pair of 14-mer triplexes immersed in an ionic solution shows that divalent ions could induce an attractive force between the triplexes. Furthermore, we investigate how the protonated cytosines in the triplexes affect the stability of the triplex helices.     

Predicting RNA-Metal Ion Binding with Ion Dehydration Effects

Metal ions play essential roles in nucleic acids folding and stability. The interaction between metal ions and nucleic acids can be highly complicated because of the interplay between various effects such as ion correlation, fluctuation, and dehydration. These effects may be particularly important for multivalent ions such as Mg²⁺ ions. Previous efforts to model ion correlation and fluctuation effects led to the development of the Monte Carlo tightly bound ion model. Here, by incorporating ion hydration/dehydration effects into the Monte Carlo tightly bound ion model, we develop a, to our knowledge, new approach to predict ion binding. The new model enables predictions for not only the number of bound ions but also the three-dimensional spatial distribution of the bound ions. Furthermore, the new model reveals several intriguing features for the bound ions such as the mutual enhancement/inhibition in ion binding between the fully hydrated (diffuse) ions, the outer-shell dehydrated ions, and the inner-shell dehydrated ions and novel features for the monovalent-divalent ion interplay due to the hydration effect.     

Zinc coordination geometry and ligand binding affinity: the structural and kinetic analysis of the second-shell serine 228 residue and the methionine 180 residue of the aminopeptidase from Vibrio proteolyticus

The chemical properties of zinc make it an ideal metal to study the role of coordination strain in enzymatic rate enhancement. The zinc ion and the protein residues that are bound directly to the zinc ion represent a functional charge/dipole complex, and polarization of this complex, which translates to coordination distortion, may tune electrophilicity, and hence, reactivity. Conserved protein residues outside of the charge/dipole complex, such as second-shell residues, may play a role in supporting the electronic strain produced as a consequence of functional polarization. To test the correlation between charge/dipole polarity and ligand binding affinity, structure-function studies were carried out on the dizinc aminopeptidase from Vibrio proteolyticus. Alanine substitutions of S228 and M180 resulted in catalytically diminished enzymes whose crystal structures show very little change in the positions of the metal ions and the protein residues. However, more detailed inspections of the crystal structures show small positional changes that account for differences in the zinc ion coordination geometry. Measurements of the binding affinity of leucine phosphonic acid, a transition state analogue, and leucine, a product, show a correlation between coordination geometry and ligand binding affinity. These results suggest that the coordination number and polarity may tune the electrophilicity of zinc. This may have provided the evolving enzyme with the ability to discriminate between reaction coordinate species.     

Comparison of polyelectrolyte theories of the binding of cations to DNA.

Predictions of the binding of counterions to DNA made using the counterion condensation theory developed by Manning are compared with those made using the Poisson-Boltzmann equation, solved numerically by the Runge-Kutta procedure. Ions are defined as territorially or atmospherically bound if they fall within a given distance, defined by counterion condensation theory, from the DNA surface. Two types of experimental situations are considered. The first is the delocalized binding of a single type of counterion to DNA. In this case the Poisson-Boltzmann treatment predicts somewhat lower extents of binding TO DNA, modeled as a 10-A radius cylinder, than does Manning theory. The two theories converge as the radius decreases. The second type of experiment is the competition of ions of different valence for binding to DNA. The theories are compared with literature values of binding constants of divalent ions in the presence of monovalent ions, and of spermidine 3+ in the presence of Na+ or Mg2+. Both predict with fair accuracy the salt dependence of the equilibrium constants.     

Thermodynamic properties of the specific binding between Ag+ ions and C:C mismatched base pairs in duplex DNA

Metal-mediated base pairs formed by the interaction between metal ions and artificial bases in oligonucleotides have been developed for potential applications in nanotechnology. We recently found that a natural C:C mismatched base pair bound to an Ag(+) ion to generate a novel metal-mediated base pair in duplex DNA. Preparation of the novel C-Ag-C base pair involving natural bases is more convenient than that of metal-mediated base pairs involving artificial bases because time-consuming base synthesis is not required. Here, we examined the thermodynamic properties of the binding between the Ag(+) ion and each of single and double C:C mismatched base pair in duplex DNA by isothermal titration calorimetry. The Ag(+) ion specifically bound to the C:C mismatched base pair at a 1:1 molar ratio with 10(6) M(-1) binding constant, which was significantly larger than those for nonspecific metal ion-DNA interactions. The specific binding between the Ag(+) ion and the single C:C mismatched base pair was mainly driven by the positive dehydration entropy change and the negative binding enthalpy change. In the interaction between the Ag(+) ion and each of the consecutive and interrupted double C:C mismatched base pairs, stoichiometric binding at a 1:1 molar ratio was achieved in each step of the first and second Ag(+) binding. The binding affinity for the second Ag(+) binding was similar to that for the first Ag(+) binding. Stoichiometric binding without interference and negative cooperativity may be favorable for aligning multiple Ag(+) ions in duplex DNA for applications of the metal-mediated base pairs in nanotechnology.     

Poisson-Boltzmann calculations of nonspecific salt effects on protein-protein binding free energies

The salt dependence of the binding free energy of five protein-protein hetero-dimers and two homo-dimers/tetramers was calculated from numerical solutions to the Poisson-Boltzmann equation. Overall, the agreement with experimental values is very good. In all cases except one involving the highly charged lactoglobulin homo-dimer, increasing the salt concentration is found both experimentally and theoretically to decrease the binding affinity. To clarify the source of salt effects, the salt-dependent free energy of binding is partitioned into screening terms and to self-energy terms that involve the interaction of the charge distribution of a monomer with its own ion atmosphere. In six of the seven complexes studied, screening makes the largest contribution but self-energy effects can also be significant. The calculated salt effects are found to be insensitive to force-field parameters and to the internal dielectric constant assigned to the monomers. Nonlinearities due to high charge densities, which are extremely important in the binding of proteins to negatively charged membrane surfaces and to nucleic acids, make much smaller contributions to the protein-protein complexes studied here, with the exception of highly charged lactoglobulin dimers. Our results indicate that the Poisson-Boltzmann equation captures much of the physical basis of the nonspecific salt dependence of protein-protein complexation.     

Density functional theory for the nonspecific binding of salt to polyelectrolytes: thermodynamic properties

The thermodynamics of the nonspecific binding of salt to a polyelectrolyte molecule is studied using a density functional approach. The polyelectrolyte molecule is modeled as an infinite, inflexible, and impenetrable charged cylinder and the counterions and co-ions are modeled as charged hard spheres of equal diameter. The density functional theory is based on a hybrid approach where the hard-sphere contribution to the one-particle correlation function is evaluated nonperturbatively and the ionic contribution to the one-particle correlation function is evaluated perturbatively. The advantage of the approach is that analytical expressions are available for all the correlation functions. The calculated single ion preferential interaction coefficients, excess free energy, and activity coefficients show a nonmonotonic variation as a function of polyion charge in the presence of divalent ions. These properties display considerable departure from the predictions of the nonlinear Poisson-Boltzmann (NLPB) equation, with qualitative differences in some cases, which may be attributed to correlation effects neglected in the NLPB theory.     

Structural and energetic analysis of metal ions essential to SRP signal recognition domain assembly

The signal recognition particle (SRP) targets proteins to the endoplasmic reticulum in eukaryotes or to the inner membrane in prokaryotes by binding to hydrophobic signal sequences. Signal peptide recognition occurs within the highly conserved RNA-protein core of the SRP, underscoring the importance of this complex in SRP function. Structural analysis of the RNA and protein components of the prokaryotic SRP in the free and bound states revealed that the RNA undergoes a significant conformational change upon protein binding involving the uptake of several monovalent and divalent cations. To investigate the role of these metal ions in formation of the functional SRP complex, we used binding affinity assays and X-ray crystallography to analyze the specificity and energetic contributions of mono- and divalent metal ions bound in the RNA. Our results demonstrate that several metal ion binding sites important for RNA conformation can accommodate chemically distinct ions, often without affecting the structure of the complex. Thus, while these metal ions are highly ordered and essential for the formation and stability of the SRP complex, they behave like nonspecific metal ions.     

RNA helix stability in mixed Na+/Mg2+ solution

A recently developed tightly bound ion model can account for the correlation and fluctuation (i.e., different binding modes) of bound ions. However, the model cannot treat mixed ion solutions, which are physiologically relevant and biologically significant, and the model was based on B-DNA helices and thus cannot directly treat RNA helices. In the present study, we investigate the effects of ion correlation and fluctuation on the thermodynamic stability of finite length RNA helices immersed in a mixed solution of monovalent and divalent ions. Experimental comparisons demonstrate that the model gives improved predictions over the Poisson-Boltzmann theory, which has been found to underestimate the roles of multivalent ions such as Mg2+ in stabilizing DNA and RNA helices. The tightly bound ion model makes quantitative predictions on how the Na+-Mg2+ competition determines helix stability and its helix length-dependence. In addition, the model gives empirical formulas for the thermodynamic parameters as functions of Na+/Mg2+ concentrations and helix length. Such formulas can be quite useful for practical applications.     

Computational approach to site-directed ligand discovery

A computational approach, Systematic Conformational Search & Induced Fit (SCI&FI), to site-directed ligand discovery (Tethering) is presented. SCI&FI has the ability to predict the binding site, binding mode, and bound dynamics of small molecule fragments covalently tethered to a protein. The SCI&FI method was engineered with the ability to model induced fit conformational changes of the protein because of the binding of the tether. SCI&FI generates comprehensive picture of the binding preferences of the tether to the protein by elucidating potential binding sites of the tether and by describing regions of receptor space capable of conformational change because of the binding of the tether. The SCI&FI method provides a complementary approach to experimental tethering. Initial validation of the SCI&FI method is reported by predicting the 3D structure of two Interleukin-2 and an Interleukin-4 tethered-protein systems.     

Structure model of a complex between the factor for inversion stimulation (FIS) and DNA: Modeling protein-DNA complexes with dyad symmetry and known protein structures

A method is presented to predict overall conformations of protein-DNA complexes on the basis of the known three-dimensional structures of the proteins. The method is restricted to proteins with a common twofold symmetry axis, which show only minor conformational changes upon binding to DNA. The method uses a numerical finite difference solution of the linearized Poisson-Boltzmann equation and subsequent energy minimization cycles. Structural parameters—the rotation angle of the DNA relative to the protein around the common symmetry axis, the protein-DNA distance, and intermolecular hydrogen-bonding contacts—are presented for two test cases, DNA bound to CAP (catabolite gene activator protein) and to the Cro-repressor of bacteriophage 434. The DNA curvature in the starting model of the docking procedure was chosen as a smoothed approximation of the conformation found in the X-ray structures of these complexes. The method is further used to predict the unknown structure of the complex between the factor for inversion stimulation (FIS) and DNA, which is bent upon binding to FIS. In contrast to the test cases, the unknown curvature of the starting model is derived from a calibration of electrostatic precalculations for different proteins according to crystallographically observed DNA bending. The results of the modeling are in good accordance with the experimentally observed overall structure of protein-DNA complexes for the two test cases; for FIS, they correspond to several of the experimentally proposed protein-DNA contacts. © 1996 Wiley-Liss, Inc.     

Metal binding motif in the active site of the HDV ribozyme binds divalent and monovalent ions

The hepatitis delta virus (HDV) ribozyme uses both metal ion and nucleobase catalysis in its cleavage mechanism. A reverse G·U wobble was observed in a recent crystal structure of the precleaved state. This unusual base pair positions a Mg(2+) ion to participate in catalysis. Herein, we used molecular dynamics (MD) and X-ray crystallography to characterize the conformation and metal binding characteristics of this base pair in product and precleaved forms. Beginning with a crystal structure of the product form, we observed formation of the reverse G·U wobble during MD trajectories. We also demonstrated that this base pair is compatible with the diffraction data for the product-bound state. During MD trajectories of the product form, Na(+) ions interacted with the reverse G·U wobble in the RNA active site, and a Mg(2+) ion, introduced in certain trajectories, remained bound at this site. Beginning with a crystal structure of the precleaved form, the reverse G·U wobble with bound Mg(2+) remained intact during MD simulations. When we removed Mg(2+) from the starting precleaved structure, Na(+) ions interacted with the reverse G·U wobble. In support of the computational results, we observed competition between Na(+) and Mg(2+) in the precleaved ribozyme crystallographically. Nonlinear Poisson-Boltzmann calculations revealed a negatively charged patch near the reverse G·U wobble. This anionic pocket likely serves to bind metal ions and to help shift the pK(a) of the catalytic nucleobase, C75. Thus, the reverse G·U wobble motif serves to organize two catalytic elements, a metal ion and catalytic nucleobase, within the active site of the HDV ribozyme.     

Prediction of structures of zinc‐binding proteins through explicit modeling of metal coordination geometry

Metal ions play an essential role in stabilizing protein structures and contributing to protein function. Ions such as zinc have well‐defined coordination geometries, but it has not been easy to take advantage of this knowledge in protein structure prediction efforts. Here, we present a computational method to predict structures of zinc‐binding proteins given knowledge of the positions of zinc‐coordinating residues in the amino acid sequence. The method takes advantage of the “atom‐tree” representation of molecular systems and modular architecture of the Rosetta3 software suite to incorporate explicit metal ion coordination geometry into previously developed de novo prediction and loop modeling protocols. Zinc cofactors are tethered to their interacting residues based on coordination geometries observed in natural zinc‐binding proteins. The incorporation of explicit zinc atoms and their coordination geometry in both de novo structure prediction and loop modeling significantly improves sampling near the native conformation. The method can be readily extended to predict protein structures bound to other metal and/or small chemical cofactors with well‐defined coordination or ligation geometry.     

Lanthanide binding and IgG affinity construct: Potential applications in solution NMR, MRI, and luminescence microscopy

Paramagnetic lanthanide ions when bound to proteins offer great potential for structural investigations that utilize solution nuclear magnetic resonance spectroscopy, magnetic resonance imaging, or optical microscopy. However, many proteins do not have native metal ion binding sites and engineering a chimeric protein to bind an ion while retaining affinity for a protein of interest represents a significant challenge. Here we report the characterization of an immunoglobulin G-binding protein redesigned to include a lanthanide binding motif in place of a loop between two helices (Z-L2LBT). It was shown to bind Tb(3+) with 130 nM affinity. Ions such as Dy(3+) , Yb(3+) , and Ce(3+) produce paramagnetic effects on NMR spectra and the utility of these effects is illustrated by their use in determining a structural model of the metal-complexed Z-L2LBT protein and a preliminary characterization of the dynamic distribution of IgG Fc glycan positions. Furthermore, this designed protein is demonstrated to be a novel IgG-binding reagent for magnetic resonance imaging (Z-L2LBT:Gd(3+) complex) and luminescence microscopy (Z-L2LBT: Tb(3+) complex).     

Identification of the Third Na Site and the Sequence of Extracellular Binding Events in the Glutamate Transporter

The transport cycle in the glutamate transporter (GlT) is catalyzed by the cotransport of three Na⁺ ions. However, the positions of only two of these ions (Na1 and Na2 sites) along with the substrate have been captured in the crystal structures reported for both the outward-facing and the inward-facing states of Glt_ph. Characterizing the third ion binding site (Na3) is necessary for structure-function studies attempting to investigate the mechanism of transport in GlTs at an atomic level, particularly for the determination of the sequence of the binding events during the transport cycle. In this study, we report a series of molecular dynamics simulations performed on various bound states of Glt_ph (the apo state, as well as in the presence of Na⁺, the substrate, or both), which have been used to identify a putative Na3 site. The calculated trajectories have been used to determine the water accessibility of potential ion-binding residues in the protein, as a prerequisite for their ion binding. Combined with conformational analysis of the key regions in the protein in different bound states and several additional independent simulations in which a Na⁺ ion was randomly introduced to the interior of the transporter, we have been able to characterize a putative Na3 site and propose a plausible binding sequence for the substrate and the three Na⁺ ions to the transporter during the extracellular half of the transport cycle. The proposed Na3 site is formed by a set of highly conserved residues, namely, Asp³¹², Thr⁹², and Asn³¹⁰, along with a water molecule. Simulation of a fully bound state, including the substrate and the three Na⁺ ions, reveals a stable structure—showing closer agreement to the crystal structure when compared to previous models lacking an ion in the putative Na3 site. The proposed sequence of binding events is in agreement with recent experimental models suggesting that two Na⁺ ions bind before the substrate, and one after that. Our results, however, provide additional information about the sites involved in these binding events.     

Tests of continuum theories as models of ion channels. I. Poisson-Boltzmann theory versus Brownian dynamics

Continuum theories of electrolytes are widely used to describe physical processes in various biological systems. Although these are well-established theories in macroscopic situations, it is not clear from the outset that they should work in small systems whose dimensions are comparable to or smaller than the Debye length. Here, we test the validity of the mean-field approximation in Poisson-Boltzmann theory by comparing its predictions with those of Brownian dynamics simulations. For this purpose we use spherical and cylindrical boundaries and a catenary shape similar to that of the acetylcholine receptor channel. The interior region filled with electrolyte is assumed to have a high dielectric constant, and the exterior region representing protein a low one. Comparisons of the force on a test ion obtained with the two methods show that the shielding effect due to counterions is overestimated in Poisson-Boltzmann theory when the ion is within a Debye length of the boundary. As the ion gets closer to the boundary, the discrepancy in force grows rapidly. The implication for membrane channels, whose radii are typically smaller than the Debye length, is that Poisson-Boltzmann theory cannot be used to obtain reliable estimates of the electrostatic potential energy and force on an ion in the channel environment.     

Nucleic acid helix stability: effects of salt concentration, cation valence and size, and chain length

Metal ions play crucial roles in thermal stability and folding kinetics of nucleic acids. For ions (especially multivalent ions) in the close vicinity of nucleic acid surface, interion correlations and ion-binding mode fluctuations may be important. Poisson-Boltzmann theory ignores these effects whereas the recently developed tightly bound ion (TBI) theory explicitly accounts for these effects. Extensive experimental data demonstrate that the TBI theory gives improved predictions for multivalent ions (e.g., Mg2+) than the Poisson-Boltzmann theory. In this study, we use the TBI theory to investigate how the metal ions affect the folding stability of B-DNA helices. We quantitatively evaluate the effects of ion concentration, ion size and valence, and helix length on the helix stability. Moreover, we derive practically useful analytical formulas for the thermodynamic parameters as functions of finite helix length, ion type, and ion concentration. We find that the helix stability is additive for high ion concentration and long helix and nonadditive for low ion concentration and short helix. All these results are tested against and supported by extensive experimental data.