Abstract
HIV-1 protease (HIV-1 PR) is an essential enzyme for the replication process of its virus, and therefore considered an important target for the development of drugs against the acquired immunodeficiency syndrome (AIDS). Our previous study shows that the catalytic mechanism of subtype B/C-SA HIV-1 PR follows a one-step concerted acyclic hydrolysis reaction process using a two-layered ONIOM B3LYP/6-31++G(d,p) method. This present work is aimed at exploring the proposed mechanism of the proteolysis catalyzed by HIV-1 PR and to ensure our proposed mechanism is not an artefact of a single theoretical technique. Hence, we present umbrella sampling method that is suitable for calculating potential mean force (PMF) for non-covalent ligand/substrate-enzyme association/dissociation interactions which provide thermodynamic details for molecular recognition. The free activation energy results were computed in terms of PMF analysis within the hybrid QM(DFTB)/MM approach. The theoretical findings suggest that the proposed mechanism corresponds in principle with experimental data. Given our observations, we suggest that the QM/MM MD method can be used as a reliable computational technique to rationalize lead compounds against specific targets such as the HIV-1 protease.
![](/cms/asset/aae6ee9b-5e77-4c33-9eba-379b1f90e8bc/tbsd_a_1832578_uf0001_c.jpg)
Acknowledgements
The authors thank the College of Health Sciences, University of KwaZulu-Natal, Asphen Pharmacare, Medical Research Council, and the National Research Foundation (all in South Africa) for financial support. The authors also acknowledge the University of Florida Research Computing for providing computational resources and support that contributed to the research results reported in this publication. GK thank AR for hosting ZS and ML at UF during the project. We are also grateful to the Centre for High-Performance Computing (www.chpc.ac.za) for computational resources.
Disclosure statement
The authors declare that they have no competing interests.