Role of Intercellular Interactions on Single Cell and Population Level Responses: Considerations for Multicellular Bioreporter Design

Douglas Lin^*Michael Martin

• University of Louisville, Louisville, United States

Abstract Bioreporters are genetically engineered cells that produce detectable responses in the presence of specific analytes, providing a cheap, mass-producible, and accurate method of analyte detection. Most research focuses on the single cell-level, where all engineering is concentrated on the interactions within a single cell. However, intercellular communication is a well-known natural phenomenon that has been associated with sensitive responses to certain chemical stimuli, yet incorporation of intercellular communication into bioreporter design is exceedingly rare and the effect of intercellular signaling on single-cellular and population level responses has not been explicitly characterized before. In this work, a multicellular simulator implementing the Gillespie algorithm and compatible with Virtual Cell-designed networks is created and used to demonstrate nuances to multicellular stochastic simulations. The algorithm was used to simulate multiple cells in a reaction network in which a self-promoting and membrane permeable transcription factor also induces production of a cell-bound reporter protein. Simulations with 1, 10, 25, 50, 75, and 100 cells for 20 seconds of simulated time show that while final reporter protein count per cell decreased as cell count increased, aggregate final reporter protein number across all cells significantly increased. Interestingly, 50 second simulations show final reporter protein count per cell significantly increasing as the number of cells increases. The results show significant differences between the results of multicellular and single cellular simulations, which demonstrates the importance of simulating multiple cells to obtain nuanced results. A proof-of-concept bioreporter that responded to an environmental analyte while leveraging intercellular interactions for signal production was also designed and simulated for 50 seconds of simulation time in systems containing 1, 5, or 10 bioreporter cells. Greater number of bioreporter cells resulted in significantly greater amounts of signal protein produced in response to the same amount of starting analyte both on the population level and on the individual cellular level, indicating great potential for multicellular bioreporter designs to amplify response magnitude and sensitivity on the individual cellular level. These results warrant further research into the application of different simulation algorithms and multicellular bioreporter design and modeling.