Proteins Deleterious on Overexpression Are Associated with High Intrinsic Disorder, Specific Interaction Domains, and Low Abundance

Novel insights into construct toxicity, strain optimization, and primary sequence design for producing recombinant silk fibroin and elastin-like peptide in E. coli

Spider silk proteins (spidroins) are a remarkable class of biomaterials that exhibit a unique combination of high-value attributes and can be processed into numerous morphologies for targeted applications in diverse fields. Recombinant production of spidroins represents the most promising route towards establishing the industrial production of the material, however, recombinant spider silk production suffers from fundamental difficulties that includes low titers, plasmid instability, and translational inefficiencies. In this work, we sought to gain a deeper understanding of upstream bottlenecks that exist in the field through the production of a panel of systematically varied spidroin sequences in multiple E. coli strains. A restriction on basal expression and specific genetic mutations related to stress responses were identified as primary factors that facilitated higher titers of the recombinant silk constructs. Using these findings, a novel strain of E. coli was created that produces recombinant silk constructs at levels 4-33 times higher than standard BL21(DE3). However, these findings did not extend to a similar recombinant protein, an elastin-like peptide. It was found that the recombinant silk proteins, but not the elastin-like peptide, exert toxicity on the E. coli host system, possibly through their high degree of intrinsic disorder. Along with strain engineering, a bioprocess design that utilizes longer culturing times and attenuated induction was found to raise recombinant silk titers by seven-fold and mitigate toxicity. Targeted alteration to the primary sequence of the recombinant silk constructs was also found to mitigate toxicity. These findings identify multiple points of focus for future work seeking to further optimize the recombinant production of silk proteins and is the first work to identify the intrinsic disorder and subsequent toxicity of certain spidroin constructs as a primary factor related to the difficulties of production.

Natural variation in the consequences of gene overexpression and its implications for evolutionary trajectories

Copy number variation through gene or chromosome amplification provides a route for rapid phenotypic variation and supports the long-term evolution of gene functions. Although the evolutionary importance of copy-number variation is known, little is understood about how genetic background influences its tolerance. Here, we measured fitness costs of over 4000 overexpressed genes in 15 Saccharomyces cerevisiae strains representing different lineages, to explore natural variation in tolerating gene overexpression (OE). Strain-specific effects dominated the fitness costs of gene OE. We report global differences in the consequences of gene OE, independent of the amplified gene, as well as gene-specific effects that were dependent on the genetic background. Natural variation in the response to gene OE could be explained by several models, including strain-specific physiological differences, resource limitations, and regulatory sensitivities. This work provides new insight on how genetic background influences tolerance to gene amplification and the evolutionary trajectories accessible to different backgrounds.

Gene-wise variant burden and genomic characterization of nearly every gene

Aim: Current gene-level prioritization methods aim to provide information for further prioritization of ‘disease-causing’ mutations. Since, they are inherently biased toward disease genes, methods specific to pharmacogenetic (PGx) genes are required. Methods: We proposed a gene-wise variant burden (GVB) method that integrates in silico deleteriousness scores of the multitude of variants of a given gene at a personal-genome level. Results: GVB in its simplest form outperformed the two state-of-the-art methods with regard to predicting pharmacogenes and complex disease genes but not for rare Mendelian disease genes. GVB* adjusted by paralog counts robustly performed well in most of the pharmacogenetic subcategories. Seven molecular genetic features well characterized the unique genomic properties of PGx, complex, and Mendelian disease genes. Conclusion: Altogether, GVB is an individual-specific genescore, especially advantageous for PGx studies.

Faster Evolving Primate Genes Are More Likely to Duplicate

An attractive and long-standing hypothesis regarding the evolution of genes after duplication posits that the duplication event creates new evolutionary possibilities by releasing a copy of the gene from constraint. Apparent support was found in numerous analyses, particularly, the observation of higher rates of evolution in duplicated as compared with singleton genes. Could it, instead, be that more duplicable genes (owing to mutation, fixation, or retention biases) are intrinsically faster evolving? To uncouple the measurement of rates of evolution from the determination of duplicate or singleton status, we measure the rates of evolution in singleton genes in outgroup primate lineages but classify these genes as to whether they have duplicated or not in a crown group of great apes. We find that rates of evolution are higher in duplicable genes prior to the duplication event. In part this is owing to a negative correlation between coding sequence length and rate of evolution, coupled with a bias toward smaller genes being more duplicable. The effect is masked by difference in expression rate between duplicable genes and singletons. Additionally, in contradiction to the classical assumption, we find no convincing evidence for an increase in dN/dS after duplication, nor for rate asymmetry between duplicates. We conclude that high rates of evolution of duplicated genes are not solely a consequence of the duplication event, but are rather a predictor of duplicability. These results are consistent with a model in which successful gene duplication events in mammals are skewed toward events of minimal phenotypic impact.

Estimating the protein burden limit of yeast cells by measuring the expression limits of glycolytic proteins

The ultimate overexpression of a protein could cause growth defects, which are known as the protein burden. However, the expression limit at which the protein-burden effect is triggered is still unclear. To estimate this limit, we systematically measured the overexpression limits of glycolytic proteins in Saccharomyces cerevisiae. The limits of some glycolytic proteins were up to 15% of the total cellular protein. These limits were independent of the proteins’ catalytic activities, a finding that was supported by an in silico analysis. Some proteins had low expression limits that were explained by their localization and metabolic perturbations. The codon usage should be highly optimized to trigger the protein-burden effect, even under strong transcriptional induction. The S–S-bond-connected aggregation mediated by the cysteine residues of a protein might affect its expression limit. Theoretically, only non-harmful proteins could be expressed up to the protein-burden limit. Therefore, we established a framework to distinguish proteins that are harmful and non-harmful upon overexpression.

If a cell makes too much of a given protein, it can sometimes cause problems and impair the cell’s growth. Overproducing some proteins may deplete the cell’s limited resources, meaning it does not have enough to make other more essential proteins. This phenomenon is known as the protein burden effect. Theoretically, only harmless proteins can be overproduced up to a level where growth would be impaired in this way. Conversely, if an overproduced protein causes harm before it becomes a burden on resources, scientists must consider other mechanisms to explain the cell’s problems, namely that the protein itself is harmful.

Knowing the ultimate level of protein production that could cause the protein burden effect – the protein burden limit – would allow scientists to distinguish between harmful and non-harmful proteins. However, to date, this limit had not been defined for any cell.

Eguchi et al. have now tried to estimate the protein burden limit for budding yeast – one of the best-studied experimental organisms. The experiments first focused on enzymes involved in alcoholic fermentation because they were expected to be non-harmful. Some of these enzymes were overproduced to the level were the made up 15% of all the cell’s proteins before they started to cause growth defects. The same results were seen with versions of the enzymes that had been mutated to be less active, leading Eguchi et al. to conclude that this level is the protein burden limit.

In other experiments, harmful enzymes could only be overproduced to levels that were far less than this proposed protein burden limit. These enzymes caused problems for the yeast in several ways, including interfering with biochemical reactions and forming large aggregates in the cell. Lastly, Eguchi et al. looked at the yeast’s genetic code and saw that most of its genes seemed to have evolved to specifically limit the production of proteins to a level that would avoid the unwanted protein burden effect.

Together these findings establish a framework to clearly distinguish between harmful and non-harmful proteins. This framework will be useful to understand the different reasons why the overproduction of certain proteins, which is seen in neurodegenerative diseases and cancer cells, can cause problems for cells.

Intrinsic protein disorder reduces small-scale gene duplicability

Whereas the rate of gene duplication is relatively high, only certain duplications survive the filter of natural selection and can contribute to genome evolution. However, the reasons why certain genes can be retained after duplication whereas others cannot remain largely unknown. Many proteins contain intrinsically disordered regions (IDRs), whose structures fluctuate between alternative conformational states. Due to their high flexibility, IDRs often enable protein–protein interactions and are the target of post-translational modifications. Intrinsically disordered proteins (IDPs) have characteristics that might either stimulate or hamper the retention of their encoding genes after duplication. On the one hand, IDRs may enable functional diversification, thus promoting duplicate retention. On the other hand, increased IDP availability is expected to result in deleterious unspecific interactions. Here, we interrogate the proteomes of human, Drosophila melanogaster, Caenorhabditis elegans, Saccharomyces cerevisiae, Arabidopsis thaliana and Escherichia coli, in order to ascertain the impact of protein intrinsic disorder on gene duplicability. We show that, in general, proteins encoded by duplicated genes tend to be less disordered than those encoded by singletons. The only exception is proteins encoded by ohnologs, which tend to be more disordered than those encoded by singletons or genes resulting from small-scale duplications. Our results indicate that duplication of genes encoding IDPs outside the context of whole-genome duplication (WGD) is often deleterious, but that IDRs facilitate retention of duplicates in the context of WGD. We discuss the potential evolutionary implications of our results.

Transient protein-protein interactions perturb E. coli metabolome and cause gene dosage toxicity

Gene dosage toxicity (GDT) is an important factor that determines optimal levels of protein abundances, yet its molecular underpinnings remain unknown. Here, we demonstrate that overexpression of DHFR in E. coli causes a toxic metabolic imbalance triggered by interactions with several functionally related enzymes. Though deleterious in the overexpression regime, surprisingly, these interactions are beneficial at physiological concentrations, implying their functional significance in vivo. Moreover, we found that overexpression of orthologous DHFR proteins had minimal effect on all levels of cellular organization - molecular, systems, and phenotypic, in sharp contrast to E. coli DHFR. Dramatic difference of GDT between 'E. coli's self' and 'foreign' proteins suggests the crucial role of evolutionary selection in shaping protein-protein interaction (PPI) networks at the whole proteome level. This study shows how protein overexpression perturbs a dynamic metabolon of weak yet potentially functional PPI, with consequences for the metabolic state of cells and their fitness.

Quantitative nature of overexpression experiments

Overexpression experiments are sometimes considered as qualitative experiments designed to identify novel proteins and study their function. However, in order to draw conclusions regarding protein overexpression through association analyses using large-scale biological data sets, we need to recognize the quantitative nature of overexpression experiments. Here I discuss the quantitative features of two different types of overexpression experiment: absolute and relative. I also introduce the four primary mechanisms involved in growth defects caused by protein overexpression: resource overload, stoichiometric imbalance, promiscuous interactions, and pathway modulation associated with the degree of overexpression.

Evaluating the fitness cost of protein expression in Saccharomyces cerevisiae

Protein metabolism is one of the most costly processes in the cell and is therefore expected to be under the effective control of natural selection. We stimulated yeast strains to overexpress each single gene product to approximately 1% of the total protein content. Consistent with previous reports, we found that excessive expression of proteins containing disordered or membrane-protruding regions resulted in an especially high fitness cost. We estimated these costs to be nearly twice as high as for other proteins. There was a ten-fold difference in cost if, instead of entire proteins, only the disordered or membrane-embedded regions were compared with other segments. Although the cost of processing bulk protein was measurable, it could not be explained by several tested protein features, including those linked to translational efficiency or intensity of physical interactions after maturation. It most likely included a number of individually indiscernible effects arising during protein synthesis, maturation, maintenance, (mal)functioning, and disposal. When scaled to the levels normally achieved by proteins in the cell, the fitness cost of dealing with one amino acid in a standard protein appears to be generally very low. Many single amino acid additions or deletions are likely to be neutral even if the effective population size is as large as that of the budding yeast. This should also apply to substitutions. Selection is much more likely to operate if point mutations affect protein structure by, for example, extending or creating stretches that tend to unfold or interact improperly with membranes.

Identification of virulence genes in the crucifer anthracnose fungus Colletotrichum higginsianum by insertional mutagenesis

To investigate the molecular and genetic mechanisms underlying virulence of Colletotrichum higginsianum on Arabidopsis thaliana, a T-DNA insertion mutant library of C. higginsianum, the causal agent of crucifer anthracnose, was established using Agrobacterium tumefaciens-mediated transformation. Among 875 transformants tested for virulence on Arabidopsis, six mutants with altered virulence, including an appressorial melanin-deficient mutant T734, two mutants defective in penetration, T45 and B30, and three mutants, T679, T732 and T801, that cause hypersensitive reactions on host Arabidopsis, were obtained. Southern blot analysis indicated that the mutants T732 and T734 harbored single-site T-DNA integrations, while B30 harbored two T-DNA insertions. Border flanking sequences of T-DNAs from these mutants were recovered by inverse polymerase chain reaction (PCR) and thermal asymmetric interlaced PCR. Sequence analyses revealed that single T-DNA insertions in mutant T734 targeted the coding region of a gene with unknown function, and in mutant T732 targeted a gene encoding a copper amine oxidase. The two T-DNA insertion sites in mutant B30 were found in the coding region of a gene encoding an exosome component and in the upstream region of a DUF221-domain gene. None of these genes have previously been implicated in virulence of the phytopathogenic fungi. Among these avirulent mutants, T734 showed altered color in colony growth and produced melanin-deficient, albino appressoria. The T-DNA insert in T734 was detected in the coding region of a gene named C. higginsianum melanin-deficiency gene (Ch-MEL1), which is highly similar to a gene encoding a hypothetical protein in Colletotrichum gloeosporioides (GenBank ELA33048). To validate whether the Ch-MEL1 gene was associated with virulence of the mutant T734, a targeted gene disruption and complementation approach was used. The appressoria of ▵Ch-mel1 null mutants were defective in melanization and failed to penetrate the host epidermal cells. When inoculated onto the wounded leaf tissues, the ▵Ch-mel1 mutants grew on host tissues but failed to cause lesions beyond the wound site. In contrast, both the complement C▵Ch-mel1-2 and the wild type produced melanized appressoria and caused necrosis on leaves of Arabidopsis. Ch-MEL1 is required for both appressorial melanin production in C. higginsianum and post-invasive growth in host tissues. Together with identification of other avirulent mutants and their associated genes, this study provides novel insights into molecular mechanisms underlying virulence of the hemibiotroph, C. higginsianum.

Transcription initiation arising from E-cadherin/CDH1 intron2: a novel protein isoform that increases gastric cancer cell invasion and angiogenesis

Disruption of E-cadherin (CDH1 gene) expression, subcellular localization or function arises during initiation and progression of almost 90% of all epithelial carcinomas. Nevertheless, the mechanisms through which this occurs are largely unknown. Previous studies showed the importance of CDH1 intron 2 sequences for proper gene and protein expression, supporting these as E-cadherin cis-modulators. Through RACE and RT-PCR, we searched for transcription events arising from CDH1 intron 2 and discovered several new transcripts. One, named CDH1a, with high expression in spleen and absent from normal stomach, was demonstrated to be translated into a novel isoform, differing from canonical E-cadherin in its N-terminal, as determined by mass spectrometry. Quantitative and functional assays showed that when overexpressed in an E-cadherin negative context, CDH1a replaced canonical protein interactions and functions. However, when co-expressed with canonical E-cadherin, CDH1a increased cell invasion and angiogenesis. Further, interferon-induced gene IFITM1 and IFI27 levels were increased upon CDH1a overexpression. Effects on invasion and IFITM1 and IFI27 expression were reverted upon CDH1a-specific knockdown. Importantly, CDH1a was de novo expressed in gastric cancer cell lines. This study presents a new mechanism by which E-cadherin functions are impaired by cis-regulatory mechanisms possibly with the involvement of inflammatory machinery. If confirmed in other cancer models, our data enclose potential for designing targeted therapies to rescue E-cadherin function.

Intrinsically disordered proteins: regulation and disease

Intrinsically disordered proteins (IDPs) are enriched in signaling and regulatory functions because disordered segments permit interaction with several proteins and hence the re-use of the same protein in multiple pathways. Understanding IDP regulation is important because altered expression of IDPs is associated with many diseases. Recent studies show that IDPs are tightly regulated and that dosage-sensitive genes encode proteins with disordered segments. The tight regulation of IDPs may contribute to signaling fidelity by ensuring that IDPs are available in appropriate amounts and not present longer than needed. The altered availability of IDPs may result in sequestration of proteins through non-functional interactions involving disordered segments (i.e., molecular titration), thereby causing an imbalance in signaling pathways. We discuss the regulation of IDPs, address implications for signaling, disease and drug development, and outline directions for future research.