PMID- 26441619 OWN - NLM STAT- PubMed-not-MEDLINE DCOM- 20151007 LR - 20200930 IS - 1662-5188 (Print) IS - 1662-5188 (Electronic) IS - 1662-5188 (Linking) VI - 9 DP - 2015 TI - Neurovascular coupling: a parallel implementation. PG - 109 LID - 10.3389/fncom.2015.00109 [doi] LID - 109 AB - A numerical model of neurovascular coupling (NVC) is presented based on neuronal activity coupled to vasodilation/contraction models via the astrocytic mediated perivascular K(+) and the smooth muscle cell (SMC) Ca(2+) pathway termed a neurovascular unit (NVU). Luminal agonists acting on P2Y receptors on the endothelial cell (EC) surface provide a flux of inositol trisphosphate (IP3) into the endothelial cytosol. This concentration of IP3 is transported via gap junctions between EC and SMC providing a source of sarcoplasmic derived Ca(2+) in the SMC. The model is able to relate a neuronal input signal to the corresponding vessel reaction (contraction or dilation). A tissue slice consisting of blocks, each of which contain an NVU is connected to a space filling H-tree, simulating a perfusing arterial tree (vasculature) The model couples the NVUs to the vascular tree via a stretch mediated Ca(2+) channel on both the EC and SMC. The SMC is induced to oscillate by increasing an agonist flux in the EC and hence increased IP3 induced Ca(2+) from the SMC stores with the resulting calcium-induced calcium release (CICR) oscillation inhibiting NVC thereby relating blood flow to vessel contraction and dilation following neuronal activation. The coupling between the vasculature and the set of NVUs is relatively weak for the case with agonist induced where only the Ca(2+) in cells inside the activated area becomes oscillatory however, the radii of vessels both inside and outside the activated area oscillate (albeit small for those outside). In addition the oscillation profile differs between coupled and decoupled states with the time required to refill the cytosol with decreasing Ca(2+) and increasing frequency with coupling. The solution algorithm is shown to have excellent weak and strong scaling. Results have been generated for tissue slices containing up to 4096 blocks. FAU - Dormanns, Katharina AU - Dormanns K AD - UC HPC Supercomputing Centre, University of Canterbury Christchurch, New Zealand. FAU - Brown, Richard G AU - Brown RG AD - Institution of Mathematical Sciences, Massey University Palmerston North, New Zealand. FAU - David, Tim AU - David T AD - UC HPC Supercomputing Centre, University of Canterbury Christchurch, New Zealand. LA - eng PT - Journal Article DEP - 20150915 PL - Switzerland TA - Front Comput Neurosci JT - Frontiers in computational neuroscience JID - 101477956 PMC - PMC4569750 OTO - NOTNLM OT - agonistic behavior OT - computational biology OT - neurovascular coupling OT - neurovascular unit OT - parallel computing EDAT- 2015/10/07 06:00 MHDA- 2015/10/07 06:01 PMCR- 2015/01/01 CRDT- 2015/10/07 06:00 PHST- 2015/05/04 00:00 [received] PHST- 2015/08/24 00:00 [accepted] PHST- 2015/10/07 06:00 [entrez] PHST- 2015/10/07 06:00 [pubmed] PHST- 2015/10/07 06:01 [medline] PHST- 2015/01/01 00:00 [pmc-release] AID - 10.3389/fncom.2015.00109 [doi] PST - epublish SO - Front Comput Neurosci. 2015 Sep 15;9:109. doi: 10.3389/fncom.2015.00109. eCollection 2015.