Investigation of DNA dynamics and drug-DNA interaction by steady state fluorescence anisotropy.

We have used steady-state fluorescence polarization anisotropy (FPA) of ethidium probe molecules bound to DNA to investigate DNA-DNA interactions and the effect of high densities of intercalating drugs on the internal motions of DNA responsible for depolarization of the ethidium fluorescence. To calibrate the method, we examined the effect of DNA length on (FPA) using DNA varying in size from 10-150 base pair. The association of approximately 30 base pair DNA at high concentrations was then detected by its effect on (FPA). With sample concentrations approaching those commonly used in various physical experiments (NMR, Raman) significant DNA-DNA interactions are observed. With high molecular weight DNA (greater than 500 base pair), the limiting value of the (FPA) (0.23) is due to internal motions of the DNA (and bound chromophores). The (FPA) of ethidium probe molecules (1 drug/200 base pair) is unaffected by the addition of high levels (1 drug/2 base pair) proflavine. This indicates that either the elastic properties of DNA are unaffected by high densities of intercalated drug or that the depolarization of the ethidium fluorescence is due to highly localized motions of the base pairs that are unperturbed by binding of drugs at neighboring sites.     

Thermodynamics of drug-DNA interactions: entropy-driven intercalation and enthalpy-driven outside binding in the ellipticine series

Viscosimetric and kinetic results allow one to characterize three modes of DNA binding in the ellipticine series: (1) Ellipticine and its 9 methoxy derivative, which present maximal DNA lengthening properties and bind DNA through a single step mechanism, can be considered as pure intercalators. (2) Ellipticinium derivatives and short-chain substituted oxazolopyridocarbazoles, which present intermediate DNA lengthening properties, bind DNA through a two-step mechanism, one being intercalation. (3) Long-chain substituted oxazolopyridocarbazole derivatives, which display the smallest DNA lengthening properties, bind DNA through a single-step mechanism, probably resulting from an outside binding mode. The viscosimetric and kinetic results are compared with the thermodynamic results obtained from the temperature dependence of the binding constants. It appears that drugs binding on the outside of the DNA double helix tend to have large enthalpy and small entropy contributions, whereas pure intercalating drugs have contributions from both enthalpy and entropy, with entropy dominating by about 2:1. Drugs showing two binding modes exhibit a continuum between the aforementioned extremes, with no breaks in behavior. From this comparison, a correlation between thermodynamic data and DNA binding modes is proposed. Possible molecular implications of both enthalpy and entropy to DNA binding free energy are discussed.     

Conformation and dynamics of dimeric ureido-balhimycin

The conformation of the antibiotic ureido-balhimycin has been investigated by proton nmr and refined with computer simulations. The dimeric structure cj ureido-balhimycin was established by unambiguous identification of 6 nuclear Overhauser effects (NOEs) as intermonomeric, out of the total of 186 NOEs observed. Via distance geometry calculations the antiparallel orientation of the two monomers was demonstrated. Further refinement by molecular dynamics simulations provided the essential hydrogen bonds and aromatic interactions responsible for interfacial stabilization. The computational protocol illustrated here (distance geometry to define the coarse topology and molecular dynamics for refinement) should find general applicability in the study of homodimers. © 1995 John Wiley & Sons, Inc.     

ESR study of intercalation: quantitative evaluation of drug-DNA binding through competition with a spin-labeled 9-aminoacridine

Interaction of a spin-labeled 9-aminoacridine with DNA was studied by electron spin resonance spectroscopy. Accurate determination of the binding parameters, equilibrium dissociation constant (K_D), and total number of ligand-binding sites was obtained using Scatchard and Lineweaver-Burk plots. The competition between 9-aminoacridine and its spin-labeled derivative was examined by a similar analysis of the spin-label signals. The practical interest of this method lies in the fact that the precision and the simplicity of the measurements allow the quantitative determination of the binding capacity of any intercalative drug which interacts specifically with adenine/thymine bases of DNA.     

Conformation and dynamics of heparin and heparan sulfate

The glycosaminoglycans heparin and heparan sulfate contain similar structural units in varying proportions providing considerable diversity in sequence and biological function. Both compounds are alternating copolymers of glucosamine with both iduronate- and glucuronate-containing sequences bearing N-sulfate, N-acetyl, and O-sulfate substitution. Protein recognition of these structurally-diverse compounds depends upon substitution pattern, overall molecular shape, and on internal mobility. In this review particular attention is paid to the dynamic aspects of heparin/heparan sulfate conformation. The iduronate residue possesses an unusually flexible pyranose ring conformation. This extra source of internal mobility creates special problems in rationalization of experimental data for these compounds. We present herein the solution-state NMR parameters, fiber diffraction data, crystallographic data, and molecular modeling methods employed in the investigation of heparin and heparan sulfate. Heparin is a useful model compound for the sulfated, protein-binding regions of heparan sulfate. The literature contains a number of solution and solid-state studies of heparin oligo- and polysaccharides for both isolated heparin species and those bound to protein receptors. These studies indicate a diversity of iduronate ring conformations, but a limited range of glycosidic linkage geometries in the repeating disaccharides. In this sense, heparin exhibits a well-defined overall shape within which iduronate ring forms can freely interconvert. Recent work suggests that computational modeling could potentially identify heparin binding sites on protein surfaces.     

Conformation and dynamics in a Z-DNA tetramer

Two kinds of conformational variability have been reported for left-handed Z-DNA: the Z to Z′ transition, which involves a change in guanine sugar pucker from C-3′-endo to C-1′-exo, and the Z_Ito Z_II transition, which corresponds to a simple three-atom phosphate-group rotation. Neither of these motions substantially affects base stacking or helical twist, and this is because the degree of independent motion of phosphate groups, sugar molecules and base-pairs is greater in the left-handed Z helix than in right-handed B-DNA. Detailed considerations of Z helix geometry suggest that Z_I, Z_IIand Z′ are not separate species, but only samplings of the full range of conformation open to Z-DNA.     

Conformation and dynamics of an Fab'-bound peptide by isotope-edited NMR spectroscopy.

The dynamics and conformation of the peptide antigen MHKDFLEKIGGL bound to the Fab' fragment of the monoclonal antipeptide antibody B13A2, raised against a peptide from myohemerythrin, have been investigated by isotope-edited NMR techniques. The peptides were labeled with 15N (98%) or 13C (99%) at the backbone of individual amino acid residues. Well-resolved amide proton and nitrogen backbone resonances were obtained and assigned for eight of the 12 residues of this bound peptide. Significant resonance line width and chemical shift differences were observed. The 15N and 1H line width variations are attributed to differential backbone mobilities among the bound peptide residues which are consistent with the previously mapped epitope of this peptide antigen. Local structural information was obtained from isotope-directed NOE studies. The approximate distances associated with the experimental NOEs were estimated on the basis of a theoretical NOE analysis involving the relative integrated intensities of the NOE and source peaks. In this way, the sequential NH-NH NOEs obtained for seven of the Fab'-bound peptide residues were shown to correspond to interproton separations of approximately 3 A or less. Such short distances indicate that the backbone dihedral angles of these residues are in the alpha rather than the beta region of phi,psi conformational space; the peptide most likely adopts a helical conformation from F5 to G11 within the antibody combining site. The significance of these results with respect to the type and extent of conformational information obtainable from studies of high molecular weight systems is discussed.     

Atomic force microscope of drug-DNA interaction.

We have been investigated the possibility of B-Z transition in ZnTMPyP-DNA interaction based on the observation of spectroscopic data. In this study, we found drastic change in the AFM image of supercoiled plasmid DNA when it was interacted with TMPyPs indicating that the considerable amount of unwinding of double helix or B-Z transition is induced by the drug-DNA interaction. Such phenomena were not observed for other cationic drugs examined.     

Conformation and dynamics of an RNA internal loop

The conformation and the dynamics of an RNA oligonucleotide (26 nucleotides) which is a model for loop E in eukaryotic 5S RNA have been investigated by one- and two-dimensional NMR. The central portion of the oligonucleotide contains two G A oppositions, a common feature of ribosomal RNAs. The exchangeable proton spectrum indicates that an internal loop separates two stems of four and five base pairs. This observation is not consistent with structures for loop E containing mismatched G.A base pairs proposed from chemical and enzymatic studies on Xenopus laevis 5S RNA. The nonexchangeable proton spectrum has been assigned by two-dimensional NMR. Scalar couplings from correlated experiments and interproton distances from NOESY experiments at short mixing times have been used to determine glycosidic angles, sugar puckers, and other conformational features. The conformation of the stems is very close to standard A-form RNA, and extensive base stacking continues into the internal loop. This result provides a structural basis for the large favorable enthalpy of duplex formation determined in thermodynamic studies. Unusual structural and dynamic features are localized in the nucleotides connecting the loop to the stems.     

Molecular dynamics of structural transitions and intercalation in DNA

The large-scale flexibility of DNA and the intercalation of actinomycin D have been studied by computer simulation using molecular dynamics. The stretching and unwinding of B and Z forms of DNA and intercalation in B-DNA were examined through molecular dynamics simulations, and the energetics of transitions were calculated by the conformational energy minimization method. The principal results of this research are as follows: (1) A dynamic conformational pathway is presented for longitudinal stretching and unwinding of the double helix to open an intercalation site. (2) Large-scale transitions are possible in both B and Z forms of DNA through a conformationally allowed kinetic pathway. (3) The stretching and untwisting of a 5′(CG)3′ step is energetically more favorable than for a GC step in B-DNA. (4) The formation of an adjacent second cavity in B-DNA requires larger energy than the formation of the first cavity, affirming the neighbor-exclusion principle of intercalation. (5) Docking an intercalated actinomycin D in the stretched structure is shown to be geometrically and energetically feasible.     

Conformation and fluctuations of free stefin B: a molecular dynamics study.

Molecular dynamics study was performed on the cysteine proteinase inhibitor stefin B. Structure of inhibitor from the complex with papain was used as a starting point. Amino terminal "trunk" of the inhibitor which lies extended along the cleft of the enzyme in the complex, folded onto the body of inhibitor during MD simulation, thereby reducing the total and particularly hydrophobic surface exposed to the solvent. This effect counterbalances hydrophobic contribution of the "trunk" and explains why its deletion in stefin B and related inhibitors doesn't reduce the dissociation constant. The rest of stefin B conformation is conserved together with main chain hydrogen bonds. Fluctuations of C alpha atoms resembles crystallographic B factors with exception of residues in contact with enzyme.     

Conformation and aggregation of M13 coat protein studied by molecular dynamics.

Molecular dynamics (MD) simulations are performed on M13 coat protein, a small membrane protein for which both alpha- and beta-structures have been suggested. The simulations are started from initial conformations that are either monomers or dimers of alpha-helices or U-shaped beta-sheets. The lipid bilayer is represented by a hydrophobic potential. The results are analyzed in terms of stability, energy and secondary structure. The U-shaped beta-structure changes from a planar to a twisted form with larger twist for the monomer than the dimer. The beta-sheet is much more flexible than the alpha-helix as monitored by the root mean square (rms) fluctuations of the C alpha atoms. A comparison of the energies after 100 ps MD simulation shows that of the monomers, the alpha-helix has the lowest energy. The energy difference between alpha- and beta-structures decreases from 266 kJ/mol to 148 kJ/mol, when going from monomers to dimers. It is expected that this difference will decrease with higher aggregation numbers.     

Thermodynamics of drug-DNA interactions

Many anticancer, antibiotic, and antiviral drugs exert their primary biological effects by reversibly interacting with nucleic acids. Therefore, these biomolecules represent a major target in drug development strategies designed to produce next generation therapeutics for diseases such as cancer. In order to improve the clinical efficacy of existing drugs and also to design new ones it is necessary to understand the molecular basis of drug-DNA interactions in structural, thermodynamic, and kinetic detail. The past decade has witnessed an increase in the number of rigorous biophysical studies of drug-DNA systems and considerable knowledge has been gained in the energetics of these binding reactions. This is, in part, due to the increased availability of high-sensitivity calorimetric techniques, which have allowed the thermodynamics of drug-DNA interactions to be probed directly and accurately. The focus of this article is to review thermodynamic approaches to examining drug-DNA recognition. Specifically, an overview of a recently developed method of analysis that dissects the binding free energy of these reactions into five component terms is presented. The results of applying this analysis to the DNA binding interactions of both minor groove drugs and intercalators are discussed. The solvent water plays a key role in nucleic acid structure and consequently in the binding of ligands to these biomolecules. Any rational approach to DNA-targeted drug design requires an understanding of how water participates in recognition and binding events. Recent studies examining hydration changes that accompany DNA binding by intercalators will be reviewed. Finally some aspects of cooperativity in drug-DNA interactions are described and the importance of considering cooperative effects when examining these reactions is highlighted.     

Conformation of secretin in dimethyl sulfoxide solution. NMR studies and restrained molecular dynamics

The solution conformation of the 27-residue polypeptide hormone secretin in dimethyl sulfoxide has been determined on the basis of 1H-NMR measurements. The experimental data set used in the structure determination consisted of 98 nuclear-Overhauser-enhancement-derived interproton and dihedral angle restraints from coupling constants. The NH-NH and H alpha-NH NOEs were determined from build-up rates, while the remaining distances were classified in a qualitative manner. The structure calculations consisted of two phases. First, dynamical simulated annealing calculations were carried out to find conformations of the peptide which satisfy NOE and phi dihedral restraints. The convergence of ten calculated structures was good except for those regions of the molecule where NOE data were not unambiguous. From the calculated set another initial structure was built which was again minimized in several 5-ps calculations now employing the full empirical energy function. The resulting structures of secretin reveal conformationally well-defined regions, but not a single uniform secondary structure. The structure is different from the calculated structure from trifluoroethanol/water measurements.     

Conformation and dynamics of the chemotactic peptide formyl-L-methionyl-L-leucyl-L-phenylalanine in solution

Several conformational and dynamic features of the chemotactic peptide formyl-L-methionyl-L-leucyl-L-phenylalanine in solution have been delineated by investigations of NMR and IR spectroscopic parameters. Both 1D and 2D NMR experiments have been performed for detection of scalar and dipolar proton-proton connectivities, whereas 13C and 1H relaxation parameters have been interpreted in terms of molecular dynamics. The main conformation appeared to be unfolded with the three hydrophobic side chains extending in divergent directions with respect to the backbone. The existence of relatively weak intermolecular hydrogen bonds was demonstrated, involving the formamide end group, with increase in the hydrophobicity of the external surface.     

Sequence specificity of drug-DNA interactions

Methods for determining sequence specificities of anticancer drugs, carcinogens, and mutagens which interact with natural DNA's are presented. For drugs which nick or covalently bind to DNA and thus leave a permanent record of their residence position on the helix, the sequences important in drug action can be readily determined. For agents which interact with DNA in an equilibrium fashion, "footprinting" analysis, a technique used to investigate protein-DNA binding, has proved to be useful in studying drug-DNA interactions. The sequence specificities of a number of small ligands which interact with natural DNA's are also presented.     

Conformation and dynamics of bovine brain S-100a protein determined by fluorescence spectroscopy.

We have used time-resolved laser fluorescence spectroscopy to investigate the intensity and anisotropy decays of the single tryptophan residue in bovine brain S-100a (alpha beta) protein. The steady-state and acrylamide quenching results indicated that the Trp 90 of the alpha-subunit was partially buried in a relatively nonpolar environment at pH 7.5. Both Ca2+ and pH 8.5 slightly enhanced the exposure of the residue to the solvent, but the residue remained partially buried in the calcium complex at both pH values. The best representation of the intensity decays was a linear combination of three exponential terms, regardless of solvent condition and temperature. The three lifetimes (tau i) were in the range of 0.4-5 ns and insensitive to emission wavelength, but their fractional amplitudes (alpha i) shifted in favor of the shortest component (alpha 1) when the decays were measured at the blue end of the emission spectrum. These results suggest that an excited-state interaction between the indole ring and the side chain of an adjacent residue may be responsible for the observed shortest lifetime. In the presence of Ca2+, the three lifetimes remained relatively unaltered, but the values of alpha 1 decreased by a factor of 2.3 at pH 7.2 and a factor of 1.8 at pH 8.2. This Ca(2+)-induced decrease may be attributed to disruption of the putative excited-state interaction resulting from reorientations of the alpha-helical segments flanking a Ca(2+)-binding loop (residues 62-73). At both pH 7.2 and 8.4, the anisotropy decays of the apoprotein followed a biexponential decay law.(ABSTRACT TRUNCATED AT 250 WORDS)     

Conformation of a neurokinin antagonist in solution. 2D NMR and restrained molecular dynamics study.

The hexapeptide [cyclo(Leu1 psi(CH2NH2)Leu2-Gln3-Trp4-Phe5-Gly6)]+1 is a potent antagonist of neurokinin A activity in tissues of hamster urinary bladder. The solution conformation of this cyclic hexapeptide has been characterized by the combined use of two dimensional nuclear magnetic resonance spectroscopy and restrained molecular dynamics. The proton spectrum of the peptide was fully assigned by the sequential assignment procedure. Interproton distances were derived from crosspeak volumes in two dimensional Nuclear Overhauser Effect spectra, and dihedral angles were calculated from appropriate coupling constants. Temperature coefficients of the amide protons were determined. Restrained molecular dynamics simulations were carried out using the backbone interproton distances as constraints. During 210 ps of restrained molecular dynamics the peptide interconverted among three closely related families of conformations. These interconversions occurred at picosecond timescales under the simulation conditions.     

Conformation and molecular dynamics calculations on uteroglobin fragment 18-47

The conformational properties of fragment 18–47 of rabbit uteroglobin in aqueous solution containing SDS micelles were investigated by two-dimensional nmr spectroscopy and molecular dynamics calculations. The fragment comprises helices II and III and the β-turn connecting the two helices. The nmr results and nmr-restrained molecular dynamics calculations showed that in the isolated fragment the elements of secondary structure present in the intact protein are preserved only in part. Specifically, a well-defined α-helix was found in the sequence 33–44, corresponding to helix III of uteroglobin, while the regions of helix II and β-turn are characterized by high flexibility in the fragment. © 1994 John Wiley & Sons, Inc.     

Conformation of viroids.

Viroids are uncoated infectious RNA molecules (MW 107 000-127 000) known as pathogens of certain higher plants. Thermodynamic and kinetic studies were carried out on highly purified viroid preparations by applying UV-absorption melting analysis and temperature jump methods. The thermal denaturation of viroids is characterized by high thermal stability, high cooperativity and a high degree of base pairing. Two relaxation processes could be resolved; a process in the sec range could be evaluated as an independent all-or-none-transition with the following properties: reaction enthalpy= 550 kcal/mol, activation enthalpy of the dissociation = 470 kcal/mol; G : C content = 72 %. These data indicate the existence of an uninterrupted double helix of 52 base pairs. A process in the msec range involves 15 - 25 base pairs which are most probably distributed over several short double helical stretches. A tentative model for the secondary structure of viroids isproposed and the possible functional implications of their physicochemical properties are discussed.