Molecular dynamics and MM/GBSA-integrated protocol probing the correlation between biological activities and binding free energies of HIV-1 TAR RNA inhibitors.

The interaction of HIV-1 transactivator protein Tat with its cognate transactivation response (TAR) RNA has emerged as a promising target for developing antiviral compounds and treating HIV infection, since it is a crucial step for efficient transcription and replication. In the present study, molecular dynamics (MD) simulations and MM/GBSA calculations have been performed on a series of neamine derivatives in order to estimate appropriate MD simulation time for acceptable correlation between ΔGbind and experimental pIC50 values. Initially, all inhibitors were docked into the active site of HIV-1 TAR RNA. Later to explore various conformations and examine the docking results, MD simulations were carried out. Finally, binding free energies were calculated using MM/GBSA method and were correlated with experimental pIC50 values at different time scales (0-1 to 0-10 ns). From this study, it is clear that in case of neamine derivatives as simulation time increased the correlation between binding free energy and experimental pIC50 values increased correspondingly. Therefore, the binding energies which can be interpreted at longer simulation times can be used to predict the bioactivity of new neamine derivatives. Moreover, in this work, we have identified some plausible critical nucleotide interactions with neamine derivatives that are responsible for potent inhibitory activity. Furthermore, we also provide some insights into a new class of oxadiazole-based back bone cyclic peptides designed by incorporating the structural features of neamine derivatives. On the whole, this approach can provide a valuable guidance for designing new potent inhibitors and modify the existing compounds targeting HIV-1 TAR RNA.