Accurate prediction of relative binding affinities of a series of HIV-1 protease inhibitors using semi-empirical quantum mechanical charge

Designing drugs when there is low data availability: one-shot learning and other approaches to face the issues of a long-term concern

INTRODUCTION

Modern drug discovery generally is accessed by useful information from previous large databases or uncovering novel data. The lack of biological and/or chemical data tends to slow the development of scientific research and innovation. Here, approaches that may help provide solutions to generate or obtain enough relevant data or improve/accelerate existing methods within the last five years were reviewed.

AREAS COVERED

One-shot learning (OSL) approaches, structural modeling, molecular docking, scoring function space (SFS), molecular dynamics (MD), and quantum mechanics (QM) may be used to amplify the amount of available data to drug design and discovery campaigns, presenting methods, their perspectives, and discussions to be employed in the near future.

EXPERT OPINION

Recent works have successfully used these techniques to solve a range of issues in the face of data scarcity, including complex problems such as the challenging scenario of drug design aimed at intrinsically disordered proteins and the evaluation of potential adverse effects in a clinical scenario. These examples show that it is possible to improve and kickstart research from scarce available data to design and discover new potential drugs.

Ritonavir and xk263 Binding-Unbinding with HIV-1 Protease: Pathways, Energy and Comparison

Understanding non-covalent biomolecular recognition, which includes drug-protein bound states and their binding/unbinding processes, is of fundamental importance in chemistry, biology, and medicine. Fully revealing the factors that govern the binding/unbinding processes can further assist in designing drugs with desired binding kinetics. HIV protease (HIVp) plays an integral role in the HIV life cycle, so it is a prime target for drug therapy. HIVp has flexible flaps, and the binding pocket can be accessible by a ligand via various pathways. Comparing ligand association and dissociation pathways can help elucidate the ligand-protein interactions such as key residues directly involved in the interaction or specific protein conformations that determine the binding of a ligand under certain pathway(s). Here, we investigated the ligand unbinding process for a slow binder, ritonavir, and a fast binder, xk263, by using unbiased all-atom accelerated molecular dynamics (aMD) simulation with a re-seeding approach and an explicit solvent model. Using ritonavir-HIVp and xk263-HIVp ligand-protein systems as cases, we sampled multiple unbinding pathways for each ligand and observed that the two ligands preferred the same unbinding route. However, ritonavir required a greater HIVp motion to dissociate as compared with xk263, which can leave the binding pocket with little conformational change of HIVp. We also observed that ritonavir unbinding pathways involved residues which are associated with drug resistance and are distal from catalytic site. Analyzing HIVp conformations sampled during both ligand-protein binding and unbinding processes revealed significantly more overlapping HIVp conformations for ritonavir-HIVp rather than xk263-HIVp. However, many HIVp conformations are unique in xk263-HIVp unbinding processes. The findings are consistent with previous findings that xk263 prefers an induced-fit model for binding and unbinding, whereas ritonavir favors a conformation selection model. This study deepens our understanding of the dynamic process of ligand unbinding and provides insights into ligand-protein recognition mechanisms and drug discovery.