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.

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.