Document Detail

Convective washout reduces the antidiarrheal efficacy of enterocyte surface-targeted antisecretory drugs.
Jump to Full Text
MedLine Citation:
PMID:  23359285     Owner:  NLM     Status:  MEDLINE    
Secretory diarrheas such as cholera are a major cause of morbidity and mortality in developing countries. We previously introduced the concept of antisecretory therapy for diarrhea using chloride channel inhibitors targeting the cystic fibrosis transmembrane conductance regulator channel pore on the extracellular surface of enterocytes. However, a concern with this strategy is that rapid fluid secretion could cause convective drug washout that would limit the efficacy of extracellularly targeted inhibitors. Here, we developed a convection-diffusion model of washout in an anatomically accurate three-dimensional model of human intestine comprising cylindrical crypts and villi secreting fluid into a central lumen. Input parameters included initial lumen flow and inhibitor concentration, inhibitor dissociation constant (K(d)), crypt/villus secretion, and inhibitor diffusion. We modeled both membrane-impermeant and permeable inhibitors. The model predicted greatly reduced inhibitor efficacy for high crypt fluid secretion as occurs in cholera. We conclude that the antisecretory efficacy of an orally administered membrane-impermeant, surface-targeted inhibitor requires both (a) high inhibitor affinity (low nanomolar K(d)) to obtain sufficiently high luminal inhibitor concentration (>100-fold K(d)), and (b) sustained high luminal inhibitor concentration or slow inhibitor dissociation compared with oral administration frequency. Efficacy of a surface-targeted permeable inhibitor delivered from the blood requires high inhibitor permeability and blood concentration (relative to K(d)).
Byung-Ju Jin; Jay R Thiagarajah; A S Verkman
Related Documents :
23841875 - Deconstructing activation events in rhodopsin.
24900225 - Antiproliferative and differentiating activities of a novel series of histone deacetyla...
23695255 - Structure of 2-haloacid dehalogenase from pseudomonas syringae pv. tomato dc3000.
24982075 - Molecular mechanisms of substrate recognition and specificity of new delhi metallo β-l...
2684275 - Pyridoxal 5'-phosphate dependent histidine decarboxylase: overproduction, purification,...
11034785 - Distinct effect of benzalkonium chloride on the esterase and aryl acylamidase activitie...
Publication Detail:
Type:  Journal Article; Research Support, N.I.H., Extramural    
Journal Detail:
Title:  The Journal of general physiology     Volume:  141     ISSN:  1540-7748     ISO Abbreviation:  J. Gen. Physiol.     Publication Date:  2013 Feb 
Date Detail:
Created Date:  2013-01-29     Completed Date:  2013-07-11     Revised Date:  2014-05-23    
Medline Journal Info:
Nlm Unique ID:  2985110R     Medline TA:  J Gen Physiol     Country:  United States    
Other Details:
Languages:  eng     Pagination:  261-72     Citation Subset:  IM    
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Antidiarrheals / administration & dosage*,  pharmacokinetics*
Computer Simulation
Diarrhea / drug therapy*,  physiopathology*
Enterocytes / drug effects,  metabolism*
Metabolic Clearance Rate
Models, Biological*
Grant Support
Reg. No./Substance:
Comment In:
J Gen Physiol. 2013 Aug;142(2):173   [PMID:  23898009 ]
J Gen Physiol. 2013 Aug;142(2):171-2   [PMID:  23898008 ]

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): J Gen Physiol
Journal ID (iso-abbrev): J. Gen. Physiol
Journal ID (hwp): jgp
ISSN: 0022-1295
ISSN: 1540-7748
Publisher: The Rockefeller University Press
Article Information
Download PDF
© 2013 Jin et al.
Received Day: 23 Month: 8 Year: 2012
Accepted Day: 3 Month: 1 Year: 2013
Print publication date: Month: 2 Year: 2013
Volume: 141 Issue: 2
First Page: 261 Last Page: 272
PubMed Id: 23359285
ID: 3557305
Publisher Id: 201210885
DOI: 10.1085/jgp.201210885

Convective washout reduces the antidiarrheal efficacy of enterocyte surface–targeted antisecretory drugs Alternate Title:Convection–diffusion model of intestine
Byung-Ju Jin12
Jay R. Thiagarajah123
A.S. Verkman12
1Department of Medicine and 2Department of Physiology, University of California, San Francisco, San Francisco, CA 94143
3Department of Pediatrics, Massachusetts General Hospital, Boston, MA 02114
Correspondence: Correspondence to Alan S. Verkman:


Excessive intestinal fluid secretion occurs in enterotoxin-mediated secretory diarrheas caused by Vibrio cholera and enterotoxigenic Escherichia coli (Field, 1979). The rate-limiting step in fluid secretion is chloride transport from the enterocyte cytoplasm into the intestinal lumen, which creates the electroosmotic force driving sodium and water secretion (Murek et al., 2010; Venkatasubramanian et al., 2010; Thiagarajah and Verkman, 2012). Cell culture and animal models (Chao et al., 1994; Gabriel et al., 1994; Thiagarajah et al., 2004) indicate that elevation in cyclic nucleotides caused by bacterial enterotoxins activates the CFTR, a chloride channel expressed on the luminal surface of enterocytes. CFTR inhibition is thus predicted to be of clinical benefit as antisecretory therapy in diarrheas caused by bacterial enterotoxins (Al-Awqati, 2002; Zhang et al., 2012).

We identified a class of small molecules, the glycine hydrazides and the related malonic acid hydrazides (MalH), as CFTR inhibitors that target the extracellular-facing pore of CFTR (Muanprasat et al., 2004). An extracellular site-of-action was suggested by patch-clamp measurements showing outwardly rectifying whole-cell currents and rapid single-channel flicker (Muanprasat et al., 2004), and proven from CFTR inhibition by membrane-impermeant MalH–polyethylene glycol conjugates (Sonawane et al., 2006). Subsequently, multivalent membrane-impermeant conjugates of MalH with lectins (Sonawane et al., 2007) and polyethylene glycols (Sonawane et al., 2008) were synthesized with Kd < 100 nM for inhibition of CFTR chloride conductance. These membrane-impermeant, nonabsorbable conjugates showed antisecretory efficacy in closed-loop and sucking mouse models of cholera. Membrane-permeable, absorbable glycine hydrazide analogues were also synthesized for potential therapy of polycystic kidney disease, in which cyst growth is CFTR dependent (Yang et al., 2008).

Nonabsorbable CFTR inhibitors are potentially useful for antisecretory therapy because of their minimal systemic exposure. A glycine hydrazide analogue, iOWH032 (de Hostos et al., 2011), with modest CFTR inhibition potency (Kd of ∼8 µM), is currently in early-stage clinical trials for cholera therapy. A natural product, crofelemer, which is thought to act by inhibition of CFTR and calcium-activated chloride channels (Tradtrantip et al., 2010), is also in clinical trials (Cottreau et al., 2012). Crofelemer consists of large proanthocyanidin oligomers that are predicted to be membrane impermeable and thus externally acting.

A potentially important though largely ignored concern in the application of externally targeted inhibitors is washout as a result of convection caused by rapid fluid secretion. Convection effects cannot be addressed easily in animal models because of differences from humans in crypt–villus geometry and fluid secretion rates, as well as the lack of suitable animal models of secretory diarrhea. As diagrammed in Fig. 1 A, the intestine contains a dense array of long, narrow cylindrical crypts that secrete fluid into a central lumen. Because of fluid convection, crypt fluid secretion is predicted to reduce inhibitor concentration at the enterocyte surface and hence reduce inhibitor efficacy. Further, if therapeutic inhibitor concentration is not sustained in the intestinal lumen, its antisecretory effect would diminish rapidly unless inhibitor dissociation from its target at the enterocyte surface is slow compared with the time between oral doses. Here, we model CFTR inhibitor washout to predict the conditions (concentration/target affinity, dissociation rate) under which an extracellularly targeted antisecretory drug could be effective in reducing intestinal fluid secretion.

Model formulation

Fig. 1 A diagrams a cross section of a generic intestinal villus–crypt unit, as found in mid-jejunum, which consists of a long, narrow aqueous compartment bounded by a cylindrical layer of epithelial cells. The crypts project in an outward radial direction from the intestinal wall, and the villi project in an inward radial direction. The inhibitor target, CFTR, is heterogeneously distributed in the lumen-facing (apical) membrane of epithelial cells along the crypt–villus axis, with its extracellular surface exposed to the luminal fluid. Based on published anatomical data in human male adult mid-jejunum (Loehry and Creamer, 1969; Marsh and Swift, 1969; Trbojević-Stanković et al., 2010), typical crypt length is 150 µm and inner diameter (of the aqueous lumen) is 20–25 µm; typical villus length is 350 µm and inter-crypt spacing is ∼100 µm. Volume flux across the epithelium at location z along the crypt–villus, Jv(z), depends on CFTR density, the magnitude of the secretory stimulus, and the concentration of inhibitor contacting the villus–crypt wall. In cholera, typical fluid secretion rate in mid-jejunum is ∼4.5 ml/h per cm intestine (Banwell et al., 1970; Van Loon et al., 1992; Bearcroft et al., 1997), which, together with a crypt density of 12,500 cm−2 (Trbojević-Stanković et al., 2010) and a jejunal diameter of 4–5 cm, corresponds to a single-crypt fluid secretion rate of ∼4 × 10−7 cm3/min (∼7 × 10−2 µL/cm2/s at the crypt surface). Fluid viscosity in the crypt lumen is ∼10-fold greater than that of water, as measured by photobleaching and dye polarization (Naftalin et al., 1995). Material properties and baseline parameters are listed in Table 1. For a single cylindrical crypt–villus unit, the computed parameters include inhibitor concentration Ci (r,z) and velocity Vi (r,z), in which cylindrical symmetry allows specification of location by radial dimension r and axial position z.

Fig. 1 B shows the generic three-dimensional geometry of the intestine, with crypts and villi projecting radially from the cylindrical intestinal wall. The diameter of the intestinal lumen is orders of magnitude greater than crypt length. Fluid entering the intestinal lumen has a specified inhibitor concentration and flow velocity, which may vary in time (because of intermittent oral administration). Because the intestinal lumen is cylindrically symmetric and large compared with crypt dimensions, diffusion–convection for the geometry shown in Fig. 1 B can be simplified to that of a truncated radial wedge of fluid with attached crypt–villus units as diagrammed in Fig. 1 C, as described in Results.

List of parameters
Parameters for single-crypt/villus computations:

D, diffusion coefficient in crypt/villus lumen (cm2/s)

Jv(z), volume flux across the epithelium at location z along the crypt (µL/cm2/s)

Jvo(z), volume flux across the epithelium in the absence of inhibitor (µL/cm2/s)

ΣJv, total volume flux from crypt (µL/s/crypt)

Ci(r,z), inhibitor concentration at location r,z (µM)

Vi(r,z), vectorial flow velocity at location r,z (cm/s)

Pi(r,z), pressure at location r,z (Pa)

ΔP, pressure difference between Pi(r,z) and reference pressure in the intestine (Pa)

Co, inhibitor concentration at outer boundary (µM)

Jb(z), transepithelial volume flux of (permeable) inhibitor (mol/cm2/s)

Jbo(z), initial transepithelial volume flux of inhibitor (mol/cm2/s)

Pinh, transepithelial inhibitor permeability (cm/s)

Cb, inhibitor concentration in blood (µM)

Kd, inhibitor–CFTR binding dissociation constant (µM)

k-1, dissociation rate constant (min−1)

k1, association rate constant (µM−1 min−1)

ρ, fluid density (kg/m3)

η, fluid dynamic viscosity (Pa · s)

Δt, computation time step (s)

L, crypt length (µm)

d, crypt diameter (µm)

fb(z), occupied fraction of inhibitor binding site

Additional parameters for multi-crypt/villus computations:

DL, diffusion coefficient in the intestinal lumen (cm2/s)

DM, diffusion coefficient in the mucosa (crypt/villus region) (cm2/s)

ηL, fluid dynamic viscosity in the lumen (Pa · s)

ηM, fluid dynamic viscosity in the mucosa (Pa · s)

LV, villus length (µm)

LC, crypt length (µm)

dV, villus diameter (µm)

dC, crypt diameter (µm)

dL, intestinal lumen diameter (µm)

U(x,y,z), inhibitor velocity in the intestinal lumen (cm/s)

Umean, mean inhibitor velocity in the lumen (cm/s)

s, distance from the center of the lumen (cm)

so, distance from the center of the lumen to the top of the villus (cm)

Model computations

Models of convection–diffusion in single crypt–villus units (two-dimensional) and the full intestine (three-dimensional) were implemented using COMSOL Multiphysics (version 3.4; COMSOL). The model is specified by a diffusion equation describing inhibitor convection–diffusion coupled with the Navier–Stokes equation describing the fluid flow field. The convection–diffusion equation is

[Formula ID: fd1]
which contains convection and diffusion terms. The velocity field, Vi, in the convection term was computed from the Navier–Stokes equation for an incompressible fluid and the continuity equation,
[Formula ID: fd2]
[Formula ID: fd3]
where ρ is fluid density, Pi is pressure, and η is dynamic viscosity.

For computations on single-crypt/villus units, as shown in Fig. 2 A, constant inhibitor concentration was imposed at the outer boundary and an insulation boundary condition (no flux through membrane) was imposed at the crypt/villus surface. Volume flux at the crypt surface in the steady state, Jv(z), depends on volume flux in the absence of inhibitor, Jvo(z), and local inhibitor concentration, Ci(ro,z),

[Formula ID: fd4]
Jv(z)= Jvo(z)/[1+Ci(ro,z)/Kd],
where Kd is the equilibrium dissociation constant for inhibitor-target (CFTR) binding.

To model a membrane-permeable inhibitor in which inhibitor is transported from the circulation into the crypt/villus lumen, the transepithelial flux of inhibitor, Jb(z), is described by

[Formula ID: fd5]
Jb(z)= Pinh⋅[Cb−Ci(r,z)]= Jbo(z)⋅[1−Ci(ro,z)/Cb],
where Pinh describes inhibitor permeability across the intestinal epithelium, and Cb is inhibitor concentration in blood. A convective flux boundary condition was imposed at the outer boundary and at the crypt/villus surface.

In multi-crypt computations, the same flux boundary condition is imposed at the crypt surface. The parabolic velocity profile at the lumen inlet was calculated from mean lumen velocity as

[Formula ID: fd6]
U(x,y,z)= 2Umean⋅[1−(s/so)2],
where so is the distance from the lumen center to the top of villi, and s is the distance from the center of the lumen. Symmetric boundary conditions are imposed at the two side walls and a slip boundary condition at the top surface. Details on boundary conditions are provided in the Appendix (Figs. A1 and A2, and Table 2).

The computation time to obtain the steady-state solution was ∼10 min for single-crypt computations and 4–24 h for multi-crypt computations, as performed on an HP Z600 workstation (Xeon E5645 CPU and 32G RAM; Intel). The time step was automatically computed in COMSOL Multiphysics from mesh size and property value variation. Computation validation studies are provided in the supplemental text.

Inhibition of net fluid secretion was computed as the ratio of total crypt–villus fluid secretion in the presence versus absence of inhibitor,

[Formula ID: fd7]
% inhibition = 100⋅(1−∑Jv/∑Jvo).

For modeling the kinetics of inhibitor washout, the fraction of bound inhibitor, fb (0 to 1, fb = Ci/(Ci + Kd) in the steady state), is described by the differential equation,

[Formula ID: fd8]
dfb(z)/dt = k1⋅[1−fb(z)]⋅ Ci−k−1⋅fb(z),
where k−1 is inhibitor dissociation rate constant, and k1 is inhibitor bimolecular association rate constant, subject to the condition, Kd = k−1/k1.

The percent inhibition of fluid secretion along the length of the intestine was deduced from the computational results done for small segments, in which inhibitor concentration in each segment is diluted progressively because of fluid secretion. Conservation of inhibitor molecules (Eq. 9) and of fluid volume (Eq. 10) requires

[Formula ID: fd9]
UoutCout = UinCin
[Formula ID: fd10]
Uout⋅AL= Uin⋅AL+∑Jv⋅(πdL/w),
where Cin and Cout are inhibitor concentration at the inlet and the outlet in each segment, Uin and Uout are mean lumen velocity, AL is luminal cross-sectional area, dL is lumen diameter, w is segment length, and ∑Jv (equal to ∑Jvo [1 − 0.01 · % inhibition]) is total single-segment secretion rate. Percent inhibition at each segment was specified by an empirical fit of computed results at the segment.

Total percent inhibition (over the length of the intestine) is the ratio of integrated secreted fluid without versus with inhibitor. Total secreted fluid in the absence of inhibitor was computed by summation of flux from each segment,

[Formula ID: fd11]
Secreted fluid (no inhibitor)=∑Jvo⋅(πdL/w)⋅(Lint/Δx),
where ∑Jvo is initial single-segment secretion rate, dL is lumen diameter, w is segment length, Lint intestinal length, and Δx is segment length. Total secreted fluid in the presence of inhibitor was computed from the difference in luminal flow in the final segment without versus with inhibitor, so that percent inhibition becomes
[Formula ID: fd12]
total % inhibition = 100⋅[1−(UoutN−Uinl)AL/{∑Jvo⋅(πdL/w)⋅(Lint/Δx)}],
where UoutN is mean flow velocity in the final segment, and Uin1 is mean inlet flow velocity at the first segment, where UoutN is calculated iteratively applying conservation conditions (Eqs. 9 and 10).

Online supplemental material

The online supplemental material includes additional information on the mathematical modeling, including the validations of computation time and the number of mesh elements, and details of the secretion ratio (Jvo) calculation. It is available at

Single crypt–villus computations for a membrane-impermeant, surface-targeted inhibitor

We first analyzed convection–diffusion for a single crypt–villus unit to explore the main features of the model. Fig. 2 A (left) shows a cross section of the crypt–villus structure in mid-jejunum (where the majority of fluid secretion occurs in cholera), which consists of a 350-µm-long villus shown aligned with a 150-µm-long crypt. Based on reported data on CFTR distribution (Jakab et al., 2011), relative ∑Jv(ro,z) in the absence of inhibitor was taken as 1.0 (lower 50% of crypt), 0.6 (upper 50% of crypt), 0.28 (lower 50% of villus), and 0.15 (upper 50% of villus). Fig. 2 A (middle) shows mesh element in the two-dimensional computation region, which, by symmetry, solves the equivalent three-dimensional crypt–villus problem (Fig. 2 A, right).

For the two-dimensional computations, the intervillous region was modeled as a simple cylinder approximating the actual three-dimensional geometry, which is modeled below. We note that the further simplification of our model to one dimension is not possible because of the need to impose unrealistic boundary conditions, which necessitated the inclusion of a “buffer zone” at the crypt–villus outlet in the single-crypt computations. As fluid leaves the crypt–villus unit it takes with it solute and so produces a complex two-dimensional profile extending from the outlet that can be quite large for high flow velocities as found in diarrhea.

Fig. 2 B shows steady-state profiles of drug concentration (left) and fluid flow velocity (middle) for different (dimensionless) Co/Kd, in which Co/Kd is fixed at the outer computational region. The increase in pressure is negligible, as shown in the right panel of Fig. 2 B for Co/Kd = 100. Convective effects are observed in the profiles, with marked axial concentration and velocity gradients. At fixed fluid secretion rate (in the absence of inhibitor) of Jvo = 4 × 10−2 µL/cm2/s, increasing Co/Kd to 100 allowed inhibitor access deep in crypts (Fig. 2 B, left). Increased inhibitor concentration reduced fluid velocity. Fig. 2 C shows the steady-state drug concentration and fluid flow velocity profiles at fixed Co/Kd of 1 for different Jvo, and the pressure profile for Jvo = 10−1 µL/cm2/s. Increasing Jvo reduced inhibitor concentration in crypt-villus units because of convection, with a very small pressure increase of <2 Pa (compared with atmospheric pressure of 105 Pa).

Results from computations as in Fig. 2 are summarized in Fig. 3. Fig. 3 A shows the percent inhibition of secreted fluid for the crypt–villus geometry of the mid-jejunum, as deduced from volume flux exiting the crypt without inhibitor versus with inhibitor. The percent inhibition of fluid secretion increased from 0 to 100% with increasing Co/Kd, as expected, although the apparent IC50 for inhibition was substantially increased with increasing Jvo. For low, near static flow (Jvo = 10−3 µL/cm2/s), Co/Kd of 1 produced ∼50% inhibition of fluid secretion, as expected. For Jvo of 7 × 10−2 µL/cm2/s, which is typical of cholera, inhibition of net fluid secretion by 50% required Co/Kd of ∼200, and clinically relevant inhibition by 90% required Co/Kd of ∼750. Fig. 3 B shows similar computations for the colon, which consists of long, narrow crypts (length 430 µm, diameter 16 µm) without villi. A similar dependence of percent inhibition of fluid secretion on Co/Kd and Jvo was seen, in which 50% inhibition of fluid secretion required Co/Kd of ∼40 for Jvo = 7 × 10−3 µL/cm2/s, typical of cholera in colon.

Single crypt/villus computations for a membrane-permeable, surface-targeted inhibitor

Computations were also done for a membrane-permeable, surface-targeted CFTR inhibitor in which the inhibitor enters the crypt–villus lumen from the blood compartment rather than the intestinal lumen. We assume constant inhibitor concentration in blood, with entry into the crypt–villus lumen determined by the (blood-to-lumen) inhibitor concentration gradient and enterocyte permeability. The hypothesis is that high blood inhibitor concentration and enterocyte permeability are required for antisecretory efficacy of a membrane-permeant, externally targeted inhibitor. Our reasoning, as diagrammed in Fig. 4 A, is that inhibitor molecules crossing the enterocyte barrier are convected away from the luminal surface, impairing their ability to inhibit CFTR.

Fig. 4 B shows drug concentration profiles for different drug permeability coefficients for transepithelial transport from blood into the crypt–villus lumen, Pinh (left), and for different fluid secretion rates, Jvo. Convective washout of luminal inhibitor is seen as reduced Ci/Cb with increasing Jvo. Fig. 4 C summarizes the percent inhibition of net secreted fluid as a function of Pinh for different Jvo at fixed Cb = 10 µM and Kd = 10 µM. At a low Jvo of 10−3 µL/cm2/s, half-maximal inhibition required Pinh of ∼10−7 cm/s, whereas at high Jvo of 7 × 10−2 µL/cm2/s, typical of cholera, 50% of maximal inhibition required an ∼100-fold increased Pinh of ∼10−5 cm/s. Maximal inhibition was 50% for Cb = Kd, as expected. Similar computations are shown in Fig. 4 D for Cb = 1 µM and Kd = 0.1 µM, where maximal inhibition is ∼100%. At low Jvo of 10−3 µL/cm2/s, half-maximal inhibition required Pinh of ∼10−9 cm/s, whereas at higher Jvo of 7 × 10−2 µL/cm2/s, half-maximal inhibition required Pinh of ∼5 × 10−6 cm/s. These computations show that convective washout effects can occur for a membrane-permeable, surface-targeted inhibitor.

Multi-crypt computations with three-dimensional intestinal geometry

Computations were also done for the three-dimensional intestine as diagrammed in Fig. 1 B. Our approach was to first analyze percent inhibition versus Co/Kd (as done in Fig. 3 A for single crypt–villus units) for a short segment of intestine (mid-jejunum), and then extend the computation to model many meters of intestine. Fig. 5 A shows an en face view of the crypt–villus geometry in this model, as deduced from electron microscopy data (Loehry and Creamer, 1969; Marsh and Swift, 1969). Finger-shaped villi of 100–250-µm diameter extend out of the plane, with inter-crypt distance of ∼100 µm and a crypt to villus ratio of ∼3. We modeled crypts and villi as cylinders with semispherical outer surfaces. Villus length was taken as 350 µm, with an outer diameter of 120 µm and inter-villus spacing of 150 µm. Crypt length was taken as 150 µm, with an inner diameter of 20 µm and inter-crypt distance of 100 µm. The villi to crypt ratio was 3, giving a crypt density of 13,333 cm−2, similar to the measured value of 12,500 cm−2 (Trbojević-Stanković et al., 2010). Fig. 6 A (middle) shows a small segment of 1-mm length and 150-µm width used as the computational volume (chosen to allow a crypt to villus ratio of 3). The computational volume is composed of two layers, mucosa and lumen, of different viscosities. Fig. 5 A (right) shows 80,000–150,000 mesh elements in the computation volume.

Fig. 5 B shows the steady-state profiles of drug concentration. Similar to the profiles in Fig. 2 B, increasing Jvo produced greater convection, reducing inhibitor access deep in crypts. Fig. 5 C shows percent net fluid inhibition versus Co/Kd for different Jvo, analogous to the single crypt–villus analysis in Fig. 3 C. Co/Kd was ∼1 for low Jvo, as expected, increasing substantially for higher Jvo. These computations for three-dimensional crypt–villus geometry support the robustness of the conclusion that greatly increased Co/Kd is needed for effective inhibition at high fluid secretion rates.

The short segment of intestinal surface with anatomically accurate three-dimensional geometry was extended to model the human mid-jejunum with 5-mm length and 5-cm diameter (Fig. 5 D, top). Assuming a well-mixed intestinal lumen beyond the level of villus tips, the steady-state profile of fluid inhibition along the intestine was computed as described under Model computations in Materials and methods. Fig. 5 D (bottom) shows percent inhibition along the intestine for different initial lumen axial flow velocities, Umean, representing different axial starting points along the length of the intestine. The percent inhibition decreases along the length of intestine as inhibitor is diluted by secreted fluid from crypts and villi. Fig. 5 E shows deduced percent inhibition versus Co/Kd plots for different Umean for Jvo of 7 × 10−2 µL/cm2/s, typical of cholera. Because of inhibitor dilution effects, the percent inhibition versus Co/Kd profiles is right-shifted, indicating the requirement of high Co/Kd (>1,000) for efficient reduction in intestinal fluid secretion. Even greater Co/Kd would be predicted without the assumption of a well-mixed intestinal lumen.

Inhibitor washout kinetics limits antidiarrheal efficacy

The computations described above for an orally administered (membrane-impermeable) inhibitor were done for the steady state, with constant inhibitor concentration entering the intestinal lumen (at the boundary of the computation region). However, because of intermittent oral administration, inhibitor concentration varies over time, with little inhibitor present in the intestinal lumen much of the time. To model inhibitor washout, inhibitor concentration in the lumen was reduced at time 0 from Co/Kd to 0 (step-function). Fig. 6 A shows the kinetics of inhibitor washout for a single crypt–villus unit for dissociation rate constant, k−1 = 10−1 min−1. Washout occurs over time as CFTR inhibitor dissociation releases inhibitor molecules that are convected out of the crypt. In addition to crypt–villus geometry, the determinants of washout kinetics include k−1, fluid secretion rate, and Co/Kd. Fig. 6 B shows the kinetics of loss of inhibition of fluid secretion for mid-jejunum, with parameters chosen to give ∼90% inhibition of fluid secretion before washout. For a k−1 of 10−1 min−1 (well below that expected for small-molecule inhibitors with micromolar Kd) and a high Co/Kd of 10, washout occurs in tens of minutes. Very low k−1 is required to greatly slow inhibitor washout. For the three-dimensional crypt–villus geometry, predicted washout is even more rapid, occurring over a few minutes for a k−1 of 10−1 min−1 (Fig. 6 C).


CFTR offers a unique target for antisecretory therapy of enterotoxin-mediated diarrheas because it is expressed on the luminal membrane of enterocytes and is hence exposed to the aqueous-phase contents of the intestinal lumen. The luminal exposure of CFTR allows the possibility of developing inhibitors that are targeted to its external surface that could act after oral administration and without systemic absorption. High throughput screening of chemically diverse small-molecule libraries identified glycine hydrazides and related analogues that blocked the external CFTR pore directly (Muanprasat et al., 2004). The glycine hydrazides were engineered into nonabsorbable macromolecular conjugates with high affinity binding to CFTR and slow dissociation (Sonawane et al., 2007, 2008). However, because of the high volume of fluid that can be secreted in diarrhea, inhibitor washout from the surface of the intestine is a concern. Convection–diffusion, as fluid is secreted out of the long, narrow intestinal crypts, is predicted to reduce inhibitor concentration at the luminal surface of the crypt epithelial cells. Washout effects can also be important for systemically absorbed inhibitors that target an extracellular site on crypt epithelial cells, as inhibitor dissociation from the surface and subsequent convection out of crypts can reduce surface inhibitor concentration. Although orally administered nonabsorbable drugs are in use (Charmot, 2012), such as cholesterol-binding resins and antimicrobials, drug targeting to the external site of a luminal receptor on the intestinal epithelium is a new concept in drug delivery. Because there was no precedent for existing drugs using this targeting strategy, or data to evaluate the magnitude or kinetics of washout effects, we modeled the CFTR–inhibitor system in the intestine using data applicable to secretory diarrhea, although our model would apply in general to drugs that target luminal receptors on enterocytes.

Qualitatively, the effects of convection and diffusion on washout can be derived from a simple dimensional analysis of the convection and diffusion terms in Eq. 1, giving a single dimensionless parameter, the Peclet number (Pe),

washout effect(Pe)= convection termdiffusion term≈Vi⋅∇CiD∇2Ci≈ViCiLDCiL2= ViLD≈4L2JvodD,
where Vi is flow velocity, L is length of crypt or villus, D is diffusion coefficient, Jvo is initial volume flux, and d is crypt/villus diameter. The washout effect is proportional to crypt or villus length and volume flux, and inversely proportional to crypt or villus radius and diffusion coefficient. Substantial washout effect occurs for small crypt/villi radius and long crypt/villi length. The geometric term (L2/d) for villi (L = 350 µm, d = 120 µm from effective diameter) is similar to that for crypts (L = 150 µm, d = 20 µm), which gives a Pe of ∼14 and 16, respectively, for the severe diarrhea case (Jvo = 7 × 10−2 µL/cm2/s). Therefore, both crypts and villi are important determinants of washout, although the washout effect in crypts is greater because of higher fluid secretion rate. In any case, quantitative analysis required full computational analysis because of nonlinear effects and the complex geometry of crypt–villus units.

The main conclusion of the modeling is that the antisecretory efficacy of a surface-targeted channel inhibitor requires both (a) very high luminal inhibitor concentration compared with its Kd, and (b) sustained very high luminal inhibitor concentration or slow inhibitor dissociation compared with the interval between oral dosing. For a crypt secretion rate of 7 × 10−2 µL/cm2/s, typical of cholera, the luminal concentration of a membrane-impermeant inhibitor should be >200 × Kd to reduce net fluid secretion by 50% and 800 × Kd to reduce fluid secretion by 90% (Figs. 3 A and 5 C). Even higher inhibitor concentration is needed if the inhibitor diffusion coefficient is lower than that used in the model here, which may be the case because of mucus and other viscous substances in the aqueous-phase fluid. Further, even if adequate inhibitor concentration in the intestinal lumen is achieved for antisecretory efficacy, a high concentration must be sustained, as delivery of inhibitor-free fluid at the luminal surface results in rapid loss of inhibitor efficacy unless inhibitor dissociation from its luminal target is very slow. For a typical small molecule with low micromolar Kd, dissociation occurs over a few seconds or less, producing significant loss of antisecretory efficacy over minutes (Fig. 6). For oral administration of an inhibitor every few hours, the lumen is exposed to cyclical changes in inhibitor concentration, with little inhibitor present most of the time. Sustained release or other formulations, although potentially helpful, are unlikely to produce sustained therapeutic inhibitor concentrations when the fluid secretion rate is high as in cholera and traveler’s diarrhea.

Of the available externally targeted CFTR inhibitors, GlyH-101 (Sonawane et al., 2006) and small-molecule analogues such as MalH (Sonawane et al., 2007) and iOWH032 (de Hostos et al., 2011) have IC50 in the range of 2–8 µM and dissociate from CFTR within seconds. We conclude from the modeling here that such compounds would not fulfill criteria (a) or (b). MalH-lectin conjugates have higher affinity (Kd of ∼50 nM) and slow washout over several hours or longer (Sonawane et al., 2007), as the lectin moiety is bound tightly and trapped in the dense glycocalyx lining the intestinal surface. However, although MalH-lectin conjugates can thus fulfill criteria (a) and (b), lectin conjugates may not be practical for antisecretory therapy in developing countries because the lectin moiety significantly increases the drug cost and reduces its stability during storage.

The antidiarrheal efficacy of a surface-targeted, membrane-permeable inhibitor that is systemically absorbed and enters the intestinal lumen from the blood side was also considered. Fluid convection can cause washout of inhibitor molecules that enter the crypt lumen from the blood. Typical intestinal transepithelial permeability coefficients for absorbable drugs, based on Caco-2 data, are in the range of 10−6 to 10−5 cm/s, with values up to 7 × 10−5 cm/s for the most highly permeable drugs (Usansky and Sinko, 2005). The computations in Fig. 4 indicate that antidiarrheal efficacy is limited by inhibitor convective washout at high drug Kd ∼ Cb for typical drug permeability (10−6 ∼ 10−5 cm/s). A more potent drug (Kd of ∼0.1 Cb) predicted good antidiarrheal efficacy. Therefore, provided that therapeutic drug concentration in the blood can be maintained, and that drug Kd is sufficiently low and drug permeability sufficiently high, a surface-targeted, membrane-permeable inhibitor can be efficacious in secretory diarrhea.

Although our model was formulated to closely recapitulate the in vivo three-dimensional geometry of a high density of radially oriented crypt–villus units lining the intestinal lumen, several features of the intestine were not included in the model. The model does not consider differences in the detailed geometry and transport properties of the different regions of the small and large intestine, nor does it include fluid absorption or mechanical peristalsis. Although these additional complexities were not modeled explicitly, the major conclusions of our model are robust and should not be influenced by these fine details. The majority of secreted fluid comes from the jejunum, as was modeled here, and fluid secretion in secretory diarrheas such as cholera greatly exceeds absorption. Peristaltic effects would be predicted to have little time-averaged effect on inhibitor diffusion in the unstirred crypt luminal compartment.

Notwithstanding these potential limitations of our model, we conclude that washout in secretory diarrhea is a major factor for nonabsorbable drugs that target surface receptors or transporters on enterocytes in intestinal crypts. Therapeutic efficacy requires sustained high drug concentration in the intestinal lumen. For the inhibitor–CFTR system, surface-targeted small molecules with micromolar Kd are unlikely to be efficacious for therapy of secretory diarrhea, whereas macromolecular conjugates with nanomolar Kd and very slow dissociation might be effective. We conclude, therefore, that although surface-targeted glycine hydrazide and MalH originally appeared to be attractive candidates for antisecretory therapy, absorbable CFTR inhibitors, such as thiazolidinones (Thiagarajah et al., 2004) and PPQ/BPO compounds (Tradtrantip et al., 2009; Snyder et al., 2011), are better development candidates for CFTR inhibition therapy of enterotoxin-mediated secretory diarrheas.

Boundary conditions

As described in Materials and methods, the three-dimensional multi-crypt computation involves solution of the diffusion–convection equation for inhibitor concentration and the (incompressible fluid) Navier–Stokes equation for the fluid field. Fig. A1 and Table 2 provide details on boundary conditions.

Lumen inlet boundary condition

The circular lumen inlet has parabolic velocity profile with mean velocity, Umean. Because the intestinal lumen (∼5-cm diameter) is much bigger than crypt dimensions (20-µm diameter and 150-µm length), the computation was simplified to greatly reduce computation time. A small segment near the mucosa (crypt and villus), with computation volume of 1,250 µm (crypt length of 150 µm, villus length of 350 µm), was used (Fig. 1 C). The velocity profile at the lumen inlet was taken as part of the parabolic profile as shown in Fig. A2. The lumen inlet boundary condition (the magnitude of mean velocity) had negligible effect on percent inhibition (not depicted) because the viscosity in the mucosa is >10-fold greater than that in the intestinal lumen, so that the flow field produced by secretion is dominant in the crypt/villus mucosa. A constant inhibitor concentration boundary condition was used for the inhibitor concentration at the lumen inlet.

Left and right surface boundary condition

As shown in Fig. 5 A, the small segment containing crypts and villi unit is repeated symmetrically around the lumen, producing laminar flow in the lumen (viscous forces dominant over inertia), Re = UmeandL/(ν) < 0.1, where dL is lumen diameter and ν is kinematic viscosity. Thus, symmetric boundary conditions were imposed at the left and the right surface boundaries.

Top surface boundary condition

A slip boundary condition was used for the Navier–Stokes equation, and an insulation (no flux) boundary condition was used for the diffusion–convection equation.

Crypt and villus surface boundary condition

As described in Materials and methods, the rate of fluid secretion at the crypt and villus surface was determined by initial secretion rate and inhibitor concentration. For a membrane-impermeant inhibitor insulation (no flux), boundary conditions were imposed at the crypt and villus surfaces.

Lumen outlet boundary condition

Because the lumen outlet is open, a no-viscous stress and pressure boundary condition was used for the Navier–Stokes equation, and a convective flux boundary condition was used for the diffusion–convection equation.

Supplementary Material Supplemental Material


Abbreviation used in this paper:
MalH malonic acid hydrazides

This work was supported by National Institutes of Health grants DK72517, EB00415, HL73856, DK35124, and EY13574.

Edward N. Pugh Jr. served as editor.

Al-Awqati Q.. Year: 2002Alternative treatment for secretory diarrhea revealed in a new class of CFTR inhibitors. J. Clin. Invest.110:1599–160112464662
Banwell J.G.,Pierce N.F.,Mitra R.C.,Brigham K.L.,Caranasos G.J.,Keimowitz R.I.,Fedson D.S.,Thomas J.,Gorbach S.L.,Sack R.B.,Mondal A.. Year: 1970Intestinal fluid and electrolyte transport in human cholera. J. Clin. Invest.49:183–19510.1172/JCI1062175409804
Bearcroft C.P.,André E.A.,Farthing M.J.. Year: 1997In vivo effects of the 5-HT3 antagonist alosetron on basal and cholera toxin-induced secretion in the human jejunum: a segmental perfusion study. Aliment. Pharmacol. Ther.11:1109–111410.1046/j.1365-2036.1997.d01-1389.x9663837
Chao A.C.,de Sauvage F.J.,Dong Y.J.,Wagner J.A.,Goeddel D.V.,Gardner P.. Year: 1994Activation of intestinal CFTR Cl channel by heat-stable enterotoxin and guanylin via cAMP-dependent protein kinase. EMBO J.13:1065–10727510634
Charmot D.. Year: 2012Non-systemic drugs: a critical review. Curr. Pharm. Des.18:1434–144522300258
Cottreau J.,Tucker A.,Crutchley R.,Garey K.W.. Year: 2012Crofelemer for the treatment of secretory diarrhea. Expert Rev. Gastroenterol. Hepatol.6:17–2310.1586/egh.11.8722149578
de Hostos E.L.,Choy R.K.,Nguyen T.. Year: 2011Developing novel antisecretory drugs to treat infectious diarrhea. Future Med Chem.3:1317–132510.4155/fmc.11.8721859305
Field M.. Year: 1979Mechanisms of action of cholera and Escherichia coli enterotoxins. Am. J. Clin. Nutr.32:189–19632766
Gabriel S.E.,Brigman K.N.,Koller B.H.,Boucher R.C.,Stutts M.J.. Year: 1994Cystic fibrosis heterozygote resistance to cholera toxin in the cystic fibrosis mouse model. Science.266:107–10910.1126/science.75241487524148
Jakab R.L.,Collaco A.M.,Ameen N.A.. Year: 2011Physiological relevance of cell-specific distribution patterns of CFTR, NKCC1, NBCe1, and NHE3 along the crypt-villus axis in the intestine. Am. J. Physiol. Gastrointest. Liver Physiol.300:G82–G9810.1152/ajpgi.00245.201021030607
Loehry C.A.,Creamer B.. Year: 1969Three-dimensional structure of the human small intestinal mucosa in health and disease. Gut.10:6–1210.1136/gut.10.1.65784163
Marsh M.N.,Swift J.A.. Year: 1969A study of the small intestinal mucosa using the scanning electron microscope. Gut.10:940–94910.1136/gut.10.11.9405358588
Muanprasat C.,Sonawane N.D.,Salinas D.,Taddei A.,Galietta L.J.,Verkman A.S.. Year: 2004Discovery of glycine hydrazide pore-occluding CFTR inhibitors: Mechanism, structure–activity analysis, and in vivo efficacy. J. Gen. Physiol.124:125–13710.1085/jgp.20040905915277574
Murek M.,Kopic S.,Geibel J.. Year: 2010Evidence for intestinal chloride secretion. Exp. Physiol.95:471–47810.1113/expphysiol.2009.04944520233891
Naftalin R.J.,Zammit P.S.,Pedley K.C.. Year: 1995Concentration polarization of fluorescent dyes in rat descending colonic crypts: evidence of crypt fluid absorption. J. Physiol.487:479–4958558478
Snyder D.S.,Tradtrantip L.,Yao C.,Kurth M.J.,Verkman A.S.. Year: 2011Potent, metabolically stable benzopyrimido-pyrrolo-oxazine-dione (BPO) CFTR inhibitors for polycystic kidney disease. J. Med. Chem.54:5468–547710.1021/jm200505e21707078
Sonawane N.D.,Hu J.,Muanprasat C.,Verkman A.S.. Year: 2006Luminally active, nonabsorbable CFTR inhibitors as potential therapy to reduce intestinal fluid loss in cholera. FASEB J.20:130–13216317066
Sonawane N.D.,Zhao D.,Zegarra-Moran O.,Galietta L.J.,Verkman A.S.. Year: 2007Lectin conjugates as potent, nonabsorbable CFTR inhibitors for reducing intestinal fluid secretion in cholera. Gastroenterology.132:1234–124410.1053/j.gastro.2007.02.01817408659
Sonawane N.D.,Zhao D.,Zegarra-Moran O.,Galietta L.J.,Verkman A.S.. Year: 2008Nanomolar CFTR inhibition by pore-occluding divalent polyethylene glycol-malonic acid hydrazides. Chem. Biol.15:718–72810.1016/j.chembiol.2008.05.01518635008
Thiagarajah J.R.,Verkman A.S.. Year: 2012CFTR inhibitors for treating diarrheal disease. Clin. Pharmacol. Ther.92:287–29010.1038/clpt.2012.11422850599
Thiagarajah J.R.,Broadbent T.,Hsieh E.,Verkman A.S.. Year: 2004Prevention of toxin-induced intestinal ion and fluid secretion by a small-molecule CFTR inhibitor. Gastroenterology.126:511–51910.1053/j.gastro.2003.11.00514762788
Tradtrantip L.,Sonawane N.D.,Namkung W.,Verkman A.S.. Year: 2009Nanomolar potency pyrimido-pyrrolo-quinoxalinedione CFTR inhibitor reduces cyst size in a polycystic kidney disease model. J. Med. Chem.52:6447–645510.1021/jm900987319785436
Tradtrantip L.,Namkung W.,Verkman A.S.. Year: 2010Crofelemer, an antisecretory antidiarrheal proanthocyanidin oligomer extracted from Croton lechleri, targets two distinct intestinal chloride channels. Mol. Pharmacol.77:69–7810.1124/mol.109.06105119808995
Trbojević-Stanković J.B.,Milićević N.M.,Milosević D.P.,Despotović N.,Davidović M.,Erceg P.,Bojić B.,Bojić D.,Svorcan P.,Protić M.,et al. Year: 2010Morphometric study of healthy jejunal and ileal mucosa in adult and aged subjects. Histol. Histopathol.25:153–15820017102
Usansky H.H.,Sinko P.J.. Year: 2005Estimating human drug oral absorption kinetics from Caco-2 permeability using an absorption-disposition model: model development and evaluation and derivation of analytical solutions for ka and Fa. J. Pharmacol. Exp. Ther.314:391–39910.1124/jpet.104.07618215833900
Van Loon F.P.,Rabbani G.H.,Bukhave K.,Rask-Madsen J.. Year: 1992Indomethacin decreases jejunal fluid secretion in addition to luminal release of prostaglandin E2 in patients with acute cholera. Gut.33:643–64510.1136/gut.33.5.6431612480
Venkatasubramanian J.,Ao M.,Rao M.C.. Year: 2010Ion transport in the small intestine. Curr. Opin. Gastroenterol.26:123–12810.1097/MOG.0b013e3283358a4520010100
Yang B.,Sonawane N.D.,Zhao D.,Somlo S.,Verkman A.S.. Year: 2008Small-molecule CFTR inhibitors slow cyst growth in polycystic kidney disease. J. Am. Soc. Nephrol.19:1300–131010.1681/ASN.200707082818385427
Zhang W.,Fujii N.,Naren A.P.. Year: 2012Recent advances and new perspectives in targeting CFTR for therapy of cystic fibrosis and enterotoxin-induced secretory diarrheas. Future Med Chem.4:329–34510.4155/fmc.12.122393940


[Figure ID: fig1]
Figure 1. 

Model of drug convection–diffusion in intestinal crypt–villus units. (A) Schematic of epithelial cell–lined crypt–villus units (lengths, Lcrypt and Lvillus). Fluid secretion into the lumen produces convective (upward) solute transport. (B) Three-dimensional model of radially oriented crypt–villus units projecting from a cylindrical lumen. (C) Section near the intestinal wall showing crypt–villus units and a wedge of luminal aqueous-phase volume. Cylindrical symmetry reduces computation complexity. An example of pseudocolored inhibitor concentration profile is shown.

[Figure ID: fig2]
Figure 2. 

Convective inhibitor washout in single crypts. Inhibitor convection–diffusion was solved for a cylindrical crypt–villus unit exposed at its outer boundary (intestinal lumen) to a constant concentration of membrane-impermeant inhibitor, Co/Kd. (A) Cross section of a cylindrical crypt–villus unit showing relative Jvo at indicated regions (left). Two-dimensional computation volume showing mesh elements (middle). Equivalent three-dimensional crypt–villus unit showing downward inhibitor diffusion and upward convection. (B) Steady-state profiles of drug concentration (left), fluid flow velocity (middle), and pressure (right) for single-crypt diffusion–convection for indicated Co/Kd, with Jvo = 4 × 10−2 µL/cm2/s. (C) Steady-state profiles for indicated Jvo, with Co/Kd = 1.

[Figure ID: fig3]
Figure 3. 

Convective inhibitor washout requires a very high concentration of a membrane-impermeant inhibitor in the intestinal lumen for antisecretory efficacy. (A) Computations done for human mid-jejunal anatomy as in Fig. 2. Percent inhibition of net secreted fluid as a function of Co/Kd for indicated Jvo. Jvo of ∼7 × 10−2 µL/cm2/s is typical in cholera. (B) Computations done for human colonic anatomy with long, narrow crypts without villi. Percent inhibition of net secreted fluid as a function of Co/Kd for indicated Jvo. Jvo of ∼7 × 10−3 µL/cm2/s in colon is typical in cholera.

[Figure ID: fig4]
Figure 4. 

Model of antidiarrheal efficacy of a membrane-permeant inhibitor. Computations done for human mid-jejunal crypt–villus anatomy. (A) Schematic showing inhibitor permeation across enterocytes from blood into the crypt–villus lumen, which depends on inhibitor transepithelial permeability coefficient, Pinh, and blood concentration Cb. (B) Steady-state inhibitor concentration profiles for indicated Pinh for Jvo = 2 × 10−2 µL/cm2/s (left), and for indicated Jvo for Pinh = 10−5 cm/s (right), both for Cb = 10 µM and Kd = 10 µM. (C) Percent inhibition of net secreted fluid as a function of Pinh for indicated Jvo at Cb = 10 µM and Kd = 10 µM. (D) Same as in C, with Cb = 1 µM and Kd = 0.1 µM.

[Figure ID: fig5]
Figure 5. 

Inhibitor washout in a three-dimensional multi-crypt model of the intestine. (A) Crypt–villus geometry in the three-dimensional model showing a matrix of spatially distinct narrow crypts and relatively wide villi. (B) Drug concentration profiles for steady-state solution of multi-crypt diffusion–convection for indicated Jvo for Co/Kd = 10. (C) Percent inhibition of net secreted fluid as a function of Co/Kd for indicated Jvo. Jvo of ∼7 × 10−2 µL/cm2/s is typical in cholera. (D) Percent inhibition along the intestine for indicated lumen axial mean velocity, Umean. Segmental percent inhibition in C was described by an empirical regression to the equation: % inhibition = α · exp[β · log10(Cin/Kd)], with α = 3.523 and β = 1.112 for Jvo = 7 × 10−2 µL/cm2/s, typical of cholera. (E) Percent inhibition for 5-m-long intestinal segment as a function of Co/Kd for indicated Umean. Results for a short (1-mm) intestinal segment are shown for comparison.

[Figure ID: fig6]
Figure 6. 

Rapid washout of a luminally delivered membrane-impermeant inhibitor. After steady-state inhibition, inhibitor concentration in the intestinal lumen was reduced from Co/Kd to 0. (A) Drug concentrations profiles for kinetic (presteady-state) solution of single-crypt diffusion–convection for Co/Kd = 10, Jvo = 10−3 µL/cm2/s, and k−1 = 10−1 min−1. (B) Kinetics of inhibition of fluid secretion after inhibitor washout. Percent inhibition of net fluid secretion as a function of time after washout for Co/Kd = 10, Jvo = 10−3 µL/cm2/s, and indicated k−1. (C) Washout kinetics for the three-dimensional multi-crypt model for same parameters as in A.

[Figure ID: fig7]
Figure A1. 

Schematic of boundary conditions for the three-dimensional multi-crypt computation. (A) Boundary conditions for the Navier–Stokes equation. (B) Boundary conditions for the diffusion–convection equation.

[Figure ID: fig8]
Figure A2. 

Schematic of inlet boundary condition, where s is distance from the center of the lumen, so is distance from the center of the lumen to the top of the villus, and Umean is the mean inhibitor velocity in the lumen.

[TableWrap ID: tbl1] Table 1. 

Material properties and baseline parameters

Material property Baseline parameter
ρ Fluid density = 103 (kg/m3)
DL Diffusion coefficient in the intestinal lumen = 2 × 10−9 (m2/s)
DM Diffusion coefficient in the mucosa (crypt/villus region) = 2 × 10−10 (m2/s)
ηL Fluid dynamic viscosity in the lumen = 10−3 (Pa · s)
ηM Fluid dynamic viscosity in the mucosa = 10−2 (Pa · s)
LV Villus length = 350 (µm)
LC Crypt length = 150 (µm)
dV Villus diameter = 120 (µm)
dC Crypt diameter = 20 (µm)
dL Intestinal lumen diameter = 5 (cm)

[TableWrap ID: tbl2] Table 2. 

Boundary conditions for Navier–Stokes and diffusion–convection equations

Navier–Stokes equation Diffusion–convection equation
(1) Inlet BC U = 2Umean · [1−(s/so)2] (1) Inlet BC C = Co
(lumen flow velocity) Vi = −U·n (constant Co)
(2) Symmetric BC n·Vi = 0, t · [−pI + µ(∇Vi + (∇Vi)T] = 0 (2) Symmetric BC n · n = 0, N = −D · ∇C + CVi
(3) Slip BC n·Vi = 0, t · [−pI + µ(∇Vi + (∇Vi)T] = 0 (3) Insulation BC n · n = 0, N = −D · ∇C + CVi
(4) Flux BC Vi = -Jv · n, Jv = Jvo(1 − Ci/(Ci + Kd)) (4) Insulation BC n · n = 0, N = −D · ∇C + CVi
(5) Outlet BC [µ(∇Vi + (∇Vi)T]n = 0, p = p0 (5) Outlet BC n·(−D∇C) = 0
(no viscous stress) (convective flux)

BC, boundary condition; n, surface normal vector; I, unit vector; ∇, gradient operator. Other variables are defined in the main text.

Article Categories:
  • Article

Previous Document:  Interaction between residues in the Mg2+-binding site regulates BK channel activation.
Next Document:  Effects of PFM rehabilitation on PFM function and morphology in older women.