Contributions of Hyperactive Mutations in Mpro from SARS-CoV-2 to Drug Resistance.

The appearance and spread of mutations that cause drug resistance in rapidly evolving diseases, including infections by the SARS-CoV-2 virus, are major concerns for human health. Many drugs target enzymes, and resistance-conferring mutations impact inhibitor binding or enzyme activity. Nirmatrelvir, the most widely used inhibitor currently used to treat SARS-CoV-2 infections, targets the main protease (Mpro) preventing it from processing the viral polyprotein into active subunits. Our previous work systematically analyzed resistance mutations in Mpro that reduce binding to inhibitors; here, we investigate mutations that affect enzyme function. Hyperactive mutations that increase Mpro activity can contribute to drug resistance but have not been thoroughly studied. To explore how hyperactive mutations contribute to resistance, we comprehensively assessed how all possible individual mutations in Mpro affect enzyme function using a mutational scanning approach with a fluorescence resonance energy transfer (FRET)-based yeast readout. We identified hundreds of mutations that significantly increased the Mpro activity. Hyperactive mutations occurred both proximal and distal to the active site, consistent with protein stability and/or dynamics impacting activity. Hyperactive mutations were observed 3 times more than mutations which reduced apparent binding to nirmatrelvir in recent studies of laboratory-grown viruses selected for drug resistance. Hyperactive mutations were also about three times more prevalent than nirmatrelvir binding mutations in sequenced isolates from circulating SARS-CoV-2. Our findings indicate that hyperactive mutations are likely to contribute to the natural evolution of drug resistance in Mpro and provide a comprehensive list for future surveillance efforts.

Crystal Structures of Inhibitor-Bound Main Protease from Delta- and Gamma-Coronaviruses.

With the spread of SARS-CoV-2 throughout the globe causing the COVID-19 pandemic, the threat of zoonotic transmissions of coronaviruses (CoV) has become even more evident. As human infections have been caused by alpha- and beta-CoVs, structural characterization and inhibitor design mostly focused on these two genera. However, viruses from the delta and gamma genera also infect mammals and pose a potential zoonotic transmission threat. Here, we determined the inhibitor-bound crystal structures of the main protease (Mpro) from the delta-CoV porcine HKU15 and gamma-CoV SW1 from the beluga whale. A comparison with the apo structure of SW1 Mpro, which is also presented here, enabled the identification of structural arrangements upon inhibitor binding at the active site. The cocrystal structures reveal binding modes and interactions of two covalent inhibitors, PF-00835231 (active form of lufotrelvir) bound to HKU15, and GC376 bound to SW1 Mpro. These structures may be leveraged to target diverse coronaviruses and toward the structure-based design of pan-CoV inhibitors.

Allosteric quinoxaline-based inhibitors of the flavivirus NS2B/NS3 protease.

Viruses from the Flavivirus genus infect millions of people worldwide and cause severe diseases, including recent epidemics of dengue virus (DENV), and Zika virus (ZIKV). There is currently no antiviral treatment against flavivirus infections, despite considerable efforts to develop inhibitors against essential viral enzymes including NS2B/NS3 protease. Targeting the flavivirus NS2B/NS3 protease proved to be challenging because of the conformational dynamics, topology, and electrostatic properties of the active site. Here, we report the identification of quinoxaline-based allosteric inhibitors by fragment-based drug discovery approach as a promising new drug-like scaffold to target the NS2B/NS3 protease. Enzymatic assays and mutational analysis of the allosteric site in ZIKV NS2B/NS3 protease support noncompetitive inhibition mechanism as well as engineered DENV protease construct indicating the compounds likely compete with the NS2B cofactor for binding to the protease domain. Furthermore, antiviral activity confirmed the therapeutic potential of this new inhibitor scaffold.

Defining the substrate envelope of SARS-CoV-2 main protease to predict and avoid drug resistance.

Coronaviruses can evolve and spread rapidly to cause severe disease morbidity and mortality, as exemplified by SARS-CoV-2 variants of the COVID-19 pandemic. Although currently available vaccines remain mostly effective against SARS-CoV-2 variants, additional treatment strategies are needed. Inhibitors that target essential viral enzymes, such as proteases and polymerases, represent key classes of antivirals. However, clinical use of antiviral therapies inevitably leads to emergence of drug resistance. In this study we implemented a strategy to pre-emptively address drug resistance to protease inhibitors targeting the main protease (Mpro) of SARS-CoV-2, an essential enzyme that promotes viral maturation. We solved nine high-resolution cocrystal structures of SARS-CoV-2 Mpro bound to substrate peptides and six structures with cleavage products. These structures enabled us to define the substrate envelope of Mpro, map the critical recognition elements, and identify evolutionarily vulnerable sites that may be susceptible to resistance mutations that would compromise binding of the newly developed Mpro inhibitors. Our results suggest strategies for developing robust inhibitors against SARS-CoV-2 that will retain longer-lasting efficacy against this evolving viral pathogen.

Deciphering the Molecular Mechanism of HCV Protease Inhibitor Fluorination as a General Approach to Avoid 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.

Viral proteases: Structure, mechanism and inhibition.

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.

Pan-3C Protease Inhibitor Rupintrivir Binds SARS-CoV-2 Main Protease in a Unique Binding Mode.

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.

Discovery of Quinoxaline-Based P1-P3 Macrocyclic NS3/4A Protease Inhibitors with Potent Activity against Drug-Resistant Hepatitis C Virus Variants.

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.

Deciphering Antifungal Drug Resistance in Pneumocystis jirovecii DHFR with Molecular Dynamics and Machine Learning.

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.

Structural basis of substrate specificity in human cytidine deaminase family APOBEC3s.

The human cytidine deaminase family of APOBEC3s (A3s) plays critical roles in both innate immunity and the development of cancers. A3s comprise seven functionally overlapping but distinct members that can be exploited as nucleotide base editors for treating genetic diseases. Although overall structurally similar, A3s have vastly varying deamination activity and substrate preferences. Recent crystal structures of ssDNA-bound A3s together with experimental studies have provided some insights into distinct substrate specificities among the family members. However, the molecular interactions responsible for their distinct biological functions and how structure regulates substrate specificity are not clear. In this study, we identified the structural basis of substrate specificities in three catalytically active A3 domains whose crystal structures have been previously characterized: A3A, A3B- CTD, and A3G-CTD. Through molecular modeling and dynamic simulations, we found an interdependency between ssDNA substrate binding conformation and nucleotide sequence specificity. In addition to the U-shaped conformation seen in the crystal structure with the CTC0 motif, A3A can accommodate the CCC0 motif when ssDNA is in a more linear (L) conformation. A3B can also bind both U- and L-shaped ssDNA, unlike A3G, which can stably recognize only linear ssDNA. These varied conformations are stabilized by sequence-specific interactions with active site loops 1 and 7, which are highly variable among A3s. Our results explain the molecular basis of previously observed substrate specificities in A3s and have implications for designing A3-specific inhibitors for cancer therapy as well as engineering base-editing systems for gene therapy.

Deciphering Complex Mechanisms of Resistance and Loss of Potency through Coupled Molecular Dynamics and Machine Learning.

Drug resistance threatens many critical therapeutics through mutations in the drug target. The molecular mechanisms by which combinations of mutations, especially those remote 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 coupled parallel molecular dynamics simulations with experimental potency to identify specific conserved mechanisms underlying resistance. Physical features were extracted from the simulations, analyzed, and integrated into one consistent and interpretable elastic network model. To rigorously test this strategy, HIV-1 protease variants with diverse mutations were used, with potencies ranging from picomolar to micromolar to the drug darunavir. Feature reduction resulted in a model with four specific features that predicts for both the training and test sets inhibitor binding free energy within 1 kcal/mol of the experimental value over this entire range of potency. These predictive features are physically interpretable, as they vary specifically with affinity and diagonally transverse across the protease homodimer. This physics-based strategy of parallel molecular dynamics and machine learning captures mechanisms by which complex combinations of mutations confer resistance and identify critical features that serve as bellwethers of affinity, which will be critical in future drug design.

Inhibiting HTLV-1 Protease: A Viable Antiviral Target.

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 submicromolar 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.

Crystal Structure of SARS-CoV-2 Main Protease in Complex with the Non-Covalent Inhibitor ML188.

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 µ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 Å 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.

Drug Design Strategies to Avoid Resistance in Direct-Acting Antivirals and Beyond.

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.

Crystal Structure of a Soluble APOBEC3G Variant Suggests ssDNA to Bind in a Channel that Extends between the Two Domains.

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 Å 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.

A cross-reactive human IgA monoclonal antibody blocks SARS-CoV-2 spike-ACE2 interaction.

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.

Structural Analysis of Potent Hybrid HIV-1 Protease Inhibitors Containing Bis-tetrahydrofuran in a Pseudosymmetric Dipeptide Isostere.

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.

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

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.IMPORTANCE 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.

Optimizing the refinement of merohedrally twinned P61 HIV-1 protease-inhibitor cocrystal structures.

Twinning is a crystal-growth anomaly in which protein monomers exist in different orientations but are related in a specific way, causing diffraction reflections to overlap. Twinning imposes additional symmetry on the data, often leading to the assignment of a higher symmetry space group. Specifically, in merohedral twinning, reflections from each monomer overlap and require a twin law to model unique structural data from overlapping reflections. Neglecting twinning in the crystallographic analysis of quasi-rotationally symmetric homo-oligomeric protein structures can mask the degree of structural non-identity between monomers. In particular, any deviations from perfect symmetry will be lost if higher than appropriate symmetry is applied during crystallographic analysis. Such cases warrant choosing between the highest symmetry space group possible or determining whether the monomers have distinguishable structural asymmetries and thus require a lower symmetry space group and a twin law. Using hexagonal cocrystals of HIV-1 protease, a C2-symmetric homodimer whose symmetry is broken by bound ligand, it is shown that both assigning a lower symmetry space group and applying a twin law during refinement are critical to achieving a structural model that more accurately fits the electron density. By re-analyzing three recently published HIV-1 protease structures, improvements in nearly every crystallographic metric are demonstrated. Most importantly, a procedure is demonstrated where the inhibitor can be reliably modeled in a single orientation. This protocol may be applicable to many other homo-oligomers in the PDB.