Here, we propose a framework for the design of synthetic protein networks from modular protein-protein or protein-peptide interactions and provide a starter toolkit of protein building blocks. Our proof of concept experiments outline a general work flow for part-based protein systems engineering. We streamlined the iterative BioBrick cloning protocol and assembled 25 synthetic multidomain proteins each from seven standardized DNA fragments. A systematic screen revealed two main factors controlling protein expression in Escherichia coli: obstruction of translation initiation by mRNA secondary structure or toxicity of individual domains. Eventually, 13 proteins were purified for further characterization. Starting from well-established biotechnological tools, two general-purpose interaction input and two readout devices were built and characterized in vitro. Constitutive interaction input was achieved with a pair of synthetic leucine zippers. The second interaction was drug-controlled utilizing the rapamycin-induced binding of FRB(T2098L) to FKBP12. The interaction kinetics of both devices were analyzed by surface plasmon resonance. Readout was based on Förster resonance energy transfer between fluorescent proteins and was quantified for various combinations of input and output devices. Our results demonstrate the feasibility of parts-based protein synthetic biology. Additionally, we identify future challenges and limitations of modular design along with approaches to address them.

To explore the conservation of Src homology 3 (SH3) domain-mediated networks in evolution, we compared the specificity landscape of these domains among four yeast species, Saccharomyces cerevisiae, Ashbya gossypii, Candida albicans, and Schizosaccharomyces pombe, encompassing 400 million years of evolution. We first aligned and catalogued the families of SH3-containing proteins in these four species to determine the relationships between homologous domains. Then, we tagged and purified all soluble SH3 domains (82 in total) to perform a quantitative peptide assay (SPOT) for each SH3 domain. All SPOT readouts were hierarchically clustered and we observed that the organization of the SH3 specificity landscape in three distinct profile classes remains conserved across these four yeast species. We also produced a specificity profile for each SH3 domain from manually aligned top SPOT hits and compared the within-family binding motif consensus. This analysis revealed a striking example of binding motif divergence in a C. albicans Rvs167 paralog, which cannot be explained by overall SH3 sequence or interface residue divergence, and we validated this specificity change with a yeast two-hybrid (Y2H) assay. In addition, we show that position-weighted matrices (PWM) compiled from SPOT assays can be used for binding motif screening in potential binding partners and present cases where motifs are either conserved or lost among homologous SH3 interacting proteins. Finally, by comparing pairwise SH3 sequence identity to binding profile correlation we show that for ~75% of all analyzed families the SH3 specificity profile was remarkably conserved over a large evolutionary distance. Thus, a high sequence identity within an SH3 domain family predicts conserved binding specificity, whereas divergence in sequence identity often coincided with a change in binding specificity within this family. As such, our results are important for future studies aimed at unraveling complex specificity networks of peptide recognition domains in higher eukaryotes, including mammals.

This Pulmonary Perspective summarizes the content and main conclusions of an international workshop on personalized respiratory medicine coorganized by the Barcelona Respiratory Network ( www.brn.cat ) and the AJRCCM in June 2014. It discusses (1) its definition and historical, social, legal, and ethical aspects; (2) the view from different disciplines, including basic science, epidemiology, bioinformatics, and network/systems medicine; (3) the bottlenecks and opportunities identified by some currently ongoing projects; and (4) the implications for the individual, the healthcare system and the pharmaceutical industry. The authors hope that, although it is not a systematic review on the subject, this document can be a useful reference for researchers, clinicians, healthcare managers, policy-makers, and industry parties interested in personalized respiratory medicine.

Copy number alterations (CNAs) in cancer patients show a large variability in their number, length and position, but the sources of this variability are not known. CNA number and length are linked to patient survival, suggesting clinical relevance. We have identified genes that tend to be mutated in samples that have few or many CNAs, which we term CONIM genes (COpy Number Instability Modulators). CONIM proteins cluster into a densely connected subnetwork of physical interactions and many of them are epigenetic modifiers. Therefore, we investigated how the epigenome of the tissue-of-origin influences the position of CNA breakpoints and the properties of the resulting CNAs. We found that the presence of heterochromatin in the tissue-of-origin contributes to the recurrence and length of CNAs in the respective cancer type.

Retinopathies are a group of monogenetic or complex retinal diseases associated with high unmet medical need. Monogenic disorders are caused by rare genetic variation and usually arise early in life. Other diseases, such as age-related macular degeneration (AMD), develop late in life and are considered to be of complex origin as they develop from a combination of genetic, ageing, environmental and lifestyle risk factors. Here, we contrast the underlying disease networks and pathological mechanisms of monogenic as opposed to complex retinopathies, using AMD as an example of the latter. We show that, surprisingly, genes associated with the different forms of retinopathies in general do not overlap despite their overlapping retinal phenotypes. Further, AMD risk genes participate in multiple networks with interaction partners that link to different ubiquitous pathways affecting general tissue integrity and homeostasis. Thus AMD most likely represents an endophenotype with differing underlying pathogenesis in different subjects. Localising these pathomechanisms and processes within and across different retinal anatomical compartments provides a novel representation of AMD that may be extended to complex disease in general. This approach may generate improved treatment options that target multiple processes with the aim of restoring tissue homeostasis and maintaining vision.

Orchestration of signaling, photoreceptor structural integrity, and maintenance needed for mammalian vision remain enigmatic. By integrating three proteomic data sets, literature mining, computational analyses, and structural information, we have generated a multiscale signal transduction network linked to the visual G protein-coupled receptor (GPCR) rhodopsin, the major protein component of rod outer segments. This network was complemented by domain decomposition of protein-protein interactions and then qualified for mutually exclusive or mutually compatible interactions and ternary complex formation using structural data. The resulting information not only offers a comprehensive view of signal transduction induced by this GPCR but also suggests novel signaling routes to cytoskeleton dynamics and vesicular trafficking, predicting an important level of regulation through small GTPases. Further, it demonstrates a specific disease susceptibility of the core visual pathway due to the uniqueness of its components present mainly in the eye. As a comprehensive multiscale network, it can serve as a basis to elucidate the physiological principles of photoreceptor function, identify potential disease-associated genes and proteins, and guide the development of therapies that target specific branches of the signaling pathway.

Mycoplasma pneumoniae, a threatening pathogen with a minimal genome, is a model organism for bacterial systems biology for which substantial experimental information is available. With the goal of understanding the complex interactions underlying its metabolism, we analyzed and characterized the metabolic network of M. pneumoniae in great detail, integrating data from different omics analyses under a range of conditions into a constraint-based model backbone. Iterating model predictions, hypothesis generation, experimental testing, and model refinement, we accurately curated the network and quantitatively explored the energy metabolism. In contrast to other bacteria, M. pneumoniae uses most of its energy for maintenance tasks instead of growth. We show that in highly linear networks the prediction of flux distributions for different growth times allows analysis of time-dependent changes, albeit using a static model. By performing an in silico knock-out study as well as analyzing flux distributions in single and double mutant phenotypes, we demonstrated that the model accurately represents the metabolism of M. pneumoniae. The experimentally validated model provides a solid basis for understanding its metabolic regulatory mechanisms.

We developed a comprehensive resource for the genome-reduced bacterium Mycoplasma pneumoniae comprising 1748 consistently generated '-omics' data sets, and used it to quantify the power of antisense non-coding RNAs (ncRNAs), lysine acetylation, and protein phosphorylation in predicting protein abundance (11%, 24% and 8%, respectively). These factors taken together are four times more predictive of the proteome abundance than of mRNA abundance. In bacteria, post-translational modifications (PTMs) and ncRNA transcription were both found to increase with decreasing genomic GC-content and genome size. Thus, the evolutionary forces constraining genome size and GC-content modify the relative contributions of the different regulatory layers to proteome homeostasis, and impact more genomic and genetic features than previously appreciated. Indeed, these scaling principles will enable us to develop more informed approaches when engineering minimal synthetic genomes.

One of the biggest mysteries in cancer research remains why mutations in certain genes cause cancer only at specific sites in the human body. The poor correlation between the expression level of a cancer gene and the tissues in which it causes malignant transformations raises the question of which factors determine the tissue-specific effects of a mutation. Here, we explore why some cancer genes are associated only with few different cancer types (i.e., are specific), while others are found mutated in a large number of different types of cancer (i.e., are general). We do so by contrasting cellular functions of specific-cancer genes with those of general ones to identify properties that determine where in the body a gene mutation is causing malignant transformations. We identified different groups of cancer genes that did not behave as expected (i.e., DNA repair genes being tissue specific, immune response genes showing a bimodal specificity function or strong association of generally expressed genes to particular cancers). Analysis of these three groups demonstrates the importance of environmental impact for understanding why certain cancer genes are only involved in the development of some cancer types but are rarely found mutated in other types of cancer.

Protein post-translational modifications (PTMs) represent important regulatory states that when combined have been hypothesized to act as molecular codes and to generate a functional diversity beyond genome and transcriptome. We systematically investigate the interplay of protein phosphorylation with other post-transcriptional regulatory mechanisms in the genome-reduced bacterium Mycoplasma pneumoniae. Systematic perturbations by deletion of its only two protein kinases and its unique protein phosphatase identified not only the protein-specific effect on the phosphorylation network, but also a modulation of proteome abundance and lysine acetylation patterns, mostly in the absence of transcriptional changes. Reciprocally, deletion of the two putative N-acetyltransferases affects protein phosphorylation, confirming cross-talk between the two PTMs. The measured M. pneumoniae phosphoproteome and lysine acetylome revealed that both PTMs are very common, that (as in Eukaryotes) they often co-occur within the same protein and that they are frequently observed at interaction interfaces and in multifunctional proteins. The results imply previously unreported hidden layers of post-transcriptional regulation intertwining phosphorylation with lysine acetylation and other mechanisms that define the functional state of a cell.

Many aspects of cell signalling, trafficking, and targeting are governed by interactions between globular protein domains and short peptide segments. These domains often bind multiple peptides that share a common sequence pattern, or "linear motif" (e.g., SH3 binding to PxxP). Many domains are known, though comparatively few linear motifs have been discovered. Their short length (three to eight residues), and the fact that they often reside in disordered regions in proteins makes them difficult to detect through sequence comparison or experiment. Nevertheless, each new motif provides critical molecular details of how interaction networks are constructed, and can explain how one protein is able to bind to very different partners. Here we show that binding motifs can be detected using data from genome-scale interaction studies, and thus avoid the normally slow discovery process. Our approach based on motif over-representation in non-homologous sequences, rediscovers known motifs and predicts dozens of others. Direct binding experiments reveal that two predicted motifs are indeed protein-binding modules: a DxxDxxxD protein phosphatase 1 binding motif with a KD of 22 microM and a VxxxRxYS motif that binds Translin with a KD of 43 microM. We estimate that there are dozens or even hundreds of linear motifs yet to be discovered that will give molecular insight into protein networks and greatly illuminate cellular processes.

Formation of multiprotein complexes is a common theme to pattern a cell, thereby generating spatially and functionally distinct entities at specialised regions. Central components of these complexes are scaffold proteins, which contain several protein-protein interaction domains and provide a platform to recruit a variety of additional components. There is increasing evidence that protein complexes are dynamic structures and that their components can undergo various interactions depending on the cellular context. However, little is known so far about the factors regulating this behaviour. One evolutionarily conserved protein complex, which can be found both in Drosophila and mammalian epithelial cells, is composed of the transmembrane protein Crumbs/Crb3 and the scaffolding proteins Stardust/Pals1 and DPATJ/PATJ, respectively, and localises apically to the zonula adherens. Here we show by in vitro analysis that, similar as in vertebrates, the single PDZ domain of Drosophila DmPar-6 can bind to the four C-terminal amino acids (ERLI) of the transmembrane protein Crumbs. To further evaluate the binding capability of Crumbs to DmPar-6 and the MAGUK protein Stardust, analysis of the PDZ structural database and modelling of the interactions between the C-terminus of Crumbs and the PDZ domains of these two proteins were performed. The results suggest that both PDZ domains bind Crumbs with similar affinities. These data are supported by quantitative yeast two-hybrid interactions. In vivo analysis performed in cell cultures and in the Drosophila embryo show that the cytoplasmic domain of Crumbs can recruit DmPar-6 and DaPKC to the plasma membrane. The data presented here are discussed with respect to possible dynamic interactions between these proteins.

Identification of small open reading frames (smORFs) encoding small proteins (≤ 100 amino acids; SEPs) is a challenge in the fields of genome annotation and protein discovery. Here, by combining a novel bioinformatics tool (RanSEPs) with "-omics" approaches, we were able to describe 109 bacterial small ORFomes. Predictions were first validated by performing an exhaustive search of SEPs present in Mycoplasma pneumoniae proteome via mass spectrometry, which illustrated the limitations of shotgun approaches. Then, RanSEPs predictions were validated and compared with other tools using proteomic datasets from different bacterial species and SEPs from the literature. We found that up to 16 ± 9% of proteins in an organism could be classified as SEPs. Integration of RanSEPs predictions with transcriptomics data showed that some annotated non-coding RNAs could in fact encode for SEPs. A functional study of SEPs highlighted an enrichment in the membrane, translation, metabolism, and nucleotide-binding categories. Additionally, 9.7% of the SEPs included a N-terminus predicted signal peptide. We envision RanSEPs as a tool to unmask the hidden universe of small bacterial proteins.

DNA methylation is an epigenetic mechanism known to affect gene expression and aberrant DNA methylation patterns have been described in cancer. However, only a small fraction of differential methylation events target genes with a defined role in cancer, raising the question of how aberrant DNA methylation contributes to carcinogenesis. As recently a link has been suggested between methylation patterns arising in ageing and those arising in cancer, we asked which aberrations are unique to cancer and which are the product of normal ageing processes. We therefore compared the methylation patterns between ageing and cancer in multiple tissues. We observed that hypermethylation preferentially occurs in regulatory elements, while hypomethylation is associated with structural features of the chromatin. Specifically, we observed consistent hypomethylation of late-replicating, lamina-associated domains. The extent of hypomethylation was stronger in cancer, but in both ageing and cancer it was proportional to the replication timing of the region and the cell division rate of the tissue. Moreover, cancer patients who displayed more hypomethylation in late-replicating, lamina-associated domains had higher expression of cell division genes. These findings suggest that different cell division rates contribute to tissue- and cancer type-specific DNA methylation profiles.

Mycoplasmas are important model organisms for Systems and Synthetic Biology, and are pathogenic to a wide variety of species. Despite their relevance, many of the tools established for genome editing in other microorganisms are not available for Mycoplasmas. The Tn4001 transposon is the reference tool to work with these bacteria, but the transformation efficiencies (TEs) reported for the different species vary substantially. Here, we explore the mechanisms underlying these differences in four Mycoplasma species, Mycoplasma agalactiae, Mycoplasma feriruminatoris, Mycoplasma gallisepticum and Mycoplasma pneumoniae, selected for being representative members of each cluster of the Mycoplasma genus. We found that regulatory regions (RRs) driving the expression of the transposase and the antibiotic resistance marker have a major impact on the TEs. We then designed a synthetic RR termed SynMyco RR to control the expression of the key transposon vector elements. Using this synthetic RR, we were able to increase the TE for M. gallisepticum, M. feriruminatoris and M. agalactiae by 30-, 980- and 1036-fold, respectively. Finally, to illustrate the potential of this new transposon, we performed the first essentiality study in M. agalactiae, basing our study on more than 199,000 genome insertions.

Cancer is a rapidly evolving, multifactorial disease that accumulates numerous genetic and epigenetic alterations. This results in molecular and phenotypic heterogeneity within the tumor, the complexity of which is further amplified through specific interactions between cancer cells. We aimed to dissect the molecular mechanisms underlying the cooperation between different clones.

The C-terminal sequence of a protein is involved in processes such as efficiency of translation termination and protein degradation. However, the general relationship between features of this C-terminal sequence and levels of protein expression remains unknown. Here, we identified C-terminal amino acid biases that are ubiquitous across the bacterial taxonomy (1,582 genomes). We showed that the frequency is higher for positively charged amino acids (lysine, arginine), while hydrophobic amino acids and threonine are lower. We then studied the impact of C-terminal composition on protein levels in a library of Mycoplasma pneumoniae mutants, covering all possible combinations of the two last codons. We found that charged and polar residues, in particular lysine, led to higher expression, while hydrophobic and aromatic residues led to lower expression, with a difference in protein levels up to fourfold. We further showed that modulation of protein degradation rate could be one of the main mechanisms driving these differences. Our results demonstrate that the identity of the last amino acids has a strong influence on protein expression levels.

Microtia is a congenital malformation that encompasses mild hypoplasia to complete loss of the external ear, or pinna. Although the contribution of genetic variation and environmental factors to microtia remains elusive, Amerindigenous populations have the highest reported incidence. Here, using both transmission disequilibrium tests and association studies in microtia trios (parents and affected child) and microtia cohorts enrolled in Latin America, we map an ∼10-kb microtia locus (odds ratio = 4.7; P = 6.78e-18) to the intergenic region between Roundabout 1 (ROBO1) and Roundabout 2 (ROBO2) (chr3: 78546526 to 78555137). While alleles at the microtia locus significantly increase the risk of microtia, their penetrance is low (<1%). We demonstrate that the microtia locus contains a polymorphic complex repeat element that is expanded in affected individuals. The locus is located near a chromatin loop region that regulates ROBO1 and ROBO2 expression in induced pluripotent stem cell–derived neural crest cells. Furthermore, we use single nuclear RNA sequencing to demonstrate ROBO1 and ROBO2 expression in both fibroblasts and chondrocytes of the mature human pinna. Because the microtia allele is enriched in Amerindigenous populations and is shared by some East Asian subjects with craniofacial malformations, we propose that both populations share a mutation that arose in a common ancestor prior to the ancient migration of Eurasian populations into the Americas and that the high incidence of microtia among Amerindigenous populations reflects the population bottleneck that occurred during the migration out of Eurasia.

Regulation of microtubule (MT) dynamics is key for mitotic spindle assembly and faithful chromosome segregation. Here we show that polyglutamylation, a still understudied post-translational modification of spindle MTs, is essential to define their dynamics within the range required for error-free chromosome segregation. We identify TTLL11 as an enzyme driving MT polyglutamylation in mitosis and show that reducing TTLL11 levels in human cells or zebrafish embryos compromises chromosome segregation fidelity and impairs early embryonic development. Our data reveal a mechanism to ensure genome stability in normal cells that is compromised in cancer cells that systematically downregulate TTLL11. Our data suggest a direct link between MT dynamics regulation, MT polyglutamylation and two salient features of tumour cells, aneuploidy and chromosome instability (CIN).

To assist in the advancement of the large-scale production of safe Mycoplasma vaccines and other Mycoplasma-based therapies, we developed a culture medium free of animal serum and other animal components for Mycoplasma pneumoniae growth. By establishing a workflow method to systematically test different compounds and concentrations, we provide optimized formulations capable of supporting serial passaging and robust growth reaching 60 to 70% of the biomass obtained in rich medium. Global transcriptomic and proteomic analysis showed minor physiological changes upon cell culture in the animal component-free medium, supporting its suitability for the production of M. pneumoniae-based therapies. The major contributors to growth performance were found to be glucose as a carbon source, glycerol, cholesterol, and phospholipids as a source of fatty acids. Bovine serum albumin or cyclodextrin (in the animal component-free medium) were required as lipid carriers to prevent lipid toxicity. Connaught Medical Research Laboratories medium (CMRL) used to simplify medium preparation as a source of amino acids, nucleotide precursors, vitamins, and other cofactors could be substituted by cysteine. In fact, the presence of protein hydrolysates such as yeastolate or peptones was found to be essential and preferred over free amino acids, except for the cysteine. Supplementation of nucleotide precursors and vitamins is not strictly necessary in the presence of yeastolate, suggesting that this animal origin-free hydrolysate serves as an efficient source for these compounds. Finally, we adapted the serum-free medium formulation to support growth of Mycoplasma hyopneumoniae, a swine pathogen for which inactivated whole-cell vaccines are available. IMPORTANCE Mycoplasma infections have a significant negative impact on both livestock production and human health. Vaccination is often the first option to control disease and alleviate the economic impact that some Mycoplasma infections cause on milk production, weight gain, and animal health. The fastidious nutrient requirements of these bacteria, however, challenges the industrial production of attenuated or inactivated whole-cell vaccines, which depends on the use of animal serum and other animal raw materials. Apart from their clinical relevance, some Mycoplasma species have become cellular models for systems and synthetic biology, owing to the small size of their genomes and the absence of a cell wall, which offers unique opportunities for the secretion and delivery of biotherapeutics. This study proposes medium formulations free of serum and animal components with the potential of supporting large-scale production upon industrial optimization, thus contributing to the development of safe vaccines and other Mycoplasma-based therapies.