Molecular Dynamics Approaches Estimate the Binding Energy of HIV-1 Integrase Inhibitors and Correlate with In Vitro Activity

The design of novel integrase (IN) inhibitors has been aided by recent crystal structures revealing the binding mode of these compounds with a full-length prototype foamy virus (PFV) IN and synthetic viral DNA ends. Earlier docking studies relied on incomplete structures and did not include the contribution of the viral DNA to inhibitor binding. Using the structure of PFV IN as the starting point, we generated a model of the corresponding HIV-1 complex and developed a molecular dynamics (MD)-based approach that correlates with thein vitro activities of novel compounds. Four well-characterized compounds (raltegravir, elvitegravir, MK-0536, and dolutegravir) were used as a training set, and the data for theirin vitro activity against the Y143R, N155H, and G140S/Q148H mutants were used in addition to the wild-type (WT) IN data. Three additional compounds were docked into the IN-DNA complex model and subjected to MD simulations. All three gave interaction potentials within 1 standard deviation of values estimated from the training set, and the most active compound was identified. Additional MD analysis of the raltegravir- and dolutegravir-bound complexes gave internal and interaction energy values that closely match the experimental binding energy of a compound related to raltegravir that has similar activity. These approaches can be used to gain a deeper understanding of the interactions of the inhibitors with the HIV-1 intasome and to identify promising scaffolds for novel integrase inhibitors, in particular, compounds that retain activity against a range of drug-resistant mutants, making it possible to streamline synthesis and testing.