Concerted hydrolysis mechanism of HIV-1 natural substrate against subtypes B and C-SA PR: Insight through Molecular Dynamics and Hybrid QM/MM studies

It is well known that understanding the catalytic mechanism of HIV-1 PR is the rationale on which its inhibitors were developed; therefore, a better understanding of the mechanism of natural substrate hydrolysis is important. Herein, the reaction mechanism of HIV-1 natural substrates with subtypes B and common mutant in South Africa (subtype C-SA) protease was studied through transition state modelling, using a general acid-general base (GA-GB) one-step concerted process. The activation free energies of enzyme—substrate complexes were compared based on their rate of hydrolysis using a two-layered ONIOM (B3LYP/6-31++G(d,p):AMBER) method. We expanded our computational model to obtain a better understanding of the mechanism of hydrolysis as well as how the enzyme recognises or chooses the cleavage site of the scissile bonds. Using this model, a potential substrate-based inhibitor could be developed with better potency. The calculated activation energies of natural substrates in our previous study correlated well with experimental data. A similar trend was observed for the Gag and Gag-Pol natural substrates in the present work for both enzyme complexes except for the PR-RT substrate. Natural bond orbital (NBO) analysis was also applied to determine the extent of charge transfer within the QM part of both enzymes considered and the PR-RT natural substrate. The result of this study shows that the method can be utilized as a dependable computational technique to rationalize lead compounds against specific targets.