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.

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.