A hyperelastic model for simulating cells in flow.
Ontology highlight
ABSTRACT:
SUBMITTER: Muller SJ
PROVIDER: S-EPMC7979664 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Biological and synthetic materials often exhibit intrinsic variability in their elastic responses under large strains, owing to microstructural inhomogeneity or when elastic data are extracted from viscoelastic mechanical tests. For these materials, although hyperelastic models calibrated to mean data are useful, stochastic representations accounting also for data dispersion carry extra information about the variability of material properties found in practical applications. We combine finite elasticity and information theories to construct homogeneous isotropic hyperelastic models with random field parameters calibrated to discrete mean values and standard deviations of either the stress-strain function or the nonlinear shear modulus, which is a function of the deformation, estimated from experimental tests. These quantities can take on different values, corresponding to possible outcomes of the experiments. As multiple models can be derived that adequately represent the observed phenomena, we apply Occam's razor by providing an explicit criterion for model selection based on Bayesian statistics. We then employ this criterion to select a model among competing models calibrated to experimental data for rubber and brain tissue under single or multiaxial loads.
| S-EPMC5897763 | biostudies-literature
Project description:Shear stress induces directed endothelial cell (EC) migration in blood vessels leading to vessel diameter increase and induction of vascular maturation. Other factors, such as EC elongation and interaction between ECs and non-vascular areas are also important. Computational models have previously been used to study collective cell migration. These models can be used to predict EC migration and its effect on vascular remodelling during embryogenesis. We combined live time-lapse imaging of the remodelling vasculature of the quail embryo yolk sac with flow quantification using a combination of micro-Particle Image Velocimetry and computational fluid dynamics. We then used the flow and remodelling data to inform a model of EC migration during remodelling. To obtain the relation between shear stress and velocity in vitro for EC cells, we developed a flow chamber to assess how confluent sheets of ECs migrate in response to shear stress. Using these data as an input, we developed a multiphase, self-propelled particles (SPP) model where individual agents are driven to migrate based on the level of shear stress while maintaining appropriate spatial relationship to nearby agents. These agents elongate, interact with each other, and with avascular agents at each time-step of the model. We compared predicted vascular shape to real vascular shape after 4 hours from our time-lapse movies and performed sensitivity analysis on the various model parameters. Our model shows that shear stress has the largest effect on the remodelling process. Importantly, however, elongation played an especially important part in remodelling. This model provides a powerful tool to study the input of different biological processes on remodelling.
| S-EPMC7478591 | biostudies-literature
Project description:A gray lattice Boltzmann model has previously been developed by the authors of this article to simulate fluid flow in porous media that contain both resolved pores and grains as well as aggregates of unresolved smaller pores and grains. In this model, a single parameter is introduced to prescribe the amount of fluid to be bounced back at each aggregate cell. This model has been shown to recover Darcy-Brinkman flow but with effective viscosity and permeability correlated through the model parameter. In this paper, we prove that the model parameter relates to the fraction of the solid phase of a sub-pore system for a specific set of bounce-back conditions. We introduce an additional parameter to the model, and this enables flow simulation in which cases with variable effective viscosity and permeability can be specified by selecting the two parameters independently. We verify and validate the model for layered channel cases and mathematically analyze fluid momentum and energy losses for the single- and two-parameter models to explain the roles of the parameters in their conservation. We introduce a strategy to upgrade our model to an isotropic version. We discuss the fundamental differences between our model and the Brinkman body-force LBM scheme.
| S-EPMC5895624 | biostudies-literature
Project description:As drug development is extremely expensive, the identification of novel indications for in-market drugs is financially attractive. Multiple algorithms are used to support such drug repurposing, but highly reliable methods combining simulation of intracellular networks and machine learning are currently not available. We developed an algorithm that simulates drug effects on the flow of information through protein-protein interaction networks, and used support vector machine to identify potentially effective drugs in our model disease, psoriasis. Using this method, we screened about 1,500 marketed and investigational substances, identified 51 drugs that were potentially effective, and selected three of them for experimental confirmation. All drugs inhibited tumor necrosis factor alpha-induced nuclear factor kappa B activity in vitro, suggesting they might be effective for treating psoriasis in humans. Additionally, these drugs significantly inhibited imiquimod-induced ear thickening and inflammation in the mouse model of the disease. All results suggest high prediction performance for the algorithm.
| S-EPMC5836852 | biostudies-literature
Project description:Given the importance of hyperelastic constitutive models in the design of engineering components, researchers have been developing the improved and new constitutive models in search of a more accurate and even universal performance. Here, a modified hyperelastic constitutive model based on the Yeoh model is proposed to improve its prediction performance for multiaxial deformation of hyperelastic polymeric materials while retaining the advantages of the original Yeoh model. The modified constitutive model has one more correction term than the original model. The specific form of the correction term is a composite function based on a power function represented by the principal stretches, which is derived from the corresponding residual strain energy when the Yeoh model predicts the equibiaxial mode of deformation. In addition, a parameter identification method based on the cyclic genetic-pattern search algorithm is introduced to accurately obtain the parameters of the constitutive model. By applying the modified model to the experimental datasets of various rubber or rubber-like materials (including natural unfilled or filled rubber, silicone rubber, extremely soft hydrogel and human brain cortex tissue), it is confirmed that the modified model not only possesses a significantly improved ability to predict multiaxial deformation, but also has a wider range of material applicability. Meanwhile, the advantages of the modified model over most existing models in the literatures are also demonstrated. For example, when characterizing human brain tissue, which is difficult for most existing models in the literature, the modified model has comparable predictive accuracy with the third-order Ogden model, while maintaining convexity in the corresponding deformation domain. Moreover, the effective prediction ability of the modified model for untested equi-biaxial deformation of different materials has also been confirmed using only the data of uniaxial tension and pure shear from various datasets. The effective prediction for the untested equibiaxial deformation makes it more suitable for the practice situation where the equibiaxial deformation of certain polymeric materials is unavailable. Finally, compared with other parameter identification methods, the introduced parameter identification method significantly improves the predicted accuracy of the constitutive models; meanwhile, the uniform convergence of introduced parameter identification method is also better.
| S-EPMC10421088 | biostudies-literature
Project description:PurposeEstimation of brain deformation is crucial during neurosurgery. Whilst mechanical characterisation captures stress-strain relationships of tissue, biomechanical models are limited by experimental conditions. This results in variability reported in the literature. The aim of this work was to demonstrate a generative model of strain energy density functions can estimate the elastic properties of tissue using observed brain deformation.MethodsFor the generative model a Gaussian Process regression learns elastic potentials from 73 manuscripts. We evaluate the use of neo-Hookean, Mooney-Rivlin and 1-term Ogden meta-models to guarantee stability. Single and multiple tissue experiments validate the ability of our generative model to estimate tissue properties on a synthetic brain model and in eight temporal lobe resection cases where deformation is observed between pre- and post-operative images.ResultsEstimated parameters on a synthetic model are close to the known reference with a root-mean-square error (RMSE) of 0.1 mm and 0.2 mm between surface nodes for single and multiple tissue experiments. In clinical cases, we were able to recover brain deformation from pre- to post-operative images reducing RMSE of differences from 1.37 to 1.08 mm on the ventricle surface and from 5.89 to 4.84 mm on the resection cavity surface.ConclusionOur generative model can capture uncertainties related to mechanical characterisation of tissue. When fitting samples from elastography and linear studies, all meta-models performed similarly. The Ogden meta-model performed the best on hyperelastic studies. We were able to predict elastic parameters in a reference model on a synthetic phantom. However, deformation observed in clinical cases is only partly explained using our generative model.
| S-EPMC7822772 | biostudies-literature
Project description:Evidence suggests that both the interaction of so-called Merkel cells and the epidermal stress distribution play an important role in the formation of fingerprint patterns during pregnancy. To model the formation of fingerprint patterns in a biologically meaningful way these patterns have to become stationary. For the creation of synthetic fingerprints it is also very desirable that rescaling the model parameters leads to rescaled distances between the stationary fingerprint ridges. Based on these observations, as well as the model introduced by Kücken and Champod we propose a new model for the formation of fingerprint patterns during pregnancy. In this anisotropic interaction model the interaction forces not only depend on the distance vector between the cells and the model parameters, but additionally on an underlying tensor field, representing a stress field. This dependence on the tensor field leads to complex, anisotropic patterns. We study the resulting stationary patterns both analytically and numerically. In particular, we show that fingerprint patterns can be modeled as stationary solutions by choosing the underlying tensor field appropriately.
| S-EPMC6534530 | biostudies-literature
Project description:We present a multiscale model and an adaptive numerical scheme for simulating cardiac action potential propagation along a linear strand of heart muscle cells. This model couples macroscale partial differential equations posed over the tissue to microscale equations posed over discrete cellular geometry. The microscopic equations are used only near action potential wave fronts, and the macroscopic equations are used everywhere else. We study the effects of gap-junctional and ephaptic coupling on conduction in the multiscale model and its fully macroscale and fully microscale analogues. Our simulations reveal that the adaptive multiscale model accurately reproduces the action potential wave forms and wave speeds of the fully microscale model. They also demonstrate that, at low gap-junctional conductivities, the accuracy of fully macroscale simulations is sensitive to numerical grid spacing. Moreover, adaptive multiscale simulations capture the effect of ephaptic coupling, whereas fully macroscale simulations do not. We propose two ways of generalizing our multiscale model to higher dimensions, and we argue that such generalizations may be necessary to obtain accurate three-dimensional simulations of cardiac conduction in certain pathophysiological parameter regimes.
| S-EPMC2930419 | biostudies-literature
Project description:In many developmental and pathological processes, including cellular migration during normal development and invasion in cancer metastasis, cells are required to withstand severe deformations. The structural integrity of eukaryotic cells under small deformations has been known to depend on the cytoskeleton including actin filaments (F-actin), microtubules (MT), and intermediate filaments (IFs). However, it remains unclear how cells resist severe deformations since both F-actin and microtubules yield or disassemble under moderate strains. Using vimentin containing IFs (VIFs) as a model for studying the large family of IF proteins, we demonstrate that they dominate cytoplasmic mechanics and maintain cell viability at large deformations. Our results show that cytoskeletal VIFs form a stretchable, hyperelastic network in living cells. This network works synergistically with other cytoplasmic components, substantially enhancing the strength, stretchability, resilience, and toughness of cells. Moreover, we find the hyperelastic VIF network, together with other quickly recoverable cytoskeletal components, forms a mechanically robust structure which can mechanically recover after damage.
| S-EPMC6717279 | biostudies-literature
Project description:The Ovine spine is an accepted model to investigate the biomechanical behaviour of the human lumbar one. Indeed, the use of animal models for in vitro studies is necessary to investigate the mechanical behaviour of biological tissue, but needs to be reduced for ethical and social reasons. The aim of this study was to create a finite element model of the lumbar intervertebral disc of the sheep that may help to refine the understanding of parallel in vitro experiments and that can be used to predict when mechanical failure occurs. Anisotropic hyperelastic material properties were assigned to the annulus fibrosus and factorial optimization analyses were performed to find out the optimal parameters of the ground substance and of the collagen fibers. For the ground substance of the annulus fibrosus the investigation was based on experimental data taken from the literature, while for the collagen fibers tensile tests on annulus specimens were conducted. Flexibility analysis in flexion-extension, lateral bending and axial rotation were conducted. Different material properties for the anterior, lateral and posterior regions of the annulus were found. The posterior part resulted the stiffest region in compression whereas the anterior one the stiffest region in tension. Since the flexibility outcomes were in a good agreement with the literature data, we considered this model suitable to be used in conjunction with in vitro and in vivo tests to investigate the mechanical behaviour of the ovine lumbar disc.
| S-EPMC5417645 | biostudies-literature