Skip to main content

Surface aggregation patterns of LDL receptors near coated pits III: potential effects of combined retrograde membrane flow-diffusion and a polarized-insertion mechanism

Abstract

Although the process of endocytosis of the low density lipoprotein (LDL) macromolecule and its receptor have been the subject of intense experimental research and modeling, there are still conflicting hypotheses and even conflicting data regarding the way receptors are transported to coated pits, the manner by which receptors are inserted before they aggregate in coated pits, and the display of receptors on the cell surface. At first it was considered that LDL receptors in human fibroblasts are inserted at random locations and then transported by diffusion toward coated pits. But experiments have not ruled out the possibility that the true rate of accumulation of LDL receptors in coated pits might be faster than predicted on the basis of pure diffusion and uniform reinsertion over the entire cell surface. It has been claimed that recycled LDL receptors are inserted preferentially in regions where coated pits form, with display occurring predominantly as groups of loosely associated units. Another mechanism that has been proposed by experimental cell biologists which might affect the accumulation of receptors in coated pits is a retrograde membrane flow. This is essentially linked to a polarized receptor insertion mode and also to the capping phenomenon, characterized by the formation of large patches of proteins that passively flow away from the regions of membrane exocytosis. In this contribution we calculate the mean travel time of LDL receptors to coated pits as determined by the ratio of flow strength to diffusion-coefficient, as well as by polarized-receptor insertion. We also project the resulting display of unbound receptors on the cell membrane. We found forms of polarized insertion that could potentially reduce the mean capture time of LDL receptors by coated pits which is controlled by diffusion and uniform insertion. Our results show that, in spite of its efficiency as a possible device for enhancement of the rate of receptor trapping, polarized insertion nevertheless fails to induce the formation of steady-state clusters of receptor on the cell membrane. Moreover, for appropriate values of the flow strength-diffusion ratio, the predicted steady-state distribution of receptors on the surface was found to be consistent with the phenomenon of capping.

Background

Endocytosis generally refers to the process by which substances are internalized through the cell membrane. Receptor-mediated endocytosis (RME) is a highly specialized kind of endocytosis in which large protein molecules called receptors project from the cell membrane and couple selectively to ligands such as low density lipoproteins (LDL) [1]. Coated pits and their associated receptors have been studied most extensively in cells grown in culture, and the LDL receptor is the one on which a majority of experimental research has been conducted. The ligand-receptor complexes aggregate in specialized cell-membrane formations called coated pits [25]. When these close, they form vesicles which then transport the ligand-receptor complexes to the interior of the cell. Once the ligand-receptor bindings are separated, the ligands are degraded at the lysosomes. In some experimental systems, including the LDL system, the receptors are returned to the cell surface for further endocytic cycles, and the ligands are degraded at the lysosomes. Experimental results reveal that RME occurs in virtually all eukaryotic cells except the mature erythrocyte [6] and that it also provides an entrance mechanism for cell nutrients such as low density lipoproteins (LDL) [7], hormones and growth factors [8], serum transport proteins and antibodies [9], as well as toxins and lectins [10] and even viruses [11].

The low density lipoprotein (LDL) particles are cholesterol-transporting macromolecules that are produced in the liver and circulate in the plasma. Membrane reconstitution in human fibroblasts requires the assimilation of cholesterol. The efficiency of this process depends on the rate of endocytosis of the LDL macromolecule and its receptor. It is thought [12] that a severely depleted number of LDL receptors promotes high levels of circulating cholesterol because LDL internalization requires the receptor-binding stage. Deficiencies in the LDL cycle are known to be responsible for the ailment known as familial hypercholesterolemia, which promotes atherosclerosis and the incidence of strokes and coronary disease [13]; this being the reason why so much experimental research has been conducted on the internalization of the LDL receptor. A substantial portion of research aimed at the elucidation of the cholesterol-uptake process has addressed the characterization of the dynamics and display of the LDL receptor. This has produced a large pool of data and a sound conceptual framework which sustains the theoretical exploration of aspects of the dynamics mentioned above.

It is generally thought that the maintenance of receptors on the cell surface is due primarily to receptor recycling rather that to de novo synthesis. These data come from experiments using cycloheximide to block protein synthesis. In this experimental setup, it is observed that the number of LDL receptors on the cell surface remain roughly constant for at least six hours [7]. Evidence of receptor internalization and reinsertion in unblocked systems would sustain the assumption that a steady-state concentration of receptors is maintained at the cell surface. The time receptors spend in the interior of the cell is negligible. The basis of this claim is the apparently undetectable pool of receptors inside the cell during endocytosis [14]. The transit time for an LDL receptor from binding on coated pits to reappearance in the membrane, and found it to be on the order of 15 seconds [14]. Based on these ideas, we will abide by the assumption that internalization and recycling of LDL receptors in human fibroblastic cells maintains the surface concentration of LDL at a steady state. Moreover, coated pits include 1% of the cell surface (coated pits include 2% of the cell surface at 4°C [15, 16], but when the temperature is raised to 37°C the number of coated pits on the surface is halved [17]). Anderson et al. [17] assert that coated pits tend to be linearly aligned over intracellular fibers. Hence, it can be assumed that the geometrical arrangement of coated pits on the cell surface of human fibroblasts can be reasonably approximated by means of a dilute and ordered system of sinks distributed on a two dimensional surface [18]. Guadorov et al. [19] reported that coated pits tend to form repeatedly at defined sites while other regions are excluded. Therefore, we can assume that coated pits are infinitely long-lived traps distributed in a dilute and ordered form over the cell surface.

It is also known that on cultured human fibroblasts, receptors for certain ligands (e.g. insulin, epidermal growth factor and  - 2-macroglobulin) cluster in coated pits only after exposure to the ligand [20], while receptors for LDL cluster in coated pits independently of ligand binding [7]. This feature of the LDL receptor pathway makes it a particularly attractive candidate for mathematical modeling, since initially we can ignore the details of the ligand-receptor interaction and still study the recycling of the receptor and its interaction with the coated pit. Preliminary mathematical models of RME are aimed at calculating the rate at which diffusing particles (receptors), which have been inserted uniformly all over a certain 2-dimensional space (the cell surface) hit traps (coated pits). This rate is known as the diffusion-limited forward-rate constant [21] and is denoted here by k d +. This rate can be calculated as the flux of particles into a trap, divided by the mean particle concentration [22]. In the two-dimensional case, for a circular sink of radius a, k d + is defined by means of the equation

k d + = 2 πaD C C r | r = a ,
(1)

where C(r) is the steady-state radial distribution function of receptors not bound to coated pits, with D > 0 their diffusion coefficient and 〈C〉 the receptor concentration averaged over all the diffusion space [21], that is,

C = 1 π b 2 0 b 2 πrC r dr .
(2)

The constant k d + times the number of traps per unit area gives the probability per unit time that a diffusing particle hits the trap. It is linked to the mean time τ d for a particle to hit a trap (mean capture time) [21] by the relation,

τ d = 1 k d + ρ ,
(3)

where ρ is the number of coated pits distributed per unit area on the cell surface.

If LDL receptors are trapped in times less than, or comparable to the average lifetime of a coated pit (≤5 min), the mean capture time is close to what it would be if the traps were infinitely long-lived. Moreover, since we are dealing with a dilute system of traps, we can follow the method of Berg and Purcell [23] and represent the real multiple-trap problem by means of single coated pit of radius a, surrounded by a reference annulus Ω of outer radius b, i.e.

Ω = r , θ | 0 α 2 π , a r b ,
(4)

where b is calculated through

ρ = 1 π b 2
(5)

Receptors are projected into the reference annulus Ω according to an insertion rate function S(r), and diffuse afterwards until they reach the boundary of the sink at r = a and are trapped. A steady-state concentration of diffusing particles will be maintained if the number of particles lost to the trap is the same as the number inserted per unit time. Then, if C ds (r) denotes the steady-state receptor concentration at a distance r from the center of the sink we will have

D 2 C ds r + S r = 0 for a r b ,
(6)

so that the inner boundary at r = a absorbs, i.e.

C ds a = 0
(7)

and the outer boundary at r = b reflects, i.e.

C ds r = 0 at r = b ,
(8)

that is, there is no net flux of receptors across the outer boundary of Ω. Solving for C ds (r) we obtain

C ds r = a r z b uS u du Dz dz .
(9)

Denoting by means of C du (r) the form of C ds (r) obtained by setting S(r) = S, and with S being a constant, this means that whenever receptors are inserted uniformly over the entire reference annulus Ω we obtain

C du r = b 2 S 2 D ln r a - S r 2 - a 2 4 D
(10)

and if k du + denotes the associated forward rate constant we have,

k du + = 2 πaD C du C du r | r = a .
(11)

Then, for this case equations (1) through (3) yield

τ du = b 4 ln b a 2 D b 2 - a 2 - 3 b 2 - a 2 8 D .
(12)

Goldstein et al. [24] assert that whether traps are ordered or disordered makes only a small difference in the rate at which a trap captures receptors, given that the traps are distributed over the entire surface of the cell in a dilute way. Goldstein et al. [21, 25] dealt with two conflicting hypotheses about the coated-pit recycling process. In either case, the coated pit effectively has a finite lifetime. However, they found that the model of equations (6)-(12) obtained under the assumption that sinks have infinite lifetimes give good approximations for k d + in the experimental system for receptors of low density lipoproteins on human fibroblastic cells. Nevertheless, they also concluded that the experiments do not rule out the possibility that the true rate of accumulation of LDL receptors in coated pits could be faster than predicted on the basis of pure diffusion and uniform reinsertion all over the cell surface. A mechanism that has been proposed by experimental cell biologists which might affect the accumulation of receptors in coated pits is convection. Evidence of a membrane flow comes from the observation that, when the leading edge of a moving cell comes in contact with a small particle, the particle often adheres to the cell surface and is transported backward [2628]. The motion of such adherent particles is consistent with particles undergoing Brownian motion in the presence of a constant membrane flow [29, 30]. It has been hypothesized that this flow originates when membrane components like specific receptors and lipids, which had been internalized from the cell surface, are recycled to the leading edge and contribute to the extension of the cell boundary [29]. Considering that a rapid rate of membrane internalization is observed in endocytosis, a significant amount of recycled membrane should be made available for delivery at the leading edge. This delivery would be restricted to just those molecules that are internalized by coated pits [14, 31, 32], thus initiating the transport of membrane components away from the leading edge and towards the rear of the cell [2631]. In addition to the specific delivery of receptors at the leading edge of the cell, their sorting within specific targets along the endocytic pathway have been proposed as components of a mechanism for the control of the internalization of specific ligands [18, 3335]. Generally, eukaryotic cells show the ability to efficiently control the traffic of many types of macromolecules, targeting proteins and lipids for a variety of destinations, based on their types and the current needs of the cell. This traffic control is crucial for the maintenance of cellular structure and function [36, 37]. This makes it reasonable to assume that reinsertion of LDL receptors could be arranged in such a way that their trapping rate might be adapted to suit specific metabolic requirements. Moreover, it has been claimed [38] that in giant HeLa cells, receptors are not inserted uniformly into the plasma membrane but rather at the periphery of the cell. Robenek and Hesz [39] claimed that their experiments with LDL particles bound to colloidal gold provided the first clear demonstration of the sequential clustering of receptors near coated pits in plaques or loosely associated groups, and not in the form of widely dispersed individual units. They concluded that this effect is produced when recycled LDL receptors are inserted preferentially in membrane regions where coated pits form. According to Wofsy et al. [40] plaques could be conceptualized as receptor clusters surrounding coated pits. The idea of a retrograde flow is essentially connected to the concept of a polarized reinsertion of receptors. However, neither the influence of this insertion mode on the aggregation rate of LDL receptors in coated pits, nor the resultant steady-state surface aggregation patterns have been formally studied.

Methods

For the aims of the present research we extend the formal framework presented in Echavarria-Heras et al. [41] for the purpose of including diffusion aided by a retrograde flow, plus polarized receptor insertion as a combination of mechanisms that influences both the mean capture time τ du and the surface aggregation patterns of unbound LDL receptors. First, we conceptualize a partitioned receptor-insertion mode in which the cell membrane is divided into disjoint regions, each one receiving a fraction of the total number of recycled receptors. In the general form of this mode, different numbers of receptors are projected over the different regions partitioning the cell surface. Polarized insertion, moreover, corresponds to the case in which receptors are recycled into two disjoint regions having abruptly contrasting rates. This partitioned mode can be used to obtain arrangements—which can be seen as special cases—comprising uniform [7], locally uniform, peripheral [38], and plaque forms [39, 40], as well as a generalized radially symmetric mode of the form described in [41].

In order to formalize the idea of a partitioned insertion mode, let’s consider the reference annulus Ω = {(r, θ)| a ≤ r ≤ b, - π ≤ θ ≤ π}, letting I Ω be the number of particles inserted in Ω, m and a being real numbers satisfying 1 ≤ m ≤ b/a, 0 ≤ α ≤ π, we define the disjoint regions Ω c m , α , Ω p m , α and Ω q m , α through

Ω c m , α = r , θ a r b , - π + α θ π - α ,
(13)
Ω p m , α = r , θ ma r b , π - α θ π + α ,
(14)
Ω q m , α = r , θ a r ma , π - α θ π + α .
(15)

Then, we have

Ω = Ω c m , α Ω p m , α Ω q m , α .
(16)

In order to characterize partitioned receptor insertion, we conceive an insertion rate function S (c, p, q, m, α)  that distributes receptors over the disjoint regions Ω c (m, α) Ω p (m, α), and Ω q (m, α) (Figure  1) through the respective non-negative and continuous functions S c m , α , S p m , α and S q m , α . Formally,

Figure 1
figure 1

The geometry of the model. a) A circular trap of radius a (the coated pit) is encircled by an annulus of radius b (the reference region Ω associated with a coated pit). LDL receptors originally inserted at a point (r, θ) inside the reference annulus Ω, move afterwards by convection and diffusion until they are trapped in coated pits. b) Receptor insertion occurs according to a partitioned insertion rate function S (c, p, q, m, α), which sorts receptors at distinct rates S c m , α , S p m , α and S q m , α linked respectively to the disjoint regions Ω c , Ω p and Ω q .

S c , p , q , m , α = S c m , α r , θ Ω c m , α S p m , α r , θ Ω p m , α S q m , α r , θ Ω q m , α ,
(17)

with

Ω c S c m , α + Ω p S p m , α + Ω q S q m , α = I Ω .
(18)

The formal properties of the partitioned insertion mode are elaborated in Appendix A. Table  1 summarizes the possible forms that S (c, p, q, m, α) can achieve. Moreover, S (c, p, q, m, α)  can be handily characterized by using the functions δ c (m, α), δ p (m, α) and δ q (m, α), which respectively stand for the proportion of I Ω that S c m , α , S p m , α and S q m , α insert over the disjoint regions Ω c m , α , Ω p m , α and Ω q m , α . For a and m fixed, equations (A4) through (A6) define the proportions δ c (m, α), δ p (m, α), and δ q (m, α). Whenever a pairwise combination of insertion proportions is vanishing, the associated form of the partitioned mode is polarized. Otherwise, the partitioned mode would be uniform all over Ω or might be locally uniform, thus sorting receptors at dissimilarly constant rates over Ω c m , α , Ω p m , α and Ω q m , α . Although, for ease of presentation, we here partitioned Ω into three disjoint regions, the form given by equation (18) could be extended to include a greater number of disjoint regions forming Ω.

Table 1 Different forms of the partitioned insertion rate function S (c, p, p, m, α) : c , p and q stand for positive constants

Let C λmα (r, θ) denote the steady-state concentration of receptors in Ω. In Appendix B we review the formal steps which establish that if receptors are inserted in Ω in accordance with the insertion rate function S (c, p, q, m, α) and afterwards move by diffusion aided by a retrograde membrane flow, then if λ = v/2D, where v is the strength of the retrograde flow and D stands for the diffusion coefficient (cf. Eq. B1), we will have

C λmα r , θ = G λmα r , θ e λ r cos θ ,
(19)

where G λmα (r, θ)  is the solution to the boundary-value problem,

2 G λmα r , θ - λ 2 G λmα r , θ = - S c , p , q , m , α D e - λb cos θ ,
(20)
G λmα a , θ = 0 ,
(21)
G λmα r , θ = G λmα r , - θ ,
(22)
G λmα r , θ periodic in θ ,
(23)
- π π G λmα r , θ r | r = b - λ cos θ G λmα b , θ e - λb cos θ = 0 .
(24)

The associated forward rate constant is denoted by means of k  + and is given by

k λmα + = π b 2 a b - π π S c , p , q , m , α rdrdθ a b - π π C λmα r , θ rdrdθ ,
(25)

and the associated mean capture time τ λmα is given by

τ λmα = a b - π π G λmα r , θ e λr cos θ rdrdθ a b - π π S c , p , q , m , α rdrdθ .
(26)

For particular characterizations of S (c, p, q, m, α), the resulting patterns of receptors can be simulated. These representations can be obtained by means of the computer graphics technique of ray tracing [42] using grey tones corresponding to the values of C λmα (r, θ).

Results and discussion

Anderson et al. [7] proposed that LDL receptors in human fibroblasts are inserted at random locations and are transported by diffusion toward coated pits, and Goldstein et al. [21, 25] have studied theoretical aspects of the dynamics of the LDL receptor on the cell surface. One of the questions addressed in these models was whether or not the random insertion of LDL receptors into the plasma membrane, followed by pure diffusion and using the measured diffusion coefficient for LDL receptors, could give a rapid enough aggregation of receptors in coated pits to account for the observed rate of removal of LDL from the cell surface. Since an experimentally determined lower bound for the forward rate constant k du + is (2.3 × 10- 10 ± 1.6) × 10- 10 cm 2/s and the value obtained from their model is 1.9 × 10- 10 cm 2/s, the answer to that question is that the hypothesis of random insertion and pure diffusion of LDL receptors to coated pits is consistent with experimental observations but just barely. Therefore, even though the diffusion-limited forward-rate constant k du + calculated by these authors is consistent with available experimental results on the rate of LDL internalization, these do not rule out the possibility that the true value of the mean capture time of receptors is actually smaller than τ du , the value predicted on the basis of random insertion of receptors followed by pure diffusion to coated pits. Indeed, if LDL receptors on human fibroblasts are not inserted uniformly into the plasma membrane, but are inserted preferentially into specialized regions where coated pits form, and if the resulting reduction in the time required for LDL receptors to diffuse to coated pits was significant, the conclusion that diffusion of LDL receptors to coated pits is the limiting step in the interaction of cell surface LDL receptors with coated pits could be wrong. Also, diffusion could be aided by a more active transport in order to increase the rate of accumulation of receptors in coated pits. In particular, we have studied here the effects of a polarized insertion of receptors as well as those related to a transport mechanism set by diffusion aided by a retrograde membrane flow.

Experimentally determined values for coated pit radius a and the outer radius b of the reference region Ω are respectively a = 0.10 ± 0.05 μm [43] and b = 1.0 ± 0.2 μm; these were obtained by using equation (5) with the calculated value for the suitable density of coated pits ρ = 0.31 ± 0.09/μm 2 [16]. In order to typify values of the fundamental ratio λ = v/2D (cf. Eq. B1) we will use a low reference value v = v 0 for the strength v of the retrograde membrane flow. The value v 0 = 9.6 × 10- 7 cm/min represents the maximum of a hypothetical local flow that originates when coated pits invaginate to internalize the trapped ligand-receptor complexes [44]. As a second reference value for the strength v of the retrograde flow we may take v 1 = 10- 4 cm/min. This value gives the approximate rate at which objects on the surface of fibroblasts are swept backwards [26, 14]. Since v 1 is approximately equal to 100v 0, this retrograde membrane-flow rate would sustain an extremely fast convective transport in comparison with one having a strength v 0. In order to characterize a variation range for the diffusion coefficient D, we will consider a normal process for the LDL receptor. This is associated with a diffusion coefficient value of D 0 = (2.7 ± 0.09) × 10- 9 cm2/min, which was determined by Barak and Webb [45] using fluorescence photo bleaching recovery (FBR) in an experimental arrangement involving LDL receptors diffusing on the surface of human fibroblast cells. We can also contemplate a fast diffusion process corresponding to D 1 = 8.1 × 10- 8 cm2/min, this being the average of the values defining the observed range of variation for the diffusion coefficient of the dil (3)- LDL receptor complex on blebs induced in the cell membrane. And we could also consider an extreme reference value of D ext  = 6 × 10- 7 cm 2/min, that is D ext  = 7.40D 1.

The phenomenon of capping is characterized by the formation of large patches of proteins that passively flow away from the regions of membrane exocytosis. Individual proteins can overcome this convective transport when they have diffusion coefficients (D) of about D ext [14]. Moreover, a value for D of D ext /6 is comparable to the largest estimation for the diffusion coefficient of the dil (3)- LDL receptor complex on blebs [45]. For the parameter values of the LDL system on human fibroblasts, the calculated value for τ du , the mean capture time when receptors move by pure diffusion and are inserted uniformly over the entire cell surface, is τ du  = 2.92 min (cf. Eq. (11)).

In the present setting, a uniform insertion-rate function can be acquired by, for instance taking α = 0, m = b/a and S c m , 0 = c , where c is a positive constant and S p m , 0 = S q m , 0 = 0 . This representation is known here as a c-uniform insertion mode (Figure  2 and Eq. A22). A second description of a uniform insertion mode is what we call a pq-uniform rate function, in which receptors are inserted solely over the regions Ω p (m, π) and Ω q (m, π); this is realized by setting S p m , π = S q m , π = c , where c is a positive constant (Figure  3 and Eq. A24). Finally, a description that we refer to as cpq-uniform insertion is achieved by sorting receptors over the three disjoint regions Ω c (m, α), Ω p (m, α), Ω q (m, α)  with S c m , α = S p m , α = S q m , α = c , where c is a positive constant (Figure  4 and Eq. A26). Consistently, for the case where no convective flow is present (v = 0) and receptors diffuse normally (D = D 0), all these representations indistinctively yield τ λmα  = τ du .

Figure 2
figure 2

The c -uniform insertion mode. A receptor insertion mechanism symbolized here by means of S (c, 0, 0, b/a, 0), which is obtained by setting S p b / a , 0 = S q m , α = 0 and S c b / a , 0 = c , with c a positive constant.

Figure 3
figure 3

The pq -uniform insertion mode. This is denoted by means of the symbol S (0, c, c, m, π), and obtained by setting the conditions S c m , π = 0 and S p m , π = S q m , π = c , with c a positive constant.

Figure 4
figure 4

The cpq -uniform insertion mode. A device symbolized here using S (c, c, c, m, α) and is obtained by setting S c m , α = S p m , α = S q m , α = c , with c a positive constant.

For the case in which convective transport is moderate (v = v 0) and receptors diffuse normally (D = D 0), that is, the flow strength to diffusion ratio λ has a value λ = v 0/2D 0, calculating τ λmα as given by equation (B27) for S (c, p, q, m, α) where insertion is uniform over all of Ω, we obtained τ λmα  = 1.003 τ du . Hence, the present model predicts that if LDL receptors diffuse normally, a retrograde membrane flow having strength comparable to v 0 would have an insignificant effect on τ du . But for that same convection-to-diffusion ratio, we found that a locally uniform insertion mode can prompt greater receptor trapping rates than those linked to uniform insertion all over Ω. But now suppose that receptor insertion takes place over the regions Ω c (m, α), Ω p (m, α) and Ω q (m, α), at different constant rates S c m , α = c , S p m , α = p and S q m , α = q ; this results in what we term here a cpq-locally uniform insertion mode (Figure  5 and Eq. A28). Then, if for instance we choose m = 1.5, α = π/2, δ c (m, α) = 0.10, δ p (m, α) = 0.30 and δ q (m, α) = 0.60, the model yields τ λmα  = 0.49τ du , which amounts to an important reduction of τ du . For α = π we obtain what we label here as a pq-locally uniform insertion form (Figure  6 and Eq. 30) denoted by the symbol S (0, p, q, m, π). Receptors are inserted in a radially symmetric manner over each of the regions Ω q (m, π) and Ω p (m, π), but through different constant rates S p m , α = p and S q m , α = q . For m = 2.0, δ p (m, π) = 0.20 and δ q (m, π) = 0.80, the model yields τ λmα  = 0.41τ du . Now, since for a pq-locally uniform insertion form we have δ p (m, π) + δ q (m, π) = 1, then letting δ p (m, π) approach zero, receptor insertion will be gradually accommodated within the region Ω q (m, α)  so that eventually, when δ p (m, π) vanishes, we will obtain what we identify as a q-plaque form insertion mode denoted by means of the symbol S (0, 0, q, m, π) (Figure  7 and Eq. A32). This is actually the plaque-form insertion mechanism envisioned by Wofsy et al. [40] for modeling preferential insertion as conceived by Robeneck and Hesz [39]. For λ = v 0/D 0, and m = 2.0, the mode S (0, 0, q, m, π) yields τ λmα  = 0.26τ du , again a noticeable reduction in τ du . This mode is radially symmetric and can be already considered as a form of polarized insertion. Similarly, if we initially chose a locally uniform mode S (0, p, q, m, π), then letting δ q (m, π) approach zero we will force receptors to be mainly sorted over the Ω p (m, π), region and in due course when δ q (m, π) vanishes will produce what we call a p-peripheral insertion mode symbolized by means of S (0, p, 0, m, π) (Figure  8 and Eq. A34). This insertion rate function is a peripheral form of polarized insertion; receptors are inserted in an annulus contiguous to the outer boundary of Ω. For λ = v 0/D 0, peripheral insertion yields τ λmα  ≥ 1.003 τ du with the lower bound of 1.003 τ du attained in the limiting case when m approaches one. Hence, if enhancement of LDL receptor trapping rate is required, a peripheral insertion mode for m > 1 turns out to be an inefficient mechanism. For instance taking λ = v 0/2D 0, and m = 9.7 we obtain τ λmα  = 1.15τ du .

Figure 5
figure 5

The cpq -locally uniform insertion mode. This insertion rate function is labeled by means of S (c, p, q, m, α), it is linked to the conditions S c m , α = c , S p m , α = p , and S q m , α = q , with c, p and q positive constants.

Figure 6
figure 6

The pq -locally uniform insertion mode. A receptor insertion mechanism which is denoted by means of S (0, p, q, m, π). It is radially symmetric over each one of the regions Ω p and Ω q . This form is associated with the case S c m , π = 0 , S p m , π = p , and S q m , π = q , with p and q positive constants.

Figure 7
figure 7

The q -plaque form insertion mode. This mode is denoted here by means of the symbol S (0, 0, q, m, π). It is radially symmetric and polarized and associated with the conditions S c m , π = 0 , S p m , π = 0 , and S q m , π = q , with q a positive constant.

Figure 8
figure 8

The p -peripheral insertion mode. An insertion form denoted here by means of S (0, p, 0, m, π). It is a radially symmetric and polarized and linked to the case S c m , π = 0 , S p m , π = p , and S q m , π = 0 , with p a positive constant.

We can deal with three different non-radially symmetric polarized insertion forms. One is obtained if we let S c m , α = 0 and S q m , α = 0 , then all recycling receptors will be sorted over the region Ω p (m, α) and we will be dealing with a p-polarized receptor insertion device denoted by means of S (0, p, 0, m, α) (Figure  9 and Eq. A36). Furthermore, for a p-polarized mode, by reducing the area of the insertion region Ω p (m, α), we could accommodate all recycled receptors in a favorable position relative to the cross section of the coated pit perpendicular to flow streamlines. Then, we may conjecture that this will allow the effect of convection to be maximal for receptor trapping rate enhancement. But we found that this device could at best yield a narrow variation range of τ λmα relative to τ du and not a relevant reduction in this value. In fact, for a p-polarized insertion mode with 0 < α ≤ π/6 and 1.1 ≤ m ≤ b/a, we found 1.01τ du  ≤ τ λmα  ≤ 1.16τ du . Now, if we assume that 0 < α < π, and also that S c m , α = 0 hold, then receptors will be sorted over the regions Ω p (m, α) and Ω q (m, α), resulting in another form of a non-radially symmetric polarized insertion mode. This is regarded here as a pq-polarized receptor insertion form and denoted through S (0, p, q, m, α) (Figure  10 and Eq. A38). In particular, for m = 2.2, α = π/2, δ p (m, α) = 0.20 and δ q (m, π/2) = 0.80, we calculated τ λmα  = 0.45τ du . That is, a pq-polarized receptor insertion form could lead to a significant reduction of τ du . But if we now set S c m , α = 0 and S p m , α = 0 , then all recycling receptors will be delivered over the region Ω q (m, α), and we will have a third form of non-radially symmetric polarized insertion mode. This is named q-polarized insertion and is symbolized by means of S (0, 0, q, m, α) (Figure  11 and Eq. A40). Additionally, this arrangement along with a slow convective transport (v = v 0) and a normal diffusion process (D = D 0)  could potentially induce a major reduction in τ du . In fact, if we set m = 2.0 and α = π/6 we obtain τ λmα  = 0.26τ du , which coincides with the plaque form insertion mode mentioned above. Hence, for the case λ = v 0/D 0 either polarized insertion over the regions Ω p (m, α) and Ω q (m, α) or over the region Ω q (m, α), both seem to provide an efficient mechanism for the reduction of τ du .

Figure 9
figure 9

The p -polarized insertion mode. This form is symbolized here by S (0, p, 0, m, α) and obtained by setting S c m , α = 0 , S p m , α = p , S q m , α = 0 , with p a positive constant. Under the condition λ = v 1/2D 0, which determines capping, this mode was found to produce values for τ λma that are equivalent to τ du .

Figure 10
figure 10

The pq -polarized insertion mode. This insertion mechanism is denoted here through S (0, p, q, m, α). It is polarized and non-radially symmetric and linked to the conditions S c m , α = 0 , S p m , α = p ,and S q m , α = q , with p and q as positive constant. Even for λ = v 1/2D 0, which produce capping-like effects, certain characterization of this mode can give substantial reductions in τ du .

Figure 11
figure 11

The q -polarized insertion mode. A receptor insertion paradigm denoted here by means of S (0, 0, q, m, α) is linked to the case S c m , α = 0 , S p m , α = 0 , and S q m , α = q , with q a positive constant. This mode is non-radially symmetric and polarized. Setting = π , the q-polarized insertion mode gives the paradigm of Wofsy et al. [35] for insertion in plaques. Even when λ = v 1/2D 0 holds so that the capping phenomenon could create a graduated distribution of unbound receptors in the direction of flow streamlines, this mode seems to provide a highly efficient form of receptor insertion by making dramatic reductions on τ du .

If we now assume that λ = v 1/2D 0, that is, a relatively fast convective transport (v = v 1) and a normal diffusion process (D = D 0), and if we again choose the q-polarized insertion mode S (0, 0, 1, 2, π/6 ) described above, then we will get τ λmα  = 0.83 τ du . This result endorses the view that whenever convection becomes a relatively more active transport than diffusion, the process will fail to induce an effective randomization in the surface distribution of receptors [14]. The relative dominance of convection over diffusion implies that those receptors favorably distributed in an area determined by the cross section of the coated pit perpendicular to flow streamlines could be trapped, while the others would be swept away by the flow. Therefore, whenever convection is relatively fast (v = v 1) we may assume that increasing D, that is, taking D > D 0 could induce relatively greater aggregation rates in coated pits. But since we have D 0 = (2.7 ± 0.09) × 10- 9 cm2/min, then the largest feasible experimentally determined value for D 0 is only D = 2.7910- 9 cm2/min, and this along with the S (0, 0, 1, 2, π/6 ) mode and a fast convective transport will produce τ λmα  = 0.75 τ du , which only amounts to a moderate reduction of the value τ λmα  = 0.83 τ du calculated above. Further, if convection is relatively fast (v = v 1) and we still have a q-polarized insertion mode S (0, 0, 1, 2, π/6 ), then unless we have D = 1.74D 0, the calculated value τ λmα will be significantly reduced: for that arrangement we obtained τ λmα  = 0.24 τ du . Therefore, whenever convection is relatively fast (v = v 1), even for an efficient insertion mode such as the q-polarized form S (0, 0, 1, 2, π/6), larger values than D = 1.74D 0 will be required for a noticeable reduction in τ du . Indeed, for a fast retrograde flow (v = v 1) and a fast diffusion process (D = D 1), we found that S (0, 0, 1, 2, π/6 ) can produce a value of τ λmα  = 0.01τ du , which is a drastic reduction of τ du . If we assume, moreover that receptors are returned to the cell surface through a p-polarized insertion mode with 0 < α ≤ π and with 9.5 ≤ m ≤ b/a, then even if convection is extremely fast (v = 300v 1), letting D = D ext , we get 0.19τ du  ≤ τ λmα  ≤ 0.65τ du . Compared with the results obtained above for the same insertion mode, this last inequality highlights the importance of diffusion in the reduction of the mean capture time of LDL receptors by coated pits. However, comparing with the expected variation range for D 0, D values beyond D = 2.7910- 8 cm2/min seem to be extremely high, and likely not achievable in the experimental system analyzed. Thus, before a faster diffusion process can be raised as an agent that can induce an effective randomization of the surface distribution of LDL receptors in the presence of a fast retrograde membrane flow, we assume instead that there must be another factor that could be adduced as a feasible device for promoting higher aggregation rates of LDL receptors in coated pits. And so, in accordance with the results presented here, we can suggest special forms of the insertion rate function S (c, p, q, m, α) as mechanisms that potentially could boost the rate of removal of receptors by coated pits. Thus under the assumption of a fast convective transport, v = v 1, and a normal diffusion process, if we allow returning receptors to be sorted by means of a q-polarized insertion form by setting m = 1.3 and α = π/2, we would get τ λmα  = 0.23τ du , implying a noticeable reduction in τ du . This means that whenever convection is fast and receptors diffuse normally, increasing the receptor-trapping rate could depend decidedly on polarized receptor insertion at preferential regions near coated pits.

For the insertion-transport mechanism for LDL receptors under consideration, we also address here the theoretical exploration of their consequent display on the cell surface. A capping-like surface display should manifest as a graduated concentration of receptors in the direction of flow streamlines. Figure  12 displays the surface aggregation patterns of receptors that result when a polarized reinsertion mode is considered. The combination of polarized insertion, fast convective transport of rate v = v 1 and a relatively slow diffusion process such as one associated to D 0 can be observed to produce a marked gradient in the distribution of unbound receptors (Figure  12a and b). This is consistent with the conceptual model for the capping phenomenon presented by Bretscher [14]. Further, and in concurrence with that author, the present model explains this gradient as being induced by a relative dominance of convection over diffusion; that is, an arrangement which precludes an effective randomization of the surface-receptor distribution. Thus, in the presence of a retrograde membrane flow with typical strength v = v 1, even single LDL particles having a diffusion coefficient value D = D 0 —which is considered normal—would produce a capping-like cell surface distribution; thus precluding both a uniform distribution and the display of clusters of unbound receptors on the surface. Even in the presence of a fast convective transport, a suitably fast diffusion process can reverse the capping effects (Figures  12c and d).

Figure 12
figure 12

Surface aggregation patterns of unbound LDL receptors associated with different values of the fundamental ratio λ and non-radially symmetric-polarized insertion modes. a) the surface pattern made by λ = v 1/2D 0 and the q-polarized insertion mode S (0, 0, q, m, α), with m = 2.3, δ q (m, α) = 1 and α = π/4; b) surface pattern associated with the case λ = v 1/2D 1 and the p-polarized insertion mode S (0, p, 0, m, α) with α = π/6, m = 9.5, and δ p (m, α) = 1; c) the surface aggregation pattern formed by λ = v 1/2D ext and q-polarized insertion S (0, p, q, m, α) with m = 1.2, α = π/6 and δ q (m, α) = 1; d) the surface receptor pattern for λ = v 1/2D ext and p-polarized insertion S (0, p, q, m, α) with m = 9.5, α = π/4 and δ p (m, α) = 1. The patterns shown are consistent with the capping phenomena; that is, when convection is fast (v = v 1) and diffusion is normal (D = D 0) a graduated distribution in the direction of flow streamlines is observed, a) and b). But in the presence of a fast convective transport, suitable values of the diffusion coefficient (e. g. D = D ext ) can reverse the capping effect and induce an effective randomization of the distribution of unbound LDL receptors, c) and d). In any event, in the presence of a retrograde membrane flow, no receptor surface clusters are formed.

Our results are mainly consistent with the capping phenomenon, and suggest that this is indeed concomitant with the idea of a retrograde membrane flow. Capping is described as the rearward transport across the leading lamella of various materials used to mark the cell surface [46, 14]. Bretscher [14] asserts that during capping, those receptors that are excluded by coated pits, and therefore do not recycle, may be continuously swept by a retrograde membrane flow. This author also states that whether such a non-circulating protein does get swept backwards or not depends on how fast it can diffuse by Brownian motion and how fast the flow is; and moreover that for individual proteins to overcome the sweeping effect of the flow, they must have diffusion coefficients (D) of about D ext [14]. This assertion was corroborated by our results, since our model predicts that for D values of the order of D ext the capping effect can be reversed. Nevertheless, a value of D ext lies far beyond the maximum experimentally feasible value for D 0. And even though Bretscher [14] claims that FBR often underestimates the true value of the diffusion coefficients of pertinent surface receptors, a value D ext   for D in the LDL experimental system seems unfeasible. Our results suggest that this is where the relative importance of non-uniform insertion becomes evident as a factor for boosting LDL receptor trapping. Indeed, even when convection is fast and receptors diffuse normally, our findings indicate that special forms of polarized receptor insertion can be arranged so as to produce a substantial reduction in τ du . This was more obvious in the case of a q-polarized mode, which could be raised as a feasible preferential-insertion mechanism for enhancing the LDL receptor trapping rate.

The problem of the surface distribution of unbound receptors is an important question pertaining to the LDL experimental system. For example, Wofsy et al. [40] argued that in the Robenek and Hesz [39] experiments, the LDL-gold particles were highly multivalent and thus may have been bound more efficiently to aggregate than single receptors, suggesting that the aggregation of newly inserted LDL receptors in regions around coated pits was still unproven. But beyond the caveats of Wofsy et al. [40] the results of Gross and Webb [47] were considered to provide a quantitative basis in support of a surface distribution of both LDL particles and their receptors in different-sized clusters. And again, the existence of these clusters was questioned by Sanan et al. [48] who reported experimental results that detected a dispersed or scattered population of LDL receptors, in addition to and clearly distinct from, clusters that formed soon after their insertion during recycling. In view of these disparate reports, Robenek et al. [49] performed new experiments on fibroblast and hepatocyte plasma membranes which confirmed their original claim that the initial insertion and display of LDL receptors in fibroblasts occurs predominantly as plaques. But experiments using fluorescently labeled LDL and influenza virus particles bound to the surface of human fibroblasts imaged with a cooled-scan CCD camera attached to a fluorescence microscope were found to yield a reasonably accurate measure of the proportion of single particles, but large errors were encountered in the proportion of large cluster sizes [50]. The high proportion of single particles found in these data provided evidence against the clusters claimed by Robenek et al. [49] and support the random insertion model. It is still a matter of disagreement among different research groups whether newly inserted LDL receptors are dispersed in the plasma membrane of cultured fibroblasts prior to entrapment in coated pits or whether they remain after reinsertion as groups of clusters. Our findings reveal that in the presence of a fast retrograde membrane flow, the expected surface display of unbound receptors must be consistent with the capping phenomenon; thus, the formation of steady-state surface clumps of LDL receptors is unachievable in the current setting (Figure  12).

Traditionally the existence of a retrograde flow has been a source of controversy among researchers, some of whom disclaim this paradigm (e.g. [5153]) while others convincingly address it (e.g. [14, 54, 55]). In particular, our results provide a theoretical support for the claim of Bretscher [14] which suggests that, if a retrograde flow does exist, the diffusion-convection ratio could be a fundamental determinant of the surface display of LDL receptors that are not bound to coated pits. And even though biochemical or biomechanical factors have been invoked in order to put forward other paradigms for the capping phenomenon (e.g. [49, 5658]), our findings support the views of Bretscher [59] and Ishihara et al. [60] on the adequacy of the retrograde-flow model. Hence, the involvement of diffusion, polarized reinsertion and a retrograde membrane flow in the formalization of the receptor-mediated endocytic cycle may provide a trustworthy paradigm for the interpretation of observations. Moreover, the fact that some characterizations of the retrograde flow-diffusion ratio could induce slower receptor trapping rates cannot provide a rebuttal for the retrograde-flow model. Indeed, Echavarria-Heras et al. [41] suggested that receptors could be projected to the plasma membrane in a preferential mode relative to the flow streamlines. We have studied here the consequences of such an arrangement and our results show that certainly, even in the presence of a fast retrograde flow, we could consider a q-polarized form as a mechanism for the enhancement of LDL receptor trapping rates.

Conclusions

Our results seem to indicate that even in the presence of a fast retrograde membrane flow, which could influence capping-like effects, diffusing LDL receptors that are sorted through polarized insertion in the form conceptualized here, could spur faster aggregation rates in coated pits than those associated with uniform insertion over the entire cell membrane. Moreover a q-polarized form could provide a paradigm for a highly efficient form of preferential insertion. The effectiveness of this mode is mainly explained by its placement relative to the flow streamlines. Furthermore, if the cross sections of the linked insertion region Ω q (m, α) and those associated with a coated pit are comparable, then even though a relative dominance of convection over diffusion might occur, the privileged position of the insertion region Ω q (m, α) would induce a fast receptor-trapping rate. But our results also hint that in spite of the efficiency of this device, the influence of a retrograde flow would cancel out the display of the corresponding surface clusters. Moreover, if new experiments do reveal that the surface clusters exist, then it could be necessary to envision a mechanism that explains their formation in the presence of a retrograde membrane flow. Meanwhile, for the case v = v 1 and D = D 0, a q-polarized insertion mode with m = 1.20 and α = π/16 provides a convenient paradigm for explaining the coexistence of a retrograde membrane flow and receptor aggregations that might be experimentally misinterpreted as surface clusters. Indeed, considering the experimental error associated with the measurement of the coated pit radius, a region of insertion of the form Ω q   with m = 1.20 could actually be a region placed over the boundary of a coated pit. And the surface clusters observed might be formed when recycling receptors are returned near the periphery of newly-formed coated pits so as to enhance the removal of the LDL ligand, but the existence of that insertion paradigm would require experimental substantiation.

Appendix A

In accordance with Eq. (17), the function S (c, p, q, m, α) inserts S * m , α particles per unit area per unit time at (r, θ)  in Ω *. Suppose that the number of particles that S * m , α inserts in Ω * is I *(m, α). Then, we have

a b - π + α π - α S c m , α rdrdθ = I c m , α ,
(A1)
ma b π - α π + α S p m , α rdrdθ = I p m , α ,
(A2)
a ma π - α π + α S q m , α rdrdθ = I q m , α .
(A3)

Let's denote by means of I Ω the total number of particles inserted by S (c, p, q, m, α) in Ω. Assume also, that a proportion δ p (m, α) of the total number of particles inserted in Ω is sorted by S (c, p, q, m, α) in Ω p (m, α), that a proportion δ q (m, α) of I Ω is placed into the region Ω q (m, α), and that a proportion δ c (m, α) of I Ω is delivered in Ω q (m, α). Then, we must have

I c m , α = δ c m , α I Ω ,
(A4)
I p m , α = δ p m , α I Ω ,
(A5)
I q m , α = δ q m , α I Ω ,
(A6)

and the balance equation

δ p m , α + δ q m , α + δ c m , α = 1 .
(A7)

From equations (A1) through (A3) and the mean value theorem there are numbers

r c , θ c Ω c m , α ,
(A8)
r p , θ p Ω p m , α ,
(A9)
r q , θ q Ω q m , α ,
(A10)

such that

δ p m , α = α S p r p θ p m , α , b 2 - a 2 m 2 I Ω ,
(A11)
δ q m , α = α S q r q θ q m , α , a 2 m 2 - a 2 I Ω ,
(A12)
δ c m , α = α S c r c θ c m , α , π - α b 2 - a 2 I Ω .
(A13)

Therefore, we have to introduce the consistency conditions

δ p b / a , α = 0 ,
(A14)
δ p m , 0 = 0 ,
(A15)
δ p 1 , π = 1 ,
(A16)
δ q 1 , α = 0 ,
(A17)
δ q b / a , π = 1 ,
(A18)
δ q m , 0 = 0 ,
(A19)
δ c m , π = 0 ,
(A20)
δ c m , 0 = 1 .
(A21)

A uniform insertion mode is obtained by setting S c b / a , 0 = c , where c is a positive constant. This is symbolically represented by means of S (c, 0, 0, b/a, 0) and given by

S c , 0 , 0 , b / a , 0 = c r , θ Ω c b / a , 0 0 r , θ Ω p b / a , 0 0 r , θ Ω q b / a , 0 ,
(A22)

where

c = I Ω π b 2 - a 2 .
(A23)

Also, uniform insertion is associated in the case S c m , π = 0 , S p m , π = S q m , π = c . This mode is known as pq-uniform insertion mode and denoted by means of S (0, p, q, m, π),

S 0 , c , c , m , π = 0 r , θ Ω c m , π c r , θ Ω p m , π c r , θ Ω q m , π .
(A24)

Again, from equations (A11) through (A13),

c = I Ω α b 2 - a 2 .
(A25)

Moreover, the insertion rate function

S c , c , c , m , α = c r , θ Ω c m , α c r , θ Ω p m , α c r , θ Ω q m , α ,
(A26)

which is obtained for the case S c m , α = S p m , α = S q m , α = c , also inserts particles uniformly all over Ω. This is known as a cpq-uniform insertion mode. From equations (A11) through (A13) we have

c = I Ω α b 2 - a 2 m 2 + π - α b 2 - a 2 + α a 2 m 2 - a 2 .
(A27)

The form of S (c, p, q, m, α) obtained by setting S c m , α = c , S p m , α = p and S q m , α = q , where c, p and q are positive constants will be known as a cpq-locally uniform insertion mode. That is, from equation (17) we will have

S c , p , q , m , α = c r , θ Ω c m , α p r , θ Ω p m , α q r , θ Ω q m , α
(A28)

with

αp b 2 - a 2 m 2 I Ω + αq a 2 m 2 - a 2 I Ω + αc π - α b 2 - a 2 I Ω = 1 .
(A29)

A second form of locally uniform insertion denoted by means of S (0, p, q, m, π) and known as a pq-locally uniform insertion mode is obtained for the case S c m , π = 0 , S p m , π = p and S q m , π = q , that is,

S 0 , p , q , m , π = 0 r , θ Ω c m , π p r , θ Ω p m , π q r , θ Ω q m , π
(A30)

where, from equations (A11) through (A13) p and q satisfy

αp b 2 - a 2 m 2 I Ω + αq a 2 m 2 - a 2 I Ω = 1 .
(A31)

The case in which receptors are inserted into two disjoint regions of abruptly contrasting rates is defined as a polarized insertion mode. It will be typified here as a mode that sorts the total number of recycled receptors over either the region Ω p (m, α), or Ω q (m, α), or Ω p (m, α)  Ω q (m, α). Moreover, a polarized insertion mode can be obtained in five different forms, two of these forms being radially symmetric.

The insertion rate function known as a q-plaque form is denoted by means of S (0, 0, q, m, π) and obtained for S c m , π = 0 , S p m , π = 0 , and S q m , π = q , being q a positive constant, that is,

S 0 , 0 , q , m , π = 0 r , θ Ω c m , α 0 r , θ Ω p m , α q r , θ Ω q m , α
(A32)

is radially symmetric and polarized, and from equations (A11) through (A13) we have

q = I Ω π a 2 m 2 - a 2 .
(A33)

The function S (0, 0, q, m, π) moreover, coincides with the plaque form insertion mode proposed by Wofsy et al. [40]. And, by setting S (0, 0, q, m, π) = Sr - β one obtains from equation (17) the decreasing and radially symmetric insertion-rate function considered by Echavarria-Heras et al. [41].

Further, what we call here a p-peripheral insertion mode was proposed by Bretscher [38]. It can be also considered as a radially symmetric polarized insertion mode. It is denoted through S (0, p, 0, m, π) and is obtained by setting S c m , π = 0 , S p m , π = p and S q m , π = 0 , with p a positive constant, that is

S 0 , p , 0 , m , π = 0 r , θ Ω c m , π p r , θ Ω p m , π 0 r , θ Ω q m , π ,
(A34)

and from equations (A11) through (A13) we have

p = I Ω π b 2 - a 2 m 2 .
(A35)

A first form of a non-radially symmetric polarized insertion mode is known here as p-polarized. It is denoted by means of S (0, p, 0, m, α)  and arises whenever S c m , α = 0 , S p m , α = p and S q m , α = 0 , where p is a positive constant, that is,

S 0 , p , 0 , m , α = 0 r , θ Ω c m , α p r , θ Ω p m , α 0 r , θ Ω q m , α
(A36)

where from equations (A11) through (A13) we have

p = I Ω α b 2 - a 2 m 2 .
(A37)

A second form of a non-radially symmetric polarized receptor insertion mode is known here as pq-polarized. It is denoted by means of S (0, p, q, m, α), and is linked to the case S c m , α = 0 , S p m , α = p , and S q m , α = q , with p and q positive constants, that is,

S 0 , p , q , m , α = 0 r , θ Ω c m , α p r , θ Ω p m , α q r , θ Ω q m , α
(A38)

and from equations (A11) through (A13) we have

I Ω p m , α α b 2 - a 2 m 2 + I Ω q m , α α a 2 m 2 - a 2 = 1 .
(A39)

A third non-radially symmetric and polarized form called q-polarized is symbolized by means of S (0, 0, q, m, α). It is linked to the case S c m , α = 0 , S p m , α = 0 and S q m , α = q , where q is a positive constant, that is,

S 0 , 0 , q , m , α = 0 r , θ Ω c m , α 0 r , θ Ω p m , α q r , θ Ω q m , α
(A40)

where from equations (A11) through (A13) we have

q = I Ω α a 2 m 2 - a 2 .
(A41)

Appendix B

In this section we outline the adaptation of the model for receptor mediated endocytosis presented in Echavarria-Heras et al. [41] for the case of a partitioned receptor insertion mode. For detailed mathematical proofs the reader is referred to Echavarria-Heras and Leal-Ramirez [61]. In the present settings the reference region Ω = x , y | a x 2 + y 2 b coincides with the Berg-Purcell [23] geometry. Equation (5) scales the outer radius b. A retrograde membrane flow will be represented by means of the vector field V x , y , t = v i with i being the unitary vector pointing in the direction of the positive x-axis and v, the strength of the flow in units of cm/s. Since we have assumed here a convective diffusion transport and also a partitioned receptor-insertion mode (cf. Eq. (17)) we must extend the model of equations (20) through (22) in order to include these mechanisms. The symbol D will stand for the 2-dimensional diffusion coefficient, which has units of cm 2/s and will be assumed to be constant. The ratio

λ = v 2 D ,
(B1)

having units cm - 1 is known as the fundamental ratio of the transport process because it combines the parameters v and D, which represent the factors that control receptor movement on the cell surface. Also, S(m, α, x, y)  will denote the rate of particle insertion, with units: particles/cm 2 - s. At steady state, the corresponding concentration of particles C(λ, x, y) satisfies the partial differential equation in the annulus Ω [62],

- div D C + div C V = S m , α , x , y
(B2)

Equivalently,

D 2 C - v C x + S m , α , x , y = 0 .
(B3)

Due to the geometry determined by Ω, it is convenient to switch to polar coordinates. Then, if we denote by means of n the normal unit vector pointing out radially towards (r, θ) for 0 ≤ θ ≤ 2π, and if J stands for the flux vector, the projection of the flux vector over n is given by

J n = D C x + vC cos θ .
(B4)

We will consider the partitioned insertion mode of equation (17), that is

S m , α , r , θ = S c m , α r , θ Ω c m , α S p m , α r , θ Ω p m , α S q m , α r , θ Ω q m , α .
(B5)

Because of the nature of the convective vector field C λmα (r, θ) is not expected to be radially symmetric if the strength of the membrane flow is different from zero. However, C λmα (r, θ)  will be symmetric about the x-axis. This property is formally expressed by

C λmα r , θ = C λmα r , - θ .
(B6)

The boundary at r = a remains as absorbing, i.e.

C λmα r , θ | r = a = 0 .
(B7)

On the other hand, the reflecting boundary condition given by equation (8) must be modified in order to account for the influence of the retrograde flow. We will consider that this boundary condition can be replaced by a flux-vanishing boundary condition at r = b. Hence the boundary condition at r = b will be formally represented by means of the equation

Ω b J n d Ω = 0 .
(B8)

Where ∂Ωb denotes the outer boundary of the annulus Ω, since we have assumed a steady-state for the receptors. This requirement for a vanishing flux of particles across the boundary at r = b can be considered a model for a dynamics in which those receptors crossing the boundary at r = b are transported into a contiguous influence area determined by the direction of the flow streamlines.

The further formal treatment can be greatly simplified by introducing a function G λmα (r, θ) defined through,

G λmα r , θ = C λmα r , θ e - λr cos θ .
(B9)

The substitution of Eq. (B9) into Eq. (B3) and the fulfilment of the symmetry and boundary conditions for C λmα (r, θ) show that, G λmα (r, θ) must satisfy

2 G λmα r , θ - λ 2 G λmα r , θ = - S c , p , q , m , α D e - λb cos θ ,
(B10)
G λmα a , θ = 0 ,
(B11)
G λmα r , θ = G λmα r , - θ ,
(B12)
G λmα r , θ periodic in θ ,
(B13)
- π π G λmα r , θ r | r = b - λ cos θ G λmα b , θ e - λb cos θ = 0 .
(B14)

Finally, for consistency a limiting radial symmetry condition must be added when the strength of the flow approaches zero. That is, we will require that for fixed D > 0, G λmα (r, θ) satisfies

lim λ 0 G λmα r , θ = C ds r
(B15)

Where C ds (r) is the solution of equation (6) subject to boundary conditions (7) and (8).

Using the boundary conditions (B11) trough (B14) it is possible to obtain, using separation of variables, a solution G λmα (r, θ) for equation (B10). Furthermore G λmα (r, θ) can be written in the form,

G λmα r , θ = k = 0 g λmα k r cos ,
(B16)

where for k = 0, 1 …, a ≤ m ≤ b/a, 0 ≤ α ≤ π and λR, the function g λmα k r satisfies the non-homogeneous problem

L λ k u + λ 2 u = σ k D - π π S c , p , q , m , α e - λr cos θ cos kθdθ ,
(B17)

with

σ k = 1 2 π k = 0 1 π k = 1 , 2 , , ,
(B18)

and

L λ k u = 1 r - u ' ' + k 2 r u
(B19)

defined for uC 2([a, b]) and satisfying boundary conditions

u a = 0 ,
(B20)
u ' b + γ k λ u b = 0 .
(B21)

Where γ k (λ) is the function given by

γ k λ = - λ I k λb / I k λb λ 0 - k / b λ = 0 ,
(B22)

with I k (∙) being the modified Bessel function of the first class and of order k for k = 0, 1 …, and λR. Notice also that the operator L λ k u given by equation (B19), as well as, the boundary conditions (B20) and (B21) originate when eigenfunctions ϕ λ (r, θ) for the operator A λ  = - 2 are supposed to have the form ϕ λ (r, θ) = R λ (r)ψ(θ). The function γ k (λ) given by equation (B22) is motivated by the requirement of equation (B15) of radial symmetry in the limiting case λ → 0 for D > 0 fixed.

For the operator L λ k u a Green’s function H λ k r , t exists. It is given by

H λ k r , t = I k λa K k λr - K k λa I k λr I k λt I k λa r t I k λa K k λt - K k λa I k λt I k λr I k λa t r ,
(B23)

where K k (∙) stands for the modified Bessel function of the second kind and of order r. Using formulas (9.6.10) and (9.6.53) in Abranowitz and Stegun [63] we can obtain an expression for H λ k r , t depending only on I k (z). This permits the proof of the existence of the solution G λmα (r, θ) in form given by equation (B16) and which depends continuously on the boundary data, namely

C λmα r , θ = k = 0 g λmα k r e λr cos θ cos
(B24)

with g λmα k r expressed in terms of the Green’s function H λ k r , t through the expression

g λmα k r = σ k D a b - t H λ k r , t - π π S c , p , q , m , α e - λt cos θ cos kθdθ dt ,
(B25)

and σ k given by equation (B18) [61].

The symbol k λmα + denotes the pertinent forward rate constant, which can be obtained by calculating the number of receptors hitting the trap in unit time and dividing it by the average concentration (cf. Eq. (1)). We have nevertheless assumed the existence of a steady-state concentration C λmα (r, θ) in the reference annulus Ω ¯ for all λ. So the number of receptors hitting the trap in unit time coincides with the number of particles that are inserted into the reference annulus in unit time. Hence, we have equivalently

k λmα + = a b - π π S c , p , q , m , α rdrdθ 1 π b 2 a b - π π C λmα r , θ rdrdθ ,
(B26)

and by virtue of equations (1) and (3) we must have for τ λmα the associated mean capture time,

τ λmα = k = 0 a b - π π g λmα k r e λr cos θ cos kθrdrdθ a b - π π S c , p , q , m , α rdrdθ
(B27)

For the case of a uniform receptor insertion (i.e. the case m = b/a and α = 0 in equation (B5)), we will have S (c, p, q, m, α) = S, a constant for a ≤ r ≤ b. Denoting by g λs k r the corresponding form of g λmα k r and then using the result,

- π π e λt cos θ cos kθdθ = 2 π - 1 k I k λt
(B28)

(cf. Abramowitz and Stegun [63]) equation (B25) yields

g λS k = 2 πS σ k - 1 k + 1 D a b t H λ k r , t I k λt dt .
(B29)

Meanwhile, the corresponding mean capture time τ λs becomes

τ λs = k = 0 ϑ k - 1 k + 1 a b r a b t H λ k r , t I k λt dt I k λr dr D b 2 - a 2 ,
(B30)

where

ϑ k = 2 if k = 0 4 if k 1 .

It can be also shown [61] that for D > 0 and fixed, we have

lim λ τ λs = b 4 ln b a 2 D b 2 - a 2 - 3 b 2 - a 2 8 D ;
(B31)

that is, if receptor insertion is uniform all over Ω, then whenever λ approaches zero, τ λs approaches the value τ du of equation (12), calculated by Berg and Purcell [23].

References

  1. Basu Sandip K: Receptor-mediated endocytosis: an overview of a dynamic process. J Biosci. 1984, 6 (4): 535-542. 10.1007/BF02703904.

    Article  Google Scholar 

  2. Roth TF, Porter KR: Yolk protein uptake in the oocyte of the mosquito A edes aegypti.L. J Cell Biol. 1964, 20 (3): 13-32.

    Article  Google Scholar 

  3. Fawcett DW: Surface specializations of absorbing cells. J Histochem Cytochem. 1965, 13: 75-91. 10.1177/13.2.75.

    Article  CAS  PubMed  Google Scholar 

  4. Friend DS, Farquhar MG: Functions of coated vesicles during protein absorption in the rat vas deferens. J Cell Biol. 1967, 35: 357-76. 10.1083/jcb.35.2.357.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Pearse BMF: Chathrin: a unique protein associated with intracellular transfer of membrane by coated vesicles. Proc Natl Acad Sci U S A. 1980, 73: 1255-1259.

    Article  Google Scholar 

  6. Stahl P, Schwarts AL: Receptor mediated endocytosis. J Clin Investig. 1986, 7: 657-662.

    Article  Google Scholar 

  7. Anderson RGW, Goldstein JL, Brown MS: A mutation that impairs the ability of lipoprotein receptors to localize in coated pits on the cell surface of human fibroblasts. Nature. 1977, 270: 695-699. 10.1038/270695a0.

    Article  CAS  PubMed  Google Scholar 

  8. Kerkhof PV, Govers R, Santos CM A d, Strous GJ: Endocytosis and degradation of the growth hormone receptor are proteasome-dependent. The Journal of Biological Chemestry. 2002, 275 (3): 1575-1580.

    Article  Google Scholar 

  9. Molfetta R, Belleudi F, Peruzzi G, Morrone S, Leone L, Dikic I, Piccoli M, Frati L, Torrisi MR, Santoni A, Paolini R: CIN85 Regulates the ligand-dependent endocytosis of the IgE receptor: a new molecular mechanism to dampen mast cell function. J Immunol. 2005, 175: 4208-4216. 10.4049/jimmunol.175.7.4208.

    Article  CAS  PubMed  Google Scholar 

  10. Yang J, Chen H, Vlahov IR, Chen JH, Low PS: Evaluation of disulfide reduction during receptor-mediated endocytosis by usin FRET imaging. Proc Natl Acad Sci. 2006, 103: 3872-19877.

    Google Scholar 

  11. Beckebaum S, Cicinnati VR, Zhang X, Ferencik S, Frilling A, Grosse-Wilde H, Broelsch CE, Gerken G: Hepatits B virus-induced defect of monocyte-derived dendritic cells leads to impaired T helper type 1 response in vitro: mechanisms for viral immune escape. Immunology. 2003, 109: 487-495. 10.1046/j.1365-2567.2003.01699.x.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Brown MS, Goldstein JL: How LDL receptors influence cholesterol and atherosclerosis. Sci Am. 1984, 251: 52-60.

    Article  Google Scholar 

  13. Goldstein JL, Brown MS: The low-density lipoprotein pathway and its relation to atherosclerosis. Annu Rev Biochem. 1977, 46: 897-930. 10.1146/annurev.bi.46.070177.004341.

    Article  CAS  PubMed  Google Scholar 

  14. Bretscher MS: Endocytosis: relation to capping and cell locomotion. Science. 1984, 224: 681-685. 10.1126/science.6719108.

    Article  CAS  PubMed  Google Scholar 

  15. Anderson RGW, Goldstein JL, Brown MS: Localization of low density lipoprotein receptors on plasma membranes of normal human fibroblasts and their absence in cells from a familial hypercholesterolemia homozygote. Proc Natl Acad Sci U S A. 1976, 73: 2434-2438. 10.1073/pnas.73.7.2434.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Orci L, Carpentier JL, Perrelet A, Anderseon RGW, Goldstein JL, Brown MS: Occurrence of low density lipoproteins receptor within large pits on the surface of human fibroblasts as demonstrated by freeze-etching. Exp Cell Res. 1978, 113: 1-13. 10.1016/0014-4827(78)90081-2.

    Article  CAS  PubMed  Google Scholar 

  17. Anderson RGW, Vasile EMello RJ, Brown MS, Goldstein JL: Immunocytochemical visualization of coated pits and vesicles in Human Fibroblasts: relation to low density lipoprotein receptor distribution. Cell. 1978, 15: 919-993. 10.1016/0092-8674(78)90276-3.

    Article  CAS  PubMed  Google Scholar 

  18. Goldstein JL, Brown MS, Anderson RGW, Russell DW, Schneider WJ: Receptor-mediated endocytosis: concepts emerging from the LDL receptor system. Annual Review Cell Biology. 1985, 1: 1-39. 10.1146/annurev.cb.01.110185.000245.

    Article  CAS  Google Scholar 

  19. Guadorov I, Santini F, Warren RA, Keen JH: Spatial control of coated-pit dynamics in living cells. Natural Cell Biology. 1999, 1: 1-7. 10.1038/8971.

    Article  Google Scholar 

  20. Schlessinger AL, Shechter Y, Willigham MC, Pastan I: Direct visualization of the binding aggregation and internalization of insulin and epidermal growth factor on fibroblastic cell. Proc Natl Acad Sci U S A. 1978, 75: 2659-2663. 10.1073/pnas.75.6.2659.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Goldstein B, Wofsy C, Bell G: Interaction of low density lipoprotein receptors with coated pits on human fibroblasts: estimate of the forward rate constant and comparison with the diffusion limit. Proc Natl Acad Sci U S A. 1981, 78 (9): 5695-5698. 10.1073/pnas.78.9.5695.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Smoluchowski MV: Drei Vorträge über Diffusion, Brownsche Molekularbewegung und Koagulation von Kolloidteilchen. Physics. 1916, 17: 585-599.

    Google Scholar 

  23. Berg HC, Purcell EM: Physics of chemoreception. Biophys J. 1977, 20: 139-219.

    Google Scholar 

  24. Goldstein B, Wofsy C, Echavarria-Heras H: Effect of membrane flow on the capture of receptors by coated pits. Biophys J. 1988, 53: 405-414. 10.1016/S0006-3495(88)83117-5.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Goldstein B, Griego R, Wofsy C: Diffusion-limited forward rate constants in two dimensions. Applications to the trapping of cell surface receptors by coated pits. Biophys J. 1984, 46: 573-586. 10.1016/S0006-3495(84)84056-4.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Abercrombie M, Heaysman JEM, Pegrum SM: The locomotion of fibroblasts in culture. III. Movement of particles on the dorsal surfaces of the leading lamella. Exp Cell Res. 1970, 62: 389-398. 10.1016/0014-4827(70)90570-7.

    Article  CAS  PubMed  Google Scholar 

  27. Ingram YM: A side view of moving fibroblasts. Nature (London). 1969, 222: 641-644. 10.1038/222641a0.

    Article  CAS  Google Scholar 

  28. Harris AK, Dunn G: Centripetal transport of attached particles on both surfaces of moving fibroblasts. Exp Cell Res. 1972, 72: 519-523. 10.1016/0014-4827(72)90022-5.

    Article  Google Scholar 

  29. Bretscher MS: Directed lipid flow in cell membranes. Nature. 1976, 20: 21-22.

    Article  Google Scholar 

  30. Dembo M, Harris AK: Motion of particles adhering to the leading lamella of crawling cells. J Cell Biol. 1981, 91: 528-536. 10.1083/jcb.91.2.528.

    Article  CAS  PubMed  Google Scholar 

  31. Thompson CRL, Bretscher MS: Cell polarity and locomotion as well as endocytosis, depend on NSF. Development. 2002, 129: 4185-4192.

    CAS  PubMed  Google Scholar 

  32. Alberts B, Bray D, Lewis J, Raff M, Roberts K, Watson JD: Molecular biology of the cell. 1983, New York and London: Garland Publishing Inc.

    Google Scholar 

  33. Brown D, Breton S: Sorting proteins to their target membranes. Kidney Int. 2000, 57: 816-824. 10.1046/j.1523-1755.2000.00920.x.

    Article  CAS  PubMed  Google Scholar 

  34. Chu CL, Buczek-Thomas JA, Nugent MA: Heparan sulfate proteoglycans modulate fibroblast growth factor-2 binding through a lipid-raft-mediated mechanism. Biochem J. 2004, 379: 331-341. 10.1042/BJ20031082.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Marco E, Wedlich-Soldner R, Li R, Altschuler SJ, Wu LF: Endocytosis optimizes the dynamic localization of membrane proteins that regulate cortical polarity. Cell. 2007, 129: 411-422. 10.1016/j.cell.2007.02.043.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Linderman JJ, Lauffenburger DA: Receptor/ligand Sorting along the Endocytic Pathway. 1989, Berlin: Springer

    Book  Google Scholar 

  37. Puertollano R: Clathrin-mediated transport: assembly required, workshop on molecular mechanisms of vesicle selectivity. EMBO Rep. 2004, 5 (10): 942-946. 10.1038/sj.embor.7400249.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Bretscher MS: Distribution of receptors for transferrin and lowdensity lipoprotein on the surface of giant HeLa cells. Proc Natl Acad Sci U S A. 1983, 80: 454-458. 10.1073/pnas.80.2.454.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Robenek H, Hesz A: Dynamics of low density lipoprotein receptors in the plasma membrane of cultured human skin fibroblasts as visualized by colloidal golds in conjunction with surface replicas. Eur J Cell Biol. 1983, 31: 275-282.

    CAS  PubMed  Google Scholar 

  40. Wofsy C, Echavarría-Heras H, Goldstein B: Effect of preferential insertion of LDL receptors near coated pits. Cell Biophys. 1985, 7: 197-204. 10.1007/BF02790465.

    Article  CAS  PubMed  Google Scholar 

  41. Echavarria-Heras H, Solana-Arellano E, Leal-Ramirez C: Surface aggregation patterns of LDL receptors near coated pits II. The retrograde membrane flow-diffusion and generalized plaque-form insertion mechanism. Mathematical Med Biol. 2010, 29 (2): 109-30.

    Article  Google Scholar 

  42. Dowell JW: Graphical Ray Tracing: An Advanced Text. 1973, El Segundo CA: Chan StreeL

    Google Scholar 

  43. Hauser J: Three-dimentional visualization of coated vesicles formation in fibroblasts. J Cell Biol. 1980, 84: 560-583. 10.1083/jcb.84.3.560.

    Article  Google Scholar 

  44. Solana-Arellano E, Echavarria-Heras E, Leal-Ramirez C: Surface aggregation patterns of LDL receptors near coated pits I. The radially convective diffusion and generalized insertion mechanism. IMA J Math Appl Med Biol. 1998, 15: 351-366. 10.1093/imammb/15.4.351.

    Article  CAS  PubMed  Google Scholar 

  45. Barak LS, Webb WW: Diffusion of low density lipoprotein-receptor complex on human fibroblasts. J Cell Biol. 1982, 95: 846-852. 10.1083/jcb.95.3.846.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Holifield BF, Ishihara A, Jacobson K: Comparative behavior of membrane protein-antibody complexes on montile fibroblasts: implications for a mechanism of capping. J Cell Biol. 1990, 111 (6 part 1): 2499-2512.

    Article  CAS  PubMed  Google Scholar 

  47. Gross D, Webb WW: Molecular counting of low-density lipoprotein particles as individuals and small clusters on cell surfaces. Biophys J. 1986, 9: 901-911.

    Article  Google Scholar 

  48. Sanan DA, Van der Westhuyzen DR, Gevers W, Coetzee GA: The surface distribution of low density lipoprotein receptors on cultured fibroblasts and endothelial cells. Ultrastructural evidence for dispersed receptors. Histochem Cell Biol. 1987, 86 (5): 517-523.

    CAS  Google Scholar 

  49. Robenek H, Harrach B, Severs NJ: Display of low density lipoprotein receptors is clustered, not dispersed, in fibroblast and hepatocyte plasma membranes. Arterioscler Thromb Vasc Biol. 1991, 11: 261-271. 10.1161/01.ATV.11.2.261.

    Article  CAS  Google Scholar 

  50. Morrison IEG, Anderson CM, Georgiou GN, Stevenson GVW, Cherry RJ: Analysis of receptor clustering on cell surfaces by imaging fluorescent particles. Biophys J. 1994, 67: 1280-1290. 10.1016/S0006-3495(94)80600-9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Kucik DF, Elson EL, Sheetz MP: Cell migration does not produce membrane flow. J Cell Biol. 1990, 111: 1617-1622. 10.1083/jcb.111.4.1617.

    Article  CAS  PubMed  Google Scholar 

  52. Lee J, Gustafsson M, Magnusson KE, Jacobson K: The direction of membrane lipid flow in locomoting polymorphonuclear leukocytes. Science. 1990, 247: 1229-1233. 10.1126/science.2315695.

    Article  CAS  PubMed  Google Scholar 

  53. Traynor D, Kay RR: Possible roles of the endocytic cycle in cell motility. J Cell Sci. 2007, 120: 2318-2327. 10.1242/jcs.007732.

    Article  CAS  PubMed  Google Scholar 

  54. Sornette D: A physical model of cell crawling motion. J de Physique France. 1989, 50: 1759-1770. 10.1051/jphys:0198900500130175900.

    Article  Google Scholar 

  55. Vallonton P, Danuser G, Bohnet S, Meister JJ, Verkhovsky AB: Tracking retrograde flow in keratocytes: new from the front. Mol Biol Cell. 2005, 16: 1223-1231. 10.1091/mbc.E04-07-0615.

    Article  Google Scholar 

  56. De Petris S: Distribution and Mobility of Plasma Membrane Components on Lymphocytes. Dynamic Aspects of Cell Surface Organization, Cell Surface Reviews. Edited by: Poste G, Nicolson L. 1977, Amsterdam: Elsevier, 643-728. 3

    Chapter  Google Scholar 

  57. Michitson TJ, Cramer LP: Actin-based cell motility and cell locomotion. Cell. 1996, 84: 371-379. 10.1016/S0092-8674(00)81281-7.

    Article  Google Scholar 

  58. Coutsias EA, Wester MJ, Perelson AS: A nucleation theory of cell surface capping. J Stat Phys. 1996, 87: 1179-1203.

    Article  Google Scholar 

  59. Bretscher MS: Fibroblasts on the move. J Cell Biol. 1988, 106: 235-237. 10.1083/jcb.106.2.235.

    Article  CAS  PubMed  Google Scholar 

  60. Ishihara A, Holifield B, Jacobson K: Analysis of lateral redistribution of a monoclonal antibody complex plasma membrane glycoprotein which occurs during cell locomotion. J Cell Biol. 1988, 106 (2): 329-343. 10.1083/jcb.106.2.329.

    Article  CAS  PubMed  Google Scholar 

  61. Echavarria-Heras HA, Leal-Ramirez C: A theoretical assessment of the combined effect of a transversely convective diffusion transport and a plaque-form reinsertion mode in the LDL endocytic system. Aportaciones Matemáticas, Serie Comunicaciones. 2007, 37: 11-58.

    Google Scholar 

  62. Stakgold I: Green’s Functions and Boundary Value Problems. 1979, New York: John Wiley and Sons Inc.

    Google Scholar 

  63. Abramowitz M, Stegun IA: Handbook of Mathematical Functions. 1965, Washington DC: National Bureau of Standards

    Google Scholar 

Download references

Acknowledgments

We are indebted to Dr. Ernesto Franco Vizcaino for language editing. We are also grateful to Jose Maria Dominguez and Francisco Ponce for the art work.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Héctor Echavarria-Heras.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors&amp;#x2019; contributions

HEH conceived, designed and performed the research. CLR performed the required mathematical proofs and numerical analysis procedures. OC supervised the research and revised the manuscript critically at a formal level. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Echavarria-Heras, H., Leal-Ramirez, C. & Castillo, O. Surface aggregation patterns of LDL receptors near coated pits III: potential effects of combined retrograde membrane flow-diffusion and a polarized-insertion mechanism. Theor Biol Med Model 11, 23 (2014). https://doi.org/10.1186/1742-4682-11-23

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1742-4682-11-23

Keywords