Mechanistic models of how single cells respond to different perturbations can help integrate disparate big data sets or predict response to varied drug combinations. However, the construction and simulation of such models have proved challenging. Here, we developed a python-based model creation and simulation pipeline that converts a few structured text files into an SBML standard and is high-performance- and cloud-computing ready. We applied this pipeline to our large-scale, mechanistic pan-cancer signaling model (named SPARCED) and demonstrate it by adding an IFNγ pathway submodel. We then investigated whether a putative crosstalk mechanism could be consistent with experimental observations from the LINCS MCF10A Data Cube that IFNγ acts as an anti-proliferative factor. The analyses suggested this observation can be explained by IFNγ-induced SOCS1 sequestering activated EGF receptors. This work forms a foundational recipe for increased mechanistic model-based data integration on a single-cell level, an important building block for clinically-predictive mechanistic models.

The current era of big biomedical data accumulation and availability brings data integration opportunities for leveraging its totality to make new discoveries and/or clinically predictive models. Black-box statistical and machine learning methods are powerful for such integration, but often cannot provide mechanistic reasoning, particularly on the single-cell level. While single-cell mechanistic models clearly enable such reasoning, they are predominantly “small-scale”, and struggle with the scalability and reusability required for meaningful data integration. Here, we present an open-source pipeline for scalable, single-cell mechanistic modeling from simple, annotated input files that can serve as a foundation for mechanistic data integration. As a test case, we convert one of the largest existing single-cell mechanistic models to this format, demonstrating robustness and reproducibility of the approach. We show that the model cell line context can be changed with simple replacement of input file parameter values. We next use this new model to test alternative mechanistic hypotheses for the experimental observations that interferon-gamma (IFNG) inhibits epidermal growth factor (EGF)-induced cell proliferation. Model- based analysis suggested, and experiments support that these observations are better explained by IFNG-induced SOCS1 expression sequestering activated EGF receptors, thereby downregulating AKT activity, as opposed to direct IFNG-induced upregulation of p21 expression. Overall, this new pipeline enables large-scale, single-cell, and mechanistically-transparent modeling as a data integration modality complementary to machine learning.

Western blotting is a widely used technique for molecular-weight-resolved analysis of proteins and their posttranslational modifications, but high-throughput implementations of the standard slab gel arrangement are scarce. The previously developed Microwestern requires a piezoelectric pipetting instrument, which is not available for many labs. Here, we report the Mesowestern blot, which uses a 3D-printable gel casting mold to enable high-throughput Western blotting without piezoelectric pipetting and is compatible with the standard sample preparation and small (∼1 μL) sample sizes. The main tradeoffs are reduced molecular weight resolution and higher sample-to-sample CV, making it suitable for qualitative screening applications. The casted polyacrylamide gel contains 336, ∼0.5 μL micropipette-loadable sample wells arranged within a standard microplate footprint. Polyacrylamide % can be altered to change molecular weight resolution profiles. Proof-of-concept experiments using both infrared-fluorescent molecular weight protein ladder and cell lysate (RIPA buffer) demonstrate that the protein loaded in Mesowestern gels is amenable to the standard Western blotting steps. The main difference between Mesowestern and traditional Western is that semidry horizontal instead of immersed vertical gel electrophoresis is used. The linear range of detection is at least 32-fold, and at least ∼500 attomols of β-actin can be detected (∼29 ng of total protein from mammalian cell lysates: ∼100–300 cells). Because the gel mold is 3D-printable, users with access to additive manufacturing cores have significant design freedom for custom layouts. We expect that the technique could be easily adopted by any typical cell and molecular biology laboratory already performing Western blots.

Systematic, genome-scale genetic screens have been instrumental for elucidating genotype-phenotype relationships, but approaches for probing genetic interactions have been limited to at most ~100 pre-selected gene combinations in mammalian cells. Here, we introduce a theory for high-throughput genetic interaction screens. The theory extends our recently developed Multiplexing using Spectral Imaging and Combinatorics (MuSIC) approach to propose ~105 spectrally unique, genetically-encoded MuSIC barcodes from 18 currently available fluorescent proteins. Simulation studies based on constraints imposed by spectral flow cytometry equipment suggest that genetic interaction screens at the human genome-scale may be possible if MuSIC barcodes can be paired to guide RNAs. While experimental testing of this theory awaits, it offers transformative potential for genetic perturbation technology and knowledge of genetic function. More broadly, the availability of a genome-scale spectral barcode library for non-destructive identification of single-cells could find more widespread applications such as traditional genetic screening and high-dimensional lineage tracing.

Western blotting is a widely-used technique for molecular-weight-resolved analysis of proteins and their post-translational modifications, but has been refractory to affordable scale-up. Here, we report the Mesowestern blot, which uses a 3D-printable gel-casting mold to enable affordable, high-throughput Western blotting with standard sample preparation and small (<1 uL) sample sizes. The casted polyacrylamide gel contains 336, 0.5 uL micropipette-loadable sample wells arranged within a standard microplate footprint. Polyacrylamide % can be altered to change molecular weight resolution range. Proof-of-concept experiments using both infrared-fluorescent molecular weight protein ladder as well as cell lysate (RIPA buffer) demonstrate protein loaded in Mesowestern gels is amenable to the standard Western blotting steps. The main difference between Mesowestern and traditional Western is that semi-dry horizontal instead of immersed vertical gel electrophoresis is used. The linear range of detection is approximately 2 orders of magnitude, with a limit of detection (for β-actin) of around 30 ng of total protein from mammalian cell lysates (~30-3000 cells). Because the gel mold is 3D-printable, users have significant design freedom for custom layouts, and there are few barriers to adoption by the typical cell and molecular biology laboratory already performing Western blots.