Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Unprocessed Viral DNA Could Be the Primary Target of the HIV-1 Integrase Inhibitor Raltegravir

  • Farah F. Ammar,

    Affiliations LBPA, UMR8113 du CNRS, Ecole Normale Supérieure de Cachan, Cedex, Cachan, France, Unité de Biochimie, Département SVT, Faculté des Sciences, Université Saint-Joseph, CST-Mar Roukoz, Beyrouth, Liban

  • Safwat Abdel-Azeim,

    Affiliation LBPA, UMR8113 du CNRS, Ecole Normale Supérieure de Cachan, Cedex, Cachan, France

  • Loussinée Zargarian,

    Affiliation LBPA, UMR8113 du CNRS, Ecole Normale Supérieure de Cachan, Cedex, Cachan, France

  • Zeina Hobaika,

    Affiliation Unité de Biochimie, Département SVT, Faculté des Sciences, Université Saint-Joseph, CST-Mar Roukoz, Beyrouth, Liban

  • Richard G. Maroun,

    Affiliation Unité de Biochimie, Département SVT, Faculté des Sciences, Université Saint-Joseph, CST-Mar Roukoz, Beyrouth, Liban

  • Serge Fermandjian

    serge.fermandjian@lbpa.ens-cachan.fr

    Affiliation LBPA, UMR8113 du CNRS, Ecole Normale Supérieure de Cachan, Cedex, Cachan, France

Abstract

Integration of HIV DNA into host chromosome requires a 3′-processing (3′-P) and a strand transfer (ST) reactions catalyzed by virus integrase (IN). Raltegravir (RAL), commonly used in AIDS therapy, belongs to the family of IN ST inhibitors (INSTIs) acting on IN-viral DNA complexes (intasomes). However, studies show that RAL fails to bind IN alone, but nothing has been reported on the behaviour of RAL toward free viral DNA. Here, we assessed whether free viral DNA could be a primary target for RAL, assuming that the DNA molecule is a receptor for a huge number of pharmacological agents. Optical spectroscopy, molecular dynamics and free energy calculations, showed that RAL is a tight binder of both processed and unprocessed LTR (long terminal repeat) ends. Complex formation involved mainly van der Waals forces and was enthalpy driven. Dissociation constants (Kds) revealed that RAL affinity for unbound LTRs was stronger than for bound LTRs. Moreover, Kd value for binding of RAL to LTRs and IC50 value (half concentration for inhibition) were in same range, suggesting that RAL binding to DNA and ST inhibition are correlated events. Accommodation of RAL into terminal base-pairs of unprocessed LTR is facilitated by an extensive end fraying that lowers the RAL binding energy barrier. The RAL binding entails a weak damping of fraying and correlatively of 3′-P inhibition. Noteworthy, present calculated RAL structures bound to free viral DNA resemble those found in RAL-intasome crystals, especially concerning the contacts between the fluorobenzyl group and the conserved 5′C4pA33′ step. We propose that RAL inhibits IN, in binding first unprocessed DNA. Similarly to anticancer drug poisons acting on topoisomerases, its interaction with DNA does not alter the cut, but blocks the subsequent joining reaction. We also speculate that INSTIs having viral DNA rather IN as main target could induce less resistance.

Introduction

Integration of the HIV-1 DNA into the host chromosome leads to the viral infection at the origin of the AIDS pandemic. Integration is catalysed by the retroviral enzyme integrase (IN) [1][4]. The whole integration involves the 3′-processing (3′-P) and the strand transfer (ST), which occur in the cytoplasm and in the nucleus, respectively. The integration is finalized by the cell enzymes which cleave the viral DNA 5′-overhang and fill the room left between the viral and cellular DNA [3][5]. A huge effort from both the public research and the pharmaceutical industry was made during the last decade to discover IN inhibitors. Only those acting on the ST step have emerged as interesting antiretroviral drugs [5][7]. Thus, Merck and Co has recently developed the raltegravir (RAL, MK-0518), a potent INSTI (IN ST Inhibitor) that derives from DKAs (Diketo Acids) and which is now widely used in AIDS therapy [8], [9]. However, RAL induces mutations located mainly into the loop 140 (Y143H/R/C, Q148H/R/K and G140S-Q148H) and the α4 helix (N155H) [10][12], entailing significant clinical resistance [11]. Positions of these mutations in the protein are consistent with the determining role hold by the α4 helix (residues 150 to 166) [13], [14] and the loop 140 (residues 140 to 149) [11], [15] of the catalytic core domain (CCD) in the IN activity and also as sites of inhibitors [15][18]. The crystal structures recently resolved by Hare et al., 2010 confirm the coordination of two Mg2+ to the three coplanar oxygen or nitrogen atoms of the metal binding motif, while the halogenated aromatic ring penetrates more or less deeply into the space made available by the opening of the conserved 5′C4pA33′ step at the end of the processed strand. The oxadiazole moiety of the RAL molecule further interacts through stacking interactions with the Tyr 143 phenolic group. However, this “π-π” stacking does not appear indispensable for the inhibitory activity of INSTIs as it can be replaced by other interactions with the loop 140 amino acid residues, as shown by the oxadiazole-lacking compounds ie EVG (elvitegravir) [15], [19], [20] and the so-called new generation of INSTIs such as MK-0536 [18], [21], MK-2048 [18], [22] and DTG (dolutegravir) [23], [24]. These compounds not only are fully active against WT-IN but they also remain effective against the RAL-resistant Y143R mutant of IN [20]. All display intermolecular contacts (“π-π”stacking) between their halogenated aromatic ring and the cytosine base of the conserved 5′C4pA33′step. The subsequent spatial displacement of the adjacent 3′-adenine (A3) bearing the functional hydroxyl group from its initial position is considered as the main event promoting the inhibition of integration [15].

Actually, it is the large number of interactions occurring between INSTIs and the viral DNA, but also the inability of RAL to bind tightly to IN taken alone [25], [26], which motivated us to examine whether the LTR ends could be the primary targets of the drugs. To this end we used RAL and several oligonucleotides mimicking or deriving from the U5 LTR extremity of viral DNA (Fig. 1). Analysis of the drug-DNA complexes was performed by UV-absorbance [27], circular dichroism (CD) [27], fluorescence [13], molecular dynamics simulations (GROMACS 4.5.3/Amber99SB-ILDN) [28], [29] and free energy calculations using the Molecular Mechanics-Poisson Boltzmann Surface Area method (MMPBSA) [30][33]. Results indicate that one molecule of RAL binds tightly to 3′-processed LTR (Kd≈6 nM) and more weakly to unprocessed LTR (Kd≈20 nM). Binding of RAL to processed LTR requires several key nucleotides including the 5′A−1C−23′ overhang, known for its strong implication in the binding and activity of IN [34]. The binding of RAL to this small dinucleotide strand in unprocessed LTR is facilitated by a major fraying in terminal base pairs that lowers the energy barrier for drug insertion [35], [36]. The insertion of RAL into the terminal base pairs, affects slightly their fraying and similarly the 3′-P reaction, the latter reaction being closely correlated to the motions at the LTR ends [35]. After the deletion of the 5′G2T13′ dinucleotide, RAL blocks the ST reaction in adopting a new position at the LTR end more conducive to binding with the two divalent cations and the cytosine C4 of the conserved malleable C4pA3 step [15], [17], [18]. Remarkably, in our two modeled RAL-LTR32 structures, the adenine A3 nucleotide that bears the essential 3′OH group has conformations similar the ones found in the crystal structures of the RAL-intasome complexes (PDB codes: 3L2T [15] and 3OYA [18]). All together, present results bring greater clarity on the inhibitory mechanism of INSTIs, especially in showing that the drugs may bind specifically to both the unprocessed and processed LTR ends. The binding of RAL to unprocessed LTR does not or little impair the 3′-P reaction because the end fraying required for the cleavage of the phosphodiester backbone by IN is only weakly altered by RAL [37]. The 3′-P reaction produces a change of the complex conformation, responsible for the blocking of the joining reaction [38], [39]. We also propose that anti-AIDS drugs having an increased number of interactions with the substrate viral DNA, at detriment of the protein active site, could induce less resistance mutations.

thumbnail
Figure 1. Molecules used in this study.

(A) The oligonucleotides are designed to adopt a hairpin structure, folded around a three thymine loop whose central thymine bears the fluorescein reporter. LTR34: unprocessed U5-LTR end with a 17 base pair stem. The numbering of the four last base pairs (+1 to +4 and −1 to −4) starts from the ultimate 3′nucleotide on the upper (+) strand and from the ultimate 5′ nucleotide on the lower (−) strand; LTR32, processed U5-LTR, with 5′A−1C−23′ as overhang on the (−) strand; LTR32-I, inversed LTR32, with 5′G2T13′ as overhang in the (+) strand; and LTR30, doubly deleted LTR34 (blunt-ended DNA). (B) Chemical structure of the here studied RAL (MK-0518). (C) Chemical structure of MK-2048. This compound is given as an example of INSTI of second generation, inducing less resistance mutations in IN. (D) PFV LTR sequences used in calculations. The numbering of the four last base pairs (+1 to +4 and −1 to −4) is the same as in (A).

https://doi.org/10.1371/journal.pone.0040223.g001

Materials and Methods

Oligonucleotides and RAL

The oligonucleotides LTR34, LTR32, LTR32-I and LTR30 (Fig. 1 A) were purchased from Eurogentec (Belgium). They were designed to adopt a monomolecular hairpin structure that remains stable at the low concentrations used in fluorescence and CD experiments (10−9 to 10−5 M). LTR34 reproduces the unprocessed version of the U5 LTR end, LTR32 is the 3′-processed version obtained by deletion of G2T13′ and LTR32-I (LTR32-inverted) is obtained by deletion of C−2A−15′; blunt-ended LTR30 is obtained by deletion of both G2T13′ and C−2A−15′. The thymine at the centre of the three thymine loop bears the fluorescein reporter for fluorescence studies. Unlabeled oligonucleotides were also prepared for UV-absorbance and CD measurements. RAL (Fig. 1 B) was purchased from CacheSyn while MK-2048 (Fig. 1 C) is given as an example of an INSTI of the new generation. PFV LTR sequences used in calculations are indicated in Fig. 1D.

UV-absorbance measurements

UV-spectrometry experiments were recorded by using an Uvikon spectrophotometer model 941 (Kontron Instruments). DNA and RAL samples were dissolved in phosphate buffer (Na/Na2 phosphate, 10 mM, pH 6, I = 0.05) in the presence of 5 mM MgCl2. Titrations were performed using the DNA as the titrant in 10 mm and 2 mm path length quartz cells and scanning the spectrum after each aliquot addition. The RAL concentration was generally maintained at 20 µM and that of DNAs was varied from 1 to 20 µM. Difference spectra between 200 nm and 380 nm were obtained by subtraction of DNA spectra from the DNA-RAL complex spectra, after subtraction of the buffer contribution.

CD measurements

CD spectra were recorded on a Jobin-Yvon CD6 dichrograph. Measurements were calibrated with (+)-10-camphorsulfonic acid. Samples were dissolved in phosphate buffer. The concentration of oligonucleotides was 10 µM while that of the titrant was varied from 10 µM to 80 µM. Samples were placed in thermally jacketed cuvettes with 1–5 mm path lengths. Spectra, recorded with 1-nm steps and corrected for the base line, were averaged over 10 scans. Before spectral recording, samples were incubated 10 min at the chosen temperature to allow the solutions to reach their equilibrium state. Spectra of DNAs and complexes of RAL-DNA were presented as CD per residue, Δε (M−1.cm−1), as a function of wavelength, λ (nm), between 200 and 330 nm. As RAL lacks chirality it does not directly contribute to the spectrum. Yet, RAL bears several chromophores which once placed in an asymmetric environment, i.e. in the vicinity of a deoxyribose ring, can acquire chirality and generate a signal in the absorption region.

Fluorescence measurements

These included both fluorescence intensity and anisotropy titration studies. Measurements were carried on a Jobin-Yvon Fluoromax II instrument. In fluorescence intensity titrations, RAL was maintained at a constant concentration while the DNA was used as the titrant. For RAL the wavelengths for maximum excitation and emission were λ = 313 nm and λ = 413 nm, respectively. For fluorescence anisotropy (A = (IIII)/(III+2I)) the parallel (III) and perpendicular (I) emission components were measured in L-format. The denominator of A was simply the total light that would be observed if no polarizers were used. With fluorescein as fluorophore grafted on DNA, the excitation from the xenon lamp (150-watt ozone-free) was performed at 488 nm with a 4-nm slit width. The emission was recorded at 516 nm with a 5-nm slit width in the case of LTR34, at 515 nm with a 4-nm slit width in cases of LTR32 and LTR32-I and at 514 nm with a 4-nm slit width in the case of LTR30. The fluorescein-labeled oligonucleotides were diluted to the desired concentration in 800 µl of phosphate buffer at the selected temperature (generally 5°C). Samples were placed in thermally jacketed 1-cm×0.5-cm quartz cuvettes, and measurements, at least 10 data points for each titration point, were recorded with an integration time of 1 s. For each fluorescence anisotropy measurement, the parallel (III) and the perpendicular (I) intensities of the background solution (i.e. buffer and RAL contributions) were subtracted from the sample value. The validity of fluorescence anisotropy measurements was controlled in measuring the total fluorescence intensity in parallel to fluorescence anisotropy. Variations of fluorescence intensity during these experiments were very weak, so we considered that the anisotropy signal contained the desired information on the complex formation. Kds (equilibrium dissociation constants) were calculated by fitting the sigmoidal curves, using GraphPad Prism 5 applying either the linear regression or non-linear regression (curve fit) “Least Squares” procedure. Binding stoichiometries were determined using the Bujalowsky and Lohman procedure [40]. The reverse experiment consisting in the analysis of the oligonucleotide binding to the drug was also carried out.

Molecular dynamics simulations of viral DNA-RAL complexes

Molecular dynamics (MD) simulations were performed using the GROMACS software (version 4.5.3) [28] with the Amber99SB-ILDN force field [29]. RAL parameters were constructed using the ACPYPE (AnteChamber PYthon Parser interfacE) [41], [42], the General Amber Force Field (GAFF) [41] and the restrained electrostatic potential (RESP) charges. RESP charges were calculated using the Antechamber program [42] using the ESP charges calculated at Hartree-Fock [43] 6–31G* level [44] using G03(Gaussian03, Revision C.02). Simulations on RAL-DNA complexes were carried out on both processed and unprocessed PFV DNAs (i.e. LTR32 and LTR34) with a single Mg2+ ion. Initial coordinates of the LTR32-RAL complex were extracted from the 3OYA [18] PDB structure, after IN removal. Unprocessed DNA, LTR34 of PFV (Fig. 1D), was constructed by addition of the 5′AT dinucleotide (corresponding to 5′G2T13′ in HIV) at the 3′-end of the processed strand, using the VMD program [45]. Complexes were slightly relaxed using 50 steps of steepest descent minimizer. Subsequently, the system was immersed in an explicit water box of TIP3P model [46] which extended at least 16 Å away in each direction from any DNA or RAL atom. The systems include 14993 water molecules and 35 Na ions in the case of LTR34 (46274 atoms) and 13474 water molecules and 33 Na ions in the case of LTR32 (41651 atoms).

Sodium ions were added to neutralize the system as needed for the Particle Mesh Ewald (PME) [47] calculation of the long-range electrostatic interactions, while cut-off of 10 Å was used for van der Waals (VDW) and short-range electrostatic interactions. The system was exposed to 500 step of steepest descent minimization to remove the bad contacts with the solvent. All bonds involving hydrogen atoms were constrained by LINCS algorithm [48]. Equilibration of the solvent and ions around the complexes with position constraints of the heavy atoms, were performed for two nanoseconds in the constant Number of particles, Volume, and Temperature (NVT) ensemble and in Constant Number of particles, Pressure and Temperature (NPT) thermodynamic ensemble respectively. NVT simulations were carried out using the velocity rescaling thermostat (V-rescale) [49] and the NPT using Parrinello-Rahman barostat [50] MD production simulations were performed for a total of 400 ns duration in the NPT ensemble. Moreover, we also performed a calculation of the distance (nm) evolution between the center of mass of RAL and terminal bases (5′A and 3′T) of the PFV LTR34, in the calculation course from t = 0 ns to 100 ns.

Binding free energy calculations

RAL binding free energies were estimated using the end point Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) method [30][33], [51]. The total free energy (G) in MMPBSA analysis for a given species (complex, receptor, and ligand) was determined using Eq. (1); the overall change in free energy for complex formation (ΔGbind) for a non-covalent binding event was calculated according to Eq. (2).(1)(2)The polar solvation energies (Gpolar) were computed in continuum solvent using Poisson-Boltzmann (PB) and ionic force of 0.05 as used in the experiments. The non-polar terms (Gnonpolar+γSASA+β) were estimated using solvent accessible surface areas (SASA in Å [29]) with typical values for γ = 0.00542 kcal/mol Å2 [29] and β = 0.92 kcal/mol. The EMM term represented the sum of the electrostatic (Coulombic), VDW (Lennard-Jones), and internal energies (bonds, angles, and dihedrals). The remaining term represented temperature (T) and solute entropy (S), which can be estimated from normal-mode analysis of energy-minimized structures or quasi-harmonic (QH) modes over stabilized region of MD trajectory. 500 snapshots were selected from the last 5 ns, by keeping the snapshots every 5 ps. The entropy contributions were estimated by QH mode using 5000 frames from the last 5 ns. The free energy analysis was carried out using MMPBSA.py script from Amber11 program.

Results and Discussion

DNA is a target for a wide diversity of ligands. Due to the complexity and variety of DNA structures, the binding modes of these ligands, including anticancer agents, are them also complex and varied (intercalation, groove binding, insertion in breaks…) [27], [38], [39], [52][54]. Often, a same molecule (for instance ellipticine) can act as both an intercalator and a groove binder [27]. Moreover, anticancer agents such as camptothecin and derivatives, poisons of topoisomerase I, operate through insertion into a break of the DNA double helix created by the enzyme in one DNA strand [55]. Most of results stipulate that base pairs, double helix grooves and strand breaks in DNAs can be primary binding sites for topoisomerases inhibiting drugs [39]. As IN, similarly to topoisomerases, is also a DNA cutting and joining enzyme, we decided to investigate the binding of its best known inhibitor, RAL, to DNA. To this end, we used UV absorbance, CD, fluorescence and molecular dynamics with oligonucleotides to assess the binding of the drug to viral DNA LTR ends with respect to its INSTI activity.

UV-absorbance measurements

The first evidence of an interaction of RAL with the terminal part of viral DNA either unprocessed or processed is provided by UV-absorption titrations (Fig. 2 and Fig. S1 A and B). UV spectra of selected oligonucleotides, LTR34, LTR32 and LTR30 (Fig. 1 A), between 200 and 380 nm, display a main signal centred at about 260 nm and an additional peak around 200 nm, characteristics of B-DNA (Fig. S1 A). The UV spectrum of RAL (Fig. 1 B) between 200 nm and 380 nm consists in two peaks at 210 nm and 313 nm and a shoulder at 245 nm (Fig. S1 A), corresponding to the contributions of the aromatic chromophores making up the molecule. Fig. 2, shows the variations of the RAL spectrum resulting from addition of 20 µM LTR32, LTR34 and LTR30 to 20 µM RAL (stoichiometry 1∶1) after subtracting the spectrum of free DNA at the same concentration. Completely different effects are observed. With addition of LTR32 there is an emergence of a signal at about 260 nm, exactly where the DNA contributes, consistent with a change of conformation in LTR32 in response to drug-DNA complex formation (Fig. S1 B). The increase of intensity of the band at ≈205–210 nm, could arise from changes in the contributions of both RAL and DNA, due to their interaction. In contrast, the broad signal of RAL centred at 313 nm did not manifest any change, suggesting that the RAL chromophore generating this signal remains free of interactions in the complex. Noteworthy, addition of unprocessed LTR34 to RAL produces the same effects as shown by LTR32, but, however, less intense, while addition of blunt ended LTR30 is without effects. Taken together the UV-absorption results shows that RAL binds to the terminal part of LTRs, with a preference for the 3′ processed DNA that carries the 5′A−1C−23′dinucleotide overhang. In unprocessed LTR34 the 5′A−1C−23′ dinucleotide is involved in a duplex with the undeleted 5′ G2T13′ dinucleotide, but this does not constitute a rigid barrier to the binding of RAL. This can be explained by the important end fraying in unprocessed LTR [35], [36] that facilitates the RAL accommodation into the terminal base pairs. The inability of RAL to bind the blunt ended oligonucleotide LTR30, which is further devoid of both IN binding and ST activities [34], confirms the functional importance of the terminal 5′A−1C−23′ step within either the unprocessed or the processed viral DNA, regarding the capture of ligands.

thumbnail
Figure 2. UV-absorption analysis of oligonucleotides.

Spectra of RAL 20 µM (black) together with 20 µM LTR32 (green), 20 µM LTR34 (blue), 20 µM LTR30 (red), in phosphate buffer pH 6, I = 0.05, and MgCl2 5 mM final concentration.

https://doi.org/10.1371/journal.pone.0040223.g002

CD spectroscopy measurements

The CD technique is widely used to determine the secondary structures of proteins and nucleic acids and to follow the conformational changes induced by their mutual association or their binding to any type of ligand [56], [57]. Here, CD was applied to study the binding of RAL to LTR34 and LTR32, which mimic the unprocessed and processed viral U5 DNA ends, respectively (Fig. 1 A). The spectra of LTR34 and LTR32 presented in Fig. 3 show two main signals (negative at ≈250 nm and positive at ≈280 nm) characteristics of B-DNA. In titration experiments, RAL was added (concentrations from 10 µM to 80 µm) to oligonucleotides maintained at 10 µM concentration. The spectra recorded at drug: DNA ratios of 1, 2, 4, 6 and 8 showed a gradual variation of intensity of the two B-DNA signals at ≈250 nm and ≈280 nm. Similarly to UV-absorption experiments, effects were larger with LTR32 than LTR34. However, the signals were not shifted and no new signal induced by drug chromophores buried in a chiral environment was detected. The CD changes were assigned to a rearrangement of the nucleotide bases at the RAL binding site or contiguous to it [27], [58]. We will subsequently observe that it is a single RAL molecule that is inserted at the LTR extremity, so that one cannot expect the generation of a large signal.

thumbnail
Figure 3. Circular dichroism analysis of oligonucleotides-drug complexes.

Spectra of LTR34 (A) and LTR32 (B) at 10 µM (black) and difference spectra [LTR32/34 (10 µM)+RAL (10 µM, red; 20 µM, green; 40 µM, blue; 60 µM, orange; and 80 µM, purple)−LTR32/34 (10 µM)], in phosphate buffer pH 6, I = 0.05, and MgCl2 5 mM final concentration.

https://doi.org/10.1371/journal.pone.0040223.g003

Fluorescence measurements

The quantification of interactions stabilizing partner molecules is essential to the understanding of the complex formation. Fluorescence intensity and anisotropy measurements are well suited to a quantitative analysis of complexes as long as one of the binding partners is fluorescent. The fluorescence anisotropy gives information on both the stoichiometry of the complex and the binding constant: Kd = [L]×[R]/[LR] (Kd: dissociation constant; L: ligand and R: receptor) [13], [26], [59], [60]. The Kds provided by fluorescence anisotropy are confirmed by fluorescence intensity experiments, where the fluorescent ligand molecule is titrated with increasing concentrations of oligonucleotides.

The equilibrium saturation curves of the four fluorescein labeled oligonucleotides and increasing concentrations of RAL recorded by fluorescent anisotropy are shown in Fig. 4 and 5. The results reported in Fig. 4 confirm the above UV-absorbance and CD experiments, which indicated that RAL was able to interact with processed LTR32 and unprocessed LTR34, but fail to interact with the 3′ dangling ended LTR32-I and the blunt ended LTR30. Kds determined at midpoints are of ≈6 nM for LTR32 and ≈20 nM for LTR34. Rather similar values were obtained with the fluorescence intensity approach, using the oligonucleotides to titrate RAL (Fig. S2). The stoichiometry for the binding of RAL to LTR32 was determined by titrations at three different DNA concentrations (9 nM, 20 nM and 30 nM) (Fig. 5). The monophasic curves reached a same plateau after a variation of anisotropy of ΔA≈0.015. The three Kd values were quite similar and provided a mean value of ≈6 nM. Application of the Bujalwski and Lohman procedure [40], showed that a single RAL molecule was bound to LTR32 (1∶1 stoichiometry), in agreement with the crystal structure results [15]. Noteworthy, the experimental Kd for the binding of RAL to intasome [61] is higher the one found for LTR32 alone (19 nM vs 6 nM).

thumbnail
Figure 4. Quantitative analysis of RAL binding to oligonucleotides.

Fluorescence anisotropy titration of the four oligonucleotides LTR32 (black), LTR34 (red), LTR32-I (blue) and LTR30 (green) at 20 nM by increasing concentrations of RAL (from 10−12 M to 10−4 M). Kds obtained from titrations of LTR32 and LTR34 at 20 nM are indicated near the corresponding curves.

https://doi.org/10.1371/journal.pone.0040223.g004

thumbnail
Figure 5. Thermodynamic parameters for the binding of RAL to LTR32.

Titration of LTR32, at three different concentrations: 9 nM (black), 20 nM (red), and 30 nM (blue). Curve treatment provided a 1∶1 stoichiometry for the complex formation and an average Kd of ≈6 nM for the binding affinity. Samples were in phosphate buffer pH 6, I = 0.05, at 5°C, MgCl2 5 mM final concentration.

https://doi.org/10.1371/journal.pone.0040223.g005

The above results have several implications. First, the fact that RAL binds to LTR34, is consistent with an end fraying lowering the energy barrier for RAL accommodation into the terminal bases. In retroviral DNAs, the fraying of terminal base pairs seems amplified by the small 5′C4A3∶3′G−4T−3 duplex just before in the sequence. This dinucleotide duplex, invariant in retroviral DNAs, is considered as one of the less stacked and most malleable dinucleotide [62], [63]. It could contribute to increase the motions and disruption of the connected base pairs [35], [36], [62], [63], which according to several authors facilitate the IN binding and the scissile bond cleavage [35], [36]. Present study shows further that the end fraying in LTR may contributes to the capture of INSTIs, but the latter has no significant impact on the cut of the scissile bond, as INSTIs are generally weak 3′-P inhibitors. Second, the ability of RAL to bind both LTR32 and LTR34 and its inability to bind both the blunt ended LTR30 and the 3′ dangling ended LTR32-I, suggests that the terminal 5′A−1C−23′ dinucleotide is required for the drug-DNA complex formation, either as a dinucleotide overhang at the end of processed LTR or as a dinucleotide within a duplex subject to a large fraying at the end of unprocessed LTR. Actually, a fair amount of data has been reported on the stabilizing role of the 5′A−1C−23′ overhang in the complex of IN with LTR [34], [64]. Here, we understand that the overhang can be also involved in the LTR-RAL complex stability.

Molecular modelling

Our UV-absorption, CD and fluorescence experiments provide insight on the binding of RAL to the processed (LTR32) and unprocessed (LTR34) DNAs. However, they do not give information on the binding events and the type of interactions stabilizing molecular complexes. On the other hand, MD simulations can provide atomic details on structural and dynamic events governing the complex formation, while MMPBSA provides the binding free energies and allows the quantification of the complex stability [30], [51], [65]. We performed MD simulations (400 ns in total) and MMPBSA calculations in order to unravel the role of viral DNA ends as possible primary targets in the mechanism of IN inhibition by RAL. Indeed, MMPBSA overestimate the binding values, which is not surprising since this method is known for this defect [51], [66]. However, the MMPBSA values presented hereafter provide the same ranking as that given by fluorescence, in showing that the complexes of RAL with processed LTR are more stable than those with unprocessed LTR.

Equilibration of the MD simulations

We applied MD simulations to the analysis of RAL complexes with LTR34 and LTR32 ends. MD trajectories monitored by the root-mean-square displacement (RMSD) values of heavy atoms with respect to the X-ray structure (PDB code: 3OYA) are shown in Fig. 6A and B. Similar RMSD values were obtained when the sugar C4′ atoms (green curve in Figure 6A for LTR34) or the phosphorus atoms (not shown) were monitored. The trajectories show that the system is stabilized after 2 ns and conserves the same RMSD value till the end of the simulation (100 ns) for both LTR34 and LTR32. A duplicate simulation also of 100 ns gives similar results. LTR34 displays higher RMSD values than LTR32 (the purple curve in Fig. 6A is very similar to the black and the blue in Fig. 6B). The important fraying (base pair disruption or impairing) of the two terminal base pairs of LTR34, is the main reason of the greater flexibility of LTR34 in comparison with LTR32. This is particularly obvious in the root-mean-square fluctuation (RMSF) of the sugar C4′ atoms (see Fig. 6C and D).

thumbnail
Figure 6. MD simulations of the RAL-LTR34 and RAL-LTR32 complex systems (PFV oligonucleotides), using GROMACS with the AMBER force field.

(A) Time evolution of RMSD (root mean square deviation) values based on all the heavy atoms for the two LTR34 trajectories (black: LTR34-1 and blue: LTR34-2). RMSD calculations for a single trajectory were also performed using the sugar C4′ atoms (green: LTR34-1) and repeated for LTR34 devoid of 3′-AT (purple). (B) Time evolution of RMSD values of LTR32 for two trajectories (black: LTR32-1 and blue: LTR32-2). (C) RMSF (root mean square fluctuation) variations of sugar C4′ atoms for LTR34 and (D) RMSF variations of sugar C4′ atoms for LTR32.

https://doi.org/10.1371/journal.pone.0040223.g006

RAL binding to DNA

Each oligonucleotide yields two complexes: RAL-LTR34-1, RAL-LTR34-2 and RAL-LTR32-1, RAL-LTR32-2 (Fig. 7). Table 1 summarizes the free energy values for the binding of RAL to the LTR34 and LTR32 ends, calculated with the MMPBSA method. Values agree with a favorable binding of RAL to both unprocessed LTR34 (RAL-LTR34-1, RAL-LTR34-2) and processed LTR32 (RAL-LTR32-1, RAL-LTR32-2). The binding energy characterizing RAL-LTR34-1 is less favorable compared with RAL-LTR34-2, and also compared with the two other complexes, RAL-LTR32-1and LTR32-2. Compared with RAL-LTR34-2 (and the other two complexes), RAL-LTR34-1 displays a very distinct binding mode of RAL. In fact, it is the only case where the drug uses its oxadiazole moiety to bind DNA. The ring intercalates in between the terminal base pairs, while the remaining of the molecule is solvent exposed. In RAL-LTR34-2, RAL uses its fluorobenzyl moiety to interact with the G−4 and the T−3 bases, and its pyrimidine ring to interact with the A3 of the conserved C4pA3step. Interactions of the fluorobenzyl moiety with of G−4 and T−3 bases are also found in both RAL-LTR32-1 and RAL-LTR32-2. In RAL-LTR32-1 the fluorobenzyl group interacts with the C4 base, while in RAL-LTR32-2 the pyrimidine ring interacts with the A3 base.

thumbnail
Figure 7. Snapshots from the two 100 ns trajectories of RAL in complex with unprocessed LTR (LTR34-1 and 2, top) and processed LTR (LTR32-1 and 2, bottom).

RAL is colored in slime green and bases in sandy brown, except for atoms at interacting distances which are colored using the usual code (hydrogen in white; nitrogen in blue; and oxygen in red) except for carbons, while the Mg2+ ion is represented by magenta ball. Selected snapshots are 0 ns (the initial structure), 50 ns and 100 ns.

https://doi.org/10.1371/journal.pone.0040223.g007

thumbnail
Table 1. Calculated binding parameters for the complexes of RAL with LTR32 and LTR34.

https://doi.org/10.1371/journal.pone.0040223.t001

It is worth noting, that the terminal adenine A3 bearing the recessed 3′-hydroxyl, samples two different conformations in the calculated RAL-LTR32 complexes. These two conformations are found in the crystal structures of RAL bound to the PFV intasome (PDB codes: 3L2T [15] and 3OYA [18]). In the corresponding Fig. 8, the role of the fluorobenzyl ring containing moiety appears essential to stabilization of the RAL-DNA complex and consequently to induction of inhibition. In three complexes out of four provided by calculations, the ring partakes in key interactions. Actually, the stacking of the aromatic ring on the cytosine C4 base and the interactions of the RAL amide group with the adenine A3 sugar in the invariant 5′C4A33′step, as well as the interaction of the fluorine atom with the guanine G−4 at the back of the cavity, are all found in the crystal structures of INSTI-intasome [18]. Moreover, the fluorobenzyl moiety in RAL-LTR34-2 displays the same stabilizing interactions than in RAL-LTR32-1, or 2, and in the crystal structures of RAL-intasome [15].

thumbnail
Figure 8. Details of the interactions of RAL with its surrounding amino acids and nucleotides as observed in the two X-ray structures of RAL bound to the PFV intasome (pdb codes: 3L2T and 3OYA).

RAL is shown in slime green and IN and LTR residues in sandy brown. The amino acids and nucleotides giving interactions are shown in sticks with hydrogens in white, nitrogens in blue and oxygens in red. Mg2+ ions are represented by light green balls. The main difference between 3L2T and 3OYA structures concerns the orientation of the adenine A3.

https://doi.org/10.1371/journal.pone.0040223.g008

Complex stabilities

The calculated binding energies are shown in table 1. We note that when the fluorobenzyl moiety is not involved in interactions (RAL-LTR34-1), the binding energy is the highest, which confirms the key role of the halogenated moiety in the complex stabilization. In the four RAL-DNA complexes, the van der Waals (VDW) interactions (Evdw) are predominating which is a relatively common feature in complexes of organic ligands with nucleic acids [67]. In the two complexes of RAL with unprocessed LTR34, the fluorobenzyl ring has much more favorable VDW interactions compared with the oxadiazole group: ([LTR34-2]vdw−[LTR34-1]vdw = −18 kcal/mol). The solvent (ΔGPB+ΔGSASA) also contributes to a favorable binding, while both the electrostatics and the entropy are not favorable to the binding. The latter affects apparently equally the binding of RAL to LTR34 and to LTR32. Since the electrostatics and the solvation forces neutralize each other, the VDW interactions become the main binding contribution. Finally, the binding of RAL to viral DNA ends appears as an enthalpy driven process, rather than an entropy driven one.

Binding free energies: comparisons of simulated with experimental structures

The free energy, ΔG, values for the binding of RAL to LTR34 and LTR32, provided by the MMPBSA methodology are overestimated compared with the experimental ΔG values determined by fluorescence anisotropy (ΔGLTR32≈−10.5 Kcal/mol and ΔGLTR34≈−9.8 Kcal/mol). This result is not unexpected as MMPBSA is known for this defect, especially when a non-polarizable force field is used. Yet, the fluorescence titration experiments and the molecular dynamics approaches lead to similar conclusions. Both show that RAL can form stable complexes with LTR34 and LTR32. The MMPBSA calculations predict a binding of RAL to viral DNA which is enthalpy driven and assign an important role to the VDW forces in the complex stabilization. The finding of an enthalpy driven binding is not so surprising, as RAL, similarly to intercalators, also used π-π stackings to interact with the bases of DNA (Fig. 7). A dataset consisting of 26 binding interactions has shown that intercalating molecules bind to DNA with a favorable enthalpy contribution [68], while the binding of groove-binders is due more to a favorable entropy [68]. Actually, the binding of RAL to LTRs is characterized by a mean ΔH/ΔG ratio of about 1.50, while a ΔH/ΔG ratio ranging from 0.83 to 1.97 is considered as a clear signature of enthalpy driven binding [68]. Although RAL has a structure reminiscent of some intercalators, our hydrodynamic studies (not shown) indicate that it is unable to insert into base pairs of the DNA double helix. Indeed, the anchoring of RAL at the end of unprocessed LTR needs of the end fraying. In processed LTR the room opened at the DNA end by the release of the 5′G2T13′ (5′A2T13′ for IN-PFV) dinucleotide facilitates the RAL interaction. In both processed and unprocessed LTRs we find that the 5′A−1C−23′ (5′A−1T−23′ for PFV IN) dinucleotide partakes in the DNA-drug complex stabilization. Yet, in the crystal structure of the ternary complex of RAL with the PFV intasome, the 5′AT3′ overhang rather interacts with the protein [15] as it was already the case in the intasome [15]. Here, either the interactions of the 5′A−1T−23′ overhang could be stronger with the protein than with RAL, and be prioritized, or the crystal packing could prevent RAL from reaching its target. Noteworthy, RAL makes van der Waals contacts with both the invariant 5′C4pA33′step known for its functional importance at the LTR end [34], [61] and the guanine G−4 facing the cytosine C4 of this step. These interactions are also found in calculated and crystal structures of DTG, a structural analogue showing activity against RAL-resistant mutants [69]. In that case it has been suggested that DNA could make the greatest energetic contribution to DTG binding.

The particular properties of LTR ends contribute to the molecular fitting of RAL to DNA prior 3′-P. Owing to the high flexibility prevailing at the unprocessed DNA end the energy loss during complex formation is weak. This is especially true with regard to the drug intercalation which is an entropy costly mechanism requiring both base pair destacking and DNA distortion to create a suited site for the binding, and resulting in a damping of motions and stiffening of the duplex.

Conclusion

Our experimental and theoretical studies underline the particular role of the LTR terminal nucleotides in the binding of RAL to viral DNA. In previous reports, the antiviral drug RAL has been described as an INSTI acting at the interface of DNA-enzyme according to an important lattice of interactions with the processed LTR end, two metal ions and IN [3], [15], [18]. The recently published X-ray crystallography results describing the binding of INSTIs to the PFV-intasome, have confirmed most of the previous biochemical observations, including the two metal binding, and provided outstanding information on the inhibition mechanism used by inhibitors [15], [17], [18]. Yet, so far reported studies had not addressed the possible binding of IN inhibitors to the viral DNA end, especially prior to 3′-P. Actually, like other small organic molecules with chemotherapeutic activities, RAL binds directly and selectively to DNA. However, in contrast to the anticancer agents, as for example anthracyclines and ellipticines, RAL occupies selective binding sites on DNA. Its binding to the unprocessed LTR end could be at the basis of the observed 3′-P inhibition. What appears is that the capture of the drug by unprocessed LTR entails a small damping of motions (fraying) in the terminal base pairs (Fig. 9). Initial studies [35], [36], [70][73] have shown that any restriction of fraying in the terminal base pairs by either extension of the duplex [35], [71] or chemical linkage of the duplex ends [73] impairs the 3′-P reaction. Remarkably, RAL keeps the same conformation within the complex with unprocessed and processed DNAs, with its halogenated ring in face to face contact with the cytosine base of the conserved 5′C4pA33′ step, and the adenine A3 bearing the recessed 3′-hydroxyl group moved from its operative position.

thumbnail
Figure 9. Effects of RAL on the distance and the fraying of unprocessed LTR ends.

Time evolution of the distance (nm) between the mass of RAL (D) and that of the terminal bases (T and A) of LTR34-1 (black:D-T and red: D-A) and LTR34-2 (purple: D-T and blue: D-A). Interaction of RAL with the terminal bases decreases the distance between the drug and the bases and also reduces moderately the end fraying.

https://doi.org/10.1371/journal.pone.0040223.g009

At the end our work could help to better understand some of the factors contributing to the mechanism of action of IN inhibitors. Above all, we have shown that RAL binds to the LTR end prior 3′-P reaction, but the impact on this reaction is small. The fact that the contacts established by RAL with the processed LTR end mainly concern the nucleotides C4, A3 and G−4, pushes us to speculate that new drugs giving an increased number of interactions with these highly conserved bases, at expense of interactions with amino acid side chains of the protein active site, will be better INSTIs and weaker resistance inducers.

Supporting Information

Figure S1.

UV-absorption analysis of oligonucleotides and raltegravir free and in complexes. (A) Spectra of RAL (80 µM, in red) and LTR32 (10 µM, in black) in phosphate buffer pH 6, I = 0.05, MgCl2 5 mM final concentration. (B) Spectra of RAL 20 µM, (black), in complex with LTR32 (1 µM, red), LTR32 (5 µM, blue) and LTR32 (10 µM, green).

https://doi.org/10.1371/journal.pone.0040223.s001

(TIF)

Figure S2.

Binding of oligonucleotides to raltegravir. Titration data for LTR32 (black) and LTR34 (red) and corresponding Kds are obtained from fluorescence intensity in so called reverse experiments. The spectra of raltegravir recorded at different LTR34 concentrations are given in insert.

https://doi.org/10.1371/journal.pone.0040223.s002

(TIF)

Acknowledgments

We thank Dr. M. Buckle for support of these studies.

Author Contributions

Conceived and designed the experiments: SF RGM. Performed the experiments: FFA SA LZ. Analyzed the data: FFA SA LZ ZH RGM SF. Contributed reagents/materials/analysis tools: LZ. Wrote the paper: SF FFA SA RGM.

References

  1. 1. Goff SP (1992) Genetics of Retroviral Integration. Annual Review of Genetics 26: 527–544.
  2. 2. Katz RA, Skalka AM (1994) The retroviral enzymes. Annu Rev Biochem 63: 133–173. 10.1146/annurev.bi.63.070194.001025 [doi].
  3. 3. Pommier Y, Johnson AA, Marchand C (2005) Integrase inhibitors to treat HIV/AIDS. Nat Rev Drug Discov 4: 236–248. nrd1660 [pii];10.1038/nrd1660 [doi].
  4. 4. Li X, Krishnan L, Cherepanov P, Engelman A (2011) Structural biology of retroviral DNA integration. Virology 411: 194–205. S0042-6822(10)00752-X [pii];10.1016/j.virol.2010.12.008 [doi].
  5. 5. Grobler JA, Stillmock K, Hu B, Witmer M, Felock P, et al. (2002) Diketo acid inhibitor mechanism and HIV-1 integrase: implications for metal binding in the active site of phosphotransferase enzymes. Proc Natl Acad Sci U S A 99: 6661–6666. 10.1073/pnas.092056199 [doi];092056199 [pii].
  6. 6. Marchand C, Maddali K, Metifiot M, Pommier Y (2009) HIV-1 IN inhibitors: 2010 update and perspectives. Curr Top Med Chem 9: 1016–1037. Abs-010-9-11 [pii].
  7. 7. McColl DJ, Chen X (2010) Strand transfer inhibitors of HIV-1 integrase: bringing IN a new era of antiretroviral therapy. Antiviral Res 85: 101–118. S0166-3542(09)00533-6 [pii];10.1016/j.antiviral.2009.11.004 [doi].
  8. 8. Grinsztejn B, Nguyen BY, Katlama C, Gatell JM, Lazzarin A, et al. (2007) Safety and efficacy of the HIV-1 integrase inhibitor raltegravir (MK-0518) in treatment-experienced patients with multidrug-resistant virus: a phase II randomised controlled trial. Lancet 369: 1261–1269. S0140-6736(07)60597-2 [pii];10.1016/S0140-6736(07)60597-2 [doi].
  9. 9. Koelsch KK, Cooper DA (2009) Integrase inhibitors in salvage therapy regimens for HIV-1 infection. Curr Opin HIV AIDS 4: 518–523. 10.1097/COH.0b013e328331b526 [doi];01222929-200911000-00011 [pii].
  10. 10. Marinello J, Marchand C, Mott BT, Bain A, Thomas CJ, et al. (2008) Comparison of raltegravir and elvitegravir on HIV-1 integrase catalytic reactions and on a series of drug-resistant integrase mutants. Biochemistry 47: 9345–9354. 10.1021/bi800791q [doi].
  11. 11. Metifiot M, Marchand C, Maddali K, Pommier Y (2010) Resistance to integrase inhibitors. Viruses 2: 1347–1366.
  12. 12. Reigadas S, Anies G, Masquelier B, Calmels C, Stuyver LJ, et al. (2010) The HIV-1 integrase mutations Y143C/R are an alternative pathway for resistance to Raltegravir and impact the enzyme functions. PLoS One 5: e10311. 10.1371/journal.pone.0010311 [doi].
  13. 13. Zargarian L, Benleumi MS, Renisio JG, Merad H, Maroun RG, et al. (2003) Strategy to discriminate between high and low affinity bindings of human immunodeficiency virus, type 1 integrase to viral DNA. J Biol Chem 278: 19966–19973.
  14. 14. Hobaika Z, Zargarian L, Maroun RG, Mauffret O, Burke TR, et al. (2010) HIV-1 integrase and virus and cell DNAs: complex formation and perturbation by inhibitors of integration. Neurochem Res 35: 888–893. 10.1007/s11064-009-0098-2 [doi].
  15. 15. Hare S, Gupta SS, Valkov E, Engelman A, Cherepanov P (2010) Retroviral intasome assembly and inhibition of DNA strand transfer. Nature 464: 232–236. nature08784 [pii];10.1038/nature08784 [doi].
  16. 16. Goldgur Y, Craigie R, Cohen GH, Fujiwara T, Yoshinaga T, et al. (1999) Structure of the HIV-1 integrase catalytic domain complexed with an inhibitor: a platform for antiviral drug design. Proc Natl Acad Sci U S A 96: 13040–13043.
  17. 17. Krishnan L, Li X, Naraharisetty HL, Hare S, Cherepanov P, et al. (2010) Structure-based modeling of the functional HIV-1 intasome and its inhibition. Proc Natl Acad Sci U S A 107: 15910–15915. 1002346107 [pii];10.1073/pnas.1002346107 [doi].
  18. 18. Hare S, Vos AM, Clayton RF, Thuring JW, Cummings MD, et al. (2010) Molecular mechanisms of retroviral integrase inhibition and the evolution of viral resistance. Proc Natl Acad Sci U S A 107: 20057–20062. 1010246107 [pii];10.1073/pnas.1010246107 [doi].
  19. 19. Sato M, Motomura T, Aramaki H, Matsuda T, Yamashita M, et al. (2006) Novel HIV-1 integrase inhibitors derived from quinolone antibiotics. J Med Chem 49: 1506–1508. 10.1021/jm0600139 [doi].
  20. 20. Da Silva D, Van WL, Breilh D, Reigadas S, Anies G, et al. (2010) HIV-1 resistance patterns to integrase inhibitors in antiretroviral-experienced patients with virological failure on raltegravir-containing regimens. J Antimicrob Chemother 65: 1262–1269. dkq099 [pii];10.1093/jac/dkq099 [doi].
  21. 21. Metifiot M, Johnson B, Smith S, Zhao XZ, Marchand C, et al. (2011) MK-0536 inhibits HIV-1 integrases resistant to raltegravir. Antimicrob Agents Chemother 55: 5127–5133. AAC.05288-11 [pii];10.1128/AAC.05288-11 [doi].
  22. 22. Vacca J, Wai J, Fisher T, Embrey M, Hazuda D, et al. (2007) Discovery of MK-2048 – subtle changes confer unique resistance properties to a series of tricyclic hydroxypyrrole integrase strand transfer inhibitors. Poster discussion: 4th IAS Conference on HIV Pathogenesis, Treatment and Prevention: Abstract no WEPEA088.
  23. 23. Johns B, Kawasuiji T, Taishi T, Yoshida H, Garvey E, et al. (2010) The discovery of S/GSK1349572: a once-daily next generation integrase inhibitor with a superior resistance profile. in 17th Conference on Retroviruses and Opportunistic Infections, SanFrancisco, CA.
  24. 24. Hare S, Smith SJ, Metifiot M, Jaxa-Chamiec A, Pommier Y, et al. (2011) Structural and functional analyses of the second-generation integrase strand transfer inhibitor dolutegravir (S/GSK1349572). Mol Pharmacol 80: 565–572. mol.111.073189 [pii];10.1124/mol.111.073189 [doi].
  25. 25. Espeseth AS, Felock P, Wolfe A, Witmer M, Grobler J, et al. (2000) HIV-1 integrase inhibitors that compete with the target DNA substrate define a unique strand transfer conformation for integrase. Proc Natl Acad Sci U S A 97: 11244–11249. 10.1073/pnas.200139397 [doi];200139397 [pii].
  26. 26. Fitzkee NC, Masse JE, Shen Y, Davies DR, Bax A (2010) Solution conformation and dynamics of the HIV-1 integrase core domain. J Biol Chem 285: 18072–18084. M110.113407 [pii];10.1074/jbc.M110.113407 [doi].
  27. 27. Monnot M, Mauffret O, Simon V, Lescot E, Psaume B, et al. (1991) DNA-drug recognition and effects on topoisomerase II-mediated cytotoxicity. A three-mode binding model for ellipticine derivatives. J Biol Chem 266: 1820–1829.
  28. 28. Hess B, Kutzner C, Van Der Spoel D, Lindahl E (2008) GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J Chem Theory Comput 4: 435–447.
  29. 29. Lindorff-Larsen K, Piana S, Palmo K, Maragakis P, Klepeis JL, et al. (2010) Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins 78: 1950–1958. 10.1002/prot.22711 [doi].
  30. 30. Srinivasan J, Cheatham TE, Cieplak P, Kollman PA, Case DA (1998) Continuum Solvent Studies of the Stability of DNA, RNA, and Phosphoramidate-DNA Helices. J Am Chem Soc 120: 9401–9409.
  31. 31. Massova I, Kollman PA (1999) Computational Alanine Scanning To Probe Protein-Protein Interactions: A Novel Approach To Evaluate Binding Free Energies. J Am Chem Soc 121: 8133–8143.
  32. 32. Hou T, Wang J, Li Y, Wang W (2011) Assessing the performance of the MM/PBSA and MM/GBSA methods. 1. The accuracy of binding free energy calculations based on molecular dynamics simulations. J Chem Inf Model 51: 69–82.
  33. 33. Xue W, Liu H, Yao X (2012) Molecular mechanism of HIV-1 integrase-vDNA interactions and strand transfer inhibitor action: A molecular modeling perspective. J Comput Chem 33: 527–536. 10.1002/jcc.22887 [doi].
  34. 34. Dicker IB, Samanta HK, Li Z, Hong Y, Tian Y, et al. (2007) Changes to the HIV long terminal repeat and to HIV integrase differentially impact HIV integrase assembly, activity, and the binding of strand transfer inhibitors. J Biol Chem 282: 31186–31196. M704935200 [pii];10.1074/jbc.M704935200 [doi].
  35. 35. Scottoline BP, Chow S, Ellison V, Brown PO (1997) Disruption of the terminal base pairs of retroviral DNA during integration. Genes Dev 11: 371–382.
  36. 36. Katz RA, Merkel G, Andrake MD, Roder H, Skalka AM (2011) Retroviral integrases promote fraying of viral DNA ends. J Biol Chem 286: 25710–25718. M111.229179 [pii];10.1074/jbc.M111.229179 [doi].
  37. 37. Bujacz G, Jaskolski M, Alexandratos J, Wlodawer A, Merkel G, et al. (1996) The catalytic domain of avian sarcoma virus integrase: conformation of the active-site residues in the presence of divalent cations. Structure 4: 89–96.
  38. 38. Liu LF (1989) DNA topoisomerase poisons as antitumor drugs. Annu Rev Biochem 58: 351–375. 10.1146/annurev.bi.58.070189.002031 [doi].
  39. 39. Pommier Y, Leo E, Zhang H, Marchand C (2010) DNA topoisomerases and their poisoning by anticancer and antibacterial drugs. Chem Biol 17: 421–433. S1074-5521(10)00161-4 [pii];10.1016/j.chembiol.2010.04.012 [doi].
  40. 40. Bujalowski W, Lohman TM (1987) A general method of analysis of ligand-macromolecule equilibria using a spectroscopic signal from the ligand to monitor binding. Application to Escherichia coli single-strand binding protein-nucleic acid interactions. Biochemistry 26: 3099–3106.
  41. 41. Wang J, Wolf RM, Caldwell JW, Kollman PA, Case DA (2004) Development and testing of a general amber force field. J Comput Chem 25: 1157–1174. 10.1002/jcc.20035 [doi].
  42. 42. Wang J, Wang W, Kollman PA, Case DA (2006) Automatic atom type and bond type perception in molecular mechanical calculations. J Mol Graph Model 25: 247–260. S1093-3263(05)00173-7 [pii];10.1016/j.jmgm.2005.12.005 [doi].
  43. 43. Roothaan CCJ (1951) New Developments in Molecular Orbital Theory. Rev Mod Phys 23: 69–89.
  44. 44. Rassolov VA, Ratner MA, Pople JA, Redfern PC, Curtiss LA (2001) 6-31G* Basis Set for Third-Row Atoms. J Comput Chem 22: 976–984.
  45. 45. Humphrey W, Dalke A, Schulten K (1996) VMD: visual molecular dynamics. J Mol Graph 14: 33–38. 0263785596000185 [pii].
  46. 46. Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML (1983) Comparison of Simple Potential Functions for Simulating Liquid Water. J Chem Phys 79: 926–935.
  47. 47. Darden T, Darrin Y, Pedersen L (1993) Particle mesh Ewald: An N. log(N) method for Ewald sums in large systems. J Chem Phys 98: 10089–10092.
  48. 48. Hess B, Bekker H, Berendsen HJ, Fraaije JG (1997) LINCS: A linear constraint solver for molecular simulations. J Comput Chem 18: 1463–1472.
  49. 49. Bussi G, Donadio D, Parrinello M (2007) Canonical sampling through velocity rescaling. J Chem Phys 126: 014101–014107. 10.1063/1.2408420 [doi].
  50. 50. Parrinello M, Rahman A (1981) Polymorphic transitions in single crystals: A new molecular dynamics method. J App Phys 52: 7182–7190.
  51. 51. Kollman PA, Massova I, Reyes C, Kuhn B, Huo S, et al. (2000) Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. Acc Chem Res 33: 889–897. ar000033j [pii].
  52. 52. Hurley LH (2002) DNA and its associated processes as targets for cancer therapy. Nat Rev Cancer 2: 188–200. 10.1038/nrc749 [doi].
  53. 53. Dougherty G, Pigram WJ (1982) Spectroscopic analysis of drug-nucleic acid interactions. CRC Crit Rev Biochem 12: 103–132.
  54. 54. Wilson WD, Tanious FA, Barton HJ, Wydra RL, Jones RL, et al. (1990) The interaction of unfused polyaromatic heterocycles with DNA: intercalation, groove-binding and bleomycin amplification. Anticancer Drug Des 5: 31–42.
  55. 55. Hertzberg RP, Caranfa MJ, Hecht SM (1989) On the mechanism of topoisomerase I inhibition by camptothecin: evidence for binding to an enzyme-DNA complex. Biochemistry 28: 4629–4638.
  56. 56. Fasman GD (1996) Circular dichroism and the conformationnal analysis of biomolecules. New York: Plenum Press.
  57. 57. Berova N, Nakanishi K, Woody RW (2000) Circular dichroism: principles and applications. Wiley-VCH, New York.
  58. 58. Monnot M, Mauffret O, Lescot E, Fermandjian S (1992) Probing intercalation and conformational effects of the anticancer drug 2-methyl-9-hydroxyellipticinium acetate in DNA fragments with circular dichroism. Eur J Biochem 204: 1035–1039.
  59. 59. Heyduk T, Lee JC (1990) Application of fluorescence energy transfer and polarization to monitor Escherichia coli cAMP receptor protein and lac promoter interaction. Proc Natl Acad Sci U S A 87: 1744–1748.
  60. 60. Hill JJ, Royer CA (1997) Fluorescence approaches to study of protein-nucleic acid complexation. Methods Enzymol 278: 390–416.
  61. 61. Langley DR, Samanta HK, Lin Z, Walker MA, Krystal MR, et al. (2008) The terminal (catalytic) adenosine of the HIV LTR controls the kinetics of binding and dissociation of HIV integrase strand transfer inhibitors. Biochemistry 47: 13481–13488. 10.1021/bi801372d [doi].
  62. 62. Renisio JG, Cosquer S, Cherrak I, El AS, Mauffret O, et al. (2005) Pre-organized structure of viral DNA at the binding-processing site of HIV-1 integrase. Nucleic Acids Res 33: 1970–1981. 33/6/1970 [pii];10.1093/nar/gki346 [doi].
  63. 63. Calladine CR, Drew HR (1992) Understanding DNA: The Molecule and How It Works. Academic Press, London.
  64. 64. Ellison V, Brown PO (1994) A stable complex between integrase and viral DNA ends mediates human immunodeficiency virus integration in vitro. Proc Natl Acad Sci U S A 91: 7316–7320.
  65. 65. Harris SA, Gavathiotis E, Searle MS, Orozco M, Laughton CA (2001) Cooperativity in drug-DNA recognition: a molecular dynamics study. J Am Chem Soc 123: 12658–12663. ja016233n [pii].
  66. 66. Treesuwan W, Wittayanarakul K, Anthony NG, Huchet G, Alniss H, et al. (2009) A detailed binding free energy study of 2∶1 ligand-DNA complex formation by experiment and simulation. Phys Chem Chem Phys 11: 10682–10693. 10.1039/b910574c [doi].
  67. 67. Smith GF (2009) Medicinal chemistry by the numbers: the physicochemistry, thermodynamics and kinetics of modern drug design. Prog Med Chem 48: 1–29.
  68. 68. Chaires JB (2006) A thermodynamic signature for drug-DNA binding mode. Arch Biochem Biophys 453: 26–31. S0003-9861(06)00132-9 [pii];10.1016/j.abb.2006.03.027 [doi].
  69. 69. Johnson BC, Metifiot M, Pommier Y, Hughes SH (2012) Molecular dynamics approaches estimate the binding energy of HIV-1 integrase inhibitors and correlate with in vitro activity. Antimicrob Agents Chemother 56: 411–419. AAC.05292-11 [pii];10.1128/AAC.05292-11 [doi].
  70. 70. Vink C, van Gent DC, Elgersma Y, Plasterk RH (1991) Human immunodeficiency virus integrase protein requires a subterminal position of its viral DNA recognition sequence for efficient cleavage. J Virol 65: 4636–4644.
  71. 71. Engelman A, Mizuuchi K, Craigie R (1991) HIV-1 DNA integration: mechanism of viral DNA cleavage and DNA strand transfer. Cell 67: 1211–1221. 0092-8674(91)90297-C [pii].
  72. 72. Vink C, Yeheskiely E, van der Marel GA, van Boom JH, Plasterk RH (1991) Site-specific hydrolysis and alcoholysis of human immunodeficiency virus DNA termini mediated by the viral integrase protein. Nucleic Acids Res 19: 6691–6698.
  73. 73. Agapkina J, Smolov M, Zubin E, Mouscadet JF, Gottikh M (2004) HIV-1 integrase can process a 3′-end crosslinked substrate. Eur J Biochem 271: 205–211.