Comprehensive bioinformatics analysis of cell-free protein synthesis: identification of multiple protein properties that correlate with successful expression

Protein pI and Intracellular Localization

The protein isoelectric point (pI) can be calculated from an amino acid sequence using computational analysis in a good agreement with experimental data. Availability of whole-genome sequences empowers comparative studies of proteome-wide pI distributions. It was found that the whole-proteome distributions of protein pI values are multimodal in different species. It was further hypothesized that the observed multimodality is associated with subcellular localization-specific differences in local pI distributions. Here, we overview the multimodality of proteome-wide pI distributions in different organisms focusing on the relationships between protein pI and subcellular localization. We also discuss the probable factors responsible for variation of the intracellular localization-specific pI profiles.

Cell-Free Protein Synthesis of Small Intrinsically Disordered Proteins for NMR Spectroscopy

Cell-free protein synthesis (CFPS) is an established method to produce recombinant proteins and has been used in a wide variety of applications. The use of CFPS has almost from the onset been favorably linked to the production of isotopically labelled proteins for NMR spectroscopy as the resulting labelling of the produced protein is defined by the chosen amino acids during reaction setup. Here we describe how to set up production and isotopic labelling of small intrinsically disordered proteins (IDPs) for NMR spectroscopy applications using an E. coli-based CFPS system in batch mode.

Novel molecular aspects of the CRISPR backbone protein ‘Cas7’ from cyanobacteria

The cyanobacterium Anabaena PCC 7120 shows the presence of Type I-D CRISPR system that can potentially confer adaptive immunity. The Cas7 protein (Alr1562), which forms the backbone of the type I-D surveillance complex, was characterized from Anabaena. Alr1562, showed the presence of the non-canonical RNA recognition motif and two intrinsically disordered regions (IDRs). When overexpressed in E. coli, the Alr1562 protein was soluble and could be purified by affinity chromatography, however, deletion of IDRs rendered Alr1562 completely insoluble. The purified Alr1562 was present in the dimeric or a RNA-associated higher oligomeric form, which appeared as spiral structures under electron microscope. With RNaseA and NaCl treatment, the higher oligomeric form converted to the lower oligomeric form, indicating that oligomerization occurred due to the association of Alr1562 with RNA. The secondary structure of both these forms was largely similar, resembling that of a partially folded protein. The dimeric Alr1562 was more prone to temperature-dependent aggregation than the higher oligomeric form. In vitro, the Alr1562 bound more specifically to a minimal CRISPR unit than to the non-specific RNA. Residues required for binding of Alr1562 to RNA, identified by protein modeling-based approaches, were mutated for functional validation. Interestingly, these mutant proteins, showing reduced ability to bind RNA were predominantly present in dimeric form. Alr1562 was detected with specific antiserum in Anabaena, suggesting that the type I-D system is expressed and may be functional in vivo. This is the first report that describes the characterization of a Cas protein from any photosynthetic organism.

Macromolecular crowding effects on transcription and translation are regulated by free magnesium ion

Cell-free metabolic engineering is an emerging and promising alternative platform for the production of fuels and chemicals. In recent years, macromolecular crowding effect, which is an important function in living cells but ignored in cell-free systems, has been transferred to cell-free protein synthesis (CFPS). However, inhibitory effects of crowding agents on CFPS were frequently observed, and the mechanism is unclear. In this study, free Mg²⁺ was found to be a key factor that can regulate the macromolecular crowding effect on in vitro transcription, in vitro translation, and coupled transcript/translation. Addition of crowding agents (20% of Ficoll-70 or Ficoll-400) enhanced in vitro transcription at an index of free Mg²⁺ concentration (IFMC) below 2 mM but inhibited the transcription when the IFMC was higher than 2 mM. Similarly, Ficoll-400 enhanced in vitro translation and coupled transcription/translation at a lower IFMC (0.1-2 mM) and inhibited the reactions at higher IFMC (>2 mM). Based on the results, CFPS systems could be further optimized by adjusting the content of crowding agents and the IFMC. Besides, the results also indicate that macromolecular crowding effect is important for maintaining the efficiency of in vivo transcription and translation which occur at a low intracellular IFMC (<1 mM).

Localization-specific distributions of protein pI in human proteome are governed by local pH and membrane charge

Background

Whole-proteome distributions of protein isoelectric point (pI) values in different organisms are bi- or trimodal with some variations. It was suggested that the observed multimodality of the proteome-wide pI distributions is associated with subcellular localization-specific differences in the local pI distributions. However, the factors responsible for variation of the intracellular localization-specific pI profiles have not been investigated in detail.

Results

In this work, we explored proteome-wide pI distributions of 32,138 human proteins predicted to reside in 10 subcellular compartments, as well as the pI distributions of experimentally observed lysosomal and Golgi proteins. The distributions were found to differ significantly, although all of them adhered to the major recurrent bimodal pattern. Grossly, acid-biased and alkaline-biased patterns with various minor statistical features were observed at different subcellular locations. Bioinformatics analysis revealed the existence of strong statistically significant correlations between protein pI and subcellular localization. Most markedly, protein pI was found to correlate positively with nuclear and mitochondrial locations and negatively with cytoskeletal, cytoplasmic, lysosomal and peroxisomal environment. Further analysis demonstrated that subcellular compartment-specific pI distributions are greatly influenced by local pH and organelle membrane charge. Multiple nonlinear regression analysis identified a polynomial function of the two variables that best fitted the mean pI values of the localization-specific pI distributions. A high coefficient of determination calculated for this regression (R² = 0.98) suggests that local pH and organelle membrane charge are the major factors responsible for variation of the intracellular localization-specific pI profiles.

Conclusions

Our study demonstrates that strong correlations exist between protein pI and subcellular localization. The specific pI distributions at different subcellular locations are defined by local environment. Predominantly, it is the local pH and membrane charge that shape the organelle-specific protein pI patterns. These findings expand our understanding of spatial organization of the human proteome.

Electronic supplementary material

The online version of this article (10.1186/s12860-019-0221-4) contains supplementary material, which is available to authorized users.

Troubleshooting Guide to Expressing Intrinsically Disordered Proteins for Use in NMR Experiments

Intrinsically disordered proteins (IDPs) represent a structural class of proteins that do not have a well-defined, 3D fold in solution, and often have little secondary structure. To characterize their function and molecular mechanism, it is helpful to examine their structure using nuclear magnetic resonance (NMR), which can report on properties, such as residual structure (at both the secondary and tertiary levels), ligand binding affinity, and the effect of ligand binding on IDP structure, all on a per residue basis. This brief review reports on the common problems and decisions that are involved when preparing a disordered protein for NMR studies. The paper covers gene design, expression host choice, protein purification, and the initial NMR experiments that are performed. While many of these steps are essentially identical to those for ordered proteins, a few key differences are highlighted, including the extreme sensitivity of IDPs to proteolytic cleavage, the ability to use denaturing conditions without having to refold the protein, the optimal chromatographic system choice, and the challenges of quantifying an IDP. After successful purification, characterization by NMR can be done using the standard ¹⁵N-heteronuclear single quantum coherence (¹⁵N-HSQC) experiment, or the newer CON series of experiments that are superior for disordered proteins.

Content of intrinsic disorder influences the outcome of cell-free protein synthesis

Cell-free protein synthesis is used to produce proteins with various structural traits. Recent bioinformatics analyses indicate that more than half of eukaryotic proteins possess long intrinsically disordered regions. However, no systematic study concerning the connection between intrinsic disorder and expression success of cell-free protein synthesis has been presented until now. To address this issue, we examined correlations of the experimentally observed cell-free protein expression yields with the contents of intrinsic disorder bioinformatically predicted in the expressed sequences. This analysis revealed strong relationships between intrinsic disorder and protein amenability to heterologous cell-free expression. On the one hand, elevated disorder content was associated with the increased ratio of soluble expression. On the other hand, overall propensity for detectable protein expression decreased with disorder content. We further demonstrated that these tendencies are rooted in some distinct features of intrinsically disordered regions, such as low hydrophobicity, elevated surface accessibility and high abundance of sequence motifs for proteolytic degradation, including sites of ubiquitination and PEST sequences. Our findings suggest that identification of intrinsically disordered regions in the expressed amino acid sequences can be of practical use for predicting expression success and optimizing cell-free protein synthesis.

Functional advantages of dynamic protein disorder

Intrinsically disordered proteins participate in many important cellular regulatory processes. The absence of a well-defined structure in the free state of a disordered domain, and even on occasion when it is bound to physiological partners, is fundamental to its function. Disordered domains are frequently the location of multiple sites for post-translational modification, the key element of metabolic control in the cell. When a disordered domain folds upon binding to a partner, the resulting complex buries a far greater surface area than in an interaction of comparably-sized folded proteins, thus maximizing specificity at modest protein size. Disorder also maintains accessibility of sites for post-translational modification. Because of their inherent plasticity, disordered domains frequently adopt entirely different structures when bound to different partners, increasing the repertoire of available interactions without the necessity for expression of many different proteins. This feature also adds to the faithfulness of cellular regulation, as the availability of a given disordered domain depends on competition between various partners relevant to different cellular processes.

Plant-PrAS: a database of physicochemical and structural properties and novel functional regions in plant proteomes

Arabidopsis thaliana is an important model species for studies of plant gene functions. Research on Arabidopsis has resulted in the generation of high-quality genome sequences, annotations and related post-genomic studies. The amount of annotation, such as gene-coding regions and structures, is steadily growing in the field of plant research. In contrast to the genomics resource of animals and microorganisms, there are still some difficulties with characterization of some gene functions in plant genomics studies. The acquisition of information on protein structure can help elucidate the corresponding gene function because proteins encoded in the genome possess highly specific structures and functions. In this study, we calculated multiple physicochemical and secondary structural parameters of protein sequences, including length, hydrophobicity, the amount of secondary structure, the number of intrinsically disordered regions (IDRs) and the predicted presence of transmembrane helices and signal peptides, using a total of 208,333 protein sequences from the genomes of six representative plant species, Arabidopsis thaliana, Glycine max (soybean), Populus trichocarpa (poplar), Oryza sativa (rice), Physcomitrella patens (moss) and Cyanidioschyzon merolae (alga). Using the PASS tool and the Rosetta Stone method, we annotated the presence of novel functional regions in 1,732 protein sequences that included unannotated sequences from the Arabidopsis and rice proteomes. These results were organized into the Plant Protein Annotation Suite database (Plant-PrAS), which can be freely accessed online at http://plant-pras.riken.jp/.

Identification of multiple physicochemical and structural properties associated with soluble expression of eukaryotic proteins in cell-free bacterial extracts

Bacterial extracts are widely used to synthesize recombinant proteins. Vast data volumes have been accumulated in cell-free expression databases, covering a whole range of existing proteins. It makes possible comprehensive bioinformatics analysis and identification of multiple features associated with protein solubility and aggregation. In the present paper, an approach to identify the multiple physicochemical and structural properties of amino acid sequences associated with soluble expression of eukaryotic proteins in cell-free bacterial extracts is presented. The method includes: (1) categorical assessment of expression data; (2) calculation and prediction of multiple properties of expressed sequences; (3) correlation of the individual properties with the expression scores; and (4) evaluation of statistical significance of the observed correlations. Using this method, a number of significant correlations between calculated and predicted properties of amino acid sequences and their propensity for soluble cell-free expression have been revealed.

Bioinformatics analysis and optimization of cell-free protein synthesis

Cell-free protein synthesis offers substantial advantages over cell-based expression, allowing direct access to the protein synthetic reaction and meticulous control over the reaction conditions. Recently, we identified a number of statistically significant correlations between calculated and predicted properties of amino acid sequences and their amenability to heterologous cell-free expression. These correlations can be of practical use for predicting expression success and optimizing cell-free protein synthesis. In this chapter, we describe our approach and demonstrate how computational and predictive bioinformatics can be used to analyze and optimize cell-free protein expression.

Electrostatic mis-interactions cause overexpression toxicity of proteins in E. coli

A majority of E. coli proteins when overexpressed inhibit its growth, but the reasons behind overexpression toxicity of proteins remain unknown. Understanding the mechanism of overexpression toxicity is important from evolutionary, biotechnological and possibly clinical perspectives. Here we study sequence and functional features of cytosolic proteins of E. coli associated with overexpression toxicity to understand its mechanism. We find that number of positively charged residues is significantly higher in proteins showing overexpression toxicity. Very long proteins also show high overexpression toxicity. Among the functional classes, transcription factors and regulatory proteins are enriched in toxic proteins, while catalytic proteins are depleted. Overexpression toxicity could be predicted with reasonable accuracy using these few properties. The importance of charged residues in overexpression toxicity indicates that nonspecific electrostatic interactions resulting from protein overexpression cause toxicity of these proteins and suggests ways to improve the expression level of native and foreign proteins in E. coli for basic research and biotechnology. These results might also be applicable to other bacterial species.

Highly efficient NMR assignment of intrinsically disordered proteins: application to B- and T cell receptor domains

We present an integrated approach for efficient characterization of intrinsically disordered proteins. Batch cell-free expression, fast data acquisition, automated analysis, and statistical validation with data resampling have been combined for achieving cost-effective protein expression, and rapid automated backbone assignment. The new methodology is applied for characterization of five cytosolic domains from T- and B-cell receptors in solution.

Multiple post-translational modifications affect heterologous protein synthesis

Post-translational modifications (PTMs) are required for proper folding of many proteins. The low capacity for PTMs hinders the production of heterologous proteins in the widely used prokaryotic systems of protein synthesis. Until now, a systematic and comprehensive study concerning the specific effects of individual PTMs on heterologous protein synthesis has not been presented. To address this issue, we expressed 1488 human proteins and their domains in a bacterial cell-free system, and we examined the correlation of the expression yields with the presence of multiple PTM sites bioinformatically predicted in these proteins. This approach revealed a number of previously unknown statistically significant correlations. Prediction of some PTMs, such as myristoylation, glycosylation, palmitoylation, and disulfide bond formation, was found to significantly worsen protein amenability to soluble expression. The presence of other PTMs, such as aspartyl hydroxylation, C-terminal amidation, and Tyr sulfation, did not correlate with the yield of heterologous protein expression. Surprisingly, the predicted presence of several PTMs, such as phosphorylation, ubiquitination, SUMOylation, and prenylation, was associated with the increased production of properly folded soluble proteins. The plausible rationales for the existence of the observed correlations are presented. Our findings suggest that identification of potential PTMs in polypeptide sequences can be of practical use for predicting expression success and optimizing heterologous protein synthesis. In sum, this study provides the most compelling evidence so far for the role of multiple PTMs in the stability and solubility of heterologously expressed recombinant proteins.

Cell-free protein expression under macromolecular crowding conditions

Background

Cell-free protein expression (CFPE) comprised of in vitro transcription and translation is currently manipulated in relatively dilute solutions, in which the macromolecular crowding effects present in living cells are largely ignored. This may not only affect the efficiency of protein synthesis in vitro, but also limit our understanding of the functions and interactions of biomolecules involved in this fundamental biological process.

Methodology/Principal Findings

Using cell-free synthesis of Renilla luciferase in wheat germ extract as a model system, we investigated the CFPE under macromolecular crowding environments emulated with three different crowding agents: PEG-8000, Ficoll-70 and Ficoll-400, which vary in chemical properties and molecular size. We found that transcription was substantially enhanced in the macromolecular crowding solutions; up to 4-fold increase in the mRNA production was detected in the presence of 20% (w/v) of Ficoll-70. In contrast, translation was generally inhibited by the addition of each of the three crowding agents. This might be due to PEG-induced protein precipitation and non-specific binding of translation factors to Ficoll molecules. We further explored a two-stage CFPE in which transcription and translation was carried out under high then low macromolecular crowding conditions, respectively. It produced 2.2-fold higher protein yield than the coupled CFPE control. The macromolecular crowding effects on CFPE were subsequently confirmed by cell-free synthesis of an approximately two-fold larger protein, Firefly luciferase, under macromolecular crowding environments.

Conclusions/Significance

Three macromolecular crowding agents used in this research had opposite effects on transcription and translation. The results of this study should aid researchers in their choice of macromolecular crowding agents and shows that two-stage CFPE is more efficient than coupled CFPE.

Gene optimization mechanisms: a multi-gene study reveals a high success rate of full-length human proteins expressed in Escherichia coli

The genetic code is universal, but recombinant protein expression in heterologous systems is often hampered by divergent codon usage. Here, we demonstrate that reprogramming by standardized multi-parameter gene optimization software and de novo gene synthesis is a suitable general strategy to improve heterologous protein expression. This study compares expression levels of 94 full-length human wt and sequence-optimized genes coding for pharmaceutically important proteins such as kinases and membrane proteins in E. coli. Fluorescence-based quantification revealed increased protein yields for 70% of in vivo expressed optimized genes compared to the wt DNA sequences and also resulted in increased amounts of protein that can be purified. The improvement in transgene expression correlated with higher mRNA levels in our analyzed examples. In all cases tested, expression levels using wt genes in tRNA-supplemented bacterial strains were outperformed by optimized genes expressed in non-supplemented host cells.