Chromatin reconstitution on small DNA rings. II. DNA supercoiling on the nucleosome.

DNA supercoiling on the nucleosome was investigated by relaxing with topoisomerase I mono- and dinucleosomes reconstituted on small DNA rings. Besides 359 base-pair (bp) rings whose linking differences were integers, two additional series of rings with fractional differences, 341 and 354 bp in size, were used. Mononucleosomes reconstituted on 359 bp rings were found to relax into a single mononucleosome form. In contrast, 341 and 354 bp mononucleosomes relaxed into a mixture of two forms, corresponding to two adjacent topoisomers. The observation that the ratio between these two forms was, within each ring series, virtually independent of the initial linking number of the topoisomer used for the reconstitution suggested that each partition reflected an equilibrium. Comparison with the equilibria observed for the same rings in the absence of histones showed that the formation of a single nucleosome is associated with a linking number change of -1.1(+/-0.1) turn. Dinucleosomes, in contrast, were not relaxed to completion and do not reach equilibria. The corresponding linking number change per nucleosome was, however, estimated to be similar to the above figure, in agreement with previous data from the literature obtained with circular chromatins containing larger numbers of nucleosomes. DNA structure in mononucleosomes was subsequently investigated by means of high-resolution electron microscopy and gel electrophoresis. It was found that the above linking number reduction could be ascribed to a particle with a large open extranucleosomal DNA loop and with no more than 1.5 turns of a superhelix around the histone core. A theoretical model of a nucleosome on a small ring was constructed in which one part of the DNA was wrapped around a cylinder and the other part was free to vary both in torsion and flexion. The linking number reduction predicted was found to be most consistent with experimental data when the twist of the DNA in the superhelix was between 10.5 and 10.65 pb per turn, suggesting that wrapping on the nucleosome does not alter the twist of the DNA significantly. A lower estimate of the linking number reduction associated with a two-turn nucleosome was also derived, based on an analysis of recent data obtained upon treatment of reconstituted minichromosomes with gyrase. The value, 1.6 turns, set a lower limit of 10.44 bp per turn for the twist of nucleosomal DNA, in agreement with the above estimate.(ABSTRACT TRUNCATED AT 400 WORDS)

Solution structure of an oncogenic DNA duplex, the K-ras gene and the sequence containing a central C.A or A.G mismatch as a function of pH: nuclear magnetic resonance and molecular dynamics studies.

The DNA duplex 5' d(GCCACCAGCTC)-d(GAGCTGGTGGC) corresponds to the sequence 29 to 39 of the K-ras gene, which contains a hot spot for mutations. This has been studied by one and two-dimensional nuclear magnetic resonance, energy minimization and molecular dynamics. The results show that it adopts a globally B-DNA type structure. We have introduced, at the central base-pair, the mismatches C.A and A.G. The mismatch position is that of the first base of the Gly12 codon, the hot spot. For the C.A mismatch we observe a structural change as a function of pH with an apparent pKa of 7.2. At low pH, the mismatch pair adopts a structure close to a classic wobble conformation with the cytidine residue displaced into the major groove. It is stabilised by two hydrogen bonds in which the adenosine residue is protonated and the cytidine residue has a significant C3'-endo population. At high pH, the mispair structure is in equilibrium between wobble and reverse wobble conformations. Similar studies are reported on the A.G mismatch, which also undergoes a transition as a function of pH. 31P spectra have been recorded on all systems and as a function of pH. No evidence for BII phosphodiester backbone conformations was found. The NMR results are well corroborated by molecular dynamics calculations performed with or without distance constraints. The dynamics at the mismatch sites have been examined. Although the overall structures are close to B-DNA, helical parameters fluctuate differently at these sites. Different hydrogen bonding alternatives in dynamic equilibrium that can involve three-centred hydrogen bonds are observed.

About the large fluctuations observed using gas-phase molecular dynamics in the K-ras gene containing a mismatch.

We describe how we can reduce the periodic bending motions in the simulation in vacuo of the molecular dynamics of a short DNA fragment containing the Gly 12 hot spot of the K-ras oncogene and having at its center a mismatch CA+.

Structure-activity relationships and binding mode of styrylquinolines as potent inhibitors of HIV-1 integrase and replication of HIV-1 in cell culture.

Our prior studies showed that polyhydroxylated styrylquinolines are potent HIV-1 integrase (IN) inhibitors that block the replication of HIV-1 in cell culture at nontoxic concentrations. To explore the mechanism of action of these inhibitors, various novel styrylquinoline derivatives were synthesized and tested against HIV-1 IN and in cell-based assays. Regarding the in vitro experiments, the structural requirements for biological activity are a carboxyl group at C-7, a hydroxyl group at C-8 in the quinoline subunit, and an ancillary phenyl ring. However the in vitro inhibitory profile tolerates deep alterations of this ring, e.g. by the introduction of various substituents or its replacement by heteroatomic nuclei. Regarding the ex vivo assays, the structural requirements for activity are more stringent than for in vitro inhibition. Thus, in addition to an o-hydroxy acid group in the quinoline, the presence of one ortho pair of substituents at C-3' and C-4', particularly two hydroxyl groups, in the ancillary phenyl ring is imperatively required for inhibitory potency. Starting from literature data and the SARs developed in this work, a putative binding mode of styrylquinoline inhibitors to HIV-1 IN was derived.

Modeling of the inhibition of retroviral integrases by styrylquinoline derivatives.

Styrylquinoline derivatives, known to be potent inhibitors of HIV-1 integrase, have been experimentally tested for their inhibitory effect on the disintegration reaction catalyzed by catalytic cores of HIV-1 and Rous sarcoma virus (RSV) integrases. A modified docking protocol, consisting of coupling a grid search method with full energy minimization, has been specially designed to study the interaction between the inhibitors and the integrases. The inhibitors consist of two moieties that have hydroxyl and/or carboxyl substituents: the first moiety is either benzene, phenol, catechol, resorcinol, or salicycilic acid; the hydroxyl substituents on the second (quinoline) moiety may be in the keto or in the enol forms. Several tautomeric forms of the drugs have been docked to the crystallographic structure of the RSV catalytic core. The computed binding energy of the keto forms correlates best with the measured inhibitory effect. The docking procedure shows that the inhibitors bind closely to the crystallographic catalytic Mg(2+) dication. Additional quantum chemistry computations show that there is no direct correlation between the binding energy of the drugs with the Mg(2+) dication and their in vitro inhibitory effect. The designed method is a leading way for identification of potent integrase inhibitors using in silico experiments.

Tautomers of styrylquinoline derivatives containing a methoxy substituent: computation of their population in aqueous solution and their interaction with RSV integrase catalytic core.

8-Hydroxy-2-[2-(3-hydroxy-4-methoxyphenyl)ethenyl]-7-quinoline carboxylic acid and 8-hydroxy-2-[2-(3-methoxy-4-hydroxyphenyl)ethenyl]-7-quinoline carboxylic acid inhibit the processing and strand transfer reactions catalyzed by HIV-1 integrase with an IC50 of 2 microM. Some of their spectral properties are briefly reported. Their fluorescence is so weak that it is of no use in an experimental determination of the binding to the protein and we resorted to computer simulation. Both styrylquinoline derivatives, in their monoanionic form, have several dozens of tautomers and each of these forms has four planar rotamers. In this work computer simulations have been performed to determine which tautomer is the most abundant in aqueous solution and which binds to the Rous sarcoma virus (RSV) integrase catalytic core. As the substituents on the quinoline moiety are the same as on salicylic acid, the energies of hydroxy benzoic acid tautomers were also computed both in vacuo and embedded in a continuous medium which had the dielectric constant of bulk water, using the recent CPCM technique. The CPCM method was then applied to the two integrase inhibitors to estimate the tautomer population in water. The binding site of the compounds on the RSV integrase catalytic core was determined through a docking protocol, consisting of coupling a grid search method with full energy minimization. The designed method is a way leading to identification of potent integrase inhibitors using in silico experiments.

A database of 32 DNA triplets to study triple helices by molecular mechanics and dynamics.

We present here a database of 32 deoxyribonucleotide triplets, that can be used as building blocks of triple helix forming deoxyribonucleotides on a computer. This database is made of all the pairing schemes of the triplets ATT, GCC+, ATA and GCG where the third base forms two hydrogen bonds with the purine of the first two Watson-Crick strands. The essential features of the known triple helices were preserved in the resulting structures. A triple helix can be easily built from any combination of these basic triplets. Four homogeneous and alternate triple helices thus obtained were studied by molecular mechanics and dynamics in vacuo. The results are in agreement with known experimental observations for ATT and suggest a possible structure for the GCG triple helix. In order to characterize the geometry of the structures obtained, the definitions of nucleic acid structure parameters (R.E. Dickerson et al., EMBO J. 8 (1989) 1-4) have been extended to triple helical polynucleotides.

Counter-ion condensation and system dimensionality.

Condensation of the counter-ions around a highly charged infinitely long cylindrical molecule, such as DNA, can be described in detail in terms of the solutions of the Poisson-Boltzmann (Gouy-Chapman) equation. By using the Alfrey-Berg-Morawetz (1951) solution of this equation one can show that a certain fraction of the counter-ions remain within finite distances of the poly-ion even when the volume of the system is expanded indefinitely; these ions can be appropriately called "condensed". The fraction of the macromolecule's charge represented by these ions is just 1-1/xi, where xi is the linear charge-density parameter of the macromolecule; this is also the value given by Manning's theory. The question arises: Is this property unique to the infinite cylinder? Using the same PB equation, we can consider the infinite charged plane and a large finite charged sphere for comparison. In the case of the plane all of the counter-ions are condensed in the above sense, not just a fraction, for any surface charge density of the plane. These ions form the classical Gouy double layer. On the other hand, none of the counter-ions of the charged sphere are condensed in the above sense, no matter how high the surface charge density. Thus the cylinder is a unique intermediate case in which a fraction of the counter-ions are condensed if the linear charge density is higher than the critical value of unity.

Molecular mechanics and dynamics of an abasic frameshift in DNA and comparison to NMR data.

In a previous publication (Ph. Cuniasse, L.C. Sowers, R. Eritja, B. Kaplan, M.F. Goodman, J.A.H. Cognet, M. Le Bret, W. Guschlbauer and G.V. Fazakerley, Biochemistry 28, 2018 (1989), we determined by two dimensional NMR studies and molecular mechanics calculations the three-dimensional structure of a non-selfcomplementary oligonucleotide: [sequence; see text] where dr, at the center of the first strand, is a model abasic site. In order to explain all the results arising from NMR measurements, we found that an equilibrium between two conformations was necessary. These conformations differ mainly by the sugar pucker of G5 which is C2' endo or C3' endo. The latter is stabilized by addition of counterions between phosphate residues P3 and P4. In this paper, we have constructed systematically, all possible structures as a function of torsion angles delta of dr4 and of G5 by molecular mechanics in the presence or absence of counterions. Since these conformations were not forced with NMR distance measurements, this method allows detailed comparisons between all possible conformations and NMR data. Maps of contour lines of the potential energy, of fits to NMR distance measurements, and of helical twist as a function of torsion angles delta of dr4 and of G5 unravel the difficulties associated with the study of the G5 sugar pucker conformation equilibrium. Sugar puckers and proton distances are very sensitive criteria to monitor molecular dynamics. Relying on these experimental criteria, we have tested many molecular dynamics preparation phases and we propose a new warm-up and equilibration procedure for molecular dynamics. Thus we show with a 290 ps molecular dynamic run that G5 is in conformational equilibrium and that all NMR data are well reproduced.

A molecular mechanics and dynamics study of alternate triple-helices involving the integrase-binding site of the HIV-1 virus and oligonucleotides having a 3'-3' internucleotide junction.

Triple helix formation by oligonucleotides can be extended beyond polypurine tracts with the help of specially designed linkers. In this paper we focus our attention on the integrase-binding site of the HIV-1 virus located on the U5 LTR end which contains two adjacent purine tracts on opposite strands. Two alternate triple helices with a 3'-3' junction in the third strand are considered: 5'-GGTTTTp3'-3'pTGTGT-5' and 5'-GGAAAAp3'-3'pAGAGA-5' The structural plausibility of these triplexes is investigated using molecular mechanics and dynamics simulations, both in vacuo and in aqua. The non-isomorphism of the triplets in the GpT steps in the first sequence, gives rise to non canonical conformations in the torsion angles, hydration appears to be crucial for this triplex. Sugar puckers are predominantly South during in vacuo simulations while they turn East in aqua. In the simulation in aqua the triplexes are shrouded by an hydration shell, however, we have not been able to detect any permanent hydrogen bond bridge between DNA and water. The solvation of ions as well as their radial distribution, appear to be relatively well behaved despite the artifacts known to be generated by the simulation procedure. The experimental feasibility of these structures is discussed.

Reassessment of structural characteristics of the d(CGCG)2:actinomycin D complex from complete 1H and 31P NMR.

Complexes formed between Actinomycin D (ActD) and the tetranucleotides d(AGCT)2 and d(CGCG)2 were studied in detail by one and two-dimensional 1H and 31P NMR. The 31P two dimensional chemical exchange experiment, at room temperature on saturated complexes (1:1), showed unambiguously that the asymmetrical phenoxazone ring binds to the unique GC site under the two possible orientations in the d(AGCT)2 tetranucleotide but adopts a single orientation in the d(CGCG)2 tetranucleotide. For the d(CGCG)2:Act D saturated complex, complete assignments of all protons and phosphorus signals of the two-nucleotide strands, as well as of the two cyclic pentapeptide chains has allowed us to study in details the conformational features of the complex from NOE and coupling constants analysis. The tetranucleotide remains in a right-handed duplex, but the sugar puckers are modified for residues at the intercalation site. A uniform C2' endo pucker is observed for residues on the strand facing the quinoid side of the phenoxazone ring while a C2' endo-C3-endo equilibrium about 60% of C2' endo is proposed for the two residues on the strand facing the benzenoid side of the phenoxazone ring. In contrast to previous studies on ActD-DNA interactions, we have been able to measure the 3J phosphorus-proton coupling constants at the intercalation site but also adjacent to it, showing that 31P chemical shifts are not simply related to the backbone conformation. Molecular mechanics calculations, using empirical distances deduced from NOE effects as restrained distances during minimizations, led to a model differing mainly from those previously published by orientation of the N methyl groups of both N-Methyl-Valines.

Binding of Net-Fla, a netropsin-flavin hybrid molecule, to DNA: molecular mechanics and dynamics studies in vacuo and in water solution.

We have studied the binding of the hybrid netropsin-flavin (Net-Fla) molecule onto four sequences containing four A. T base pairs. Molecular mechanics minimizations in vacuo show numerous minimal conformations separated by one base pair. 400 ps molecular dynamics simulations in vacuo have been performed using the lowest minima as the starting conformations. During these simulations, the flavin moiety of the drug makes two hydrogen bonds with an amino group of a neighboring guanine. A 200 ps molecular dynamics simulation in explicit water solution suggests that the binding of Net-Fla upon the DNA substrate is enhanced by water bridges. A water molecule bridging the amidinium of Net-Fla to the N3 atom of an adenine seems to be stuck in the drug-DNA complex during the whole simulation. The fluctuations of the DNA helical parameters and of the torsion angles of the sugar-phosphate backbone are very similar in the simulations in vacuo and in water. The time auto-correlation functions for the DNA helical parameters decrease rapidly in the picosecond range in vacuo. The same functions computed from the water solution molecular dynamics simulations seem to have two modes: the rapid mode is similar to the behavior in vacuo, and is followed by a slower mode in the 10 ps range.

Quenching of DNA-ethidium fluorescence by amsacrine and other antitumor agents: a possible electron-transfer effect.

The antitumor agent amsacrine, 4'-(9-acridinylamino)methanesulfon-m-anisidide (m-AMSA), when bound to double-stranded DNA, particularly poly(deoxyadenylicthymidylic acid), reduced the fluorescence of bound ethidium without physically displacing it from DNA. Fluorescence lifetime measurements showed that the reduction of fluorescence was not due to reduction of the lifetime of the excited state of ethidium. Rather, a proportion of the DNA-bound ethidium changed to a state where the fluorescence was highly quenched. Several other 9-anilinoacridine derivatives, and also 9-hydroxyellipticine, caused quenching of ethidium-DNA fluorescence, whereas 9-aminoacridine, proflavin, and ellipticine had no effect. Resonance energy transfer (Förster transfer) is not responsible for the effect since there is no spectral overlap between the absorption spectrum of any of the agents and the fluorescence emission spectrum of ethidium. It is suggested that quenching may occur as a result of reversible formation of electron-transfer complexes between the intercalating drug and the excited state of ethidium.

Object Command Language: a formalism to build molecular models and to analyze structural parameters in macromolecules, with applications to nucleic acids.

We have written a programming language OCL (Object Command Language) to solve, in a general way, two recurring problems that arise during the construction of molecular models and during the geometrical characterization of macromolecules: how to move precisely and reproducibly any part of a molecular model in any user-defined local reference axes; and how to calculate standard or user-defined structural parameters that characterize the complex geometries of any macromolecule. OCL endows the user with three main capabilities: the definition of subsets of the macromolecule, called objects in OCL, with a formalism from elementary set theory or lexical analysis; the definition of sequences of elementary geometrical operations, called procedures in OCL, enabling one to build arbitrary three-dimensional (3D) orthonormal reference frames, to be associated with previously defined objects; and the transmission of these definitions to programs that allow one to display, to modify and to analyze interactively the molecular structure, or to programs that perform energy minimizations or molecular dynamics. Several applications to nucleic acids are presented.

Accessibility to bacterial nucleic acids of the intercalating drug aliphatic amino acid ellipticinium derivatives in Escherichia coli and Salmonella typhimurium.

The fluorescence of the aliphatic (amino acido)-N2-methyl-9-hydroxyellipticinium (AA-NMHE) derivatives [Auclair, C., Voisin, E., Banoun, H., Bernardou, J., Meunier, B., & Paoletti, C. (1984) J. Med. Chem. 27, 1161-1166], namely, dehydroglycino-NMHE, dehydroalanino-NMHE, dehydrovalino-NMHE, and dehydroleucino-NMHE, has been characterized. The changes in the fluorescence properties of the drugs, including increase in quantum yields, increase in fluorescence lifetimes, and occurrence of energy transfer upon binding to DNA in vitro, have been further investigated. The measurement of the fluorescence increment of AA-NMHE when bound to fluorescent sites inside intact bacteria has been found to be suitable for the determination of the accessibility of the drugs to bacterial nucleic acids according to the method of Lambert and Le Pecq [Lambert, B., & Le Pecq, J.B. (1984) Biochemistry 23, 166-176]. With this methodology, the kinetics of drug uptake, the ability of the drug to reach the bacterial nucleic acids at equilibrium, and the nature of the ligand binding model have been determined in two AA-NMHE-sensitive strains, Escherichia coli BL 101 (Lambert & Le Pecq, 1984) and Salmonella typhimurium TA 98 [Ames, B.N., Lee, F.D., & Durston, W.E. (1973) Proc. Natl. Acad. Sci. U.S.A. 70, 782-786]. The main results obtained are the following: At nonsaturating concentrations, each AA-NMHE exhibits a marked difference in its ability to reach the bacterial nucleic acids. This parameter seems to be correlated with the antibacterial efficiency of the drugs.(ABSTRACT TRUNCATED AT 250 WORDS)

Stability of G,A triple helices.

In this work we selected double-stranded DNA sequences capable of forming stable triplexes at 20 or 50 degrees C with corresponding 13mer purine oligonucleotides. This selection was obtained by a double aptamer approach where both the starting sequences of the oligonucleotides and the target DNA duplex were random. The results of selection were confirmed by a cold exchange method and the influence of the position of a 'mismatch' on the stability of the triplex was documented in several cases. The selected sequences obey two rules: (i) they have a high G content; (ii) for a given G content the stability of the resulting triplex is higher if the G residues lie in stretches. The computer simulation of the Mg2+, Na+and Cl-environment around three triplexes by a density scaled Monte Carlo method provides an interpretation of the experimental observations. The Mg2+cations are statistically close to the G N7 and relatively far from the A N7. The presence of an A repels the Mg2+from adjacent G residues. Therefore, the triplexes are stabilized when the Mg2+can form a continuous spine on G N7.

Solution structure of two mismatches A.A and T.T in the K-ras gene context by nuclear magnetic resonance and molecular dynamics.

Two mismatches, one homopurine (A.A) and the other homopyrimidine (T.T), have been incorporated at the central position N of: 5'd(GCCACNAGCTC).d(GAGCTNGTGGC) in order to study nuclear magnetic resonance spectra and molecular dynamics. These duplexes constitute the sequence 29-39 of the K-ras gene coding for Gly12, a hot spot for mutation. The NMR spectra show that the duplexes are not greatly distorted by the introduction of the mismatches and their global conformation is that of a canonical B-form double helix. For both systems, no structural change is observed in the pH range 4.7-9. For the duplex containing the homopurine A.A mismatch, we propose a type of pairing involving one hydrogen bond between the amino group of one central adenine and the nitrogen N1 of the opposite adenine. For the duplex containing the mispaired T.T bases, NMR spectra recorded in H2O at 282 K indicate that these central bases are engaged in wobble pairing, involving two imino-carbonyl hydrogen bonds. For both systems two conformations with the same donor and acceptor pattern can coexist, one being obtained from the other by a 180 degrees rotation about the pseudodyadic axis. Exchange between the two forms is observed by NMR at low temperature for the T.T mispair and also inferred from NMR measurements on the A.A system. The presence of this exchange and its pathway has been investigated by molecular dynamics calculations on both systems. Distance restrained and unrestrained molecular dynamics are in very good agreement with the NMR data. The average structure for either mispair shows only small conformational change from normal B DNA. For each, a systematic pathway is observed for exchange between the two conformations.

Sequence requirements for mammalian topoisomerase II mediated DNA cleavage stimulated by an ellipticine derivative.

Various antitumor drugs stabilize DNA topoisomerase II-DNA transient covalent complexes. The complexes distribution along pBR322 DNA was shown previously to depend upon the nature of the drug (Tewey et al. (1984) Science 226, 466-468). The position in pBR322 of DNA cleavage by calf DNA topoisomerase II for 115 such sites stabilized by an ellipticine derivative and the relative frequency of cleavage at most of these sites were determined. The nucleotide sequence surrounding the 25 strongest sites was analyzed and the following ellipticine specific consensus sequence was deduced: 5'-ANCNT(A/G)T.NN(G/C)N(A/G)-3' where cleavage occurs at the indicated mark. A thymine is always present at the 3' end of at least one strand of the strong cleavage sites, and the dinucleotide AT or GT at the 3' end of the break plays a major role in the complex stabilisation. The predictive value of cleavage of the consensus was tested for two regions of SV40 DNA and cleavage was indeed detected at the majority of the sites matching the consensus. Some complexes stabilized by ellipticine are resistant to salt dissociation and this property seems to be correlated with the presence of symmetrical sequences in the cleavage site with a center of symmetry staggered relatively to the center of symmetry of cleavage.

Electron microscopic measurement of chain flexibility of poly(dG-dC).poly(dG-dC) modified by cis-diamminedichloroplatinum(II).

The antitumor drug cis-diamminedichloroplatinum (II) (cis-Pt) forms bidentate adducts with guanine residues of poly(dG-dC).poly(dG-dC). The secondary structure of the polymer is altered. In this work, high resolution pictures of naked molecules, obtained by dark field electron microscopy reveal DNA chain distortions with radii as small as 30 A. The extent of distortion increases with the drug/nucleotide ratio (rb). These alterations of the secondary structure are responsible for the apparent shortening of the molecules. Measurements of the persistence lengths of the polymer as well as the end-to-end distances of elementary segments of various lengths, are obtained from digitized electron micrographs. The measurements are used to monitor and quantify the observed modifications of polymer structure upon cis-Pt binding at various rb or incubation times. Poly(dG-m5dC).poly(dG-m5dC) in the B and Z forms have different persistence lengths. In the B form, this polymer is more altered by cis-Pt than in the Z one.