Graph Theory for Modeling and Analysis of the Human Lymphatic System

Computational fluid dynamic modeling of the lymphatic system: a review of existing models and future directions

Historically, research into the lymphatic system has been overlooked due to both a lack of knowledge and limited recognition of its importance. In the last decade however, lymphatic research has gained substantial momentum and has included the development of a variety of computational models to aid understanding of this complex system. This article reviews existing computational fluid dynamic models of the lymphatics covering each structural component including the initial lymphatics, pre-collecting and collecting vessels, and lymph nodes. This is followed by a summary of limitations and gaps in existing computational models and reasons that development in this field has been hindered to date. Over the next decade, efforts to further characterize lymphatic anatomy and physiology are anticipated to provide key data to further inform and validate lymphatic fluid dynamic models. Development of more comprehensive multiscale- and multi-physics computational models has the potential to significantly enhance the understanding of lymphatic function in both health and disease.

Graph theory applications in congenital heart disease

Graph theory can be used to address problems with complex network structures. Congenital heart diseases (CHDs) involve complex abnormal connections between chambers, vessels, and organs. We proposed a new method to represent CHDs based on graph theory, wherein vertices were defined as the spaces through which blood flows and edges were defined by the blood flow between the spaces and direction of the blood flow. The CHDs of tetralogy of Fallot (TOF) and transposition of the great arteries (TGA) were selected as examples for constructing directed graphs and binary adjacency matrices. Patients with totally repaired TOF, surgically corrected d-TGA, and Fontan circulation undergoing four-dimensional (4D) flow magnetic resonance imaging (MRI) were included as examples for constructing the weighted adjacency matrices. The directed graphs and binary adjacency matrices of the normal heart, extreme TOF undergoing a right modified Blalock-Taussig shunt, and d-TGA with a ventricular septal defect were constructed. The weighted adjacency matrix of totally repaired TOF was constructed using the peak velocities obtained from 4D flow MRI. The developed method is promising for representing CHDs and may be helpful in developing artificial intelligence and conducting future research on CHD.

Mathematical Modeling of Lymph Node Drainage Function by Neural Network

The lymph node (LN) represents a key structural component of the lymphatic system network responsible for the fluid balance in tissues and the immune system functioning. Playing an important role in providing the immune defense of the host organism, LNs can also contribute to the progression of pathological processes, e.g., the spreading of cancer cells. To gain a deeper understanding of the transport function of LNs, experimental approaches are used. Mathematical modeling of the fluid transport through the LN represents a complementary tool for studying the LN functioning under broadly varying physiological conditions. We developed an artificial neural network (NN) model to describe the lymph node drainage function. The NN model predicts the flow characteristics through the LN, including the exchange with the blood vascular systems in relation to the boundary and lymphodynamic conditions, such as the afferent lymph flow, Darcy’s law constants and Starling’s equation parameters. The model is formulated as a feedforward NN with one hidden layer. The NN complements the computational physics-based model of a stationary fluid flow through the LN and the fluid transport across the blood vessel system of the LN. The physical model is specified as a system of boundary integral equations (IEs) equivalent to the original partial differential equations (PDEs; Darcy’s Law and Starling’s equation) formulations. The IE model has been used to generate the training dataset for identifying the NN model architecture and parameters. The computation of the output LN drainage function characteristics (the fluid flow parameters and the exchange with blood) with the trained NN model required about 1000-fold less central processing unit (CPU) time than computationally tracing the flow characteristics of interest with the physics-based IE model. The use of the presented computational models will allow for a more realistic description and prediction of the immune cell circulation, cytokine distribution and drug pharmacokinetics in humans under various health and disease states as well as assisting in the development of artificial LN-on-a-chip technologies.