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.

APOBEC3 enzymes are key enzymes in our innate immune system regulating antiviral response in HIV and unfortunately adding diversity in cancer as they deaminate cytosine. Seven unique single and double domain APOBEC3s provide them with unique activity and specificity profiles for this deamination. Recent crystal and NMR structures of APOBEC3 complexes are unraveling the variety of epitopes involved in binding nucleic acids, including at the catalytic site, elsewhere on the catalytic domain and in the inactive N-terminal domain. The interplay between these diverse interactions is critical to uncovering the mechanisms by which APOBEC3s recognize and process their substrates.

Drug resistance is prevalent across many diseases, rendering therapies ineffective with severe financial and health consequences. Rather than accepting resistance after the fact, proactive strategies need to be incorporated into the drug design and development process to minimize the impact of drug resistance. These strategies can be derived from our experience with viral disease targets where multiple generations of drugs had to be developed to combat resistance and avoid antiviral failure. Significant efforts including experimental and computational structural biology, medicinal chemistry, and machine learning have focused on understanding the mechanisms and structural basis of resistance against direct-acting antiviral (DAA) drugs. Integrated methods show promise for being predictive of resistance and potency. In this review, we give an overview of this research for human immunodeficiency virus type 1, hepatitis C virus, and influenza virus and the lessons learned from resistance mechanisms of DAAs. These lessons translate into rational strategies to avoid resistance in drug design, which can be generalized and applied beyond viral targets. While resistance may not be completely avoidable, rational drug design can and should incorporate strategies at the outset of drug development to decrease the prevalence of drug resistance.

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.

COVID-19 caused by SARS-CoV-2 has become a global pandemic requiring the development of interventions for the prevention or treatment to curtail mortality and morbidity. No vaccine to boost mucosal immunity, or as a therapeutic, has yet been developed to SARS-CoV-2. In this study, we discover and characterize a cross-reactive human IgA monoclonal antibody, MAb362. MAb362 binds to both SARS-CoV and SARS-CoV-2 spike proteins and competitively blocks ACE2 receptor binding, by overlapping the ACE2 structural binding epitope. Furthermore, MAb362 IgA neutralizes both pseudotyped SARS-CoV and SARS-CoV-2 in 293 cells expressing ACE2. When converted to secretory IgA, MAb326 also neutralizes authentic SARS-CoV-2 virus while the IgG isotype shows no neutralization. Our results suggest that SARS-CoV-2 specific IgA antibodies, such as MAb362, may provide effective immunity against SARS-CoV-2 by inducing mucosal immunity within the respiratory system, a potentially critical feature of an effective vaccine.

Type 1 diabetes (T1D) is caused by the autoimmune destruction of pancreatic beta cells. Pluripotent stem cells can now be differentiated into beta cells, thus raising the prospect of a cell replacement therapy for T1D. However, autoimmunity would rapidly destroy newly transplanted beta cells. Using a genome-scale CRISPR screen in a mouse model for T1D, we show that deleting RNLS, a genome-wide association study candidate gene for T1D, made beta cells resistant to autoimmune killing. Structure-based modelling identified the U.S. Food and Drug Administration-approved drug pargyline as a potential RNLS inhibitor. Oral pargyline treatment protected transplanted beta cells in diabetic mice, thus leading to disease reversal. Furthermore, pargyline prevented or delayed diabetes onset in several mouse models for T1D. Our results identify RNLS as a modifier of beta cell vulnerability and as a potential therapeutic target to avert beta cell loss in T1D.

Over the last few decades, protein NMR isotope labeling methods using E. coli based expression have revolutionized the information accessible from biomolecular NMR experiments. Selective labeling of a protein of interest in a multi-protein complex can significantly reduce the number of cross-peaks and allow for study of large protein complexes. However, limitations still remain since some proteins are not stable independently and cannot be separately labeled in either NMR active isotope enriched or unenriched media and reconstituted into a multimeric complex. To overcome this limitation, the LEGO-NMR method was previously developed using protein expression plasmids containing T7 or araBAD promoters to separately express proteins in the same E. coli after changing between labeled and unlabeled media. Building on this, we developed a method to label the Human Immunodeficiency Virus type 1 viral infectivity factor (HIV-1 Vif), a monomerically unstable protein, in complex with CBFβ, it’s host binding partner. We designed a dual promoter plasmid containing both T7 and araBAD promoters to independently control the expression of HIV-1 Vif in NMR active isotope enriched media and CBFβ in unenriched media. Using this method, we assigned the backbone resonance and directly observed the binding of HIV-1 Vif with APOBEC3G, a host restriction factor to HIV-1.

SummaryInfluenza viruses pose severe public health threats; they cause millions of infections and tens of thousands of deaths annually in the US. Influenza viruses are extensively pleomorphic, in both shape and size as well as organization of viral structural proteins. Analysis of influenza morphology and ultrastructure can help elucidate viral structure-function relationships as well as aid in therapeutics and vaccine development. While cryo-electron tomography (cryoET) can depict the 3D organization of pleomorphic influenza, the low signal-to-noise ratio inherent to cryoET and extensive viral heterogeneity have precluded detailed characterization of influenza viruses. In this report, we developed a cryoET processing pipeline leveraging convolutional neural networks (CNNs) to characterize the morphological architecture of the A/Puerto Rico/8/34 (H1N1) influenza strain. Our pipeline improved the throughput of cryoET analysis and accurately identified viral components within tomograms. Using this approach, we successfully characterized influenza viral morphology, glycoprotein density, and conduct subtomogram averaging of HA glycoproteins. Application of this processing pipeline can aid in the structural characterization of not only influenza viruses, but other pleomorphic viruses and infected cells.Graphical abstract

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.

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.

Antibodies elicited by infection accumulate somatic mutations in germinal centers that can increase affinity for cognate antigens. We analyzed 6 independent groups of clonally related severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) Spike receptor-binding domain (RBD)-specific antibodies from 5 individuals shortly after infection and later in convalescence to determine the impact of maturation over months. In addition to increased affinity and neutralization potency, antibody evolution changed the mutational pathways for the acquisition of viral resistance and restricted neutralization escape options. For some antibodies, maturation imposed a requirement for multiple substitutions to enable escape. For certain antibodies, affinity maturation enabled the neutralization of circulating SARS-CoV-2 variants of concern and heterologous sarbecoviruses. Antibody-antigen structures revealed that these properties resulted from substitutions that allowed additional variability at the interface with the RBD. These findings suggest that increasing antibody diversity through prolonged or repeated antigen exposure may improve protection against diversifying SARS-CoV-2 populations, and perhaps against other pandemic threat coronaviruses.

Drug resistance impacts the effectiveness of many new therapeutics. Mutations in the therapeutic target confer resistance, however deciphering which mutations, often remote from the enzyme active site, drive resistance is challenging. In a series of Pneumocystis Jirovecii dihydrofolate reductase variants we elucidate which interactions are key bellwethers to confer resistance to trimethoprim using homology modeling, molecular dynamics and machine learning. Six molecular features involving mainly residues that did not vary, were the best indicators of resistance.

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.

Discovery and optimization of small molecule inhibitors as therapeutic drugs have immensely benefited from rational structure-based drug design. With recent advances in high-resolution structure determination, computational power, and machine learning methodology, it is becoming more tractable to elucidate the structural basis of drug potency. However, the applicability of machine learning models to drug design is limited by the interpretability of the resulting models in terms of feature importance. Here, we take advantage of the large number of available inhibitor-bound HIV-1 protease structures and associated potencies to evaluate inhibitor diversity and machine learning models to predict ligand affinity. First, using a hierarchical clustering approach, we grouped HIV-1 protease inhibitors and identified distinct core structures. Explicit features including protein-ligand interactions were extracted from high-resolution cocrystal structures as 3D-based fingerprints. We found that a gradient boosting machine learning model with this explicit feature attribution can predict binding affinity with high accuracy. Finally, Shapley values were derived to explain local feature importance. We found specific van der Waals (vdW) interactions of key protein residues are pivotal for the predicted potency. Protein-specific and interpretable prediction models can guide the optimization of many small molecule drugs for improved potency.

Direct acting antivirals have dramatically increased the efficacy and tolerability of hepatitis C treatment, but drug resistance has emerged with some of these inhibitors, including non-structural protein 3/4A protease inhibitors (PIs). Although many co-crystal structures of PIs with the NS3/4A protease have been reported, a systematic review of these crystal structures in the context of the rapidly emerging drug resistance especially for early PIs has not been performed. To provide a framework for designing better inhibitors with higher barriers to resistance, we performed a quantitative structural analysis using co-crystal structures and models of HCV NS3/4A protease in complex with natural substrates and inhibitors. By comparing substrate structural motifs and active site interactions with inhibitor recognition, we observed that the selection of drug resistance mutations correlates with how inhibitors deviate from viral substrates in molecular recognition. Based on this observation, we conclude that guiding the design process with native substrate recognition features is likely to lead to more robust small molecule inhibitors with decreased susceptibility to resistance.

Retroviral proteases (PR) have a unique specificity that allows cleavage of sites with or without a P1’ proline. A P1’ proline is required at the MA/CA cleavage site due to its role in a post-cleavage conformational change in the capsid protein. However, the HIV-1 PR prefers to have large hydrophobic amino acids flanking the scissile bond, suggesting PR recognizes two different classes of substrate sequences. We analyzed the cleavage rate of over 150 iterations of six different HIV-1 cleavage sites to explore rate determinants of cleavage. We found that cleavage rates are strongly influenced by the two amino acids flanking the amino acids at the scissile bond (P2-P1/P1’-P2’), with two complementary sets of rules. When P1’ is proline, the P2 side chain interacts with a polar region in the S2 subsite of the PR, while the P2’ amino acid interacts with a hydrophobic region of the S2’ subsite. When P1’ is not proline, the orientations of the P2 and P2’ side chains with respect to the scissile bond are reversed; P2 residues interact with a hydrophobic face of the S2 subsite while the P2’ amino acid usually engages hydrophilic amino acids in the S2’ subsite. These results reveal that the HIV-1 PR has evolved bi-functional S2 and S2’ subsites to accommodate the steric effects imposed by a P1’ proline on the orientation of P2 and P2’ substrate side chains. These results also suggest a new strategy for inhibitor design to engage the multiple specificities in these subsites.

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.

Antibodies elicited by infection accumulate somatic mutations in germinal centers that can increase affinity for cognate antigens. We analyzed 6 independent groups of clonally related severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) Spike receptor-binding domain (RBD)-specific antibodies from 5 individuals shortly after infection and later in convalescence to determine the impact of maturation over months. In addition to increased affinity and neutralization potency, antibody evolution changed the mutational pathways for the acquisition of viral resistance and restricted neutralization escape options. For some antibodies, maturation imposed a requirement for multiple substitutions to enable escape. For certain antibodies, affinity maturation enabled the neutralization of circulating SARS-CoV-2 variants of concern and heterologous sarbecoviruses. Antibody-antigen structures revealed that these properties resulted from substitutions that allowed additional variability at the interface with the RBD. These findings suggest that increasing antibody diversity through prolonged or repeated antigen exposure may improve protection against diversifying SARS-CoV-2 populations, and perhaps against other pandemic threat coronaviruses., Graphical abstract, SARS-CoV-2 RBD-binding antibodies accumulate mutations over months of convalescence, but how maturation affects antibody properties is unclear. Muecksch et al. analyze 6 groups of clonally related neutralizing antibodies from convalescent donors and reveal that antibody maturation increases affinity, potency, and breadth, and restricts viral escape options. They also identify the structural basis for these improved antibody properties during maturation.

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.

APOBEC3G (A3G) is a single-stranded DNA (ssDNA) cytosine deaminase that can restrict HIV-1 infection by mutating the viral genome. A3G consists of a non-catalytic N-terminal domain (NTD) and a catalytic C-terminal domain (CTD) connected by a short linker. While the CTD catalyzes cytosine deamination, the NTD is believed to provide additional affinity for ssDNA. Structures of both A3G domains have been solved individually; however, a full-length A3G structure has been challenging. Recently, crystal structures of full-length rhesus macaque A3G variants were solved which suggested dimerization mechanisms and RNA binding surfaces, whereas the dimerization appeared to compromise catalytic activity. We determined the crystal structure of a soluble variant of human A3G (sA3G) at 2.5 A and from these data generated a model structure of wild-type A3G. This model demonstrated that the NTD was rotated 90° relative to the CTD along the major axis of the molecule, an orientation that forms a positively charged channel connected to the CTD catalytic site, consisting of NTD loop-1 and CTD loop-3. Structure-based mutations, in vitro deamination and DNA binding assays, and HIV-1 restriction assays identify R24, located in the NTD loop-1, as essential to a critical interaction with ssDNA. Furthermore, sA3G was shown to bind a deoxy-cytidine dinucleotide near the catalytic Zn2+, yet not in the catalytic position, where the interactions between deoxy-cytidines and CTD loop-1 and loop-7 residues were different from those formed with substrate. These new interactions suggest a mechanism explaining why A3G exhibits a 3′ to 5′ directional preference in processive deamination.

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.

COVID-19 caused by SARS-CoV-2 has become a global pandemic requiring the development of interventions for the prevention or treatment to curtail mortality and morbidity. No vaccine to boost mucosal immunity, or as a therapeutic, has yet been developed to SARS-CoV-2. In this study, we discover and characterize a cross-reactive human IgA monoclonal antibody, MAb362. MAb362 binds to both SARS-CoV and SARS-CoV-2 spike proteins and competitively blocks ACE2 receptor binding, by overlapping the ACE2 structural binding epitope. Furthermore, MAb362 IgA neutralizes both pseudotyped SARS-CoV and SARS-CoV-2 in 293 cells expressing ACE2. When converted to secretory IgA, MAb326 also neutralizes authentic SARS-CoV-2 virus while the IgG isotype shows no neutralization. Our results suggest that SARS-CoV-2 specific IgA antibodies, such as MAb362, may provide effective immunity against SARS-CoV-2 by inducing mucosal immunity within the respiratory system, a potentially critical feature of an effective vaccine., Here, Ejemel et al. report the identification and characterization of a cross-neutralizing human IgA monoclonal antibody, named MAb362, that binds the receptor-binding domain of SARS-CoV-2 Spike, blocking its interaction with the ACE2 host receptor.

Broadly neutralizing antibodies (bnAbs) targeting conserved influenza A virus (IAV) hemagglutinin (HA) epitopes can provide valuable information for accelerating universal vaccine designs. Here, we report structural details for heterosubtypic recognition of HA from circulating and emerging IAVs by the human antibody 3I14. Somatic hypermutations play a critical role in shaping the HCDR3, which alone and uniquely among VH3-30 derived antibodies, forms contacts with five sub-pockets within the HA-stem hydrophobic groove. 3I14 light-chain interactions are also key for binding HA and contribute a large buried surface area spanning two HA protomers. Comparison of 3I14 to bnAbs from several defined classes provide insights to the bias selection of VH3-30 antibodies and reveals that 3I14 represents a novel structural solution within the VH3-30 repertoire. The structures reported here improve our understanding of cross-group heterosubtypic binding activity, providing the basis for advancing immunogen designs aimed at eliciting a broadly protective response to IAV., Previously, a broadly neutralizing antibody, 3I14, active against groups 1 and 2 influenza A viruses was isolated from human memory B cells and showed protection in mice from lethal viral challenge. Here, Harshbarger and Deming et al. provide the crystal structure of 3I14 Fab in complex with H3, H6, and H10.

Over the past several decades, atomistic simulations of biomolecules, whether carried out using molecular dynamics or Monte Carlo techniques, have provided detailed insights into their function. Comparing the results of such simulations for a few closely related systems has guided our understanding of the mechanisms by which changes such as ligand binding or mutation can alter the function. The general problem of detecting and interpreting such mechanisms from simulations of many related systems, however, remains a challenge. This problem is addressed here by applying supervised and unsupervised machine learning techniques to a variety of thermodynamic observables extracted from molecular dynamics simulations of different systems. As an important test case, these methods are applied to understand the evasion by human immunodeficiency virus type-1 (HIV-1) protease of darunavir, a potent inhibitor to which resistance can develop via the simultaneous mutation of multiple amino acids. Complex mutational patterns have been observed among resistant strains, presenting a challenge to developing a mechanistic picture of resistance in the protease. In order to dissect these patterns and gain mechanistic insight into the role of specific mutations, molecular dynamics simulations were carried out on a collection of HIV-1 protease variants, chosen to include highly resistant strains and susceptible controls, in complex with darunavir. Using a machine learning approach that takes advantage of the hierarchical nature in the relationships among the sequence, structure, and function, an integrative analysis of these trajectories reveals key details of the resistance mechanism, including changes in the protein structure, hydrogen bonding, and protein-ligand contacts.

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.

The human APOBEC3G protein is a cytidine deaminase that generates cytidine to deoxy-uridine mutations in single-stranded DNA (ssDNA), and capable of restricting replication of HIV-1 by generating mutations in viral genome. The mechanism by which APOBEC3G specifically deaminates 5′-CC motifs has remained elusive since structural studies have been hampered due to apparently weak ssDNA binding of the catalytic domain of APOBEC3G. We overcame the problem by generating a highly active variant with higher ssDNA affinity. Here, we present the crystal structure of this variant complexed with a ssDNA substrate at 1.86 Å resolution. This structure reveals atomic-level interactions by which APOBEC3G recognizes a functionally-relevant 5′-TCCCA sequence. This complex also reveals a key role of W211 in substrate recognition, implicating a similar recognition in activation-induced cytidine deaminase (AID) with a conserved tryptophan., APOBEC3G (A3G) is a single-stranded DNA (ssDNA) cytidine deaminase that restricts HIV-1. Here the authors provide molecular insights into A3G substrate recognition by determining the 1.86 Å resolution crystal structure of its catalytic domain bound to ssDNA.

The APOBEC3 (A3) family of human cytidine deaminases is renowned for providing a first line of defense against many exogenous and endogenous retroviruses. However, the ability of these proteins to deaminate deoxycytidines in ssDNA makes A3s a double-edged sword. When overexpressed, A3s can mutate endogenous genomic DNA resulting in a variety of cancers. Although the sequence context for mutating DNA varies among A3s, the mechanism for substrate sequence specificity is not well understood. To characterize substrate specificity of A3A, a systematic approach was used to quantify the affinity for substrate as a function of sequence context, length, secondary structure, and solution pH. We identified the A3A ssDNA binding motif as (T/C)TC(A/G), which correlated with enzymatic activity. We also validated that A3A binds RNA in a sequence specific manner. A3A bound tighter to substrate binding motif within a hairpin loop compared to linear oligonucleotide, suggesting A3A affinity is modulated by substrate structure. Based on these findings and previously published A3A–ssDNA co-crystal structures, we propose a new model with intra-DNA interactions for the molecular mechanism underlying A3A sequence preference. Overall, the sequence and structural preferences identified for A3A leads to a new paradigm for identifying A3A’s involvement in mutation of endogenous or exogenous DNA.

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.

In the era of state-of-the-art inhibitor design and high-resolution structural studies, detection of significant but small protein structural differences in the inhibitor-bound forms is critical to further developing the inhibitor. Here, we probed differences in HIV-1 protease (PR) conformation among darunavir and four analogous inhibitor-bound forms and compared them with a drug-resistant mutant using nuclear magnetic resonance chemical shifts. Changes in amide chemical shifts of wild-type (WT) PR among these inhibitor-bound forms, ΔCSP, were subtle but detectable and extended10 Å from the inhibitor-binding site, asymmetrically between the two subunits of PR. Molecular dynamics simulations revealed differential local hydrogen bonding as the molecular basis of this remote asymmetric change. Inhibitor-bound forms of the drug-resistant mutant also showed a similar long-range ΔCSP pattern. Differences in ΔCSP values of the WT and the mutant (ΔΔCSPs) were observed at the inhibitor-binding site and in the surrounding region. Comparing chemical shift changes among highly analogous inhibitors and ΔΔCSPs effectively eliminated local environmental effects stemming from different chemical groups and enabled exploitation of these sensitive parameters to detect subtle protein conformational changes and to elucidate asymmetric and remote conformational effects upon inhibitor interaction.

Summary The heavy chain IGHV1-69 germline gene exhibits a high level of polymorphism and shows biased use in protective antibody (Ab) responses to infections and vaccines. It is also highly expressed in several B cell malignancies and autoimmune diseases. G6 is an anti-idiotypic monoclonal Ab that selectively binds to IGHV1-69 heavy chain germline gene 51p1 alleles that have been implicated in these Ab responses and disease processes. Here, we determine the co-crystal structure of humanized G6 (hG6.3) in complex with anti-influenza hemagglutinin stem-directed broadly neutralizing Ab D80. The core of the hG6.3 idiotope is a continuous string of CDR-H2 residues starting with M53 and ending with N58. G6 binding studies demonstrate the remarkable breadth of binding to 51p1 IGHV1-69 Abs with diverse CDR-H3, light chain, and antigen binding specificities. These studies detail the broad expression of the G6 cross-reactive idiotype (CRI) that further define its potential role in precision medicine., Graphical Abstract, Highlights • G6 binds to a subset of IGHV1-69 germline-based anti-influenza Abs • The structure of humanized G6 with a IGHV1-69 anti-influenza Ab is reported • Various binding assays further define the G6 cross-reactive binding idiotope • The core binding idiotope of G6 is deduced, G6 is an exceptional anti-idiotypic antibody that binds to antibodies encoded by the immunoglobulin heavy chain germline gene IGHV1-69. Avnir et al. describe how G6 binds to its cross-reactive idiotope by a set of binding assays and by solving the structure of humanized-G6 complexed with an IGHV1-69-anti-influenza Ab.

HIV-1 protease is responsible for the cleavage of 12 non-homologous sites within the Gag and Gag-Pro-Pol polyproteins in the viral genome. Under the selective pressure of protease inhibition, the virus evolves mutations within (primary) and outside of (secondary) the active site allowing the protease to process substrates while simultaneously countering inhibition. The primary protease mutations impede inhibitor binding directly, while the secondary mutations are considered accessory mutations that compensate for a loss in fitness. However, the role of secondary mutations in conferring drug resistance remains a largely unresolved topic. We have shown previously that mutations distal to the active site are able to perturb binding of darunavir (DRV) via the protein’s internal hydrogen-bonding network. In this study we show that mutations distal to the active site, regardless of context, can play an interdependent role in drug resistance. Applying eigenvalue decomposition to collections of hydrogen bonding and van der Waals interactions from a series of molecular dynamics simulations of 15 diverse HIV-1 protease variants, we identify sites in the protease where amino acid substitutions lead to perturbations in non-bonded interactions with DRV and/or the hydrogen-bonding network of the protease itself. While primary mutations are known to drive resistance in HIV-1 protease, these findings delineate the significant contributions of accessory mutations to resistance. Identifying the variable positions in the protease that have the greatest impact on drug resistance may aid in future structure-based design of inhibitors.

A substrate envelope-guided design strategy is reported for improving the resistance profile of HCV NS3/4A protease inhibitors. Analogues of 5172-mcP1P3 were designed by incorporating diverse quinoxalines at the P2 position that predominantly interact with the invariant catalytic triad of the protease. Exploration of structure-activity relationships showed that inhibitors with small hydrophobic substituents at the 3-position of P2 quinoxaline maintain better potency against drug resistant variants, likely due to reduced interactions with residues in the S2 subsite. In contrast, inhibitors with larger groups at this position were highly susceptible to mutations at Arg155, Ala156 and Asp168. Excitingly, several inhibitors exhibited exceptional potency profiles with EC50 values ≤ 5 nM against major drug resistant HCV variants. These findings support that inhibitors designed to interact with evolutionarily constrained regions of the protease, while avoiding interactions with residues not essential for substrate recognition, are less likely to be susceptible to drug 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.

C-terminal binding proteins 1 and 2 (CtBP1 and CtBP2) are transcriptional regulators that activate or repress many genes involved in cellular development, apoptosis, and metastasis. NADH-dependent CtBP activation has been implicated in multiple types of cancer and poor patient prognosis. Central to understanding activation of CtBP in oncogenesis is uncovering how NADH triggers protein assembly, what level of assembly occurs, and if oncogenic activity depends upon such assembly. Here, we present the cryoelectron microscopic structures of two different constructs of CtBP2 corroborating that the native state of CtBP2 in the presence of NADH is tetrameric. The physiological relevance of the observed tetramer was demonstrated in cell culture, showing that CtBP tetramer-destabilizing mutants are defective for cell migration, transcriptional repression of E-cadherin, and activation of TIAM1. Together with our cryoelectron microscopy studies, these results highlight the tetramer as the functional oligomeric form of CtBP2.

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.

This Editorial discusses the state of research on drug resistance in the fields of cancer, infectious disease, and agriculture. Reaching across the aisle for a more cross-collaborative approach may lead to exciting breakthroughs toward tackling the challenges of drug resistance in each field.

APOBEC3s proteins (A3s), a family of human cytidine deaminases, protect the host cell from endogenous retro-elements and exogenous viral infections by introducing hypermutations. However, the ability to mutate genomic DNA makes A3s a potential cancer source. Of the 7 human A3s, A3B has been implicated as an endogenous cause for multiple cancers. Despite overall similarity, A3s have distinct deamination activity with A3B among the least catalytically active. Over the past few years, several structures of apo as well as DNA-bound A3 proteins have been determined. These structures revealed the molecular determinants of nucleotide specificity and the importance of the loops around the active site in DNA binding. However, for A3B, the structural basis for regulation of deamination activity and the role of active site loops in coordinating DNA had remained unknown. In this study, using a combination of advanced molecular modelling followed by experimental mutational analysis and dynamics simulations, we investigated molecular mechanism of A3B regulating activity and DNA binding. We identified a unique auto-inhibited conformation of A3B that restricts access and binding of DNA to the active site, mainly due to the extra PLV residues in loop 1. We modelled DNA binding to fully native A3B and found that Arg211 in the arginine patch of loop1 is the gatekeeper while Arg212 stabilizes the bound DNA. This model also identified the critical residues for substrate specificity, especially at the -1 position. Our results reveal the structural basis for relatively lower catalytic activity of A3B and provide opportunities for rational design of inhibitors that specifically target A3B to benefit cancer therapeutics.

Though the steps of human immunodeficiency virus type 1 (HIV-1) virion maturation are well documented, the mechanisms regulating the proteolysis of the Gag and Gag-Pro-Pol polyproteins by the HIV-1 protease (PR) remain obscure. One proposed mechanism argues that the maturation intermediate p15NC must interact with RNA for efficient cleavage by the PR. We investigated this phenomenon and found that processing of multiple substrates by the HIV-1 PR was enhanced in the presence of RNA. The acceleration of proteolysis occurred independently from the substrate's ability to interact with nucleic acid, indicating that a direct interaction between substrate and RNA is not necessary for enhancement. Gel-shift assays demonstrated the HIV-1 PR is capable of interacting with nucleic acids, suggesting that RNA accelerates processing reactions by interacting with the PR rather than the substrate. All HIV-1 PRs examined have this ability; however, the HIV-2 PR does not interact with RNA and does not exhibit enhanced catalytic activity in the presence of RNA. No specific sequence or structure was required in the RNA for a productive interaction with the HIV-1 PR, which appears to be principally, though not exclusively, driven by electrostatic forces. For a peptide substrate, RNA increased the kinetic efficiency of the HIV-1 PR by an order of magnitude, affecting both turnover rate ( k cat ) and substrate affinity ( K m ). These results suggest that an allosteric binding site exists on the HIV-1 PR and that HIV-1 PR activity during maturation could be regulated in part by the juxtaposition of the enzyme with virion-packaged RNA.

The human APOBEC3G (A3G) DNA cytosine deaminase restricts and hypermutates DNA-based parasites including HIV-1. The viral infectivity factor (Vif) prevents restriction by triggering A3G degradation. Although the structure of the A3G catalytic domain is known, the structure of the N-terminal Vif-binding domain has proven more elusive. Here, we used evolution- and structure-guided mutagenesis to solubilize the Vif-binding domain of A3G, thus permitting structural determination by NMR spectroscopy. A smaller zinc-coordinating pocket and altered helical packing distinguish the structure from previous catalytic-domain structures and help to explain the reported inactivity of this domain. This soluble A3G N-terminal domain is bound by Vif; this enabled mutagenesis and biochemical experiments, which identified a unique Vif-interacting surface formed by the α1-β1, β2-α2 and β4-α4 loops. This structure sheds new light on the Vif-A3G interaction and provides critical information for future drug development.

Nucleic acid editing enzymes are essential components of the immune system that lethally mutate viral pathogens and somatically mutate immunoglobulins, and contribute to the diversification and lethality of cancers. Among these enzymes are the seven human APOBEC3 deoxycytidine deaminases, each with unique target sequence specificity and subcellular localization. While the enzymology and biological consequences have been extensively studied, the mechanism by which APOBEC3s recognize and edit DNA remains elusive. Here we present the crystal structure of a complex of a cytidine deaminase with ssDNA bound in the active site at 2.2 Å. This structure not only visualizes the active site poised for catalysis of APOBEC3A, but pinpoints the residues that confer specificity towards CC/TC motifs. The APOBEC3A–ssDNA complex defines the 5′–3′ directionality and subtle conformational changes that clench the ssDNA within the binding groove, revealing the architecture and mechanism of ssDNA recognition that is likely conserved among all polynucleotide deaminases, thereby opening the door for the design of mechanistic-based therapeutics., Cytidine deaminases are evolutionarily conserved enzymes that edit genomes by deaminating cytidine to uridine. Here the authors present the crystal structure of APOBEC3A with a single-stranded DNA substrate bound in the active site to shed light on the mechanism and specificity of substrate recognition.

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 are part of the highly active antiretroviral therapy effectively used in the treatment of HIV infection and AIDS. Darunavir (DRV) is the most potent of these inhibitors, soliciting drug resistance only when a complex combination of mutations occur both inside and outside the protease active site. With few exceptions, the role of mutations outside the active site in conferring resistance remains largely elusive. Through a series of DRV-protease complex crystal structures, inhibition assays, and molecular dynamics simulations, we find that single and double site mutations outside the active site often associated with DRV resistance alter the structure and dynamic ensemble of HIV-1 protease active site. These alterations correlate with the observed inhibitor binding affinities for the mutants, and suggest a network hypothesis on how the effect of distal mutations are propagated to pivotal residues at the active site and may contribute to conferring drug 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 oncogenic corepressors C-terminal Binding Protein (CtBP) 1 and 2 harbor regulatory d-isomer specific 2-hydroxyacid dehydrogenase (d2-HDH) domains. 4-Methylthio 2-oxobutyric acid (MTOB) exhibits substrate inhibition and can interfere with CtBP oncogenic activity in cell culture and mice. Crystal structures of human CtBP1 and CtBP2 in complex with MTOB and NAD+ revealed two key features: a conserved tryptophan that likely contributes to substrate specificity and a hydrophilic cavity that links MTOB with an NAD+ phosphate. Neither feature is present in other d2-HDH enzymes. These structures thus offer key opportunities for the development of highly selective anti-neoplastic CtBP inhibitors.

Here we present a novel, end-point method using the dead-end-elimination and A* algorithms to efficiently and accurately calculate the change in free energy, enthalpy, and configurational entropy of binding for ligand− receptor association reactions. We apply the new approach to the binding of a series of human immunodeficiency virus (HIV-1) protease inhibitors to examine the effect ensemble reranking has on relative accuracy as well as to evaluate the role of the absolute and relative ligand configurational entropy losses upon binding in affinity differences for structurally related inhibitors. Our results suggest that most thermodynamic parameters can be estimated using only a small fraction of the full configurational space, and we see significant improvement in relative accuracy when using an ensemble versus single-conformer approach to ligand ranking. We also find that using approximate metrics based on the single-conformation enthalpy differences between the global minimum energy configuration in the bound as well as unbound states also correlates well with experiment. Using a novel, additive entropy expansion based on conditional mutual information, we also analyze the source of ligand configurational entropy loss upon binding in terms of both uncoupled per degree of freedom losses as well as changes in coupling between inhibitor degrees of freedom. We estimate entropic free energy losses of approximately +24 kcal/mol, 12 kcal/mol of which stems from loss of translational and rotational entropy. Coupling effects contribute only a small fraction to the overall entropy change (1−2 kcal/mol) but suggest differences in how inhibitor dihedral angles couple to each other in the bound versus unbound states. The importance of accounting for flexibility in drug optimization and design is also discussed.

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)