A minimal self-organisation model of the Golgi apparatus
Figures

Model of golgi structure and transport.
(A) Top: classical representation of the structure of the Golgi, showing some of the important protein actors. The three main membrane identities are shown in different colors (cis: blue, medial: orange, trans: green): ER = Endoplasmic Reticulum, ERGIC = ER Golgi intermediate Compartment, TGN = Trans Golgi-Network, PM = Plasma Membrane. Bottom: Sketches of the four main models of Golgi transport (see text). (B) Our quantitative model of Golgi self-organisation. The left boundary is the ER, composed of a cis-membrane identity, and the right boundary is the TGN, composed of a trans-membrane identity. Golgi compartments self-organise via three stochastic mechanisms: Fusion: (1) All compartments can aggregate using homotypic fusion mechanisms: the fusion rate is higher between compartment of similar identities. (2) Each compartment can exit the system by fusing homotypically with the boundaries. Budding: (3) Each compartment larger than a vesicle can create a vesicle by losing a patch of membrane. Biochemical conversion: (4) Each patch of membrane undergoes a conversion from a cis to a trans-identity. New cis-vesicles (0) bud from the ER at a constant rate. In the sketch, the boundaries also contain neutral (gray) membrane species that dilute their identity (impact of this dilution in Appendix 7).

Steady-state of the self-organized model of Golgi apparatus.
(A) Size distribution of compartments for a biochemical conversion rate and different values of the budding rate . The red bar shows the characteristic compartment size. (B) Phase diagram of the size of the compartments as the function of and . Black lines: theoretical prediction (Equations 4,6) (C) Phase diagram of the purity of the system as a function of and . Dashed line is , black lines are a theoretical prediction (see Appendix 3). (D) System snapshots for showing the mixed regime (square - ), the sorted regime (triangle - ), and the vesicular regime (circle - ) - see text. (E) Average compartment purity as a function of their average size, obtained by varying the budding rate for (black dots: simulation results – gray line: guide for the eyes). For all panels, . See also Appendix 5 for further characterizations of the steady-state organisation, Appendix 7 for the role of composition of the exit compartment, and Appendix 8 for results with alternative budding and fusion kinetics.

Directionality of vesicular transport, displayed varying the budding rate for a biochemical conversion rate (equal amount of all species in the system).
(A) System purity as a function of the budding rate – a cut through the purity phase diagram shown Figure 2C. (B) Normalized enrichment in cis (blue), medial (orange) and trans (green) identities between the acceptor and donor compartments during vesicular exchange (solid lines are guides for the eyes). See Appendix 6 for non-normalized fluxes. (C) Sketches showing the direction of the dominant vesicular fluxes – consistent with data of A and B and data from Appendix 6—figure 1 – omitting less contributive transports. Before the purity transition (low values of - purity ~ 0.5) the vesicular flux is retrograde. After the transition (high values of - purity ~ 1) the vesicular flux is anterograde. Around the crossover (), the vesicular flux is centripetal and oriented toward medial-compartments. The centripetal flux disappears if inter-compartments fusion is prohibited, see Appendix 8.

Relationship between the structure of the system and the vesicular fluxes.
Distribution of the total size of the system (truncated scale, maximum value in brackets) and number of compartments, and the vesicular flux between them, shown in a triangular composition space. Each point represents a given composition of compartments and the triangle vertices are compartments of pure cis, medial and trans identity. Arrows show the vesicular flux. The base of each arrow is the composition of the donor compartment and the tip is the average composition of receiving compartments (ignoring back fusion). The opacity of the arrows is proportional to the flux of transported cargo going through this path per unit time (), normalized by the total number of cargo-proteins. The dashed line on the top-left triangle is the theoretical prediction given by Equation 9. See also Appendix 6 for further characterizations of the vesicular transport, and Appendix 8 to discuss the impact of other budding and fusion implementations on the vesicular fluxes.

Link between the total size of the Golgi and its steady-state organisation.
(A) Size dependence of the steady-state organisation of the Golgi (to be compared with Figure 2 in the main text). We look at two different systems that only differ in term of the average total size , that equals 30 or 1000. Typical size of compartments and average system purity are shown as a function of and . In both cases (see Appendix 1). (B) Average steady-state total size , as a function of and , for simulation in which the influx is set to have (using Equation 6).

Graphical representation of the size and purity of compartments.
Each compartment is represented as a sphere with an area equal to the area of a vesicle (the smallest elements in this snapshot) times the number of vesicles that fused together to create the compartment. The color of the compartment reflects its composition, according to the color code on the triangular composition phase space: 100% cis at the top, medial at the bottom right, trans at the bottom left. The purity of each of these compartments is its distance from a perfectly mixed compartment at the center of the triangle (middle arrow). This distance is normalized so that it equals 1 for a perfectly pure compartment (bottom arrow).

Analytical description of the cis-compartment purity.
Related to Figure 2 - main text. (A) Analytical model of the steady-state contamination by a medial-species (orange) of a cis-compartment (blue). The fraction occupied by the medial-species is named . The eight events (1–8) impacting are listed in the text. (B-C). Time variation as a function of , with and for different (0.9 dashed line, 1.05 solid line, and dotted line). has a stable value when and or (black arrows). It is potentially unstable when and (red arrows). (B) For (dashed) exhibits a stable fixed-point for but with as (solid) and (dotted) exhibit a stable concentration for . (C) For (dashed) exhibits a stable fixed-point for (solid) exhibits an unstable fixed-point near (dotted) does not display this unstable fixed-point as (see text).

Transient regime towards the steady-state.
System’s size (total number of membrane patches in the system) as a function of time (in unit ) for three different values of the budding rate (A: , B: , C: ). The top row shows the full simulations, with a vertical dashed bar when the number of simulation steps reaches 106. The steady-state data analysed and discussed in this paper are recorded after this time, and the full simulation typically lasted 107 time-steps. The bottom row shows the evolution at early time (de-novo formation). The red lines are the analytical solution given by Appendix 4 Equation 3.

Exit kinetics of passive cargo.
Simulated amount of cargo in the Golgi following a pulse-chase like simulation in which 1000 cargo molecules are trapped in the ER and allowed to be packaged into ER exiting vesicles at . This is expected to mimick the iFRAP experimental setup used in Patterson et al., 2008 and to the RUSH setup used in Boncompain et al., 2012. The decay of the amount of cargo is close to a single-time exponential decay, is agreement with experimental observations for transiting cargo (Figure 2i of Patterson et al., 2008 and Figure 3e in Boncompain et al., 2012). Parameters are , corresponding to the ‘sorted regime’ (the triangle of Figure 2) which we expect to be the most reasonable for mammalian Golgi.

Steady-state composition and fluctuation of compartments.
Related to Figure 2 - main text. (A–B) Characterization of the temporal fluctuations of the system as a function of the budding rate . (A) Average values of the system total size (top panel) and standard deviation of the size (bottom panel). (B) Average values of the system purity (top panel) and standard deviation of the purity (bottom panel). (C) Distribution of the total size, number of compartments and their average size, as a function of the composition of the compartments, shown in a triangle plot (see Appendix 2) for different values of the (normalized) budding rate . The average size for a given composition is calculated as the average size divided by the average number of compartments. Plots are shown for four different values of , with . The dashed black line on the top, left panel (low regime) corresponds to the theoretical prediction of Equation 9. (D) The evolution of the compartment’s size as a function of time, in the low limit, agrees with the linear prediction of Equation 7.

Detailed characterizations of the vesicular transport.
Related to Figures 3–4 - main text. for all panels. (A) Characterization of vesicle back-fusion (fusion of vesicles with the same compartment or with a compartment of comparable composition). Top panel: fraction of vesicular flux undergoing back-fusion, as a function of the composition of the budding compartments, for different values of the budding rate . Bottom panel: non-normalized enrichment in cis (blue), medial (orange) and trans-species (green) seen by a passive cargo, between a budding and the next fusion event, as a function of . (B) Vesicular fluxes of cis (blue), medial (orange) and trans (green) vesicles for different values of the budding rate . The base of each arrow gives the composition of the donor compartment and the tip gives the average composition of receiving compartments - ignoring back-fusion - as in Figure 4, main text. The opacity of the arrows represents the normalized number of cargo (divided by the total number of cargo in the system) transported per unit time (normalized by the fusion rate). (C) Net vesicular flux, excluding vesicle biochemical conversion, for different values of the budding rate . Top panel: vesicular flux is represented as a vector field as in Figure 4, main text. Bottom panel: vesicular flux is represented as the enrichment between the donor and the acceptor compartment, excluding back-fusion, as in Figure 3, main text.

Impact of the composition of the boundaries of the system (ER and TGN) on the steady-state organisation.
Related to Figure 2 - main text. (A–B). The fraction of active species driving homotypic fusion with the boundaries are identical for the ER and the TGN. A. (same value as in the main text) B. . For both cases, the top panel is the purity diagram as a function of the normalized biochemical conversion rate and budding rate . Dashed lines correspond to and solid lines to . For both cases, the bottom panel shows some examples of the size distribution for and different values of . (C) Purity diagram as a function of and when fusion with the ER is abolished: whereas . Solid lines are the predictions of the analytical model developed in Appendix 3.

Impact of different budding and fusion kernels on the system's organisation and vesicular transport.
Related to Figures 2–4 - main text. . (A–C) Linear budding scheme. (D–F) Inter-compartment fusion is prohibited. (A and D) Size distribution of compartment for equals 102 (top), 1 (middle) and 102 (bottom). (B and E) Steady-state purity as a function of . (C and F) Enrichments in cis (blue), medial (orange) and trans-species (green) during vesicular transport (details of computations of this enrichment presented in Appendix 2) as a function of .

Possible experimental strategies to test the predicted correlation between size and purity of Golgi cisternae.
The left column reproduces Figure 2B–C, showing the phase diagrams for the purity (top) and size (bottom) of Golgi cisternae upon varying the ratio of budding to fusion rates and the ratio of biochemical conversion to fusion rate . The right column shows a zoom of the dashed square area, showing the observed result of physiological perturbations of cisternae purity and size (solid arrows) and prediction for combined perturbations (dashed arrows). Red arrow: observed result of a decrease of the budding rate by impairing COP I (Papanikou et al., 2015). Dashed red arrow: prediction for the combined decrease of the biochemical conversion rate by Ypt1 down-regulation. Green arrow: observed result of an increase of the biochemical conversion rate by Ypt1 over-expression (Kim et al., 2016). Dashed green arrow: prediction for the combined increase of the budding rate by COP I up-regulation.
Tables
System’s parameters.
Budding rate per patches of membrane in a compartment | |
---|---|
biochemical conversion rate for each patch of membrane (cis→medial→trans) | |
Maximum fusion rate between two compartments | |
Injection rate of cis-vesicles from the ER | |
Budding rate normalized by the fusion rate () | |
biochemical conversion rate normalized by the fusion rate () | |
Injection rate normalized by the fusion rate () | |
Fraction of active species (driving homotypic fusion) in the ER (cis) or the TGN (trans) | |
When different between boundaries, in the boundary ( equals ER or TGN) |
Compartments’ description.
Fraction of a given species i (cis, medial or trans) in a compartment | |
---|---|
Composition of a compartment, whose components are the fractions | |
Size of a compartment, (i.e. number of membrane patches in the compartment) | |
Number of membrane patches of the identity, in a compartment | |
Number of different identities (from 1 to 3) in a single a compartment | |
Number of vesicles surrounding a compartment that share its composition | |
Purity of a compartment | |
Budding flux of vesicles decorated by the identity , from a compartment |
Steady-state description.
Size distribution of compartments in the system | |
---|---|
Cutoff in the size distribution, namely the typical size of big compartments | |
Total number of membrane patches, in all compartments, in the steady-state system | |
Total number of patches of the i species (cis, medial or trans) |
Additional files
-
Source code 1
Simulation source code.
- https://cdn.elifesciences.org/articles/47318/elife-47318-code1-v1.zip
-
Transparent reporting form
- https://cdn.elifesciences.org/articles/47318/elife-47318-transrepform-v1.pdf