  The dynamics of cellular pattern formation is crucial for understanding
embryonic development and tissue morphogenesis. Recent studies have shown that
human dermal fibroblasts cultured on liquid crystal elastomers can exhibit an
increase in orientational alignment over time, accompanied by cell
proliferation, under the influence of the weak guidance of a molecularly
aligned substrate. However, a comprehensive understanding of how this order
arises remains largely unknown. This knowledge gap may be attributed, in part,
to a scarcity of mechanistic models that can capture the temporal progression
of the complex nonequilibrium dynamics during the cellular alignment process.
The orientational alignment occurs primarily when cells reach a high density
near confluence. Therefore, for accurate modeling, it is crucial to take into
account both the cell-cell interaction term and the influence from the
substrate, acting as a one-body external potential term. To fill in this gap,
we develop a hybrid procedure that utilizes statistical learning approaches to
extend the state-of-the-art physics models for quantifying both effects. We
develop a more efficient way to perform feature selection that avoids testing
all feature combinations through simulation. The maximum likelihood estimator
of the model was derived and implemented in computationally scalable algorithms
for model calibration and simulation. By including these features, such as the
non-Gaussian, anisotropic fluctuations, and limiting alignment interaction only
to neighboring cells with the same velocity direction, this model
quantitatively reproduce the key system-level parameters--the temporal
progression of the velocity orientational order parameters and the variability
of velocity vectors, whereas models missing any of the features fail to capture
these temporally dependent parameters.
