Cited by 222Cellular transport and secretionLipid Membrane Structure and BehaviorFungal and yeast genetics research

The DOI matches its Crossref record; 59 of the 59 references it cites were found in a registry or an index.

The Mechanochemistry of Endocytosis

Open access

Abstract

Endocytic vesicle formation is a complex process that couples sequential protein recruitment and lipid modifications with dramatic shape transformations of the plasma membrane. Although individual molecular players have been studied intensively, how they all fit into a coherent picture of endocytosis remains unclear. That is, how the proper temporal and spatial coordination of endocytic events is achieved and what drives vesicle scission are not known. Drawing upon detailed knowledge from experiments in yeast, we develop the first integrated mechanochemical model that quantitatively recapitulates the temporal and spatial progression of endocytic events leading to vesicle scission. The central idea is that membrane curvature is coupled to the accompanying biochemical reactions. This coupling ensures that the process is robust and culminates in an interfacial force that pinches off the vesicle. Calculated phase diagrams reproduce endocytic mutant phenotypes observed in experiments and predict unique testable endocytic phenotypes in yeast and mammalian cells. The combination of experiments and theory in this work suggest a unified mechanism for endocytic vesicle formation across eukaryotes.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1371/journal.pbio.1000204

Cites

59 of 59 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1 Endocytic dynamics in budding yeast. (A) Timelines for endocytic protein recruitment as determined by multicolor fluorescence microscopy analysis. Sla1p, which is an endocytic adaptor protein, represents the endocytic coat. Abp1p is an actin-binding protein and faithfully reports on actin dynamics. Sjl2p is the yeast synaptojanin that hydrolyzes PIP2. PIP2 represents the lipid module and is believed to be the recruitment signal for many endocytic proteins. Rvs167p, yeast Amphiphysin, contains a BAR domain capable of sensing/binding curved membranes and deforming membranes. (Sla1 and Abp1 data are from [8], Sjl2 data are from [11], Rvs167 data are determined in this work from six individual patches in cells expressing Rvs167-GFP and aligned to the relative timing of Sjl2 appearance.) (B) Spatial profiles of endocytic membrane and the key endocytic proteins as revealed by EM [15].
Figure 1 Endocytic dynamics in budding yeast. (A) Timelines for endocytic protein recruitment as determined by multicolor fluorescence microscopy analysis. Sla1p, which is an endocytic adaptor protein, represents the endocytic coat. Abp1p is an actin-binding protein and faithfully reports on actin dynamics. Sjl2p is the yeast synaptojanin that hydrolyzes PIP2. PIP2 represents the lipid module and is believed to be the recruitment signal for many endocytic proteins. Rvs167p, yeast Amphiphysin, contains a BAR domain capable of sensing/binding curved membranes and deforming membranes. (Sla1 and Abp1 data are from [8], Sjl2 data are from [11], Rvs167 data are determined in this work from six individual patches in cells expressing Rvs167-GFP and aligned to the relative timing of Sjl2 appearance.) (B) Spatial profiles of endocytic membrane and the key endocytic proteins as revealed by EM [15].CC BY · PubMed Central
Figure 2 Mechanochemical feedback mechanism for endocytosis in budding yeast. The thin arrows represent activation effects, and the bar ends represent inhibition effects. The local spatial coordinate along the membrane surface is the arc length s with unit length 1 nm. The bud region is defined by the arc length 0≤s≤s 1, the lipid phase boundary is at s = s 1, and the tubule region starts from s = s 1+1, where s 1 is chosen to be 100. We assume that membrane shape is cylindrically symmetric. ϕ(s) is the membrane tangential angle and r(s) is the radius of the tubule. The mean curvature, Ω(s), is the average of the curvatures in the radial and tangential directions; it measures the overall extent of the PIP2 head group exposure.
Figure 2 Mechanochemical feedback mechanism for endocytosis in budding yeast. The thin arrows represent activation effects, and the bar ends represent inhibition effects. The local spatial coordinate along the membrane surface is the arc length s with unit length 1 nm. The bud region is defined by the arc length 0≤s≤s 1, the lipid phase boundary is at s = s 1, and the tubule region starts from s = s 1+1, where s 1 is chosen to be 100. We assume that membrane shape is cylindrically symmetric. ϕ(s) is the membrane tangential angle and r(s) is the radius of the tubule. The mean curvature, Ω(s), is the average of the curvatures in the radial and tangential directions; it measures the overall extent of the PIP2 head group exposure.CC BY · PubMed Central
Figure 3 Fitting of the results calculated from the model to experimental results. (A) Timelines of functional modules during endocytosis in budding yeast (continuous lines represent calculated values, and the discontinuous lines are experimental measurements—same as Figure 1A—with standard deviation). In the model, the instantaneous total levels for the individual modules (except for actin) at the endocytic site were obtained by summing their local levels over their respective locations on the membrane surface. The instantaneous total level of actin was obtained by summing over the entire bud region the product of the local actin level and its distance from cell cortex (proportional to the length of the actin filaments). To obtain the intensity plot for each of the modules, we normalized the curve for its total levels over time in accordance to its respective peak value. We then scaled the resulting curve by setting its peak value to be the same as that of the peak intensity measured experimentally. We thus can compare the computed time-lapse curve for each module to those from experimental observations. (B) Calculated endocytic membrane shape changes. The calculation of membrane shape was carried out in 3-D. Membrane shape is shown in 2-D for clarity. The parameters in the model used for curve fitting are listed in Table 1 in Protocol S1. If not stated otherwise, the parameters are fixed throughout this paper.
Figure 3 Fitting of the results calculated from the model to experimental results. (A) Timelines of functional modules during endocytosis in budding yeast (continuous lines represent calculated values, and the discontinuous lines are experimental measurements—same as Figure 1A—with standard deviation). In the model, the instantaneous total levels for the individual modules (except for actin) at the endocytic site were obtained by summing their local levels over their respective locations on the membrane surface. The instantaneous total level of actin was obtained by summing over the entire bud region the product of the local actin level and its distance from cell cortex (proportional to the length of the actin filaments). To obtain the intensity plot for each of the modules, we normalized the curve for its total levels over time in accordance to its respective peak value. We then scaled the resulting curve by setting its peak value to be the same as that of the peak intensity measured experimentally. We thus can compare the computed time-lapse curve for each module to those from experimental observations. (B) Calculated endocytic membrane shape changes. The calculation of membrane shape was carried out in 3-D. Membrane shape is shown in 2-D for clarity. The parameters in the model used for curve fitting are listed in Table 1 in Protocol S1. If not stated otherwise, the parameters are fixed throughout this paper.CC BY · PubMed Central
Figure 4 Two positive mechanochemical feedback loops between membrane shape changes and local chemical reactions. (A) Membrane tubulation by BDPs. (B) Development of interfacial forces that drive vesicle scission.
Figure 4 Two positive mechanochemical feedback loops between membrane shape changes and local chemical reactions. (A) Membrane tubulation by BDPs. (B) Development of interfacial forces that drive vesicle scission.CC BY · PubMed Central
Figure 5 Development of the interfacial force during endocytosis. (A) Schematics of interfacial forces that consist of two components. The first is the line tension. Because less PIP2 is hydrolyzed on the tubule, a higher hydrogen bond density is created adjacent to the bud. The imbalance in electrostatic attraction from hydrogen bonds between the two adjacent regions (bud and tubule) results in a line tension encircling the neck. The second force is the lateral pressure in the cytoplasmic leaflet of the bud membrane. The average area per PIP2 in the membrane is determined by the force balance between steric repulsion (i.e., arising from both the hydrocarbon chain and the polar head groups) and attractive electrostatic interactions (e.g., hydrogen bonds). The net effect of PIP2 hydrolysis is to decrease the electrostatic attraction more than the steric repulsion, causing the PIP2 leaflet to expand [56],[57]. The osmotic pressure in the cell inhibits the expansion in the normal direction, and so the cytoplasmic leaflet expands tangentially. (B) The calculated time course of the interfacial force. The threshold value for the interfacial force was determined by a force-balance calculation similar to [24]. (C) The computed time course for PIP2 levels around the lipid phase boundary (at the arc length s = 100).
Figure 5 Development of the interfacial force during endocytosis. (A) Schematics of interfacial forces that consist of two components. The first is the line tension. Because less PIP2 is hydrolyzed on the tubule, a higher hydrogen bond density is created adjacent to the bud. The imbalance in electrostatic attraction from hydrogen bonds between the two adjacent regions (bud and tubule) results in a line tension encircling the neck. The second force is the lateral pressure in the cytoplasmic leaflet of the bud membrane. The average area per PIP2 in the membrane is determined by the force balance between steric repulsion (i.e., arising from both the hydrocarbon chain and the polar head groups) and attractive electrostatic interactions (e.g., hydrogen bonds). The net effect of PIP2 hydrolysis is to decrease the electrostatic attraction more than the steric repulsion, causing the PIP2 leaflet to expand [56],[57]. The osmotic pressure in the cell inhibits the expansion in the normal direction, and so the cytoplasmic leaflet expands tangentially. (B) The calculated time course of the interfacial force. The threshold value for the interfacial force was determined by a force-balance calculation similar to [24]. (C) The computed time course for PIP2 levels around the lipid phase boundary (at the arc length s = 100).CC BY · PubMed Central
Figure 6 Phase diagrams for endocytic dynamics. The shaded areas represent the parameter regions for successful endocytosis; the star in each phase diagram represents the parameter set used in the fitting plot in Figure 3. (A) Strength of BDP PIP2 protection: K 2 versus curvature-dependent PIP2 hydrolysis rate k 2; (B) Curvature-dependent factor for phosphatase recruitment rate, α versus phosphatase recruitment rate k 3; (C) Relative rate of BDP dynamics versus actin polymerization rate k 7; (D) Curvature-dependent factor of BDP recruitment rate χ versus interfacial force constant λ 0. Each phenotype is characterized by: (a) time-lapse plot for the coat proteins (red), actin (blue), BDP (green), phosphatase (orange), and the membrane tip position (black); (b) the time course for interfacial force development (purple); (c) the time course for membrane shape change (black). The intensity of each functional module in the phenotype plots is normalized relative to the corresponding wild-type normalized intensity shown in Figure 3, thus representing the relative abundance. Phenotype 1: Without PIP2 hydrolysis [k2 reduces from 20 (nm) per second to 0]. Phenotype 2: Increased protection strength of PIP2 hydrolysis at the tubule region [ increases from to ]. Phenotype 3: Increased phosphatase recruitment rate [α increases from 100 nm to 500 nm]. Phenotype 4: BDP recruitment does not occur.
Figure 6 Phase diagrams for endocytic dynamics. The shaded areas represent the parameter regions for successful endocytosis; the star in each phase diagram represents the parameter set used in the fitting plot in Figure 3. (A) Strength of BDP PIP2 protection: K 2 versus curvature-dependent PIP2 hydrolysis rate k 2; (B) Curvature-dependent factor for phosphatase recruitment rate, α versus phosphatase recruitment rate k 3; (C) Relative rate of BDP dynamics versus actin polymerization rate k 7; (D) Curvature-dependent factor of BDP recruitment rate χ versus interfacial force constant λ 0. Each phenotype is characterized by: (a) time-lapse plot for the coat proteins (red), actin (blue), BDP (green), phosphatase (orange), and the membrane tip position (black); (b) the time course for interfacial force development (purple); (c) the time course for membrane shape change (black). The intensity of each functional module in the phenotype plots is normalized relative to the corresponding wild-type normalized intensity shown in Figure 3, thus representing the relative abundance. Phenotype 1: Without PIP2 hydrolysis [k2 reduces from 20 (nm) per second to 0]. Phenotype 2: Increased protection strength of PIP2 hydrolysis at the tubule region [ increases from to ]. Phenotype 3: Increased phosphatase recruitment rate [α increases from 100 nm to 500 nm]. Phenotype 4: BDP recruitment does not occur.CC BY · PubMed Central
Figure 7 Schematics comparing endocytosis in yeast and mammalian cells. (A) Model for yeast endocytosis. (B) Model for mammalian endocytosis (see text).
Figure 7 Schematics comparing endocytosis in yeast and mammalian cells. (A) Model for yeast endocytosis. (B) Model for mammalian endocytosis (see text).CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 59 references

Cited by

Sorted by citations

Metadata

Authors
Liu, J., Sun, Y., Drubin, D. G., & Oster, G. F.
Affiliations
University of California, Berkeley (US)
Corresponding author
  • Drubin, David G.
  • Oster, George F.
Venue
PLoS Biology
Published
2009-09-29
Volume
7
Issue
9
Pages
e1000204
ISSN
  • 1544-9173 (ISSN-L)
  • 1545-7885
OpenAlex
W2019040974
Cited by
222
References
59
Citations by year
2026: 1 · 2025: 7 · 2024: 8 · 2023: 6 · 2022: 9 · 2021: 11 · 2020: 6 · 2019: 9 · 2018: 15 · 2017: 15 · 2016: 10 · 2015: 16 · 2014: 21 · 2013: 17 · 2012: 24
FWCI
10.53
Type
journal article
Language
en
Open access
Locations
  • PLoS Biology (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion)
  • Figshare (submittedVersion, cc-by)
Funders
  • National Science Foundation (0414039)
  • National Institute of General Medical Sciences (1 R01 GM 50399, R01GM042759, R01GM050399, 1 R01 GM 42759)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Cell Biology › Cellular transport and secretion
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Lipid Membrane Structure and Behavior
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Fungal and yeast genetics research
Keywords
endocytosis, mechanochemistry, vesicle scission, membrane curvature, protein recruitment, lipid modification, interfacial forces, mechanochemical model, phase diagram, vesicle formation
MeSH terms
Actins, Animals, Biomechanical Phenomena, Cell Membrane, Endocytosis, Saccharomycetales, Hydrolysis, Lipids, Mammals, Models, Biological, Phosphoric Monoester Hydrolases, Time Factors, Protein Structure, Tertiary, Phosphatidylinositol 4,5-Diphosphate, Transport Vesicles, Feedback, Physiological, Saccharomyces cerevisiae Proteins

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review