Conformational Studies of Nucleic Acids: IV. The Conformational Energetics of Oligonucleotides: d(ApApApA) and ApApApA

Abstract

Utilizing a new method for modeling furanose pseudorotation (D. A Pearlman and S.-H. Kim, J. Biomol. Struct. Dyn. 3, 85 (1985)) and the empirical multiple correlations between nucleic acid torsion angles we derived in the previous report (D. A Pearlman and S.-H. Kim, previous paper in this issue), we have made an energetic examination of the entire conformational spaces available to two nucleic acid oligonucleotides: d(ApApApA) and ApApApA The energies are calculated using a semi-empirical potential function. From the resulting body of data, energy contour map pairs (one for the DNA molecule, one for the RNA structure) have been created for each of the 21 possible torsion angle pairs in a nucleotide repeating unit. Of the 21 pairs, 15 have not been reported previously. The contour plots are different from those made earlier in that for each point in a particular angle-angle plot, the remaining five variable torsion angles are rotated to the values which give a minimum energy at this point. The contour maps are overall quite consistent with the experimental distribution of oligonucleotide data. A number of these maps are of particular interest: δ (C5′-C4′-C3′-03′)χ (04′-C1′-N9- C4), where the energetic basis for an approximately linear δ-χ correlation can be seen; ζ (C3′- 03′-P-05′)-δ, in which the experimentally observed linear correlation between ζ and δ in DNA (220° < ζ <280°) is clearly predicted; ζ-ε (C4′-C3′-03′-P), which shows that e increases with decreasing ζ <260°; α (03′-P-05′-C5′)-γ (05′-C5′-C4′-C3′) where a clear linear correlation between these angles is also apparent, consistent with experiment; and several others. For the DNA molecule studied here, the sugar torsion Ô is predicted to be the most flexible, while for the RNA molecule, the greatest amount of flexibility is expected to reside in a and y. Both the DNA and RNA molecules are predicted to be highly polymorphic. Complete energy minimization has been performed on each of the minima found in the energy searches and the results further support this prediction. Possible pathways for B-form to A-form DNA interconversion suggested by the results of this study are discussed. The results of these calculations support use of the new sugar modeling technique and torsion angle correlations in future conformational studies of nucleic acids.     

Conformational studies of nucleic acids: IV. The conformational energetics of oligonucleotides: d(ApApApA) and ApApApA

Utilizing a new method for modeling furanose pseudorotation (D. A. Pearlman and S.-H. Kim, J. Biomol. Struct. Dyn. 3, 85 (1985)) and the empirical multiple correlations between nucleic acid torsion angles we derived in the previous report (D. A. Pearlman and S.-H. Kim, previous paper in this issue), we have made an energetic examination of the entire conformational spaces available to two nucleic acid oligonucleotides: d(ApApApA) and ApApApA. The energies are calculated using a semi-empirical potential function. From the resulting body of data, energy contour map pairs (one for the DNA molecule, one for the RNA structure) have been created for each of the 21 possible torsion angle pairs in a nucleotide repeating unit. Of the 21 pairs, 15 have not been reported previously. The contour plots are different from those made earlier in that for each point in a particular angle-angle plot, the remaining five variable torsion angles are rotated to the values which give a minimum energy at this point. The contour maps are overall quite consistent with the experimental distribution of oligonucleotide data. A number of these maps are of particular interest: delta (C5'-C4'-C3'-O3')-chi (O4'-C1'-N9-C4), where the energetic basis for an approximately linear delta-chi correlation can be seen: zeta (C3'-O3'-P-O5')-delta, in which the experimentally observed linear correlation between zeta and delta in DNA(220 degrees less than zeta less than 280 degrees) is clearly predicted; zeta-epsilon (C4'-C3'-O3'-P), which shows that epsilon increases with decreasing zeta less than 260 degrees; alpha (O3'-P-O5'-C5')-gamma (O5'-C5'-C4'-C3') where a clear linear correlation between these angles is also apparent, consistent with experiment; and several others. For the DNA molecule studied here, the sugar torsion delta is predicted to be the most flexible, while for the RNA molecule, the greatest amount of flexibility is expected to reside in alpha and gamma. Both the DNA and RNA molecules are predicted to be highly polymorphic. Complete energy minimization has been performed on each of the minima found in the energy searches and the results further support this prediction. Possible pathways for B-form to A-form DNA interconversion suggested by the results of this study are discussed. The results of these calculations support use of the new sugar modeling technique and torsion angle correlations in future conformational studies of nucleic acids.     

Structural organization of guanosine derivatives in dilute solutions: small angle neutron scattering analysis

Guanosine derivatives, dissolved in water, can form cholesteric and hexagonal mesophases. The common structural unit is a chiral rod-shaped aggregate consisting of a stack of Hoogsten-bonded guanosine tetrameric disks. In order to elucidate the self-association process, we decided to investigate, by small-angle neutron scattering, the structural properties of d(pG), d(GpG), d(GpGpG), d(GpGpGpG) and d(GpGpGpGpGpG) derivatives in very dilute solutions. Under our experimental conditions only d(pG) seems not to form detectable particles. On the other hand, the results for the other derivatives indicate that cylindrical aggregates, having a 10 cross-section gyration radius and a length of about 70 Å, exist in the isotropic phase. According to the structure of the hexagonal and cholesteric phases, we fitted the experimental data by using a model of rod-shaped aggregates formed by stacking about 18 to 20 guanosine tetramers. Moreover, from the measurement of the concentration of scattering particles, we deduced that guanosine derivatives are only partially aggregated, depending on their ability to form mesophases. Correspondence to: P. Mariani     

Structure and Dynamics of Double Helical DNA in Torsion Angle Hyperspace: A Molecular Mechanics Approach

Abstract

Analysis of the conformational space populated by the torsion angles and the correlation between the conformational energy and the sequence of DNA are important for fully understanding DNA structure and function. Presence of seven variable torsion angles about single covalent bonds in DNA main chain puts a big challenge for such analysis. We have carried out restrained energy minimization studies for four representative dinucleosides, namely d(ApA):d(TpT), d(CpG):d(CpG), d(GpC):d(GpC) and d(CpA):d(TpG) to determine the energy hyperspace of DNA in context to the values of the torsion angles and the structural properties of the DNA conformations populating the favorable regions of this energy hyperspace. The torsion angles were manipulated by constraining their values at the reference points and then performing energy minimization. The energy minima obtained on the potential energy contour plots mostly correspond to the conformations populated in crystal structures of DNA. Some novel favorable conformations that are not present in crystal structure data are also found. The plots also suggest few low energy routes for conformational transitions or the associated energy barrier heights. Analyses of base pairing and stacking possibility reveal structural changes accompanying these transitions as well as the flexibility of different base steps towards variations in different torsion angles.     

Model for the conformational analysis of hydrated peptides. Effect of hydration on the conformational stability of the terminally blocked residues of the 20 naturally occurring amino acids

The effects of hydration are included in empirical conformational energy computations on oligopeptides by means of a modified hydration-shell model. Free energy terms are introduced to account for “specific hydration” due to water–solute hydrogen bonding and for “nonspecific hydration” describing the interaction of the solute with water molecules in a first-neighbor shell. The dielectric constant has been doubled (over the value used for calculations in the absence of water) to take into account the presence of solvent. Computations were carried out for the N-acetyl-N′-methylamides of the 20 naturally occurring amino acids. Conformational energy maps are compared with similar maps calculated in the absence of hydration. Minimum-energy conformations are located and compared with the corresponding minima for unhydrated peptides in terms of ordering with respect to potential energy, the dihedral angles at the minima, and the presence of intramolecular hydrogen bonds. The Boltzmann factors for various conformational regions are altered significantly on hydration in some cases. These changes can be explained in terms of differences in the hydration free energy terms for various conformations.     

Influence of the base sequence on the conformational behaviour of DNA polynucleotides in solution

NMR studies were carried out on samples of the non-self-complementary tetramers d(C-A-C-A), d(T-G-T-G), d(G-A-G-A) and d(T-C-T-C) and of 1:1 mixtures of the complementary tetramers d(C-A-C-A) X d(T-G-T-G) and d(G-A-G-A) X d(T-C-T-C) at two DNA concentrations and of the self-complementary octamers d(C-A-C-A-T-G-T-G) and d(G-A-G-A-T-C-T-C). Assignments, based upon one-dimensional NOE and homonuclear-decoupling and two-dimensional correlated and NOE spectroscopies are given of the resonances of most of the base and sugar protons. Chemical shift vs temperature profiles, constructed for all samples, yielded insight into the temperature- and concentration-dependent conformational behaviour of the compounds and were used to obtain thermodynamic parameters pertaining to the stacked-single-strand----random-coil and duplex----random-coil equilibria. Vicinal proton-proton couplings were analyzed in terms of the conformation of the deoxyribose rings in the single-stranded tetramers and duplexed octamers. The NOE patterns, chemical shift profiles, imino-proton resonances and coupling data revealed that the compounds adopt B-DNA-like structures. The ratio duplexed/stacked-single-strand/random coil depends upon external conditions as well as upon base sequence. The thermodynamic data indicate that: in terms of single-helical stacking, the R-R steps (Tm 321-328 K) appear more stable than the Y-R or R-Y steps (Tm 308-316 K) and the Y-Y steps score least (Tm 290-300 K), and the duplexes consisting of alternating, d(Y-R)n, strands are more stable, in terms of delta H degrees, compared to the d(R-R)n X d(Y-Y)n duplexes. The analyses of the couplings demonstrated that the sugars of the single-stranded tetramers and duplexed octamers occur as a blend of N- and S-type conformers, with a preference for the S-type (C2'-endo) sugar conformation: upon duplex formation, no significant shift in the N-type/S-type ratio was observed. The fraction S-type sugar conformation of a given residue, %S, in the stacked-single strands was found to depend upon the nature of its own base and that of the adjacent residues: sugars in an R-R stretch display high values of %S (90-100), whereas those in Y-Y stretches show relatively low values (approximately equal to 65).(ABSTRACT TRUNCATED AT 400 WORDS)     

Thermodynamics of the first transition in writhe of a small circular DNA by Monte Carlo simulation

Monte Carlo simulations are employed to investigate the thermodynamics of the first transition in writhe of a circular model filament corresponding to a 468 base-pair DNA. Parameters employed in these simulations are the torsional rigidity, C = 2.0 × 10⁻¹⁹ dyne cm², and persistence length, P = 500 Å. Intersubunit interactions are modeled by a screened Coulomb potential. For a straight line of subunits this accurately approximates the nonlinear Poisson-Boltzmann potential of a cylinder with the linear charge density of DNA. Curves of relative free energy vs writhe at fixed linking difference (Δ1) exhibit two minima, one corresponding to slightly writhed circles and one to slightly underwrithed figure-8's, whenever Δ1 lies in the transition region. The free energies of the two minima are equal when Δ1_c = 1.35, which defines the midpoint of the transition. At this midpoint, the free energy barrier between the two minima is found to be ΔG_bar = (0.20) k_BT at 298 K. Curves of mean potential energy vs writhe at fixed linking difference similarly exhibit two minima for Δ1 values in the transition region, and the two minimum mean potential energies are equal when Δ1 = 1.50. At the midpoint writhe, Δ1_c = 1.35, the difference in mean potential energy between the minimum free energy figure-8 and circle states is (1.3) k_BT, and the difference in their entropies is 1.3 k_B. Thus, the entropy of the minimum free energy figure-8 state significantly exceeds that of the circle at the midpoint of the transition. The first transition in writhe is found to occur over a rather broad range of Δ1 values from 0.85 to 1.85. The twist energy parameter (E_T), which governs the overall free energy of supercoiling, undergoes a sigmoidal decrease, while the translational diffusion coefficient undergoes a sigmoidal increase, over this same range. The static structure factor exhibits an increase, which reflects a decrease in radius of gyration associated with the circle to figure-8 transition. © 1996 John Wiley & Sons, Inc.     

Modeling DNA Hydration: Comparison of Calculated and Experimental Hydration Properties of Nuclic Acid Bases

Abstract

Hydration properties of individual nucleic acid bases were calculated and compared with the available experimental data. Three sets of classical potential functions (PF) used in simulations of nucleic acid hydration were juxtaposed: (i) the PF developed by Poltev and Malenkov (PM), (ii) the PF of Weiner and Kollman (WK), which together with Jorgensen's TIP3P water model are widely used in the AMBER program, and (HI) OPLS (optimized potentials for liquid simulations) developed by Jorgensen (J). The global minima of interaction energy of single water molecules with all the natural nucleic acid bases correspond to the formation of two water-base hydrogen bonds (water bridging of two hydrophilic atoms of the base). The energy values of these minima calculated via PM potentials are in somewhat better conformity with mass-spectrometric data than the values calculated via WK PF. OPLS gave much weaker water-base interactions for all compounds considered, thus these PF were not used in further computations. Monte Carlo simulations of the hydration of 9- methyladenine, 1-methyluracil and 1-methylthymine were performed in systems with 400 water molecules and periodic boundary conditions. Results of simulations with PM potentials give better agreement with experimental data on hydration energies than WK PF. Computations with PM PF of the hydration energy of keto and enol tautomers of 9-methyl- guanine can account for the shift in the tautomeric equilibrium of guanine in aqueous media to a dominance of the keto form in spite of nearly equal intrinsic stability of keto and enol tautomers. The results of guanine hydration computations are discussed in relation to mechanisms of base mispairing errors in nucleic acid biosynthesis. The data presented in this paper along with previous results on simulation of hydration shell structures in DNA duplex grooves provide ample evidence for the advantages of PM PF in studies of nucleic-acid hydration.     

Solvent effects on the conformational preferences of model peptoids. MP2 study

The influence of aqueous environment on the main‐chain conformation (ω₀, ?, and ψ dihedral angles) of two model peptoids: N‐acetyl‐N‐methylglycine N’‐methylamide (Ac‐N(Me)‐Gly‐NHMe) ( 1 ) and N‐acetyl‐N‐methylglycine N’,N’‐dimethylamide (Ac‐N(Me)‐Gly‐NMe₂) ( 2 ) was investigated by MP2/6‐311++G(d,p) method. The Ramachandran maps of both studied molecules with cis and trans configuration of the N‐terminal amide bond in the gas phase and in water environment were obtained and all energy minima localized. The polarizable continuum model was applied to estimate the solvation effect on conformation. Energy minima of the Ac‐N(Me)‐Gly‐NHMe and Ac‐N(Me)‐Gly‐NMe₂ have been analyzed in terms of the possible hydrogen bonds and C = O dipole attraction. To validate the theoretical results obtained, conformations of the similar structures gathered in the Cambridge Crystallographic Data Centre were analyzed. Obtained results indicate that aqueous environment in model peptoids 1 and 2 favors the conformation F (? and ψ = ?70º, 180º), and additionally significantly increases the percentage of structures with cis configuration of N‐terminal amide bond in studied compounds. Copyright © 2014 European Peptide Society and John Wiley & Sons, Ltd.     

Simulation of Interactions between Nucleic Acid Bases by Refined Atom-Atom Potential Functions

Abstract

Energy of interaction between nitrogen bases of nucleic acids has been calculated as a function of parameters determining the mutual position of two bases. Refined atom-atom potential functions are suggested. These functions contain terms proportional to the first (electrostatics), sixth (or tenth for the atoms forming a hydrogen bond) and twelfth (repulsion of all atoms) powers of interatomic distance. Calculations have shown that there are two groups of minima of the base interaction energy. The minima of the first group correspond to coplanar arrangement of the base pairs and hydrogen bond formation. The minima of the second group correspond to the position of bases one above the other in almost parallel planes. There are 28 energy minima corresponding to the formation of coplanar pairs with two (three for the G:C pair) almost linear N-H … O and (or) N-H … N hydrogen bonds. The position of nitrogen bases paired by two such H-bonds in any crystal of nucleic acid component, in polynucleotide complexes and in tRNA is close to the position in one of these minima. Besides, for each pair there are energy minima corresponding to the formation of a single N-H … O or N-H … N and one C-H … O or C-H … N hydrogen bond. The form of potential surface in the vicinity of minima has been characterized. The results of calculations agree with the experimental data and with more rigorous calculations based on quantum- mechanical approach.     

Accurate ab initio calculation of the Ar–CF4 intermolecular potential energy surface

Ar–CF₄ intermolecular interaction potential is studied by ab initio calculations at the MP2 and CCSD(T) levels of theory containing the so-called bond functions ({3s3p2d1f} basis set was chosen) both with and without a correction for the basis-set superposition error. The calculations were performed with Dunning's correlation consistent basis sets (aug-cc-pVXZ, X = D, T, Q, 5) to extrapolate the Ar–CF₄ potential energy minimum and intermolecular distance to their complete basis set (CBS) limits. It is shown that the addition of bond functions results in a dramatic improvement in the convergence of the calculated interaction energies at the MP2/aug-cc-pVTZ level. The MP2/{3s3p2d1f}-aug-cc-pVTZ potential energy surface even approaches the CCSD(T)/aug-cc-pVQZ potential energy surface. The potential energy minima and the intermolecular distances are both significantly closer to the CBS limit when using the bond functions, and it implies that adding bond functions in the calculation has a great effect on the interaction energies. We also find that with bond functions included in the CCSD(T)/aug-cc-pVDZ model chemistry, the potential energy minima are extremely close to the CBS limit and are better than the CCSD(T)/aug-cc-pVQZ values. Several levels of theory described in the text were used to determine pairwise analytic potential energy surfaces for Ar+CF₄. The analytic potential energy surfaces are in very good agreement with the ab initio values.     

Influence of local interactions on protein structure. I. Conformational energy studies of N-acetyl-N′-methylamides of pro-X and X-pro dipeptides

Conformational energy calculations using an Empirical Conformational Energy Program for Peptides (ECEPP) were carried out on the N-acetyl-N′-methylamides of Pro-X, where X = Ala, Asn, Asp, Gly, Leu, Phe, Ser, and Val, and of X-Pro, where X = Ala, Asn, Gly, and Pro. The conformational energy was minimized from starting conformations which included all combinations of low-energy single-residue minima and several standard bend structures. It was found that almost all resulting minima are combinations of low-energy single-residue minima, suggesting that intra residue interactions predominate in determining conformation. The calculations also indicate, however, that inter residue interactions can be important. In addition, librational entropy was found to influence the relative stabilities of some minima. Because of the existence of 10–100 low-energy minima for each dipeptide, the normalized statistical weight of an individual minimum rarely exceeds 0.3, suggesting that these dipeptides have considerable conformational flexibility and exist as statistical ensembles of low-energy structures. The propensity of each dipeptide to form bend conformations was calculated, and the results were compared with available experimental data. It was found that bends are favored in Pro-X dipeptides because ?_Pro is fixed by the pyrrolidine ring in a conformation which is frequently found in bends, but that bends are not favored in X-Pro dipeptides because interactions between the X residue and the pyrrolidine ring restrict the X residue to conformations which are not usually found in bends.     

Mass Spectrometry Study and Infrared Spectroscopy of the Complex Between Camphor and the Two Enantiomers of Protonated Alanine: The Role of Higher‐Energy Conformers in the Enantioselectivity of the Dissociation Rate Constants

The properties of the protonated complexes built from S camphor and R or S alanine were studied in a Paul ion trap at room temperature by collision‐induced dissociation (CID) and infrared multiple‐photon dissociation spectroscopy (IRMPD), as well as molecular dynamics and ab initio calculations. While the two diastereomer complexes display very similar vibrational spectra in the fingerprint region, in line with similar structures, and almost identical calculated binding energies, their collision‐induced dissociation rates are different. Comparison of the IRMPD results to computed spectra shows that the SS and SR complexes both contain protonated alanine strongly hydrogen‐bonded to the keto group of camphor. The floppiness of this structure around the NH⁺…O = C hydrogen bond results in a complex potential energy surface showing multiple minima. Calculating the dissociation rate constant within the frame of the transition state theory shows that the fragmentation rate larger for the heterochiral SR complex than the homochiral SS complex can be explained in terms of two almost isoenergetic low‐energy conformers in the latter that are not present for the former. Chirality 25:436‐443, 2013. © 2013 Wiley Periodicals, Inc.     

Hybrid MC/QC simulations of water-assisted proton transfer in nucleosides. Guanosine and its analog acyclovir

To provide an in-depth insight into the molecular basis of spontaneous tautomerism in DNA and RNA base pairs, a hybrid Monte Carlo (MC)–quantum chemical (QC) methodology is implemented to map two-dimensional potential energy surfaces along the reaction coordinates of solvent-assisted proton transfer processes in guanosine and its analog acyclovir in aqueous solution. The solvent effects were simulated by explicit inclusion of water molecules that model the relevant part of the first hydration shell around the solute. The position of these water molecules was estimated by carrying out a classical Metropolis Monte Carlo simulation of dilute water solutions of the guanosine (Gs) and acyclovir (ACV) and subsequently analyzing solute–solvent intermolecular interactions in the statistically-independent MC-generated configurations. The solvent-assisted proton transfer processes were further investigated using two different ab initio MP2 quantum chemical approaches. In the first one, potential energy surfaces of the ‘bare’ finite solute–solvent clusters containing Gs/ACV and four water molecules (MP2/6-31+G(d,p) level) were explored, while within the second approach, these clusters were embedded in ‘bulk’ solvent treated as polarizable continuum (C-PCM/MP2/6-31+G(d,p) level of theory). It was found that in the gas phase and in water solution, the most stable tautomer for guanosine and acyclovir is the 1H-2-amino-6-oxo form followed by the 2-amino-6-(sZ)-hydroxy form. The energy barriers of the water-assisted proton transfer reaction in guanosine and in acyclovir are found to be very similar – 11.74 kcal mol^?1 for guanosine and 11.16 kcal mol^?1 for acyclovir, and the respective rate constants (k = 1.5?×?10¹ s^?1, guanosine and k = 4.09?×?10¹ s^?1, acyclovir), are sufficiently large to generate the 2-amino-6-(sZ)-hydroxy tautomer. The analysis of the reaction profiles in both compounds shows that the proton transfer processes occur through the asynchronous concerted mechanism.     

Conformational studies on polynucleotide chains. III. Intramolecular energy maps and comparison with experiments

Conformational energy maps for the four combinations of two consecutive torsional angles of the backbone structure of polydeoxyribonucleotides are presented. Both the C(2′)-endo and the C(3′)-endo conformation of sugar rings were considered. The energies were evaluated with an analytical expression representing the best fit to ab initio energies computed in the Hartree-Fock approximation, and consisting of a contribution from nonbonded interactions of the Lennard-Jones 6-12 type and an intrinsic torsional potential. It is shown that the minima of these maps are in excellent agreement with the most stable conformations as obtained from x-ray crystallographic analysis of nucleic acids and polynucleotides.     

Theoretical ab Initio Study of the Effects of Methylation on Structure and Stability of G:C Watson-Crick Base Pair

Abstract

Methylation of DNA occurs most readily at N(3), N(7), and O(6) of purine bases and N(3) and O(2) of pyrimidines. Methylated bases are continuously formed through endogenous and exogenous mechanisms. The results of a theoretical ab initio study on the methylation of G:C base pair components are reported. The geometries of the local minima were optimized without symmetry restrictions by the gradient procedure at DFT level of theory and were verified by energy second derivative calculations. The standard 6–31G(d) basis set was used. The single-point calculations have been performed at the MP2/6–31G(d,p), MP2/6–31₊₊G(d,p), and MP2/6–311₊₊G(2d,2p) levels of theory. The geometrical parameters, relative stability and counterpoise corrected interaction energies are reported. Also, using a variation-perturbation energy decomposition scheme we have found the vital contributions to the total interaction energy.     

Computer simulation of the conformational properties of retro–inverso peptides. I. Empirical force field calculations of rigid and flexible geometries of N-acetylglycine-N′- methylamide,bis(acetamido) methane,and N,N′- dimethylmalonamide and their corresponding Cα-methylated analogs

Rigid and flexible geometry calculations are described for N-acetylglycine-N′-methylamide, N-acetylalanine-N′-methylamide, and their retro-inverso analogs, bis(acetamido) methane, 1,1-bis(acetamido) ethane, N,N′-dimethylmalonamide, and N,N′-dimethyl-2-methyl-malonamide. The significance of relaxing all degrees of freedom, especially angular flexibility is demonstrated. The flexible geometry approach yields energy maps similar to those from rigid geometry, but the energy barriers between minima are substantially reduced, leading in general, to more probable transitions and a higher volume of accessible conformational space. Whereas the glycine and alanine derivatives exhibit their lowest energy minima in the C region, the gem-diaminoalkyl and malonyl residues show their lowest minima in the “α-helical” regions. With respect to the effect of side chains (H versus CH₃), the greatest conformational influence appears with the gem-diaminoalkyl residues. These results indicate significantly different conformational behavior of retro peptides and the implications of these pairwise incorporations of retro-inverso residues in peptide chains, are discussed.     

An optimized potential function for the calculation of nucleic acid interaction energies I. Base stacking

An optimized potential function for base-stacking interaction is constructed. Stacking energies between the complementary pairs of a dimer are calculated as a function of the rotational angle and separation distance. Using several different sets of atomic charges, the electrostatic component in the monopole-monopole approximation (MMA) is compared to the more refined segmented multipole–multipole representation (SMMA); the general features of the stacking minima are found to be correctly reproduced with IEHT or CNDO atomic charges. The electrostatic component is observed to control the location of stacking minima. The MMA, in general, is not a reliable approximation of the SMMA in regions away from minima; however, the MMA is reliable in predicting the location and nature of stacking minima. The attractive part of the Lennard-Jones 6–12 potential is compared to and parameterized against the expression for the second-order interaction terms composed of multipole-bond polarizability for the polarization energy and transition-dipole bond polarizabilities for approximation of the dispersion energy. The repulsive part of the Lennard-Jones potential is compared to a Kitaygorodski-type repulsive function; changing the exponent from its usual value of 12 to 11.7 gives significantly better agreement with the more refined repulsive function. Stacking minima calculated with the optimized potential method are compared with various perturbation-type treatments. The optimized potential method yields results that compare as well with melting data as do any of the more recent and expensive perturbation methods.     

Binding of Hoechst 33258 and its Derivatives to DNA

Abstract

In the present work, we employed UV-VIS spectroscopy, fluorescence methods, and circular dichroism spectroscopy (CD) to study the interaction of dye Hoechst 33258, Hoechst 33342, and their derivatives to poly[d(AT)]·poly[d(AT)], poly(dA)·poly(dT), and DNA dodecamer with the sequence 5′-CGTATATATACG-3′. We identified three types of complexes formed by Hoechst 33258, Hoechst 33342, and methylproamine with DNA, corresponding to the binding of each drug in monomer, dimer, and tetramer forms. In a dimer complex, two dye molecules are sandwiched in the same place of the minor DNA groove. Our data show that Hoechst 33258, Hoechst 33342, and methylproamine also form complexes of the third type that reflects binding of dye associates (probably tetramers) to DNA. Substitution of a hydrogen atom in the ortho position of the phenyl ring by a methyl group has a little effect on binding of monomers to DNA. However it reduces strength of binding of tetramers to DNA. In contrast, a Hoechst derivative containing the ortho-isopropyl group in the phenyl ring exhibits a low affinity to poly(dA)·poly(dT) and poly[d(AT)]·poly[d(AT)] and binds to DNA only in the monomer form. This can be attributed to a sterical hindrance caused by the ortho-isopropyl group for side-by-side accommodation of two dye molecules in the minor groove. Our experiments show that mode of binding of Hoechst 33258 derivatives and their affinity for DNA depend on substituents in the ortho position of the phenyl ring of the dye molecule. A statistical mechanical treatment of binding of Hoechst 33258 and its derivatives to a polynucleotide lattice is described and used for determination of binding parameters of Hoechst 33258 and its derivatives to poly[d(AT)]·poly[d(AT)] and poly(dA)·poly(dT).