Dataset Information

Smeared Multiscale Finite Element Models for Mass Transport and Electrophysiology Coupled to Muscle Mechanics.

ABSTRACT: Mass transport represents the most fundamental process in living organisms. It includes delivery of nutrients, oxygen, drugs, and other substances from the vascular system to tissue and transport of waste and other products from cells back to vascular and lymphatic network and organs. Furthermore, movement is achieved by mechanical forces generated by muscles in coordination with the nervous system. The signals coming from the brain, which have the character of electrical waves, produce activation within muscle cells. Therefore, from a physics perspective, there exist a number of physical fields within the body, such as velocities of transport, pressures, concentrations of substances, and electrical potential, which is directly coupled to biochemical processes of transforming the chemical into mechanical energy and further internal forces for motion. The overall problems of mass transport and electrophysiology coupled to mechanics can be investigated theoretically by developing appropriate computational models. Due to the enormous complexity of the biological system, it would be almost impossible to establish a detailed computational model for the physical fields related to mass transport, electrophysiology, and coupled fields. To make computational models feasible for applications, we here summarize a concept of smeared physical fields, with coupling among them, and muscle mechanics, which includes dependence on the electrical potential. Accuracy of the smeared computational models, also with coupling to muscle mechanics, is illustrated with simple example, while their applicability is demonstrated on a liver model with tumors present. The last example shows that the introduced methodology is applicable to large biological systems.

SUBMITTER: Kojic M

PROVIDER: S-EPMC6914730 | biostudies-literature | 2019

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Smeared Multiscale Finite Element Models for Mass Transport and Electrophysiology Coupled to Muscle Mechanics.

Kojic Milos M Milosevic Miljan M Simic Vladimir V Milicevic Bogdan B Geroski Vladimir V Nizzero Sara S Ziemys Arturas A Filipovic Nenad N Ferrari Mauro M

Frontiers in bioengineering and biotechnology 20191210

Mass transport represents the most fundamental process in living organisms. It includes delivery of nutrients, oxygen, drugs, and other substances from the vascular system to tissue and transport of waste and other products from cells back to vascular and lymphatic network and organs. Furthermore, movement is achieved by mechanical forces generated by muscles in coordination with the nervous system. The signals coming from the brain, which have the character of electrical waves, produce activati ...[more]

PMID: 31921800

Similar Datasets

Project description:BACKGROUND: Biological mass transport processes determine the behavior and function of cells, regulate interactions between synthetic agents and recipient targets, and are key elements in the design and use of biosensors. Accurately predicting the outcomes of such processes is crucial to both enhancing our understanding of how these systems function, enabling the design of effective strategies to control their function, and verifying that engineered solutions perform according to plan. METHODS: A Galerkin-based finite element model was developed and implemented to solve a system of two coupled partial differential equations governing biomolecule transport and reaction in live cells. The simulator was coupled, in the framework of an inverse modeling strategy, with an optimization algorithm and an experimental time series, obtained by the Fluorescence Recovery after Photobleaching (FRAP) technique, to estimate biomolecule mass transport and reaction rate parameters. In the inverse algorithm, an adaptive method was implemented to calculate sensitivity matrix. A multi-criteria termination rule was developed to stop the inverse code at the solution. The applicability of the model was illustrated by simulating the mobility and binding of GFP-tagged glucocorticoid receptor in the nucleoplasm of mouse adenocarcinoma. RESULTS: The numerical simulator shows excellent agreement with the analytic solutions and experimental FRAP data. Detailed residual analysis indicates that residuals have zero mean and constant variance and are normally distributed and uncorrelated. Therefore, the necessary and sufficient criteria for least square parameter optimization, which was used in this study, were met. CONCLUSION: The developed strategy is an efficient approach to extract as much physiochemical information from the FRAP protocol as possible. Well-posedness analysis of the inverse problem, however, indicates that the FRAP protocol provides insufficient information for unique simultaneous estimation of diffusion coefficient and binding rate parameters. Care should be exercised in drawing inferences, from FRAP data, regarding concentrations of free and bound proteins, average binding and diffusion times, and protein mobility unless they are confirmed by long-range Markov Chain-Monte Carlo (MCMC) methods and experimental observations.

Project description:This work presents a few variational multiscale models for charge transport in complex physical, chemical and biological systems and engineering devices, such as fuel cells, solar cells, battery cells, nanofluidics, transistors and ion channels. An essential ingredient of the present models, introduced in an earlier paper (Bulletin of Mathematical Biology, 72, 1562-1622, 2010), is the use of differential geometry theory of surfaces as a natural means to geometrically separate the macroscopic domain from the microscopic domain, meanwhile, dynamically couple discrete and continuum descriptions. Our main strategy is to construct the total energy functional of a charge transport system to encompass the polar and nonpolar free energies of solvation, and chemical potential related energy. By using the Euler-Lagrange variation, coupled Laplace-Beltrami and Poisson-Nernst-Planck (LB-PNP) equations are derived. The solution of the LB-PNP equations leads to the minimization of the total free energy, and explicit profiles of electrostatic potential and densities of charge species. To further reduce the computational complexity, the Boltzmann distribution obtained from the Poisson-Boltzmann (PB) equation is utilized to represent the densities of certain charge species so as to avoid the computationally expensive solution of some Nernst-Planck (NP) equations. Consequently, the coupled Laplace-Beltrami and Poisson-Boltzmann-Nernst-Planck (LB-PBNP) equations are proposed for charge transport in heterogeneous systems. A major emphasis of the present formulation is the consistency between equilibrium LB-PB theory and non-equilibrium LB-PNP theory at equilibrium. Another major emphasis is the capability of the reduced LB-PBNP model to fully recover the prediction of the LB-PNP model at non-equilibrium settings. To account for the fluid impact on the charge transport, we derive coupled Laplace-Beltrami, Poisson-Nernst-Planck and Navier-Stokes equations from the variational principle for chemo-electro-fluid systems. A number of computational algorithms is developed to implement the proposed new variational multiscale models in an efficient manner. A set of ten protein molecules and a realistic ion channel, Gramicidin A, are employed to confirm the consistency and verify the capability. Extensive numerical experiment is designed to validate the proposed variational multiscale models. A good quantitative agreement between our model prediction and the experimental measurement of current-voltage curves is observed for the Gramicidin A channel transport. This paper also provides a brief review of the field.

Project description:Over the past two decades, the increase in prevalence of cardiovascular diseases and the limited availability of autologous blood vessels and saphenous vein grafts have motivated the development of tissue-engineered vascular grafts (TEVGs). However, compliance mismatch and poor mechanical properties of the TEVGs remain as two major issues that need to be addressed. Researchers have investigated the role of various culture conditions and mechanical conditioning in deposition and orientation of collagen fibers, which are the key structural components in the vascular wall; however, the intrinsic complexity of mechanobiological interactions demands implementing new engineering approaches that allow researchers to investigate various scenarios more efficiently. In this study, we utilized a coupled agent-based finite element analysis (AB-FEA) modeling approach to study the effect of various loading modes (uniaxial, biaxial, and equibiaxial), boundary conditions, stretch magnitudes, and growth factor concentrations on growth and remodeling of smooth muscle cell-populated TEVGs, with specific focus on collagen deposition and orientation. Our simulations (12 weeks of culture) showed that biaxial cyclic loading (and not uniaxial or equibiaxial) leads to alignment of collagen fibers in the physiological directions. Moreover, axial boundary conditions of the TEVG act as determinants of fiber orientations. Decreasing the serum concentration, from 10% to 5% or 1%, significantly decreased the growth and remodeling speed, but only affected the fiber orientation in the 1% serum case. In conclusion, in silico tissue engineering has the potential to evolve the future of tissue engineering, as it will allow researchers to conceptualize various interactions and investigate numerous scenarios with great speed. In this study, we were able to predict the orientation of collagen fibers in TEVGs using a coupled AB-FEA model in less than 8?h. Impact Statement Tissue-engineered vascular grafts (TEVGs) hold potential to replace the current gold standard of vascular grafting, saphenous vein grafts. However, developing TEVGs that mimic the mechanical performance of the native tissue remains a challenging task. We developed a computational model of the grafts' remodeling processes and studied the effects of various loading mechanisms and culture conditions on collagen fiber orientation, which is a key factor in mechanical performance of the grafts. We were able to predict the fiber orientations accurately and show that biaxial loading and axial boundary conditions are important factors in collagen fiber organization.

Dataset Information

Smeared Multiscale Finite Element Models for Mass Transport and Electrophysiology Coupled to Muscle Mechanics.

Publications

Smeared Multiscale Finite Element Models for Mass Transport and Electrophysiology Coupled to Muscle Mechanics.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets