Flow-induced ATP release in patient-specific arterial geometries--a comparative study of computational models

Genetic polymorphisms of purinergic P2Y2 receptor were associated with the susceptibility to essential hypertension in Chinese postmenopausal women

Hypertension has become a prominent public health concern. Essential hypertension (EH) is a polygenic disorder caused by multiple susceptibility genes. It has been previously shown that the purinergic P2Y2 receptor (P2Y2R) regulates blood pressure; however, whether P2Y2R genetic polymorphisms correlate with EH has not been investigated in Chinese. Our study included 500 EH cases and 504 controls who are Chinese postmenopausal women. We used allele-specific polymerase chain reaction (ASPCR) to genotype five single-nucleotide polymorphism (SNPs) in the P2Y2R gene, i.e., rs4944831, rs12366239, rs1783596, rs4382936, and rs10898909. We assessed the association of P2Y2R genetic polymorphisms with EH susceptibility. The results demonstrated that P2Y2R rs4382936A was correlated with a high risk of EH; particularly, the participants with the rs4382936A allele and CA/AA/(CA+AA) genotypes were at higher risks to EH compared to the subjects with the rs4382936C allele and CC genotype. Moreover, haplotype CAG combined by rs1783596-rs4382936-rs10898909 was a susceptible haplotype for EH, whereas haplotype CCG was a protective haplotype for EH. These results may provide new evidence for applying P2Y2R genetic polymorphisms as useful markers in clinic screening or monitoring potential EH cases in a population of Chinese postmenopausal women.

Hemodynamics of vascular shunts: trends, challenges, and prospects

Vascular bypass surgery takes a significant place in the treatment of vascular disease. According to various assessments, this type of surgery is associated with almost 20 % of all vascular surgery episodes (up to 23 % according to the Federal Neurosurgical Center of Novosibirsk). Even though the problem of using of vascular grafts is obvious and natural, many problems associated with them are not still elucidated. From the mechanics' point of view, a vascular bypass is a converging or diverging tee, and the functioning of such structures still does not have strict mathematical formulations and proofs in the general case, which forces many researchers to solve specific engineering problems associated with shunting. Mathematical modeling, which is the gold standard for virtual simulations of industrial and medical problems, faces great difficulties and limitations in solving problems for vascular bypasses. Complications in the treatment of the vascular disease may follow the difficulties in mathematical modeling, and the price can be a cardiac arrest or a stroke. This work is devoted to the main aspects of the medical application of vascular bypasses and their functioning as a mechanical system, as well the mathematical aspects of their possible setup.

Modeling a Nociceptive Neuro-Immune Synapse Activated by ATP and 5-HT in Meninges: Novel Clues on Transduction of Chemical Signals Into Persistent or Rhythmic Neuronal Firing

Extracellular ATP and serotonin (5-HT) are powerful triggers of nociceptive firing in the meninges, a process supporting headache and whose cellular mechanisms are incompletely understood. The current study aimed to develop, with the neurosimulator NEURON, a novel approach to explore in silico the molecular determinants of the long-lasting, pulsatile nature of migraine attacks. The present model included ATP and 5-HT release, ATP diffusion and hydrolysis, 5-HT uptake, differential activation of ATP P2X or 5-HT3 receptors, and receptor subtype-specific desensitization. The model also tested the role of branched meningeal fibers with multiple release sites. Spike generation and propagation were simulated using variable contribution by potassium and sodium channels in a multi-compartment fiber environment. Multiple factors appeared important to ensure prolonged nociceptive firing potentially relevant to long-lasting pain. Crucial roles were observed in: (i) co-expression of ATP P2X2 and P2X3 receptor subunits; (ii) intrinsic activation/inactivation properties of sodium Nav1.8 channels; and (iii) temporal and spatial distribution of ATP/5-HT release sites along the branches of trigeminal nerve fibers. Based on these factors we could obtain either persistent activation of nociceptive firing or its periodic bursting mimicking the pulsating nature of pain. In summary, our model proposes a novel tool for the exploration of peripheral nociception to test the contribution of clinically relevant factors to headache including migraine pain.

3D time-varying simulations of Ca2+ dynamics in arterial coupled cells: A massively parallel implementation

Preferential locations of atherosclerotic plaque are strongly associated with the areas of low wall shear stress and disturbed haemodynamic characteristics such as flow detachment, flow recirculation and oscillatory flow. The areas of low wall shear stress are also associated with the reduced production of adenosine triphosphate in the endothelial layer, as well as the resulting reduced production of inositol trisphosphate (IP³ ). The subsequent variation in Ca²⁺ signalling and nitric oxide synthesis could lead to the impairment of the atheroprotective function played by nitric oxide. In previous studies, it has been suggested that the reduced IP³ and Ca²⁺ signalling can explain the correlation of atherosclerosis with induced low WSS and disturbed flow characteristics. The massively parallel implementation described in this article provides insight into the dynamics of coupled smooth muscle cells and endothelial cells mapped onto the surface of an idealised arterial bifurcation. We show that variations in coupling parameters, which model normal and pathological conditions, provide vastly different smooth muscle cell Ca²⁺ dynamics and wave propagation profiles. The extensibility of the coupled cells model and scalability of the implementation provide a solid framework for in silico investigations of the interaction between complex cellular chemistry and the macro-scale processes determined by fluid dynamics. © 2016 The Authors. International Journal for Numerical Methods in Biomedical Engineering published by John Wiley & Sons Ltd.

Patient-specific fracture risk assessment of vertebrae: A multiscale approach coupling X-ray physics and continuum micromechanics

While in clinical settings, bone mineral density measured by computed tomography (CT) remains the key indicator for bone fracture risk, there is an ongoing quest for more engineering mechanics-based approaches for safety analyses of the skeleton. This calls for determination of suitable material properties from respective CT data, where the traditional approach consists of regression analyses between attenuation-related grey values and mechanical properties. We here present a physics-oriented approach, considering that elasticity and strength of bone tissue originate from the material microstructure and the mechanical properties of its elementary components. Firstly, we reconstruct the linear relation between the clinically accessible grey values making up a CT, and the X-ray attenuation coefficients quantifying the intensity losses from which the image is actually reconstructed. Therefore, we combine X-ray attenuation averaging at different length scales and over different tissues, with recently identified 'universal' composition characteristics of the latter. This gives access to both the normally non-disclosed X-ray energy employed in the CT-device and to in vivo patient-specific and location-specific bone composition variables, such as voxel-specific mass density, as well as collagen and mineral contents. The latter feed an experimentally validated multiscale elastoplastic model based on the hierarchical organization of bone. Corresponding elasticity maps across the organ enter a finite element simulation of a typical load case, and the resulting stress states are increased in a proportional fashion, so as to check the safety against ultimate material failure. In the young patient investigated, even normal physiological loading is probable to already imply plastic events associated with the hydrated mineral crystals in the bone ultrastructure, while the safety factor against failure is still as high as five. Copyright © 2016 John Wiley & Sons, Ltd.

Mathematical model for shear stress dependent NO and adenine nucleotide production from endothelial cells

We developed a mass transport model for a parallel-plate flow chamber apparatus to predict the concentrations of nitric oxide (NO) and adenine nucleotides (ATP, ADP) produced by cultured endothelial cells (ECs) and investigated how the net rates of production, degradation, and mass transport for these three chemical species vary with changes in wall shear stress (τw). These simulations provide an improved understanding of experimental results obtained with parallel-plate flow chambers and allows quantitative analysis of the relationship between τw, adenine nucleotide concentrations, and NO produced by ECs. Experimental data obtained after altering ATP and ADP concentrations with apyrase were analyzed to quantify changes in the rate of NO production (RNO). The effects of different isoforms of apyrase on ATP and ADP concentrations and nucleotide-dependent changes in RNO could be predicted with the model. A decrease in ATP was predicted with apyrase, but an increase in ADP was simulated due to degradation of ATP. We found that a simple proportional relationship relating a component of RNO to the sum of ATP and ADP provided a close match to the fitted curve for experimentally measured changes in RNO with apyrase. Estimates for the proportionality constant ranged from 0.0067 to 0.0321 μM/s increase in RNO per nM nucleotide concentration, depending on which isoform of apyrase was modeled, with the largest effect of nucleotides on RNO at low τw (<6 dyn/cm(2)).

A reduced-dimensional model for near-wall transport in cardiovascular flows

Near-wall mass transport plays an important role in many cardiovascular processes, including the initiation of atherosclerosis, endothelial cell vasoregulation, and thrombogenesis. These problems are characterized by large Péclet and Schmidt numbers as well as a wide range of spatial and temporal scales, all of which impose computational difficulties. In this work, we develop an analytical relationship between the flow field and near-wall mass transport for high-Schmidt-number flows. This allows for the development of a wall-shear-stress-driven transport equation that lies on a codimension-one vessel-wall surface, significantly reducing computational cost in solving the transport problem. Separate versions of this equation are developed for the reaction-rate-limited and transport-limited cases, and numerical results in an idealized abdominal aortic aneurysm are compared to those obtained by solving the full transport equations over the entire domain. The reaction-rate-limited model matches the expected results well. The transport-limited model is accurate in the developed flow regions, but overpredicts wall flux at entry regions and reattachment points in the flow.

Multidimensional persistence in biomolecular data

Persistent homology has emerged as a popular technique for the topological simplification of big data, including biomolecular data. Multidimensional persistence bears considerable promise to bridge the gap between geometry and topology. However, its practical and robust construction has been a challenge. We introduce two families of multidimensional persistence, namely pseudomultidimensional persistence and multiscale multidimensional persistence. The former is generated via the repeated applications of persistent homology filtration to high-dimensional data, such as results from molecular dynamics or partial differential equations. The latter is constructed via isotropic and anisotropic scales that create new simiplicial complexes and associated topological spaces. The utility, robustness, and efficiency of the proposed topological methods are demonstrated via protein folding, protein flexibility analysis, the topological denoising of cryoelectron microscopy data, and the scale dependence of nanoparticles. Topological transition between partial folded and unfolded proteins has been observed in multidimensional persistence. The separation between noise topological signatures and molecular topological fingerprints is achieved by the Laplace-Beltrami flow. The multiscale multidimensional persistent homology reveals relative local features in Betti-0 invariants and the relatively global characteristics of Betti-1 and Betti-2 invariants.