1
|
Zhang X, Kang R, Liu Y, Yan Z, Xu Y, Yue T. From reversible to irreversible: When the membrane nanotube pearling is coupled with phase separation. Colloids Surf B Biointerfaces 2021; 209:112160. [PMID: 34736219 DOI: 10.1016/j.colsurfb.2021.112160] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2021] [Revised: 09/17/2021] [Accepted: 10/07/2021] [Indexed: 12/11/2022]
Abstract
Membrane nanotubes, which are ubiquitous in biology and act as channels maintaining transport between different cells and organelles, readily undergo pearling in response to external stimuli. Membrane nanotube pearling involves generation of heterogeneous curvature coupled with redistribution of membrane components that may interfere with the shape recovery of pearled nanotubes. However, the mechanism underlying such delicate process remains unclear and difficult to study at the molecular scale in vivo. By means of molecular dynamics simulation, here we investigate pearling of multi-component membrane nanotubes and reversibility through manipulating system temperature and osmotic pressure. With the equilibrium shape of membrane nanotubes controlled by the osmotic pressure, our results demonstrate that the process of membrane nanotube pearling can be reversible or irreversible, depending on the phase segregation state. For the pearled nanotube releasing high surface energy, different lipid components redistribute along the tube axial direction. Lipids with unsaturated tails prefer gathering at the high-curvature shrinking region, whereas the swelling region is constituted by saturated lipids forming the liquid-ordered phase of a higher bending rigidity. Such curvature sensitive phase segregation minimizes the system free energy by reducing both the membrane bending energy and line tension at the phase boundary. As such, the pearled nanotube fails to recover its shape upon retracting stimuli, suggesting irreversibility of the membrane nanotube pearling coupled with phase separation. Given importance of membrane nanotube pearling in various cellular activities, these results provide a new mechanism of controlling equilibrium shapes of membrane nanotubes in complex cellular environment.
Collapse
Affiliation(s)
- Xiaoyang Zhang
- College of Chemical Engineering, China University of Petroleum (East China), Qingdao 266580, China
| | - Runshan Kang
- College of Chemical Engineering, China University of Petroleum (East China), Qingdao 266580, China
| | - Yingjie Liu
- Institute of Coastal Environmental Pollution Control, Key Laboratory of Marine Environment and Ecology, Ministry of Education, Frontiers Science Center for Deep Ocean Multispheres and Earth System, Ocean University of China, Qingdao 266100, China
| | - Zengshuai Yan
- National Laboratory of Solid State Microstructures, School of Physics, Nanjing University, Nanjing 210093, China
| | - Yan Xu
- College of Electronic Engineering and Automation, Shandong University of Science and Technology, Qingdao 266590, China
| | - Tongtao Yue
- Institute of Coastal Environmental Pollution Control, Key Laboratory of Marine Environment and Ecology, Ministry of Education, Frontiers Science Center for Deep Ocean Multispheres and Earth System, Ocean University of China, Qingdao 266100, China; Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266237, China.
| |
Collapse
|
2
|
Knyazev E, Nersisyan S, Tonevitsky A. Endocytosis and Transcytosis of SARS-CoV-2 Across the Intestinal Epithelium and Other Tissue Barriers. Front Immunol 2021; 12:636966. [PMID: 34557180 PMCID: PMC8452982 DOI: 10.3389/fimmu.2021.636966] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2020] [Accepted: 08/20/2021] [Indexed: 12/14/2022] Open
Abstract
Since 2003, the world has been confronted with three new betacoronaviruses that cause human respiratory infections: SARS-CoV, which causes severe acute respiratory syndrome (SARS), MERS-CoV, which causes Middle East respiratory syndrome (MERS), and SARS-CoV-2, which causes Coronavirus Disease 2019 (COVID-19). The mechanisms of coronavirus transmission and dissemination in the human body determine the diagnostic and therapeutic strategies. An important problem is the possibility that viral particles overcome tissue barriers such as the intestine, respiratory tract, blood-brain barrier, and placenta. In this work, we will 1) consider the issue of endocytosis and the possibility of transcytosis and paracellular trafficking of coronaviruses across tissue barriers with an emphasis on the intestinal epithelium; 2) discuss the possibility of antibody-mediated transcytosis of opsonized viruses due to complexes of immunoglobulins with their receptors; 3) assess the possibility of the virus transfer into extracellular vesicles during intracellular transport; and 4) describe the clinical significance of these processes. Models of the intestinal epithelium and other barrier tissues for in vitro transcytosis studies will also be briefly characterized.
Collapse
Affiliation(s)
- Evgeny Knyazev
- Laboratory of Microfluidic Technologies for Biomedicine, Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry of the Russian Academy of Sciences, Moscow, Russia.,Faculty of Biology and Biotechnology, National Research University Higher School of Economics (HSE), Moscow, Russia
| | - Stepan Nersisyan
- Faculty of Biology and Biotechnology, National Research University Higher School of Economics (HSE), Moscow, Russia
| | - Alexander Tonevitsky
- Laboratory of Microfluidic Technologies for Biomedicine, Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry of the Russian Academy of Sciences, Moscow, Russia.,Faculty of Biology and Biotechnology, National Research University Higher School of Economics (HSE), Moscow, Russia
| |
Collapse
|
3
|
van Hilten N, Stroh KS, Risselada HJ. Membrane Thinning Induces Sorting of Lipids and the Amphipathic Lipid Packing Sensor (ALPS) Protein Motif. Front Physiol 2020; 11:250. [PMID: 32372966 PMCID: PMC7177014 DOI: 10.3389/fphys.2020.00250] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2020] [Accepted: 03/05/2020] [Indexed: 12/16/2022] Open
Abstract
Heterogeneities (e.g., membrane proteins and lipid domains) and deformations (e.g., highly curved membrane regions) in biological lipid membranes cause lipid packing defects that may trigger functional sorting of lipids and membrane-associated proteins. To study these phenomena in a controlled and efficient way within molecular simulations, we developed an external field protocol that artificially enhances packing defects in lipid membranes by enforcing local thinning of a flat membrane region. For varying lipid compositions, we observed strong thinning-induced depletion or enrichment, depending on the lipid's intrinsic shape and its effect on a membrane's elastic modulus. In particular, polyunsaturated and lysolipids are strongly attracted to regions high in packing defects, whereas phosphatidylethanolamine (PE) lipids and cholesterol are strongly repelled from it. Our results indicate that externally imposed changes in membrane thickness, area, and curvature are underpinned by shared membrane elastic principles. The observed sorting toward the thinner membrane region is in line with the sorting expected for a positively curved membrane region. Furthermore, we have demonstrated that the amphipathic lipid packing sensor (ALPS) protein motif, a known curvature and packing defect sensor, is effectively attracted to thinner membrane regions. By extracting the force that drives amphipathic molecules toward the thinner region, our thinning protocol can directly quantify and score the lipid packing sensing of different amphipathic molecules. In this way, our protocol paves the way toward high-throughput exploration of potential defect- and curvature-sensing motifs, making it a valuable addition to the molecular simulation toolbox.
Collapse
Affiliation(s)
- Niek van Hilten
- Leiden Institute of Chemistry, Leiden University, Leiden, Netherlands
| | - Kai Steffen Stroh
- Institute for Theoretical Physics, Georg August University Göttingen, Göttingen, Germany
| | - Herre Jelger Risselada
- Leiden Institute of Chemistry, Leiden University, Leiden, Netherlands.,Institute for Theoretical Physics, Georg August University Göttingen, Göttingen, Germany
| |
Collapse
|
4
|
Fonda P, Rinaldin M, Kraft DJ, Giomi L. Thermodynamic equilibrium of binary mixtures on curved surfaces. Phys Rev E 2019; 100:032604. [PMID: 31639923 DOI: 10.1103/physreve.100.032604] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2019] [Indexed: 06/10/2023]
Abstract
We study the global influence of curvature on the free energy landscape of two-dimensional binary mixtures confined on closed surfaces. Starting from a generic effective free energy, constructed on the basis of symmetry considerations and conservation laws, we identify several model-independent phenomena, such as a curvature-dependent line tension and local shifts in the binodal concentrations. To shed light on the origin of the phenomenological parameters appearing in the effective free energy, we further construct a lattice-gas model of binary mixtures on nontrivial substrates, based on the curved-space generalization of the two-dimensional Ising model. This allows us to decompose the interaction between the local concentration of the mixture and the substrate curvature into four distinct contributions, as a result of which the phase diagram splits into critical subdiagrams. The resulting free energy landscape can admit, as stable equilibria, strongly inhomogeneous mixed phases, which we refer to as "antimixed" states below the critical temperature. We corroborate our semianalytical findings with phase-field numerical simulations on realistic curved lattices. Despite this work being primarily motivated by recent experimental observations of multicomponent lipid vesicles supported by colloidal scaffolds, our results are applicable to any binary mixture confined on closed surfaces of arbitrary geometry.
Collapse
Affiliation(s)
- Piermarco Fonda
- Instituut-Lorentz, Universiteit Leiden, P.O. Box 9506, 2300 RA Leiden, Netherlands
| | - Melissa Rinaldin
- Instituut-Lorentz, Universiteit Leiden, P.O. Box 9506, 2300 RA Leiden, Netherlands
- Huygens-Kamerlingh Onnes Laboratory, Universiteit Leiden, P.O. Box 9504, 2300 RA Leiden, Netherlands
| | - Daniela J Kraft
- Huygens-Kamerlingh Onnes Laboratory, Universiteit Leiden, P.O. Box 9504, 2300 RA Leiden, Netherlands
| | - Luca Giomi
- Instituut-Lorentz, Universiteit Leiden, P.O. Box 9506, 2300 RA Leiden, Netherlands
| |
Collapse
|
5
|
Beltrán-Heredia E, Tsai FC, Salinas-Almaguer S, Cao FJ, Bassereau P, Monroy F. Membrane curvature induces cardiolipin sorting. Commun Biol 2019; 2:225. [PMID: 31240263 PMCID: PMC6586900 DOI: 10.1038/s42003-019-0471-x] [Citation(s) in RCA: 86] [Impact Index Per Article: 17.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2018] [Accepted: 05/23/2019] [Indexed: 11/17/2022] Open
Abstract
Cardiolipin is a cone-shaped lipid predominantly localized in curved membrane sites of bacteria and in the mitochondrial cristae. This specific localization has been argued to be geometry-driven, since the CL's conical shape relaxes curvature frustration. Although previous evidence suggests a coupling between CL concentration and membrane shape in vivo, no precise experimental data are available for curvature-based CL sorting in vitro. Here, we test this hypothesis in experiments that isolate the effects of membrane curvature in lipid-bilayer nanotubes. CL sorting is observed with increasing tube curvature, reaching a maximum at optimal CL concentrations, a fact compatible with self-associative clustering. Observations are compatible with a model of membrane elasticity including van der Waals entropy, from which a negative intrinsic curvature of -1.1 nm-1 is predicted for CL. The results contribute to understanding the physicochemical interplay between membrane curvature and composition, providing key insights into mitochondrial and bacterial membrane organization and dynamics.
Collapse
Affiliation(s)
- Elena Beltrán-Heredia
- Departamento de Estructura de la Materia, Física Térmica y Electrónica, Universidad Complutense de Madrid, Plaza de Ciencias, 1, 28040 Madrid, Spain
- Departamento de Química Física, Universidad Complutense de Madrid, Avda. Complutense, s/n, 28040 Madrid, Spain
| | - Feng-Ching Tsai
- Laboratoire Physico Chimie Curie, Institut Curie, PSL Research University, CNRS UMR168, 75005 Paris, France
| | - Samuel Salinas-Almaguer
- Departamento de Química Física, Universidad Complutense de Madrid, Avda. Complutense, s/n, 28040 Madrid, Spain
| | - Francisco J. Cao
- Departamento de Estructura de la Materia, Física Térmica y Electrónica, Universidad Complutense de Madrid, Plaza de Ciencias, 1, 28040 Madrid, Spain
- Instituto Madrileño de Estudios Avanzados en Nanociencia, IMDEA Nanociencia, Calle Faraday, 9, 28049 Madrid, Spain
| | - Patricia Bassereau
- Laboratoire Physico Chimie Curie, Institut Curie, PSL Research University, CNRS UMR168, 75005 Paris, France
- Sorbonne Université, UPMC Univ Paris 06, 75005 Paris, France
| | - Francisco Monroy
- Departamento de Química Física, Universidad Complutense de Madrid, Avda. Complutense, s/n, 28040 Madrid, Spain
- Unit of Translational Biophysics, Instituto de Investigación Sanitaria Hospital Doce de Octubre (imas12), Avda. de Córdoba, s/n, 28041 Madrid, Spain
- Institute for Quantitative Biosciences-QB3, University of California at Berkeley, Berkeley, CA 94720 USA
| |
Collapse
|
6
|
Wu Z, Yuan H, Zhang X, Yi X. Sidewall contact regulating the nanorod packing inside vesicles with relative volumes. SOFT MATTER 2019; 15:2552-2559. [PMID: 30839980 DOI: 10.1039/c8sm01656a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Intracellular packing of one-dimensional and rodlike materials plays an important role in many biological processes such as cell mimicking, microtubule protrusion, cell division, frustrated phagocytosis, and pathogenicity. To understand the mechanical interplay between cells/intracellular membranous organelles and encapsulated rodlike materials, we perform theoretical analyses to investigate how the morphologies and mechanical behaviors of lipid vesicles of given relative volumes are regulated by encapsulated rigid nanorods of finite diameters and selected geometries, including a cylindrical nanorod, a nanorod with one widened end, and a cone-shaped nanorod. The contact between the vesicle protrusion and the sidewall of the rod, neglected in most theoretical studies, is shown to play an important role in regulating vesicle tubulation, membrane tension, and axial contact force on the nanorod. As the nanorod length increases, the confining vesicle evolves from a prolate into different shapes, such as a lemon, a conga drum, a cherry, and a bowling pin, depending on the radical size of the nanorod and the relative vesicle volume. The corresponding morphological phase diagrams are determined. Moreover, phase diagrams of the buckling of the encapsulated nanorods are determined based on the classical Euler buckling theory. It is shown that there exists an optimal filament number at which the encapsulated weakly cross-linked filament bundle maintains the largest length in a mechanically stable state. Similarities and differences between the nanorod packing in vesicles at a given pressure difference and a relative volume are discussed. Our results provide valuable insight into the biophysics underlying cell interactions with one-dimensional and rodlike materials.
Collapse
Affiliation(s)
- Zeming Wu
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing 100871, China.
| | | | | | | |
Collapse
|
7
|
Luo Y, Maibaum L. Relating the structure factors of two-dimensional materials in planar and spherical geometries. SOFT MATTER 2018; 14:5686-5692. [PMID: 29947410 DOI: 10.1039/c8sm00978c] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Scattering structure factors provide essential insight into material properties and are routinely obtained in experiments, computer simulations, and theoretical analyses. Different approaches favor different geometries of the material. In case of lipid bilayers, scattering experiments can be performed on spherical vesicles, while simulations and theory often consider planar membrane patches. We derive an approximate relationship between the structure functions of such a material in planar and spherical geometries. We illustrate the usefulness of this relationship in a case study of a Gaussian material that supports both homogeneous and microemulsion phases. Within its range of applicability, this relationship enables a model-free comparison of structure factors of the same material in different geometries.
Collapse
Affiliation(s)
- Yongtian Luo
- Department of Chemistry, University of Washington, Seattle, WA 98195, USA.
| | | |
Collapse
|
8
|
Tian F, Yue T, Dong W, Yi X, Zhang X. Size-dependent formation of membrane nanotubes: continuum modeling and molecular dynamics simulations. Phys Chem Chem Phys 2018; 20:3474-3483. [DOI: 10.1039/c7cp06212e] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
With continuum theory and molecular dynamics simulations we demonstrated that the lipid membrane upon extraction exhibits size- and tension-dependent mechanical behaviors, and different structural lipid rearrangements in different leaflets.
Collapse
Affiliation(s)
- Falin Tian
- State Key Laboratory of Organic-Inorganic Composites
- Beijing University of Chemical Technology
- Beijing 100029
- China
- Laboratoire de Chimie
| | - Tongtao Yue
- State Key Laboratory of Heavy Oil Processing
- Center for Bioengineering and Biotechnology
- China University of Petroleum (East China)
- Qingdao 266580
- China
| | - Wei Dong
- Laboratoire de Chimie
- Ecole Normale Superieure de Lyon
- 69364 Lyon Cedex 07
- France
| | - Xin Yi
- Department of Mechanics and Engineering Science
- College of Engineering
- Peking University
- Beijing 100871
- China
| | - Xianren Zhang
- State Key Laboratory of Organic-Inorganic Composites
- Beijing University of Chemical Technology
- Beijing 100029
- China
| |
Collapse
|
9
|
Yu S, Wang H, Ni Y, He L, Huang M, Lin Y, Qian J, Jiang H. Tuning interfacial patterns of molecular bonds via surface morphology. SOFT MATTER 2017; 13:5970-5976. [PMID: 28869265 DOI: 10.1039/c7sm01278k] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Many studies have demonstrated that the mechanical properties of the extracellular matrix can significantly influence the morphology, strength and lifetime of focal adhesions. However, how the morphology of the contact surface affects the pattern formation of the molecular bonds still remains largely unknown. Here, by simplifying the cell and extracellular matrix to two opposing elastic bodies and considering the lateral diffusion as well as the bonding/debonding of molecular bonds, we study the clustering behavior of receptor-ligand bonds between curved surfaces and the phase diagrams of cluster patterns. We reveal the important role of surface morphology and bond kinetics in regulating the patterns of bond clusters. We further investigate the segregation dynamics of the interfacial bonds under various loading speeds, and we show that the average interfacial stress is rate-dependent while the rupture stress is rate-independent. Finally, we demonstrate that programmable patterning of bond clusters can be achieved through the designed surface morphology.
Collapse
Affiliation(s)
- Sai Yu
- CAS Key Laboratory of Mechanical Behavior and Design of Materials, Department of Modern Mechanics, University of Science and Technology of China, Hefei, Anhui 230026, China.
| | | | | | | | | | | | | | | |
Collapse
|
10
|
Zhu Q, Zheng F, Liu AP, Qian J, Fu C, Lin Y. Shape Transformation of the Nuclear Envelope during Closed Mitosis. Biophys J 2017; 111:2309-2316. [PMID: 27851952 DOI: 10.1016/j.bpj.2016.10.004] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2016] [Revised: 09/28/2016] [Accepted: 10/07/2016] [Indexed: 11/18/2022] Open
Abstract
The nuclear envelope (NE) in lower eukaryotes such as Schizosaccharomyces pombe undergoes large morphology changes during closed mitosis. However, which physical parameters are important in governing the shape evolution of the NE, and how defects in the dividing chromosomes/microtubules are reflected in those parameters, are fundamental questions that remain unresolved. In this study, we show that improper separation of chromosomes in genetically deficient cells leads to membrane tethering or asymmetric division in contrast to the formation of two equal-sized daughter nuclei in wild-type cells. We hypothesize that the poleward force is transmitted to the nuclear membrane through its physical contact with the separated sister chromatids at the two spindle poles. A theoretical model is developed to predict the morphology evolution of the NE where key factors such as the work done by the poleward force and bending and surface energies stored in the membrane have been taken into account. Interestingly, the predicted phase diagram, summarizing the dependence of nuclear shape on the size of the load transmission regions, and the pole-to-pole distance versus surface area relationship all quantitatively agree well with our experimental observations, suggesting that this model captures the essential physics involved in closed mitosis.
Collapse
Affiliation(s)
- Qian Zhu
- Department of Mechanical Engineering, The University of Hong Kong, Hong Kong SAR, China
| | - Fan Zheng
- Chinese Academy of Sciences Center for Excellence in Molecular Cell Biology, School of Life Sciences, University of Science and Technology of China, Hefei, Anhui, China
| | - Allen P Liu
- Department of Mechanical Engineering, University of Michigan, Ann Arbor, Michigan
| | - Jin Qian
- Department of Engineering Mechanics, Zhejiang University, Hangzhou, Zhejiang, China
| | - Chuanhai Fu
- Chinese Academy of Sciences Center for Excellence in Molecular Cell Biology, School of Life Sciences, University of Science and Technology of China, Hefei, Anhui, China.
| | - Yuan Lin
- Department of Mechanical Engineering, The University of Hong Kong, Hong Kong SAR, China.
| |
Collapse
|
11
|
Brummel BE, Braun AR, Sachs JN. Polyunsaturated chains in asymmetric lipids disorder raft mixtures and preferentially associate with α-Synuclein. BIOCHIMICA ET BIOPHYSICA ACTA. BIOMEMBRANES 2017; 1859:529-536. [PMID: 27742354 PMCID: PMC5749632 DOI: 10.1016/j.bbamem.2016.10.006] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/30/2016] [Revised: 09/30/2016] [Accepted: 10/05/2016] [Indexed: 01/12/2023]
Abstract
Using molecular dynamics simulations, we have explored the effect of asymmetric lipids-specifically those that contain one polyunsaturated (PUFA) and one saturated fatty acid chain-on phase separation in heterogeneous membranes. These lipids are prevalent in neuronal membranes, particularly in synaptic membranes, where the Parkinson's Disease protein α-Synuclein (αS) is found. We have therefore explored the relationship between asymmetric, PUFA-containing lipids, and αS. The simulations show that asymmetric lipids partition to the liquid disordered (Ld) phase of canonical raft mixtures because of the highly disordered PUFA chain. In the case of a membrane built to mimic the lipid composition of a synaptic vesicle, the PUFA-containing asymmetric lipids completely disrupt phase separation. Because αS is positively charged, we show that it partitions with negatively charged lipids, regardless of the saturation state of the chains. Additionally, αS preferentially associates with the polyunsaturated fatty acid tails of both charged and neutral lipids. This is a consequence of those chains' ability to accommodate the void beneath the amphipathic helix. This article is part of a Special Issue entitled: Lipid order/lipid defects and lipid-control of protein activity edited by Dirk Schneider.
Collapse
Affiliation(s)
- Benjamin E Brummel
- Department of Biomedical Engineering, University of Minnesota, 7-124 Nils Hasselmo Hall, 312 Church St SE, Minneapolis, MN 55455, USA.
| | - Anthony R Braun
- Department of Biomedical Engineering, University of Minnesota, 7-124 Nils Hasselmo Hall, 312 Church St SE, Minneapolis, MN 55455, USA.
| | - Jonathan N Sachs
- Department of Biomedical Engineering, University of Minnesota, 7-124 Nils Hasselmo Hall, 312 Church St SE, Minneapolis, MN 55455, USA.
| |
Collapse
|
12
|
Klaus CJS, Raghunathan K, DiBenedetto E, Kenworthy AK. Analysis of diffusion in curved surfaces and its application to tubular membranes. Mol Biol Cell 2016; 27:3937-3946. [PMID: 27733625 PMCID: PMC5170615 DOI: 10.1091/mbc.e16-06-0445] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2016] [Revised: 09/08/2016] [Accepted: 10/04/2016] [Indexed: 11/11/2022] Open
Abstract
Diffusion of particles in curved surfaces is inherently complex compared with diffusion in a flat membrane, owing to the nonplanarity of the surface. The consequence of such nonplanar geometry on diffusion is poorly understood but is highly relevant in the case of cell membranes, which often adopt complex geometries. To address this question, we developed a new finite element approach to model diffusion on curved membrane surfaces based on solutions to Fick's law of diffusion and used this to study the effects of geometry on the entry of surface-bound particles into tubules by diffusion. We show that variations in tubule radius and length can distinctly alter diffusion gradients in tubules over biologically relevant timescales. In addition, we show that tubular structures tend to retain concentration gradients for a longer time compared with a comparable flat surface. These findings indicate that sorting of particles along the surfaces of tubules can arise simply as a geometric consequence of the curvature without any specific contribution from the membrane environment. Our studies provide a framework for modeling diffusion in curved surfaces and suggest that biological regulation can emerge purely from membrane geometry.
Collapse
Affiliation(s)
| | - Krishnan Raghunathan
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN 37232
| | | | - Anne K Kenworthy
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN 37232 .,Department of Cell and Developmental Biology, Vanderbilt University School of Medicine, Nashville, TN 37232
| |
Collapse
|
13
|
Su YC, Chen JZY. A model of vesicle tubulation and pearling induced by adsorbing particles. SOFT MATTER 2015; 11:4054-4060. [PMID: 25907594 DOI: 10.1039/c5sm00565e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
We study the basic theoretical model of a deformable vesicle immersed in a solution of particles that can adsorb onto one of the two surfaces of a membrane. The model consists of an adsorption energy gain for the adsorbing particles and the Canham-Helfrich membrane bending energy, in which the spontaneous curvature is coupled with the adsorption area. We demonstrate that bud, pearling, and tube conformations can be stabilized after minimizing the free energy and that the pearling-tubulation transition has the characteristics of an abrupt structural transition.
Collapse
Affiliation(s)
- Yu-Cheng Su
- Department of Physics and Astronomy, University of Waterloo, Waterloo, Ontario, CanadaN2L 3G1.
| | | |
Collapse
|
14
|
Jiang H, Qian J, Lin Y, Ni Y, He L. Aggregation dynamics of molecular bonds between compliant materials. SOFT MATTER 2015; 11:2812-2820. [PMID: 25706682 DOI: 10.1039/c4sm02903h] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
In this paper, we develop a mechanochemical modeling framework in which the spatial-temporal evolution of receptor-ligand bonds takes place at the interface between two compliant media in the presence of an externally applied tensile load. Bond translocation, dissociation and association occur simultaneously, resulting in dynamic aggregation of molecular bonds that is regulated by mechanical factors such as material compliance and applied stress. The results show that bond aggregation is energetically favorable in the out-of-equilibrium process with convoluted time scales from bond diffusion and reaction. Material stiffness is predicted to contribute to adhesion growth and an optimal level of applied stress leads to the maximized size of bond clusters for integrin-based adhesion, consistent with related experimental observations on focal adhesions of cell-matrix interactions. The stress distribution within bond clusters is generally non-uniform and governed by the stress concentration index.
Collapse
Affiliation(s)
- Hongyuan Jiang
- Department of Modern Mechanics, CAS Key Laboratory of Mechanical Behavior and Design of Materials, University of Science and Technology of China, Hefei, Anhui 230026, China
| | | | | | | | | |
Collapse
|
15
|
Ramakrishnan N, Sunil Kumar PB, Radhakrishnan R. Mesoscale computational studies of membrane bilayer remodeling by curvature-inducing proteins. PHYSICS REPORTS 2014; 543:1-60. [PMID: 25484487 PMCID: PMC4251917 DOI: 10.1016/j.physrep.2014.05.001] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
Abstract
Biological membranes constitute boundaries of cells and cell organelles. These membranes are soft fluid interfaces whose thermodynamic states are dictated by bending moduli, induced curvature fields, and thermal fluctuations. Recently, there has been a flood of experimental evidence highlighting active roles for these structures in many cellular processes ranging from trafficking of cargo to cell motility. It is believed that the local membrane curvature, which is continuously altered due to its interactions with myriad proteins and other macromolecules attached to its surface, holds the key to the emergent functionality in these cellular processes. Mechanisms at the atomic scale are dictated by protein-lipid interaction strength, lipid composition, lipid distribution in the vicinity of the protein, shape and amino acid composition of the protein, and its amino acid contents. The specificity of molecular interactions together with the cooperativity of multiple proteins induce and stabilize complex membrane shapes at the mesoscale. These shapes span a wide spectrum ranging from the spherical plasma membrane to the complex cisternae of the Golgi apparatus. Mapping the relation between the protein-induced deformations at the molecular scale and the resulting mesoscale morphologies is key to bridging cellular experiments across the various length scales. In this review, we focus on the theoretical and computational methods used to understand the phenomenology underlying protein-driven membrane remodeling. Interactions at the molecular scale can be computationally probed by all atom and coarse grained molecular dynamics (MD, CGMD), as well as dissipative particle dynamics (DPD) simulations, which we only describe in passing. We choose to focus on several continuum approaches extending the Canham - Helfrich elastic energy model for membranes to include the effect of curvature-inducing proteins and explore the conformational phase space of such systems. In this description, the protein is expressed in the form of a spontaneous curvature field. The approaches include field theoretical methods limited to the small deformation regime, triangulated surfaces and particle-based computational models to investigate the large-deformation regimes observed in the natural state of many biological membranes. Applications of these methods to understand the properties of biological membranes in homogeneous and inhomogeneous environments of proteins, whose underlying curvature fields are either isotropic or anisotropic, are discussed. The diversity in the curvature fields elicits a rich variety of morphological states, including tubes, discs, branched tubes, and caveola. Mapping the thermodynamic stability of these states as a function of tuning parameters such as concentration and strength of curvature induction of the proteins is discussed. The relative stabilities of these self-organized shapes are examined through free-energy calculations. The suite of methods discussed here can be tailored to applications in specific cellular settings such as endocytosis during cargo trafficking and tubulation of filopodial structures in migrating cells, which makes these methods a powerful complement to experimental studies.
Collapse
Affiliation(s)
- N. Ramakrishnan
- Department of Chemical and Biomolecular Engineering, Department of Bioengineering, Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, PA-19104
| | - P. B. Sunil Kumar
- Department of Physics, Indian Institute of Technology Madras, Chennai, India - 600036
| | - Ravi Radhakrishnan
- Department of Chemical and Biomolecular Engineering, Department of Bioengineering, Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, PA-19104
| |
Collapse
|
16
|
Tourdot RW, Ramakrishnan N, Radhakrishnan R. Defining the free-energy landscape of curvature-inducing proteins on membrane bilayers. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:022717. [PMID: 25215768 PMCID: PMC4336182 DOI: 10.1103/physreve.90.022717] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/05/2014] [Indexed: 05/08/2023]
Abstract
Curvature-sensing and curvature-remodeling proteins, such as Amphiphysin, Epsin, and Exo70, are known to reshape cell membranes, and this remodeling event is essential for key biophysical processes such as tubulation, exocytosis, and endocytosis. Curvature-inducing proteins can act as curvature sensors; they aggregate to membrane regions matching their intrinsic curvature; as well as induce curvature in cell membranes to stabilize emergent high curvature, nonspherical, structures such as tubules, discs, and caveolae. A definitive understanding of the interplay between protein recruitment and migration, the evolution of membrane curvature, and membrane morphological transitions is emerging but remains incomplete. Here, within a continuum framework and using the machinery of Monte Carlo simulations, we introduce and compare three free-energy methods to delineate the free-energy landscape of curvature-inducing proteins on bilayer membranes. We demonstrate the utility of the Widom test particle (or field) insertion methodology in computing the excess chemical potentials associated with curvature-inducing proteins on the membrane-in particular, we use this method to track the onset of morphological transitions in the membrane at elevated protein densities. We validate this approach by comparing the results from the Widom method with those of thermodynamic integration and Bennett acceptance ratio methods. Furthermore, the predictions from the Widom method have been tested against analytical calculations of the excess chemical potential at infinite dilution. Our results are useful in precisely quantifying the free-energy landscape, and also in determining the phase boundaries associated with curvature-induction, curvature-sensing, and morphological transitions. This approach can be extended to studies exploring the role of thermal fluctuations and other external (control) variables, such as membrane excess area, in shaping curvature-mediated interactions on bilayer membranes.
Collapse
Affiliation(s)
- Richard W. Tourdot
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - N. Ramakrishnan
- Department of Bioengineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Ravi Radhakrishnan
- Department of Chemical and Biomolecular Engineering, Department of Bioengineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
17
|
Weflen AW, Baier N, Tang QJ, Van den Hof M, Blumberg RS, Lencer WI, Massol RH. Multivalent immune complexes divert FcRn to lysosomes by exclusion from recycling sorting tubules. Mol Biol Cell 2013; 24:2398-405. [PMID: 23741050 PMCID: PMC3727932 DOI: 10.1091/mbc.e13-04-0174] [Citation(s) in RCA: 65] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
Study of receptor sorting between recycling and degradative pathways shows that sorting into the recycling pathway depends not only on recognition of sorting motifs by cytosolic adaptors, but also on the physical properties of the endosomal luminal complexes, as shown by the neonatal receptor for IgG FcRn. The neonatal receptor for immunoglobulin G (IgG; FcRn) prevents IgG degradation by efficiently sorting IgG into recycling endosomes and away from lysosomes. When bound to IgG-opsonized antigen complexes, however, FcRn traffics cargo into lysosomes, where antigen processing can occur. Here we address the mechanism of sorting when FcRn is bound to multivalent IgG-opsonized antigens. We find that only the unbound receptor or FcRn bound to monomeric IgG is sorted into recycling tubules emerging from early endosomes. Cross-linked FcRn is never visualized in tubules containing the unbound receptor. Similar results are found for transferrin receptor, suggesting a general mechanism of action. Deletion or replacement of the FcRn cytoplasmic tail does not prevent diversion of trafficking to lysosomes upon cross-linking. Thus physical properties of the lumenal ligand–receptor complex appear to act as key determinants for sorting between the recycling and lysosomal pathways by regulating FcRn entry into recycling tubules.
Collapse
Affiliation(s)
- Andrew W Weflen
- Division of Gastroenterology, Boston Children's Hospital, Boston, MA 02115, USA
| | | | | | | | | | | | | |
Collapse
|
18
|
Jiang H. Dynamic sorting of lipids and proteins in multicomponent membranes. PHYSICAL REVIEW LETTERS 2012; 109:198101. [PMID: 23215429 DOI: 10.1103/physrevlett.109.198101] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/18/2012] [Indexed: 06/01/2023]
Abstract
Dynamic sorting of lipids and proteins in cellular membranes plays a critical role in establishing and maintaining distinct compositions in various organelles. Recent experiments found that the lipid sorting in a membrane tube highly depends on the pulling speed at the tip. However, the mechanism of this velocity dependence has not yet been revealed. In this Letter, we found that when a membrane is deformed rapidly, the lipid flow induced by fast membrane shape change will significantly affect the sorting results. The competition between the curvature-driven lipid sorting and the pulling-induced lipid flow leads to novel behaviors. When a membrane tube is pulled out from a liquid ordered (L(o)) domain at a constant speed, slow pulling leads to the formation of a liquid disordered (L(d)) tube, while fast pulling results in a L(o) tube. Interestingly, in a membrane tube pulled at an intermediate speed, alternate L(d) and L(o) domains appear in the tube. The sorting dynamics and the corresponding pulling force were systematically studied. The results of this study could lead to a better understanding of the dynamic sorting and traffic of lipids and proteins in living cells.
Collapse
Affiliation(s)
- Hongyuan Jiang
- Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA.
| |
Collapse
|
19
|
Idema T, Storm C. Analytical expressions for the shape of axisymmetric membranes with multiple domains. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2011; 34:67. [PMID: 21751093 DOI: 10.1140/epje/i2011-11067-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/05/2011] [Revised: 06/09/2011] [Accepted: 06/21/2011] [Indexed: 05/31/2023]
Abstract
Based on the Canham-Helfrich free energy, we derive analytical expressions for the shapes of axisymmetric membranes consisting of multiple domains. We give explicit equations for both closed vesicles and almost cylindrical tubes. Using these expressions, we also find the shape of a tube attached to a spherical vesicle. The resulting shapes compare well to numerical data, and our expressions can be used to easily determine membrane parameters from experimentally obtained shapes.
Collapse
Affiliation(s)
- T Idema
- Department of Physics and Astronomy, University of Pennsylvania, 209 S 33rd street, 19104 Philadelphia, Pennsylvania, USA.
| | | |
Collapse
|
20
|
Callan-Jones A, Sorre B, Bassereau P. Curvature-driven lipid sorting in biomembranes. Cold Spring Harb Perspect Biol 2011; 3:cshperspect.a004648. [PMID: 21421916 DOI: 10.1101/cshperspect.a004648] [Citation(s) in RCA: 97] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
It has often been suggested that the high curvature of transport intermediates in cells may be a sufficient means to segregate different lipid populations based on the relative energy costs of forming bent membranes. In this review, we present in vitro experiments that highlight the essential physics of lipid sorting at thermal equilibrium: It is driven by a trade-off between bending energy, mixing entropy, and interactions between species. We collect evidence that lipid sorting depends strongly on lipid-lipid and protein-lipid interactions, and hence on the underlying composition of the membrane and on the presence of bound proteins.
Collapse
Affiliation(s)
- Andrew Callan-Jones
- CNRS UMR 5207, Laboratoire de Physique Théorique et Astroparticules, Université Montpellier II, 34095 Montpellier Cedex 05, France
| | | | | |
Collapse
|
21
|
Dynamic sorting of lipids and proteins in membrane tubes with a moving phase boundary. Proc Natl Acad Sci U S A 2010; 107:7208-13. [PMID: 20368457 DOI: 10.1073/pnas.0913997107] [Citation(s) in RCA: 106] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Cellular organelle membranes maintain their integrity, global shape, and composition despite vigorous exchange among compartments of lipids and proteins during trafficking and signaling. Organelle homeostasis involves dynamic molecular sorting mechanisms that are far from being understood. In contrast, equilibrium thermodynamics of membrane mixing and sorting, particularly the phase behavior of binary and ternary model membrane mixtures and its coupling to membrane mechanics, is relatively well characterized. Elucidating the continuous turnover of live cell membranes, however, calls for experimental and theoretical membrane models enabling manipulation and investigation of directional mass transport. Here we introduce the phenomenon of curvature-induced domain nucleation and growth in membrane mixtures with fluid phase coexistence. Membrane domains were consistently observed to nucleate precisely at the junction between a strongly curved cylindrical (tube) membrane and a pipette-aspirated giant unilamellar vesicle. This experimental geometry mimics intracellular sorting compartments, because they often show tubular-vesicular membrane regions. Nucleated domains at tube necks were observed to present diffusion barriers to the transport of lipids and proteins. We find that curvature-nucleated domains grow with characteristic parabolic time dependence that is strongly curvature-dependent. We derive an analytical model that reflects the observed growth dynamics. Numerically calculated membrane shapes furthermore allow us to elucidate mechanical details underlying curvature-dependent directed lipid transport. Our observations suggest a novel dynamic membrane sorting principle that may contribute to intracellular protein and lipid sorting and trafficking.
Collapse
|
22
|
Zanghellini J, Wodlei F, von Grünberg HH. Phospholipid demixing and the birth of a lipid droplet. J Theor Biol 2010; 264:952-61. [PMID: 20184900 DOI: 10.1016/j.jtbi.2010.02.025] [Citation(s) in RCA: 54] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2009] [Revised: 01/25/2010] [Accepted: 02/13/2010] [Indexed: 11/26/2022]
Abstract
The biogenesis of lipid droplets (LD) in the yeast Saccharomyces cerevisiae was theoretically investigated on basis of a biophysical model. In accordance with the prevailing model of LD formation, we assumed that neutral lipids oil-out between the membrane leaflets of the endoplasmic reticulum (ER), resulting in LD that bud-off when a critical size is reached. Mathematically, LD were modeled as spherical protuberances in an otherwise planar ER membrane. We estimated the local phospholipid composition, and calculated the change in elastic free energy of the membrane caused by nascent LD. Based on this model calculation, we found a gradual demixing of lipids in the membrane leaflet that goes along with an increase in surface curvature at the site of LD formation. During demixing, the phospholipid monolayer was able to gain energy during LD growth, which suggested that the formation of curved interfaces was supported by or even driven by lipid demixing. In addition, we show that demixing is thermodynamically necessary as LD cannot bud-off otherwise. In the case of Saccharomyces cerevisiae our model predicts a LD bud-off diameter of about 12 nm. This diameter is far below the experimentally determined size of typical yeast LD. Thus, we concluded that if the standard model of LD formation is valid, LD biogenesis is a two step process. Small LD are produced from the ER, which subsequently ripe within the cytosol through a series of fusions.
Collapse
Affiliation(s)
- J Zanghellini
- School of Mathematics and Statistics, University of South Australia, Mawson Lakes, SA 5095, Australia.
| | | | | |
Collapse
|
23
|
Tian A, Capraro BR, Esposito C, Baumgart T. Bending stiffness depends on curvature of ternary lipid mixture tubular membranes. Biophys J 2009; 97:1636-46. [PMID: 19751668 DOI: 10.1016/j.bpj.2009.07.012] [Citation(s) in RCA: 64] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2009] [Revised: 07/01/2009] [Accepted: 07/13/2009] [Indexed: 12/14/2022] Open
Abstract
Lipid and protein sorting and trafficking in intracellular pathways maintain cellular function and contribute to organelle homeostasis. Biophysical aspects of membrane shape coupled to sorting have recently received increasing attention. Here we determine membrane tube bending stiffness through measurements of tube radii, and demonstrate that the stiffness of ternary lipid mixtures depends on membrane curvature for a large range of lipid compositions. This observation indicates amplification by curvature of cooperative lipid demixing. We show that curvature-induced demixing increases upon approaching the critical region of a ternary lipid mixture, with qualitative differences along two roughly orthogonal compositional trajectories. Adapting a thermodynamic theory earlier developed by M. Kozlov, we derive an expression that shows the renormalized bending stiffness of an amphiphile mixture membrane tube in contact with a flat reservoir to be a quadratic function of curvature. In this analytical model, the degree of sorting is determined by the ratio of two thermodynamic derivatives. These derivatives are individually interpreted as a driving force and a resistance to curvature sorting. We experimentally show this ratio to vary with composition, and compare the model to sorting by spontaneous curvature. Our results are likely to be relevant to the molecular sorting of membrane components in vivo.
Collapse
Affiliation(s)
- Aiwei Tian
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| | | | | | | |
Collapse
|
24
|
Curvature-driven lipid sorting needs proximity to a demixing point and is aided by proteins. Proc Natl Acad Sci U S A 2009; 106:5622-6. [PMID: 19304798 DOI: 10.1073/pnas.0811243106] [Citation(s) in RCA: 297] [Impact Index Per Article: 19.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
Sorting of lipids and proteins is a key process allowing eukaryotic cells to execute efficient and accurate intracellular transport and to maintain membrane homeostasis. It occurs during the formation of highly curved transport intermediates that shuttle between cell compartments. Protein sorting is reasonably well described, but lipid sorting is much less understood. Lipid sorting has been proposed to be mediated by a physical mechanism based on the coupling between membrane composition and high curvature of the transport intermediates. To test this hypothesis, we have performed a combination of fluorescence and force measurements on membrane tubes of controlled diameters pulled from giant unilamellar vesicles. A model based on membrane elasticity and nonideal solution theory has also been developed to explain our results. We quantitatively show, using 2 independent approaches, that a difference in lipid composition can build up between a curved and a noncurved membrane. Importantly, and consistent with our theory, lipid sorting occurs only if the system is close to a demixing point. Remarkably, this process is amplified when even a low fraction of lipids is clustered upon cholera toxin binding. This can be explained by the reduction of the entropic penalty of lipid sorting when some lipids are bound together by the toxin. Our results show that curvature-induced lipid sorting results from the collective behavior of lipids and is even amplified in the presence of lipid-clustering proteins. In addition, they suggest a generic mechanism by which proteins can facilitate lipid segregation in vivo.
Collapse
|