Next Article in Journal
Impact of Machine Perfusion on Biliary Complications after Liver Transplantation
Next Article in Special Issue
Auxin Controlled by Ethylene Steers Root Development
Previous Article in Journal
Regulation of Leptin Methylation Not via Apoptosis by Melatonin in the Rescue of Chronic Programming Liver Steatosis
Previous Article in Special Issue
Advances in Understanding the Mechanism of Action of the Auxin Permease AUX1
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relative Contribution of PIN-Containing Secretory Vesicles and Plasma Membrane PINs to the Directed Auxin Transport: Theoretical Estimation

1
Mathematical Institute, Faculty of Science, Leiden University, 2333 CA Leiden, The Netherlands
2
Institute of Science and Technology (IST) Austria, Am Campus 1, 3400 Klosterneuburg, Austria
3
Department Experimental Plant Biology, Faculty of Science, Charles University, 12843 Prague, Czech Republic
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2018, 19(11), 3566; https://doi.org/10.3390/ijms19113566
Submission received: 17 September 2018 / Revised: 2 November 2018 / Accepted: 4 November 2018 / Published: 12 November 2018
(This article belongs to the Special Issue Auxins and Cytokinins in Plant Development)

Abstract

:
The intercellular transport of auxin is driven by PIN-formed (PIN) auxin efflux carriers. PINs are localized at the plasma membrane (PM) and on constitutively recycling endomembrane vesicles. Therefore, PINs can mediate auxin transport either by direct translocation across the PM or by pumping auxin into secretory vesicles (SVs), leading to its secretory release upon fusion with the PM. Which of these two mechanisms dominates is a matter of debate. Here, we addressed the issue with a mathematical modeling approach. We demonstrate that the efficiency of secretory transport depends on SV size, half-life of PINs on the PM, pH, exocytosis frequency and PIN density. 3D structured illumination microscopy (SIM) was used to determine PIN density on the PM. Combining this data with published values of the other parameters, we show that the transport activity of PINs in SVs would have to be at least 1000× greater than on the PM in order to produce a comparable macroscopic auxin transport. If both transport mechanisms operated simultaneously and PINs were equally active on SVs and PM, the contribution of secretion to the total auxin flux would be negligible. In conclusion, while secretory vesicle-mediated transport of auxin is an intriguing and theoretically possible model, it is unlikely to be a major mechanism of auxin transport in planta.

1. Introduction

The plant hormone auxin (indole-3-acetic acid; IAA) is subject to intercellular polar auxin transport (PAT) mediated by diffusion and the action of efflux and influx carriers. PAT generates local auxin maxima that are crucial for a plethora of developmental processes [1]. Therefore, studies into the mechanism of auxin transport and its regulation have had a prominent place in the auxin field. Auxin maxima-driven developmental events depend on the activity of the PIN-formed (PIN) IAA efflux carriers, which provide directionality to intercellular IAA transport through their asymmetric subcellular localization [1,2,3]. The PIN proteins constitutively cycle between the plasma membrane (PM) and endosomal compartments [4,5,6], and the developmental importance of this energetically demanding phenomenon has not been unequivocally explained. Three hypotheses have been suggested [7] to explain the requirement of PIN recycling: (1) recycling enables rapid relocation of PINs and thereby the rapid redirection of auxin transport; (2) PINs serve as IAA transceptors (transporters and receptors at the same time) and their recycling is important for the process of signal transduction; and (3) the PIN-containing secretory vesicles (SVs) are filled with auxin, which is released into the apoplast upon fusion of the SV with the PM. This process is analogous to neurotransmitter release in animals, and is important for the transport of auxin as suggested by previous reports [7,8] (Figure 1).
Several studies claimed to have proven the vesicular transport of IAA [9,10,11]. However, the validity of the evidence presented in most of these studies has been questioned by many researchers in the field, and thus the hypothesis of vesicular transport of auxin remains controversial [12,13]. A major reason this controversy exists is the inability to uncouple vesicular trafficking from PIN occurrence at the PM using existing cell biology methods. For example, if one could genetically or pharmacologically completely and specifically block the movement of PIN-containing secretory vesicles, it would be impossible to conclude whether the resulting reduction in auxin transport was caused by the lack of IAA secreting vesicles or by the lack of PINs at the PM.
Here, we have constructed a simple mathematical model to estimate the parameters under which auxin transport by intracellular vesicles could explain the measured net flux values of PAT. We show that a vesicular transport model would require at least 1000 times greater PIN activity than the conventional PM transport model to generate the same net flux values.

2. Results

2.1. Model Assumptions: PINs Can Be Active on the Endomembrane Vesicles Exclusively or in Addition to Their Activity on the Plasma Membrane (PM)

PAT is characterized by the sustained velocity of auxin over long distances (several millimeters, which is much greater than a typical cell length). The magnitude of this velocity is in the range of 1 10 μ m s for different species and different types of plant tissue [14]. Arabidopsis root epidermal cells serve as a useful example because many of the physiological parameters that influence the PAT, and thus used in our calculations, have been experimentally measured in this cell type. Furthermore, vesicular auxin transport has been proposed to play a role in the transition zone of the root [9]. Thus, we chose to apply our model to this tissue. In the roots, auxin is transported towards the root tip inside the central cells (stele), whereas in the outer cell layer (lateral root cap cells and epidermal cells) it is transported in the reverse direction: from the root tip towards the shoot [15]. The velocity of this directed transport in Arabidopsis roots was found to be in a range 2 3 μ m s [16]. We therefore assume that v P A T 1 μ m s in the epidermal cells in the root transition zone and estimated the permeability of PINs needed to yield this value.
In our model we consider only secretory/recycling vesicles (SVs) that fuse with the PM, and consider these secretory vesicles as the method of PIN delivery from the endosome to the PM. For such auxin-transporting vesicles we considered two hypothetical scenarios: (1) PINs transport auxin exclusively in these endomembrane vesicles, which we compare to the case when PINs are active only on the PM (published results) [17,18,19]; and (2) PINs transport auxin both in the vesicles and on the PM. The auxin permeability of a PIN-containing membrane (PPIN) remains undetermined. However, PPIN has been calculated previously for a scenario when PINs are active solely on the PM ( P P I N P M ) [17,18]. Therefore, the permeability of PINs in vesicles P P I N v is a readout rather than input to our calculations. We calculated (1) the relative values of permeability P P I N v and the individual activity of PINs proteins needed for the vesicular transport model to generate the same net flux as the PM model and (2) the relative contribution of the vesicular transport to the net flux, assuming that PINs activity is the same in the vesicles and in the PM. A plan of the model and performed calculations is presented in Figure 2A.

2.2. Short Summary of Model Results

(1) For the case when only vesicular PINs are actively transporting auxin, we derived equations for auxin concentration inside the vesicles, which depends on the permeability of the vesicle membrane to auxin due to PINs ( P P I N v   [ μ m s ] , the measure of how many auxin molecules per second are transported through 1   μ m 2 area of a membrane when it is unit auxin concentration ( 1 mol μ m 2 ) outside the vesicle). At the same time, intravesicular auxin concentration determines the amount of auxin released into the cell wall. Therefore, we related macroscopic auxin transport velocity ( v P A T   [ μ m s ] ) to P P I N v and calculated the coefficient of proportionality between these variables for the physiological values of the other parameters: P P I N v - o n l y 42 · v P A T   .
However, if PINs are only active on the PM, then PINs permeability P P I N PM - o n l y v P A T   [17]. Thus, for the minimal velocity v P A T = 1 μ m s , as observed in epidermal root cells [14], the lowest possible value of permeability is P P I N v - o n l y   42 μ m s for vesicular transport, which is much higher than the permeability required for PM transport P P I N P M - o n l y   1 μ m s .
As the permeability is proportional to the density of PINs (the number of PIN molecules per unit area of the membrane), we generalized our equations via dividing permeability values by the density of PINs. The resulting quantity is the activity of individual PINs. We show that on the vesicle membrane PINs’ activity has to be 3 orders of magnitude higher than in the scenario where they are active only on the PM, in order to produce the same PAT velocity:
PINs   activity   on   SV   only PINs   activity   on   PM   only >   4200 .
(2) For the case when PINs are active both in vesicles and on the PM, our model provides the ratio of auxin amount passing through these two mechanisms depending on individual activities of PIN molecules in these domains:
Auxin   flux   through   PM Auxin   flux   through   SV = 4200   PINs   activity   on   PM PINs   activity   on   SV .
Thus, if activities are equal, auxin flux through SVs contributes no more than 0.02% of the total flux.
Conclusion: The transport of auxin by vesicles is ~1000 times less efficient than through PINs active on the PM. Permeability values of PINs (   P P I N v and P P I N P M ) still await direct measurements, but values as high as P P I N v = 42 μ m s , as estimated in our study, are unlikely. Therefore, this argues against vesicular transport of auxin in SVs as a major mechanism of directional auxin transport.

2.3. Detailed Model Description

2.3.1. PINs Active Only in Vesicles. Comparison to the Case When PINs are Active Only on PM

List of parameters used in the model is provided in the Table 1. Assumptions taken throughout the model description are listed in Table 2.

2.3.1.1. How Much Auxin Should a Vesicle with Active PIN Contain to Produce the PAT?

(1) Size of a Vesicle
The size of exocytotic vesicles that are considered to perform auxin transport is well known; they are spherical in shape and their diameter was reported to be d = 0.06   0.08   μ m   [20] (see Table 1 for corresponding surface area and volume).
(2) Biochemical Constituents and Transport Processes, Equation for Auxin Concentration inside the Vesicle
The total concentration of auxin (IAA) is the mixture of anions IAA and protonated form (IAAH), and the ratio between them depends on pH of the solution:
p H c = 7 in cytoplasm, ~99% of IAA, 1% of IAAH
p H v = 5.5   t o   6.5 in vesicle, ~83–98% of IAA, 17–2% of IAAH
p K a = 4.8 for IAA.
The fractions of IAA in anion form in cytosol and vesicle are denoted by f a c and f a v respectively. The values were computed using:
f a v = 1 1 + 10 p K a p H v ,
f a c = 1 1 + 10 p K a p H c .
Thus, ( 1 f a v ) and ( 1 f a c ) are fractions of IAAH in vesicle and cytosol, respectively (see Figure 2B and Table 3 for IAA fractions depending on pH).
We denote:
A c —cytoplasmic concentration of auxin [mol/L],
A v —concentration of auxin inside the vesicle [mol/L].
Fast diffusion in the cytoplasm ( D 670   μ m 2 / s [17,25]) ensures fast mixing of IAA inside the small vesicle, that takes τ d i f f 10 6 s in vesicles of diameter d 100   nm . Thus, IAA concentration inside the vesicle can be considered homogenous and increases proportionally to the flux through the membrane:
d A v d t = 1 V v J i n n e t ,
where J i n n e t is the net influx of IAA into vesicle (units [mol/s]), which has two components: transport of anions IAA by PINs and transport of IAA and IAAH by diffusion through the membrane:
  J i n n e t = J P I N I A A + J d i f f .
PINs transport anions IAA and are presumed to be permeable only in one direction (into vesicle) resulting in the unknown permeability   P P I N v of the vesicular membrane. We neglect dependence of PIN-mediated influx on intravesicular auxin concentration and assume that reverse permeability is zero. By doing so, we make PINs more efficient in our considerations than they can possibly be in reality. This approach is acceptable because our aim is to find the lower bound of the permeability of PINs which is able to produce the physiologically observed PAT transport.
Thus, PIN-mediated transport is simplified to:
  J P I N I A A = P P I N v   S v   f a c   A c .
Diffusive flux across the membrane is governed by Fick’s law:
  J d i f f = J d i f f I A A H + J d i f f I A A = P d i f f I A A H S v   ( ( 1   f a c   ) A c ( 1   f a v   ) A v ) + P d i f f I A A S v   ( f a c A c f a v   A v ) .
Diffusional permeability for IAAH is known: P d i f f I A A H =   0.5   μ m / s   [18,24]. IAAH is assumed to have the same diffusional permeability P d i f f I A A H in both directions.
Anion diffusion permeability is estimated to be much lower: P d i f f I A A P d i f f ( IAAH ) 100 [24]. However, at p H v = 7 , concentration of [ I A A ] 165 · [ I A A H ] and thus diffusional flux of I A A . can be of the same order of magnitude as IAAH flux. Nevertheless, to simplify the equations we omit the term for   I A A diffusion. As diffusion of I A A acts against PINs-mediated transport, this simplification also favors auxin accumulation into vesicles.
Combining Equations (7) and (8) and neglecting J d i f f I A A   result in the following expression for the net total flux:
J i n n e t =   P P I N v   S v   f a c   A c + P d i f f I A A H S v   ( ( 1   f a c   ) A c ( 1   f a v   ) A v ) .  
Thus, Equation (5) can be rewritten as:
d A v d t = A c · α λ · A v ,
where,
α = P P I N v   S v V v   f a c + P d i f f I A A H S v V v ( 1   f a c   )   ,
λ = P d i f f I A A H S v V v ( 1   f a v ) .
(3) Maximum Loading of Vesicles with Auxin is Proportional to the Permeability of PINs
We assume that vesicles exist in the cytoplasm long enough that the internal vesicular auxin concentration reaches its maximum: steady state concentration A v . The analytical solution of Equation (10) reads:
A v ( t ) = A c α λ ( 1 e λ t ) + e λ t A v ( 0 ) .
Thus, characteristic time to load vesicle from zero concentration till A v = A c α λ ( 1 1 e ) is τ l o a d = 1 λ = 4.6   s for p H v = 7 (and faster for p H v < 7 ). For times t τ l o a d vesicles will be fully filled. We consider fully filled vesicles because our aim is to find the minimal requirements for the vesicular transport, and partly filled vesicles would require higher permeability values to produce the PAT.
At steady state, when d A v d t = 0 , Equation (10) simplifies to:
A c · α = λ · A v .
Thus, the number of IAA molecules (in moles) in one vesicle:
N v m a x = A v m a x V v = A c V v α λ .
Substituting α and λ by their expressions gives:
N v m a x = A c V v [ P P I N v P d i f f I A A H f a c 1 f a v + 1 f a c 1 f a v ] ,
that shows that the number of auxin molecules loaded inside the vesicle is proportional to P P I N v .
We can also rewrite this equation in the form:
N v m a x = A c V v R ,
where,
  R = A v m a x A c = P P I N v P d i f f I A A H f a c 1 f a v + 1 f a c 1 f a v  
is the accumulation ratio—the ratio of intravesicular IAA concentration to that in the cytoplasm surrounding the vesicle. Expressions for accumulation ratio and maximum loading of vesicles for different p H v are presented in Table 3.
(4) Lower Bound for Accumulation Ratio Necessary to Produce the PAT
In this section we estimated the auxin accumulation ratio R in vesicles, which is necessary to produce v P A T if auxin is transported via vesicles only. Auxin flux density (molecules passed per unit area per second) on the polar domain of the PM that corresponds to v P A T velocity is given by:
  Φ c e l l a p o p l a s t = v P A T · A c     [ m o l   μ m 2 s ] .
Maximum flux of auxin, that vesicles can carry through the PM is:
Φ v = N v m a x · F m a x +     [ m o l   μ m 2 s ] ,
where F m a x + is the maximum exocytosis frequency (vesicles per second fusing with the unit area of cell face). This flux should be not lower than the yielded flux:
Φ v Φ c e l l a p o p l a s t .
This equation takes into account that diffusion diminishes the directional transport (see Box 4 for explanation).
Thus,
  N v m a x · F m a x + v P A T · A c .
which gives, using Equation (17):
  R · A c V v · F m a x + v P A T · A c .
Consequently,
R v P A T V v · F m a x + .
Substituting F m a x + with its expression (see Equation (1.7) in Box 1) F m a x + = l n 2 S v   ·   τ 1 / 2 ( ρ P M ρ v ) m a x , we can find the lower bound of the accumulation ratio, that is required for v P A T   :
  R v P A T   · S v     ·     τ 1 / 2 V v · l n 2 ( ρ v ρ P M ) m i n .
Note, that intracellular auxin concentration A c cancels (Equation (23)), which means that the equations are valid for any A c . Nevertheless, we checked that vesicles will contain at least one molecule of auxin for the physiological value of A c , a condition required for auxin transport to be theoretically possible (see Box 3). We also calculated the minimum accumulation ratio required for v P A T   = 1   μ m / s (Box 3).
Box 1. Expression for F + , frequency of secretory vesicles fusing with PM, and estimation of its maximum value F m a x + .
In the following text we will consider only the polar domain of PM, which is the front membrane in the direction of PAT. It contains more PINs than the neighboring sides of the PM, and is visible by fluorescence microscopy [23]. The density of PINs depends on the rate of their delivery to PM by vesicles and on the rate of their removal. We assume that within the polar domain PINs are homogeneously distributed. Mass conservation for the number of PINs on a polar domain of the PM of area S P M then reads:
d n d t = n v F + S P M β · n ,
where n   —number of PINs on PM, n v —number of PINs in one vesicle, F + is frequency of exocytosis. PINs come to PM via F + S P M vesicles per second. Removal of PINs is proportional to n with the decay coefficient β , which describes any possible decay mechanism (through endocytosis, diffusion to other PM domains, degradation etc.).
The decay coefficient β can be found from the following arguments. The equation for PINs that are already delivered to the PM (let us call them “old PINs”) and are only removed from the PM is the same as Equation (1.1) without exocytosis term:
d n o l d d t = β · n o l d .
The solution of this equation reads:
n o l d ( t ) = n o l d ( 0 ) e β · t .
Introducing half-life τ 1 / 2 , time when 50% of PINs has been removed ( n o l d ( τ 1 / 2 ) = n o l d ( 0 ) / 2 ), gives β :
β = l n 2 τ 1 / 2 .
τ 1 / 2 = 1.3 × 10 4   s , which has been measured by Jásik et al. ([23]) for functional PIN2-Dendra translational fusion protein in epidermal cells of Arabidopsis root. In their experiments, Jásik et al. ([23]) visualized “old PINs” separately from newly arriving PINs by use of photoconvertion of PIN2-Dendra.
Assuming that β is constant and does not depend on the method of PINs visualization, and assuming a steady state of PINs concentration on the PM, d n d t = 0 , Equation (1.1) transforms to:
n v F + S P M = β · n .
Equations (1.4) and (1.5) give us expression for F + , the frequency of vesicle fusion with the PM:
F + = l n 2 S P M     ·     τ 1 / 2 n n v ,
We rewrite it using ρ P M = n S P M , ρ v = n v S v —densities of PINs on the PM and vesicle membrane respectively:
F + = l n 2 S P M · τ 1 / 2 n n v = l n 2 S P M   ·   τ 1 / 2 ρ P M ρ v S P M S v = l n 2 S v · τ 1 / 2 ρ P M ρ v .
Finally, introducing values for S v and τ 1 / 2   gives:
F + = 0.69 1.1   × 10 2   μ m 2   · 1.3 × 10 4   s   ρ P M ρ v = ρ P M ρ v 0.005 1   μ m 2   s   ,  
Introducing an estimation of minimum value for densities ratio ( ρ v ρ P M ) m i n = 0.01 (see Box 2) into Equation (1.8) gives the upper bound for F + :
F m a x + = 0.5 v e s i c l e s   μ m 2   s .
The experimentally measured endocytosis rate gives an alternative estimate of the exocytosis rate.
In accordance with the above result, the measured rate of endocytosis is ~ 0.5 v e s i c l e s   μ m 2   s (Table 4). The total area of vesicles fusing with the PM in one cell has to be balanced by the area of endocytosed vesicles. Given that sizes of exocytotic and endocytotic vesicles are the same, the maximal rate of endocytosis gives a rough estimate of the maximal possible rate of exocytosis, which is in the same order of magnitude as our maximal theoretical estimate.
It is worth noting that the half-life of PINs τ 1 / 2 was measured for “old” PINs and does not account for the possibility of “old” PINs being removed and brought back to the membrane via vesicles. This process would effectively reduce τ 1 / 2 and allow for a higher F m a x + . However, based on our comparison with experimental values of endocytosis, we argue that F m a x + cannot be much higher than 0.5 v e s i c l e s   μ m 2   s . This restriction comes from the physiological in vivo rate of endocytic vesicle formation (taking ~18–22 s per vesicle, see Table 4 [26]) and limited area of the cell membrane. Thus, as ( ρ v ρ P M ) m a x = 1 , τ 1 / 2 can be reduced not more than by 100.
Box 2. Estimation of minimal PINs density ratio on vesicles and the PM: ( ρ v ρ P M ) m i n .
The lowest possible density of PINs on the vesicular membrane is one PIN molecule per vesicle:   ( ρ v ) m i n =   1 S v = 50 1 μ m 2 .
The maximum PIN density can be calculated as follows: diameter of the globular 60 kDa protein is ~8 nm. If PINs on the PM are densely packed, their density is at most ( ρ P M   ) m a x   2 × 10 4 1 μ m 2 , which is the maximum value of PIN density on the PM.
Thus,
( ρ v ρ P M ) m i n = ( ρ v ) m i n ( ρ P M ) m a x = 0.002 .
However, PINs probably cannot reach this maximum density in vivo, because numerous other proteins occupy space in the membrane. Also, quantification of PIN2-GFP in live epidermal cells of Arabidopsis root cells, observed using 3D structured illumination microscopy (SIM) with an x and y resolution of 106 nm and 104 nm respectively (Figure 3a), and an average acquisition speed faster than the average rate of endocytosis (see Methods & Table 4), showed that separated source-spots of GFP signal are always resolved with maximum density of 5.6 spots μ m 2 (Figure 3b). Each of the spot contains at least one PIN-GFP protein. However, most spots are large: mean area of the spots is 0.034 ± 0.002   μ m 2 , max area 0.098   μ m 2 , min area 0.01   μ m 2 (corresponds to resolution limit). Thus, most spots are likely to contain many PINs, which together fill at most 1/5 part of the membrane area (as 5.6 1 μ m 2 · 0.034   μ m 2 = 0.19 ). This result argues that upper bound for the PIN density on PM is less than close-packing: ( ρ P M   ) m a x 4 × 10 3 1 μ m 2 . Correspondingly,
( ρ v ρ P M ) m i n = 0.01
is a more realistic value of minimal PINs density ratio for epidermal root cells.
Box 3. How many IAA molecules should be inside of one vesicle?
Rewriting Equation (25) as inequality gives the accumulation ratio R, which is necessary for minimal PAT ( v P A T = 1 μ m s ) :
R v m i n   S v     ·     τ 1 2 V v · l n 2 ( ρ v ρ P M ) m i n = 1 μ m s   2 ×   10 2 μ m 2   1.3   ×   10 4   s 2.7   ×   10 4 μ m 3   · 0.69 0.002 > 2800 .
This minimal required accumulation ratio is valid for any A c . Nevertheless, we have to check, that vesicles will contain at least one molecule. In the volume of a vesicle (maximum 2.7 × 105 nm3) there will be initially N v t = 0 = A c V v = 0.005 molecules if cytoplasmic auxin concentration A c = 30   nM , which is a plausible estimate for the root epidermis [29]. For the lowest possible accumulation ratio R = 2800 , the number of molecules in one vesicle will be:
N v = A c V v · R 14 ,
which is a minimal requirement for vesicular transport.

2.3.1.2. Minimal Permeability of PINs in Vesicles Necessary for Vesicular Auxin Transport

By combining Equations (18) and (25) we can find the relation between vesicular PINs’ permeability P P I N v - o n l y and the transport velocity v P A T   in the case that PINs are only active in vesicles and thus all PAT is produced by vesicles:
  P P I N v - o n l y P d i f f I A A H f a c 1 f a v + 1 f a c 1 f a v v P A T   S v · τ 1 / 2 V v · l n 2 ( ρ v ρ P M ) m i n .
Consequently,
  P P I N v - o n l y P d i f f I A A H ( v P A T   · S v · τ 1 / 2 V v · l n 2 ( ρ v ρ P M ) m i n 1 f a c 1 f a v ) / ( f a c 1 f a v ) .
Thus, because 1 f a c 1 f a v 1 can be neglected being much smaller than the first term v P A T   · S v   ·   τ 1 / 2 V v · l n 2 ( ρ v ρ P M ) m i n 10 4 ,
P P I N v - o n l y P d i f f I A A H · v P A T   · S v · τ 1 / 2 V v · l n 2 ( 1 f a v ) f a c ( ρ v ρ P M ) m i n .
Grouping all parameters except v P A T   and ( ρ v ρ P M ) m i n into constant Κ gives:
P P I N v - o n l y Κ · v P A T   · ( ρ v ρ P M ) m i n .  
From Table 3 we find m i n ( 1 f a v f a c ) = 0.006 for p H v = 7.0 , and using parameter values from Table 1, d = 0.08 μm, we calculate value:
Κ = ( P d i f f I A A H · S v · τ 1 2 V v · l n 2 · 1 f a v f a c ) m i n = 10 2 · 6 · 1.3 · 10 4 8 · 0.69 · 6 · 10 3 2 4240   .
Using ρ v ρ P M = 0.01 from Box 2, gives an estimate of the lower bound of PIN permeability in vesicles:
P P I N v - o n l y 42 · v P A T   .
For p H v < 7.0 bounds for P P I N v - o n l y P d i f f I A A H are provided in Table 3. For a pH value that we consider realistic, p H v = 6.2 permeability has to be higher: P P I N v - o n l y 260 · v P A T   .

2.3.1.3. Comparison to the Case when PINs are Active Only on the PM: Permeability due to PINs has to be Much Greater on the Vesicles than on PM to Produce the Same Auxin Transport Velocity

It would be informative to compare our estimate of the permeability value (Equation (30)) P P I N v - o n l y 42 · v P A T   with experimental measurements. Unfortunately, the permeability of PINs has never been measured directly in the intact tissues. Also measurements of the influx carrier permeability are not possible in explants like protoplasts, as PINs do not stay on the PM, but are instead internalized [18], resulting in no contribution to the efflux of auxin as in in vivo systems.
Fortunately, analogous estimates for PIN permeability have been done for the case when PINs are active only on the PM (see Box 4, [18]), which have shown that:
  P P I N P M - o n l y v P A T .
This relation was proved in theoretical studies using simple mathematical models and confirmed by computational models of multicellular tissues (see Box 4). It is determined under the assumption that PINs are active only on the PMs and facilitate auxin transport without any action from vesicles. Consequently, in epidermal cells of Arabidopsis root P P I N 2 P M 1   μ m / s   .
We conclude from Equations (31) and (32) that if all other parameters are held in physiological range, to yield the same PAT velocity as PM PINs, vesicular PINs permeability has to be much higher: P P I N v - o n l y P P I N P M - o n l y 42 . However, P P I N v - o n l y and P P I N P M - o n l y also depend on density of PINs on the membranes. To clarify this issue, in the next section we calculated the ratio of PINs activity, which is characteristic of individual transporters and which does not depend on the densities of PINs.
Box 4. The permeability of PINs on the plasma membrane P P I N P M equals the directional auxin transport velocity v P A T .
P P I N P M can be derived from the following considerations, analogous to Equations (19)–(21). Auxin flux density due to PINs through the PM has the following expression because PINs are transporting anions IAA and depend only on its intracellular concentration f a c · A C :
Φ P M = P P I N P M · f a c · A C ,
that coincides with Equation (39).
It should be not lower than the PAT flux Equation (19): Φ P M Φ c e l l a p o p l a s t , because other transport mechanisms, diffusion and non-polar transporters, are diminishing polar auxin transport velocity. The relationship between PAT velocity and the polar permeability accounting for the non-polar permeability was derived by Mitchison [17] (Equation (2)), which is discussed below in this box. Thus, an underestimate for P P I N P M is obtained by equating the directional transport to transport by PM-PINs:
P P I N P M · f a c · A C v P A T · A c .
Moreover, non-polar term is at least 100 less than P P I N P M . Consequently, we can consider
P P I N P M · f a c = v P A T ,
or, as f a c 1 , for pH = 7
P P I N P M v P A T .
Thus, for v P A T = 1   μ m / s , P P I N P M 1   μ m / s .
In fact, References [17,19] have shown that for a file of cells transport speed is comparable to the efflux permeability, and this conclusion was confirmed by computer simulations of multilayered tissues [17,18,19]. In these classical works, it has been proven that the “advection” of auxin can be just a result of combined polar membrane transport and cytoplasmic diffusion. Thus, macroscopic advection velocity is limited either by polar membrane transport or by the rate of auxin’s transfer along the cell length; whichever value is lower. One can calculate that diffusion along the longest cell length (~   100   μ m ) is faster than the measured velocity, proving that cytoplasmic transport does not limit the “macroscopic advection velocity v P A T ” [14]. In this case, polar membrane transport governs macroscopic velocity [17,19]. To show this, we provide derivation done by Mitchison [17].
Derivation of P P I N P M that would yield v P A T =   1   μ m / s from theory developed by Mitchison [17]:
Equation (2) from [17] reads:
1 v P A T = 1 p + ( 1 + 2 q p ) L 2 D ,
where v P A T —macroscopic auxin velocity, p —polar efflux permeability, q —non-polar permeability (by diffusion), L—cell length, D—diffusion coefficient of auxin inside the cell. Because L 2 D 1 ×   10 4 m 2 · 6.7   ×   10 10 m 2 s = 7.46 × 10 4   s / m , v P A T = 1 × 10 6   m /s:
p = 1 + 2 q L 2 D 1   v P A T L 2 D 1.1 · 10 6 m / s .
This equation shows that the efflux permeability p   v P A T , because 2 q L 2 D 1 for realistic maximum value of q = 10 6   m / s . Note, that condition 1   v P A T L 2 D 0 must hold, so that the value of permeability is positive, which is true for v P A T 2 D L = 13.4 × 10 6   m / s . The polar permeability value estimated in [17] is also in the order of 10 6 m / s .
In principle, auxin-containing vesicles can also contribute to directional transport within the cell. However, as noted above, diffusion is already sufficient to transport auxin inside the cell, so such additional “acceleration” is not relevant for the macroscopic transport rate.

2.3.1.4. Individual Activity of PINs Has to Be Much Higher on the Vesicles than on the PM to Produce the Same Auxin Transport Velocity

To compare the efficiency of transport, it is necessary to normalize permeabilities P P I N v - o n l y and   P P I N P M - o n l y to the corresponding density of PINs. Membrane permeability due to PINs is proportional to the PINs’ density and can be expressed as a product of 1) density of transporters on the membrane ρ P I N   [ mol μ m 2 ] and (2) individual activity of one transporter protein p P I N [ μ m / s mol / μ m 2 ], which depends on affinity to auxin, electrical potential across the membrane ([19]), phosphorylation status, and any other parameters.
  P P I N = ρ P I N · p P I N .
The ratio of normalized permeabilities equals the ratio of the individual activity of PIN transporters situated in vesicles ( p v - o n l y ) and on the PM ( p P M - o n l y ): p v - o n l y p P M - o n l y = P P I N v - o n l y ρ v   ρ P M P P I N P M - o n l y .
Combining this equation and Equation (29) (equality form):
  P P I N v - o n l y = Κ · v P A T   ρ v ρ P M  
gives:
  p v - o n l y p P M - o n l y = v P A T   · Κ P P I N P M - o n l y ρ v ρ P M   ·   ρ P M ρ v .
Note that density ratio ρ v ρ P M cancels in this equation:
  p v - o n l y p P M - o n l y = v P A T   · Κ P P I N P M - o n l y .
By calculating the minimum of the right-hand side and using the fact that v P A T   = P P I N P M - o n l y · f a c (see Box 4), we can find the lower bound for the PINs activity ratio:
  p v - o n l y p P M - o n l y Κ · f a c  
Introducing the minimal value of K (Equation (30)) and f a c = 0.994 :
p v - o n l y p P M - o n l y >   4200 .
For p H v = 6.2 , (see Table 3), the ratio is even higher: p v - o n l y p P M - o n l y >   26600 .
Conclusion I: For vesicle transport to be able to produce all of the observed IAA flux, the activity of PINs on the vesicle membrane has to be at least three orders of magnitude greater than that estimated for the case when PINs are only active on the PM. Our calculations show that transporting auxin directly through the PM is 1000 times more effective than by means of SVs. This, in our opinion, is an argument against SVs-mediated transport of auxin as a major mechanism of directional auxin transport.

2.3.2. PINs Active Both in Vesicles and on the PM

Next, we derived the expression for the ratio of auxin fluxes through both mechanisms in the case that PINs are active on vesicles and the PM. Auxin efflux density caused by (active) PINs on the PM is:
Φ P M = P P I N P M · f a c · A C .
Auxin efflux density caused by arriving PIN delivery vesicles (from Equation (20) using Equation (16) and omitting the term 1 f a c 1 f a v as in Equation (28)):
    Φ v   = N v · F + P P I N v P d i f f I A A H f a c 1 f a v V v · F + · A c .  
Both effluxes are proportional to intracellular auxin concentration A c , which is assumed to be the same at the PM and around SVs.
Dividing Equation (39) by Equation (40) gives the flux ratio:
  Φ P M   Φ v = P P I N P M · f a c A c P P I N v P d i f f I A A H f a c 1 f a v V v · F + · A c =   P P I N P M P P I N v · P d i f f I A A H ( 1 f a v ) V v · 1 F + .  
Using Equation (33),
Φ P M   Φ v = P d i f f I A A H ( 1 f a v ) V v · 1 F +   ρ P M ρ v p P M p v .
Substituting F + by its expression from Equation (1.7) gives:
Φ P M   Φ v =   P d i f f I A A H · ( 1 f a v ) V v ρ v ρ P M S v   · τ 1 / 2 l n 2 ρ P M ρ v p P M p v ,
where ρ v ρ P M cancel. Substituting with the constant K (Equation (30)) we reduce it further:
Φ P M   Φ v = K   · f a c p P M p v .
The coefficient here is the same as in Equation (38):
Φ P M Φ v 4200   p P M p v ,
where p P M p v is the ratio of individual PINs activity on the PM and vesicles if PINs are active on both membranes. Note that the densities of PINs cancel out and flux distributions between the two mechanisms depend only on the activity of PINs. Also Equation (38) can be derived directly from Equation (45) by assuming Φ P M   Φ v = 1 .
If the individual activity of PINs on the vesicle and on the PM are equal ( p v = p P M ), then,
Φ P M Φ v 4200 ,
which means that the flux through the vesicle mechanism does not contribute more than 0.02% to the total flux.
For p H v = 6.2 the lower bound of the ratio is higher: Φ P M   Φ v 26600 , which implies that the contribution of the vesicular mechanism is less than 0.004%, which is negligibly small compared to PM flux.
The flux ratio Φ P M   Φ v would be higher (favoring PM flux over SV-mediated flux) if:
  • PINs half-life on the PM is higher;
  • pH in the vesicle is lower;
  • vesicles do not stay in the cytoplasm for long enough before fusing with the PM to be fully filled with IAA;
  • the size of vesicles is smaller.
The opposite changes of parameter values would increase the SV-mediated contribution to the auxin flux (see Equation (43)).
Conclusion II: The activity of PIN transporters in the vesicle has to be at least 1000 times greater than on the PM to make a substantial contribution to the total directional auxin transport. When realistic physiological parameter values were used in the model, a factor greater than 10 4 was determined.

3. Discussion

The hypothesis of vesicular secretion of auxin was postulated more than 15 years ago [7,8]. While several studies have claimed to provide experimental evidence to support this concept since [9,10,11], their conclusions are currently greatly debated in the field [12,13]. Hence the question whether vesicular secretion significantly contributes to intercellular auxin transport remains unresolved. Here, we took a modeling approach to estimate whether such mode of auxin transport is even theoretically possible.
First, we compared individual activities of PINs active exclusively on vesicles or on the PM that are necessary to yield the same PAT velocity (Equations (37) and (38)):
p v - o n l y p P M - o n l y P d i f f I A A H ( 6 · τ 1 / 2 d · l n 2 · 1 f a v f a c ) m i n .
This ratio is valid for any PAT velocity and depends on the following measurable parameters: diameter of vesicles d , half-life time of PINs on the PM τ 1 / 2 , pH-dependent fraction of IAAH in a vesicle ( 1 f a v ) , fraction of IAA in cytoplasm f a c . Experimental parameter values (Table 1) gave us p v - o n l y p P M - o n l y >   4200 , meaning that the individual activity of PINs on SVs ( p v ) needs to be at least three orders of magnitude greater than the activity of PINs on the PM ( p P M ) in order to produce the same auxin transport.
We also provide an estimate for the permeability of PINs in the vesicles required to yield v P A T = 1 μ m s in epidermal root cells: P P I N v 42 μ m s , which is not the measure of individual transporters, but characteristic of the unit area of the membrane. This value is much higher than any measured permeabilities to date [18,24]. The measurements of permeability due to PINs either on the PM, in vesicles or in both domains will allow one to draw precise conclusions from our model.
Calculation of permeability P P I N v requires the knowledge of relative density of PINs proteins on vesicles and the PM. We provide measurement of PM density of PINs using SIM microscopy: ρ P M 4 × 10 3 1 μ m 2 , which gives an estimate of minimal possible density ratio: ( ρ v ρ P M ) m i n = 0.01 (Box 2).
Finally, we derived an expression that relates auxin flux driven by PINs on the polar domain of PM and auxin flux driven by PINs at the vesicles, if they are both active (Equation (43), using Equation (30)):
Φ P M   Φ v = P d i f f I A A H · ( 1 f a v ) d · 6   · τ 1 / 2 l n 2   p P M p v .
This shows that relative auxin fluxes through both mechanisms depend on: diameter of vesicles d , half-life time of PINs on the PM τ 1 / 2 , pH-dependent fraction of IAAH in a vesicle ( 1 f a v ) and individual activities of PIN proteins p P M and p v .
The main parameter, namely PIN activity, has not been experimentally determined on the PM or in vesicles, which makes it impossible to find the actual ratio of fluxes in vivo. Nevertheless, by estimating the extreme values of the rest of parameters from available experimental data (Table 1), we were able to show that Φ P M   Φ v > 4200   p P M p v . This expression suggests that if PINs were equally active both at the PM and on the SVs, the contribution of secreted IAA to the overall net IAA flux would be less than 0.02% at pH = 7 inside vesicles of maximal size, which is the maximal possible contribution in that case.
Moreover, instead of the half-life time of PINs on the PM τ 1 / 2 and vesicle diameter, other measurable parameters can be put in this equation (using Equation (1.7), as done in Equation (42)): exocytosis frequency F + , densities of PINs ρ v , P M , pH-dependent fraction of IAAH in a vesicle ( 1 f a v ) and vesicular volume V v :
Φ P M   Φ v = P d i f f I A A H ( 1 f a v ) V v · 1 F + ρ P M ρ v p P M p v .
Interestingly, our theoretical estimate of maximum exocytosis rate Equation (1.7) coincides with experimentally measured maximum endocytosis rate: F m a x + = 0.5 v e s i c l e s   μ m 2   s (see end of Box 1), which provides additional confirmation that the chosen parameter ranges are plausible. Equation (49) shows that the restrictions on frequency of exocytosis is one of the factors that limits the capacity of vesicles to release auxin at a high rate. Also limiting for SV-mediated transport is the fact that the density of PINs on the vesicle membrane cannot be higher than the density of PINs on the PM, which is almost at its possible maximum.
Our assumptions and calculations are based on the model of a PIN2-expressing root epidermal cell. We cannot rule out that vesicular transport of auxin could prevail in different cell types or physiological contexts. However, all parameters in Equations (48) and (49) except PINs activities cannot have zero values, thus, they would have to be changed by orders of magnitude to make vesicular auxin transport the dominant mechanism, causing massive changes in cell homeostasis. Therefore, we consider this unlikely and assume that our conclusions can be extrapolated to any cell type. If experimental values of parameters in any other cell type become available in the future, our equations can be used to resolve this question.
If it was found that PINs transporters are active exclusively in vesicles and are inactive on the PM ( p P M = 0 ), Equations (48) and (49) would become meaningless. Then the testable relation between the individual activity of PINs and the measurable parameters is (from Equation (28)):
p v - o n l y =   P P I N v - o n l y ρ v = P d i f f I A A H ρ P M · v P A T   · 6   · τ 1 / 2 d   · l n 2 ( 1 f a v ) f a c .
However, the individual activity of PINs remains a technically demanding value to determine experimentally. Moreover, SV localization of a member of the AGCVIII kinase family, the action of which is required for PIN activity [30], has not been reported, arguing that PINs in the SVs are possibly inactive [13]. An additional argument that PINs are active at the PM is that increasing PINs incidence at the cell surface, e.g., by inhibiting their endocytosis leads to increased auxin efflux [31,32]. Overall, methods for separate measurement or perturbation of parameters in Equations (48) and (49) are needed to unambiguously distinguish between vesicular-mediated and PM-mediated auxin fluxes.
In conclusion, we have created a simple mathematical model to calculate the efficiency of PIN-mediated vesicular secretion of auxin compared to transport across the PM. Even under the most “vesicular transport-favoring” values of the parameters the vesicular transport of auxin was still determined to be several orders of magnitude less efficient compared to the membrane transport. Our calculations showed that PINs on the PM can produce auxin transport having much less individual activity than required for PINs in the secretory vesicles (SVs). Therefore, we consider it unlikely that PINs are active only in the secretory vesicles; and in case they are active on both PM and SVs, vesicular transport would play a negligible role in PAT.

4. Materials and Methods

Measurement of PIN Density Using Structured Illumination Microscopy (SIM)

Live Arabidopsis roots samples were prepared as previously described by Johnson and Vert [26], where the imaging media was supplemented with 30% opti-prep and the coverslips were fixed on to the microscope slide. Cells in the elongation zone were imaged using an OMX BLAZE v4 3D SIM (Applied Precision). A 60 × 1.42 NA Oil Immersion objective and a 100 mw 488 nm laser was used. The microscope was aligned and calibrated using 100 nm tetra spark beads (Thermo Scientific) and found to have a working resolution of 106 nm and 104 nm in the x and y dimension and 336.9 nm in the z dimension. This SIM system was used due to its enhanced speed of acquisition, when compared to other systems, in order to minimize the effect of rapidly recycling proteins during acquisition. To overcome the curved and non-uniform shape of the roots, Z-stacks (with an optimal z spacing of 125 nm) were created to image the entire lateral membrane of the PIN2 polar domain. This is particularly important when imaging the lateral membrane next to the apical face of the cell, as it curves away from the coverslip. Depending on the shape of the cell, 1.625–3 µm stacks were made of each cell to ensure the entire lateral membrane region was captured, with an average acquisition time for the total z-stack of 17.7 s, with each z position having an average dwell time of 1.08 s. Each z-position image is based on 15 images generated from 3 different angles and 5 different SIM patterns using SOFTWORX software (Applied Precision). A maximum projection of the stacks was used for analysis. Images were made binary and subjected to watershed segmentation using Fiji [33]. PIN2-GFP spots were then detected using TrackMate [34]. The density of PIN2-GFP, mean, and maximum and minimum area of the PIN2-GFP spots were calculated from 12,343 spots in 9 different cells from 3 independent roots.

Author Contributions

Conceptualization, S.H., M.G. and J.F.; Methodology, S.H., M.A., A.J.; Formal Analysis, S.H., M.A.; Investigation, S.H., M.A., M.G., A.J.; Visualization, S.H., A.J.; Funding Acquisition, S.H., J.F.; Writing—Original Draft Preparation, M.A., M.G.; Writing—Review and Editing, S.H., M.A., M.G., A.J. and J.F.

Funding

European Research Council (ERC): 742985 to Jiri Friml; M.A. was supported by the Austrian Science Fund (FWF) (M2379-B28); AJ was supported by the Austria Science Fund (FWF): I03630 to Jiri Friml.

Acknowledgments

We cordially thank František Baluška, Kees Boot, Bert van Duijn, Geneviève Hines, Kees Libbenga and Eric Kramer for their valuable input to the initial discussions that led to this work. These took place at the “Quantitative Biology of Auxin Transport” workshop at the Lorentz Center, Leiden, Netherlands, 12–16 January 2015, and we would therefore also like to acknowledge the organizers of this workshop. We acknowledge the Advanced Microscopy Facility (advMICRO) of the Vienna Biocenter Core Facilities (VBCF), member of the Vienna Biocenter (VBC) Austria, for use of the OMX v4 3D SIM microscope and Lijuan Zhang for her technical assistance. We thank Ivan Kulik for comments on the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

PMPlasma membrane
SVSecretory vesicle
IAAIndole-3-acetic acid = auxin
PINsPIN-formed transporters
TGN/EETrans-Golgi network/ early endosome

References

  1. Adamowski, M.; Friml, J. PIN-Dependent Auxin Transport: Action, Regulation, and Evolution. Plant Cell 2015, 27, 20–32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Kleine-Vehn, J.; Langowski, L.; Wisniewska, J.; Dhonukshe, P.; Brewer, P.B.; Friml, J. Cellular and Molecular Requirements for Polar PIN Targeting and Transcytosis in Plants. Mol. Plant 2008, 1, 1056–1066. [Google Scholar] [CrossRef] [PubMed]
  3. Zhou, J.-J.; Luo, J. The PIN-FORMED Auxin Efflux Carriers in Plants. Int. J. Mol. Sci. 2018, 19, 2759. [Google Scholar] [CrossRef] [PubMed]
  4. Geldner, N.; Friml, J.; Stierhof, Y.-D.D.; Jürgens, G.; Palme, K. Auxin Transport Inhibitors Block PIN1 Cycling and Vesicle Trafficking. Nature 2001, 413, 425–428. [Google Scholar] [CrossRef] [PubMed]
  5. Dhonukshe, P.; Aniento, F.; Hwang, I.; Robinson, D.G.; Mravec, J.; Stierhof, Y.-D.D.; Friml, J.J. Clathrin-Mediated Constitutive Endocytosis of PIN Auxin Efflux Carriers in Arabidopsis. Curr. Biol. 2007, 17, 520–527. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Adamowski, M.; Narasimhan, M.; Kania, U.; Glanc, M.; De Jaeger, G.; Friml, J. A Functional Study of AUXILIN-LIKE1 and 2, Two Putative Clathrin Uncoating Factors in Arabidopsis. Plant Cell 2018. [Google Scholar] [CrossRef] [PubMed]
  7. Friml, J.; Palme, K. Polar Auxin Transport—Old Questions and New Concepts? Plant Mol. Biol. 2002, 49, 273–284. [Google Scholar] [CrossRef] [PubMed]
  8. Baluška, F.; Šamaj, J.; Menzel, D. Polar Transport of Auxin: Carrier-Mediated Flux across the Plasma Membrane or Neurotransmitter-like Secretion? Trends Cell Biol. 2003, 13, 282–285. [Google Scholar] [CrossRef]
  9. Mancuso, S.; Marras, A.M.; Mugnai, S.; Schlicht, M.; Zársky, V.; Li, G.; Song, L.; Xue, H.-W.; Baluska, F. Phospholipase Dzeta2 Drives Vesicular Secretion of Auxin for Its Polar Cell-Cell Transport in the Transition Zone of the Root Apex. Plant Signal. Behav. 2007, 2, 240–244. [Google Scholar] [CrossRef] [PubMed]
  10. Mettbach, U.; Strnad, M.; Mancuso, S.; Baluška, F. Immunogold-EM Analysis Reveal Brefeldin a-Sensitive Clusters of Auxin in Arabidopsis Root Apex Cells. Commun. Integr. Biol. 2017, 10, e1327105. [Google Scholar] [CrossRef] [PubMed]
  11. Schlicht, M.; Strnad, M.; Scanlon, M.J.; Mancuso, S.; Hochholdinger, F.; Palme, K.; Volkmann, D.; Menzel, D.; Baluska, F. Auxin Immunolocalization Implicates Vesicular Neurotransmitter-like Mode of Polar Auxin Transport in Root Apices. Plant Signal. Behav. 2006, 1, 122–133. [Google Scholar] [CrossRef] [PubMed]
  12. Alpi, A.; Amrhein, N.; Bertl, A.; Blatt, M.R.; Blumwald, E.; Cervone, F.; Dainty, J.; De Michelis, M.I.; Epstein, E.; Galston, A.W.; et al. Plant Neurobiology: No Brain, No Gain? Trends Plant Sci. 2007, 12, 135–136. [Google Scholar] [CrossRef] [PubMed]
  13. Robinson, D.G.; Hawes, C.; Hillmer, S.; Jürgens, G.; Schwechheimer, C.; Stierhof, Y.-D.; Viotti, C. Auxin and Vesicle Traffic. Plant Physiol. 2018, 176, 1884–1888. [Google Scholar] [CrossRef] [PubMed]
  14. Kramer, E.M.; Rutschow, H.L.; Mabie, S.S. AuxV: A Database of Auxin Transport Velocities. Trends Plant Sci. 2011, 16, 461–463. [Google Scholar] [CrossRef] [PubMed]
  15. Petrášek, J.; Friml, J. Auxin Transport Routes in Plant Development. Development 2009, 136, 2675–2688. [Google Scholar] [CrossRef] [PubMed]
  16. Rashotte, A.M.; Poupart, J.; Waddell, C.S.; Muday, G.K. Transport of the Two Natural Auxins, Indole-3-Butyric Acid and Indole-3-Acetic Acid, in Arabidopsis. Plant Physiol. 2003, 133, 761–772. [Google Scholar] [CrossRef] [PubMed]
  17. Mitchison, G.J. The Dynamics of Auxin Transport. Proc. R. Soc. B Biol. Sci. 1980, 209, 489–511. [Google Scholar] [CrossRef]
  18. Rutschow, H.L.; Baskin, T.I.; Kramer, E.M. The Carrier AUXIN RESISTANT (AUX1) Dominates Auxin Flux into Arabidopsis Protoplasts. New Phytol. 2014, 204, 536–544. [Google Scholar] [CrossRef] [PubMed]
  19. Goldsmith, M.H.; Goldsmith, T.H.; Martin, M.H. Mathematical Analysis of the Chemosmotic Polar Diffusion of Auxin through Plant Tissues. Proc. Natl. Acad. Sci. USA 1981, 78, 976–980. [Google Scholar] [CrossRef] [PubMed]
  20. Pimpl, P.; Movafeghi, A.; Coughlan, S.; Denecke, J.; Hillmer, S.; Robinson, D.G. In Situ Localization and in Vitro Induction of Plant COPI-Coated Vesicles. Plant Cell 2000, 12, 2219–2236. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Luo, Y.; Scholl, S.; Doering, A.; Zhang, Y.; Irani, N.G.; Di Rubbo, S.; Neumetzler, L.; Krishnamoorthy, P.; Van Houtte, I.; Mylle, E.; et al. V-ATPase Activity in the TGN/EE Is Required for Exocytosis and Recycling in Arabidopsis. Nat. Plants 2015, 1, 15094. [Google Scholar] [CrossRef] [PubMed]
  22. Shen, J.; Zeng, Y.; Zhuang, X.; Sun, L.; Yao, X.; Pimpl, P.; Jiang, L. Organelle PH in the Arabidopsis Endomembrane System. Mol. Plant 2013, 6, 1419–1437. [Google Scholar] [CrossRef] [PubMed]
  23. Jásik, J.; Boggetti, B.; Baluška, F.; Volkmann, D.; Gensch, T.; Rutten, T.; Altmann, T.; Schmelzer, E. PIN2 Turnover in Arabidopsis Root Epidermal Cells Explored by the Photoconvertible Protein Dendra2. PLoS ONE 2013, 8, e61403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Delbarre, A.; Muller, P.; Imhoff, V.; Guern, J. Comparison of Mechanisms Controlling Uptake and Accumulation of 2,4-Dichlorophenoxy Acetic Acid, Naphthalene-1-Acetic Acid, and Indole-3-Acetic Acid in Suspension-Cultured Tobacco Cells. Planta 1996, 198, 532–541. [Google Scholar] [CrossRef] [PubMed]
  25. Swarup, R.; Kramer, E.M.; Perry, P.; Knox, K.; Leyser, H.M.O.; Haseloff, J.; Beemster, G.T.S.S.; Bhalerao, R.; Bennett, M.J.; Leyser, O.; et al. Root Gravitropism Requires Lateral Root Cap and Epidermal Cells for Transport and Response to a Mobile Auxin Signal. Nat. Cell Biol. 2005, 7, 1057–1065. [Google Scholar] [CrossRef] [PubMed]
  26. Johnson, A.; Vert, G. Single Event Resolution of Plant Plasma Membrane Protein Endocytosis by TIRF Microscopy. Front. Plant Sci. 2017, 8, 612. [Google Scholar] [CrossRef] [PubMed]
  27. Konopka, C.A.; Bednarek, S.Y. Comparison of the Dynamics and Functional Redundancy of the Arabidopsis Dynamin-Related Isoforms DRP1A and DRP1C during Plant Development. Plant Physiol. 2008, 147, 1590–1602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Konopka, C.A.; Backues, S.K.; Bednarek, S.Y. Dynamics of Arabidopsis Dynamin-Related Protein 1C and a Clathrin Light Chain at the Plasma Membrane. Plant Cell Online 2008, 20, 1363–1380. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Fendrych, M.; Akhmanova, M.; Merrin, J.; Glanc, M.; Hagihara, S.; Takahashi, K.; Uchida, N.; Torii, K.U.; Friml, J. Rapid and Reversible Root Growth Inhibition by TIR1 Auxin Signalling. Nat. Plants 2018, 4, 453–459. [Google Scholar] [CrossRef] [PubMed]
  30. Zourelidou, M.; Absmanner, B.; Weller, B.; Barbosa, I.C.; Willige, B.C.; Fastner, A.; Streit, V.; Port, S.A.; Colcombet, J.; de la Fuente van Bentem, S.; et al. Auxin Efflux by PIN-FORMED Proteins Is Activated by Two Different Protein Kinases, D6 PROTEIN KINASE and PINOID. Elife 2014, 3, e02860. [Google Scholar] [CrossRef] [PubMed]
  31. Paciorek, T.; Zazímalová, E.; Ruthardt, N.; Petrásek, J.; Stierhof, Y.-D.; Kleine-Vehn, J.; Morris, D.A.; Emans, N.; Jürgens, G.; Geldner, N.; et al. Auxin Inhibits Endocytosis and Promotes Its Own Efflux from Cells. Nature 2005, 435, 1251–1256. [Google Scholar] [CrossRef] [PubMed]
  32. Simon, S.; Kubeš, M.; Baster, P.; Robert, S.; Dobrev, P.I.; Friml, J.; Petrášek, J.; Zažímalová, E. Defining the Selectivity of Processes along the Auxin Response Chain: A Study Using Auxin Analogues. New Phytol. 2013, 200, 1034–1048. [Google Scholar] [CrossRef] [PubMed]
  33. Schindelin, J.; Arganda-Carreras, I.; Frise, E.; Kaynig, V.; Longair, M.; Pietzsch, T.; Preibisch, S.; Rueden, C.; Saalfeld, S.; Schmid, B.; et al. Fiji: An Open-Source Platform for Biological-Image Analysis. Nat. Methods 2012, 9, 676–682. [Google Scholar] [CrossRef] [PubMed]
  34. Tinevez, J.-Y.; Perry, N.; Schindelin, J.; Hoopes, G.M.; Reynolds, G.D.; Laplantine, E.; Bednarek, S.Y.; Shorte, S.L.; Eliceiri, K.W. TrackMate: An Open and Extensible Platform for Single-Particle Tracking. Methods 2017, 115, 80–90. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (AC) Schematics of PIN-formed (PIN) endomembrane trafficking and possible auxin transport mechanisms.
Figure 1. (AC) Schematics of PIN-formed (PIN) endomembrane trafficking and possible auxin transport mechanisms.
Ijms 19 03566 g001
Figure 2. (A) Map of the model steps and usage of results from side-calculations presented in boxes; (B) scheme of the auxin transport into the secretory vesicle (SV) and illustration of variables used in the model.
Figure 2. (A) Map of the model steps and usage of results from side-calculations presented in boxes; (B) scheme of the auxin transport into the secretory vesicle (SV) and illustration of variables used in the model.
Ijms 19 03566 g002
Figure 3. PIN2 translational fusion to Green Fluorescent Protein (PIN2-GFP) density measured by 3D structured illumination microscopy (SIM). (a) Example max projection of a 3D SIM image of the lateral membrane of root epidermal cells; (b) Left panel; the magnified yellow rectangle in (a). Middle panel; the image is made binary and subjected to watershed segmentation. Right panel; pink circles denote detected PIN2-GFP spots. Scale bars; (a), 5 µm, (b), 2 µm.
Figure 3. PIN2 translational fusion to Green Fluorescent Protein (PIN2-GFP) density measured by 3D structured illumination microscopy (SIM). (a) Example max projection of a 3D SIM image of the lateral membrane of root epidermal cells; (b) Left panel; the magnified yellow rectangle in (a). Middle panel; the image is made binary and subjected to watershed segmentation. Right panel; pink circles denote detected PIN2-GFP spots. Scale bars; (a), 5 µm, (b), 2 µm.
Ijms 19 03566 g003
Table 1. Values of parameters used in this study.
Table 1. Values of parameters used in this study.
ParameterDescriptionValue/RangeRef.
d Typical diameter of a vesicle 0.06 0.08   μ m [20]
S v = πd2Surface area of a vesicle 1.1 2.0   × 10 2   μ m 2
V v = ( 1 / 6 ) πd3Volume of a vesicle 1.1 2.7   × 10 4   μ m 3
pHvpH in the lumen of the trans-Golgi network/early endosome (TGN/EE)-derived vesicle5.5–6.5[21,22]
pHcpH in the cytoplasm7[21,22]
v P A T Speed of directional plant hormone auxin (IAA) transport in epidermal cells of Arabidopsis root tip 2 3   μ m / s [14,16]
τ 1 / 2 Half-life of PINs on the plasma membrane (PM) 1.3 × 10 4   s [23]
  ( ρ v ) m i n Minimum PIN density in a vesicle 50 1 μ m 2 One molecule per vesicle
  ( ρ P M   ) m a x Maximum PIN density on a plasma membrane 2 × 10 4 1 μ m 2    
4 × 10 3 1 μ m 2
Close-packing of equal spheres of 4 nm radius.Estimate based on experimental data, this study.
P d i f f ( IAAH ) Diffusional permeability of IAA through membrane 0.5   μ m / s [18,24]
Table 2. List of assumptions *.
Table 2. List of assumptions *.
AssumptionEffect on Efficiency of SVs to Transport AuxinEffect on Efficiency of PM to Transport Auxin
1Steady state/homogeneity: auxin fluxes and parameters that influence auxin fluxes are constant in time. Note, that experimentally measured values of parameters used in the model were also assumed to be constant when measured (Table 1): size of SV, pH in the cell and in SV, and listed below:
1.1PINs density on the PM equals an average density of PINs, being constant and homogenous along the whole polar side of the PM.
1.2PINs density on SVs is constant, all SV are identical and contain the same density of PINs.
1.3Intracellular auxin concentration near the PM is constant and homogenous.
1.4Intravesicular auxin concentration is homogenous.
1.5Removal of PINs from the PM is proportional to the PINs concentration on the PM and is a constant process in time, characterized by constant half-life of PINs on the PM. It doesn’t depend on the exocytosis rate.Less efficientNo effect
2SVs fuse with the PM and deliver all PIN proteins that they contain to the PM.Less efficientNo effect
3SV fills with auxin to its maximum concentration and then all auxin inside the SV is released outside the cell.More efficientNo effect
4Diffusion of anion form (IAA) through the membrane (PM and SV) is negligible.More efficient-
5PINs transport auxin (IAA) only in one direction (inside the vesicle, outside the cell on the PM).More efficientMore efficient
6Non-polar auxin transport is neglected when equating auxin flux through PM or via SVs and directional auxin transport rateMore efficient, magnitude of correction is the same for PM and SVs
* Validation of each assumption is provided in the text; most assumptions favor vesicular transport (increase auxin efflux mediated by SVs) and are taken to estimate minimal requirements for vesicular auxin transport.
Table 3. Dependence of the accumulation ratio R (Equation (18)) and maximum number of molecules in the vesicle N v m a x (Equation (15)) on pH in the vesicle. The last column provides the lower bounds for P P I N v - o n l y P d i f f I A A H , required to yield directional transport of auxin v P A T   depending on p H v . f H v / c —fraction of IAA in the vesicle/cytoplasm, ( 1 f a v / c ) —fraction of IAAH in the vesicle/cytoplasm.
Table 3. Dependence of the accumulation ratio R (Equation (18)) and maximum number of molecules in the vesicle N v m a x (Equation (15)) on pH in the vesicle. The last column provides the lower bounds for P P I N v - o n l y P d i f f I A A H , required to yield directional transport of auxin v P A T   depending on p H v . f H v / c —fraction of IAA in the vesicle/cytoplasm, ( 1 f a v / c ) —fraction of IAAH in the vesicle/cytoplasm.
  p H v     f a v     1 f a v   R *   N v m a x     P P I N v - o n l y P d i f f I A A H  
5.50.8330.167 6.0 · P P I N v P d i f f + 0.036 ** 0.096 · P P I N v P d i f f 2330 s μ m   · v P A T  
6.20.9620.038 26.2 · P P I N v P d i f f + 0.16 ** 0.42 · P P I N v P d i f f 535 s μ m   · v P A T  
6.50.9800.020 49.7 · P P I N v P d i f f + 0.3 ** 0.80 · P P I N v P d i f f 280 s μ m   · v P A T  
7.00.994 = f a c 0.006 = 1 f a c 166 · P P I N v P d i f f + 1 ** 2.7 · P P I N v P d i f f 84 s μ m   · v P A T  
* For lower p H v accumulation ratios R have lower coefficient in front of P P I N v - o n l y P d i f f I A A H because diffusion from the vesicle, that acts against PIN-mediated influx, is higher for lower p H v . ** the right terms are always much smaller than the left terms and can be neglected, because v P A T 1 μ m s and thus P P I N v P d i f f > 84 .
Table 4. Calculation of rate of endocytosis (based on values from [27,28]).
Table 4. Calculation of rate of endocytosis (based on values from [27,28]).
Endocytic MarkerDRP1C-GFP (Dynamin-Related Protein 1C-GFP)CLC-GFP (Clathrin Light Chain-GFP)
foci per µm23.543.48
standard deviation (SD)0.620.55
average lifetime (s)17.719.7
SD8.86.8
foci per model cell (15 × 15 micron)796.5783
SD139.5123.75
Endocytosis events per cell per secondDRP1C-GFPCLC-GFP
Average45.039.7
max105.270.3
min24.824.9
Maximum rate of endocytosis per second per µm20.4670.312

Share and Cite

MDPI and ACS Style

Hille, S.; Akhmanova, M.; Glanc, M.; Johnson, A.; Friml, J. Relative Contribution of PIN-Containing Secretory Vesicles and Plasma Membrane PINs to the Directed Auxin Transport: Theoretical Estimation. Int. J. Mol. Sci. 2018, 19, 3566. https://doi.org/10.3390/ijms19113566

AMA Style

Hille S, Akhmanova M, Glanc M, Johnson A, Friml J. Relative Contribution of PIN-Containing Secretory Vesicles and Plasma Membrane PINs to the Directed Auxin Transport: Theoretical Estimation. International Journal of Molecular Sciences. 2018; 19(11):3566. https://doi.org/10.3390/ijms19113566

Chicago/Turabian Style

Hille, Sander, Maria Akhmanova, Matouš Glanc, Alexander Johnson, and Jiří Friml. 2018. "Relative Contribution of PIN-Containing Secretory Vesicles and Plasma Membrane PINs to the Directed Auxin Transport: Theoretical Estimation" International Journal of Molecular Sciences 19, no. 11: 3566. https://doi.org/10.3390/ijms19113566

APA Style

Hille, S., Akhmanova, M., Glanc, M., Johnson, A., & Friml, J. (2018). Relative Contribution of PIN-Containing Secretory Vesicles and Plasma Membrane PINs to the Directed Auxin Transport: Theoretical Estimation. International Journal of Molecular Sciences, 19(11), 3566. https://doi.org/10.3390/ijms19113566

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop