Skip to main content

OPINION article

Front. Pediatr., 23 March 2022
Sec. Neonatology
This article is part of the Research Topic Cardio-Circulatory Support of Neonatal Transition View all 5 articles

Multiscale Modeling Is Required for the Patent Ductus Arteriosus in Preterm Infants

  • 1Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
  • 2Department of Ultrasound, Chongqing General Hospital, University of Chinese Academy of Sciences, Chongqing, China

Introduction

The ductus arteriosus (DA) in the fetus diverts most of the deoxygenated blood returning from the head, upper extremities and coronary sinus into the descending aorta, bypassing the pulmonary circulation (1, 2). The DA in full-term newborns closes within 24 to 48 h after delivery to facilitate the blood perfusion of lung tissues. However, DA often fails to close in preterm neonates. Indeed, 70% of preterm infants delivered before 28 weeks of gestation require surgical closure or pharmaceutical treatment (1). Patent DA (PDA) in preterm infants causes severe left-to-right shunting. The overloaded pulmonary perfusion causes pulmonary edema, bronchopulmonary dysplasia, pulmonary hypertension and other conditions such as hyperactive precordium and cardiomegaly (3, 4). From a hemodynamic perspective, the high resistance in the pulmonary bed and higher pressure in the pulmonary artery drive the blood flow through the DA into systemic circulation during fetal life. After delivery, the pressure gradient between pulmonary and systemic circulations reverses and exposes the pulmonary microvasculature to systemic blood pressure and increased pulmonary blood flow from PDAs (1, 5).

As a conduit for blood flow, it is critical to determine the so-called “hemodynamically significant” PDA (hsPDA), defined as the shunting of left-to-right blood flow that increases the pulmonary blood flow and decreases the systemic blood due to the amount of blood running from the descending aorta into the lung circulation (4), and imposes a threat to the survival of a preterm infant (6). Multi-center clinical trials, e.g., that described in (6), and animal (sheep) models (7) have investigated the hemodynamic implication of PDAs. The blood flow in DAs depends upon the pressure gradient between the pulmonary and systemic circulations. In the sheep model, reversed flow through the DA contributed up to 50% of total pulmonary blood flow at 30 min after the onset of pulmonary ventilation (7).

Underlying the hemodynamic phenomena in PDA and the pulmonary and systemic circulations are the fundamental changes in PDA's walls. In term infants, the smooth muscle layer of DA develops ischemic hypoxia, which drives the cell death and inflammatory cascade, leading to the remodeling of DA, and transforming the DA into a non-contractile ligament (8). In preterm infants, there is no vulnerable region of the wall that is at risk for loss of vasa vasorum flow, therefore the PDA is less likely to develop the severe degree of hypoxia that is necessary for ductus remodeling (8).

Hemodynamics in the PDA Examined From Echocardiography

Echocardiography is the gold standard for the diagnosis of PDA (4, 5). Echocardiographic assessments include measuring the ductal size, ductal flow patterns, peak systolic and end-diastolic flow velocities, and the presence of retrograde diastolic flow (5, 9). The cut-off value for hsPDA is the ductal diameter ≥ 1.5 mm (4), which should be measured from 2D ultrasonic imaging rather than color Doppler as the latter may overestimate the ductal diameter (5). Besides the ductal diameter, the flow direction and pulsatile flow pattern are also valuable for assessing psPDA. Specifically, a pulsatile, non-restrictive left-to-right Doppler pattern is most sensitive (93.5%) and specific (100%) for predicting its development (5). Furthermore, several other parameters, such as the left ventricular output, volume and pressure, can be used as a surrogate for the degree of shunting at the ductal level (5). In addition, parameters for systemic hypo-perfusion e.g., the absent or retrograde diastolic flow in the abdominal aorta (4), may also be measured from echocardiography (5).

Pertinent to hemodynamic modeling of PDAs, contemporary echocardiography can reveal the pressure gradient between the systemic and pulmonary circulations (Figure 1B) (4), which is of great value for providing boundary conditions for flow equations.

FIGURE 1
www.frontiersin.org

Figure 1. Patent ductus arteriosus (PDA) in preterm infants: (A) the left-to-right or reversed flow in PDA after birth. Note the constriction where the DA joins aorta introduces high pressure gradient according to the Bernoulli's equation; (B) echocardiography of PDA also reveals the pressure gradient between systemic and pulmonary circulations that aids hemodynamic modeling (4); (C) in silico models may be developed for the remodeling of different layers of the DA; (D) a typical pharmacokinetics model where blood flow rates are used to derive the drug amount in an organ. PA, pulmonary artery; Ao, Aorta; EC, endothelial cell.

Modeling for the Hemodynamics in PDA

In silico models have been used to simulate multiple aspects of the cardiovascular system. However, few have been developed for the circulation transitions from fetuses to neonates. Even less is the in silico models developed for PDAs. The PDA model in (10) uses two semi-cylinders (representing the pulmonary artery and aorta) joined by a plate (representing the DA), which is highly simplified and idealistic in its geometry. For in silico models to contribute to unraveling the key mechanisms underlying PDA, thoughts must be placed on the intimate connections between remodeling of the DA, the stimuli for vessel tone regulation, and the hemodynamics in it. Toward that end, we suggest a paradigm i.e., multiscale in silico modeling that closely links the hemodynamics in PDA with vessel constituents and vascular tones.

Hemodynamic models have been developed for the cardiovascular system of fetuses. For example, van den Wijingaard et al. designed a distributed model that consists of 13 arterial segments and nine vascular beds for the abnormal arterial flows in fetuses (11). Myers and Capper used a transmission line model to simulate the pulsatile arterial flow in the thoracic aorta and other arteries in fetuses (12). However, neither study included the DA and pulmonary artery in the arterial model. Blood flow simulations have been performed at a single-vessel level for the umbilical artery (13) and the ductus venosus (14). In a fetal circulation model, the DA is designed as an electronic component connecting the descending aorta and the pulmonary artery (15). Similar treatment is used where the model for PDA is based on PHYSBE and is modified in a way to include an additional flow derived from the aorta toward lungs (16).

Concerning modeling techniques, the in silico blood flow models mentioned above range from lumped parameter (or 0D) models to the transmission line and 1D models where the vessel geometry (length and diameter) is incorporated. If the purpose of the simulation is to reproduce the velocity waveform observed in echocardiography and the impedance of the vascular bed, then a coupled 1D and lumped parameter model, e.g., of (17), or coupled 1D and transmission line model, e.g., of (18) could serve the purpose. On the other hand, if the modeling goal is to investigate the detailed 3D flow patterns, e.g., the formation of vortex and flow separations after constriction in PDA, and the wall shear stress, then a 3D flow model should be applied (19, 20). A recent model simulated the 3D flow in the PDA to investigate the effects of the PDA on the flow features of the modified Blalock-Taussig shunt (21). Another study investigated the 3D flow in the PDA with respect to three typical pulmonary artery (PA) dysplasia structures and different sizes of PDA (22). Similarly, Hsia et al. developed a multi-dimensional model for the Norwood operation with a right ventricle-to-pulmonary artery shunt (23).

Hence, various modeling methods can be adopted for the left-to-right flow in the PDA in preterm infants, depending on the level of details required for the hemodynamics in the PDA. Szpinda et al. found that the reversed flow in the PDA was less than that expected from the pressure gradient (2). This is because the constriction at the joint between the PDA and the aorta introduces a higher pressure gradient according to the Bernoulli's equation (Figure 1A). Thus, the rapid flow at the constriction suits a flow model of a higher dimension (3D or 1D), whereas the resistance from the pulmonary bed can be lumped into a 0D model. For example, with a model that has a 3D simulation for the aortic arch and the shunt, 0D models for pulmonary, coronary, upper and lower body circulations imposed from the inlet outlet or inlet boundaries at the pulmonary or systemic circulation bed (22, 23).

While insights could be gained from the hemodynamics in PDA and its impacts on the pulmonary and systemic circulations from the above mentioned models, the closure of PDA is an event associated with the remodeling of the arterial wall of the DA, which none of the hemodynamic models address directly. To that end a different category of models, i.e., for wall remodeling, is required.

Multiscale Modeling for the Vascular Evolution

The DA in term infants is closed within a few hours after birth, because the increased arterial PaO2 and decreased circulating prostaglandin allow the smooth muscle media of the ductus to constrict (8). The constriction leads to ischemic hypoxia in the inner muscle wall of the DA, which consequently develops vascular endothelial growth factor that transform the ductus into a non-contractile ligament (8). In preterm infants, the intrinsic tone of the ductus is lower than that in full-term infants and has less capacity to constrict (1). Low fetal systemic arterial oxygen tension and circulating prostaglandin keep the lumen of DA patent.

Mechano-transduction plays an important role in the closure of DA: the wall shear stress acts on endothelial cells (ECs), triggering a cascade of intracellular events (24, 25). In addition, shear stress induced by transmural interstitial flow acts on the arterial media layer (Figure 1C) and induces prostaglandin production (26). In the in vitro smooth muscle cell gel model of (26), the production of prostaglandin E1 was enhanced with a 1 dyne/cm2 shear stress, but not by a 0.15 dyne/cm2 shear stress at any time.

The modeling framework for arterial remodeling has been proposed in several works (27, 28). For example, in such a model for aneurysms where the arterial wall remodels to accommodate complex arterial flows, the medial wall degradation is driven by an inflammatory response, followed by the changes in the distribution of collagen fibers. From mathematical modeling's perspective, the fibroblast-mediated collagen growth is governed by a system of ordinary differential equations to regulate the Transforming Growth Factor (TGF)-β, a key promoter of matrix deposition (28). However, efforts are required to assess whether these models are applicable to PDA and its closure.

Discussion and Conclusion

The unique anatomical structure of the DA is a “double-edged-sword”. While it results in cardiac decompensation and pulmonary distress syndrome, it may also provide the only life-sustaining conduit to preserve pulmonary or systemic blood flow (3). While many preterm infants receive surgical procedures or pharmaceutical treatments for closing PDAs, it is still under debate whether a more invasive approach (e.g., surgeries) should be adopted, or a more conservative approach (e.g., machine taking) be adopted PDA be left alone.

Given the challenges in designing clinical or experimental protocols, in silico models are of great value to the studies of many subtle phenomena of PDAs. For example, the ductus in preterm infants is more sensitive to the vasodilating effects of nitric oxide and prostaglandin than full-term infants (1). In silico models could therefore be developed to check whether the arterial constituents and blood flow contribute to this phenomenon. This requires a multiscale approach where several in silico models are linked: the flow in the PDA is solved through a 0D-1D model, whereby the results are supplied to an arterial wall model of PDA for its modeling. The former model occurs at the vessel (millimeter) level, while the latter model occurs at cellular (e.g., collagen, micrometer) level.

An in-depth understanding of the underlying vascular physiology will also aid drug development for PDA. One example is that the drug for PDA, indomethacin, reduces renal, mesenteric and cerebral blood flows due to its role as the inhibitor of vasodilators, e.g., prostaglandin. Renal and cerebral flows occur in organ-level scales, while drug molecules access celluar membrane receptors at nanometer scales (3). A flow-limited pharmacokinetics model requires the data of blood flow rate to induce the amount of drug distributed to the organ (Figure 1D). Therefore, thoughts could be placed on coupling the blood flow arriving at the PDA with the pharmacodynamics effect of the drug, and the remodeling of the DA wall that spans tissue, cell and molecule levels (29). The circulation of prostaglandin and its signaling for the smooth muscle of the PDA provides an excellent example for such a multiscale model. These events are at the very core of the pathophysiology and clinical treatment of PDAs, where the application of in silico models is yet to be tested.

Multiscale and multi-modality modeling is a promising research direction in PDA and other arterial diseases in neonates and infants. With such a strategy, biomechanical simulations are connected with imaging and histological studies of the arterial wall. A class of mechano-biological models for growth and remodeling of the arterial wall and their intimate interaction with hemodynamics, cell activity, and arterial wall mechanics are described in (30). Specifically, in the work of (20), the three-dimensional flow simulations for a subject-specific carotid artery are compared with the endarterectomy specimen's histological maps, which could provide clues for cellular level models linking to the wall shear stress (WSS) that triggers arterial plaque. A roadmap for addressing the technical challenges for incorporating the proposed growth and remodeling into existing hemodynamic simulators are discussed in (29, 30).

In conclusion, we have provided a light literature review of the pathophysiology of PDA, and some in silico models used for blood flow and arterial wall modeling. We stress that a multiscale paradigm needs to be adopted when developing future in silico models for PDA and its treatment.

Author Contributions

HH drafted the paper. XJ provided the medical context. HH and XJ reviewed and agreed the article. Both authors contributed to the article and approved the submitted version.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

HH acknowledges the financial support from the Li Ka Shing Foundation.

References

1. Clyman RI. Ibuprofen and patent ductus arteriosus. N Engl J Med. (2000) 343:728–30. doi: 10.1056/NEJM200009073431009

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Szpinda M, Szwesta A, Szpinda E. Morphometric study of the ductus arteriosus during human development. Ann Anat. (2007) 189:47–52. doi: 10.1016/j.aanat.2006.06.008

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Friedman WF, Hirschklau MJ, Printz MP, Pitlick PT, Kirkpatrick SE. Pharmacologic closure of patent ductus arteriosus in the premature infant. N Engl J Med. (1976) 295:526–9. doi: 10.1056/NEJM197609022951003

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Arlettaz R. Echocardiographic evaluation of patent ductus arteriosus in preterm infants. Front Pediatr. (2017) 5: 147. doi: 10.3389/fped.2017.00147

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Babla K, Shetty S, Kulkarni A. A clinical and echocardiographic approach to evaluation of patent ductus arteriosus in preterm infants. J Paediatr Child Health. (2020) 30: 129–34. doi: 10.1016/j.paed.2020.01.003

CrossRef Full Text | Google Scholar

6. Ellison RC, Peckham GJ, Lang P, Talner NS, Lerer TJ, Lin L, Dooley KJ, Nadas AS. Evaluation of the preterm infant for patent ductus arteriosus. Pediatr. (1983) 71:364–372. doi: 10.1542/peds.71.3.364

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Crossley KJ, Allison BJ, Polglase GR, Morley CJ, Davis PG, Hooper SB. Dynamic changes in the direction of blood flow through the ductus arteriosus at birth. J Physiol. (2009) 587:4695–704. doi: 10.1113/jphysiol.2009.174870

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Hermes-DeSantis ER, Clyman RI. Patent ductus arteriosus: pathophysiology and management. J Perinatol. (2006) 26:S14–8. doi: 10.1038/sj.jp.7211465

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Hsu K-H, Nguyen J, Dekom S, Ramanathan R, Noori S. Effects of patent ductus arteriosus on organ blood flow in infants born very preterm: a prospective study with serial echocardiography. Eur J Pediatr. (2020) 216:95–100.e2. doi: 10.1016/j.jpeds.2019.08.057

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Setchi A, Mestel AJ, Siggers JH, Parker KH, Tan MW, Wong K. Mathematical model of flow through the patent ductus arteriosus. J Math Biol. (2013) 67:1487–506. doi: 10.1007/s00285-012-0596-8

PubMed Abstract | CrossRef Full Text | Google Scholar

11. van den Wijngaard JPHM, Westerhof BE, Faber DJ, Ramsay MM, Westerhof N, van Gemert MJC. Abnormal arterial flows by a distributed model of the fetal circulation. Am J Physiol Regul. (2006) 291:R1222–33. doi: 10.1152/ajpregu.00212.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Myers LJ, Capper WL. A transmission line model of the human foetal circulatory system. Med Eng Phys. (2002) 24:285–94. doi: 10.1016/S1350-4533(02)00019-X

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Wilke DJ, Denier JP, Khong TY, Mattner TW. Pressure and flow in the umbilical cord. J Biomech. (2018) 79:78–87. doi: 10.1016/j.jbiomech.2018.07.044

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Leinan PR, Degroote J, Kiserud T, Skallerud B, Vierendeels J, Hellevik LR. Velocity profiles in the human ductus venosus: a numerical fluid structure interaction study. Biomech Model Mechanobiol. (2013) 12:1019–35. doi: 10.1007/s10237-012-0460-1

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Garcia AJ, Koschnitzky JE, Ramirez J-M. The physiological determinants of sudden infant death syndrome. Respir Physiol Neurobiol. (2013) 189:288–300. doi: 10.1016/j.resp.2013.05.032

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Bakhshali MA, Mafi M, Daneshvar S. Mathematical modelling of the patent ductus arteriosus (PDA). Math Comput Model Dyn Syst. (2013) 19:238–49. doi: 10.1080/13873954.2012.727187

CrossRef Full Text | Google Scholar

17. Formaggia L, Lamponi D, Tuveri M, Veneziani A. Numerical modeling of 1D arterial networks coupled with a lumped parameters description of the heart. Comput Methods Biomech Biomed Engin. (2006) 9:273–88. doi: 10.1080/10255840600857767

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Baker N, Clarke R, Ho H, A. coupled one dimension and transmission line model for arterial flow simulation. Int J Numer Method Biomed Eng. (2020) 36:e3327. doi: 10.1002/cnm.3327

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Formaggia L, Gerbeau JF, Nobile F, Quarteroni A. On the coupling of 3D and 1D Navier-Stokes equations for flow problems in compliant vessels. Comput Methods Appl Mech Eng. (2001) 191:561–82. doi: 10.1016/S0045-7825(01)00302-4

CrossRef Full Text | Google Scholar

20. Soloperto G, Keenan NG, Sheppard MN, Ohayon J, Wood NB, Pennell DJ, et al. Combined imaging, computational and histological analysis of a ruptured carotid plaque: a patient-specific analysis. Artery Res. (2010) 4:59–65. doi: 10.1016/j.artres.2010.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Xiong J, Sun Q, Qian Y, Hu L, Tong Z, Liu J, Liu J. Effects of patent ductus arteriosus on the hemodynamics of modified blalock–taussig shunt based on patient-specific simulation. Front Physiol. (2021) 12:707128. doi: 10.3389/fphys.2021.707128

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Xu P, Yuan H, Zhuang J, Zhang N, Jia Q, Dong Y, et al. The hemodynamics of patent ductus arteriosus in patients after central shunt operation. Comput Math Methods Med. (2021) 2021:e6675613. doi: 10.1155/2021/6675613

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Hsia T-Y, Cosentino D, Corsini C, Pennati G, Dubini G, Migliavacca F. Use of mathematical modeling to compare and predict hemodynamic effects between hybrid and surgical norwood palliations for hypoplastic left heart syndrome circulation. (2011) 124:S204–10. doi: 10.1161/CIRCULATIONAHA.110.010769

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Chien S. Mechanotransduction and endothelial cell homeostasis: the wisdom of the cell. Am J Physiol Heart Circ Physiol. (2007) 292:H1209–1224. doi: 10.1152/ajpheart.01047.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Tarbell JM, Shi Z-D, Dunn J, Jo H. Fluid mechanics, arterial disease, and gene expression. Annu Rev Fluid Mech. (2014) 46:591–614. doi: 10.1146/annurev-fluid-010313-141309

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Wang Su, Tarbell John M. Effect of fluid flow on smooth muscle cells in a 3-dimensional collagen gel model. Arterioscler Thromb Vasc Biol. (2000) 20:2220–25. doi: 10.1161/01.ATV.20.10.2220

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Watton PN, Selimovic A, Raberger NB, Huang P, Holzapfel GA, Ventikos Y. Modelling evolution and the evolving mechanical environment of saccular cerebral aneurysms. Biomech Model Mechanobiol. (2010) 10:109–32. doi: 10.1007/s10237-010-0221-y

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Aparício P, Thompson MS, Watton PN. A novel chemo-mechano-biological model of arterial tissue growth and remodelling. J Biomech. (2016) 49:2321–30. doi: 10.1016/j.jbiomech.2016.04.037

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Ho H, Suresh V, Kang W, Cooling MT, Watton PN, Hunter PJ. Multiscale modeling of intracranial aneurysms: cell signaling, hemodynamics, and remodeling. IEEE Trans Biomed. (2011) 58:2974–7. doi: 10.1109/TBME.2011.2160638

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Cardamone L, Humphrey JD. “Arterial growth and remodelling is driven by hemodynamics,” In: Ambrosi D, Quarteroni A, Rozza G, editors. Modeling of Physiological Flows. MS&A — Modeling, Simulation and Applications. Milano: Springer Milan (2012). p. 187–203.

Google Scholar

Keywords: ductus arteriosus, fetus, neonates, in silico modeling, circulation, hemodynamics

Citation: Ho H and Ji X (2022) Multiscale Modeling Is Required for the Patent Ductus Arteriosus in Preterm Infants. Front. Pediatr. 10:857434. doi: 10.3389/fped.2022.857434

Received: 18 January 2022; Accepted: 28 February 2022;
Published: 23 March 2022.

Edited by:

Stefan Kurath-Koller, Medical University Graz, Austria

Reviewed by:

Weiguang Yang, Stanford University, United States

Copyright © 2022 Ho and Ji. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Harvey Ho, aGFydmV5LmhvJiN4MDAwNDA7YXVja2xhbmQuYWMubno=; Xiaojuan Ji, aml4aWFvanVhbjIwMDMmI3gwMDA0MDsxNjMuY29t

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.