in-silico-clinical-trials-for-cardiovascular-disease

<em>In Silico</em> Clinical Trials for Cardiovascular Disease

The SILICOFCM project mainly aims to develop a computational platform for in silico clinical trials of familial cardiomyopathies (FCMs). The unique ...

Bioengineering Research and Development Center BioIRC

Serbian Academy of Science and Arts

Circulating tumor cells (CTCs) are known to be a harbinger of cancer metastasis. The CTCs are known to circulate as individual cells or as a group of interconnected cells called CTC clusters. Since both single CTCs and CTC clusters have been detected in venous blood samples of cancer patients, they needed to traverse at least one capillary bed when crossing from arterial to venous circulation. The diameter of a typical capillary is about 7µm, whereas the size of an individual CTC or CTC clusters can be greater than 20µm and thus size exclusion is believed to be an important factor in the capillary arrest of CTCs - a key early event in metastasis. To examine the biophysical conditions needed for capillary arrest, we have developed a custom-built viscoelastic solid-fluid 3D computational model that enables us to calculate, under physiological conditions, the maximal CTC diameter that will pass through the capillary. We show that large CTCs and CTC clusters can successfully cross capillaries if their stiffness is relatively small. Specifically, under physiological conditions, a 13µm diameter CTC passes through a 7µm capillary only if its stiffness is less than 500Pa and conversely, for a stiffness of 10Pa the maximal passing diameter can be as high as 140µm, such as for a cluster of CTCs. By exploring the parameter space, a relationship between the capillary blood pressure gradient and the CTC mechanical properties (size and stiffness) was determined. The presented computational platform and the resulting pressure-size-stiffness relationship can be employed as a tool to help study the biomechanical conditions needed for capillary arrest of CTCs and CTC clusters, provide predictive capabilities in disease progression based on biophysical CTC parameters, and aid in the rational design of size-based CTC isolation technologies where CTCs can experience large deformations due to high pressure gradients.

A computational study of circulating large tumor cells traversing microvessels.

We have recently introduced a composite smeared finite element (CSFE) to model gradient-driven mass transport in biological tissue. The transport from capillary system is smeared in a way to transform 1D transport to a continuum, while the tissue is considered as a continuum. Coupling between the smeared pressure and concentration field is achieved through 1D connectivity elements assigned at each FE node. Here we extend our smeared model to include the lymphatic system. The lymphatic vessels are treated in a way analogous to the capillaries, by introducing the corresponding Darcy and diffusion tensors. New connectivity elements are added. In the numerical examples we demonstrate accuracy of the smeared model and the effects of the lymph on the pressure and concentration within extracellular space are evaluated, assuming that there is no transport to the cell space.

Extension of the composite smeared finite element (CSFE) to include lymphatic system in modeling mass transport in capillary systems and biological tissue.

One of the basic and vital processes in living organisms is mass exchange, which occurs on several levels: it goes from blood vessels to cells and organelles within cells. On that path, molecules, as oxygen, metabolic products, drugs, etc. Traverse different macro and micro environments - blood, extracellular/intracellular space, and interior of organelles; and also biological barriers such as walls of blood vessels and membranes of cells and organelles. Many aspects of this mass transport remain unknown, particularly the biophysical mechanisms governing drug delivery. The main research approach relies on laboratory and clinical investigations. In parallel, considerable efforts have been directed to develop computational tools for additional insight into the intricate process of mass exchange and transport. Along these lines, we have recently formulated a composite smeared finite element (CSFE) which is composed of the smeared continuum pressure and concentration fields of the capillary and lymphatic system, and of these fields within tissue. The element offers an elegant and simple procedure which opens up new lines of inquiry and can be applied to large systems such as organs and tumors models. Here, we extend this concept to a multiscale scheme which concurrently couples domains that span from large blood vessels, capillaries and lymph, to cell cytosol and further to organelles of nanometer size. These spatial physical domains are coupled by the appropriate connectivity elements representing biological barriers. The composite finite element has "degrees of freedom" which include pressures and concentrations of all compartments of the vessels-tissue assemblage. The overall model uses the standard, measurable material properties of the continuum biological environments and biological barriers. It can be considered as a framework into which we can incorporate various additional effects (such as electrical or biochemical) for transport through membranes or within cells. This concept and the developed FE software within our package PAK offers a computational tool that can be applied to whole-organ systems, while also including specific domains such as tumors. The solved examples demonstrate the accuracy of this model and its applicability to large biological systems.

Multiscale smeared finite element model for mass transport in biological tissue: From blood vessels to cells and cellular organelles.

Due to the relative ease of producing nanofibers with a core⁻shell structure, emulsion electrospinning has been investigated intensively in making nanofibrous drug delivery systems for controlled and sustained release. Predictions of drug release rates from the poly (d,l-lactic-co-glycolic acid) (PLGA) produced via emulsion electrospinning can be a very difficult task due to the complexity of the system. A computational finite element methodology was used to calculate the diffusion mass transport of Rhodamine B (fluorescent drug model). Degradation effects and hydrophobicity (partitioning phenomenon) at the fiber/surrounding interface were included in the models. The results are validated by experiments where electrospun PLGA nanofiber mats with different contents were used. A new approach to three-dimensional (3D) modeling of nanofibers is presented in this work. The authors have introduced two original models for diffusive drug release from nanofibers to the 3D surrounding medium discretized by continuum 3D finite elements: (1) A model with simple radial one-dimensional (1D) finite elements, and (2) a model consisting of composite smeared finite elements (CSFEs). Numerical solutions, compared to experiments, demonstrate that both computational models provide accurate predictions of the diffusion process and can therefore serve as efficient tools for describing transport inside a polymer fiber network and drug release to the surrounding porous medium.

A Computational Model for Drug Release from PLGA Implant.

We couple a tumor growth model embedded in a microenvironment, with a bio distribution model able to simulate a whole organ. The growth model yields the evolution of tumor cell population, of the differential pressure between cell populations, of porosity of ECM, of consumption of nutrients due to tumor growth, of angiogenesis, and related growth factors as function of the locally available nutrient. The bio distribution model on the other hand operates on a frozen geometry but yields a much refined distribution of nutrient and other molecules. The combination of both models will enable simulating the growth of a tumor in a whole organ, including a realistic distribution of therapeutic agents and allow hence to evaluate the efficacy of these agents.

Coupling tumor growth and bio distribution models.

Basic functions of living organisms are governed by the nervous system through bidirectional signals transmitted from the brain to neural networks. These signals are similar to electrical waves. In electrophysiology the goal is to study the electrical properties of biological cells and tissues, and the transmission of signals. From a physics perspective, there exists a field of electrical potential within the living body, the nervous system, extracellular space and cells. Electrophysiological problems can be investigated experimentally and also theoretically by developing appropriate mathematical or computational models. Due to the enormous complexity of biological systems, it would be almost impossible to establish a detailed computational model of the electrical field, even for only a single organ (e.g. heart), including the entirety of cells comprising the neural network. In order to make computational models feasible for practical applications, we here introduce the concept of smeared fields, which represents a generalization of the previously formulated multiscale smeared methodology for mass transport in blood vessels, lymph, and tissue. We demonstrate the accuracy of the smeared finite element computational models for the electric field in numerical examples. The electrical field is further coupled with ionic mass transport within tissue composed of interstitial spaces extracellularly and by cytoplasm and organelles intracellularly. The proposed methodology, which couples electrophysiology and molecular ionic transport, is applicable to a variety of biological systems.

Smeared multiscale finite element model for electrophysiology and ionic transport in biological tissue.

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.

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

The authors present the preparation procedure and a computational model of a three-layered fibrous scaffold for prolonged drug release. The scaffold, produced by emulsion/sequential electrospinning, consists of a poly(D,L-lactic-co-glycolic acid) (PLGA) fiber layer sandwiched between two poly(ε-caprolactone) (PCL) layers. Experimental results of drug release rates from the scaffold are compared with the results of the recently introduced computational finite element (FE) models for diffusive drug release from nanofibers to the three-dimensional (3D) surrounding medium. Two different FE models are used: (1) a 3D discretized continuum and fibers represented by a simple radial one-dimensional (1D) finite elements, and (2) a 3D continuum discretized by composite smeared finite elements (CSFEs) containing the fiber smeared and surrounding domains. Both models include the effects of polymer degradation and hydrophobicity (as partitioning) of the drug at the fiber/surrounding interface. The CSFE model includes a volumetric fraction of fibers and diameter distribution, and is additionally enhanced by using correction function to improve the accuracy of the model. The computational results are validated on Rhodamine B (fluorescent drug l) and other hydrophilic drugs. Agreement with experimental results proves that numerical models can serve as efficient tools for drug release to the surrounding porous medium or biological tissue. It is demonstrated that the introduced three-layered scaffold delays the drug release process and can be used for the time-controlled release of drugs in postoperative therapy.

Preparation and modeling of three-layered PCL/PLGA/PCL fibrous scaffolds for prolonged drug release.

Metastatic cancer disease is the major cause of death in cancer patients. Because those small secondary tumors are clinically hardly detectable in their early stages, little is known about drug biodistribution and permeation into those metastatic tumors potentially contributing to insufficient clinical success against metastatic disease. Our recent studies indicated that breast cancer liver metastases may have compromised perfusion of intratumoral capillaries hindering the delivery of therapeutics for yet unknown reasons. To understand the microcirculation of small liver metastases, we have utilized computational simulations to study perfusion and oxygen concentration fields in and around the metastases smaller than 700 µm in size at the locations of portal vessels, central vein, and liver lobule acinus. Despite tumor vascularization, the results show that blood flow in those tumors can be substantially reduced indicating the presence of inadequate blood pressure gradients across tumors. A low blood pressure may contribute to the collapsed intratumoral capillary lumen limiting tumor perfusion that phenomenologically corroborates with our previously published in vivo studies. Tumors that are smaller than the liver lobule size and originating at different lobule locations may possess a different microcirculation environment and tumor perfusion. The acinus and portal vessel locations in the lobule were found to be the most beneficial to tumor growth based on tumor access to blood flow and intratumoral oxygen. These findings suggest that microcirculation states of small metastatic tumors can potentially contribute to physiological barriers preventing efficient delivery of therapeutic substances into small tumors.

Milos Kojic

Bioengineering Research and Development Center BioIRC

In Silico Clinical Trials for Cardiovascular Disease