HIV-1 integrase catalytic core: molecular dynamics and simulated fluorescence decays.

Two molecular dynamics simulations have been carried out on the HIV-1 integrase catalytic core starting from fully determined crystal structures. During the first one, performed in the absence of divalent cation (6-ns long), the catalytic core took on two main conformations. The conformational transition occurs at approximately 3.4 ns. In contrast, during the second one, in the presence of Mg(2+) (4-ns long), there were no such changes. The molecular dynamics simulations were used to compute the fluorescence intensity decays emitted by the four tryptophan residues considered as the only chromophores. The decay was computed by following, frame by frame, the amount of chromophores that remained excited at a certain time after light absorption. The simulation took into account the quenching through electron transfer to the peptide bond and the fluorescence resonance energy transfer between the chromophores. The fit to the experimental intensity decays obtained at 5 degrees C and at 30 degrees C is very good. The fluorescence anisotropy decays were also simulated. Interestingly, the fit to the experimental anisotropy decay was excellent at 5 degrees C and rather poor at 30 degrees C. Various hypotheses such as dimerization and abnormal increase of uncorrelated internal motions are discussed.

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.

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.

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.

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.

Hydration of the dTn.dAn x dTn parallel triple helix: a Fourier transform infrared and gravimetric study correlated with molecular dynamics simulations.

We present a comparative analysis of the water organization around the dTn.dAn x dTn triple helix and the Watson-Crick double helix dTn.dAn respectively by means of gravimetric measurements, infrared spectroscopy and molecular dynamics simulations. The hydration per nucleotide determined by gravimetric and spectroscopic methods correlated with the molecular dynamics simulations shows that at high relative humidity (98% RH) the triple helix is less solvated than the duplex (17 +/- 2 water molecules per nucleotide instead of 21 +/-1). The experimental desorption curves are different for both structures and indicate that below 81% RH the triplex becomes more hydrated than the duplex. At this RH the FTIR spectra show the emergence of N-type sugars in the adenosine strand of the triplex. When the third strand is bound in the major groove of the Watson-Crick duplex molecular dynamics simulations show the formation of a spine of water molecules between the two thymidine strands.

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.

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.

Sequence-directed single strand cleavage of DNA by a netropsin-flavin hybrid molecule.

In an attempt to obtain sequence specific DNA-cleaving molecules, we have synthesized a series of hybrid minor groove binders composed of a photoactiveable isoalloxazine (flavin) chromophore linked through a polymethylenic chain to a bis-pyrrolecarboxamide moiety related to netropsin. Like netropsin, the hybrid derivatives preferentially bind to A+T-rich sequences. Activation of the flavin chromophore by visible light results in the appearance of single strand breaks in the vicinity of the DNA binding site. We have further investigated the cleavage affinity properties of one of these compounds referred to as netropsin-flavin (Net-Fla) and considered as representative of the series. Net-Fla cleaves only one strand at a specific locus downstream of 5'-AAAT-3', upstream of 5'-TAAA-3' and on either side of a 5'-AAAA-3' sequence. Net-Fla cleaves both strands downstream to 5'-AATT-3'. This makes the properties of Net-Fla similar to that of a restriction endonuclease and provides additional insight into establishing the rules for the readout of B-DNA helix by non-nucleotidic compounds. Using molecular modeling, we show that Net-Fla binds to an asymmetric site in one orientation. The values of the energetic minima lie in the same order as expected from the cleavage patterns, which suggests that the oriented cleavage is a consequence of a sequence-oriented binding of Net-Fla in the DNA minor groove.

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.

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+.

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.

Triple helix formation with short oligonucleotide-intercalator conjugates matching the HIV-1 U3 LTR end sequence.

In an attempt to target short purine sequences in view of pharmacological application, we have synthesized three new TFO (triple-helix-forming oligonucleotide) conjugates in which an intercalating oxazolopyridocarbazole (OPC) chromophore is linked by a pentamethylene linker to a 7-mer oligonucleotide matching the polypurine/polypyrimidine sequence located in the HIV-1 U3 LTR end region. The TFO moiety of conjugates are 5'CCTTCCC, 5'GGGAAGG, and 5'GGGTTGG. Their ability to bind to double-stranded DNA targets was examined. This binding is demonstrated by a footprinting technique using DNase I as a cleaving agent. The complex involved intermolecular pyr-pur*pyr or pur-pur*pyr triple helix. Pyrimidine TFO-OPC binds in a pH-dependent manner, whereas the others do not. The formation of the complex has been investigated at neutral pH and increasing temperature. We observed that the protection due to the purine and mixed TFO-OPC was pH independent and remained identical up to 40 degrees C. To determine the position of the OPC chromophore, molecular modeling was undertaken on the purine-conjugate/target complex. It has been suggested that the complex involved the intercalation of the OPC at the triplex-duplex junction with a small unwinding at the next excluded site.

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.

The pH dependent configurations of the C.A mispair in DNA.

The structure of the cytosine-adenine mispair in a 7 base pair duplex has been investigated by proton NMR spectroscopy. At low pH, the predominant structure is protonated on the A residue and assumes a wobble conformation consistent with previous findings. The C residue of the mispair is found in a C2'-C3' endo equilibrium. This is confirmed by molecular dynamics calculations which suggest that the conformation of the protonated wobble is flexible and not as rigid as a normal base pair. As the solution pH is raised, a structural transition is observed with an apparent pK of 7.54 at 23 degrees C. At higher pH the predominant structure is one in which both the C and A residues are intrahelical. Evidence is presented that this structure corresponds to a reverse wobble in which the two bases are held together by one hydrogen bond. This structure is much less stable than the protonated form and even at low temperature single strands are observed in slow exchange with the neutral duplex form.

Solution conformation of an oligonucleotide containing a G.G mismatch determined by nuclear magnetic resonance and molecular mechanics.

We have determined by two-dimensional nuclear magnetic resonance studies and molecular mechanics calculations the three dimensional solution structure of the non-selfcomplementary oligonucleotide, d(GAGGAGGCACG). d(CGTGCGTCCTC) in which the central base pair is G.G. This is the first structural determination of a G.G mismatch in a oligonucleotide. Two dimensional nuclear magnetic resonance spectra show that the bases of the mismatched pair are stacked into the helix and that the helix adopts a classical B-DNA form. Spectra of the exchangeable protons show that the two guanosines are base paired via their imino protons. For the non-exchangeable protons and for some of the exchangeable protons nuclear Overhauser enhancement build up curves at short mixing times have been measured. These give 84 proton-proton distances which are sensitive to the helix conformation. One of the guanosines adopts a normal anti conformation while the other is syn or close to syn. All non-terminal sugars are C2' endo. These data sets were incorporated into the refinement of the oligonucleotide structure by molecular mechanics calculations. The G.G mismatch shows a symmetrical base pairing structure. Although the mismatch is very bulky many of its features are close to that of normal B-DNA. The mismatch induces a small lateral shift in the helix axis and the sum of the helical twist above and below the mismatch is close to that of B-DNA.

The solution structure of a DNA hairpin containing a loop of three thymidines determined by nuclear magnetic resonance and molecular mechanics.

We have determined by two-dimensional nuclear magnetic resonance studies and molecular mechanics calculations the three-dimensional solution structure of a 21 residue oligonucleotide capable of forming a hairpin structure with a loop of three thymidine residues. This structure is in equilibrium with a duplex form. At 33 degrees C, low ionic strength and in the presence of MgCl2 the hairpin form dominates in solution. Six Watson-Crick base pairs are formed topped by the loop structure. The residues 1-3 and 18-21 are not complementary and form dangling ends. Distance constraints have been derived from nuclear Overhauser enhancement measurements. These, together with molecular mechanics calculations, have been used to determine the structure. We do not observe stacking of thymidine residues either over the 3' or the 5' end of the stem.

Comparison of the bis-intercalating complexes formed between either ditercalinium or a flexible analogue and d(CpGpCpG)2 or d(TpTpCpGpCpGpApA)2 minihelices: 1H- and 31P-NMR analyses.

The 400-MHz 1H- and 162-MHz 31P-nmr have been used to study complexes constituted by (a) the d(TpTpCpGpCpGpApA)2 or the d(CpGpCpG)2 self-complementary oligonucleotides and (b) two bifunctional 7H-pyrido [4,3-c] carbazole dimer drugs, the antitumoral ditercalinium (NSC 366241), a dimer with a rigid bis-piperidine linking chain and its pharmacologically inactive analogue, a dimer with a flexible spermine-like linking chain. Nearly all proton and phosphorus signals have been assigned by two-dimensional (2D) nmr (correlated spectroscopy, homonuclear Hartmann-Hahn, nuclear Overhauser enhancement spectroscopy, 2D 31P (1H) heteronuclear correlated spectroscopy and 31P-31P chemical exchange experiments). Both drugs bis-intercalate into the two CpG sites. The complexes show small differences in the position of the 7H-pyrido [4,3-c] carbazole ring into the intercalation site and possibly in the ribose-phosphate backbone deformation. However, the inactive analogue exhibits a longer residence lifetime in octanucleotide than the ditercalinium does. All these results are discussed in terms of differences in dimer activities.

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.