RESUMO
Although Venezuelan equine encephalitis virus (VEEV) is a life-threatening pathogen with a capacity for epidemic outbreaks, there are no FDA-approved VEEV antivirals for humans. VEEV cytotoxicity is partially attributed to the formation of a tetrameric complex between the VEEV capsid protein, the nuclear import proteins importin-α and importin-ß, and the nuclear export protein CRM1, which together block trafficking through the nuclear pore complex. Experimental studies have identified small molecules from the CL6662 scaffold as potential inhibitors of the viral nuclear localization signal (NLS) sequence binding to importin-α. However, little is known about the molecular mechanism of CL6662 inhibition. To address this issue, we employed all-atom replica exchange molecular dynamics simulations to probe, in atomistic detail, the binding mechanism of CL6662 ligands to importin-α. Three ligands, including G281-1485 and two congeners with varying hydrophobicities, were considered. We investigated the distribution of ligand binding poses, their locations, and ligand specificities measured by the strength of binding interactions. We found that G281-1485 binds nonspecifically without forming well-defined binding poses throughout the NLS binding site. Binding of the less hydrophobic congener becomes strongly on-target with respect to the NLS binding site but remains nonspecific. However, a more hydrophobic congener is a strongly specific binder and the only ligand out of three to form a well-defined binding pose, while partially overlapping with the NLS binding site. On the basis of free energy estimates, we argue that all three ligands weakly compete with the viral NLS sequence for binding to importin-α in an apparent compromise to preserve host NLS binding. We further show that all-atom replica exchange binding simulations are a viable tool for studying ligands binding nonspecifically without forming well-defined binding poses.