Native-state conformational dynamics of GART: a regulatory pH-dependent coil-helix transition examined by electrostatic calculations.

Glycinamide ribonucleotide transformylase (GART) undergoes a pH-dependent coil-helix transition with pK(a) approximately 7. An alpha-helix is formed at high pH spanning 8 residues of a 21-residue-long loop, comprising the segment Thr120-His121-Arg122-Gln123-Ala124-Leu125-Glu126-Asn127. To understand the electrostatic nature of this loop-helix, called the activation loop-helix, which leads to the formation and stability of the alpha-helix, pK(a) values of all ionizable residues of GART have been calculated, using Poisson-Boltzmann electrostatic calculations and crystallographic data. Crystallographic structures of high and low pH E70A GART have been used in our analysis. Low pK(a) values of 5.3, 5.3, 3.9, 1.7, and 4.7 have been calculated for five functionally important histidines, His108, His119, His121, His132, and His137, respectively, using the high pH E70A GART structure. Ten theoretical single and double mutants of the high pH E70A structure have been constructed to identify pairwise interactions of ionizable residues, which have aided in elucidating the multiplicity of electrostatic interactions of the activation loop-helix, and the impact of the activation helix on the catalytic site. Based on our pK(a) calculations and structural data, we propose that: (1) His121 forms a molecular switch for the coil-helix transition of the activation helix, depending on its protonation state; (2) a strong electrostatic interaction between His132 and His121 is observed, which can be of stabilizing or destabilizing nature for the activation helix, depending on the relative orientation and protonation states of the rings of His121 and His132; (3) electrostatic interactions involving His119 and Arg122 play a role in the stability of the activation helix; and (4) the activation helix contains the helix-promoting sequence Arg122-Gln123-Ala124-Leu125-Glu126, but its alignment relative to the N and C termini of the helix is not optimal, and is possibly of a destabilizing nature. Finally, we provide electrostatic evidence that the formation and closure of the activation helix create a hydrophobic environment for catalytic-site residue His108, to facilitate catalysis.

Proton transfer dynamics of GART: the pH-dependent catalytic mechanism examined by electrostatic calculations.

The enzyme glycinamide ribonucleotide transformylase (GART) catalyzes the transfer of a formyl group from formyl tetrahydrofolate (fTHF) to glycinamide ribonucleotide (GAR), a process that is pH-dependent with pK(a) of approximately 8. Experimental studies of pH-rate profiles of wild-type and site-directed mutants of GART have led to the proposal that His108, Asp144, and GAR are involved in catalysis, with His108 being an acid catalyst, while forming a salt bridge with Asp144, and GAR being a nucleophile to attack the formyl group of fTHF. This model implied a protonated histidine with pK(a) of 9.7 and a neutral GAR with pK(a) of 6.8. These proposed unusual pK(a)s have led us to investigate the electrostatic environment of the active site of GART. We have used Poisson-Boltzmann-based electrostatic methods to calculate the pK(a)s of all ionizable groups, using the crystallographic structure of a ternary complex of GART involving the pseudosubstrate 5-deaza-5,6,7,8-THF (5dTHF) and substrate GAR. Theoretical mutation and deletion analogs have been constructed to elucidate pairwise electrostatic interactions between key ionizable sites within the catalytic site. Also, a construct of a more realistic catalytic site including a reconstructed pseudocofactor with an attached formyl group, in an environment with optimal local van der Waals interactions (locally minimized) that imitates closely the catalytic reactants, has been used for pK(a) calculations. Strong electrostatic coupling among catalytic residues His108, Asp144, and substrate GAR was observed, which is extremely sensitive to the initial protonation and imidazole ring flip state of His108 and small structural changes. We show that a proton can be exchanged between GAR and His108, depending on their relative geometry and their distance to Asp144, and when the proton is attached on His108, catalysis could be possible. Using the formylated locally minimized construct of GART, a high pK(a) for His108 was calculated, indicating a protonated histidine, and a low pK(a) for GAR(NH(2)) was calculated, indicating that GAR is in neutral form. Our results are in qualitative agreement with the current mechanistic picture of the catalytic process of GART deduced from the experimental data, but they do not reproduce the absolute magnitude of the pK(a)s extracted from fits of k(cat)-pH profiles, possibly because the static time-averaged crystallographic structure does not describe adequately the dynamic nature of the catalytic site during binding and catalysis. In addition, a strong effect on the pK(a) of GAR(NH(2)) is produced by the theoretical mutations of His108Ala and Asp144Ala, which is not in agreement with the observed insensitivity of the pK(a) of GAR(NH(2)) modeled from the experimental data using similar mutations. Finally, we show that important three-way electrostatic interactions between highly conserved His137, with His108 and Asp144, are responsible for stabilizing the electrostatic microenvironment of the catalytic site. In conclusion, our data suggest that further detailed computational and experimental work is necessary.

Prediction of functionally important residues based solely on the computed energetics of protein structure.

Catalytic and other functionally important residues in proteins can often be mutated to yield more stable proteins. Many of these residues are charged residues that are located in electrostatically unfavorable environments. Here it is demonstrated that because continuum electrostatics methods can identify these destabilizing residues, the same methods can also be used to identify functionally important residues in otherwise uncharacterized proteins. To establish this point, detailed calculations are performed on six proteins for which good structural and mutational data are available from experiments. In all cases it is shown that functionally important residues known to be destabilizing experimentally are among the most destabilizing residues found in the calculations. A larger scale analysis performed on 216 different proteins demonstrates the existence of a general relationship between the calculated electrostatic energy of a charged residue and its degree of evolutionary conservation. This relationship becomes obscured when electrostatic energies are calculated using Coulomb's law instead of the more complete continuum electrostatics method. Finally, in a first predictive application of the method, calculations are performed on three proteins whose structures have recently been reported by a structural genomics consortium.

Identification of protein oligomerization states by analysis of interface conservation.

The discrimination of true oligomeric protein-protein contacts from nonspecific crystal contacts remains problematic. Criteria that have been used previously base the assignment of oligomeric state on consideration of the area of the interface and/or the results of scoring functions based on statistical potentials. Both techniques have a high success rate but fail in more than 10% of cases. More importantly, the oligomeric states of several proteins are incorrectly assigned by both methods. Here we test the hypothesis that true oligomeric contacts should be identifiable on the basis of an increased degree of conservation of the residues involved in the interface. By quantifying the degree of conservation of the interface and comparing it with that of the remainder of the protein surface, we develop a new criterion that provides a highly effective complement to existing methods.

Calculation of weak protein-protein interactions: the pH dependence of the second virial coefficient.

Interactions between proteins are often sufficiently weak that their study through the use of conventional structural techniques becomes problematic. Of the few techniques capable of providing experimental measures of weak protein-protein interactions, perhaps the most useful is the second virial coefficient, B(22), which quantifies a protein solution's deviations from ideal behavior. It has long been known that B(22) can in principle be computed, but only very recently has it been demonstrated that such calculations can be performed using protein models of true atomic detail (Biophys. J. 1998, 75:2469-2477). The work reported here extends these previous efforts in an attempt to develop a transferable energetic model capable of reproducing the experimental trends obtained for two different proteins over a range of pH and ionic strengths. We describe protein-protein interaction energies by a combination of three separate terms: (i) an electrostatic interaction term based on the use of effective charges, (ii) a term describing the electrostatic desolvation that occurs when charged groups are buried by an approaching protein partner, and (iii) a solvent-accessible surface area term that is used to describe contributions from van der Waals and hydrophobic interactions. The magnitude of the third term is governed by an adjustable, empirical parameter, gamma, that is altered to optimize agreement between calculated and experimental values of B(22). The model is applied separately to the proteins lysozyme and chymotrypsinogen, yielding optimal values of gamma that are almost identical. There are, however, clear difficulties in reproducing B(22) values at the extremes of pH. Explicit calculation of the protonation states of ionizable amino acids in the 200 most energetically favorable protein-protein structures suggest that these difficulties are due to a neglect of the protonation state changes that can accompany complexation. Proper reproduction of the pH dependence of B(22) will, therefore, almost certainly require that account be taken of these protonation state changes. Despite this problem, the fact that almost identical gamma values are obtained from two different proteins suggests that the basic energetic formulation used here, which can be evaluated very rapidly, might find use in dynamical simulations of weak protein-protein interactions at intermediate pH values.

Realistic modeling of the denatured states of proteins allows accurate calculations of the pH dependence of protein stability.

Computational techniques based on continuum electrostatics treatments have been successful in predicting and interpreting the pKa values of ionizable amino acids in folded proteins. Despite this progress, efforts to reproduce the pH-dependence of protein stability have met with only limited success: agreement with experimental results has been only qualitative. It has been argued previously that the most likely reason for discrepancies is the presence of residual electrostatic interactions in the unfolded state, which cause pKa values to be shifted from their model compound values. Here we show that by constructing atomistic models of the unfolded state with a simple molecular mechanics protocol that uses the native state as a starting point, much improved reproduction of pH effects on protein stability can be obtained. In contrast, when a fully extended model of the unfolded state is used, no such improvement is obtained, a result that suggests that local interactions with residues nearby in the sequence are not sufficient to properly account for the pKa shifts in the unfolded state. In comparison to model compound values, the pKa values of acidic residues in "native-like" unfolded states are typically found to be shifted downwards by approximately 0.3 pH unit, in good agreement with the average downward shift deduced from experimental measurements. Given its success in the present situation, the protocol employed here for developing simple models of the unfolded state may prove useful in other computer simulation applications.

Computer simulations of actin polymerization can explain the barbed-pointed end asymmetry.

Computer simulations of actin polymerization were performed to investigate the role of electrostatic interactions in determining polymerization rates. Atomically detailed models of actin monomers and filaments were used in conjunction with a Brownian dynamics method. The simulations were able to reproduce the measured barbed end association rates over a range of ionic strengths and predicted a slower growing pointed end, in agreement with experiment. Similar simulations neglecting electrostatic interactions indicate that configurational and entropic factors may actually favor polymerization at the pointed end, but electrostatic interactions remove this trend. This result would indicate that polymerization at the pointed end is not only limited by diffusion, but faces electrostatic forces that oppose binding. The binding of the actin depolymerizing factor (ADF) and G-actin complex to the end of a filament was also simulated. In this case, electrostatic steering effects lead to an increase in the simulated association rate. Together, the results indicate that simulations provide a realistic description of both polymerization and the binding of more complex structures to actin filaments.

Computer simulation of protein-protein association kinetics: acetylcholinesterase-fasciculin.

Computer simulations were performed to investigate the role of electrostatic interactions in promoting fast association of acetylcholinesterase with its peptidic inhibitor, the neurotoxin fasciculin. The encounter of the two macromolecules was simulated with the technique of Brownian dynamics (BD), using atomically detailed structures, and association rate constants were calculated for the wild-type and a number of mutant proteins. In a first set of simulations, the ordering of the experimental rate constants for the mutant proteins was correctly reproduced, although the absolute values of the rate constants were overestimated by a factor of around 30. Rigorous calculations of the full electrostatic interaction energy between the two proteins indicate that this overestimation of association rates results at least in part from approximations made in the description of interaction energetics in the BD simulations. In particular, the initial BD simulations neglect the unfavourable electrostatic desolvation effects that result from the exclusion of high dielectric solvent that accompanies the approach of the two low dielectric proteins. This electrostatic desolvation component is so large that the overall contribution of electrostatics to the binding energy of the complex is unlikely to be strongly favourable. Nevertheless, electrostatic interactions are still responsible for increased association rates, because even if they are unfavourable in the fully formed complex, they are still favourable at intermediate protein-protein separation distances. It therefore appears possible for electrostatic interactions to promote the kinetics of binding even if they do not make a strongly favourable contribution to the thermodynamics of binding. When an approximate description of these electrostatic desolvation effects is included in a second set of BD simulations, the relative ordering of the mutant proteins is again correctly reproduced, but now association rate constants that are much closer in magnitude to the experimental values are obtained. Inclusion of electrostatic desolvation effects also improves reproduction of the experimental ionic strength dependence of the wild-type association rate.

Rapid binding of a cationic active site inhibitor to wild type and mutant mouse acetylcholinesterase: Brownian dynamics simulation including diffusion in the active site gorge.

It is known that anionic surface residues play a role in the long-range electrostatic attraction between acetylcholinesterase and cationic ligands. In our current investigation, we show that anionic residues also play an important role in the behavior of the ligand within the active site gorge of acetylcholinesterase. Negatively charged residues near the gorge opening not only attract positively charged ligands from solution to the enzyme, but can also restrict the motion of the ligand once it is inside of the gorge. We use Brownian dynamics techniques to calculate the rate constant kon, for wild type and mutant acetylcholinesterase with a positively charged ligand. These calculations are performed by allowing the ligand to diffuse within the active site gorge. This is an extension of previously reported work in which a ligand was allowed to diffuse only to the enzyme surface. By setting the reaction criteria for the ligand closer to the active site, better agreement with experimental data is obtained. Although a number of residues influence the movement of the ligand within the gorge, Asp74 is shown to play a particularly important role in this function. Asp74 traps the ligand within the gorge, and in this way helps to ensure a reaction.

The stability of salt bridges at high temperatures: implications for hyperthermophilic proteins.

Salt bridges have been proposed to play a crucial role in promoting hyperthermostability in proteins, yet they appear to make little contribution to protein stability at room temperature. The latter point has been rationalized previously on the basis that the association of two charged molecules to form a salt bridge incurs a substantial desolvation penalty, which is seldom completely compensated by favourable interactions within the salt bridge and with the rest of the protein. Here a continuum solvation model is used to investigate how this same argument applies at temperatures more appropriate to hyperthermophiles. The solvation model employed was previously parameterised to reproduce the hydration free energies of neutral and charged amino acid side-chains in the temperature range from 5-100 degreesC. A key result of the previous work was that the hydration free energies of charged side-chains are more adversely affected by increasing temperature than are the hydration free energies of hydrophobic side-chains of identical size and shape (isosteres). As is shown here, a direct consequence of the temperature dependence of the hydration free energies is that at high temperatures the desolvation penalty for formation of a salt bridge is markedly reduced in magnitude. As a result, the argument that relative to hydrophobic isosteres, salt bridges destabilise proteins, may no longer be true at high temperatures. We demonstrate this point first in the setting of a small model system, but then also show that the same argument is likely to carry over to real proteins. We compare three hyperthermophilic proteins with their mesophilic homologues and find that hydration effects preferentially stabilise the hyperthermophiles at high temperatures. When the hydration effects are incorporated into a model for the free energy of folding of the proteins, it is found that in each case, the hyperthermophile is predicted to remain stable to a temperature 20-25 deg.C higher than the corresponding mesophile. Higher thermal stability for the hyperthermophile is obtained even if the mesophile is more stable at room temperature. The results obtained therefore suggest one possible way in which the apparently destabilising effects of salt bridges at room temperature can be reconciled with their increased abundance in hyperthermophilic proteins.

Electrostatic contributions to the stability of halophilic proteins.

Examination of the first crystal structures of proteins from a halophilic organism suggests that an abundance of acidic residues distributed over the protein surface is a key determinant of adaptation to high-salt conditions. Although one extant theory suggests that acidic residues are favored because of their superior water-binding capacity, it is clear that extensive repulsive electrostatic interactions will also be present in such proteins at physiological pH. To investigate the magnitude and importance of such electrostatic interactions, we conducted a theoretical analysis of their contributions to the salt and pH-dependence of stability of two halophilic proteins. Our approach centers on use of the Poisson-Boltzmann equation of classical electrostatics, applied at an atomic level of detail to crystal structures of the proteins. We first show that in using the method, it is important to account for the fact that the dielectric constant of water decreases at high salt concentrations, in order to reproduce experimental changes in pKa values of small acids and bases. We then conduct a comparison of salt and pH effects on the stability of 2Fe-2S ferredoxins from the halophile Haloarcula marismortui and the non-halophile anabaena. In both proteins, substantial upward shifts in pKa accompany protein folding, though shifts are considerably larger, on average, in the halophile. Upward shifts for basic residues occur because of favorable salt-bridge interactions, whilst upward shifts for acidic residues result from unfavorable electrostatic interactions with other acidic groups. Our calculations suggest that at pH 7 the stability of the halophilic protein is decreased by 18.2 kcal/mol on lowering the salt concentration from 5 M to 100 mM, a result that is in line with the fact that halophilic proteins generally unfold at low salt concentrations. For comparison, the non-halophilic ferredoxin is calculated to be destabilized by only 5.1 kcal/mol over the same range. Analysis of the pH stability curve suggests that lowering the pH should increase the intrinsic stability of the halophilic protein at low salt concentrations, although in practice this is not observed because of aggregation effects. We report the results of a similar analysis carried out on the tetrameric malate dehydrogenase from H. marismortui. In this case, we investigated the salt and pH dependence of the various monomer-monomer interactions present in the tetramer. All monomer-monomer interactions are found to make substantial contributions to the salt-dependence of stability of the tetramer. Excellent agreement is obtained between our calculated results for the stability of the tetramer and experimental results. In particular, the finding that at 4 M NaCl, the tetramer is stable only between pH 4.8 and 10 is accurately reproduced. Taken together, our results suggest that repulsive electrostatic interactions between acidic residues are a major factor in the destabilization of halophilic proteins in low-salt conditions, and that these interactions remain destabilizing even at high salt concentrations. As a consequence, the role of acidic residues in halophilic proteins may be more to prevent aggregation than to make a positive contribution to intrinsic protein stability.

Electrostatic effects in homeodomain-DNA interactions.

We report here an investigation of the role of electrostatics in homeodomain-DNA interactions using techniques based around the use of the Poisson-Boltzmann equation. In the present case such a study is of particular interest, since in contrast to other proteins previously studied with this method, the homeodomain is a small, highly charged protein that forms extensive ion pairs upon binding DNA. We have investigated the salt dependence of the binding constant for specific association and for a variety of models for non-specific association. The results indicate that, in line with the models proposed by Manning and Record, the entropy of counterion release accounts for a significant fraction of the salt dependence of the binding free energy, though this is perhaps due to fortuitous cancellation of other contributing terms. The thermodynamic effects of a number of specific homeodomain mutants were also investigated, and partly rationalized in terms of favorable electrostatic interactions in the major goove of DNA. Investigation of the temperature-dependence of the free energy of association indicates that the electrostatic contributions become increasingly favorable as the temperature rises. For this particular system, however, there appears to be no significant electrostatic contribution to the delta(delta C(p)) of association. Finally, an analysis of the free energy of interaction when the homeodomain is moved ca one Debye length from the DNA suggests that pure electrostatic forces are able to steer the homeodomain into a partially correct orientation for binding to the DNA.

Electrostatic channeling of substrates between enzyme active sites: comparison of simulation and experiment.

Recent simulation work has indicated that channeling of charged substrates between the active sites of bifunctional enzymes or bienzyme complexes can be significantly enhanced by favorable interactions with the electrostatic field of the enzymes. The results of such simulations are expressed in terms of transfer efficiencies, which describe the probability that substrate leaving the active site of the first enzyme will reach the active site of the second enzyme before escaping out into bulk solution. The experimental indicators of channeling, on the other hand, are factors such as a decrease in the transient (lag) time for appearance of the final product of the coupled enzyme reaction or a decrease in the susceptibility of the overall reaction rate to the presence of competing enzymes or competitive inhibitors. The work reported here aims to establish a connection between the transfer efficiencies obtained from simulation, with the above-mentioned experimental observables. This is accomplished by extending previously reported analytical approaches to combine the simulated transfer efficiency with the Michaelis-Menten kinetic parameters Km and Vmax of the enzymes involved; expressions are derived to allow both transient times and steady state rates to be calculated. These results are applied to the two systems that have been studied both theoretically and experimentally. In the first case, that of the bifunctional enzyme dihydrofolate reductase-thymidylate synthase (DHFR-TS), the experimentally observed decrease in transient times is found to be consistent with a transfer efficiency of >/=80%. In the second case, that of a citrate synthase-malate dehydrogenase fusion protein, a transfer efficiency of 73% is consistent with the experimental transient time measurements. Separate and independent analysis of the effects of adding the competing enzyme aspartate aminotransferase gives a transfer efficiency of 69%, in excellent agreement with the transient time results. The transfer efficiencies thus obtained from experimental results are in both cases in good agreement with those obtained from simulations that include electrostatic interactions. One important discrepancy between simulation and experiment, is however, found in the reported effects of adding a competitive inhibitor in the DHFR-TS system: qualitatively different results are expected from the theoretical analysis. A possible reason for this apparent contradiction is discussed.

Evidence for electrostatic channeling in a fusion protein of malate dehydrogenase and citrate synthase.

Brownian dynamics simulations were performed to investigate a possible role for electrostatic channeling in transferring substrate between two of the enzymes of the citric acid cycle. The diffusion of oxaloacetate from one of the active sites of malate dehydrogenase (MDH) to the active sites of citrate synthase (CS) was simulated in the presence and absence of electrostatic forces using a modeled structure for a MDH-CS fusion protein. In the absence of electrostatic forces, fewer than 1% of substrate molecules leaving the MDH active site are transferred to CS. When electrostatic forces are present at zero ionic strength however, around 45% of substrate molecules are successfully channeled. As expected for an electrostatic mechanism of transfer, increasing the ionic strength in the simulations reduces the calculated transfer efficiency. Even at 150 mM however, the inclusion of electrostatic forces results in an increase in transfer efficiency of more than 1 order of magnitude. The simulations therefore provide evidence for the involvement of electrostatic channeling in guiding substrate transfer between two of the enzymes of the citric acid cycle. Similar effects may operate between other members of the citric acid metabolon.

Electrostatic channeling in the bifunctional enzyme dihydrofolate reductase-thymidylate synthase.

The bifunctional enzyme dihydrofolate reductase-thymidylate synthase (DHFR-TS) carries out two distinct reactions, with the dihydrofolate produced by the TS-catalyzed reaction acting as the substrate for the DHFR-catalyzed reaction. Brownian dynamics simulation techniques were used to investigate the possible role of electrostatics in determining efficient channeling of the substrate, by explicitly simulating substrate diffusion between the two active sites. With a substrate charge of -2, almost all (> 95%) substrate molecules leaving the TS active site reached the DHFR active site at zero ionic strength. Under the same conditions, but in the absence of electrostatic effects, successful channeling was reduced to only around 6%: electrostatic effects therefore appear essential to explain the efficient channeling observed experimentally. The importance of substrate charge, the relative contributions of specific basic residues in the protein, the role played by the second monomer of the dimer in channeling and the effects of changing ionic strength were all investigated. Simulations performed for substrate transfer in the opposite direction suggest that channeling in DHFR-TS is not strongly directional and that the role of electrostatics is perhaps more one of restricting diffusion of the substrate than one of actively guiding it from the TS to the DHFR active site. The results demonstrate that electrostatic channeling can be a highly efficient means of transferring charged substrates between active sites in solvent-exposed environments.

Theoretical studies of the intercalation of 9-hydroxyellipticine in DNA.

Extensive molecular dynamics (MD) simulations have been used to investigate the intercalative binding of 9-hydroxyellipticine to the DNA oligonucleotide d(ATATATATATAT)2. Four independent simulations differing in the initial orientation of the drug at the intercalation site were carried out, and compared both with each other and a control simulation of the free DNA sequence. The structure of the latter was compared with structures obtained from x-ray crystallography and nmr spectroscopy, as well as the theoretically derived "alternating B-DNA" model [A. Klug et al. (1979), Journal of Molecular Biology, Vol. 131, p. 669]. The alternation of twist angles observed in experimental structures was reproduced in the simulation. All four independent simulations of the drug-DNA intercalation complex converged in placing the pyridine ring of the ellipticine chromophore in the major groove; in one case this involved a 180 degrees rotation of the drug at the intercalation site. At a more detailed level, the drug is seen to be capable of adopting several distinct orientations, each stable over a period of hundreds of pico-seconds. Despite the presence of several polar groups in the drug, however, no direct hydrogen bonding to the DNA occurs; instead, interactions between the methyl groups of the drug and the thymine bases at the intercalation site appear important in determining the orientational preferences of the drug. Comparison of the intercalation complexes with the free DNA sequence shows a degree of unwinding resulting from intercalation, in good agreement with experimental results, but spread over the three central base-pair steps, not confined to the intercalation site itself. Measurements of torsional rigidity indicate only a slight stiffening of the DNA restricted to the immediate site of intercalation. The structures obtained from the MD simulations were used to calculate theoretical CD spectra, with separate simulations giving very different results. This appears to indicate that given an accurate assignment of the main electronic transition dipole moment of the ellipticine chromophore, discrimination of the more realistic binding geometries may be possible. The relative merits of the various drug orientations observed in the simulations are discussed and a perpendicular orientation of the drug at the intercalation site is considered to be the most consistent with experimental data. While the simulations themselves represent a total of over 2 ns, however, the differences apparent between independent runs indicate that longer simulation times will be required before a complete, unequivocal view of DNA intercalation is obtained.

SUBCUR: visualization of structural differences between DNA duplexes.

A computer program, SUBCUR, is described which permits analysis and rapid identification of geometrical differences and patterns of variance between two DNA duplexes. The program is compatible with the CURVES 3.1 package and allows graphical visualization of the structural differences. Examples are provided which illustrate the applicability of the program in analyzing the different backbone conformations of two helices and the different curvatures of two helices.

A binding mode of lambda-[tris(1,10-phenanthroline)ruthenium(II)]2+ exhibiting preference for purine-3',5'-pyrimidine sites of DNA.

Molecular mechanics calculations and molecular dynamics simulations have been used to study the binding of the partially inserted major groove complex of Lambda-[Ru(1,10-phenanthroline)3]2+ with DNA. Energy refinements of this complex showed a clear preference for binding at purine-3',5'-pyrimidine sites over pyrimidine-3',5'-purine sites. The basis for this difference is shown to be a slight change in the binding orientation induced by interchanging the purine and pyrimidine bases. This in turn provides for a better secondary interaction with the helix backbone at a point beyond the immediate binding site. It is this secondary interaction that provides the additional energetic stabilisation for complexes formed at purine-3',5'-pyrimidine sites. Molecular dynamics simulations including explicit representation of solvent support these conclusions and provide an insight into the positional stability of the ligand at a particular site. Repuckering of specific deoxyribose rings to the C3'-endo conformation seems to be an important feature of the DNA/ligand complex.

Sequence selective binding to the DNA major groove: tris(1,10-phenanthroline) metal complexes binding to poly(dG-dC) and poly(dA-dT).

Molecular modelling and energy minimisation calculations that incorporate solvent effects have been used to investigate the complexation of delta and lambda-[Ru(1,10-phenanthroline]2+ to DNA. The most stable binding geometry for both enantiomers is one in which a phenanthroline chelate is positioned in the major groove. The chelate is partially inserted between neighbouring base pairs, but is not intercalated. For delta, though not for lambda, a geometry with two chelates in the major groove is only slightly less favourable. Minor groove binding is shown to be no more favourable than external electrostatic binding. The optimised geometries of the DNA/[Ru(1,10-phenanthroline]2+ complexes enable published linear dichroism spectra to be used to determine the percentage of each enantiomer in the two most favourable major groove sites. For delta 57 +/- 15% and for lambda 82 +/- 7% of bound molecules are in the partially inserted site.