Computational model for the regulation of extracellular ATP and adenosine in airway epithelia

Physiology and pathophysiology of human airway mucus

The mucus clearance system is the dominant mechanical host defense system of the human lung. Mucus is cleared from the lung by cilia and airflow, including both two-phase gas-liquid pumping and cough-dependent mechanisms, and mucus transport rates are heavily dependent on mucus concentration. Importantly, mucus transport rates are accurately predicted by the gel-on-brush model of the mucociliary apparatus from the relative osmotic moduli of the mucus and periciliary-glycocalyceal (PCL-G) layers. The fluid available to hydrate mucus is generated by transepithelial fluid transport. Feedback interactions between mucus concentrations and cilia beating, via purinergic signaling, coordinate Na+ absorptive vs Cl- secretory rates to maintain mucus hydration in health. In disease, mucus becomes hyperconcentrated (dehydrated). Multiple mechanisms derange the ion transport pathways that normally hydrate mucus in muco-obstructive lung diseases, e.g., cystic fibrosis (CF), chronic obstructive pulmonary disease (COPD), non-CF bronchiectasis (NCFB), and primary ciliary dyskinesia (PCD). A key step in muco-obstructive disease pathogenesis is the osmotic compression of the mucus layer onto the airway surface with the formation of adherent mucus plaques and plugs, particularly in distal airways. Mucus plaques create locally hypoxic conditions and produce airflow obstruction, inflammation, infection, and, ultimately, airway wall damage. Therapies to clear adherent mucus with hydrating and mucolytic agents are rational, and strategies to develop these agents are reviewed.

Co-localization and confinement of ecto-nucleotidases modulate extracellular adenosine nucleotide distributions

Nucleotides comprise small molecules that perform critical signaling roles in biological systems. Adenosine-based nucleotides, including adenosine tri-, di-, and mono-phosphate, are controlled through their rapid degradation by diphosphohydrolases and ecto-nucleotidases (NDAs). The interplay between nucleotide signaling and degradation is especially important in synapses formed between cells, which create signaling 'nanodomains'. Within these 'nanodomains', charged nucleotides interact with densely-packed membranes and biomolecules. While the contributions of electrostatic and steric interactions within such nanodomains are known to shape diffusion-limited reaction rates, less is understood about how these factors control the kinetics of nucleotidase activity. To quantify these factors, we utilized reaction-diffusion numerical simulations of 1) adenosine triphosphate (ATP) hydrolysis into adenosine monophosphate (AMP) and 2) AMP into adenosine (Ado) via two representative nucleotidases, CD39 and CD73. We evaluate these sequentially-coupled reactions in nanodomain geometries representative of extracellular synapses, within which we localize the nucleotidases. With this model, we find that 1) nucleotidase confinement reduces reaction rates relative to an open (bulk) system, 2) the rates of AMP and ADO formation are accelerated by restricting the diffusion of substrates away from the enzymes, and 3) nucleotidase co-localization and the presence of complementary (positive) charges to ATP enhance reaction rates, though the impact of these contributions on nucleotide pools depends on the degree to which the membrane competes for substrates. As a result, these contributions integratively control the relative concentrations and distributions of ATP and its metabolites within the junctional space. Altogether, our studies suggest that CD39 and CD73 nucleotidase activity within junctional spaces can exploit their confinement and favorable electrostatic interactions to finely control nucleotide signaling.

A physiologically-motivated model of cystic fibrosis liquid and solute transport dynamics across primary human nasal epithelia

Cystic fibrosis (CF) disease is caused by mutations affecting the gene coding for the cystic fibrosis transmembrane conductance regulator (CFTR), an anion channel expressed in the mucosal side of epithelial tissue. In the airway, dysfunctional CFTR results in a transepithelial osmotic imbalance leading to hyperabsorption of airway surface liquid mucostasis, chronic inflammation, and eventual respiratory failure. Human nasal epithelial cell cultures from healthy and CF donors were used to perform studies of liquid and solute transport dynamics at an air/liquid interface in order to emulate the in vivo airway. Then, these results were used to inform a quantitative systems pharmacology model of airway epithelium describing electrically and chemically driven transcellular ionic transport, contributions of both convective and diffusive paracellular solute transport, and osmotically driven transepithelial water dynamics. Model predictions showed CF cultures, relative to non-CF ones, have increased apical and basolateral water permeabilities, and increase paracellular permeability and transepithelial chemical driving force for a radiolabeled tracer used to track small molecule absorption. These results provide a computational platform to better understand and probe the mechanisms behind the liquid hyperabsorption and small molecule retention profiles observed in the CF airway.

Imaging extracellular ATP with a genetically-encoded, ratiometric fluorescent sensor

Extracellular adenosine triphosphate (ATP) is a key purinergic signal that mediates cell-to-cell communication both within and between organ systems. We address the need for a robust and minimally invasive approach to measuring extracellular ATP by re-engineering the ATeam ATP sensor to be expressed on the cell surface. Using this approach, we image real-time changes in extracellular ATP levels with a sensor that is fully genetically-encoded and does not require an exogenous substrate. In addition, the sensor is ratiometric to allow for reliable quantitation of extracellular ATP fluxes. Using live-cell microscopy, we characterize sensor performance when expressed on cultured Neuro2A cells, and we measure both stimulated release of ATP and its clearance by ectonucleotidases. Thus, this proof-of-principle demonstrates a first-generation sensor to report extracellular ATP dynamics that may be useful for studying purinergic signaling in living specimens.

Mathematical model reveals role of nucleotide signaling in airway surface liquid homeostasis and its dysregulation in cystic fibrosis

Mucociliary clearance is composed of three components (i.e., mucin secretion, airway surface hydration, and ciliary-activity) which function coordinately to clear inhaled microbes and other foreign particles from airway surfaces. Airway surface hydration is maintained by water fluxes driven predominantly by active chloride and sodium ion transport. The ion channels that mediate electrogenic ion transport are regulated by extracellular purinergic signals that signal through G protein-coupled receptors. These purinoreceptors and the signaling pathways they activate have been identified as possible therapeutic targets for treating lung disease. A systems-level description of airway surface liquid (ASL) homeostasis could accelerate development of such therapies. Accordingly, we developed a mathematical model to describe the dynamic coupling of ion and water transport to extracellular purinergic signaling. We trained our model from steady-state and time-dependent experimental measurements made using normal and cystic fibrosis (CF) cultured human airway epithelium. To reproduce CF conditions, reduced chloride secretion, increased potassium secretion, and increased sodium absorption were required. The model accurately predicted ASL height under basal normal and CF conditions and the collapse of surface hydration due to the accelerated nucleotide metabolism associated with CF exacerbations. Finally, the model predicted a therapeutic strategy to deliver nucleotide receptor agonists to effectively rehydrate the ASL of CF airways.

Effects of airway surface liquid height on the kinetics of extracellular nucleotides in airway epithelia

Experimental techniques aimed at measuring the concentration of signaling molecules in the airway surface liquid (ASL) often require an unrealistically large ASL volume to facilitate sampling. This experimental limitation, prompted by the difficulty of pipetting liquid from a very shallow layer (~15 μm), leads to dilution and the under-prediction of physiologic concentrations of signaling molecules that are vital to the regulation of mucociliary clearance. Here, we use a computational model to describe the effect of liquid height on the kinetics of extracellular nucleotides in the airway surface liquid coating respiratory epithelia. The model consists of a reaction-diffusion equation with boundary conditions that represent the enzymatic reactions occurring on the epithelial surface. The simulations reproduce successfully the kinetics of extracellular ATP following hypotonic challenge for ASL volumes ranging from 25 μl to 500 μl in a 12-mm diameter cell culture. The model reveals that [ATP] and [ADO] reach 1200 nM and 2200 nM at the epithelial surface, respectively, while their volumetric averages remain less than 200 nM at all times in experiments with a large ASL volume (500 μl). These findings imply that activation of P2Y2 and A2B receptors is robust after hypotonic challenge, in contrast to what could be concluded based on experimental measurements of volumetric concentrations in large ASL volumes. Finally, given the central role that ATP and ADO play in regulating mucociliary clearance, we investigated which enzymes, when inhibited, provide the greatest increase in ATP and ADO concentrations. Our findings suggest that inhibition of NTPDase1/highTNAP would cause the greatest increase in [ATP] after hypotonic challenge, while inhibition of the transporter CNT3 would provide the greatest increase in [ADO].

A mechanochemical model for auto-regulation of lung airway surface layer volume

We develop a proof-of-principle model for auto-regulation of water volume in the lung airway surface layer (ASL) by coupling biochemical kinetics, transient ASL volume, and homeostatic mechanical stresses. The model is based on the hypothesis that ASL volume is sensed through soluble mediators and phasic stresses generated by beating cilia and air drag forces. Model parameters are fit based on the available data on human bronchial epithelial cell cultures. Simulations then demonstrate that homeostatic volume regulation is a natural consequence of the processes described. The model maintains ASL volume within a physiological range that modulates with phasic stress frequency and amplitude. Next, we show that the model successfully reproduces the responses of cell cultures to significant isotonic and hypotonic challenges, and to hypertonic saline, an effective therapy for mucus hydration in cystic fibrosis patients. These results compel an advanced airway hydration model with therapeutic value that will necessitate detailed kinetics of multiple molecular pathways, feedback to ASL viscoelasticity properties, and stress signaling from the ASL to the cilia and epithelial cells.