Interdependence of Inhibitor Recognition in HIV-1 Protease

Abstract

Molecular recognition is a highly interdependent process. Subsite couplings within the active site of proteases are most often revealed through conditional amino acid preferences in substrate recognition. However, the potential effect of these couplings on inhibition and thus inhibitor design is largely unexplored. The present study examines the interdependency of subsites in HIV-1 protease using a focused library of protease inhibitors, to aid in future inhibitor design. Previously a series of darunavir (DRV) analogs was designed to systematically probe the S1′ and S2′ subsites. Co-crystal structures of these analogs with HIV-1 protease provide the ideal opportunity to probe subsite interdependency. All-atom molecular dynamics simulations starting from these structures were performed and systematically analyzed in terms of atomic fluctuations, intermolecular interactions, and water structure. These analyses reveal that the S1′ subsite highly influences other subsites: the extension of the hydrophobic P1′ moiety results in 1) reduced van der Waals contacts in the P2′ subsite, 2) more variability in the hydrogen bond frequencies with catalytic residues and the flap water, and 3) changes in the occupancy of conserved water sites both proximal and distal to the active site. In addition, one of the monomers in this homodimeric enzyme has atomic fluctuations more highly correlated with DRV than the other monomer. These relationships intricately link the HIV-1 protease subsites and are critical to understanding molecular recognition and inhibitor binding. More broadly, the interdependency of subsite recognition within an active site requires consideration in the selection of chemical moieties in drug design; this strategy is in contrast to what is traditionally done with independent optimization of chemical moieties of an inhibitor.

Exploring molecular mechanism of allosteric inhibitor to relieve drug resistance of multiple mutations in HIV‐1 protease by enhanced conformational sampling

“…Except for the experimental studies, different simulation methods were also employed to partially probe drug resistance induced by the mutated PR . It is worth noting that Schiffer's group used crystallography and simulations to study bindings of inhibitors to PR, and their results clarified drug‐resistant mechanism of mutations in PR and binding modes of inhibitors to PR . Although more traditional approaches, namely design and development of inhibitors targeting the active site, are adopted by many research groups to overcome drug resistance of PR, by now it is still a big challenge to solve this issue.…”

Avoiding Drug Resistance by Substrate Envelope Guided Design: Toward Potent and Robust HCV NS3/4A Protease Inhibitors

“…where r is the distance within 6 Å between atom pairs i of the protease and j of the inhibitor, is the well depth, and is the van der Waals radius (62). Data availability.…”

Selection of HIV-1 for Resistance to Fifth Generation Protease Inhibitors Reveals Two Independent Pathways to High-Level Resistance

“…Selections were performed with 11 analogous inhibitors derived from a common 233 scaffold ( Table 1) , Paulsen et al, 2017. We were interested to see if subtle 234 chemical differences between the inhibitors could impact selection for different resistance 235 pathways.…”