Reverse engineering morphogenesis through Bayesian optimization of physics-based models.

Morphogenetic programs coordinate cell signaling and mechanical interactions to shape organs. In systems and synthetic biology, a key challenge is determining optimal cellular interactions for predicting organ shape, size, and function. Physics-based models defining the subcellular force distribution facilitate this, but it is challenging to calibrate parameters in these models from data. To solve this inverse problem, we created a Bayesian optimization framework to determine the optimal cellular force distribution such that the predicted organ shapes match the experimentally observed organ shapes. This integrative framework employs Gaussian Process Regression, a non-parametric kernel-based probabilistic machine learning modeling paradigm, to learn the mapping functions relating to the morphogenetic programs that maintain the final organ shape. We calibrated and tested the method on Drosophila wing imaginal discs to study mechanisms that regulate epithelial processes ranging from development to cancer. The parameter estimation framework successfully infers the underlying changes in core parameters needed to match simulation data with imaging data of wing discs perturbed with collagenase. The computational pipeline identifies distinct parameter sets mimicking wild-type shapes. It enables a global sensitivity analysis to support the regulation of actomyosin contractility and basal ECM stiffness to generate and maintain the curved shape of the wing imaginal disc. The optimization framework, combined with experimental imaging, identified that Piezo, a mechanosensitive ion channel, impacts fold formation by regulating the apical-basal balance of actomyosin contractility and elasticity of ECM. This workflow is extensible toward reverse-engineering morphogenesis across organ systems and for real-time control of complex multicellular systems.

Balancing competing effects of tissue growth and cytoskeletal regulation during Drosophila wing disc development.

How a developing organ robustly coordinates the cellular mechanics and growth to reach a final size and shape remains poorly understood. Through iterations between experiments and model simulations that include a mechanistic description of interkinetic nuclear migration, we show that the local curvature, height, and nuclear positioning of cells in the Drosophila wing imaginal disc are defined by the concurrent patterning of actomyosin contractility, cell-ECM adhesion, ECM stiffness, and interfacial membrane tension. We show that increasing cell proliferation via different growth-promoting pathways results in two distinct phenotypes. Triggering proliferation through insulin signaling increases basal curvature, but an increase in growth through Dpp signaling and Myc causes tissue flattening. These distinct phenotypic outcomes arise from differences in how each growth pathway regulates the cellular cytoskeleton, including contractility and cell-ECM adhesion. The coupled regulation of proliferation and cytoskeletal regulators is a general strategy to meet the multiple context-dependent criteria defining tissue morphogenesis.

Optimal Performance Objectives in the Highly Conserved Bone Morphogenetic Protein Signaling Pathway.

Throughout development, complex networks of cell signaling pathways drive cellular decision-making across different tissues and contexts. The transforming growth factor ß (TGF-ß) pathways, including the BMP/Smad pathway, play crucial roles in these cellular responses. However, as the Smad pathway is used reiteratively throughout the life cycle of all animals, its systems-level behavior varies from one context to another, despite the pathway connectivity remaining nearly constant. For instance, some cellular systems require a rapid response, while others require high noise filtering. In this paper, we examine how the BMP- Smad pathway balances trade-offs among three such systems-level behaviors, or "Performance Objectives (POs)": response speed, noise amplification, and the sensitivity of pathway output to receptor input. Using a Smad pathway model fit to human cell data, we show that varying non-conserved parameters (NCPs) such as protein concentrations, the Smad pathway can be tuned to emphasize any of the three POs and that the concentration of nuclear phosphatase has the greatest effect on tuning the POs. However, due to competition among the POs, the pathway cannot simultaneously optimize all three, but at best must balance trade-offs among the POs. We applied the multi-objective optimization concept of the Pareto Front, a widely used concept in economics to identify optimal trade-offs among various requirements. We show that the BMP pathway efficiently balances competing POs across species and is largely Pareto optimal. Our findings reveal that varying the concentration of NCPs allows the Smad signaling pathway to generate a diverse range of POs. This insight identifies how signaling pathways can be optimally tuned for each context.

Reverse engineering morphogenesis through Bayesian optimization of physics-based models.

Morphogenetic programs direct the cell signaling and nonlinear mechanical interactions between multiple cell types and tissue layers to define organ shape and size. A key challenge for systems and synthetic biology is determining optimal combinations of intra- and inter-cellular interactions to predict an organ's shape, size, and function. Physics-based mechanistic models that define the subcellular force distribution facilitate this, but it is extremely challenging to calibrate parameters in these models from data. To solve this inverse problem, we created a Bayesian optimization framework to determine the optimal cellular force distribution such that the predicted organ shapes match the desired organ shapes observed within the experimental imaging data. This integrative framework employs Gaussian Process Regression (GPR), a non-parametric kernel-based probabilistic machine learning modeling paradigm, to learn the mapping functions relating to the morphogenetic programs that generate and maintain the final organ shape. We calibrated and tested the method on cross-sections of Drosophila wing imaginal discs, a highly informative model organ system, to study mechanisms that regulate epithelial processes that range from development to cancer. As a specific test case, the parameter estimation framework successfully infers the underlying changes in core parameters needed to match simulation data with time series imaging data of wing discs perturbed with collagenase. Unexpectedly, the framework also identifies multiple distinct parameter sets that generate shapes similar to wild-type organ shapes. This platform enables an efficient, global sensitivity analysis to support the necessity of both actomyosin contractility and basal ECM stiffness to generate and maintain the curved shape of the wing imaginal disc. The optimization framework, combined with fixed tissue imaging, identified that Piezo, a mechanosensitive ion channel, impacts fold formation by regulating the apical-basal balance of actomyosin contractility and elasticity of ECM. This framework is extensible toward reverse-engineering the morphogenesis of any organ system and can be utilized in real-time control of complex multicellular systems.

Quantitative insights in tissue growth and morphogenesis with optogenetics.

Cells communicate with each other to jointly regulate cellular processes during cellular differentiation and tissue morphogenesis. This multiscale coordination arises through the spatiotemporal activity of morphogens to pattern cell signaling and transcriptional factor activity. This coded information controls cell mechanics, proliferation, and differentiation to shape the growth and morphogenesis of organs. While many of the molecular components and physical interactions have been identified in key model developmental systems, there are still many unresolved questions related to the dynamics involved due to challenges in precisely perturbing and quantitatively measuring signaling dynamics. Recently, a broad range of synthetic optogenetic tools have been developed and employed to quantitatively define relationships between signal transduction and downstream cellular responses. These optogenetic tools can control intracellular activities at the single cell or whole tissue scale to direct subsequent biological processes. In this brief review, we highlight a selected set of studies that develop and implement optogenetic tools to unravel quantitative biophysical mechanisms for tissue growth and morphogenesis across a broad range of biological systems through the manipulation of morphogens, signal transduction cascades, and cell mechanics. More generally, we discuss how optogenetic tools have emerged as a powerful platform for probing and controlling multicellular development.

Vestigial-dependent induction contributes to robust patterning but is not essential for wing-fate recruitment in Drosophila.

Cell recruitment is a process by which a differentiated cell induces neighboring cells to adopt its same cell fate. In Drosophila, cells expressing the protein encoded by the wing selector gene, vestigial (vg), drive a feed-forward recruitment signal that expands the Vg pattern as a wave front. However, previous studies on Vg pattern formation do not reveal these dynamics. Here, we use live imaging to show that multiple cells at the periphery of the wing disc simultaneously activate a fluorescent reporter of the recruitment signal, suggesting that cells may be recruited without the need for their contact neighbors be recruited in advance. In support of this observation, when Vg expression is inhibited either at the dorsal-ventral boundary or away from it, the activation of the recruitment signal still occurs at a distance, suggesting that Vg expression is not absolutely required to send or propagate the recruitment signal. However, the strength and extent of the recruitment signal is clearly compromised. We conclude that a feed-forward, contact-dependent cell recruitment process is not essential for Vg patterning, but it is necessary for robustness. Overall, our findings reveal a previously unidentified role of cell recruitment as a robustness-conferring cell differentiation mechanism.

The multimodal action of G alpha q in coordinating growth and homeostasis in the Drosophila wing imaginal disc.

Background: G proteins mediate cell responses to various ligands and play key roles in organ development. Dysregulation of G-proteins or Ca 2+ signaling impacts many human diseases and results in birth defects. However, the downstream effectors of specific G proteins in developmental regulatory networks are still poorly understood. Methods: We employed the Gal4/UAS binary system to inhibit or overexpress Gαq in the wing disc, followed by phenotypic analysis. Immunohistochemistry and next-gen RNA sequencing identified the downstream effectors and the signaling cascades affected by the disruption of Gαq homeostasis. Results: Here, we characterized how the G protein subunit Gαq tunes the size and shape of the wing in the larval and adult stages of development. Downregulation of Gαq in the wing disc reduced wing growth and delayed larval development. Gαq overexpression is sufficient to promote global Ca 2+ waves in the wing disc with a concomitant reduction in the Drosophila final wing size and a delay in pupariation. The reduced wing size phenotype is further enhanced when downregulating downstream components of the core Ca 2+ signaling toolkit, suggesting that downstream Ca 2+ signaling partially ameliorates the reduction in wing size. In contrast, Gαq -mediated pupariation delay is rescued by inhibition of IP 3 R, a key regulator of Ca 2+ signaling. This suggests that Gαq regulates developmental phenotypes through both Ca 2+ -dependent and Ca 2+ -independent mechanisms. RNA seq analysis shows that disruption of Gαq homeostasis affects nuclear hormone receptors, JAK/STAT pathway, and immune response genes. Notably, disruption of Gαq homeostasis increases expression levels of Dilp8, a key regulator of growth and pupariation timing. Conclusion: Gαq activity contributes to cell size regulation and wing metamorphosis. Disruption to Gαq homeostasis in the peripheral wing disc organ delays larval development through ecdysone signaling inhibition. Overall, Gαq signaling mediates key modules of organ size regulation and epithelial homeostasis through the dual action of Ca 2+ -dependent and independent mechanisms.

MAPPER: An Open-Source, High-Dimensional Image Analysis Pipeline Unmasks Differential Regulation of Drosophila Wing Features.

Phenomics requires quantification of large volumes of image data, necessitating high throughput image processing approaches. Existing image processing pipelines for Drosophila wings, a powerful genetic model for studying the underlying genetics for a broad range of cellular and developmental processes, are limited in speed, precision, and functional versatility. To expand on the utility of the wing as a phenotypic screening system, we developed MAPPER, an automated machine learning-based pipeline that quantifies high-dimensional phenotypic signatures, with each dimension quantifying a unique morphological feature of the Drosophila wing. MAPPER magnifies the power of Drosophila phenomics by rapidly quantifying subtle phenotypic differences in sample populations. We benchmarked MAPPER's accuracy and precision in replicating manual measurements to demonstrate its widespread utility. The morphological features extracted using MAPPER reveal variable sexual dimorphism across Drosophila species and unique underlying sex-specific differences in morphogen signaling in male and female wings. Moreover, the length of the proximal-distal axis across the species and sexes shows a conserved scaling relationship with respect to the wing size. In sum, MAPPER is an open-source tool for rapid, high-dimensional analysis of large imaging datasets. These high-content phenomic capabilities enable rigorous and systematic identification of genotype-to-phenotype relationships in a broad range of screening and drug testing applications and amplify the potential power of multimodal genomic approaches.

Principles for the design of multicellular engineered living systems.

Remarkable progress in bioengineering over the past two decades has enabled the formulation of fundamental design principles for a variety of medical and non-medical applications. These advancements have laid the foundation for building multicellular engineered living systems (M-CELS) from biological parts, forming functional modules integrated into living machines. These cognizant design principles for living systems encompass novel genetic circuit manipulation, self-assembly, cell-cell/matrix communication, and artificial tissues/organs enabled through systems biology, bioinformatics, computational biology, genetic engineering, and microfluidics. Here, we introduce design principles and a blueprint for forward production of robust and standardized M-CELS, which may undergo variable reiterations through the classic design-build-test-debug cycle. This Review provides practical and theoretical frameworks to forward-design, control, and optimize novel M-CELS. Potential applications include biopharmaceuticals, bioreactor factories, biofuels, environmental bioremediation, cellular computing, biohybrid digital technology, and experimental investigations into mechanisms of multicellular organisms normally hidden inside the "black box" of living cells.

Rational Design and Identification of Harmine-Inspired, N-Heterocyclic DYRK1A Inhibitors Employing a Functional Genomic In Vivo Drosophila Model System.

Deregulation of dual-specificity tyrosine phosphorylation-regulated kinase 1A (DYRK1A) plays a significant role in developmental brain defects, early-onset neurodegeneration, neuronal cell loss, dementia, and several types of cancer. Herein, we report the discovery of three new classes of N-heterocyclic DYRK1A inhibitors based on the potent, yet toxic kinase inhibitors, harmine and harmol. An initial in vitro evaluation of the small molecule library assembled revealed that the core heterocyclic motifs benzofuranones, oxindoles, and pyrrolones, showed statistically significant DYRK1A inhibition. Further, the utilization of a low cost, high-throughput functional genomic in vivo model system to identify small molecule inhibitors that normalize DYRK1A overexpression phenotypes is described. This in vivo assay substantiated the in vitro results, and the resulting correspondence validates generated classes as architectural motifs that serve as potential DYRK1A inhibitors. Further expansion and analysis of these core compound structures will allow discovery of safe, more effective chemical inhibitors of DYRK1A to ameliorate phenotypes caused by DYRK1A overexpression.

From spikes to intercellular waves: Tuning intercellular calcium signaling dynamics modulates organ size control.

Information flow within and between cells depends significantly on calcium (Ca2+) signaling dynamics. However, the biophysical mechanisms that govern emergent patterns of Ca2+ signaling dynamics at the organ level remain elusive. Recent experimental studies in developing Drosophila wing imaginal discs demonstrate the emergence of four distinct patterns of Ca2+ activity: Ca2+ spikes, intercellular Ca2+ transients, tissue-level Ca2+ waves, and a global "fluttering" state. Here, we used a combination of computational modeling and experimental approaches to identify two different populations of cells within tissues that are connected by gap junction proteins. We term these two subpopulations "initiator cells," defined by elevated levels of Phospholipase C (PLC) activity, and "standby cells," which exhibit baseline activity. We found that the type and strength of hormonal stimulation and extent of gap junctional communication jointly determine the predominate class of Ca2+ signaling activity. Further, single-cell Ca2+ spikes are stimulated by insulin, while intercellular Ca2+ waves depend on Gαq activity. Our computational model successfully reproduces how the dynamics of Ca2+ transients varies during organ growth. Phenotypic analysis of perturbations to Gαq and insulin signaling support an integrated model of cytoplasmic Ca2+ as a dynamic reporter of overall tissue growth. Further, we show that perturbations to Ca2+ signaling tune the final size of organs. This work provides a platform to further study how organ size regulation emerges from the crosstalk between biochemical growth signals and heterogeneous cell signaling states.

Pinching and pushing: fold formation in the Drosophila dorsal epidermis.

Epithelial folding is a fundamental morphogenetic process that shapes planar epithelial sheets into complex three-dimensional structures. Multiple mechanisms can generate epithelial folds, including apical constriction, which acts locally at the cellular level, differential growth on the tissue scale, or buckling because of compression from neighboring tissues. Here, we investigate the formation of dorsally located epithelial folds at segment boundaries during the late stages of Drosophila embryogenesis. We found that the fold formation at the segment boundaries occurs through the juxtaposition of two key morphogenetic processes: local apical constriction and tissue-level compressive forces from posterior segments. Further, we found that epidermal spreading and fold formation are accompanied by spatiotemporal pulses of Hedgehog (Hh) signaling. A computational model that incorporates the local forces generated from the differential tensions of the apical, basal, and lateral sides of the cell and active forces generated within the whole tissue recapitulates the overall fold formation process in wild-type and Hh overexpression conditions. In sum, this work demonstrates how epithelial folding depends on multiple, separable physical mechanisms to generate the final morphology of the dorsal epidermis. This work illustrates the modularity of morphogenetic unit operations that occur during epithelial morphogenesis.

Epithelial organ shape is generated by patterned actomyosin contractility and maintained by the extracellular matrix.

Epithelial sheets define organ architecture during development. Here, we employed an iterative multiscale computational modeling and quantitative experimental approach to decouple direct and indirect effects of actomyosin-generated forces, nuclear positioning, extracellular matrix, and cell-cell adhesion in shaping Drosophila wing imaginal discs. Basally generated actomyosin forces generate epithelial bending of the wing disc pouch. Surprisingly, acute pharmacological inhibition of ROCK-driven actomyosin contractility does not impact the maintenance of tissue height or curved shape. Computational simulations show that ECM tautness provides only a minor contribution to modulating tissue shape. Instead, passive ECM pre-strain serves to maintain the shape independent from actomyosin contractility. These results provide general insight into how the subcellular forces are generated and maintained within individual cells to induce tissue curvature. Thus, the results suggest an important design principle of separable contributions from ECM prestrain and actomyosin tension during epithelial organogenesis and homeostasis.

Rab11b-mediated integrin recycling promotes brain metastatic adaptation and outgrowth.

Breast cancer brain metastases (BCBM) have a 5-20 year latency and account for 30% of mortality; however, mechanisms governing adaptation to the brain microenvironment remain poorly defined. We combine time-course RNA-sequencing of BCBM development with a Drosophila melanogaster genetic screen, and identify Rab11b as a functional mediator of metastatic adaptation. Proteomic analysis reveals that Rab11b controls the cell surface proteome, recycling proteins required for successful interaction with the microenvironment, including integrin ß1. Rab11b-mediated control of integrin ß1 surface expression allows efficient engagement with the brain ECM, activating mechanotransduction signaling to promote survival. Lipophilic statins prevent membrane association and activity of Rab11b, and we provide proof-of principle that these drugs prevent breast cancer adaptation to the brain microenvironment. Our results identify Rab11b-mediated recycling of integrin ß1 as regulating BCBM, and suggest that the recycleome, recycling-based control of the cell surface proteome, is a previously unknown driver of metastatic adaptation and outgrowth.

Interplay between morphogen-directed positional information systems and physiological signaling.

The development of an organism from an undifferentiated single cell into a spatially complex structure requires spatial patterning of cell fates across tissues. Positional information, proposed by Lewis Wolpert in 1969, has led to the characterization of many components involved in regulating morphogen signaling activity. However, how morphogen gradients are established, maintained, and interpreted by cells still is not fully understood. Quantitative and systems-based approaches are increasingly needed to define general biological design rules that govern positional information systems in developing organisms. This short review highlights a selective set of studies that have investigated the roles of physiological signaling in modulating and mediating morphogen-based pattern formation. Similarities between neural transmission and morphogen-based pattern formation mechanisms suggest underlying shared principles of active cell-based communication. Within larger tissues, neural networks provide directed information, via physiological signaling, that supplements positional information through diffusion. Further, mounting evidence demonstrates that physiological signaling plays a role in ensuring robustness of morphogen-based signaling. We conclude by highlighting several outstanding questions regarding the role of physiological signaling in morphogen-based pattern formation. Elucidating how physiological signaling impacts positional information is critical for understanding the close coupling of developmental and cellular processes in the context of development, disease, and regeneration.

Rapid Fabrication of Custom Microfluidic Devices for Research and Educational Applications.

Microfluidic devices allow for the manipulation of fluids, particles, cells, micro-sized organs or organisms in channels ranging from the nano to submillimeter scales. A rapid increase in the use of this technology in the biological sciences has prompted a need for methods that are accessible to a wide range of research groups. Current fabrication standards, such as PDMS bonding, require expensive and time consuming lithographic and bonding techniques. A viable alternative is the use of equipment and materials that are easily affordable, require minimal expertise and allow for the rapid iteration of designs. In this work we describe a protocol for designing and producing PET-laminates (PETLs), microfluidic devices that are inexpensive, easy to fabricate, and consume significantly less time to generate than other approaches to microfluidics technology. They consist of thermally bonded film sheets, in which channels and other features are defined using a craft cutter. PETLs solve field-specific technical challenges while dramatically reducing obstacles to adoption. This approach facilitates the accessibility of microfluidics devices in both research and educational settings, providing a reliable platform for new methods of inquiry.

Multi-cellular engineered living systems: building a community around responsible research on emergence.

Ranging from miniaturized biological robots to organoids, multi-cellular engineered living systems (M-CELS) pose complex ethical and societal challenges. Some of these challenges, such as how to best distribute risks and benefits, are likely to arise in the development of any new technology. Other challenges arise specifically because of the particular characteristics of M-CELS. For example, as an engineered living system becomes increasingly complex, it may provoke societal debate about its moral considerability, perhaps necessitating protection from harm or recognition of positive moral and legal rights, particularly if derived from cells of human origin. The use of emergence-based principles in M-CELS development may also create unique challenges, making the technology difficult to fully control or predict in the laboratory as well as in applied medical or environmental settings. In response to these challenges, we argue that the M-CELS community has an obligation to systematically address the ethical and societal aspects of research and to seek input from and accountability to a broad range of stakeholders and publics. As a newly developing field, M-CELS has a significant opportunity to integrate ethically responsible norms and standards into its research and development practices from the start. With the aim of seizing this opportunity, we identify two general kinds of salient ethical issues arising from M-CELS research, and then present a set of commitments to and strategies for addressing these issues. If adopted, these commitments and strategies would help define M-CELS as not only an innovative field, but also as a model for responsible research and engineering.

Combined Scaffold Evaluation and Systems-Level Transcriptome-Based Analysis for Accelerated Lead Optimization Reveals Ribosomal Targeting Spirooxindole Cyclopropanes.

With evolutionary drug resistance impacting efforts to treat disease, the need for small molecules that exhibit novel molecular mechanisms of action is paramount. In this study, we combined scaffold-directed synthesis with a hybrid experimental and transcriptome analysis to identify bis-spirooxindole cyclopropanes that inhibit cancer cell proliferation through disruption of ribosomal function. These findings demonstrate the value of an integrated, biologically inspired synthesis and assay strategy for the accelerated identification of first-in-class cancer therapeutic candidates.

Tools to reverse-engineer multicellular systems: case studies using the fruit fly.

Reverse-engineering how complex multicellular systems develop and function is a grand challenge for systems bioengineers. This challenge has motivated the creation of a suite of bioengineering tools to develop increasingly quantitative descriptions of multicellular systems. Here, we survey a selection of these tools including microfluidic devices, imaging and computer vision techniques. We provide a selected overview of the emerging cross-talk between engineering methods and quantitative investigations within developmental biology. In particular, the review highlights selected recent examples from the Drosophila system, an excellent platform for understanding the interplay between genetics and biophysics. In sum, the integrative approaches that combine multiple advances in these fields are increasingly necessary to enable a deeper understanding of how to analyze both natural and synthetic multicellular systems.

Microfluidics on the fly: Inexpensive rapid fabrication of thermally laminated microfluidic devices for live imaging and multimodal perturbations of multicellular systems.

Microfluidic devices provide a platform for analyzing both natural and synthetic multicellular systems. Currently, substantial capital investment and expertise are required for creating microfluidic devices using standard soft-lithography. These requirements present barriers to entry for many nontraditional users of microfluidics, including developmental biology laboratories. Therefore, fabrication methodologies that enable rapid device iteration and work "out-of-the-box" can accelerate the integration of microfluidics with developmental biology. Here, we have created and characterized low-cost hybrid polyethylene terephthalate laminate (PETL) microfluidic devices that are suitable for cell and micro-organ culture assays. These devices were validated with mammalian cell lines and the Drosophila wing imaginal disc as a model micro-organ. First, we developed and tested PETLs that are compatible with both long-term cultures and high-resolution imaging of cells and organs. Further, we achieved spatiotemporal control of chemical gradients across the wing discs with a multilayered microfluidic device. Finally, we created a multilayered device that enables controllable mechanical loading of micro-organs. This mechanical actuation assay was used to characterize the response of larval wing discs at different developmental stages. Interestingly, increased deformation of the older wing discs for the same mechanical loading suggests that the compliance of the organ is increased in preparation for subsequent morphogenesis. Together, these results demonstrate the applicability of hybrid PETL devices for biochemical and mechanobiology studies on micro-organs and provide new insights into the mechanics of organ development.