An inverse modelling study on the local volume changes during early morphoelastic growth of the fetal human brain.

We take a data-driven approach to deducing the local volume changes accompanying early development of the fetal human brain. Our approach uses fetal brain atlas MRI data for the geometric changes in representative cases. Using a nonlinear continuum mechanics model of morphoelastic growth, we invert the deformation obtained from MRI registration to arrive at a field for the growth deformation gradient tensor. Our field inversion uses a combination of direct and adjoint methods for computing gradients of the objective function while constraining the optimization by the physics of morphoelastic growth. We thus infer a growth deformation gradient field that obeys the laws of morphoelastic growth. The errors between the MRI data and the forward displacement solution driven by the inverted growth deformation gradient field are found to be smaller than the reference displacement by well over an order of magnitude, and can be driven even lower. The results thus reproduce the three-dimensional growth during the early development of the fetal brain with controllable error. Our findings confirm that early growth is dominated by in plane cortical expansion rather than thickness increase.

Correction to: System inference for the spatio-temporal evolution of infectious diseases: Michigan in the time of COVID-19.

[This corrects the article DOI: 10.1007/s00466-020-01894-2.].

System inference for the spatio-temporal evolution of infectious diseases: Michigan in the time of COVID-19.

We extend the classical SIR model of infectious disease spread to account for time dependence in the parameters, which also include diffusivities. The temporal dependence accounts for the changing characteristics of testing, quarantine and treatment protocols, while diffusivity incorporates a mobile population. This model has been applied to data on the evolution of the COVID-19 pandemic in the US state of Michigan. For system inference, we use recent advances; specifically our framework for Variational System Identification (Wang et al. in Comput Methods Appl Mech Eng 356:44-74, 2019; arXiv:2001.04816 [cs.CE]) as well as Bayesian machine learning methods.

A mechanical model reveals that non-axisymmetric buckling lowers the energy barrier associated with membrane neck constriction.

Membrane neck formation is essential for scission, which, as recent experiments on tubules have demonstrated, can be location dependent. The diversity of biological machinery that can constrict a neck such as dynamin, actin, ESCRTs and BAR proteins, and the range of forces and deflection over which they operate, suggest that the constriction process is functionally mechanical and robust to changes in biological environment. In this study, we used a mechanical model of the lipid bilayer to systematically investigate the influence of location, symmetry constraints, and helical forces on membrane neck constriction. Simulations from our model demonstrated that the energy barriers associated with constriction of a membrane neck are location-dependent. Importantly, if symmetry restrictions are relaxed, then the energy barrier for constriction is dramatically lowered and the membrane buckles at lower values of forcing parameters. Our simulations also show that constriction due to helical proteins further reduces the energy barrier for neck formation when compared to cylindrical proteins. These studies establish that despite different molecular mechanisms of neck formation in cells, the mechanics of constriction naturally leads to a loss of symmetry that can lower the energy barrier to constriction.

A Diffuse Interface Framework for Modeling the Evolution of Multi-cell Aggregates as a Soft Packing Problem Driven by the Growth and Division of Cells.

We present a model for cell growth, division and packing under soft constraints that arise from the deformability of the cells as well as of a membrane that encloses them. Our treatment falls within the framework of diffuse interface methods, under which each cell is represented by a scalar phase field and the zero level set of the phase field represents the cell membrane. One crucial element in the treatment is the definition of a free energy density function that penalizes cell overlap, thus giving rise to a simple model of cell-cell contact. In order to properly represent cell packing and the associated free energy, we include a simplified representation of the anisotropic mechanical response of the underlying cytoskeleton and cell membrane through penalization of the cell shape change. Numerical examples demonstrate the evolution of multi-cell clusters and of the total free energy of the clusters as a consequence of growth, division and packing.

A computational study of stress fiber-focal adhesion dynamics governing cell contractility.

We apply a recently developed model of cytoskeletal force generation to study a cell's intrinsic contractility, as well as its response to external loading. The model is based on a nonequilibrium thermodynamic treatment of the mechanochemistry governing force in the stress fiber-focal adhesion system. Our computational study suggests that the mechanical coupling between the stress fibers and focal adhesions leads to a complex, dynamic, mechanochemical response. We collect the results in response maps whose regimes are distinguished by the initial geometry of the stress fiber-focal adhesion system, and by the external load on the cell. The results from our model connect qualitatively with recent studies on the force response of smooth muscle cells on arrays of polymeric microposts.

Perspectives on biological growth and remodeling.

The continuum mechanical treatment of biological growth and remodeling has attracted considerable attention over the past fifteen years. Many aspects of these problems are now well-understood, yet there remain areas in need of significant development from the standpoint of experiments, theory, and computation. In this perspective paper we review the state of the field and highlight open questions, challenges, and avenues for further development.

In silico estimates of the free energy rates in growing tumor spheroids.

The physics of solid tumor growth can be considered at three distinct size scales: the tumor scale, the cell-extracellular matrix (ECM) scale and the sub-cellular scale. In this paper we consider the tumor scale in the interest of eventually developing a system-level understanding of the progression of cancer. At this scale, cell populations and chemical species are best treated as concentration fields that vary with time and space. The cells have chemo-mechanical interactions with each other and with the ECM, consume glucose and oxygen that are transported through the tumor, and create chemical by-products. We present a continuum mathematical model for the biochemical dynamics and mechanics that govern tumor growth. The biochemical dynamics and mechanics also engender free energy changes that serve as universal measures for comparison of these processes. Within our mathematical framework we therefore consider the free energy inequality, which arises from the first and second laws of thermodynamics. With the model we compute preliminary estimates of the free energy rates of a growing tumor in its pre-vascular stage by using currently available data from single cells and multicellular tumor spheroids.

The micromechanics of fluid-solid interactions during growth in porous soft biological tissue.

In this paper, we address some modelling issues related to biological growth. Our treatment is based on a formulation for growth that was proposed within the context of mixture theory (J Mech Phys Solids 52:1595-1625, 2004). We aim to make this treatment more appropriate for the physics of porous soft tissues, paying particular attention to the nature of fluid transport, and mechanics of fluid and solid phases. The interactions between transport and mechanics have significant implications for growth and swelling. We also reformulate the governing differential equations for reaction-transport of solutes to represent the incompressibility constraint on the fluid phase of the tissue. This revision enables a straightforward implementation of numerical stabilisation for the advection-dominated limit of these equations. A finite element implementation with operator splitting is used to solve the coupled, non-linear partial differential equations that arise from the theory. We carry out a numerical and analytic study of the convergence of the operator splitting scheme subject to strain- and stress-homogenisation of the mechanics of fluid-solid interactions. A few computations are presented to demonstrate aspects of the physical mechanisms, and the numerical performance of the formulation.