Search

The three pan-genotypic HCV NS3/4A protease inhibitors (PIs) currently in clinical use-grazoprevir, glecaprevir, and voxilaprevir-are quinoxaline-based P2-P4 macrocycles and thus exhibit similar resistance profiles. Using our quinoxaline-based P1-P3 macrocyclic lead compounds as an alternative chemical scaffold, we explored structure-activity relationships (SARs) at the P2 and P4 positions to develop pan-genotypic PIs that avoid drug resistance. A structure-guided strategy was used to design and synthesize two series of compounds with different P2 quinoxalines in combination with diverse P4 groups of varying sizes and shapes, with and without fluorine substitutions. Our SAR data and cocrystal structures revealed the interplay between the P2 and P4 groups, which influenced inhibitor binding and the overall resistance profile. Optimizing inhibitor interactions in the S4 pocket led to PIs with excellent antiviral activity against clinically relevant PI-resistant HCV variants and genotype 3, providing potential pan-genotypic inhibitors with improved resistance profiles.

Third generation Hepatitis C virus (HCV) NS3/4A protease inhibitors (PIs), glecaprevir and voxilaprevir, are highly effective across genotypes and against many resistant variants. Unlike earlier PIs, these compounds have fluorine substitutions on the P2-P4 macrocycle and P1 moieties. Fluorination has long been used in medicinal chemistry as a strategy to improve physicochemical properties and potency. However, the molecular basis by which fluorination improves potency and resistance profile of HCV NS3/4A PIs is not well understood. To systematically analyze the contribution of fluorine substitutions to inhibitor potency and resistance profile, we used a multi-disciplinary approach involving inhibitor design and synthesis, enzyme inhibition assays, co-crystallography, and structural analysis. A panel of inhibitors in matched pairs were designed with and without P4 cap fluorination, tested against WT protease and the D168A resistant variant, and a total of 22 high-resolution co-crystal structures were determined. While fluorination did not significantly improve potency against the WT protease, PIs with fluorinated P4 caps retained much better potency against the D168A protease variant. Detailed analysis of the co-crystal structures revealed that PIs with fluorinated P4 caps can sample alternate binding conformations that enable adapting to structural changes induced by the D168A substitution. Our results elucidate molecular mechanisms of fluorine-specific inhibitor interactions that can be leveraged in avoiding drug resistance.

Rupintrivir targets the 3C cysteine proteases of the picornaviridae family, which includes rhinoviruses and enteroviruses that cause a range of human diseases. Despite being a pan-3C protease inhibitor, rupintrivir activity is extremely weak against the homologous 3C-like protease of SARS-CoV-2. In this study, the crystal structures of rupintrivir were determined bound to enterovirus 68 (EV68) 3C protease and the 3C-like main protease (Mpro) from SARS-CoV-2. While the EV68 3C protease-rupintrivir structure was similar to previously determined complexes with other picornavirus 3C proteases, rupintrivir bound in a unique conformation to the active site of SARS-CoV-2 Mpro splitting the catalytic cysteine and histidine residues. This bifurcation of the catalytic dyad may provide a novel approach for inhibiting cysteine proteases.

A structure-guided design strategy was used to improve the resistance profile of HIV-1 protease inhibitors by optimizing hydrogen bonding and van der Waals interactions with the protease while staying within the substrate envelope. Stereoisomers of 4-(1-hydroxyethyl)benzene and 4-(1,2-dihydroxyethyl)benzene moieties were explored as P2' ligands providing pairs of diastereoisomers epimeric at P2', which exhibited distinct potency profiles depending on the configuration of the hydroxyl group and size of the P1' group. While compounds with the 4-(1-hydroxyethyl)benzene P2' moiety maintained excellent antiviral potency against a panel of multidrug-resistant HIV-1 strains, analogues with the polar 4-(1,2-dihydroxyethyl)benzene moiety were less potent, and only the (R)-epimer incorporating a larger 2-ethylbutyl P1' group showed improved potency. Crystal structures of protease-inhibitor complexes revealed strong hydrogen bonding interactions of both (R)- and (S)-stereoisomers of the hydroxyethyl group with Asp30'. Notably, the (R)-dihydroxyethyl group was involved in a unique pattern of direct hydrogen bonding interactions with the backbone amides of Asp29' and Asp30'. The SAR data and analysis of crystal structures provide insights for optimizing these promising HIV-1 protease inhibitors.

Drug resistance threatens many critical therapeutics through mutations in the drug target. The molecular mechanisms by which combinations of mutations, especially involving those distal from the active site, alter drug binding to confer resistance are poorly understood and thus difficult to counteract. A machine learning strategy was developed that couples parallel molecular dynamics simulations and experimental potency to identify specific conserved mechanisms underlying resistance. A series of 28 HIV-1 protease variants with 0-24 substitutions each were used as a rigorous model of this strategy. Many of the mutations were distal from the active site and the potency of variants to a drug (darunavir) varied from low picomolar to near micromolar. With features extracted from the simulations, elastic network machine learning was applied to correlate physical interactions with loss of potency and succeeded to within 1 kcal/mol of experimental affinity for both the training and test sets, outperforming MM/GBSA calculations. Feature reduction resulted in a model with 4 specific features that describe interactions critical for potency for all 28 variants. These predictive features, that specifically vary with potency, occur throughout the enzyme and would not have been identified without dynamics and machine learning. This strategy thus captures the conserved dynamic mechanisms by which complex combinations of mutations confer resistance and identifies critical features that serve as bellwethers of loss of inhibitor potency. Machine learning models leveraging molecular dynamics can thus elucidate mechanisms of drug resistance that confer loss of affinity and will serve as predictive tools in future drug design.

Human T-cell lymphotropic virus type 1 (HTLV-1) is a retrovirus that can cause severe paralytic neurologic disease and immune disorders as well as cancer. An estimated 20 million people worldwide are infected with HTLV-1, with prevalence reaching 30% in some parts of the world. In stark contrast to HIV-1, no direct acting antivirals (DAAs) exist against HTLV-1. The aspartyl protease of HTLV-1 is a dimer similar to that of HIV-1 and processes the viral polyprotein to permit viral maturation. We report that the FDA-approved HIV-1 protease inhibitor darunavir (DRV) inhibits the enzyme with 0.8 µM potency and provides a scaffold for drug design against HTLV-1. Analogs of DRV that we designed and synthesized achieved sub-micromolar inhibition against HTLV-1 protease and inhibited Gag processing in viral maturation assays and in a chronically HTLV-1 infected cell line. Cocrystal structures of these inhibitors with HTLV-1 protease highlight opportunities for future inhibitor design. Our results show promise toward developing highly potent HTLV-1 protease inhibitors as therapeutic agents against HTLV-1 infections.

Viral proteases are critical enzymes for the maturation of many human pathogenic viruses and thus are key targets for direct acting antivirals (DAAs). The current viral pandemic caused by SARS-CoV-2 is in dire need of DAAs. The Main protease (Mpro) is the focus of extensive structure-based drug design efforts which are mostly covalent inhibitors targeting the catalytic cysteine. ML188 is a non-covalent inhibitor designed to target SARS-CoV-1 Mpro, and provides an initial scaffold for the creation of effective pan-coronavirus inhibitors. In the current study, we found that ML188 inhibits SARS-CoV-2 Mpro at 2.5 &micro, M, which is more potent than against SAR-CoV-1 Mpro. We determined the crystal structure of ML188 in complex with SARS-CoV-2 Mpro to 2.39 &Aring, resolution. Sharing 96% sequence identity, structural comparison of the two complexes only shows subtle differences. Non-covalent protease inhibitors complement the design of covalent inhibitors against SARS-CoV-2 main protease and are critical initial steps in the design of DAAs to treat CoVID 19.

Viral proteases are diverse in structure, oligomeric state, catalytic mechanism, and substrate specificity. This chapter focuses on proteases from viruses that are relevant to human health: human immunodeficiency virus subtype 1 (HIV-1), hepatitis C (HCV), human T-cell leukemia virus type 1 (HTLV-1), flaviviruses, enteroviruses, and coronaviruses. The proteases of HIV-1 and HCV have been successfully targeted for therapeutics, with picomolar FDA-approved drugs currently used in the clinic. The proteases of HTLV-1 and the other virus families remain emerging therapeutic targets at different stages of the drug development process. This chapter provides an overview of the current knowledge on viral protease structure, mechanism, substrate recognition, and inhibition. Particular focus is placed on recent advances in understanding the molecular basis of diverse substrate recognition and resistance, which is essential toward designing novel protease inhibitors as antivirals.

The design, synthesis, and X-ray structural analysis of hybrid HIV-1 protease inhibitors (PIs) containing bis-tetrahydrofuran (bis-THF) in a pseudo-C2-symmetric dipeptide isostere are described. A series of PIs were synthesized by incorporating bis-THF of darunavir on either side of the Phe-Phe isostere of lopinavir in combination with hydrophobic amino acids on the opposite P2/P2' position. Structure-activity relationship studies indicated that the bis-THF moiety can be attached at either the P2 or P2' position without significantly affecting potency. However, the group on the opposite P2/P2' position had a dramatic effect on potency depending on the size and shape of the side chain. Cocrystal structures of inhibitors with wild-type HIV-1 protease revealed that the bis-THF moiety retained similar interactions as observed in the darunavir-protease complex regardless of the position on the Phe-Phe isostere. Analyses of cocrystal structures and molecular dynamics simulations provide insights into optimizing HIV-1 PIs containing bis-THF in non-sulfonamide dipeptide isosteres.

Despite significant progress, hepatitis C virus (HCV) continues to be a major health problem with millions of people infected worldwide and thousands dying annually due to resulting complications. Recent antiviral combinations can achieve >95% cure, but late diagnosis, low access to treatment, and treatment failure due to drug resistance continue to be roadblocks against eradication of the virus. We report the rational design of two series of HCV NS3/4A protease inhibitors with improved resistance profiles by exploiting evolutionarily constrained regions of the active site using the substrate envelope model. Optimally filling the S4 pocket is critical to avoid resistance and improve potency. Our results provide drug design strategies to avoid resistance that are applicable to other quickly evolving viral drug targets., Hepatitis C virus (HCV) infects millions of people worldwide, causing chronic liver disease that can lead to cirrhosis, hepatocellular carcinoma, and liver transplant. In the last several years, the advent of direct-acting antivirals, including NS3/4A protease inhibitors (PIs), has remarkably improved treatment outcomes of HCV-infected patients. However, selection of resistance-associated substitutions and polymorphisms among genotypes can lead to drug resistance and in some cases treatment failure. A proactive strategy to combat resistance is to constrain PIs within evolutionarily conserved regions in the protease active site. Designing PIs using the substrate envelope is a rational strategy to decrease the susceptibility to resistance by using the constraints of substrate recognition. We successfully designed two series of HCV NS3/4A PIs to leverage unexploited areas in the substrate envelope to improve potency, specifically against resistance-associated substitutions at D168. Our design strategy achieved better resistance profiles over both the FDA-approved NS3/4A PI grazoprevir and the parent compound against the clinically relevant D168A substitution. Crystallographic structural analysis and inhibition assays confirmed that optimally filling the substrate envelope is critical to improve inhibitor potency while avoiding resistance. Specifically, inhibitors that enhanced hydrophobic packing in the S4 pocket and avoided an energetically frustrated pocket performed the best. Thus, the HCV substrate envelope proved to be a powerful tool to design robust PIs, offering a strategy that can be translated to other targets for rational design of inhibitors with improved potency and resistance profiles.

Hepatitis C virus (HCV), causative agent of chronic viral hepatitis, infects 71 million people worldwide and is divided into seven genotypes and multiple subtypes with sequence identities between 68 to 82%. While older generation direct-acting antivirals (DAAs) had varying effectiveness against different genotypes, the newest NS3/4A protease inhibitors including glecaprevir (GLE) have pan-genotypic activity. The structural basis for pan-genotypic inhibition and effects of polymorphisms on inhibitor potency were not well known due to lack of crystal structures of GLE-bound NS3/4A or genotypes other than 1. In this study, we determined the crystal structures of NS3/4A from genotypes 1a, 3a, 4a and 5a in complex with GLE. Comparison with the highly similar grazoprevir (GZR) indicated the mechanism of GLE’s drastic improvement in potency. We found that while GLE is highly potent against wild type NS3/4A of all genotypes, specific resistance-associated substitutions (RASs) confer orders of magnitude loss in inhibition. Our crystal structures reveal molecular mechanisms behind pan-genotypic activity of GLE, including potency loss due to RASs at D168. Our structures permit for the first time analysis of changes due to polymorphisms among genotypes, providing insights into design principles that can aid future drug development and potentially can be extended to other proteins.

Summary Well-designed viral protease inhibitors (PIs) potently inhibit replication as well as create a high genetic barrier for resistance. Through in vivo selective pressure, we have generated high-level resistance against ten HIV-1 PIs and their precursor, the FDA-approved drug darunavir (DRV), achieving 1,000-fold resistance over the starting EC50. The accumulation of mutations revealed two pathways to high-level resistance, resulting in protease variants with up to 14 mutations in and outside of the active site. The two pathways demonstrate the interplay between drug resistance and viral fitness. Replicate selections showed that one inhibitor could select for resistance through either pathway, although subtle changes in chemical structure of the inhibitors led to preferential use of one pathway over the other. Viral variants from the two pathways showed differential selection of compensatory mutations in Gag cleavage sites. These results reveal the high-level of selective pressure that is attainable with these fourth-generation protease inhibitors, and the interplay between selection of mutations to confer resistance while maintaining viral fitness.

Water is an essential in many biological processes, and the hydration structure plays a critical role in facilitating protein folding, dynamics and ligand binding. A variety of biophysical spectroscopic techniques have been used to probe the water solvating proteins, often complemented with molecular dynamics (MD) simulations to resolve the spatial and dynamic features of the hydration shell, but comparing relative water structure is challenging. In this study 1 microsecond MD simulations were performed to identify and characterize hydration sites around HIV-1 protease bound to an inhibitor, darunavir (DRV). The water density, hydration site occupancy, extent and anisotropy of fluctuations, coordinated water molecules, and hydrogen bonds were characterized and compared to the properties of bulk water. The water density of the principal hydration shell was found to be higher than bulk, dependent on the topology and physio-chemical identity of the biomolecular surface. The dynamics of water molecules occupying principal hydration sites was highly dependent on the number of water-water interactions, and inversely correlated with hydrogen bonds to the protein–inhibitor complex. While many waters were conserved following the symmetry of homodimeric HIV protease, the asymmetry induced by DRV resulted in asymmetric lower-occupancy hydration sites at the concave surface of the active site. Key interactions between water molecules and the protease, that stabilize the protein in the inhibited form, were altered in a drug resistant variant of the protease indicating that modulation of solvent–solute interactions might play a key role in conveying drug resistance. Our analysis provides insights into the interplay between an enzyme inhibitor complex and the hydration shell and has implications in elucidating water structure in a variety of biological processes and applications including ligand binding, inhibitor design and resistance.

Protease inhibitors have the highest potency among antiviral therapies against HIV-1 infections, yet the virus can evolve resistance. Darunavir (DRV), currently the most potent Food and Drug Administration-approved protease inhibitor, retains potency against single-site mutations. However, complex combinations of mutations can confer resistance to DRV. While the interdependence between mutations within HIV-1 protease is key for inhibitor potency, the molecular mechanisms that underlie this control remain largely unknown. In this study, we investigated the interdependence between the L89V and L90M mutations and their effects on DRV binding. These two mutations have been reported to be positively correlated with one another in HIV-1 patient-derived protease isolates, with the presence of one mutation making the probability of the occurrence of the second mutation more likely. The focus of our investigation is a patient-derived isolate, with 24 mutations that we call "KY"; this variant includes the L89V and L90M mutations. Three additional KY variants with back-mutations, KY(V89L), KY(M90L), and the KY(V89L/M90L) double mutation, were used to experimentally assess the individual and combined effects of these mutations on DRV inhibition and substrate processing. The enzymatic assays revealed that the KY(V89L) variant, with methionine at residue 90, is highly resistant, but its catalytic function is compromised. When a leucine to valine mutation at residue 89 is present simultaneously with the L90M mutation, a rescue of catalytic efficiency is observed. Molecular dynamics simulations of these DRV-bound protease variants reveal how the L90M mutation induces structural changes throughout the enzyme that undermine the binding interactions.

Over the last two decades, both the sensitivity of NMR and the time scale of molecular dynamics (MD) simulation have increased tremendously and have advanced the field of protein dynamics. HIV-1 protease has been extensively studied using these two methods, and has presented a framework for cross-evaluation of structural ensembles and internal dynamics by integrating the two methods. Here, we review studies from our laboratories over the last several years, to understand the mechanistic basis of protease drug-resistance mutations and inhibitor responses, using NMR and crystal structure-based parallel MD simulations. Our studies demonstrate that NMR relaxation experiments, together with crystal structures and MD simulations, significantly contributed to the current understanding of structural/dynamic changes due to HIV-1 protease drug resistance mutations.

Dengue virus (DENV), transmitted predominantly in tropical and subtropical regions by the mosquito Aedes aegypti, infects millions of people and leads to dengue fever and thousands of deaths each year. There are no direct-acting antivirals to combat DENV, and molecular and structural knowledge is required to develop such compounds. The dengue NS2B/NS3 protease is a promising target for direct-acting antivirals, as viral polyprotein cleavage during replication is required for the maturation of the viral particle. The NS2B/NS3 protease processes 8 of the 13 viral polyprotein cleavage sites to allow viral maturation. Although these sites share little sequence homology beyond the P1 and P2 positions, most are well conserved among the serotypes. How the other substrate residues, especially at the P′ side, affect substrate recognition remains unclear. We exploited the tight-binding general serine protease inhibitor aprotinin to investigate protease–substrate interactions at the molecular level. We engineered aprotinin’s binding loop with sequences mimicking the P′ side of DENV substrates. P′ residues significantly modulate substrate affinity to protease, with inhibition constants varying from nanomolar to sub-millimolar. Structural and dynamic analysis revealed the molecular basis of this modulation and allowed identifying optimal residues for each of the P′ positions. In addition, isothermal titration calorimetry showed binding to be solely entropy driven for all constructs. Potential flaviviral P′ side inhibitors could benefit from mimicking the optimal residues at P′ positions and incorporate hydrophobicity and rigidity to maintain entropic advantage for potency.

Drug resistance is a major problem in health care, undermining therapy outcomes and necessitating novel approaches to drug design. Extensive studies on resistance to viral protease inhibitors, particularly those of HIV-1 and hepatitis C virus (HCV) protease, revealed a plethora of information on the structural and molecular mechanisms underlying resistance. These insights led to several strategies to improve viral protease inhibitors to counter resistance, such as exploiting the essential biological function and leveraging evolutionary constraints. Incorporation of these strategies into structure-based drug design can minimize vulnerability to resistance, not only for viral proteases but for other quickly evolving drug targets as well, toward designing inhibitors one step ahead of evolution to counter resistance with more intelligent and rational design.

HIV-1 protease is one of the prime targets of agents used in antiretroviral therapy against HIV. However, under selective pressure of protease inhibitors, primary mutations at the active site weaken inhibitor binding to confer resistance. Darunavir (DRV) is the most potent HIV-1 protease inhibitor in clinic; resistance is limited, as DRV fits well within the substrate envelope. Nevertheless, resistance is observed due to hydrophobic changes at residues including I50, V82 and I84 that line the S1/S1’ pocket within the active site. Through enzyme inhibition assays and a series of 12 crystal structures, we interrogated susceptibility of DRV and two potent analogs to primary S1’ mutations. The analogs had modifications at the hydrophobic P1’ moiety to better occupy the unexploited space in the S1’ pocket where the primary mutations were located. Considerable losses of potency were observed against protease variants with I84V and I50V mutations for all three inhibitors. The crystal structures revealed an unexpected conformational change in the flap region of I50V protease bound to the analog with the largest P1’ moiety, indicating interdependency between the S1’ subsite and the flap region. Collective analysis of protease-inhibitor van der Waals (vdW) interactions in the crystal structures using principle component analysis indicated I84V mutation underlying the largest variation in the vdW contacts. Interestingly, the principle components were able to distinguish inhibitor identity and relative potency solely based on vdW interactions of active site residues in the crystal structures. Our results reveal the interplay between inhibitor P1’ moiety and primary S1’ mutations, as well as suggesting a novel method for distinguishing the interdependence of resistance through principle component analyses.

[Image: see text] A series of linear HCV NS3/4A protease inhibitors was designed by eliminating the P2–P4 macrocyclic linker in grazoprevir, which, in addition to conferring conformational flexibility, allowed structure–activity relationship (SAR) exploration of diverse quinoxalines at the P2 position. Biochemical and replicon data indicated preference for small hydrophobic groups at the 3-position of P2 quinoxaline for maintaining potency against resistant variants R155K, A156T, and D168A/V. The linear inhibitors, though generally less potent than the corresponding macrocyclic analogues, were relatively easier to synthesize and less susceptible to drug resistance. Three inhibitor cocrystal structures bound to wild-type NS3/4A protease revealed a conformation with subtle changes in the binding of P2 quinoxaline, depending on the 3-position substituent, likely impacting both inhibitor potency and resistance profile. The SAR and structural analysis highlight inhibitor features that strengthen interactions of the P2 moiety with the catalytic triad residues, providing valuable insights to improve potency against resistant variants.

The mosquito-transmitted dengue virus (DENV) infects millions of people in tropical and subtropical regions. Maturation of DENV particles requires proper cleavage of the viral polyprotein, including processing of 8 of the 13 substrate cleavage sites by dengue virus NS2B/NS3 protease. With no available direct-acting antiviral targeting DENV, NS2/NS3 protease is a promising target for inhibitor design. Current design efforts focus on the nonprime side of the DENV protease active site, resulting in highly hydrophilic and nonspecific scaffolds. However, the prime side also significantly modulates DENV protease binding affinity, as revealed by engineering the binding loop of aprotinin, a small protein with high affinity for DENV protease. In this study, we designed a series of cyclic peptides interacting with both sides of the active site as inhibitors of dengue virus protease. The design was based on two aprotinin loops and aimed to leverage both key specific interactions of substrate sequences and the entropic advantage driving aprotinin's high affinity. By optimizing the cyclization linker, length, and amino acid sequence, the tightest cyclic peptide achieved a K i value of 2.9 μM against DENV3 wild-type (WT) protease. These inhibitors provide proof of concept that both sides of DENV protease active site can be exploited to potentially achieve specificity and lower hydrophilicity in the design of inhibitors targeting DENV. IMPORTANCE Viruses of the flaviviral family, including DENV and Zika virus transmitted by Aedes aegypti , continue to be a threat to global health by causing major outbreaks in tropical and subtropical regions, with no available direct-acting antivirals for treatment. A better understanding of the molecular requirements for the design of potent and specific inhibitors against flaviviral proteins will contribute to the development of targeted therapies for infections by these viruses. The cyclic peptides reported here as DENV protease inhibitors provide novel scaffolds that enable exploiting the prime side of the protease active site, with the aim of achieving better specificity and lower hydrophilicity than those of current scaffolds in the design of antiflaviviral inhibitors.

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.

HIV-1 protease inhibitors (PIs) are competitive active-site inhibitors that mimic the transition state of the enzyme’s substrate, and are the most potent antiretroviral drugs against HIV infection. HIV-1 protease processes the viral polyproteins at specific cleavage sites and allows infectious mature virions and hence spread of the virus. Unfortunately rapid viral evolution combined with selective pressure of therapy causes selection of many drug-resistant variants that are no longer efficiently inhibited by the PIs. HIV-1 protease can tolerate extensive mutations, with close to half of the 99-residues making up each of the chains in the homodimeric protease and residues at substrate cleavage sites mutating to escape PI pressure. Structural and biophysical studies of many drug-resistant HIV-1 protease variants revealed insights into how mutations at and outside of the protease active site are able to confer PI resistance while still allowing recognition and processing of substrates, and why substrate mutations coevolve with primary protease mutations. We summarize the main molecular mechanisms underlying PI resistance due to primary, secondary, and substrate coevolved mutations and how this knowledge may guide the design of robust inhibitors to avoid resistance.

Under the selective pressure of therapy, HIV-1 protease mutants resistant to inhibitors evolve to confer drug resistance. Such mutations can impact both the dynamics and structures of the bound and unbound forms of the enzyme. Flap+ is a multidrug-resistant variant of HIV-1 protease with a combination of primary and secondary resistance mutations (L10I, G48V, I54V, V82A) and a strikingly altered thermodynamic profile for darunavir (DRV) binding relative to the wild-type protease. We elucidated the impact of these mutations on protein dynamics in the DRV-bound state using molecular dynamics simulations and NMR relaxation experiments. Both methods concur in that the conformational ensemble and dynamics of protease are impacted by the drug resistance mutations in Flap+ variant. Surprisingly this change in ensemble dynamics is different from that observed in the unliganded form of the same variant (Cai, Y. et al. J. Chem. Theory Comput.2012, 8, 3452–3462). Our comparative analysis of both inhibitor-free and bound states presents a comprehensive picture of the altered dynamics in drug-resistant mutant HIV-1 protease and underlies the importance of incorporating dynamic analysis of the whole system, including the unliganded state, into revealing drug resistance mechanisms.

The rapid evolution of HIV under selective drug pressure has led to multidrug resistant (MDR) strains that evade standard therapies. We designed highly potent HIV-1 protease inhibitors (PIs) using the substrate envelope model, which confines inhibitors within the consensus volume of natural substrates, providing inhibitors less susceptible to resistance because a mutation affecting such inhibitors will simultaneously affect viral substrate processing. The designed PIs share a common chemical scaffold but utilize various moieties that optimally fill the substrate envelope, as confirmed by crystal structures. The designed PIs retain robust binding to MDR protease variants and display exceptional antiviral potencies against different clades of HIV as well as a panel of 12 drug-resistant viral strains. The substrate envelope model proves to be a powerful strategy to develop potent and robust inhibitors that avoid drug resistance., National Institute of General Medical Sciences (U.S.) (Grant P01-GM66524), National Institute of General Medical Sciences (U.S.) (Grant AI41404), National Institute of General Medical Sciences (U.S.) (Grant AI43198), National Institute of General Medical Sciences (U.S.) (Grant GM065418), National Institute of General Medical Sciences (U.S.) (Grant GM082209), United States. American Recovery and Reinvestment Act of 2009 (Supplement P01GM066524-08S1)

Drug resistance occurs through a series of subtle changes that maintain substrate recognition but no longer permit inhibitor binding. In HIV-1 protease, mutations at I50 are associated with such subtle changes that confer differential resistance to specific inhibitors. Residue I50 is located at the protease flap tips, closing the active site upon ligand binding. Under selective drug pressure, I50V/L substitutions emerge in patients, compromising drug susceptibility and leading to treatment failure. The I50V substitution is often associated with amprenavir (APV) and darunavir (DRV) resistance, while the I50L substitution is observed in patients failing atazanavir (ATV) therapy. To explain how APV, DRV, and ATV susceptibility are influenced by mutations at residue 50 in HIV-1 protease, structural and binding thermodynamics studies were carried out on I50V/L-substituted protease variants in the compensatory mutation A71V background. Reduced affinity to both I50V/A71V and I50L/A71V double mutants is largely due to decreased binding entropy, which is compensated for by enhanced enthalpy for ATV binding to I50V variants and APV binding to I50L variants, leading to hypersusceptibility in these two cases. Analysis of the crystal structures showed that the substitutions at residue 50 affect how APV, DRV, and ATV bind the protease with altered van der Waals interactions and that the selection of I50V versus I50L is greatly influenced by the chemical moieties at the P1 position for APV/DRV and the P2 position for ATV. Thus, the varied inhibitor susceptibilities of I50V/L protease variants are largely a direct consequence of the interdependent changes in protease inhibitor interactions.

Understanding the interdependence of multiple mutations in conferring drug resistance is crucial to the development of novel and robust inhibitors. As HIV-1 protease continues to adapt and evade inhibitors while still maintaining the ability to specifically recognize and efficiently cleave its substrates, the problem of drug resistance has become more complicated. Under the selective pressure of therapy, correlated mutations accumulate throughout the enzyme to compromise inhibitor binding, but characterizing their energetic interdependency is not straightforward. A particular drug resistant variant (L10I/G48V/I54V/V82A) displays extreme entropy-enthalpy compensation relative to wild-type enzyme but a similar variant (L10I/G48V/I54A/V82A) does not. Individual mutations of sites in the flaps (residues 48 and 54) of the enzyme reveal that the thermodynamic effects are not additive. Rather, the thermodynamic profile of the variants is interdependent on the cooperative effects exerted by a particular combination of mutations simultaneously present.

Hepatitis C virus (HCV), affecting an estimated 150 million people worldwide, is the leading cause of viral hepatitis, cirrhosis and hepatocellular carcinoma. HCV is genetically diverse with six genotypes (GTs) and multiple subtypes of different global distribution and prevalence. Recent development of direct-acting antivirals against HCV including NS3/4A protease inhibitors (PIs) has greatly improved treatment outcomes for GT-1. However, all current PIs exhibit significantly lower potency against GT-3. Lack of structural data on GT-3 protease has limited our ability to understand PI failure in GT-3. In this study the molecular basis for reduced potency of current inhibitors against GT-3 NS3/4A protease is elucidated with structure determination, molecular dynamics simulations and inhibition assays. A chimeric GT-1a3a NS3/4A protease amenable to crystallization was engineered to recapitulate decreased sensitivity of GT-3 protease to PIs. High-resolution crystal structures of this GT-1a3a bound to 3 PIs, asunaprevir, danoprevir and vaniprevir, had only subtle differences relative to GT-1 despite orders of magnitude loss in affinity. In contrast, hydrogen-bonding interactions within and with the protease active site and dynamic fluctuations of the PIs were drastically altered. The correlation between loss of intermolecular dynamics and inhibitor potency suggests a mechanism where polymorphisms between genotypes (or selected mutations) in the drug target confer resistance through altering the intermolecular dynamics of the protein-inhibitor complex.

Recent advances in direct-acting antivirals against Hepatitis C Virus (HCV) have led to the development of potent inhibitors, including MK-5172, that target the viral NS3/4A protease with relatively low susceptibility to resistance. MK-5172 has a P2–P4 macrocycle and a unique binding mode among current protease inhibitors where the P2 quinoxaline packs against the catalytic residues H57 and D81. However, the effect of macrocyclization on this binding mode is not clear, as is the relation between macrocyclization, thermodynamic stabilization, and susceptibility to the resistance mutation A156T. We have determined high-resolution crystal structures of linear and P1–P3 macrocyclic analogs of MK-5172 bound to WT and A156T protease and compared these structures, their molecular dynamics, and experimental binding thermodynamics to the parent compound. We find that the “unique” binding mode of MK-5172 is conserved even when the P2–P4 macrocycle is removed or replaced with a P1–P3 macrocycle. While beneficial to decreasing the entropic penalty associated with binding, the constraint exerted by the P2–P4 macrocycle prevents efficient rearrangement to accommodate the A156T mutation, a deficit alleviated in the linear and P1–P3 analogs. Design of macrocyclic inhibitors against NS3/4A needs to achieve the best balance between exerting optimal conformational constraint for enhancing potency, fitting within the substrate envelope and allowing adaptability to be robust against resistance mutations.

HIV-1 protease is a major drug target against AIDS as it permits viral maturation by processing the gag and pol polyproteins of the virus. The cleavage sites in these polyproteins do not have obvious sequence homology or a binding motif and the specificity of the protease is not easily determined. We used various threading approaches, together with the crystal structures of substrate complexes which served as template structures, to study the substrate specificity of HIV-1 protease with the aim of obtaining a better differentiation between binding and nonbinding sequences. The predictions from threading improved when distance-dependent interaction energy functions were used instead of contact matrices. To rank the peptides and properly account for the peptide's conformation in the total energy, the results from using short-range potentials on multiple template structures were averaged. Finally, a dynamic threading approach is introduced which is potentially useful for cases when there is only one template structure available. The conformational energy of the peptide—especially the term accounting for the side chains—was found to be important in differentiating between binding and nonbinding sequences. Hence, the substrate specificity, and thus the ability of the virus to mature, is affected by the compatibility of the substrate peptide to fit within the limited conformational space of the active site groove.

The dynamics of HIV-1 protease, both in unliganded and substrate-bound forms have been analyzed by using an analytical method, Gaussian network model (GNM). The method is applied to different conformations accessible to the protein backbone in the native state, observed in crystal structures and snapshots from fully atomistic molecular dynamics (MD) simulation trajectories. The modes of motion obtained from GNM on different conformations of HIV-1 protease are conserved throughout the MD simulations. The flaps and 40's loop of the unliganded HIV-1 protease structure are identified as the most mobile regions. However, in the liganded structure these flaps lose mobility, and terminal regions of the monomers become more flexible. Analysis of the fast modes shows that residues important for stability are in the same regions of all the structures examined. Among these, Gly86 appears to be a key residue for stability. The contribution of residues in the active site region and flaps to the stability is more pronounced in the substrate-bound form than in the unliganded form. The convergence of modes in GNM to similar regions of HIV-1 protease, regardless of the conformation of the protein, supports the robustness of GNM as a potentially useful and predictive tool.

Drug resistance mutations in response to HIV-1 protease inhibitors are selected not only in the drug target but elsewhere in the viral genome, especially at the protease cleavage sites in the precursor protein Gag. To understand the molecular basis of this protease–substrate coevolution, we solved the crystal structures of drug resistant I50V/A71V HIV-1 protease with p1-p6 substrates bearing coevolved mutations. Analyses of the protease–substrate interactions reveal that compensatory coevolved mutations in the substrate do not restore interactions lost due to protease mutations, but instead establish other interactions that are not restricted to the site of mutation. Mutation of a substrate residue has distal effects on other residues’ interactions as well, including through the induction of a conformational change in the protease. Additionally, molecular dynamics simulations suggest that restoration of active site dynamics is an additional constraint in the selection of coevolved mutations. Hence, protease–substrate coevolution permits mutational, structural, and dynamic changes via molecular mechanisms that involve distal effects contributing to drug resistance.

The development of HIV-1 protease inhibitors has been the historic paradigm of rational structure-based drug design, where structural and thermodynamic analyses have assisted in the discovery of novel inhibitors. While the total enthalpy and entropy change upon binding determine the affinity, often the thermodynamics are considered in terms of inhibitor properties only. In the current study, profound changes are observed in the binding thermodynamics of a drug-resistant variant compared to wild-type HIV-1 protease, irrespective of the inhibitor bound. This variant (Flap+) has a combination of flap and active site mutations and exhibits extremely large entropy-enthalpy compensation compared to wild-type protease, 5-15 kcal/mol, while losing only 1-3 kcal/mol in total binding free energy for any of six FDA-approved inhibitors. Although entropy-enthalpy compensation has been previously observed for a variety of systems, never have changes of this magnitude been reported. The co-crystal structures of Flap+ protease with four of the inhibitors were determined and compared with complexes of both the wild-type protease and another drug-resistant variant that does not exhibit this energetic compensation. Structural changes conserved across the Flap+ complexes, which are more pronounced for the flaps covering the active site, likely contribute to the thermodynamic compensation. The finding that drug-resistant mutations can profoundly modulate the relative thermodynamic properties of a therapeutic target independent of the inhibitor presents a new challenge for rational drug design.

In the rapidly evolving disease of HIV drug resistance readily emerges, nullifying the effectiveness of therapy. Drug resistance has been extensively studied in HIV-1 protease where resistance occurs when the balance between enzyme inhibition and substrate recognition and turn-over is perturbed to favor catalytic activity. Mutations which confer drug resistance can impact the dynamics and structure of both the bound and unbound forms of the enzyme. Flap+ is a multi-drug-resistant variant of HIV-1 protease with a combination of mutations at the edge of the active site, within the active site, and in the flaps (L10I, G48V, I54V, V82A). The impact of these mutations on the dynamics in the unliganded form in comparison with the wild-type protease was elucidated with Molecular Dynamic simulations and NMR relaxation experiments. The comparative analyses from both methods concur in showing that the enzyme's dynamics are impacted by the drug resistance mutations in Flap+ protease. These alterations in the enzyme dynamics, particularly within the flaps, likely modulate the balance between substrate turn-over and drug binding, thereby conferring drug resistance.

WOS: 000455572800006, This study was to determine the severity of the tandoor burns in a city in Eastern Anatolia and to attract attention for the prevention of tandoor burns. Subjects and Methods: A retrospective evaluation was made of the registration data of burns hospitalized. A record was made of patient variables including burn type, age, gender, socioeconomic status, and amputation site and level. Statistical Analysis Used: Descriptive statistics were used to summarize the data, and results were reported as number and percentage. Categorical variables were evaluated with the Chi-square test. P