Barcode fusion genetics-protein-fragment complementation assay (BFG-PCA): tools and resources that expand the potential for binary protein interaction discovery.

Barcode fusion genetics (BFG) utilizes deep sequencing to improve the throughput of protein-protein interaction (PPI) screening in pools. BFG has been implemented in Yeast two-hybrid (Y2H) screens (BFG-Y2H). While Y2H requires test protein pairs to localize in the nucleus for reporter reconstruction, dihydrofolate reductase protein-fragment complementation assay (DHFR-PCA) allows proteins to localize in broader subcellular contexts and proves to be largely orthogonal to Y2H. Here, we implemented BFG to DHFR-PCA (BFG-PCA). This plasmid-based system can leverage ORF collections across model organisms to perform comparative analysis, unlike the original DHFR-PCA that requires yeast genomic integration. The scalability and quality of BFG-PCA were demonstrated by screening human and yeast interactions for >11 000 bait-prey pairs. BFG-PCA showed high-sensitivity and high-specificity for capturing known interactions for both species. BFG-Y2H and BFG-PCA capture distinct sets of PPIs, which can partially be explained based on the domain orientation of the reporter tags. BFG-PCA is a high-throughput protein interaction technology to interrogate binary PPIs that exploits clone collections from any species of interest, expanding the scope of PPI assays.

Genome Scale-Differential Flux Analysis reveals deregulation of lung cell metabolism on SARS-CoV-2 infection.

The COVID-19 pandemic is posing an unprecedented threat to the whole world. In this regard, it is absolutely imperative to understand the mechanism of metabolic reprogramming of host human cells by SARS-CoV-2. A better understanding of the metabolic alterations would aid in design of better therapeutics to deal with COVID-19 pandemic. We developed an integrated genome-scale metabolic model of normal human bronchial epithelial cells (NHBE) infected with SARS-CoV-2 using gene-expression and macromolecular make-up of the virus. The reconstructed model predicts growth rates of the virus in high agreement with the experimental measured values. Furthermore, we report a method for conducting genome-scale differential flux analysis (GS-DFA) in context-specific metabolic models. We apply the method to the context-specific model and identify severely affected metabolic modules predominantly comprising of lipid metabolism. We conduct an integrated analysis of the flux-altered reactions, host-virus protein-protein interaction network and phospho-proteomics data to understand the mechanism of flux alteration in host cells. We show that several enzymes driving the altered reactions inferred by our method to be directly interacting with viral proteins and also undergoing differential phosphorylation under diseased state. In case of SARS-CoV-2 infection, lipid metabolism particularly fatty acid oxidation, cholesterol biosynthesis and beta-oxidation cycle along with arachidonic acid metabolism are predicted to be most affected which confirms with clinical metabolomics studies. GS-DFA can be applied to existing repertoire of high-throughput proteomic or transcriptomic data in diseased condition to understand metabolic deregulation at the level of flux.

Probing the molecular mechanism of aggressive infection by antimony resistant Leishmania donovani.

The disease visceral leishmaniasis (VL) or kala azar is caused by the protozoan parasite, Leishmania donovani (LD). For many decades the pentavalent antimonial drugs countered the successive epidemics of the disease in the Indian sub-continent and elsewhere. With time, antimony resistant LD (LDR) developed and the drug in turn lost its efficacy. Infection of mammals with LDR gives rise to aggressive infection as compared to its sensitive counterpart (LDS) coupled with higher surge of IL-10 and TGF-ß. The IL-10 causes upregulation of multidrug resistant protein-1 which causes efflux of antimonials from LDR infected cells. This is believed to be a key mechanism of antimony resistance. MicroRNAs (miRNAs) are tiny post-transcriptional regulators of gene expression in mammalian cells and in macrophage play a pivotal role in controlling the expression of cytokines involved in infection process. Therefore, a change in miRNA profiles of macrophages infected with LDS or LDR could explain the differential cytokine response observed. Interestingly, the outcome of LD infection is also governed by the critical balance of pro- and anti-inflammatory cytokines which is inturn regulated by miRNA-Ago2 or miRNP complex and its antagonist RNA binding protein HuR. Here Ago2 plays the fulcrum whose phosphorylation and de-phosphorylation dictates the process; which in turn is controlled by PP2A and HuR. LDS and LDR upregulate PP2A and downregulate HuR at different magnitude leading to various levels of anti-inflammatory to proinflammatory cytokine production and resulting pathology in the host. While ectopic HuR expression alone is sufficient to clear LDS infection, simultaneous upregulation of HuR and inhibition of PP2A is required to inhibit LDR mediated infection. Therefore, tampering with miRNA pathway could be a new strategy to control infection caused by LDR parasite.

A dynamic network model predicts the phenotypes of multicellular clusters from cellular properties.

Cell division without cell separation produces multicellular clusters in budding yeast. Two fundamental characteristics of these clusters are their size (the number of cells per cluster) and cellular composition: the fractions of cells with different phenotypes. Using cells as nodes and links between mother and daughter cells as edges, we model cluster growth and breakage by varying three parameters: the cell division rate, the rate at which intercellular connections break, and the kissing number (the maximum number of connections to one cell). We find that the kissing number sets the maximum possible cluster size. Below this limit, the ratio of the cell division rate to the connection breaking rate determines the cluster size. If links have a constant probability of breaking per unit time, the probability that a link survives decreases exponentially with its age. Modeling this behavior recapitulates experimental data. We then use this framework to examine synthetic, differentiating clusters with two cell types, faster-growing germ cells and their somatic derivatives. The fraction of clusters that contain both cell types increases as either of two parameters increase: the kissing number and difference between the growth rate of germ and somatic cells. In a population of clusters, the variation in cellular composition is inversely correlated (r2 = 0.87) with the average fraction of somatic cells in clusters. Our results show how a small number of cellular features can control the phenotypes of multicellular clusters that were potentially the ancestors of more complex forms of multicellular development, organization, and reproduction.

Alternating selection for dispersal and multicellularity favors regulated life cycles.

The evolution of complex multicellularity opened paths to increased morphological diversity and organizational novelty. This transition involved three processes: cells remained attached to one another to form groups, cells within these groups differentiated to perform different tasks, and the groups evolved new reproductive strategies.1,2,3,4,5 Recent experiments identified selective pressures and mutations that can drive the emergence of simple multicellularity and cell differentiation,6,7,8,9,10,11 but the evolution of life cycles, particularly how simple multicellular forms reproduce, has been understudied. The selective pressure and mechanisms that produced a regular alternation between single cells and multicellular collectives are still unclear.12 To probe the factors regulating simple multicellular life cycles, we examined a collection of wild isolates of the budding yeast S. cerevisiae.12,13 We found that all these strains can exist as multicellular clusters, a phenotype that is controlled by the mating-type locus and strongly influenced by the nutritional environment. Inspired by this variation, we engineered inducible dispersal in a multicellular laboratory strain and demonstrated that a regulated life cycle has an advantage over constitutively single-celled or constitutively multicellular life cycles when the environment alternates between favoring intercellular cooperation (a low sucrose concentration) and dispersal (a patchy environment generated by emulsion). Our results suggest that the separation of mother and daughter cells is under selection in wild isolates and is regulated by their genetic composition and the environments they encounter and that alternating patterns of resource availability may have played a role in the evolution of life cycles.

Multicellular growth as a dynamic network of cells.

Cell division without cell separation produces multicellular clusters in budding yeast. Two fundamental characteristics of these clusters are their size (the number of cells per cluster) and cellular composition: the fractions of cells with different phenotypes. However, we do not understand how different cellular features quantitatively influence these two phenotypes. Using cells as nodes and links between mother and daughter cells as edges, we model cluster growth and breakage by varying three parameters: the cell division rate, the rate at which intercellular connections break, and the kissing number (the maximum number of connections to one cell). We find that the kissing number sets the maximum possible cluster size. Below this limit, the ratio of the cell division rate to the connection breaking rate determines the cluster size. If links have a constant probability of breaking per unit time, the probability that a link survives decreases exponentially with its age. Modeling this behavior recapitulates experimental data. We then use this framework to examine synthetic, differentiating clusters with two cell types, faster-growing germ cells and their somatic derivatives. The fraction of clusters that contain both cell types increases as either of two parameters increase: the kissing number and difference between the growth rate of germ and somatic cells. In a population of clusters, the variation in cellular composition is inversely correlated (r2=0.87) with the average fraction of somatic cells in clusters. Our results show how a small number of cellular features can control the phenotypes of multicellular clusters that were potentially the ancestors of more complex forms of multicellular development, organization, and reproduction.

Reconstruction and analysis of genome-scale metabolic model of weak Crabtree positive yeast Lachancea kluyveri.

Lachancea kluyveri, a weak Crabtree positive yeast, has been extensively studied for its unique URC pyrimidine catabolism pathway. It produces more biomass than Saccharomyces cerevisiae due to the underlying weak Crabtree effect and resorts to fermentation only in oxygen limiting conditions that renders it as a suitable industrial host. The yeast also produces ethyl acetate as a major overflow metabolite in aerobic conditions. Here, we report the first genome-scale metabolic model, iPN730, of L. kluyveri comprising of 1235 reactions, 1179 metabolites, and 730 genes distributed in 8 compartments. The in silico viability in different media conditions and the growth characteristics in various carbon sources show good agreement with experimental data. Dynamic flux balance analysis describes the growth dynamics, substrate utilization and product formation kinetics in various oxygen-limited conditions. We have also demonstrated the effect of switching carbon sources on the production of ethyl acetate under varying oxygen uptake rates. A phenotypic phase plane analysis described the energetic cost penalty of ethyl acetate and ethanol production on the specific growth rate of L. kluyveri. We generated the context specific models of L. kluyveri growing on uracil or ammonium salts as the sole nitrogen source. Differential flux calculated using flux variability analysis helped us in highlighting pathways like purine, histidine, riboflavin and pyrimidine metabolism associated with uracil degradation. The genome-scale metabolic construction of L. kluyveri will provide a better understanding of metabolism behind ethyl acetate production as well as uracil catabolism (pyrimidine degradation) pathway. iPN730 is an addition to genome-scale metabolic models of non-conventional yeasts that will facilitate system-wide omics analysis to understand fungal metabolic diversity.