Premium
Efficiency of a second‐generation HIV‐1 protease inhibitor studied by molecular dynamics and absolute binding free energy calculations
Author(s) -
Lepšík M.,
Kříž Z.,
Havlas Z.
Publication year - 2004
Publication title -
proteins: structure, function, and bioinformatics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.699
H-Index - 191
eISSN - 1097-0134
pISSN - 0887-3585
DOI - 10.1002/prot.20192
Subject(s) - indinavir , molecular dynamics , hiv 1 protease , protease , chemistry , saquinavir , hydrogen bond , side chain , solvation , van der waals force , binding energy , computational chemistry , stereochemistry , human immunodeficiency virus (hiv) , enzyme , biochemistry , molecule , physics , organic chemistry , virology , biology , sida , viral disease , antiretroviral therapy , viral load , nuclear physics , polymer
A subnanomolar inhibitor of human immunodeficiency virus type 1 (HIV‐1) protease, designated QF34, potently inhibits the wild‐type and drug‐resistant enzyme. To explain its broad activity, the binding of QF34 to the wild‐type HIV‐1 protease is investigated by molecular dynamics simulations and compared to the binding of two inhibitors that are used clinically, saquinavir (SQV) and indinavir (IDV). Analysis of the flexibility of protease residues and inhibitor segments in the complex reveals that segments of QF34 were more mobile during the dynamics studies than the segments of SQV and IDV. The dynamics of hydrogen bonding show that QF34 forms a larger number of stable hydrogen bonds than the two inhibitors that are used clinically. Absolute binding free energies were calculated with molecular mechanics–generalized Born surface area (MM–GBSA) methodology using three protocols. The most consistent results were obtained using the single‐trajectory approach, due to cancellation of errors and inadequate sampling in the separate‐trajectory protocols. For all three inhibitors, energy components in favor of binding include van der Waals and electrostatic terms, whereas polar solvation and entropy terms oppose binding. Decomposition of binding energies reveals that more protease residues contribute significantly to the binding of QF34 than to the binding of SQV and IDV. Moreover, contributions from protease main chains and side chains are balanced in the case of QF34 (52:48 ratio, respectively), whereas side chain contributions prevail in both SQV and IDV (main‐chain:side‐chain ratios of 41:59 and 45:55, respectively). The presented results help explain the ability of QF34 to inhibit multiple resistant mutants and should be considered in the design of broad‐specificity second‐generation HIV‐1 protease inhibitors. Proteins 2004. © 2004 Wiley‐Liss, Inc.