Evolutionary druggability for low-dimensional fitness landscapes toward new metrics for antimicrobial applications.

The term 'druggability' describes the molecular properties of drugs or targets in pharmacological interventions and is commonly used in work involving drug development for clinical applications. There are no current analogues for this notion that quantify the drug-target interaction with respect to a given target variant's sensitivity across a breadth of drugs in a panel, or a given drug's range of effectiveness across alleles of a target protein. Using data from low-dimensional empirical fitness landscapes composed of 16 ß-lactamase alleles and 7 ß-lactam drugs, we introduce two metrics that capture (i) the average susceptibility of an allelic variant of a drug target to any available drug in a given panel ('variant vulnerability'), and (ii) the average applicability of a drug (or mixture) across allelic variants of a drug target ('drug applicability'). Finally, we (iii) disentangle the quality and magnitude of interactions between loci in the drug target and the seven drug environments in terms of their mutation by mutation by environment (G x G x E) interactions, offering mechanistic insight into the variant variability and drug applicability metrics. Summarizing, we propose that our framework can be applied to other datasets and pathogen-drug systems to understand which pathogen variants in a clinical setting are the most concerning (low variant vulnerability), and which drugs in a panel are most likely to be effective in an infection defined by standing genetic variation in the pathogen drug target (high drug applicability).

The immune-evasive proline-283 substitution in influenza nucleoprotein increases aggregation propensity without altering the native structure.

Nucleoprotein (NP) is a key structural protein of influenza ribonucleoprotein complexes and is central to viral RNA packing and trafficking. NP also determines the sensitivity of influenza to myxovirus resistance protein 1 (MxA), an innate immunity factor that restricts influenza replication. A few critical MxA-resistant mutations have been identified in NP, including the highly conserved proline-283 substitution. This essential proline-283 substitution impairs influenza growth, a fitness defect that becomes particularly prominent at febrile temperature (39°C) when host chaperones are depleted. Here, we biophysically characterize proline-283 NP and serine-283 NP to test whether the fitness defect is caused by the proline-283 substitution introducing folding defects. We show that the proline-283 substitution changes the folding pathway of NP, making NP more aggregation prone during folding, but does not alter the native structure of the protein. These findings suggest that influenza has evolved to hijack host chaperones to promote the folding of otherwise biophysically incompetent viral proteins that enable innate immune system escape.

ER procollagen storage defect without coupled unfolded protein response drives precocious arthritis.

Collagenopathies are a group of clinically diverse disorders caused by defects in collagen folding and secretion. For example, mutations in the gene encoding collagen type-II, the primary collagen in cartilage, can lead to diverse chondrodysplasias. One example is the Gly1170Ser substitution in procollagen-II, which causes precocious osteoarthritis. Here, we biochemically and mechanistically characterize an induced pluripotent stem cell-based cartilage model of this disease, including both hetero- and homozygous genotypes. We show that Gly1170Ser procollagen-II is notably slow to fold and secrete. Instead, procollagen-II accumulates intracellularly, consistent with an endoplasmic reticulum (ER) storage disorder. Owing to unique features of the collagen triple helix, this accumulation is not recognized by the unfolded protein response. Gly1170Ser procollagen-II interacts to a greater extent than wild-type with specific proteostasis network components, consistent with its slow folding. These findings provide mechanistic elucidation into the etiology of this disease. Moreover, the cartilage model will enable rapid testing of therapeutic strategies to restore proteostasis in the collagenopathies.

Epistasis and pleiotropy shape biophysical protein subspaces associated with drug resistance.

Protein space is a rich analogy for genotype-phenotype maps, where amino acid sequence is organized into a high-dimensional space that highlights the connectivity between protein variants. It is a useful abstraction for understanding the process of evolution, and for efforts to engineer proteins towards desirable phenotypes. Few mentions of protein space consider how protein phenotypes can be described in terms of their biophysical components, nor do they rigorously interrogate how forces like epistasis-describing the nonlinear interaction between mutations and their phenotypic consequences-manifest across these components. In this study, we deconstruct a low-dimensional protein space of a bacterial enzyme (dihydrofolate reductase; DHFR) into "subspaces" corresponding to a set of kinetic and thermodynamic traits [k_{cat}, K_{M}, K_{i}, and T_{m} (melting temperature)]. We then examine how combinations of three mutations (eight alleles in total) display pleiotropy, or unique effects on individual subspace traits. We examine protein spaces across three orthologous DHFR enzymes (Escherichia coli, Listeria grayi, and Chlamydia muridarum), adding a genotypic context dimension through which epistasis occurs across subspaces. In doing so, we reveal that protein space is a deceptively complex notion, and that future applications to bioengineering should consider how interactions between amino acid substitutions manifest across different phenotypic subspaces.

The Immune-Evasive Proline 283 Substitution in Influenza Nucleoprotein Increases Aggregation Propensity Without Altering the Native Structure.

Nucleoprotein (NP) is a key structural protein of influenza ribonucleoprotein complexes and is central to viral RNA packing and trafficking. In human cells, the interferon induced Myxovirus resistance protein 1 (MxA) binds to NP and restricts influenza replication. This selection pressure has caused NP to evolve a few critical MxA-resistant mutations, particularly the highly conserved Pro283 substitution. Previous work showed that this essential Pro283 substitution impairs influenza growth, and the fitness defect becomes particularly prominent at febrile temperature (39 °C) when host chaperones are depleted. Here, we biophysically characterize Pro283 NP and Ser283 NP to test if the fitness defect is owing to Pro283 substitution introducing folding defects. We show that the Pro283 substitution changes the folding pathway of NP without altering the native structure, making NP more aggregation prone during folding. These findings suggest that influenza has evolved to hijack host chaperones to promote the folding of otherwise biophysically incompetent viral proteins that enable innate immune system escape. Teaser: Pro283 substitution in flu nucleoprotein introduces folding defects, and makes influenza uniquely dependent on host chaperones.

Epistasis meets pleiotropy in shaping biophysical protein subspaces associated with antimicrobial resistance.

Protein space is a rich analogy for genotype-phenotype maps, where amino acid sequence is organized into a high-dimensional space that highlights the connectivity between protein variants. It is a useful abstraction for understanding the process of evolution, and for efforts to engineer proteins towards desirable phenotypes. Few framings of protein space consider how higher-level protein phenotypes can be described in terms of their biophysical dimensions, nor do they rigorously interrogate how forces like epistasis-describing the nonlinear interaction between mutations and their phenotypic consequences-manifest across these dimensions. In this study, we deconstruct a low-dimensional protein space of a bacterial enzyme (dihydrofolate reductase; DHFR) into "subspaces" corresponding to a set of kinetic and thermodynamic traits [(kcat, KM, Ki, and Tm (melting temperature)]. We then examine how three mutations (eight alleles in total) display pleiotropy in their interactions across these subspaces. We extend this approach to examine protein spaces across three orthologous DHFR enzymes (Escherichia coli, Listeria grayi, and Chlamydia muridarum), adding a genotypic context dimension through which epistasis occurs across subspaces. In doing so, we reveal that protein space is a deceptively complex notion, and that the process of protein evolution and engineering should consider how interactions between amino acid substitutions manifest across different phenotypic subspaces.

Evolutionary druggability: leveraging low-dimensional fitness landscapes towards new metrics for antimicrobial applications.

The term "druggability" describes the molecular properties of drugs or targets in pharmacological interventions and is commonly used in work involving drug development for clinical applications. There are no current analogues for this notion that quantify the drug-target interaction with respect to a given target variant's sensitivity across a breadth of drugs in a panel, or a given drug's range of effectiveness across alleles of a target protein. Using data from low-dimensional empirical fitness landscapes composed of 16 ß-lactamase alleles and seven ß-lactam drugs, we introduce two metrics that capture (i) the average susceptibility of an allelic variant of a drug target to any available drug in a given panel ("variant vulnerability"), and (ii) the average applicability of a drug (or mixture) across allelic variants of a drug target ("drug applicability"). Finally, we (iii) disentangle the quality and magnitude of interactions between loci in the drug target and the seven drug environments in terms of their mutation by mutation by environment (G × G × E) interactions, offering mechanistic insight into the variant variability and drug applicability metrics. Summarizing, we propose that our framework can be applied to other datasets and pathogen-drug systems to understand which pathogen variants in a clinical setting are the most concerning (low variant vulnerability), and which drugs in a panel are most likely to be effective in an infection defined by standing genetic variation in the pathogen drug target (high drug applicability).

Viral Evolution Shaped by Host Proteostasis Networks.

Understanding the factors that shape viral evolution is critical for developing effective antiviral strategies, accurately predicting viral evolution, and preventing pandemics. One fundamental determinant of viral evolution is the interplay between viral protein biophysics and the host machineries that regulate protein folding and quality control. Most adaptive mutations in viruses are biophysically deleterious, resulting in a viral protein product with folding defects. In cells, protein folding is assisted by a dynamic system of chaperones and quality control processes known as the proteostasis network. Host proteostasis networks can determine the fates of viral proteins with biophysical defects, either by assisting with folding or by targeting them for degradation. In this review, we discuss and analyze new discoveries revealing that host proteostasis factors can profoundly shape the sequence space accessible to evolving viral proteins. We also discuss the many opportunities for research progress proffered by the proteostasis perspective on viral evolution and adaptation.

Using CRISPR/Cas9 to generate a heterozygous COL2A1 p.R719C iPSC line (MCRIi019-A-6) model of human precocious osteoarthritis.

The human iPSC line MCRIi019-A-6 was generated using CRISPR/Cas9-mediated gene editing to introduce a heterozygous COL2A1 exon 33 c.2155C>T (p.R719C) mutation into the control human iPSC line MCRIi019-A. Both the edited and parental lines display typical iPSC characteristics, including the expression of pluripotency markers, the ability to be differentiated into the three germ lines, and a normal karyotype. This cell line, along with the isogenic control line, can be used to study the molecular pathology of precocious osteoarthritis in a human model, more broadly understand type II collagenopathies, and explore novel therapeutic targets for this class of diseases.

Expanded MutaT7 toolkit efficiently and simultaneously accesses all possible transition mutations in bacteria.

Targeted mutagenesis mediated by nucleotide base deaminase-T7 RNA polymerase fusions has recently emerged as a novel and broadly useful strategy to power genetic diversification in the context of in vivo directed evolution campaigns. Here, we expand the utility of this approach by introducing a highly active adenosine deaminase-T7 RNA polymerase fusion protein (eMutaT7A→G), resulting in higher mutation frequencies to enable more rapid directed evolution. We also assess the benefits and potential downsides of using this more active mutator. We go on to show in Escherichia coli that adenosine deaminase-bearing mutators (MutaT7A→G or eMutaT7A→G) can be employed in tandem with a cytidine deaminase-bearing mutator (MutaT7C→T) to introduce all possible transition mutations simultaneously. We illustrate the efficacy of this in vivo mutagenesis approach by exploring mutational routes to antibacterial drug resistance. This work sets the stage for general application of optimized MutaT7 tools able to induce all types of transition mutations during in vivo directed evolution campaigns across diverse organisms.

Loss of heat shock factor 1 promotes hepatic stellate cell activation and drives liver fibrosis.

Liver fibrosis is an aberrant wound healing response that results from chronic injury and is mediated by hepatocellular death and activation of hepatic stellate cells (HSCs). While induction of oxidative stress is well established in fibrotic livers, there is limited information on stress-mediated mechanisms of HSC activation. Cellular stress triggers an adaptive defense mechanism via master protein homeostasis regulator, heat shock factor 1 (HSF1), which induces heat shock proteins to respond to proteotoxic stress. Although the importance of HSF1 in restoring cellular homeostasis is well-established, its potential role in liver fibrosis is unknown. Here, we show that HSF1 messenger RNA is induced in human cirrhotic and murine fibrotic livers. Hepatocytes exhibit nuclear HSF1, whereas stellate cells expressing alpha smooth muscle actin do not express nuclear HSF1 in human cirrhosis. Interestingly, despite nuclear HSF1, murine fibrotic livers did not show induction of HSF1 DNA binding activity compared with controls. HSF1-deficient mice exhibit augmented HSC activation and fibrosis despite limited pro-inflammatory cytokine response and display delayed fibrosis resolution. Stellate cell and hepatocyte-specific HSF1 knockout mice exhibit higher induction of profibrogenic response, suggesting an important role for HSF1 in HSC activation and fibrosis. Stable expression of dominant negative HSF1 promotes fibrogenic activation of HSCs. Overactivation of HSF1 decreased phosphorylation of JNK and prevented HSC activation, supporting a protective role for HSF1. Our findings identify an unconventional role for HSF1 in liver fibrosis. Conclusion: Our results show that deficiency of HSF1 is associated with exacerbated HSC activation promoting liver fibrosis, whereas activation of HSF1 prevents profibrogenic HSC activation.

Thermal Proteome Profiling Reveals the O-GlcNAc-Dependent Meltome.

Posttranslational modifications alter the biophysical properties of proteins and thereby influence cellular physiology. One emerging manner by which such modifications regulate protein functions is through their ability to perturb protein stability. Despite the increasing interest in this phenomenon, there are few methods that enable global interrogation of the biophysical effects of posttranslational modifications on the proteome. Here, we describe an unbiased proteome-wide approach to explore the influence of protein modifications on the thermodynamic stability of thousands of proteins in parallel. We apply this profiling strategy to study the effects of O-linked N-acetylglucosamine (O-GlcNAc), an abundant modification found on hundreds of proteins in mammals that has been shown in select cases to stabilize proteins. Using this thermal proteomic profiling strategy, we identify a set of 72 proteins displaying O-GlcNAc-dependent thermostability and validate this approach using orthogonal methods targeting specific proteins. These collective observations reveal that the majority of proteins influenced by O-GlcNAc are, surprisingly, destabilized by O-GlcNAc and cluster into distinct macromolecular complexes. These results establish O-GlcNAc as a bidirectional regulator of protein stability and provide a blueprint for exploring the impact of any protein modification on the meltome of, in principle, any organism.

Collagen misfolding mutations: the contribution of the unfolded protein response to the molecular pathology.

Mutations in collagen genes cause a broad range of connective tissue pathologies. Structural mutations that impact procollagen assembly or triple helix formation and stability are a common and important mutation class. How misfolded procollagens engage with the cellular proteostasis machinery and whether they can elicit a cytotoxic unfolded protein response (UPR) is a topic of considerable research interest. Such interest is well justified since modulating the UPR could offer a new approach to treat collagenopathies for which there are no current disease mechanism-targeting therapies. This review scrutinizes the evidence underpinning the view that endoplasmic reticulum stress and chronic UPR activation contributes significantly to the pathophysiology of the collagenopathies. While there is strong evidence that the UPR contributes to the pathology for collagen X misfolding mutations, the evidence that misfolding mutations in other collagen types induce a canonical, cytotoxic UPR is incomplete. To gain a more comprehensive understanding about how the UPR amplifies to pathology, and thus what types of manipulations of the UPR might have therapeutic relevance, much more information is needed about how specific misfolding mutation types engage differentially with the UPR and downstream signaling responses. Most importantly, since the capacity of the proteostasis machinery to respond to collagen misfolding is likely to vary between cell types, reflecting their functional roles in collagen and extracellular matrix biosynthesis, detailed studies on the UPR should focus as much as possible on the actual target cells involved in the collagen pathologies.

The endoplasmic reticulum proteostasis network profoundly shapes the protein sequence space accessible to HIV envelope.

The sequence space accessible to evolving proteins can be enhanced by cellular chaperones that assist biophysically defective clients in navigating complex folding landscapes. It is also possible, at least in theory, for proteostasis mechanisms that promote strict quality control to greatly constrain accessible protein sequence space. Unfortunately, most efforts to understand how proteostasis mechanisms influence evolution rely on artificial inhibition or genetic knockdown of specific chaperones. The few experiments that perturb quality control pathways also generally modulate the levels of only individual quality control factors. Here, we use chemical genetic strategies to tune proteostasis networks via natural stress response pathways that regulate the levels of entire suites of chaperones and quality control mechanisms. Specifically, we upregulate the unfolded protein response (UPR) to test the hypothesis that the host endoplasmic reticulum (ER) proteostasis network shapes the sequence space accessible to human immunodeficiency virus-1 (HIV-1) envelope (Env) protein. Elucidating factors that enhance or constrain Env sequence space is critical because Env evolves extremely rapidly, yielding HIV strains with antibody- and drug-escape mutations. We find that UPR-mediated upregulation of ER proteostasis factors, particularly those controlled by the IRE1-XBP1s UPR arm, globally reduces Env mutational tolerance. Conserved, functionally important Env regions exhibit the largest decreases in mutational tolerance upon XBP1s induction. Our data indicate that this phenomenon likely reflects strict quality control endowed by XBP1s-mediated remodeling of the ER proteostasis environment. Intriguingly, and in contrast, specific regions of Env, including regions targeted by broadly neutralizing antibodies, display enhanced mutational tolerance when XBP1s is induced, hinting at a role for host proteostasis network hijacking in potentiating antibody escape. These observations reveal a key function for proteostasis networks in decreasing instead of expanding the sequence space accessible to client proteins, while also demonstrating that the host ER proteostasis network profoundly shapes the mutational tolerance of Env in ways that could have important consequences for HIV adaptation.

Directed evolution in mammalian cells.

Directed evolution experiments are typically carried out using in vitro systems, bacteria, or yeast-even when the goal is to probe or modulate mammalian biology. Performing directed evolution in systems that do not match the intended mammalian environment severely constrains the scope and functionality of the targets that can be evolved. We review new platforms that are now making it possible to use the mammalian cell itself as the setting for directed evolution and present an overview of frontier challenges and high-impact targets for this approach.

Collagen's enigmatic, highly conserved N-glycan has an essential proteostatic function.

Intracellular procollagen folding begins at the protein's C-terminal propeptide (C-Pro) domain, which initiates triple-helix assembly and defines the composition and chain register of fibrillar collagen trimers. The C-Pro domain is later proteolytically cleaved and excreted from the body, while the mature triple helix is incorporated into the extracellular matrix. The procollagen C-Pro domain possesses a single N-glycosylation site that is widely conserved in all the fibrillar procollagens across humans and diverse other species. Given that the C-Pro domain is removed once procollagen folding is complete, the N-glycan might be presumed to be important for folding. Surprisingly, however, there is no difference in the folding and secretion of N-glycosylated versus non-N-glycosylated collagen type-I, leaving the function of the N-glycan unclear. We hypothesized that the collagen N-glycan might have a context-dependent function, specifically, that it could be required to promote procollagen folding only when proteostasis is challenged. We show that removal of the N-glycan from misfolding-prone C-Pro domain variants does indeed cause serious procollagen and ER proteostasis defects. The N-glycan promotes folding and secretion of destabilized C-Pro variants by providing access to the ER's lectin-based chaperone machinery. Finally, we show that the C-Pro N-glycan is actually critical for the folding and secretion of even wild-type procollagen under ER stress conditions. Such stress is commonly incurred during development, wound healing, and other processes in which collagen production plays a key role. Collectively, these results establish an essential, context-dependent function for procollagen's previously enigmatic N-glycan, wherein the carbohydrate moiety buffers procollagen folding against proteostatic challenge.

CRISPR/Cas9 editing to generate a heterozygous COL2A1 p.G1170S human chondrodysplasia iPSC line, MCRIi019-A-2, in a control iPSC line, MCRIi019-A.

To develop an in vitro disease model of a human chondrodysplasia, we used CRISPR/Cas9 gene editing to generate a heterozygous COL2A1 exon 50 c.3508 GGT > TCA (p.G1170S) mutation in a control human iPSC line. Both the control and COL2A1 mutant lines displayed typical iPSC characteristics, including normal cell morphology, expression of pluripotency markers, the ability to differentiate into endoderm, ectoderm and mesoderm lineages and normal karyotype. These chondrodysplasia mutant and isogenic control cell lines can be used to explore disease mechanisms underlying type II collagenopathies and aid in the discovery of new therapeutic strategies.

Elucidation of proteostasis defects caused by osteogenesis imperfecta mutations in the collagen-α2(I) C-propeptide domain.

Intracellular collagen assembly begins with the oxidative folding of ∼30-kDa C-terminal propeptide (C-Pro) domains. Folded C-Pro domains then template the formation of triple helices between appropriate partner strands. Numerous C-Pro missense variants that disrupt or delay triple-helix formation are known to cause disease, but our understanding of the specific proteostasis defects introduced by these variants remains immature. Moreover, it is unclear whether or not recognition and quality control of misfolded C-Pro domains is mediated by recognizing stalled assembly of triple-helical domains or by direct engagement of the C-Pro itself. Here, we integrate biochemical and cellular approaches to illuminate the proteostasis defects associated with osteogenesis imperfecta-causing mutations within the collagen-α2(I) C-Pro domain. We first show that "C-Pro-only" constructs recapitulate key aspects of the behavior of full-length Colα2(I) constructs. Of the variants studied, perhaps the most severe assembly defects are associated with C1163R C-Proα2(I), which is incapable of forming stable trimers and is retained within cells. We find that the presence or absence of an unassembled triple-helical domain is not the key feature driving cellular retention versus secretion. Rather, the proteostasis network directly engages the misfolded C-Pro domain itself to prevent secretion and initiate clearance. Using MS-based proteomics, we elucidate how the endoplasmic reticulum (ER) proteostasis network differentially engages misfolded C1163R C-Proα2(I) and targets it for ER-associated degradation. These results provide insights into collagen folding and quality control with the potential to inform the design of proteostasis network-targeted strategies for managing collagenopathies.

HSF1 Activation Can Restrict HIV Replication.

Host protein folding stress responses can play important roles in RNA virus replication and evolution. Prior work suggested a complicated interplay between the cytosolic proteostasis stress response, controlled by the transcriptional master regulator heat shock factor 1 (HSF1), and human immunodeficiency virus-1 (HIV-1). We sought to uncouple HSF1 transcription factor activity from cytotoxic proteostasis stress and thereby better elucidate the proposed role(s) of HSF1 in the HIV-1 lifecycle. To achieve this objective, we used chemical genetic, stress-independent control of HSF1 activity to establish whether and how HSF1 influences HIV-1 replication. Stress-independent HSF1 induction decreased both the total quantity and infectivity of HIV-1 virions. Moreover, HIV-1 was unable to escape HSF1-mediated restriction over the course of several serial passages. These results clarify the interplay between the host's heat shock response and HIV-1 infection and motivate continued investigation of chaperones as potential antiviral therapeutic targets.