Leveraging Conformational Diversity for Enhanced Structure-Based Virtual Screening: Insights from Molecular Dynamics Simulations of HIV-1 Protease-Ligand Complexes

TL;DR

This study shows that including multiple conformations of proteins and ligands, generated via molecular dynamics, improves the accuracy of structure-based virtual screening for HIV-1 protease inhibitors.

Contribution

It introduces a protocol combining molecular dynamics, clustering, and energy calculations to account for conformational diversity in virtual screening.

Findings

Native pairs have favorable binding energies

Clustering affects computational efficiency and accuracy

Conformational diversity enhances screening reliability

Abstract

Structure-based virtual screening aims to identify high-affinity ligands by estimating binding free energies between proteins and small molecules. However, the conformational flexibility of both proteins and ligands challenges conventional rigid docking methods that assume a fixed receptor structure. In this study, we examined the impact of conformational diversity on binding energy calculations using 79 HIV-1 protease-ligand complexes. Molecular dynamics simulations were employed to generate structural ensembles for both proteins and ligands in aqueous environments. RMSD-based clustering was applied to reduce redundancy while preserving structural diversity. Binding energies were computed using van der Waals and electrostatic interactions. The results demonstrated that native protein-ligand pairs consistently yielded favorable binding energies, whereas non-native pairings often failed to reproduce binding. Furthermore, clustering thresholds influenced the balance between computational cost and interaction accuracy. These findings underscore the importance of incorporating multiple protein and ligand conformations in SBVS protocols to improve prediction reliability and support more effective drug discovery strategies.