Comprehensive fitness landscape of SARS-CoV-2 Mpro reveals insights into viral resistance mechanisms.

With the continual evolution of new strains of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) that are more virulent, transmissible, and able to evade current vaccines, there is an urgent need for effective anti-viral drugs. The SARS-CoV-2 main protease (Mpro) is a leading target for drug design due to its conserved and indispensable role in the viral life cycle. Drugs targeting Mpro appear promising but will elicit selection pressure for resistance. To understand resistance potential in Mpro, we performed a comprehensive mutational scan of the protease that analyzed the function of all possible single amino acid changes. We developed three separate high throughput assays of Mpro function in yeast, based on either the ability of Mpro variants to cleave at a defined cut-site or on the toxicity of their expression to yeast. We used deep sequencing to quantify the functional effects of each variant in each screen. The protein fitness landscapes from all three screens were strongly correlated, indicating that they captured the biophysical properties critical to Mpro function. The fitness landscapes revealed a non-active site location on the surface that is extremely sensitive to mutation, making it a favorable location to target with inhibitors. In addition, we found a network of critical amino acids that physically bridge the two active sites of the Mpro dimer. The clinical variants of Mpro were predominantly functional in our screens, indicating that Mpro is under strong selection pressure in the human population. Our results provide predictions of mutations that will be readily accessible to Mpro evolution and that are likely to contribute to drug resistance. This complete mutational guide of Mpro can be used in the design of inhibitors with reduced potential of evolving viral resistance.

Donor-Acceptor Pyridinium Salts for Photo-Induced Electron-Transfer-Driven Modification of Tryptophan in Peptides, Proteins, and Proteomes Using Visible Light.

Tryptophan (Trp) plays a variety of critical functional roles in protein biochemistry; however, owing to its low natural frequency and poor nucleophilicity, the design of effective methods for both single protein bioconjugation at Trp as well as for in situ chemoproteomic profiling remains a challenge. Here, we report a method for covalent Trp modification that is suitable for both scenarios by invoking photo-induced electron transfer (PET) as a means of driving efficient reactivity. We have engineered biaryl N-carbamoyl pyridinium salts that possess a donor-acceptor relationship that enables optical triggering with visible light whilst simultaneously attenuating the probe's photo-oxidation potential in order to prevent photodegradation. This probe was assayed against a small bank of eight peptides and proteins, where it was found that micromolar concentrations of the probe and short irradiation times (10-60 min) with violet light enabled efficient reactivity toward surface exposed Trp residues. The carbamate transferring group can be used to transfer useful functional groups to proteins including affinity tags and click handles. DFT calculations and other mechanistic analyses reveal correlations between excited state lifetimes, relative fluorescence quantum yields, and chemical reactivity. Biotinylated and azide-functionalized pyridinium salts were used for Trp profiling in HEK293T lysates and in situ in HEK293T cells using 440 nm LED irradiation. Peptide-level enrichment from live cell labeling experiments identified 290 Trp modifications, with 82% selectivity for Trp modification over other π-amino acids, demonstrating the ability of this method to identify and quantify reactive Trp residues from live cells.

Morphological Characterization of Antibiotic Combinations.

Combination therapies are common in many therapeutic contexts, including infectious diseases and cancer. A common approach for evaluating combinations in vitro is to assess effects on cell growth as synergistic, antagonistic, or neutral using "checkerboard" experiments to systematically sample combinations of agents in multiple doses. To further understand the effects of antibiotic combinations, we employed high-content imaging to study the morphological changes caused by combination treatments in checkerboard experiments. Using an automated, unsupervised image analysis approach to group morphologies, and an "expert-in-the-loop" to annotate them, we attributed the heterogeneous morphological changes ofEscherichia coli cells to varying doses of both single-agent and combination treatments. We identified patterns of morphological change, including morphological potentiation, competition, and the emergence of unexpected morphologies. We found these frequently did not correlate with synergistic or antagonistic effects on viability, suggesting morphological approaches may provide a distinctive signature of the biological interaction between compounds over a range of conditions. Among the unexpected morphologies we observed, there were transitional changes associated with intermediate doses of compounds and uncharacterized phenotypes associated with well-studied antibiotics. Our approach exemplifies how unsupervised image analysis and expert knowledge can be combined to reckon with complex phenotypic changes arising from combination screening, dose titrations, or polypharmacology. In this way, quantification of morphological diversity across treatment conditions could aid in analysis and prioritization of complementary pairings of antibiotic agents by more closely capturing the true signature of the integrated cellular response.

NITD-688, a pan-serotype inhibitor of the dengue virus NS4B protein, shows favorable pharmacokinetics and efficacy in preclinical animal models.

Dengue virus (DENV) is a mosquito-borne flavivirus that poses a threat to public health, yet no antiviral drug is available. We performed a high-throughput phenotypic screen using the Novartis compound library and identified candidate chemical inhibitors of DENV. This chemical series was optimized to improve properties such as anti-DENV potency and solubility. The lead compound, NITD-688, showed strong potency against all four serotypes of DENV and demonstrated excellent oral efficacy in infected AG129 mice. There was a 1.44-log reduction in viremia when mice were treated orally at 30 milligrams per kilogram twice daily for 3 days starting at the time of infection. NITD-688 treatment also resulted in a 1.16-log reduction in viremia when mice were treated 48 hours after infection. Selection of resistance mutations and binding studies with recombinant proteins indicated that the nonstructural protein 4B is the target of NITD-688. Pharmacokinetic studies in rats and dogs showed a long elimination half-life and good oral bioavailability. Extensive in vitro safety profiling along with exploratory rat and dog toxicology studies showed that NITD-688 was well tolerated after 7-day repeat dosing, demonstrating that NITD-688 may be a promising preclinical candidate for the treatment of dengue.

Mapping the Heterogeneity of Histone Modifications on Hepatitis B Virus DNA Using Liver Needle Biopsies Obtained from Chronically Infected Patients.

Covalently closed circular DNA (cccDNA) forms the basis for replication and persistence of hepatitis B virus (HBV) in the chronically infected liver. We have previously shown that viral transcription is subject to regulation by posttranslational modifications (PTMs) of histone proteins bound to cccDNA through analysis of de novo HBV-infected cell lines. We now report the successful adaptation of this chromatin immunoprecipitation sequencing (ChIPseq) approach for analysis of fine-needle patient liver biopsy specimens to investigate the role of histone PTMs in chronically HBV-infected patients. Using 18 specimens from patients in different stages of chronic HBV infection, our work shows that the profile of histone PTMs in chronic infection is more nuanced than previously observed in in vitro models of acute infection. In line with our previous findings, we find that the majority of HBV-derived sequences are associated with the activating histone PTM H3K4me3. However, we show a striking interpatient variability of its deposition in this patient cohort correlated with viral transcription and patient HBV early antigen (HBeAg) status. Unexpectedly, we detected deposition of the classical inhibitory histone PTM H3K9me3 on HBV-DNA in around half of the patient biopsy specimens, which could not be linked to reduced levels of viral transcripts. Our results show that current in vitro models are unable to fully recapitulate the complex epigenetic landscape of chronic HBV infection observed in vivo and demonstrate that fine-needle liver biopsy specimens can provide sufficient material to further investigate the interaction of viral and host proteins on HBV-DNA.IMPORTANCE Hepatitis B virus (HBV) is a major global health concern, chronically infecting millions of patients and contributing to a rising burden of liver disease. The viral genome forms the basis for chronic infection and has been shown to be subject to regulation by epigenetic mechanisms, such as posttranslational modification of histone proteins. Here, we confirm and expand on previous results by adapting a high-resolution technique for analysis of histone modifications for use with patient-derived fine-needle liver biopsy specimens. Our work highlights that the situation in vivo is more complex than predicted by current in vitro models, for example, by suggesting a novel, noncanonical role of the histone modification H3K9me3 in the HBV life cycle. Importantly, enabling the use of fine-needle liver biopsy specimens for such high-resolution analyses may facilitate further research into the epigenetic regulation of the HBV genome.

The Safety and Antiviral Activity of BZF961 with or without Ritonavir in Patients Infected with Hepatitis C Virus: A Randomized, Multicenter Trial.

PURPOSE: Infection with hepatitis C virus is the leading cause of infectious disease mortality in the United States. BZF961 is a novel small molecule inhibitor of the hepatitis C virus NS3-4A protease. Here we present the results of a randomized, double-blinded, placebo-controlled, multicentered study in treatment-naïve patients with chronic hepatitis C virus genotype-1 infection. METHODS: Patients were enrolled sequentially in 2 parts and treated for 3days. BZF961 was administered as monotherapy (500mg BID for 3 days) or in combination with the cytochrome P450 3A4 inhibitor ritonavir to boost its exposure (BZF961 10, 20, or 50mg QD or BID). FINDINGS: BZF961 was safe and well tolerated in the patients studied with no serious adverse events. There were no appreciable differences in adverse events among patients who received BZF961, BZF961 with ritonavir, or placebo. There was a significant, clinically meaningful reduction in viral load from baseline in patients treated either with BZF961 500mg every 12hours alone or BZF961 50mg every 12hours in combination with ritonavir. Activity against the hepatitis C virus of the lower-dose regimens was apparent but more modest. There were no relevant changes from baseline viral loads in placebo-treated patients. IMPLICATIONS: Coadministration of ritonavir with BZF961 boosted BZF961 exposure (including Cmin, which is the clinically relevant parameter associated with antiviral activity) in a therapeutic range with less variability compared with BZF961 alone. For strategic reasons, BZF961 is no longer under development.

Mode of Action of the Monobactam LYS228 and Mechanisms Decreasing In Vitro Susceptibility in Escherichia coli and Klebsiella pneumoniae.

The monobactam scaffold is attractive for the development of new agents to treat infections caused by drug-resistant Gram-negative bacteria because it is stable to metallo-ß-lactamases (MBLs). However, the clinically used monobactam aztreonam lacks stability to serine ß-lactamases (SBLs) that are often coexpressed with MBLs. LYS228 is stable to MBLs and most SBLs. LYS228 bound purified Escherichia coli penicillin binding protein 3 (PBP3) similarly to aztreonam (derived acylation rate/equilibrium dissociation constant [k2/Kd ] of 367,504 s-1 M-1 and 409,229 s-1 M-1, respectively) according to stopped-flow fluorimetry. A gel-based assay showed that LYS228 bound mainly to E. coli PBP3, with weaker binding to PBP1a and PBP1b. Exposing E. coli cells to LYS228 caused filamentation consistent with impaired cell division. No single-step mutants were selected from 12 Enterobacteriaceae strains expressing different classes of ß-lactamases at 8× the MIC of LYS228 (frequency, <2.5 × 10-9). At 4× the MIC, mutants were selected from 2 of 12 strains at frequencies of 1.8 × 10-7 and 4.2 × 10-9 LYS228 MICs were ≤2 µg/ml against all mutants. These frequencies compared favorably to those for meropenem and tigecycline. Mutations decreasing LYS228 susceptibility occurred in ramR and cpxA (Klebsiella pneumoniae) and baeS (E. coli and K. pneumoniae). Susceptibility of E. coli ATCC 25922 to LYS228 decreased 256-fold (MIC, 0.125 to 32 µg/ml) after 20 serial passages. Mutants accumulated mutations in ftsI (encoding the target, PBP3), baeR, acrD, envZ, sucB, and rfaI These results support the continued development of LYS228, which is currently undergoing phase II clinical trials for complicated intraabdominal infection and complicated urinary tract infection (registered at ClinicalTrials.gov under identifiers NCT03377426 and NCT03354754).

Clustering of disulfide-rich peptides provides scaffolds for hit discovery by phage display: application to interleukin-23.

BACKGROUND: Disulfide-rich peptides (DRPs) are found throughout nature. They are suitable scaffolds for drug development due to their small cores, whose disulfide bonds impart extraordinary chemical and biological stability. A challenge in developing a DRP therapeutic is to engineer binding to a specific target. This challenge can be overcome by (i) sampling the large sequence space of a given scaffold through a phage display library and by (ii) panning multiple libraries encoding structurally distinct scaffolds. Here, we implement a protocol for defining these diverse scaffolds, based on clustering structurally defined DRPs according to their conformational similarity. RESULTS: We developed and applied a hierarchical clustering protocol based on DRP structural similarity, followed by two post-processing steps, to classify 806 unique DRP structures into 81 clusters. The 20 most populated clusters comprised 85% of all DRPs. Representative scaffolds were selected from each of these clusters; the representatives were structurally distinct from one another, but similar to other DRPs in their respective clusters. To demonstrate the utility of the clusters, phage libraries were constructed for three of the representative scaffolds and panned against interleukin-23. One library produced a peptide that bound to this target with an IC50 of 3.3 µM. CONCLUSIONS: Most DRP clusters contained members that were diverse in sequence, host organism, and interacting proteins, indicating that cluster members were functionally diverse despite having similar structure. Only 20 peptide scaffolds accounted for most of the natural DRP structural diversity, providing suitable starting points for seeding phage display experiments. Through selection of the scaffold surface to vary in phage display, libraries can be designed that present sequence diversity in architecturally distinct, biologically relevant combinations of secondary structures. We supported this hypothesis with a proof-of-concept experiment in which three phage libraries were constructed and panned against the IL-23 target, resulting in a single-digit µM hit and suggesting that a collection of libraries based on the full set of 20 scaffolds increases the potential to identify efficiently peptide binders to a protein target in a drug discovery program.

In Vitro Characterization of Human Cytomegalovirus-Targeting Therapeutic Monoclonal Antibodies LJP538 and LJP539.

Human cytomegalovirus (HCMV) infection is usually benign in healthy individuals but can cause life-threatening disease in those with compromised immune systems. Approved drugs available to treat HCMV disease, including ganciclovir, cidofovir, and foscarnet, have significant toxicities that limit their use in certain patient populations. LJP538 and LJP539 are human monoclonal antibodies that are being evaluated as immunoglobulin therapeutics. The antibodies target glycoproteins gB and the gH/gL/UL128/UL130/UL131a pentameric complex, respectively. Here we present an in vitro characterization of these antibodies. We show that LJP538 and LJP539 are more potent than a marketed immunoglobulin at inhibiting HCMV infection of various cell lines relevant to pathogenesis. We find that LJP538 and LJP539 are active against a panel of clinical isolates in vitro and demonstrate minor-to-moderate synergy in combination. Passage of HCMV in the presence of LJP538 or LJP539 alone resulted in resistance-associated mutations that mapped to the target genes. However, no loss of susceptibility to the combination of antibodies was observed for >400 days in culture. Finally, the binding regions of LJP538 and LJP539 are conserved among clinical isolates. Taken together, these data support the use of LJP538 and LJP539 in combination for clinical trials in HCMV patients.

Global identification and characterization of both O-GlcNAcylation and phosphorylation at the murine synapse.

O-linked N-acetylglucosamine (O-GlcNAc) is a dynamic, reversible monosaccharide modifier of serine and threonine residues on intracellular protein domains. Crosstalk between O-GlcNAcylation and phosphorylation has been hypothesized. Here, we identified over 1750 and 16,500 sites of O-GlcNAcylation and phosphorylation from murine synaptosomes, respectively. In total, 135 (7%) of all O-GlcNAcylation sites were also found to be sites of phosphorylation. Although many proteins were extensively phosphorylated and minimally O-GlcNAcylated, proteins found to be extensively O-GlcNAcylated were almost always phosphorylated to a similar or greater extent, indicating the O-GlcNAcylation system is specifically targeting a subset of the proteome that is also phosphorylated. Both PTMs usually occur on disordered regions of protein structure, within which, the location of O-GlcNAcylation and phosphorylation is virtually random with respect to each other, suggesting that negative crosstalk at the structural level is not a common phenomenon. As a class, protein kinases are found to be more extensively O-GlcNAcylated than proteins in general, indicating the potential for crosstalk of phosphorylation with O-GlcNAcylation via regulation of enzymatic activity.

ModBase, a database of annotated comparative protein structure models, and associated resources.

ModBase (http://salilab.org/modbase) is a database of annotated comparative protein structure models. The models are calculated by ModPipe, an automated modeling pipeline that relies primarily on Modeller for fold assignment, sequence-structure alignment, model building and model assessment (http://salilab.org/modeller/). ModBase currently contains 10,355,444 reliable models for domains in 2,421,920 unique protein sequences. ModBase allows users to update comparative models on demand, and request modeling of additional sequences through an interface to the ModWeb modeling server (http://salilab.org/modweb). ModBase models are available through the ModBase interface as well as the Protein Model Portal (http://www.proteinmodelportal.org/). Recently developed associated resources include the SALIGN server for multiple sequence and structure alignment (http://salilab.org/salign), the ModEval server for predicting the accuracy of protein structure models (http://salilab.org/modeval), the PCSS server for predicting which peptides bind to a given protein (http://salilab.org/pcss) and the FoXS server for calculating and fitting Small Angle X-ray Scattering profiles (http://salilab.org/foxs).

Prediction of protease substrates using sequence and structure features.

MOTIVATION: Granzyme B (GrB) and caspases cleave specific protein substrates to induce apoptosis in virally infected and neoplastic cells. While substrates for both types of proteases have been determined experimentally, there are many more yet to be discovered in humans and other metazoans. Here, we present a bioinformatics method based on support vector machine (SVM) learning that identifies sequence and structural features important for protease recognition of substrate peptides and then uses these features to predict novel substrates. Our approach can act as a convenient hypothesis generator, guiding future experiments by high-confidence identification of peptide-protein partners. RESULTS: The method is benchmarked on the known substrates of both protease types, including our literature-curated GrB substrate set (GrBah). On these benchmark sets, the method outperforms a number of other methods that consider sequence only, predicting at a 0.87 true positive rate (TPR) and a 0.13 false positive rate (FPR) for caspase substrates, and a 0.79 TPR and a 0.21 FPR for GrB substrates. The method is then applied to approximately 25 000 proteins in the human proteome to generate a ranked list of predicted substrates of each protease type. Two of these predictions, AIF-1 and SMN1, were selected for further experimental analysis, and each was validated as a GrB substrate. AVAILABILITY: All predictions for both protease types are publically available at http://salilab.org/peptide. A web server is at the same site that allows a user to train new SVM models to make predictions for any protein that recognizes specific oligopeptide ligands.

Regulatory elements within the prodomain of Falcipain-2, a cysteine protease of the malaria parasite Plasmodium falciparum.

Falcipain-2, a papain family cysteine protease of the malaria parasite Plasmodium falciparum, plays a key role in parasite hydrolysis of hemoglobin and is a potential chemotherapeutic target. As with many proteases, falcipain-2 is synthesized as a zymogen, and the prodomain inhibits activity of the mature enzyme. To investigate the mechanism of regulation of falcipain-2 by its prodomain, we expressed constructs encoding different portions of the prodomain and tested their ability to inhibit recombinant mature falcipain-2. We identified a C-terminal segment (Leu(155)-Asp(243)) of the prodomain, including two motifs (ERFNIN and GNFD) that are conserved in cathepsin L sub-family papain family proteases, as the mediator of prodomain inhibitory activity. Circular dichroism analysis showed that the prodomain including the C-terminal segment, but not constructs lacking this segment, was rich in secondary structure, suggesting that the segment plays a crucial role in protein folding. The falcipain-2 prodomain also efficiently inhibited other papain family proteases, including cathepsin K, cathepsin L, cathepsin B, and cruzain, but it did not inhibit cathepsin C or tested proteases of other classes. A structural model of pro-falcipain-2 was constructed by homology modeling based on crystallographic structures of mature falcipain-2, procathepsin K, procathepsin L, and procaricain, offering insights into the nature of the interaction between the prodomain and mature domain of falcipain-2 as well as into the broad specificity of inhibitory activity of the falcipain-2 prodomain.

MODBASE, a database of annotated comparative protein structure models and associated resources.

MODBASE (http://salilab.org/modbase) is a database of annotated comparative protein structure models. The models are calculated by MODPIPE, an automated modeling pipeline that relies primarily on MODELLER for fold assignment, sequence-structure alignment, model building and model assessment (http:/salilab.org/modeller). MODBASE currently contains 5,152,695 reliable models for domains in 1,593,209 unique protein sequences; only models based on statistically significant alignments and/or models assessed to have the correct fold are included. MODBASE also allows users to calculate comparative models on demand, through an interface to the MODWEB modeling server (http://salilab.org/modweb). Other resources integrated with MODBASE include databases of multiple protein structure alignments (DBAli), structurally defined ligand binding sites (LIGBASE), predicted ligand binding sites (AnnoLyze), structurally defined binary domain interfaces (PIBASE) and annotated single nucleotide polymorphisms and somatic mutations found in human proteins (LS-SNP, LS-Mut). MODBASE models are also available through the Protein Model Portal (http://www.proteinmodelportal.org/).

Global sequencing of proteolytic cleavage sites in apoptosis by specific labeling of protein N termini.

The nearly 600 proteases in the human genome regulate a diversity of biological processes, including programmed cell death. Comprehensive characterization of protease signaling in complex biological samples is limited by available proteomic methods. We have developed a general approach for global identification of proteolytic cleavage sites using an engineered enzyme to selectively biotinylate free protein N termini for positive enrichment of corresponding N-terminal peptides. Using this method to study apoptosis, we have sequenced 333 caspase-like cleavage sites distributed among 292 protein substrates. These sites are generally not predicted by in vitro caspase substrate specificity but can be used to predict other physiological caspase cleavage sites. Structural bioinformatic studies show that caspase cleavage sites often appear in surface-accessible loops and even occasionally in helical regions. Strikingly, we also find that a disproportionate number of caspase substrates physically interact, suggesting that these dimeric proteases target protein complexes and networks to elicit apoptosis.

Host pathogen protein interactions predicted by comparative modeling.

Pathogens have evolved numerous strategies to infect their hosts, while hosts have evolved immune responses and other defenses to these foreign challenges. The vast majority of host-pathogen interactions involve protein-protein recognition, yet our current understanding of these interactions is limited. Here, we present and apply a computational whole-genome protocol that generates testable predictions of host-pathogen protein interactions. The protocol first scans the host and pathogen genomes for proteins with similarity to known protein complexes, then assesses these putative interactions, using structure if available, and, finally, filters the remaining interactions using biological context, such as the stage-specific expression of pathogen proteins and tissue expression of host proteins. The technique was applied to 10 pathogens, including species of Mycobacterium, apicomplexa, and kinetoplastida, responsible for "neglected" human diseases. The method was assessed by (1) comparison to a set of known host-pathogen interactions, (2) comparison to gene expression and essentiality data describing host and pathogen genes involved in infection, and (3) analysis of the functional properties of the human proteins predicted to interact with pathogen proteins, demonstrating an enrichment for functionally relevant host-pathogen interactions. We present several specific predictions that warrant experimental follow-up, including interactions from previously characterized mechanisms, such as cytoadhesion and protease inhibition, as well as suspected interactions in hypothesized networks, such as apoptotic pathways. Our computational method provides a means to mine whole-genome data and is complementary to experimental efforts in elucidating networks of host-pathogen protein interactions.