Molecular mechanical studies of proflavine and acridine orange intercalation

Differential staining of peripheral nuclear chromatin with Acridine orange implies an A-form epichromatin conformation of the DNA

The chromatin observed by conventional electron microscopy under the nuclear envelope constitutes a single layer of dense 30–35 nm granules, while ∼30 nm fibrils laterally attached to them, form large patches of lamin-associated domains (LADs). This particular surface “epichromatin” can be discerned by specific (H2A+H2B+DNA) conformational antibody at the inner nuclear envelope and around mitotic chromosomes. In order to differentiate the DNA conformation of the peripheral chromatin we applied an Acridine orange (AO) DNA structural test involving RNAse treatment and the addition of AO after acid pre-treatment. MCF-7 cells treated in this way revealed yellow/red patches of LADs attached to a thin green nuclear rim and with mitotic chromosomes outlined in green, topologically corresponding to epichromatin epitope staining by immunofluorescence. Differentially from LADs, the epichromatin was unable to provide metachromatic staining by AO, unless thermally denatured at 94^oC. DNA enrichment in GC stretches has been recently reported for immunoprecipitated ∼ 1Kb epichromatin domains. Together these data suggest that certain epichromatin segments assume the relatively hydrophobic DNA A-conformation at the nuclear envelope and surface of mitotic chromosomes, preventing AO side dimerisation.  We hypothesize that epichromatin domains form nucleosome superbeads. Hydrophobic interactions stack these superbeads and align them at the nuclear envelope, while repulsing the hydrophilic LADs. The hydrophobicity of epichromatin explains its location at the surface of mitotic chromosomes and its function in mediating chromosome attachment to the restituting nuclear envelope during telophase.

Structure and dynamics of proflavine association around DNA

Proflavine is a small molecule that intercalates into DNA and, thereby, acts as an anticancer agent. Intercalation of proflavine is shown to be a two-step process in which the first step is believed to be the formation of a pre-intercalative outside bound state. Experimental studies so far have been unable to capture the nature of the outside bound state. However, the sub-millisecond timescale observed in fluorescence kinetic experiments is often attributed to the binding of proflavine outside of DNA. Here, we have performed molecular dynamics simulations with multiple proflavine molecules to study the structure and dynamics of the formation of the outside bound state of DNA at different ion concentrations. We observed that the timescale of the outside bound state formation is, at least, five orders of magnitude faster (in nanoseconds) than the experimentally reported timescale (sub-milliseconds) attributed to binding outside DNA. Moreover, we also observed the stacked arrangement of proflavine all around DNA, which is different from the experimentally predicted stacking arrangement perpendicular to the helical axis of DNA in the close vicinity of the phosphate groups. This study, therefore, provides insight into the molecular structure and dynamics of the pre-intercalative outside bound state and will help in understanding the overall intercalation mechanism.

Drug-DNA intercalation: from discovery to the molecular mechanism

The ability of small molecules to perturb the natural structure and dynamics of nucleic acids is intriguing and has potential applications in cancer therapeutics. Intercalation is a special binding mode where the planar aromatic moiety of a small molecule is inserted between a pair of base pairs, causing structural changes in the DNA and leading to its functional arrest. Enormous progress has been made to understand the nature of the intercalation process since its idealistic conception five decades ago. However, the biological functions were detected even earlier. In this review, we focus mainly on the acridine and anthracycline types of drugs and provide a brief overview of the development in the field through various experimental methods that led to our present understanding of the subject. Subsequently, we discuss the molecular mechanism of the intercalation process, free-energy landscapes, and kinetics that was revealed recently through detailed and rigorous computational studies.

Solution structure and thermodynamics of 2',5' RNA intercalation

As a means to explore the influence of the nucleic acid backbone on the intercalative binding of ligands to DNA and RNA, we have determined the solution structure of a proflavine-bound 2',5'-linked octamer duplex with the sequence GCCGCGGC. This structure represents the first NMR structure of an intercalated RNA duplex, of either backbone structural isomer. By comparison with X-ray crystal structures, we have identified similarities and differences between intercalated 3',5' and 2',5'-linked RNA duplexes. First, the two forms of RNA have different sugar pucker geometries at the intercalated nucleotide steps, yet have the same interphosphate distances. Second, as in intercalated 3',5' RNA, the phosphate backbone angle zeta at the 2',5' RNA intercalation site prefers to be in the trans conformation, whereas unintercalated 2',5' and 3',5' RNA prefer the -gauche conformation. These observations provide new insights regarding the transitions required for intercalation of a phosphodiester-ribose backbone and suggest a possible contribution of the backbone to the origin of the nearest-neighbor exclusion principle. Thermodynamic studies presented for intercalation of both structural RNA isomers also reveal a surprising sensitivity of intercalator binding enthalpy and entropy to the details of RNA backbone structure.

Intercalation complex of proflavine with DNA: structure and dynamics by solid-state NMR

The structure of the complex formed between the intercalating agent proflavine and fibrous native DNA was studied by one- and two-dimensional high-resolution solid-state nuclear magnetic resonance (NMR). Carbon-13-labeled proflavine was used to show that the drug is stacked with the aromatic ring plane perpendicular to the fiber axis and that it is essentially immobile. Natural abundance carbon-13 NMR of the DNA itself shows that proflavine binding does not change the puckering of the deoxyribose ring. However, phosphorus-31 NMR spectra show profound changes in the orientation of the phosphodiester grouping on proflavine binding, with some of the phosphodiesters tilting almost parallel to the helix axis, and a second set almost perpendicular. The first group to the phosphodiesters probably spans the intercalation sites, whereas the tilting of the second set likely compensates for the unwinding of the DNA by the intercalator.

Free energy calculation on base specificity of drug--DNA interactions: application to daunomycin and acridine intercalation into DNA

We present the results of free energy perturbation/molecular dynamics studies on B-DNA.daunomycin and B-DNA.9-aminoacridine complexes as well as on B-DNA itself in order to calculate the free energy differences between complexes having different base pair sequences. The results generally reproduce the trends observed experimentally, i.e., preferences of acridine and daunomycin to bind to a specific base sequence in the DNA. This is encouraging, given the simplicity of the molecular mechanical/dynamical model in which solvent is not explicitly included.

Strategies for the use of computational SAR methods in assessing genotoxicity

The relationship between computational SAR studies and relevant data gathering and generation activities is complex. First, the chemical class to be studied is selected on the basis of information requirements for hazard identification and assessment. Membership in the class is determined by consideration of chemical structure and reactivity. Compilation of the existing bioassay data for this chemical class follows immediately from the specification of the class. Bioassay data, qualitative knowledge of general chemical reactivities in this class, and knowledge concerning potential interactions with biomolecular targets all contribute to the derivation of possible mechanisms for biological activity. Computational studies based on modeling the proposed mechanism of action and/or the existing data base can provide a quantitative basis for the differentiation between chemicals. There is the opportunity for continuing feedback between the quantitative computational studies and the development of a relevant bioassay data base for this chemical class. The qualitative and quantitative information on the potential biological responses obtained will provide a rational basis for extrapolation from the extant data base to the chemicals of interest, and to biological responses significant to the assessment for which complete data are unavailable. Knowledge concerning possible mechanisms of action and preexisting data determine the type of computational study that will be most useful.

Harmonic dynamics of a DNA hexamer in the absence and presence of the intercalator ethidium

Vibrational normal mode calculations are presented for a DNA hexanucleoside pentaphosphate, d(CpGpCpGpCpG)2, and for its complex with the cationic intercalator ethidium. Two intercalation sites are modeled that differ in DNA backbone torsion angles. Normal mode frequencies for the DNA fragment itself are significantly lower than those reported earlier using different force fields, but an analysis of "effective" frequencies suggests that somewhat higher frequencies are more appropriate. Intercalation leads to significant lowering of mobility for the base pairs adjacent to the drug; in this sequence, the ethidium binding affects the guanosine atoms more strongly than the cytosine atoms. Motions of the bases and the intercalator are analyzed in terms of "twist" about the local helix axis and a "tilt" angle relative to this axis, and the results are compared to fluorescence studies of ethidium-DNA complexes.

Computer-aided molecular modelling: research study or research tool?

Developments in computational methods and equipment have produced a new type of research chemist, who prefers to calculate properties as well as measure them, either to gain a better understanding of microscopic molecular behaviour per se, or to guide a broader scientific study using a so-called 'rational' approach. While there is good reason to believe that significant results can be obtained this way, it is clear that only some of the 'tools of the trade' are sufficiently robust to present to those who are not experts in the field. This paper discusses the underlying basis for molecular modelling techniques, describing their history, state of maturity and prospects for future development. The intention is to outline the scope that these play in an industrial research environment, and to examine how they can successfully be incorporated as routine research tools. Computer-based molecular modelling is not a unique discipline. Some comparisons are drawn between it and other computer-based techniques that have reached a greater degree of maturity, in order to highlight the points made.

Conformations of duplex structures formed by oligodeoxynucleotides covalently linked to the intercalator 2-methoxy-6-chloro-9-aminoacridine

A family of covalent complexes between oligonucleotides and derivatives of the intercalating agent 9-amino acridine has been synthesized (Asseline, U., Thuong, N.T. and Helene, C. (1983) C.R.Acad. Sci. (Paris) 297 (III), 369-372) and studied (Lancelot, G., Asseline, U., Thuong, N.T., and Helene, C. (1985) Biochemistry 24, 2521-2529; Lancelot, G., Asseline, U., Thuong, N.T., and Helene, C. (1985) J. Biomol. Str. Dyn. 3, 913-921) with a view to understand nucleic acid-nucleic acid recognition. In order to understand the nature of interactions between the intercalator and the oligonucleotides in such complexes and the sensitivity of such interactions to the polymorphic form of the DNA, we have carried out molecular mechanics simulations on duplex deoxyoligonucleotides d(A)6.d(T)6 (A and B forms) and d(TATC).d(GATA) (B form) covalently bound to 2-methoxy-6-chloro-9-aminoacridine through a pentamethylene linker chain. Structures in which the acridine derivative is end stacked (at the 3' and 5' ends) and in which the dye is intercalated between the terminal base pairs (at both the ends) and between second and third base pairs from the 3' end are all of reasonably low energy in both A and B forms of DNA. Our studies on 3' end complexes find that in the B form, intercalation of the dye between the second and third base pairs is preferred over the other two modes of binding, while in the A form, intercalation between the terminal base pairs is preferred. In the 5' end A and B form complexes, outside stacking and intercalation between the terminal base pairs are preferred, respectively. Our calculations suggest the possibility that the presence of the dye attached covalently to the DNA can induce conformational transitions in the DNA. For example, intercalation of the dye two base pairs from the end could induce an A----B transition.

Recognition and interactions controlling the assemblies of beta barrel domains

We present a qualitative computer graphics approach to the characterization of forces important to the assembly of beta domains that should have general utility for examining protein interactions and assembly. In our approach, the nature of the molecular surface buried by the domain contacts, the specificity of the residue-to-residue interactions, and the identity of electrostatic, hydrophobic, and hydrophilic interactions are elucidated. These techniques are applied to the beta barrel domains of Cu, Zn superoxide dismutase (SOD), immunoglobulin Fab, and tomato bushy stunt virus coat protein (TBSV), a plant viral capsid protein. By looking at a set of proteins having different numbers of interacting beta domains, we have been able to see some of the variety and also some of the patterns common to these assembled domains. Strong beta domain interactions (identified by their biochemical integrity) are apparently due to chemical, electrostatic, and shape complementarity of the molecular surfaces buried from interaction with solvent molecules. Although the amount of hydrophobic buried surface area appears to correlate with the strength of the interaction, electrostatic forces appear to be important in both stabilizing and destabilizing specific contacts. In TBSV, analysis of electrostatic interactions may help explain mechanisms of subunit accommodation to different environments, particle expansion, and pathways of assembly. The possible molecular basis for observed differences in the stability and flexibility of the domain complexes is discussed.

Molecular electrostatic potentials: an effective tool for the elucidation of biochemical phenomena

The electrostatic potential V(r) that is created in the space around a molecule by its nuclei and electrons (treated as static distributions of charge) is a very useful property for analyzing and predicting molecular reactive behavior. It is rigorously defined and can be determined experimentally as well as computationally. The potential has been particularly useful as an indicator of the sites or regions of a molecule to which an approaching electrophile is initially attracted, and it has also been applied successfully to the study of interactions that involve a certain optimum relative orientation of the reactants, such as between a drug and its cellular receptor. A variety of methods for calculating V(r) is available, at different levels of rigor. For large biologically active molecules, multipole expansions and superposition of potentials computed for subunits have been found to be effective. A large number of chemical and biochemical systems and processes have now been studied in terms of electrostatic potentials. Three examples of such applications are surveyed in this paper. These deal with: (a) reactive properties of nucleic acids, including their component bases; (b) biological recognition processes, including drug-receptors and enzyme-substrate interactions; and (c) chemical carcinogenesis, referring specifically to the polycyclic aromatic hydrocarbons and halogenated olefins and their epoxides. For each of these areas, examples of the use of electrostatic potentials in elucidating structure-activity patterns are given.

Nucleic acid binding drugs. Part XIII. Molecular motion in a drug-nucleic acid model system: thermal motion analysis of a proflavine-dinucleoside crystal structure

The high-resolution crystal structure of the intercalation complex between proflavine and cytidylyl-3',5'-guanosine (CpG) has been studied by thermalmotion analysis. This has provided information on the translational and librational motions of individual groups in the complex. Many of these motions are similar to, though of larger magnitude than in uncomplexed dinucleosides. Pronounced librational effects were observed along the base pairs and in the plane of the drug chromophore.

500 MHz 1H-NMR study of the interaction of daunomycin with B and Z helices of d(CGm5CGCG)

The interaction of daunomycin with B and Z helices of a self-complementary DNA fragment d(CGm5CGCG) in solution was studied by 1H-NMR spectroscopy at 500 MHz. The results show that the B-Z transition kinetics is not affected by addition of daunomycin. Daunomycin binds exclusively to the B form of d(CGm5CGCG). Z exchanges with B while the latter also exchanges with the B duplex-daunomycin complexes.

1H-NMR study of the interaction of daunomycin with B-DNA helices of methylated oligodeoxynucleotides

The interaction of daunomycin with B-DNA double helices of several methylated deoxynucleotides, d(C-G-m5C-G), d(m5C-G-C-G), d(C-G-m5C-G-C-G) and d(m5C-G-C-G-m5C-G) in solution was investigated by 1H-NMR spectroscopy at 500 MHz. At low temperature (t less than 20 degrees C for the tetramer and t less than 40 degrees C for the hexamers), several daunomycin-DNA complexes were observed in slow exchange with the drug-free DNA duplexes. The presence of daunomycin in a self-complementary double helix cancels the conformational symmetry of the two strands; the proton signals can split into several others owing to the difference between free and intercalated duplexes and to the many possible intercalation sites in a duplex (three for a tetramer, five for an hexamer). A model relating the chemical shifts of splitted proton signals to the various intercalated duplex conformations was given. The results show that one daunomycin molecule is associated with one duplex and that it can enter any intercalation site with equal probability; no side-effects were observed even for very short helices (of a tetramer). In the case of d(C-G-m5C-G) the association constant and the dissociation and association rates of the intercalated complex were evaluated.

Structural and sequence-dependent aspects of drug intercalation into nucleic acids

Information gained from X-ray crystallographic studies on drug-nucleic acid complexes is described, with emphasis on the intercalation process. Relevant data from NMR experiments are examined in order to highlight similarities and differences between solution and solid-state structures. Theoretical analyses of intercalation complexes are also discussed and evaluated, with respect to the structural methods, with special reference being made to nucleic acid conformation and positions of drug molecules in the binding sites.

Electrostatic recognition between superoxide and copper, zinc superoxide dismutase

Electrostatic forces have been implicated in a variety of biologically important molecular interactions including drug orientation by DNA, protein folding and assembly, substrate binding and catalysis and macromolecular complementarity with inhibitors, drugs and hormones. To examine enzyme-substrate interactions in copper, zinc superoxide dismutase (SOD), we developed a method for the visualization and analysis of an enzyme's three-dimensional electrostatic vector field that allows the contributions of specific residues to be identified. We report here that the arrangement of electrostatic charges in SOD promotes productive enzyme-substrate interaction through substrate guidance and charge complementarity: sequence-conserved residues create an extensive electrostatic field that directs the negatively charged superoxide (O-2) substrate to the highly positive catalytic binding site at the bottom of the active-site channel. Dissection of the electrostatic potential gradient indicated the relative contributions of individual charged residues: Lys 134 and Glu 131 seem to have important roles in directing the long-range approach of O-2, while Arg 141 has local orienting effects. The reported methods of analysis may have general application for the elucidation of intermolecular recognition processes.

Electrostatic potential molecular surfaces

Color-coded computer graphics representations of the electrostatic potentials of trypsin, trypsin-inhibitor, prealbumin and its thyroxine complex, fragments of double-helical DNA, and a netropsin--DNA complex illustrate the electrostatic and topographic complementarity in macromolecule-ligand interactions. This approach is powerful in revealing intermolecular specificity and shows promise of having predictive value in drug design.

A molecular mechanical study of complexes formed between 4-nitroquinoline-N-oxide and dinucleoside phosphates

Molecular mechanical calculations were done on complexes of 4-nitroquinoline-N-oxide (NQO) with various dinucleoside phosphates [(ApT)2, (CpG)2, (GpC)2, and (TpA)2]. Models built using proflavine (uniform C3' endo sugar puckers) and acridine orange (mixed C3' endo (3'-5') C2' endo sugar puckers) dinucleoside phosphate X-ray structures were used in the calculations. Relative binding energies, complex geometries, and various intercalator orientations in the complexes were studied. The results suggest qualitatively different geometries for pyr-(3'-5')-pur and pur-(3'-5')-pyr sequences. Specifically, we find marked distortion in some of the complexes (i.e. there is not a parallel coplanar relationship between the base pairs and intercalator), distortion of the NQO nitro group from planarity in the complexes and mobility of NQO in the intercalation site. We suggest that experimental studies of NQO-dinucleoside phosphate complexes may reveal intercalation complexes which deviate substantially more from a nearly parallel coplanar arrangement of bases and intercalator than has been previously observed.