Pre-organized structure of viral DNA at the binding-processing site of HIV-1 integrase

Development of Force Field Parameters for the Simulation of Single- and Double-Stranded DNA Molecules and DNA-Protein Complexes

Although molecular dynamics (MD) simulations have been used extensively to study the structural dynamics of proteins, the role of MD simulation in studies of nucleic acid based systems has been more limited. One contributing factor to this disparity is the historically lower level of accuracy of the physical models used in such simulations to describe interactions involving nucleic acids. By modifying nonbonded and torsion parameters of a force field from the Amber family of models, we recently developed force field parameters for RNA that achieve a level of accuracy comparable to that of state-of-the-art protein force fields. Here we report force field parameters for DNA, which we developed by transferring nonbonded parameters from our recently reported RNA force field and making subsequent adjustments to torsion parameters. We have also modified the backbone charges in both the RNA and DNA parameter sets to make the treatment of electrostatics compatible with our recently developed variant of the Amber protein and ion force field. We name the force field resulting from the union of these three parameter sets (the new DNA parameters, the revised RNA parameters, and the existing protein and ion parameters) DES-Amber. Extensive testing of DES-Amber indicates that it can describe the thermal stability and conformational flexibility of single- and double-stranded DNA systems with a level of accuracy comparable to or, especially for disordered systems, exceeding that of state-of-the-art nucleic acid force fields. Finally, we show that, in certain favorable cases, DES-Amber can be used for long-timescale simulations of protein–nucleic acid complexes.

A New Approach to 3D Modeling of Inhomogeneous Populations of Viral Regulatory RNA

Tertiary structure (3D) is the physical context of RNA regulatory activity. Retroviruses are RNA viruses that replicate through the proviral DNA intermediate transcribed by hosts. Proviral transcripts form inhomogeneous populations due to variable structural ensembles of overlapping regulatory RNA motifs in the 5′-untranslated region (UTR), which drive RNAs to be spliced or translated, and/or dimerized and packaged into virions. Genetic studies and structural techniques have provided fundamental input constraints to begin predicting HIV 3D conformations in silico. Using SimRNA and sets of experimentally-determined input constraints of HIV^NL4-3 trans-activation responsive sequence (TAR) and pairings of unique-5′ (U5) with dimerization (DIS) or AUG motifs, we calculated a series of 3D models that differ in proximity of 5′-Cap and the junction of TAR and PolyA helices; configuration of primer binding site (PBS)-segment; and two host cofactors binding sites. Input constraints on U5-AUG pairings were most compatible with intramolecular folding of 5′-UTR motifs in energetic minima. Introducing theoretical constraints predicted metastable PolyA region drives orientation of 5′-Cap with TAR, U5 and PBS-segment helices. SimRNA and the workflow developed herein provides viable options to predict 3D conformations of inhomogeneous populations of large RNAs that have been intractable to conventional ensemble methods.

Selectivity for strand-transfer over 3'-processing and susceptibility to clinical resistance of HIV-1 integrase inhibitors are driven by key enzyme-DNA interactions in the active site

Integrase strand transfer inhibitors (INSTIs) are highly effective against HIV infections. Co-crystal structures of the prototype foamy virus intasome have shown that all three FDA-approved drugs, raltegravir (RAL), elvitegravir and dolutegravir (DTG), act as interfacial inhibitors during the strand transfer (ST) integration step. However, these structures give only a partial sense for the limited inhibition of the 3′-processing reaction by INSTIs and how INSTIs can be modified to overcome drug resistance, notably against the G140S-Q148H double mutation. Based on biochemical experiments with modified oligonucleotides, we demonstrate that both the viral DNA +1 and −1 bases, which flank the 3′-processing site, play a critical role for 3′-processing efficiency and inhibition by RAL and DTG. In addition, the G140S-Q148H (SH) mutant integrase, which has a reduced 3′-processing activity, becomes more active and more resistant to inhibition of 3′-processing by RAL and DTG in the absence of the −1 and +1 bases. Molecular modeling of HIV-1 integrase, together with biochemical data, indicate that the conserved residue Q146 in the flexible loop of HIV-1 integrase is critical for productive viral DNA binding through specific contacts with the virus DNA ends in the 3′-processing and ST reactions. The potency of integrase inhibitors against 3′-processing and their ability to overcome resistance is discussed.

A targeted DNA substrate mechanism for the inhibition of HIV-1 integrase by inhibitors with antiretroviral activity

We recently reported that viral DNA could be the primary target of raltegravir (RAL), an efficient anti‐HIV‐1 drug, which acts by inhibiting integrase. To elucidate this mechanism, we conducted a comparative analysis of RAL and TB11, a diketoacid abandoned as an anti‐HIV‐1 drug for its weak efficiency and marked toxicity, and tested the effects of the catalytic cofactor Mg²⁺ (5 mm) on drug‐binding properties. We used circular dichroism and fluorescence to determine drug affinities for viral DNA long terminal repeats (LTRs) and peptides derived from the integrase active site and DNA retardation assays to assess drug intercalation into DNA base pairs. We found that RAL bound more tightly to LTR ends than did TB11 (a diketo acid bearing an azido group) and that Mg²⁺ significantly increased the affinity of both RAL and TB11. We also observed a good relationship between drug binding with processed LTR and strand transfer inhibition. This unusual type of inhibition was caused by Mg²⁺‐assisted binding of drugs to DNA substrate, rather than to enzyme. Notably, while RAL bound exclusively to the cleavable/cleaved site, TB11 further intercalated into DNA base pairs and interacted with the integrase‐derived peptides. These unwanted binding sites explain the weaker bioavailability and higher toxicity of TB11 compared with the more effective RAL.

The intrinsic mechanics of B-DNA in solution characterized by NMR

Experimental characterization of the structural couplings in free B-DNA in solution has been elusive, because of subtle effects that are challenging to tackle. Here, the exploitation of the NMR measurements collected on four dodecamers containing a substantial set of dinucleotide sequences provides new, consistent correlations revealing the DNA intrinsic mechanics. The difference between two successive residual dipolar couplings (ΔRDCs) involving C6/8-H6/8, C3′-H3′ and C4′-H4′ vectors are correlated to the ³¹P chemical shifts (δP), which reflect the populations of the BI and BII backbone states. The δPs are also correlated to the internucleotide distances (D_inter) involving H6/8, H2′ and H2″ protons. Calculations of NMR quantities on high resolution X-ray structures and controlled models of DNA enable to interpret these couplings: the studied ΔRDCs depend mostly on roll, while D_inter are mainly sensitive to twist or slide. Overall, these relations demonstrate how δP measurements inform on key inter base parameters, in addition to probe the BI↔BII backbone equilibrium, and shed new light into coordinated motions of phosphate groups and bases in free B-DNA in solution. Inspection of the 5′ and 3′ ends of the dodecamers also supplies new information on the fraying events, otherwise neglected.

Unprocessed viral DNA could be the primary target of the HIV-1 integrase inhibitor raltegravir

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'C(4)pA(3)3' 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.

Localization of ASV integrase-DNA contacts by site-directed crosslinking and their structural analysis

Background

We applied crosslinking techniques as a first step in preparation of stable avian sarcoma virus (ASV) integrase (IN)-DNA complexes for crystallographic investigations. These results were then compared with the crystal structures of the prototype foamy virus (PFV) intasome and with published data for other retroviral IN proteins.

Methodology/Results

Photoaffinity crosslinking and site-directed chemical crosslinking were used to localize the sites of contacts with DNA substrates on the surface of ASV IN. Sulfhydryl groups of cysteines engineered into ASV IN and amino-modified nucleotides in DNA substrates were used for attachment of photocrosslinkers. Analysis of photocrosslinking data revealed several specific DNA-protein contacts. To confirm contact sites, thiol-modified nucleotides were introduced into oligo-DNA substrates at suggested points of contact and chemically crosslinked to the cysteines via formation of disulfide bridges. Cysteines incorporated in positions 124 and 146 in the ASV IN core domain were shown to interact directly with host and viral portions of the Y-mer DNA substrate, respectively. Crosslinking of an R244C ASV IN derivative identified contacts at positions 11 and 12 on both strands of viral DNA. The most efficient disulfide crosslinking was observed for complexes of the ASV IN E157C and D64C derivatives with linear viral DNA substrate carrying a thiol-modified scissile phosphate.

Conclusion

Analysis of our crosslinking results as well as published results of retroviral IN protein from other laboratories shows good agreement with the structure of PFV IN and derived ASV, HIV, and MuLV models for the core domain, but only partial agreement for the N- and C-terminal domains. These differences might be explained by structural variations and evolutionary selection for residues at alternate positions to perform analogous functions, and by methodological differences: i.e., a static picture of a particular assembly from crystallography vs. a variety of interactions that might occur during formation of functional IN complexes in solution.

Specificity of LTR DNA recognition by a peptide mimicking the HIV-1 integrase {alpha}4 helix

HIV-1 integrase integrates retroviral DNA through 3′-processing and strand transfer reactions in the presence of a divalent cation (Mg²⁺ or Mn²⁺). The α4 helix exposed at the catalytic core surface is essential to the specific recognition of viral DNA. To define group determinants of recognition, we used a model composed of a peptide analogue of the α4 helix, oligonucleotides mimicking processed and unprocessed U5 LTR end and 5 mM Mg²⁺. Circular dichroism, fluorescence and NMR experiments confirmed the implication of the α4 helix polar/charged face in specific and non-specific bindings to LTR ends. The specific binding requires unprocessed LTR ends—i.e. an unaltered 3′-processing site CA↓GT3′—and is reinforced by Mg²⁺ (K_d decreases from 2 to 0.8 nM). The latter likely interacts with the ApG and GpT3′ steps of the 3′-processing site. With deletion of GT3′, only persists non-specific binding (K_d of 100 μM). Proton chemical shift deviations showed that specific binding need conserved amino acids in the α4 helix and conserved nucleotide bases and backbone groups at LTR ends. We suggest a conserved recognition mechanism based on both direct and indirect readout and which is subject to evolutionary pressure.

HIV-1 integrase and virus and cell DNAs: complex formation and perturbation by inhibitors of integration

HIV-1 integrase (IN) catalyzes integration of viral DNA into cell DNA through 3'-processing of viral DNA and strand transfer reactions. To learn on binding of IN to DNAs and IN inhibition we applied spectroscopy (circular dichroism, fluorescence) in a simplified model consisting in a peptide analogue (K156) of alpha4 helix involved in recognition of viral and cell DNA; an oligonucleotide corresponding to the U5' LTR DNA end; and an inhibitor (TB11) of the diketo acid (DKA) family. Results extrapolated to IN show that: the enzyme binds viral DNA with high affinity and specificity, but cell DNA with low affinity and specificity; the affinity of TB11 for IN is high enough to impair the binding of IN to cell DNA, but not to viral DNA. This explains why TB11 is an inhibitor of strand transfer but not of 3'-processing. These results can help in the search of new IN inhibitors.

HIV-1 Integrase-DNA Recognition Mechanisms

Integration of a reverse transcribed DNA copy of the HIV viral genome into the host chromosome is essential for virus replication. This process is catalyzed by the virally encoded protein integrase. The catalytic activities, which involve DNA cutting and joining steps, have been recapitulated in vitro using recombinant integrase and synthetic DNA substrates. Biochemical and biophysical studies of these model reactions have been pivotal in advancing our understanding of mechanistic details for how IN interacts with viral and target DNAs, and are the focus of the present review.

Structural and dynamic characterization of the upper part of the HIV-1 cTAR DNA hairpin

First strand transfer is essential for HIV-1 reverse transcription. During this step, the TAR RNA hairpin anneals to the cTAR DNA hairpin; this annealing reaction is promoted by the nucleocapsid protein and involves an initial loop–loop interaction between the apical loops of TAR and cTAR. Using NMR and probing methods, we investigated the structural and dynamic properties of the top half of the cTAR DNA (mini-cTAR). We show that the upper stem located between the apical and the internal loops is stable, but that the lower stem of mini-cTAR is unstable. The residues of the internal loop undergo slow motions at the NMR time-scale that are consistent with conformational exchange phenomena. In contrast, residues of the apical loop undergo fast motions. The lower stem is destabilized by the slow interconversion processes in the internal loop, and thus the internal loop is responsible for asymmetric destabilization of mini-cTAR. These findings are consistent with the functions of cTAR in first strand transfer: its apical loop is suitably exposed to interact with the apical loop of TAR RNA and its lower stem is significantly destabilized to facilitate the subsequent action of the nucleocapsid protein which promotes the annealing reaction.

Changes to the HIV long terminal repeat and to HIV integrase differentially impact HIV integrase assembly, activity, and the binding of strand transfer inhibitors

Human immunodeficiency virus (HIV) integrase enzyme is required for the integration of viral DNA into the host cell chromosome. Integrase complex assembly and subsequent strand transfer catalysis are mediated by specific interactions between integrase and bases at the end of the viral long terminal repeat (LTR). The strand transfer reaction can be blocked by the action of small molecule inhibitors, thought to bind in the vicinity of the viral LTR termini. This study examines the contributions of the terminal four bases of the nonprocessed strand (G(2)T(1)C(-1)A(-2)) of the HIV LTR on complex assembly, specific strand transfer activity, and inhibitor binding. Base substitutions and abasic replacements at the LTR terminus provided a means to probe the importance of each nucleotide on the different functions. An approach is described wherein the specific strand transfer activity for each integrase/LTR variant is derived by normalizing strand transfer activity to the concentration of active sites. The key findings of this study are as follows. 1) The G(2):C(2) base pair is necessary for efficient assembly of the complex and for maintenance of an active site architecture, which has high affinity for strand transfer inhibitors. 2) Inhibitor-resistant enzymes exhibit greatly increased sensitivity to LTR changes. 3) The strand transfer and inhibitor binding defects of a Q148R mutant are due to a decreased affinity of the complex for magnesium. 4) Gln(148) interacts with G(2), T(1), and C(-1) at the 5' end of the viral LTR, with these four determinants playing important and overlapping roles in assembly, strand transfer catalysis and high affinity inhibitor binding.

NMR assessment of the global shape of a non-labelled DNA dodecamer containing a tandem of G-T mismatches

We have carried out a solution study of two non-labelled self-complementary DNA dodecamers d(GACTGTACAGTC)2 and d(GACTGTGCAGTC)2 by NMR, the second sequence composed of two G-T mismatches. Structures were determined using distances extracted from NOE effects alone or using both NOE and RDC constraints, measured in three different liquid crystalline media. We ensured that our data on the influence of the mesogen on the DNA structures, and the way in which the RDCs were incorporated as constraints in the protocol refinement, were consistent. We also tested the influence of different sets of RDCs and the best means of optimizing the calculation of D(a) and R. Resolution and accuracy of the ten best energy final structures were compared. The addition of a small set of RDC constraints significantly improves the final determined structures. We took advantage of the specificity of the RDC, i.e. it contains orientational information, and explored the global shape of the DNA duplexes; it was found that the duplexes do not have a large curvature. For the G-T base pair, we observed, in this particular sequence (tandem of G-T mismatches), a new pattern of base pairing, which involved the formation of a bifurcated hydrogen bond.

Sequence-dependent gating of an ion channel by DNA hairpin molecules

DNA hairpins produce ionic current signatures when captured by the alpha-hemolysin nano-scale pore under conditions of single molecule electrophoresis. Gating patterns produced by individual DNA hairpins when captured can be used to distinguish differences of a single base pair or even a single nucleotide [Vercoutere,W.A. et al. (2003) Nucleic Acids Res., 31, 1311–1318]. Here we investigate the mechanism(s) that may account for the ionic current gating signatures. The ionic current resistance profile of conductance states produced by DNA hairpin molecules with 3–12 bp stems showed a plateau in resistance between 10 and 12 bp, suggesting that hairpins with 10–12 bp stems span the pore vestibule. DNA hairpins with 9–12 bp stems produced gating signatures with the same relative conductance states. Systematic comparison of the conductance state dwell times and apparent activation energies for a series of 9–10 bp DNA hairpins suggest that the 3′ and 5′ ends interact at or near the limiting aperture within the vestibule of the alpha-hemolysin pore. The model presented may be useful in predicting and interpreting DNA detection using nanopore detectors. In addition, this well-defined molecular system may prove useful for investigating models of ligand-gated channels in biological membranes.

General method of preparation of uniformly 13C, 15N-labeled DNA fragments for NMR analysis of DNA structures

(13)C, (15)N labeling of biomolecules allows easier assignments of NMR resonances and provides a larger number of NMR parameters, which greatly improves the quality of DNA structures. However, there is no general DNA-labeling procedure, like those employed for proteins and RNAs. Here, we describe a general and widely applicable approach designed for preparation of isotopically labeled DNA fragments that can be used for NMR studies. The procedure is based on the PCR amplification of oligonucleotides in the presence of labeled deoxynucleotides triphosphates. It allows great flexibility thanks to insertion of a short DNA sequence (linker) between two repeats of DNA sequence to study. Size and sequence of the linker are designed as to create restriction sites at the junctions with DNA of interest. DNA duplex with desired sequence and size is released upon enzymatic digestion of the PCR product. The suitability of the procedure is validated through the preparation of two biological relevant DNA fragments.

Crystal structures of oligonucleotides including the integrase processing site of the Moloney murine leukemia virus

In the first step of retroviral integration, integrase cleaves the linear viral DNA within its long terminal repeat (LTR) immediately 3′ to the CA dinucleotide step, resulting in a reactive 3′ OH on one strand and a 5′ two base overhang on the complementary strand. In order to investigate the structural properties of the 3′ end processing site within the Moloney murine leukemia virus (MMLV) LTR d(TCTTTCATT), a host-guest crystallographic method was employed to determine the structures of four self-complementary 16 bp oligonucleotides including LTR sequences (underlined), d(TTTCATTGCAATGAAA), d(CTTTCATTAATGAAAG), d(TCTTTCATATGAAAGA) and d(CACAATGATCATTGTG), the guests, complexed with the N-terminal fragment of MMLV reverse transcriptase, the host. The structures of the LTR-containing oligonucleotides were compared to those of non-LTR oligonucleotides crystallized in the same lattice. Properties unique to the CA dinucleotide step within the LTR sequence, independent of its position from the end of the duplex, include a positive roll angle and negative slide value. This propensity for the CA dinucleotide step within the MMLV LTR sequence to adopt only positive roll angles is likely influenced by the more rigid, invariable 3′ and 5′ flanking TT dinucleotide steps and may be important for specific recognition and/or cleavage by the MMLV integrase.

Effect of DNA modifications on DNA processing by HIV-1 integrase and inhibitor binding: role of DNA backbone flexibility and an open catalytic site

Integration of the viral cDNA into host chromosomes is required for viral replication. Human immunodeficiency virus integrase catalyzes two sequential reactions, 3'-processing (3'-P) and strand transfer (ST). The first integrase inhibitors are undergoing clinical trial, but interactions of inhibitors with integrase and DNA are not well understood in the absence of a co-crystal structure. To increase our understanding of integrase interactions with DNA, we examined integrase catalysis with oligonucleotides containing DNA backbone, base, and groove modifications placed at unique positions surrounding the 3'-processing site. 3'-Processing was blocked with substrates containing constrained sugars and alpha-anomeric residues, suggesting that integrase requires flexibility of the phosphodiester backbone at the 3'-P site. Of several benzo[a]pyrene 7,8-diol 9,10-epoxide (BaP DE) adducts tested, only the adduct in the minor groove at the 3'-P site inhibited 3'-P, suggesting the importance of the minor groove contacts for 3'-P. ST occurred in the presence of bulky BaP DE DNA adducts attached to the end of the viral DNA suggesting opening of the active site for ST. Position-specific effects of these BaP DE DNA adducts were found for inhibition of integrase by diketo acids. Together, these results demonstrate the importance of DNA structure and specific contacts with the viral DNA processing site for inhibition by integrase inhibitors.

Probing of HIV-1 integrase/DNA interactions using novel analogs of viral DNA

The specific activity of the human immunodeficiency virus, type 1 (HIV-1), integrase on the viral long terminal repeat requires the binding of the enzyme to certain sequences located in the U3 and U5 regions at the ends of viral DNA, but the determinants of this specific DNA-protein recognition are not yet completely understood. We synthesized DNA duplexes mimicking the U5 region and containing either 2'-modified nucleosides or 1,3-propanediol insertions and studied their interactions with HIV-1 integrase, using Mn2+ or Mg2+ ions as integrase cofactors. These DNA modifications had no strong effect on integrase binding to the substrate analogs but significantly affected 3'-end processing rate. The effects of nucleoside modifications at positions 5, 6, and especially 3 strongly depended on the cationic cofactor used. These effects were much more pronounced in the presence of Mg2+ than in the presence of Mn2+. Modifications of base pairs 7-9 affected 3'-end processing equally in the presence of both ions. Adenine from the 3rd bp is thought to form at least two hydrogen bonds with integrase that are crucial for specific DNA recognition. The complementary base, thymine, is not important for integrase activity. For other positions, our results suggest that integrase recognizes a fine structure of the sugar-phosphate backbone rather than heterocyclic bases. Integrase interactions with the unprocessed strand at positions 5-8 are more important than interactions with the processed strand for specific substrate recognition. Based on our results, we suggest a model for integrase interaction with the U5 substrate.

Ribonucleic guanidine demonstrates an unexpected marked preference for complementary DNA rather than RNA

Replacement of the phosphodiester linkages of DNA and RNA by guanidinium linkages provides DNG and RNG. We report here the order of stability of mixed duplexes (RNG-U5.DNA-A5>>RNA-U5.RNA-A5>RNG-U5.RNA-A5>RNA-U5.DNA-A5>DNA-T5.DNA-A5). The considerable stability of RNG.DNA compared to RNG.RNA is shown to be due to the rigid backbone of RNG existing only in B-form and therefore lowering its affinity for A-RNA. RNG oligomers are putative antigene agents which are specific for DNA and would have minimal competitive binding to ncRNA.

Integration requires a specific interaction of the donor DNA terminal 5'-cytosine with glutamine 148 of the HIV-1 integrase flexible loop

Integration is essential for retroviral replication and gene therapy using retroviral vectors. Human immunodeficiency virus, type 1 (HIV-1), integrase specifically recognizes the terminal sequences of each long terminal repeat (LTR) and cleaves the 3'-end terminal dinucleotide 5'-GT. The exposed 3'-hydroxyl is then positioned for nucleophilic attack and subsequent strand transfer into another DNA duplex (target or chromosomal DNA). We report that both the terminal cytosine at the protruding 5'-end of the long terminal repeats (5'-C) and the integrase residue Gln-148 are critical for strand transfer. Proximity of the 5'-C and Gln-148 was demonstrated by disulfide cross-linking. Cross-linking is inhibited by the inhibitor 5CITEP 1-(5-chloroindol-3-yl)-3-hydroxy-3-(2H-tetrazol-5-yl)-propenone. We propose that strand transfer requires a conformational change of the integrase-viral (donor) DNA complex with formation of an H-bond between the N-3 of the 5'-C and the amine group of Gln-148. These findings have implications for the molecular mechanisms coupling 3'-processing and strand transfer as well as for the molecular pharmacology of integrase inhibitors.