Front. Cell. Neurosci.

Frontiers in Cellular Neuroscience

Front. Cell. Neurosci.

1662-5102

Frontiers Media S.A.

10.3389/fncel.2019.00161

Cellular Neuroscience

Original Research

Drug Diffusion Along an Intact Mammalian Cochlea

Sadreev

Ildar I.

¹ ^* Burwood

George W. S.

² ^† Flaherty

Samuel M.

² Kim

Jongrae

³ Russell

Ian J.

² Abdullin

Timur I.

⁴ Lukashkin

Andrei N.

² ⁵ ^*

¹Department of Medicine, Faculty of Medicine, Imperial College, London, United Kingdom ²Sensory Neuroscience Research Group, School of Pharmacy and Biomolecular Sciences, University of Brighton, Brighton, United Kingdom ³School of Mechanical Engineering, Institute of Design, Robotics and Optimisation, Aerospace Systems Engineering, University of Leeds, Leeds, United Kingdom ⁴Department of Biochemistry, Biotechnology and Pharmacology, Institute of Fundamental Medicine and Biology, Kazan Federal University, Kazan, Russia ⁵Centre for Regenerative Medicine and Devices, University of Brighton, Brighton, United Kingdom

Edited by: Sylvain Celanire, Pragma Therapeutics, France

Reviewed by: Gunnar P. H. Dietz, University of Göttingen, Germany; Alec Nicholas Salt, Washington University in St. Louis, United States

*Correspondence: Ildar I. Sadreev i.sadreev@imperial.ac.uk Andrei N. Lukashkin a.lukashkin@brighton.ac.uk

This article was submitted to Cellular Neurophysiology, a section of the journal Frontiers in Cellular Neuroscience

†Present Address: George W. S. Burwood, Oregon Hearing Research Center (OHSU), Portland, OR, United States

26 04 2019

2019

161

28 01 2019 08 04 2019

2019

Sadreev, Burwood, Flaherty, Kim, Russell, Abdullin and Lukashkin

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.

Intratympanic drug administration depends on the ability of drugs to pass through the round window membrane (RW) at the base of the cochlea and diffuse from this location to the apex. While the RW permeability for many different drugs can be promoted, passive diffusion along the narrowing spiral of the cochlea is limited. Earlier measurements of the distribution of marker ions, corticosteroids, and antibiotics demonstrated that the concentration of substances applied to the RW was two to three orders of magnitude higher in the base compared to the apex. The measurements, however, involved perforating the cochlear bony wall and, in some cases, sampling perilymph. These manipulations can change the flow rate of perilymph and lead to intake of perilymph through the cochlear aqueduct, thereby disguising concentration gradients of the delivered substances. In this study, the suppressive effect of salicylate on cochlear amplification via block of the outer hair cell (OHC) somatic motility was utilized to assess salicylate diffusion along an intact guinea pig cochlea in vivo. Salicylate solution was applied to the RW and threshold elevation of auditory nerve responses was measured at different times and frequencies after application. Resultant concentrations of salicylate along the cochlea were calculated by fitting the experimental data using a mathematical model of the diffusion and clearing of salicylate in a tube of variable diameter combined with a model describing salicylate action on cochlear amplification. Concentrations reach a steady-state at different times for different cochlear locations and it takes longer to reach the steady-state at more apical locations. Even at the steady-state, the predicted concentration at the apex is negligible. Model predictions for the geometry of the longer human cochlea show even higher differences in the steady-state concentrations of the drugs between cochlear base and apex. Our findings confirm conclusions that achieving therapeutic drug concentrations throughout the entire cochlear duct is hardly possible when the drugs are applied to the RW and are distributed via passive diffusion. Assisted methods of drug delivery are needed to reach a more uniform distribution of drugs along the cochlea.

cochlea drug delivery salicylate cochlear amplifier cochlear round window

MR/N004299/1

Medical Research Council10.13039/501100000265

�㾩julia��߲��

Introduction

The mammalian cochlea is one of the least accessible organs for drug delivery (Salt and Plontke, 2009; Rivera et al., 2012; El Kechai et al., 2015; Hao and Li, 2019). Systemic administration of many drugs, notably the most frequently used corticosteroids, and aminoglycoside antibiotics, is severely limited by the blood-labyrinth barrier (Salt and Hirose, 2018). Local intratympanic administration (Schuknecht, 1956; Bowe and Jacob, 2010) would be a preferable option for these drugs and local delivery is the only option for many old and newly emerging classes of drugs and therapies including local anesthetics, antioxidants, apoptosis inhibitors, neurotransmitters and their antagonists, monoclonal antibodies, growth factors, signaling pathway regulators, and genetic material (see Devare et al., 2018; Hao and Li, 2019 for the latest reviews). Intratympanic administration of drugs relies on their remaining in contact with the round window membrane (RW) (a membranous opening in the bony wall of the cochlea into the middle ear) long enough to allow their diffusion into the perilymph of the scala tympani (ST). The ability of drugs to pass through the RW does not, however, guarantee their sufficient distribution along the cochlear spiral. Drug distribution in the ST is limited by the low flow rate of perilymph within the cochlea and by cochlear geometry. The longitudinal flow of perilymph in the cochlea has been shown to be relatively slow, if present at all (Ohyama et al., 1988), and drug distribution in the perilymph is dominated by passive diffusion. Passive diffusion along the ST is, however, constrained because the cochlea is a relatively long and narrow tube with a cochlear cross-section that decreases gradually from the RW at the base to the apex. It is in the cochlear apex where human speech processing is initiated (e.g., Nuttall et al., 2018) and where drug delivery to the cochlea has greatest potential therapeutic and socioeconomic impact.

However, direct measurements of the distribution of marker ions and contrasting agents (Salt and Ma, 2001; Haghpanahi et al., 2013), corticosteroids (Plontke et al., 2008; Creber et al., 2018) and antibiotics (Mynatt et al., 2006; Plontke et al., 2007a) or measurements of the physiological effects of drugs (Chen et al., 2005; Borkholder et al., 2010) demonstrated that the concentration of substances applied to the RW was much higher in the cochlear base than in the apex. These measurements, however, involved perforating the cochlear bony wall and, in some cases, sampling perilymph. These manipulations can change the flow rate of perilymph (Ohyama et al., 1988; Salt and Ma, 2001) and lead to the intake of cerebrospinal fluid through the cochlear aqueduct (Salt et al., 2003), thereby disguising concentration gradients of the delivered substances.

A few studies investigated the distribution of substances applied to the RW in the intact cochlea without breaking cochlear boundaries. This was done mainly in morphological studies investigating the distribution of dexamethasone and other substances along the cochlea after their intratympanic administration (Saijo and Kimura, 1984; Imamura and Adams, 2003; Hargunani et al., 2006; Grewal et al., 2013). While these studies confirmed the existence of base-to-apex gradients, the actual concentrations of substances along the cochlea were not measured. Borkholder et al. (2014) measured the threshold elevation of distortion product otoacoustic emissions (DPOAE) produced by primary tones of different frequencies after intratympanic application of salicylate. Salicylate affects cochlear amplification in a concentration-dependent manner but the DPOAE is a non-linear phenomenon and the dependence of DPOAE thresholds on the primary tone level and cochlear amplification is complex (Lukashkin et al., 2002). As a result, salicylate concentrations along the cochlear spiral cannot be easily derived from the DPOAE threshold elevations.

The purpose of the current study is to quantify drug diffusion from the RW along an intact guinea pig cochlea, to identify the factors that limit passive drug diffusion along the cochlea, and to analyse possible solutions to overcome these limitations. Salicylate was used as a model drug with well-characterized physiological effects. A mathematical model, which includes a diffusion component and a biophysical component describing the action of salicylate on the cochlear amplifier was validated using experimental data and used to assess the distribution of substances along the human cochlea.

Materials and Methods Animals

Pigmented guinea pigs of similar weight (350–360 g) were anesthetized with the neurolept anesthetic technique (0.06 mg/kg body weight atropine sulfate s.c., 30 mg/kg pentobarbitone i.p., 500 μl/kg Hypnorm i.m.). Additional injections of Hypnorm were given every 40 min. Additional doses of pentobarbitone were administered as needed to maintain a non-reflexive state. The heart rate was monitored with a pair of skin electrodes placed on both sides of the thorax. The animals were tracheotomized and artificially respired, and their core temperature was maintained at 38°C with a heating blanket and a heated head holder. All procedures involving animals were performed in accordance with UK Home Office regulations with approval from the local ethics committee.

Signal Generation and Recording

The middle ear cavity of the ear used for the measurements was opened to reveal the RW. Compound action potentials (CAPs) of the auditory nerve were measured from the cochlear bony ridge in the proximity of the RW membrane using Teflon-coated silver wire coupled to laboratory designed and built extracellular amplifier (James Hartley). Thresholds of the N1 peak of the CAP were estimated visually using 10 ms tone stimuli at a repetition rate of 10 Hz.

For acoustic stimulation, sound was delivered to the tympanic membrane by a closed acoustic system comprising two Bruel and Kjaer 4134 ½″ microphones for delivering tones and a single Bruel and Kjaer 4133 ½″ microphone for monitoring sound pressure at the tympanum. The microphones were coupled to the ear canal via 1 cm long, 4 mm diameter tubes to a conical speculum, the 1 mm diameter opening of which was placed about 1 mm from the tympanum. The closed sound system was calibrated in situ for frequencies between 1 and 50 kHz. Known sound pressure levels were expressed in dB SPL re 2 × 10⁻⁵ Pa.

All acoustic stimuli in this work were shaped with raised cosines of 0.5 ms duration at the beginning and at the end of stimulation. White noise for acoustical calibration and tone sequences for auditory stimulation were synthesized by a Data Translation 3010 board at 250 kHz and delivered to the microphones through low-pass filters (100 kHz cut-off frequency). Signals from the acoustic measuring amplifier (James Hartley) were digitized at 250 kHz using the same board and averaged in the time domain. Experimental control, data acquisition, and data analysis were performed using a PC with programmes written in MATLAB (The MathWorks. Inc. 2018a).

Five microliters of sodium salicylate solution (either 100 mM in experiments on salicylate diffusion in the ST or 1M in experiments with complete block of the cochlear amplifier) in Hanks' Balanced Salt Solution were placed on the RW using pipettes. The solution was removed from the RW using paper wicks to observe the wash out effect.

Model Overview Diffusion and Clearing Equation

For the purpose of modeling, the ST is approximated by a tube with a decreasing diameter similar to that described in previous models, for example by Plontke et al. (2007b) (Figure 1A). The radii of the tube, r(0) and r(l), are equal to a and b at x = 0 and x = l, respectively, where l is the ST length. All the dimensions are known (Thorne et al., 1999) and symmetry along y and z axes is assumed. Zero longitudinal perilymph flow in the compartment is assumed (Ohyama et al., 1988) and only the passive diffusion of a drug (salicylate) with diffusion coefficient k_d is considered. In addition to diffusion, there is also clearing of the drug characterized by the clearing coefficient k_c. This clearing can be represented simply as a leak through the scala boundary (e.g., loss to the vasculature and tissues, and to other cochlear compartments). The diffusion and clearing processes are assumed to be completely independent. Because the tube radius is much smaller than its length, i.e., r(x) ≪ l for all x in [0, l], only diffusion along x axis is considered and the concentration c(x, t) within each cross-section for a fixed instance t is assumed to be constant, i.e., it does not change along the y axis. If the area of the cross-section is S(x) then the diffusion can be described by the following partial differential equation (see Appendix for detailed derivation):

(1)dc(x,t)dt=1S(x)·ddx(S(x)·kd·dc(x,t)dx)-c(x,t)·2kcr(x),

with the boundary conditions

(2)c(0,t)=crw,

(3)kddc(l,t)dx=0

and initial conditions

(4)c(0,0)=crw;x=0,

(5)c(x,0)=0;x>0.

The diffusion coefficient k_d is known (Lide, 2002) but the clearing coefficient k_c is unknown. The ratio of the diffusion and clearing coefficients can, however, be found via fitting the experimental data. The physical meaning of k_d/k_c can be described as the ratio between the amount of substance that diffuses through a unit surface normal to the direction of diffusion for a unit concentration gradient and the amount of drug that is cleared through a unit surface normal to the direction of substance exit for a unit substance concentration, both for unit time duration. The diffusion/clearing equation was validated using experimental data on the physiological effect of salicylate on the CAP thresholds. Because the salicylate concentrations could not be directly inferred from the physiological effect of salicylate, a biophysical element of the model was developed allowing calculations of the salicylate concentrations along the cochlea.

Figure 1

Schematic presentation of (A) the scala tympani (ST) approximated by a tube of decreasing diameter and (B) the cochlear amplifier modeled as a system with positive feedback provided by the outer hair cells (OHCs).

Link Between Position and Frequency

The dependence between frequency of stimulation f and frequency position along the length x of the basilar membrane for the guinea pig cochlea is defined by the Greenwood equation (Greenwood, 1990)

(6)f(x~)=A·(10αx~-β),

where A = 0.35, α = 2.1/18.5, β = 0.85 and x~=l-x meaning that the starting point for x~ in Greenwood (1990) is at the apex and not the base of the cochlea, as in this study.

Cochlear Amplifier

The cochlear amplifier is represented by a positive feedback system (Figure 1B) with feedback gain H(x, c(t)) due to force generation by the OHCs (Mountain et al., 1983; Yates, 1990; Lukashkin and Russell, 1999). The following assumptions are made for a small signal, linear regime:

The CAP threshold is observed for different sound pressure P_in(x, c(t)) at the tympanum but for the same BM displacements, i.e., for the same constant pressure P_out(x) at the BM for any given frequency/place x during manipulations with the cochlear amplifier. The assumption is based on good correspondence between neural and BM thresholds at the CF (Ruggero et al., 2000; Temchin et al., 2008).

Feedback gain H(x, c(t)) is proportional to the outer hair cell (OHC) force F(x, c(t)) for any given frequency/place in the cochlea (7)H(x,c(t))=α·F(x,c(t)),

where α is the gain constant. The initial feedback gain H(x, 0) for any frequency/place before application of salicylate can be found empirically (see below).

Salicylate changes only feedback gain H(x, c(t)) through changes in F(x, c(t)).

In line with other modeling studies (e.g., Meaud and Grosh, 2014; Ni et al., 2016), it is assumed that pressure/displacement at the BM is a linear combination of the passive BM response due to acoustic stimulation and active response due to the OHC forces.

The link between local salicylate concentration c(x, t) and reduction in force FR(x, c(t)) generated by the OHCs can be described by the Hill function (Hallworth, 1997)

(8)FR(x,c(t))=Vmax·c(x,t)nkn+c(x,t)n,

where V_max = 0.71629, k = 0.101, and n = 0.983.

The reduction in force is linked to the force before F(x, 0) and after F(x, c(t)) salicylate application as

(9)FR(x,c(t))=F(x,0)-F(x,c(t))F(x,0)=1-F(x,c(t))F(x,0)

or,

(10)F(x,c(t))=F(x,0)·(1-FR(x,c(t))).

It can be written for any given frequency/place before salicylate application at t = 0 (Figure 1B)

(11)Pout(x)Pin(x,0)=G1-G·H(x,0),

where G is the open loop gain. Similarly, at time t after salicylate application

(12)Pout(x)Pin(x,c(t))=G1-G·H(x,c(t)).

Dividing (11) by (12) and taking into account (7), it could be written

(13)Pin(x,c(t))Pin(x,0)=1-G·α·F(x,c(t))1-G·α·F(x,0).

Substituting F(x, c(t)) from (10) into (13) and using (7), one can obtain

(14)Pin(x,c(t))Pin(x,0)=1-G·H(x,0)·(1-FR(x,c(t)))1-G·H(x,0).

The left part of (14) is measured in the experiment. FR(x, c(t)) is calculated using the Hill function (8) with c(x, t) in this equation being calculated using the diffusion/clearing equation (1).

An analytical form of empirical dependence H(x, 0), i.e., feedback gain before salicylate application for different frequencies/locations, can be obtained as follows. Feedback from the OHCs can be completely blocked in experiments using a high concentration of salicylate. In this case H(x, c(t)) = 0 in (12) and the transfer function of the feedback system (Figure 1B) is equal to the open loop gain G. Then similar to (11) and (12)

(15)Pout(x)PinBlock(x)=G,

where P_inBlock(x) is the sound pressure required to produce a response from the auditory nerve in preparations where the cochlear amplifier is completely blocked, and it does not depend on time. Dividing (11) by (15) and rearranging gives the following equation

(16)H(x,0)=1G·(1-Pin(x,0)PinBlock(x)),

where P_inBlock(x)/P_in(x, 0) is measured in separate experiments.

G has frequently been assumed to be constant along the cochlea (e.g., Mountain et al., 1983; Yates, 1990; Lukashkin and Russell, 1999). In spite of the special design of the cochlea, which minimized energy losses when the BM traveling wave moves from the base to apex (Jones et al., 2013), some energy dissipation is still expected during wave propagation in a viscous environment. To account for energy losses, we assumed a simple linear dependence of the open loop gain G(f(x)) on frequency

(17)G(f(x))=s·f(x)+i,

where s is the slope and i is the intercept defined as i = 1 − s · f_l, with f_l = 49.9165 kHz specifying the upper frequency limit of linear dependence for G(f(x)). Hence, G(f(x)) effectively depends only on a single parameter s, which could be found by fitting the experimental data.

Initial Model Parameters

Ratio P_inBlock(x)/P_in(x, 0) was measured as a function of frequency f. Equation (6) shows how this frequency can be converted to a coordinate. An arbitrary Hill type function

(18)PinBlockPin(f)=m1fm2m3m2+fm2+m4

was fitted to the experimental data with 20log₁₀ transformation for dB using the Genetic Algorithm (GA) tool in MATLAB (The MathWorks. Inc. 2018a) (initial local fit). The obtained m₁, m₂, m₃, and m₄ (Table 1) were then used for the later optimization procedures described below (final global fit). The feedback gain before application of salicylate H(x, 0) = H(f(x), 0) was obtained according to (16) and (18) as

(19)H(f,0)=1G·(1-1/(m1fm2m3m2+fm2+m4)).

Initial values for the all parameters used in the model before the optimization procedure are shown in Table 1.

Table 1

Model parameter values.

Parameter	Unit	Initial value	Optimized value	Source of the initial value
a	mm	0.56	Fixed	Thorne et al., 1999
b	mm	0.18	Fixed	Thorne et al., 1999
l	mm	18–19	19	Thorne et al., 1999
k_d	mm²/s	0.959e-3	Fixed	Lide, 2002
ratio = k_d/k_c	mm	1–10	1.6968	Initial guess
c_rw	mM	100	Fixed	Experiment
A	kHz	0.35	Fixed	Greenwood, 1990
α	1/mm	2.1/18.5	Fixed	Greenwood, 1990
β	–	0.85	Fixed	Greenwood, 1990
k	mM	0.101	Fixed	Hallworth, 1997
n	–	0.983	Fixed	Hallworth, 1997
V_max	–	0.71629	Fixed	Hallworth, 1997
m₁	–	1011.2	Fixed	Experiment
m₂	–	8.1406	Fixed	Experiment
m₃	kHz	3.4816	Fixed	Experiment
m₄	–	31.686	Fixed	Experiment
s	1/kHz	0–0.0261	0.00014742	Initial guess

Optimized Model Parameters

Equation (14) with 20log₁₀ transformation for dB was solved in MATLAB (The MathWorks. Inc. 2018a) using pdepe solver for partial differential equations and fitted to the entire set of experimental data for all frequencies and salicylate concentrations using Genetic Algorithm (GA) tool in MATLAB (The MathWorks. Inc. 2018a). The sum of squared errors

SE=∑i=1n(Mi-Ei)2

was used as a cost function for minimization, where M are the model predictions and E are the experimental data for points i = 1…n. It is worth noting that only three model parameters were fitted during the global fit/optimization. These parameters are the cochlear length l, k_d/k_c ratio and slope s of the open loop gain G(f(x)).

Results Cochlear Amplifier Gain

Gain of the cochlear amplifier and corresponding feedback gain of the model, H(x, 0) = H(f(x), 0) [equation (16)] was determined empirically from elevation of the CAP thresholds after application of 1M salicylate solution to the RW which caused a consistent and steady increase in threshold over the entire frequency range (Figure 2, black circles). Values for m₁, m₂, m₃, and m₄ (Table 1) were determined through fit of the experimental data points by equation (18) (Figure 2, red curve) using the Genetic Algorithm (GA) tool in MATLAB (The MathWorks. Inc. 2018a). These values were used for the general optimization procedure performed at later stages.

Figure 2

Elevation of CAP thresholds after complete block of the cochlear amplifier (left Y-axis) and corresponding value of the open loop and feedback gain (right Y-axis). Black circles show the experimental values of threshold elevation (mean ± SD, n = 3). Red curve indicates fit of the experimental data points by Equation (18). Related values of the parameters m₁, m₂, m₃ and m₄ are given in Table 1. Value of the open loop (magenta curve) and feedback (blue curve) gains after the final global optimization procedure were calculated using Equations (17) and (19), respectively, with the optimized value of parameter s (Table 1).

Distribution of Salicylate Along the Guinea Pig Cochlea

One hundred mM solution of salicylate applied to the RW caused a rapid increase followed by saturation of CAP thresholds for high frequency tones (Figure 3A). CAP threshold increase for tones of lower frequencies was observed after an initial delay and did not reach saturation during the time of observation. Any changes in CAP threshold due to application of salicylate were below the noise floor of measurements for tone frequencies lower than 5 kHz, which corresponds to approximately the apical 55% of cochlear length [(Greenwood, 1990); equation (6)]. A partial recovery of the CAP threshold was observed after salicylate solution was washed out from the RW confirming that threshold elevation during application of salicylate was due to specific action of salicylate and not because of general deterioration of preparations.

Figure 3

CAP threshold elevation (A,B) and salicylate distribution (C) in the guinea pig cochlea after application of 100 mM of salicylate solution to the RW at time = 0. (A) Representative example of the CAP threshold elevation in a single preparation. Salicylate was washed out after 80 min of application. (B) Combined best fit of the entire set of experimental data on CAP threshold elevation for five preparations (Figure S1) using the parameter optimization procedure. Labels indicate percentage of the total cochlear length from the base. (C) Salicylate concentration along the cochlear length calculated using the optimized values of the model parameters (Table 1).

Pooled data from five animals were used to find an optimized set of the model parameters via fitting the entire set of experimental data using the Genetic Algorithm (GA) tool in MATLAB (The MathWorks. Inc. 2018a) (see Materials and Methods). The combined best fit to the entire set of experimental data for the optimized set of parameters is illustrated in Figure 3B. Figure S1 shows the same plots for each of the frequencies along the corresponding experimental data. It is worth noting, that the optimization procedure was performed over the entire experimental set in order to fit the data for all the experimental frequencies simultaneously (Figure S1). The optimized set of model parameters (Table 1) found due to the general optimization procedure was used to predict cochlear responses and concentrations of salicylate (Figure 3C) along the entire cochlear length and over arbitrary time duration.

The absence of CAP threshold changes at frequencies below 5 kHz was due to poor diffusion of salicylate from the RW into the cochlear apex. It required increasingly longer times for the salicylate concentration to reach steady-state in the more apical regions of the cochlea, but at 90% of cochlear length (10% from the apex), salicylate concentration was about 12 orders of magnitude smaller than at the base even at steady-state (Figure 3C). The model suggests that this steep concentration gradient is due mainly to the fast clearing of salicylate from the ST which is reflected in the small k_d/k_c ratio found in the optimization (Table 1). Because the flux J is proportional to the concentration gradient [equation (A1)], changes in salicylate concentration at the RW will not lead to changes in the concentration gradient between the cochlear base and apex. In this case all steady-state curves for different concentrations of salicylate at the RW are scaled versions of each other (data not shown). Hence, for a specific substance (i.e., for specific diffusion (k_d) and clearing (k_c) coefficients) and for a given cochlear geometry, the ratio of steady-state concentrations at the base, and apex of the cochlea is a constant and does not depend on substance concentration at the RW. This was further assessed for the human cochlea.

Diffusion of an Arbitrary Substance in the Human Cochlea

Hence, the validity of the diffusion/clearing equation has been confirmed using the experimental data on salicylate block of the cochlear amplifier, the equation can be used to make conclusions about the distribution of arbitrary substances along the human cochlea (Figure 4). Decrease in the relative contribution of clearing into the distribution of a substance along the ST, i.e.,. increase of k_d/k_c ratio, leads to a dramatic reduction in the steady-state, base-to-apex gradient of the substance concentration (Figure 4A) calculated using the non-dimensional form of the diffusion equation (A14). This result is expected because a larger amount of the substance is available for diffusion into the cochlear apex in this case. For salicylate, however, the difference between the basal and apical concentrations is even larger in the human cochlea (red cross in Figure 4A) compared to guinea pigs and reaches 16 orders of magnitude because of the increased length of the human cochlea.

Figure 4

Theoretical distribution of an arbitrary substance in the human cochlea. (A) Dependence of the ratio of basal (c10%) and apical (c90%) steady-state concentrations on the ratio of diffusion (k_d) and clearing (k_c) coefficients. Red cross indicates the point for salicylate (k_d/k_c = 1.7). (B) Normalized time (black curve, left ordinate) required to reach steady-state concentration at the cochlear apex for substances with different ratio of the diffusion (k_d) and clearing (k_c) coefficients. Red curve shows a specific example of the absolute time (right ordinate) for a substance with the diffusion coefficient (k_d) similar to that of dexamethasone. The steady-state was defined as the normalized difference between consecutive numerical values of concentration <10⁻⁴. Jitter in the curves for small k_d/k_c ratios is due to very low apical concentrations. The following geometrical parameters for human cochlea were used for calculations using the non-dimensional form of the diffusion equation (A14) a = 0.7981, b = 0.3990, and l = 28.46 mm (Thorne et al., 1999).

Figure 4A provides theoretical estimates of the minimal gradients which can be reached along the ST due to passive diffusion, when substances are in contact with the RW long enough to establish a concentration equilibrium distribution. Reduction in the base-to-apex gradient for substances with higher k_d/k_c ratios, which are better retained in the ST, comes at the expense of the much longer substance exposure times required to reach steady-state concentration gradients (Figure 4B). For example, for a drug with the diffusion coefficient k_d similar to dexamethasone, for which the clearing coefficient is unknown, it takes days of retention at the RW when realistic k_d/k_c ratios are assumed (red curve in Figure 4B). The problem is that, while it is theoretically possible to achieve smaller base-to apex concentration gradients for a drug with high k_d/k_c ratios, in practice, if the drug is active, it will be cleared from the ST into the cochlear tissue, hence k_c cannot be arbitrarily small. In this case, the minimal theoretical difference in the base-to-apex concentrations of the drug is still a few orders of magnitude.

Discussion

The existence of a base-to-apex drug concentration gradient, when drugs are applied to the RW, has been well-established. From this point of view, this study quantifies these gradients for the intact cochlea when the flow of perilymph in the ST is very small (Ohyama et al., 1988). This study does not investigate the problem of the RW permeability which is a separate challenge and requires specific considerations for particular drugs and formulations (Salt and Plontke, 2018). Instead, sodium salicylate which easily passes through the RW was used to ensure high concentrations at the cochlear base. Though, passive proton-mediated diffusion of salicylate across biomembranes is observed at micromolar concentrations (Takagi et al., 1998), the RW diffusional barrier could, presumably, be overcome by the drug at the much higher, submolar concentrations used in this study. While the RW membrane is highly permeable to salicylate and the CAP threshold elevation at high frequencies started within seconds after salicylate application, the model assumption that salicylate concentrations on both RW sides were the same might introduce some error in the calculated absolute concentrations. We would like to emphasize, however, that an error in calculation of the absolute concentrations (note that the absolute concentrations were calculated using Hallworth's (1997) empirical dependence between salicylate concentration and the OHC force reduction) does not lead to an error in calculation of the concentration gradient which is the basis for the conclusions in this study. This is true because the flux J is proportional to the concentration gradient [equation (A1)] and gradient curves calculated for different salicylate concentrations at the RW are scaled versions of each other with the same gradients. It is worth noting that, from discoveries we made in our preliminary experiments, salicylate concentrations higher than 100 mM used to study diffusion in this work caused elevation of CAP thresholds throughout the entire frequency range. This flooding of the whole cochlea with salicylate was due apparently to overloading of the cochlear clearing and other possible mechanisms involved. In this case, the dynamic equilibrium between diffusion and clearing and steady-state salicylate concentrations cannot be reached and our model cannot be applied. From an experimental standpoint, the use of higher concentrations of salicylate also made time-dependent estimates of diffusion impossible for high frequencies because the clearing mechanism in the basal turn became almost immediately saturated following salicylate application. As a result of the clearing overload and other unidentified processes, salicylate is accumulated throughout the cochlea affecting all the frequencies as it was observed in our experiments where we applied 1M salicylate to the RW. Of course, therapeutic use of concentrated drug formulations in order to overcome the issues raised by this study could be problematic due to likely side effects and/or restricted aqueous solubility and thus is not a practical solution.

While the steady-state distribution of concentrations, which is the basis for conclusions in this study, is fitted well by the simple diffusion model, the responses for the lower frequencies became gradually slower compared to the model predictions (Figure S1). This may happen because salicylate action on the cochlear amplifier is not limited by its block of the OHC motility (e.g., Russell and Schauz, 1995; Wu et al., 2010) as it is assumed in the model. A compensatory effect from a hypothetical mechanism maintaining cochlear homeostasis and OHC sensitivity and responsible for the “bounce” phenomenon after exposure to loud sounds (Kirk et al., 1997; Drexl et al., 2014) may also explain delayed threshold elevation at subtle salicylate concentrations in the low-frequency cochlear region. Finally, salicylate concentration at the cochlear base may be diluted by the cerebrospinal fluid coming through the cochlear aqueduct into the perilymph which becomes hyperosmotic due to relatively high salicylate concentration at the base. None of these mechanisms should, however, affect our conclusion about the magnitude of the steady-state concentration gradients along the ST.

For a cochlea of given geometry, the concentration gradient along the ST depends only on the relationship between diffusion and clearing and is drug specific. In terms of the current study, it is the value of k_d/k_c, which defines the ratio between the amount of drug entering through a unit surface of the ST normal to the direction of diffusion and leaving it through a unit area of the side walls within the same time period. Salicylate, which is readily cleared from the ST (k_d/k_c = 1.6968), does not in practice diffuse into the cochlear apex and the resultant theoretical base-to-apex concentration gradient is extremely high (red cross in Figure 4A). Drugs which are better retained in the ST (i.e., have higher k_d/k_c ratio) form smaller concentration gradients, but this is traded for the considerably longer time it takes for these drugs to reach steady-state concentrations in the cochlear apex (Figure 4B). Hence, this approach may not be practical when there is only a short time window for the treatment of a specific cochlear disorder. Also, using a drug form which is better retained in the ST will lead to larger concentration differences between the ST and surrounding tissue. This may be a problem for drugs with narrow therapeutic windows unless an inactive form of the drug is used for even distribution along the ST through diffusion and it is activated only when the drug is cleared into the surrounding tissue.

Because the retention of a drug at the RW does not lead to a leveling of its concentration along the cochlear spiral (see also Plontke et al., 2007b), different strategies for drug delivery to the cochlear apex should be employed. Stable drug loaded nanocarriers (Zou et al., 2014; Li et al., 2017; Kamalov et al., 2018) which can stay in the ST long enough without being cleared into the surrounding tissue may be a feasible option. When the concentration of nanocarriers along the ST reaches a constant level, the encapsulated drug could be released from the carriers through thermal or light activation (Karimi et al., 2016, 2017; Yuan et al., 2017) to obtain sufficient drug concentrations along the entire cochlear spiral. A potential problem with this approach is the substantial increase in time required to reach the equilibrium base-to-apex gradient of nanocarrier concentrations, due to the substantially smaller diffusion coefficients of even the smallest liposomes and micelles, compared to lone drug molecules (Figure 4B) (del Amo et al., 2017).

Drug loaded nanoparticles, however, could be used to take advantage of anatomical and cellular features of the cochlea which enable drug uptake through routes and pathways other than the ST route (Glueckert et al., 2018). Disulfiram loaded nanoparticles, for example, were observed in the apical part of the spiral ganglion just 1 day after their application to the RW and elevation of auditory brainstem response thresholds, due to disulfiram induced apoptosis of the ganglion neurons, was detected for frequencies corresponding to the cochlear apex within 2 days after application (Buckiová et al., 2012). Nanoparticles can also be effectively driven and distributed along the entire cochlea. Assisted diffusion of magnetically driven, prednisolone-loaded magnetic nanoparticle along the cochlea resulted in a significant increase in the protective effect of the drug against cisplatin-induced ototoxicity compared to intratympanic injections of prednisolone (Ramaswamy et al., 2017).

This study investigates passive drug diffusion along the intact cochlea when the drug is applied to the RW and highlights intrinsic problems with this method of local drug administration into the inner ear. Retaining drugs at the RW for an arbitrarily long time does not decrease its base-to-apex concentration gradient, which, at steady state, depends solely upon the relationship between drug diffusion along and clearing from the ST. Usage of drug-loaded nanocarriers which utilize the anatomical and cellular properties of the cochlea, and which can be actively distributed along the entire length of the cochlea seems to be a more promising approach.

Ethics Statement

All procedures involving animals were performed in accordance with UK Home Office regulations with approval from the University of Brighton Animal Welfare and Ethical Review Body.

Author Contributions

AL, IR, and TA conceived and designed the study. GB, SF, and AL performed the experiments and analyzed experimental results. IS, JK, and AL developed the model. IS performed numerical simulations and fitting to the experimental data. All authors contributed to analysis and discussion of the results. IS and AL wrote the manuscript with contribution from all authors.

Conflict of Interest Statement

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.

The authors thank Vadim Biktashev for his helpful advice and comments on the diffusion part of the model.

Supplementary Material

The Supplementary Material for this article can be found online at: /articles/10.3389/fncel.2019.00161/full#supplementary-material

References Borkholder

D. A.

Zhu

Frisina

R. D.

(2014). Round window membrane intracochlear drug delivery enhanced by induced advection. J. Control. Release 174, 171–176. 10.1016/j.jconrel.2013.11.021

24291333

Borkholder

D. A.

Zhu

Hyatt

B. T.

Archilla

A. S.

Livingston

I. I. I. W. J.

III Frisina

R. D.

(2010). Murine intracochlear drug delivery: reducing concentration gradients within the cochlea. Hear. Res. 268, 2–11. 10.1016/j.heares.2010.04.014

20451593

Bowe

S. N.

Jacob

(2010). Round window perfusion dynamics: implications for intracochlear therapy. Curr. Opin. Otolaryngol. Head Neck Surg. 18, 377–385. 10.1097/MOO.0b013e32833d30f0

20808222

Buckiová

Ranjan

Newman

T. A.

Johnston

A. H.

Sood

Kinnunen

P. K.

. (2012). Minimally invasive drug delivery to the cochlea through application of nanoparticles to the round window membrane. Nanomedicine 7, pp., 1339–1354. 10.2217/nnm.12.5

22475648

Chen

Kujawa

S. G.

McKenna

M. J.

Fiering

J. O.

Mescher

M. J.

Borenstein

J. T.

. (2005). Inner ear drug delivery via a reciprocating perfusion system in the guinea pig. J. Control. Release 110, pp., 1–19. 10.1016/j.jconrel.2005.09.003

16274830

Creber

N. J.

Eastwood

H. T.

Hampson

A. J.

Tan

O'Leary

S. J.

(2018). A comparison of cochlear distribution and glucocorticoid receptor activation in local and systemic dexamethasone drug delivery regimes. Hear. Res. 368, 75–85. 10.1016/j.heares.2018.03.018

29622283

del Amo

E. M.

Rimpelä

A. K.

Heikkinen

Kari

O. K.

Ramsay

Lajunen

. (2017). Pharmacokinetic aspects of retinal drug delivery. Prog. Retin. Eye Res. 57, 134–185. 10.1016/j.preteyeres.2016.12.001

28028001

Devare

Gubbels

Raphael

(2018). Outlook and future of inner ear therapy. Hear. Res. 368, 127–135. 10.1016/j.heares.2018.05.009

29804723

Drexl

Überfuhr

Weddell

T. D.

Lukashkin

A. N.

Wiegrebe

Krause

. (2014). Multiple indices of the ‘bounce’ phenomenon obtained from the same human ears. J. Assoc. Res. Otolaryngol. 15, 57–72. 10.1007/s10162-013-0424-x

24253659

El Kechai

Agnely

Mamelle

Nguyen

Ferrary

Bochot

(2015). Recent advances in local drug delivery to the inner ear. Int. J. Pharm. 494, 83–101. 10.1016/j.ijpharm.2015.08.015

26260230

Glueckert

Chacko

L. J.

Rask-Andersen

Wei

Handschuh

Schrott-Fischer

(2018). Anatomical basis of drug delivery to the inner ear. Hear. Res. 368, 10–27. 10.1016/j.heares.2018.06.017

30442227

Greenwood

D. D.

(1990). A cochlear frequency-position function for several species−29 years later. J. Acoust. Soc. Am. 87, 2592–2605. 10.1121/1.399052

2373794

Grewal

A. S.

Nedzelski

J. M.

Chen

J. M.

Lin

V. Y.

(2013). Dexamethasone uptake in the murine organ of Corti with transtympanic versus systemic administration. J. Otolaryngol.-Head N. 42:19. 10.1186/1916-0216-42-19

23663237

Haghpanahi

Gladstone

M. B.

Zhu

Frisina

R. D.

Borkholder

D. A.

(2013). Noninvasive technique for monitoring drug transport through the murine cochlea using micro-computed tomography. Ann. Biomed. Eng. 41, 2130–2142. 10.1007/s10439-013-0816-4

23636576

Hallworth

(1997). Modulation of outer hair cell compliance and force by agents that affect hearing. Hear. Res. 114, 204–212. 10.1016/S0378-5955(97)00167-6

9447933

Hao

S. K.

(2019). Inner ear drug delivery: recent advances, challenges, and perspective. Eur. J. Pharm. Sci. 126, 82–92. 10.1016/j.ejps.2018.05.020

29792920

Hargunani

C. A.

Kempton

J. B.

DeGagne

J. M.

Trune

D. R.

(2006). Intratympanic injection of dexamethasone: time course of inner ear distribution and conversion to its active form. Otol. Neurotol. 27, 564–569. 10.1097/01.mao.0000194814.07674.4f

16691147

Imamura

S. I.

Adams

J. C.

(2003). Distribution of gentamicin in the guinea pig inner ear after local or systemic application. J. Assoc. Res. Otolaryngol. 4, 176–195. 10.1007/s10162-002-2036-8 Jones

G. P.

Lukashkina

V. A.

Russell

I. J.

Elliott

S. J.

Lukashkin

A. N.

(2013). Frequency-dependent properties of the tectorial membrane facilitate energy transmission and amplification in the cochlea. Biophys. J. 104, 1357–1366. 10.1016/j.bpj.2013.02.002

23528095

Kamalov

M. I.

Dǎng

Petrova

N. V.

Laikov

A. V.

Luong

Akhmadishina

R. A.

. (2018). Self-assembled nanoformulation of methylprednisolone succinate with carboxylated block copolymer for local glucocorticoid therapy. Colloids Surf. B Biointerfaces 164, pp. 78–88. 10.1016/j.colsurfb.2018.01.014

29413623

Karimi

Sahandi Zangabad

Baghaee-Ravari

Ghazadeh

Mirshekari

Hamblin

M. R.

(2017). Smart nanostructures for cargo delivery: uncaging and activating by light. J. Am. Chem. Soc. 139, 4584–4610. 10.1021/jacs.6b08313

28192672

Karimi

Sahandi Zangabad

Ghasemi

Amiri

Bahrami

Malekzad

. (2016). Temperature-responsive smart nanocarriers for delivery of therapeutic agents: applications and recent advances. ACS Appl. Mater. Interfaces 8, 21107–21133. 10.1021/acsami.6b00371

27349465

Kirk

D. L.

Moleirinho

Patuzzi

R. B.

(1997). Microphonic and DPOAE measurements suggest a micromechanical mechanism for the ‘bounce’ phenomenon following low-frequency tones. Hear. Res. 112, 69–86. 10.1016/S0378-5955(97)00104-4

9367230

Chao

Brant

O'Malley

Jr Tsourkas

(2017). Advances in nano-based inner ear delivery systems for the treatment of sensorineural hearing loss. Adv. Drug Deliv. Rev. 108, 2–12. 10.1016/j.addr.2016.01.004

26796230

Lide

D. R.

(2002). CRC Handbook of Chemistry and Physics, 83rd Edn. Boca Raton, FL: CRC Press. Lukashkin

A. N.

Lukashkina

V. A.

Russell

I. J.

(2002). One source for distortion product otoacoustic emissions generated by low- and high-level primaries. J. Acoust. Soc. Am. 111, 2740–2748. 10.1121/1.1479151

12083209

Lukashkin

A. N.

Russell

I. J.

(1999). Analysis of the f2– f1 and 2f1– f2 distortion components generated by the hair cell mechanoelectrical transducer: Dependence on the amplitudes of the primaries and feedback gain. J. Acoust. Soc. Am. 106, 2661–2668. 10.1121/1.428096 Meaud

Grosh

(2014). Effect of the attachment of the tectorial membrane on cochlear micromechanics and two-tone suppression. Biophys. J. 106, 1398–1405. 10.1016/j.bpj.2014.01.034

24655515

Mountain

D. C.

Hubbard

A. E.

McMullen

T. A.

(1983). Electromechanical processes in the cochlea, in Mechanics of Hearing, eds. de Boer

Viergever

M. A.

(Dordrecht: Springer), 119–126. Mynatt

Hale

S. A.

Gill

R. M.

Plontke

S. K.

Salt

A. N.

(2006). Demonstration of a longitudinal concentration gradient along scala tympani by sequential sampling of perilymph from the cochlear apex. J. Assoc. Res. Otolaryngol. 7, 182–193. 10.1007/s10162-006-0034-y

16718612

Elliott

S. J.

Baumgart

(2016). Finite-element model of the active organ of Corti. J. R. Soc. Interface 13:20150913. 10.1098/rsif.2015.0913

26888950

Nuttall

A. L.

Ricci

A. J.

Burwood

Harte

J. M.

Stenfelt

Cayé-Thomasen

. (2018). A mechanoelectrical mechanism for detection of sound envelopes in the hearing organ. Nat. Commun. 9:4175. 10.1038/s41467-018-06725-w

30302006

Ohyama

Salt

A. N.

Thalmann

(1988). Volume flow rate of perilymph in the guinea-pig cochlea. Hear. Res. 35, 119–129. 10.1016/0378-5955(88)90111-6

3198505

Plontke

S. K.

Biegner

Kammerer

Delabar

Salt

A. N.

(2008). Dexamethasone concentration gradients along scala tympani after application to the round window membrane. Otol. Neurotol. 29, 401–406. 10.1097/MAO.0b013e318161aaae

18277312

Plontke

S. K.

Mynatt

Gill

R. M.

Borgmann

Salt

A. N.

(2007a). Concentration gradient along the scala tympani after local application of gentamicin to the round window membrane. Laryngoscope, 117, 1191–1198. 10.1097/MLG.0b013e318058a06b

17603318

Plontke

S. K.

Siedow

Wegener

Zenner

H. P.

Salt

A. N.

(2007b). Cochlear pharmacokinetics with local inner ear drug delivery using a three-dimensional finite-element computer model. Audiol. Neurotol. 12, 37–48. 10.1159/000097246

17119332

Ramaswamy

Roy

Apolo

A. B.

Shapiro

Depireux

D. A.

(2017). Magnetic nanoparticle mediated steroid delivery mitigates cisplatin induced hearing loss. Front. Cell. Neurosci. 11:268. 10.3389/fncel.2017.00268

28955202

Rivera

Sanz

Camarero

Varela-Nieto

(2012). Drug delivery to the inner ear: strategies and their therapeutic implications for sensorineural hearing loss. Curr. Drug Deliv. 9, 231–242. 10.2174/156720112800389098

22283653

Ruggero

M. A.

Narayan

S. S.

Temchin

A. N.

Recio

(2000). Mechanical bases of frequency tuning and neural excitation at the base of the cochlea: comparison of basilar-membrane vibrations and auditory-nerve-fiber responses in chinchilla. Proc. Natl. Acad. Sci. U.S.A. 97, 11744–11750. 10.1073/pnas.97.22.11744

11050204

Russell

I. J.

Schauz

(1995). Salicylate ototoxicity: effects on stiffness and electromotility of outer hair cells isolated from the guinea pig cochlea. Auditory Neurosci. 1, 309–319. Saijo

Kimura

R. S.

(1984). Distribution of HRP in the inner ear after injection into the middle ear cavity. Acta Otolaryngol. 97, 593–610. 10.3109/00016488409132937

6464711

Salt

A. N.

Hirose

(2018). Communication pathways to and from the inner ear and their contributions to drug delivery. Hear. Res. 362, 25–37. 10.1016/j.heares.2017.12.010

29277248

Salt

A. N.

Kellner

Hale

(2003). Contamination of perilymph sampled from the basal cochlear turn with cerebrospinal fluid. Hear. Res. 182, 24–33. 10.1016/S0378-5955(03)00137-0

12948598

Salt

A. N.

(2001). Quantification of solute entry into cochlear perilymph through the round window membrane. Hear. Res. 154, 88–97. 10.1016/S0378-5955(01)00223-4

11423219

Salt

A. N.

Plontke

S. K.

(2009). Principles of local drug delivery to the inner ear. Audiol. Neurotol. 14, 350–360. 10.1159/000241892

19923805

Salt

A. N.

Plontke

S. K.

(2018). Pharmacokinetic principles in the inner ear: influence of drug properties on intratympanic applications. Hear. Res. 368, 28–40. 10.1016/j.heares.2018.03.002

29551306

Schuknecht

H. F.

(1956). Ablation therapy for the relief of Meniere's disease. Laryngoscope 66, 859–859. 10.1288/00005537-195607000-00005

13358249

Takagi

Taki

Sakane

Nadai

Sezaki

Oku

. (1998). A new interpretation of salicylic acid transport across the lipid bilayer: implications of pH-dependent but not carrier-mediated absorption from the gastrointestinal tract. J. Pharmacol. Exp. Ther. 285, 1175–1180. Temchin

A. N.

Rich

N. C.

Ruggero

M. A.

(2008). Threshold tuning curves of chinchilla auditory-nerve fibers. I. Dependence on characteristic frequency and relation to the magnitudes of cochlear vibrations. J. Neurophysiol. 100, 2889–2898. 10.1152/jn.90637.2008

18701751

Thorne

Salt

A. N.

DeMott

J. E.

Henson

M. M.

Henson

O. W.

Jr. Gewalt

S. L.

(1999). Cochlear fluid space dimensions for six species derived from reconstructions of three-dimensional magnetic resonance images. Laryngoscope 109, 1661–1668. 10.1097/00005537-199910000-00021

10522939

Kim

H. J.

Yamoah

E. N.

Nuttall

A. L.

(2010). Effect of salicylate on KCNQ4 of the guinea pig outer hair cell. J. Neurophysiol. 103, 1969–1977. 10.1152/jn.01057.2009

20147414

Yates

G. K.

(1990). Basilar membrane nonlinearity and its influence on auditory nerve rate-intensity functions. Hear. Res. 50, 145–162. 10.1016/0378-5955(90)90041-M

2076968

Yuan

Huan

Liu

Zhang

. (2017). NIR light-activated drug release for synergetic chemo-photothermal therapy. Mol. Pharm. 14, 242–251. 10.1021/acs.molpharmaceut.6b00820

27983855

Zou

Sood

Zhang

Kinnunen

P. K.

Pyykkö

(2014). Pathway and morphological transformation of liposome nanocarriers after release from a novel sustained inner-ear delivery system. Nanomedicine 9, 2143–2155. 10.2217/nnm.13.181

24471501

Appendix

Only diffusion along the long axis x of the tube of decreasing diameter is considered (Figure 1A). The concentration c within each cross-section for a fixed instance t is assumed to be constant, i.e., c = c(x, t) is independent to the y axis. If the area of the cross-section is S(x), then the flux J along the x axis is given by:

(A1)J(x,t)=-S(x)·kd·dc(x,t)dx,

where k_d is the diffusion coefficient.

At the same time the clearing of salicylate from the tube of length Δx, its perimeter P(x) and with an area of surface P(x)Δx, can be described as:

(A2)Cl(x,t)=c(x,t)·kc·P(x)·Δx,

where k_c is the clearing coefficient.

The balance of fluxes and clearing in the volume between x₀ and x₁ can be described as:

(A3)S(x)·Δx·dc(x,t)dt=S(x1)·kd·dc(x1,t)dx-S(x0)·kd·dc(x0,t)dx-c(x,t)·kc·P(x)·Δx,

where Δx = x₁ − x₀ is positive and x is in [x₀, x₁].

Divide both sides of Equation (A3) by Δx and rearrange it as follows:

(A4)S(x)·dc(x,t)dt=S(x1)·kd·dc(x1,t)/dx-S(x0)·kd·dc(x0,t)/dxΔx-c(x,t)·kc·P(x).

Take the limit of Δx converging to zero and divide by S(x), the following is obtained:

(A5)dc(x,t)dt=1S(x)·ddx(S(x)·kd·dc(x,t)dx)-c(x,t)·L(x),

where L(x)=kc·P(x)S(x) is an integral coefficient of clearing.

The perimeter is

(A6)P(x)=2π·(mx+a)

and the area is

(A7)S(x)=π·(mx+a)2,

where m = (b − a)/l.

The integral coefficient can be written as:

(A8)L(x)=kc·P(x)S(x)=2kcr(x).

Thus, the diffusion can be described by the following partial differential equation

(A9)dc(x,t)dt=1S(x)·ddx(S(x)·kd·dc(x,t)dx)-c(x,t)·2kcr(x),

with the boundary conditions

(A10)c(0,t)=crw,

(A11)kddc(l,t)dx=0

and initial conditions

(A12)c(0,0)=crw;x=0,

(A13)c(x,0)=0;x>0.

Equation (A9) can be rewritten in the following non-dimensional form

(A14)du(χ,τ)dτ=1((b-a)χ+a)2·ddχ(((b-a)χ+a)2·du(χ,τ)dχ)-u(χ,τ)·2l2/ratio(b-a)χ+a,

with the boundary conditions

(A15)u(0,τ)=1,

(A16)ratio·du(1,τ)dχ=0

and the initial conditions

(A17)c(0,0)=1;χ=0,

(A18)c(χ,0)=0;χ>0,

where u = c/c_rw, τ=t·kd/l2 and χ = x/l.

Funding. The research was funded by a grant from the Medical Research Council (MR/N004299/1). TA is partially supported by the Program of Competitive Growth at Kazan Federal University.

�㾩julia�������߲���

�㾩julia��߲��