Bayesian cell therapy process optimization.

Optimizing complex bioprocesses poses a significant challenge in several fields, particularly in cell therapy manufacturing. The development of customized, closed, and automated processes is crucial for their industrial translation and for addressing large patient populations at a sustainable price. Limited understanding of the underlying biological mechanisms, coupled with highly resource-intensive experimentation, are two contributing factors that make the development of these next-generation processes challenging. Bayesian optimization (BO) is an iterative experimental design methodology that addresses these challenges, but has not been extensively tested in situations that require parallel experimentation with significant experimental variability. In this study, we present an evaluation of noisy, parallel BO for increasing noise levels and parallel batch sizes on two in silico bioprocesses, and compare it to the industry state-of-the-art. As an in vitro showcase, we apply the method to the optimization of a monocyte purification unit operation. The in silico results show that BO significantly outperforms the state-of-the-art, requiring approximately 50% fewer experiments on average. This study highlights the potential of noisy, parallel BO as valuable tool for cell therapy process development and optimization.

Intercellular Adhesion Stiffness Moderates Cell Decoupling as a Function of Substrate Stiffness.

The interplay between cell-cell and cell-substrate interactions is complex yet necessary for the formation and healthy functioning of tissues. The same mechanosensing mechanisms used by the cell to sense its extracellular matrix also play a role in intercellular interactions. We used the discrete element method to develop a computational model of a deformable cell that includes subcellular components responsible for mechanosensing. We modeled a three-dimensional cell pair on a patterned (two-dimensional) substrate, a simple laboratory setup to study intercellular interactions. We explicitly modeled focal adhesions and adherens junctions. These mechanosensing adhesions matured, becoming stabilized by force. We also modeled contractile stress fibers that bind the discrete adhesions. The mechanosensing fibers strengthened upon stalling. Traction exerted on the substrate was used to generate traction maps (along the cell-substrate interface). These simulated maps are compared to experimental maps obtained via traction force microscopy. The model recreates the dependence on substrate stiffness of the tractions' spatial distribution, contractile moment of the cell pair, intercellular force, and number of focal adhesions. It also recreates the phenomenon of cell decoupling, in which cells exert forces separately when substrate stiffness increases. More importantly, the model provides viable molecular explanations for decoupling: mechanosensing mechanisms are responsible for competition between different fiber-adhesion configurations present in the cell pair. The point at which an increasing substrate stiffness becomes as high as that of the cell-cell interface is the tipping point at which configurations that favor cell-substrate adhesion dominate over those favoring cell-cell adhesion. This competition is responsible for decoupling.

Modeling of Mechanosensing Mechanisms Reveals Distinct Cell Migration Modes to Emerge From Combinations of Substrate Stiffness and Adhesion Receptor-Ligand Affinity.

Mesenchymal cell migration is an integral process in development and healing. The process is regulated by both mechanical and biochemical properties. Mechanical properties of the environment are sensed through mechanosensing, which consists of molecular responses mediated by mechanical signals. We developed a computational model of a deformable 3D cell on a flat substrate using discrete element modeling. The cell is polarized in a single direction and thus moves along the long axis of the substrate. By modeling discrete focal adhesions and stress fibers, we implement two mechanosensing mechanisms: focal adhesion stabilization by force and stress fiber strengthening upon contraction stalling. Two substrate-associated properties, substrate (ligand) stiffness and adhesion receptor-ligand affinity (in the form of focal adhesion disassembly rate), were varied for different model setups in which the mechanosensing mechanisms are set as active or inactive. Cell displacement, focal adhesion number, and cellular traction were quantified and tracked in time. We found that varying substrate stiffness (a mechanical property) and adhesion receptor-ligand affinity (a biochemical property) simultaneously dictate the mode in which cells migrate; cells either move in a smooth manner reminiscent of keratocytes or in a cyclical manner reminiscent of epithelial cells. Mechanosensing mechanisms are responsible for the range of conditions in which a cell adopts a particular migration mode. Stress fiber strengthening, specifically, is responsible for cyclical migration due to build-up of enough force to elicit rupture of focal adhesions and retraction of the cellular rear. Together, both mechanisms explain bimodal dependence of cell migration on substrate stiffness observed in the literature.

Matrix deformations around angiogenic sprouts correlate to sprout dynamics and suggest pulling activity.

Angiogenesis is the formation of new blood vessels from the pre-existing vasculature. It is essential for normal tissue growth and regeneration, and also plays a key role in many diseases [Carmeliet in Nat Med 9:653-660, 2003]. Cytoskeletal components have been shown to be important for angiogenic sprout initiation and maintenance [Kniazeva and Putnam in Am J Physiol 297:C179-C187, 2009] as well as endothelial cell shape control during invasion [Elliott et al. in Nat Cell Biol 17:137-147, 2015]. The exact nature of cytoskeleton-mediated forces for sprout initiation and progression, however, remains poorly understood. Questions on the importance of tip cell pulling versus stalk cell pushing are to a large extent unanswered, which among others has to do with the difficulty of quantifying and resolving those forces in time and space. We developed methods based on time-lapse confocal microscopy and image processing-further termed 4D displacement microscopy-to acquire detailed, spatially and temporally resolved extracellular matrix (ECM) deformations, indicative of cell-ECM mechanical interactions around invading sprouts. We demonstrate that matrix deformations dependent on actin-mediated force generation are spatio-temporally correlated with sprout morphological dynamics. Furthermore, sprout tips were found to exert radially pulling forces on the extracellular matrix, which were quantified by means of a computational model of collagen ECM mechanics. Protrusions from extending sprouts mostly increase their pulling forces, while retracting protrusions mainly reduce their pulling forces. Displacement microscopy analysis further unveiled a characteristic dipole-like deformation pattern along the sprout direction that was consistent among seemingly very different sprout shapes-with oppositely oriented displacements at sprout tip versus sprout base and a transition zone of negligible displacements in between. These results demonstrate that sprout-ECM interactions are dominated by pulling forces and underline the key role of tip cell pulling for sprouting angiogenesis.

The role of actin protrusion dynamics in cell migration through a degradable viscoelastic extracellular matrix: Insights from a computational model.

Actin protrusion dynamics plays an important role in the regulation of three-dimensional (3D) cell migration. Cells form protrusions that adhere to the surrounding extracellular matrix (ECM), mechanically probe the ECM and contract in order to displace the cell body. This results in cell migration that can be directed by the mechanical anisotropy of the ECM. However, the subcellular processes that regulate protrusion dynamics in 3D cell migration are difficult to investigate experimentally and therefore not well understood. Here, we present a computational model of cell migration through a degradable viscoelastic ECM. This model is a 2D representation of 3D cell migration. The cell is modeled as an active deformable object that captures the viscoelastic behavior of the actin cortex and the subcellular processes underlying 3D cell migration. The ECM is regarded as a viscoelastic material, with or without anisotropy due to fibrillar strain stiffening, and modeled by means of the meshless Lagrangian smoothed particle hydrodynamics (SPH) method. ECM degradation is captured by local fluidization of the material and permits cell migration through the ECM. We demonstrate that changes in ECM stiffness and cell strength affect cell migration and are accompanied by changes in number, lifetime and length of protrusions. Interestingly, directly changing the total protrusion number or the average lifetime or length of protrusions does not affect cell migration. A stochastic variability in protrusion lifetime proves to be enough to explain differences in cell migration velocity. Force-dependent adhesion disassembly does not result in faster migration, but can make migration more efficient. We also demonstrate that when a number of simultaneous protrusions is enforced, the optimal number of simultaneous protrusions is one or two, depending on ECM anisotropy. Together, the model provides non-trivial new insights in the role of protrusions in 3D cell migration and can be a valuable contribution to increase the understanding of 3D cell migration mechanics.