<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Mol. Biosci.</journal-id>
<journal-title>Frontiers in Molecular Biosciences</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Mol. Biosci.</abbrev-journal-title>
<issn pub-type="epub">2296-889X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">747601</article-id>
<article-id pub-id-type="doi">10.3389/fmolb.2021.747601</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Molecular Biosciences</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Osh6 Revisited: Control of PS Transport by the Concerted Actions of PI4P and Sac1 Phosphatase</article-title>
<alt-title alt-title-type="left-running-head">Eisenreichova et&#x20;al.</alt-title>
<alt-title alt-title-type="right-running-head">Lipid Transport Mediated by Osh6</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Eisenreichova</surname>
<given-names>Andrea</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1469792/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>R&#xf3;&#x17c;ycki</surname>
<given-names>Bartosz</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1081129/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Boura</surname>
<given-names>Evzen</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Humpolickova</surname>
<given-names>Jana</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1164108/overview"/>
</contrib>
</contrib-group>
<aff id="aff1">
<label>
<sup>1</sup>
</label>Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, <addr-line>Prague</addr-line>, <country>Czechia</country>
</aff>
<aff id="aff2">
<label>
<sup>2</sup>
</label>Institute of Physics, Polish Academy of Sciences, <addr-line>Warsaw</addr-line>, <country>Poland</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/832811/overview">Aidan D. Meade</ext-link>, Technological University Dublin, Ireland</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1442402/overview">Agata Witkowska</ext-link>, Leibniz-Institut f&#xfc;r Molekulare Pharmakologie (FMP), Germany</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1204368/overview">Guang-Kui Xu</ext-link>, Xi&#x2019;an Jiaotong University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Jana Humpolickova, <email>jana.humpolickova@uochb.cas.cz</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Biophysics, a section of the journal Frontiers in Molecular Biosciences</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>12</day>
<month>10</month>
<year>2021</year>
</pub-date>
<pub-date pub-type="collection">
<year>2021</year>
</pub-date>
<volume>8</volume>
<elocation-id>747601</elocation-id>
<history>
<date date-type="received">
<day>26</day>
<month>07</month>
<year>2021</year>
</date>
<date date-type="accepted">
<day>08</day>
<month>09</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2021 Eisenreichova, R&#xf3;&#x17c;ycki, Boura and Humpolickova.</copyright-statement>
<copyright-year>2021</copyright-year>
<copyright-holder>Eisenreichova, R&#xf3;&#x17c;ycki, Boura and Humpolickova</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these&#x20;terms.</p>
</license>
</permissions>
<abstract>
<p>Osh6, a member of the oxysterol-binding protein&#x2013;related protein (ORP) family, is a lipid transport protein that is involved in the transport of phosphatidylserine (PS) between the endoplasmic reticulum (ER) and the plasma membrane (PM). We used a biophysical approach to characterize its transport mechanism in detail. We examined the transport of all potential ligands of Osh6. PI4P and PS are the best described lipid cargo molecules; in addition, we showed that PIP2 can be transported by Osh6 as well. So far, it was the exchange between the two cargo molecules, PS and PI4P, in the lipid-binding pocket of Osh6 that was considered an essential driving force for the PS transport. However, we showed that Osh6 can efficiently transport PS along the gradient without the help of PI4P and that PI4P inhibits the PS transport along its gradient. This observation highlights that the exchange between PS and PI4P is indeed crucial, but PI4P bound to the protein rather than intensifying the PS transport suppresses it. We considered this to be important for the transport directionality as it prevents PS from returning back from the PM where its concentration is high to the ER where it is synthesized. Our results also highlighted the importance of the ER resident Sac1 phosphatase that enables the PS transport and ensures its directionality by PI4P consumption. Furthermore, we showed that the Sac1 activity is regulated by the negative charge of the membrane that can be provided by PS or PI anions in the case of the ER membrane.</p>
</abstract>
<kwd-group>
<kwd>oxysterol-binding protein&#x2013;related proteins</kwd>
<kwd>OSH6</kwd>
<kwd>phosphatidylinositol 4-phosphate</kwd>
<kwd>phosphatidylserine</kwd>
<kwd>phosphatidylinositol (4, 5)-bisphosphate</kwd>
<kwd>lipid transport</kwd>
<kwd>SAC1</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>In the past few years, the emerging interest in lipid transport proteins shed some light on the biology of the oxysterol-binding protein (OSBP)-related protein (ORP) family (<xref ref-type="bibr" rid="B28">Olkkonen, 2013</xref>; <xref ref-type="bibr" rid="B27">Olkkonen and Li, 2013</xref>). ORPs have been identified as lipid transporters that are responsible for transferring specific lipids such as phosphatidylserine (PS) or cholesterol/ergosterol from the site of their synthesis, usually the ER, to their respective acceptor membrane. The mechanism which allows such a selective lipid transport across the concentration gradient relies on phosphoinositides (PPIns), mostly phosphatidylinositol 4-phosphate (PI4P) (<xref ref-type="bibr" rid="B35">von Filseck et&#x20;al., 2015b</xref>).</p>
<p>It has been postulated that the main feature of the transport mechanism is the exchange of the cargo lipid molecule at the target membrane for the PPIn molecule which, upon delivery to the ER, is hydrolyzed, providing the PPIn gradient as a driving force (energy) for the transport against the concentration gradient. This has been first shown elegantly for the yeast ergosterol transporter Osh4 (<xref ref-type="bibr" rid="B7">de Saint-Jean et&#x20;al., 2011</xref>; <xref ref-type="bibr" rid="B35">von Filseck et&#x20;al., 2015b</xref>) and later for the yeast PS transporters Osh6 and Osh7 (<xref ref-type="bibr" rid="B24">Maeda et&#x20;al., 2013</xref>; <xref ref-type="bibr" rid="B34">von Filseck et&#x20;al., 2015a</xref>), as well as for the human cholesterol transporter OSBP (<xref ref-type="bibr" rid="B26">Mesmin et&#x20;al., 2017</xref>).</p>
<p>To maintain the phosphoinositide gradient in live cells, three processes are required: i) PI transport to the PM, ii) its conversion to PI4P by phosphatidylinositol 4-kinases, and iii) PI4P hydrolysis in the ER membrane back to PI carried out by the Sac1 phosphatase (<xref ref-type="bibr" rid="B3">Cai et&#x20;al., 2014</xref>; <xref ref-type="bibr" rid="B5">Chung et&#x20;al., 2015</xref>) that has an essential role in the transport process (<xref ref-type="bibr" rid="B35">von Filseck et&#x20;al., 2015b</xref>).</p>
<p>In this study, we focused on the lipid-binding features of the best characterized PS transporter, the yeast protein Osh6. We compared its binding preferences (relative affinities) toward PS, PI4P, and also to phosphatidylinositol (4,5)-bisphosphate (PIP2) as their relative affinities are critical for defining the directionality of transport. We clearly demonstrated that PI4P is a significantly stronger binder than PS. This protects PS from being transported along the gradient from the PM back to ER. This affinity difference is a defining feature as it keeps the PS plasma membrane level high compared to the level in the endoplasmic reticulum (ER) and prevents unproductive lipid shuffling between the two membranes. We also showed that Sac1 is indispensable for preventing PI4P rebinding to the Osh6 binding pocket and demonstrated how Sac1 allosteric regulation by PS (<xref ref-type="bibr" rid="B42">Zhong et&#x20;al., 2012</xref>) contributes to the transport regulation.</p>
</sec>
<sec sec-type="materials|methods" id="s2">
<title>Materials and Methods</title>
<sec id="s2-1">
<title>Protein Expression and Purification</title>
<p>All expression vectors were constructed using restriction cloning. The genes for C2<sub>Lact</sub>-CFP, C2<sub>granuphilin</sub>-CFP, and SidC-CFP were cloned into a modified pHIS-2 vector with an N-terminal 6x-His tag followed by a tobacco etch virus (TEV) protease cleavage site and a CFP coding sequence. The gene for Osh6 was cloned into a pRSFD (Novagen) vector with an N-terminal His6x-GB1 solubility tag and a TEV cleavage site. Mutation to produce <ext-link ext-link-type="uri" xlink:href="https://www.google.com/search?sxsrf=ALeKk00hbjn30hJTvBSlThFDPqckaQNP2Q:1622318221644&amp;q=catalytically&amp;spell=1&amp;sa=X&amp;ved=2ahUKEwiqzsCq1u_wAhWENOwKHfyOAO0QkeECKAB6BAgBEDE">catalytically</ext-link> dead Sac1(C392S) was generated by site-directed mutagenesis. The proteins, including the C-terminal His-tagged Sac1 and Sac1(C392S), were expressed in <ext-link ext-link-type="uri" xlink:href="https://www.elabprotocols.com/protocols/">ZY</ext-link>P-5052 autoinduction media in <italic>Escherichia coli</italic> BL21 Star cells using our standard protocols (<xref ref-type="bibr" rid="B4">Chalupska et&#x20;al., 2017</xref>; <xref ref-type="bibr" rid="B19">Klima et&#x20;al., 2017</xref>). Briefly, the cells were harvested by centrifugation, resuspended in a lysis buffer (50&#xa0;mM Tris pH 8, 300&#xa0;mM NaCl, 20&#xa0;mM imidazol, and 3&#xa0;mM &#x3b2;-mercaptoethanol), and lysed by sonication. After the lysate was cleared by centrifugation and incubated with a nickel-chelating resin (Machery-Nagel), the proteins were eluted with the lysis buffer supplemented with 300&#xa0;mM imidazol. The hexahistidine tag or GB1 tag of SidC-CFP, C2<sub>Lact</sub>-CFP, C2<sub>granuphilin</sub>-CFP, and Osh6 was cleaved off with TEV protease. When required, the proteins were purified by cation exchange (Osh6) and anion exchange (Sac1 and Sac1 (C392S)). The final purification step of all proteins comprised exclusion chromatography on <ext-link ext-link-type="uri" xlink:href="https://www.sigmaaldrich.com/catalog/product/sigma/ge28989333?lang=en&amp;region=US">HiLoad 16/600 Superdex 75 pg (Cytiva</ext-link>) in 20mM Tris pH &#x3d; 7.4, 300&#xa0;mM NaCl, and 3mM &#x3b2;-mercaptoethanol. Purified recombinant proteins were stored at &#x2212;80&#xb0;C.</p>
</sec>
<sec id="s2-2">
<title>Preparation of LUVs</title>
<p>Large unilamellar vesicles (LUVs) were prepared as before (<xref ref-type="bibr" rid="B12">Humpolickova et&#x20;al., 2017</xref>). Briefly, stock solutions of lipids and labeled lipid analogs in organic solvents (chloroform, methanol) were mixed in the required ratio. The organic solvents were evaporated in the stream of dry nitrogen, and the lipid films were placed into the vacuum chamber for an hour to remove all the traces of the original solvents. The films were than resuspended in the LUV buffer (10&#xa0;mM TRIS, pH &#x3d; 7.4, 150&#xa0;mM NaCl, 3&#xa0;mM &#x3b2;-mercaptoethanol, and 1&#xa0;mM EDTA). The solution of multilamellar membranes was homogenized in an extruder (Avestin, Ottawa, Canada) through the membrane with 50 or 100&#xa0;nm pores. The total concentration of lipids in the LUV buffer was 1&#xa0;mM.</p>
</sec>
<sec id="s2-3">
<title>FCCS Microscopy</title>
<p>The FCCS experiments were acquired on a confocal microscope LSM 780 (Zeiss, Hamburg, Germany). The microscope was equipped with two external SPAD detectors and a FLIM upgrade kit (Hydraharp, Picoquant, Berlin, Germany), which was used for data collecting. The excitation light was guided to the water immersion objective (&#xd7; 40, N.A. &#x3d; 1.2), and the signal produced after it had been spatially filtered by a pinhole which was split on the two SPAD detectors with 482/35 and 679/41 emission filters placed in front of them. We used 458&#xa0;nm argon laser line (continuous wave&#x2014;cw) to excite CFP and 625&#xa0;nm line of the In tune laser (pulsed, 40&#xa0;MHz) to excite DiD. The data were correlated using the homemade scripts in MATLAB (Mathworks, Nattick, MA) (<xref ref-type="bibr" rid="B36">Wahl et&#x20;al., 2003</xref>), and the different TCSPC patterns resulting from the combination of the cw and pulsed excitation were used for eliminating the bleed-through signal (<xref ref-type="bibr" rid="B32">Skerle et&#x20;al., 2020</xref>).</p>
</sec>
<sec id="s2-4">
<title>FCCS Assays</title>
<p>In both PI4P and PIP2 transport assays, 10&#xa0;&#xb5;l of PI4P/PIP2-containing LUVs were mixed with 40&#xa0;&#xb5;l of LUVs that served as the target for PI4P/PIP2 and that also contained the lipid analog DiD (DiD/lipid &#x3d; 1/10,000). A biosensor (SidC-CFP and C2<sub>granuphilin</sub>-CFP) was added so that its concentration in the 200&#xa0;&#xb5;l total volume was 100&#xa0;nM. Short (200&#xa0;s) FCCS experiment was acquired prior to the addition of Osh6 (1&#xa0;&#xb5;M final concentration). Immediately upon Osh6 addition, the measurement was set for additional 30&#xa0;min. For the transport of either PI4P or PIP2, the LUVs contained 5&#xa0;mol% of the particular phosphoinositide; for the experiments following the transport of either PI4P or PIP2 in the mixture with the other phosphoinositide, the LUVs contained 5&#xa0;mol% of each. The PS content in the acceptor LUVs scales from 0 to 20&#xa0;mol%.</p>
<p>To address the effect of the Sac1 participation, the PI4P-accepting membrane contained 5&#xa0;mol% of DGS-NTA(Ni), which assures the Sac1 attachment <italic>via</italic> its C-terminal His-tag. The concentration of either wild-type Sac1 or C392S mutant in the experiment was 1&#xa0;&#xb5;M.</p>
<p>In the assay that addresses the Sac1 phosphatase activity, LUVs contained 0.1&#xa0;mol% of PI4P (to keep the changes in the dynamic range of the biosensor), 5&#xa0;mol% of DGS-NTA(Ni), and the 20&#xa0;mol% of the negatively charged lipid (PS, PI, and PG), and the lipid analog DiD (DiD/lipid &#x3d; 1/10,000). 40&#x20;&#xb5;l of LUVs were mixed with SidC-CFP (final concentration: 100&#xa0;nM), Sac1 (final concentration: 1&#xa0;&#xb5;M), and the LUV buffer to final volume of 200&#xa0;&#xb5;l. Upon Sac1 addition, the measurement was run for 30&#xa0;min.</p>
<p>All the FCCS experiments were acquired at least in three independent experiments to make sure that the observed trends are reproducible and robust.</p>
</sec>
<sec id="s2-5">
<title>FRET Assay</title>
<p>In a PS transport assay, 100&#xa0;&#x3bc;l of donor LUVs containing 90&#xa0;mol% POPC, 5&#xa0;mol% POPS, and 5&#xa0;mol% DGS-NTA(Ni) were mixed with 100&#xa0;&#x3bc;l of acceptor LUVs. The acceptor LUVs consisted of 95% POPC and 5&#xa0;mol% ATTO488-DOPE. To monitor the effect of phosphoinositides on the PS transport, the acceptor LUVs contained additional 5&#xa0;mol% PI4P or PIP2. A fluorescence biosensor for PS, C2<sub>Lact</sub>-CFP, was added to the final concentration of 125&#xa0;nM (in the final volume of 800&#xa0;&#x3bc;l). Shortly after the start of the measurement, Osh6 (100&#xa0;nM) was injected. When Sac1 or Sac1(C392S) (200&#xa0;nM) was present in the transport, it was first incubated for 5&#xa0;min with the donor LUVs to allow its binding to DGS-NTA(Ni). PS transport was monitored by measuring the FRET between C2<sub>Lact</sub>-CFP and ATTO488-DOPE at 527&#xa0;nm upon excitation at 400&#xa0;nm. The measurement was run for 30&#xa0;min on a Fluoromax 4 spectrophotometer (Horiba, Kyoto, Japan). The FRET kinetics was measured in at least three independent sets to prove the reproducibility and robustness of the results.</p>
</sec>
<sec id="s2-6">
<title>Choice of the Assay</title>
<p>In our study, we used different assays for different cargos. In the case of PS transport, we have realized that the C2<sub>Lact</sub>-CFP response depends on the presence of PPIns in the PS-accepting membrane. The biosensor was less sensitive to PS in PPIn&#x2019;s absence. This significantly complicated the FCCS approach as it was almost blind to small amounts of PS in the accepting membrane devoid of PPIns. The FRET approach was more robust, and its readout was insensitive to PPIns. On the contrary, using the FRET approach for PPIns that are less transported than PS gives small overall readout with significant noise. Altogether, we paid attention to i) whether the readout monitors the level of the evaluated cargo independently on other lipids significant for the transport process, and ii) to the extent and quality of the readout change.</p>
</sec>
<sec id="s2-7">
<title>GUV Preparation and Imaging</title>
<p>Giant unilamellar vesicles (GUVs) were prepared as before (<xref ref-type="bibr" rid="B2">Boura and Hurley, 2012</xref>; <xref ref-type="bibr" rid="B8">Dubankova et&#x20;al., 2017</xref>). Lipids were mixed in the desired ratio in chloroform so that their final concentration was 5&#xa0;&#xb5;g/&#xb5;l. 2&#x20;&#xd7; 9&#xa0;&#xb5;l was spread on two ITO-coated glass electrodes that were plane parallely assembled into a Teflon chamber. Five ml of 600&#xa0;mM sucrose was added, and a sinusoidal voltage of 10&#xa0;kHz frequency and an amplitude of 1&#xa0;V were applied on them for 1&#xa0;hour at 60&#xb0;C. After that, 50&#xa0;&#xb5;l of donor and 50&#xa0;&#xb5;l of acceptor GUVs were mixed with 100&#xa0;&#xb5;l of isosmotic buffer (25&#xa0;mM Tris pH 8, 10&#xa0;mM MgCl<sub>2</sub>, 20&#xa0;mM imidazole, 261.5&#xa0;mM NaCl, and 2&#xa0;mM &#x3b2;ME) containing the appropriate biosensor (70&#xa0;nM) and Osh6 (100&#xa0;nM). The mixture was placed on the BSA-coated cover glass (bottom of the 4-chamber dish), and the images were&#x20;taken.</p>
<p>The images were acquired on a Leica SP8 confocal microscope (Mannheim, Germany) using &#xd7;60 water immersion objective lens. To avoid bleed-through, the images were captured in line sequentially. CFP was excited by 458&#xa0;nm argon laser line, and DiD that was added to the acceptor GUVs was excited by 633&#xa0;nm He&#x2013;Ne laser line. Emission light was split onto two HyD detectors according to the emission spectra of the fluorophores.</p>
</sec>
<sec id="s2-8">
<title>Chemicals</title>
<p>All the lipids in this study were purchased from Avanti Polar Lipids (Alabaster, AL). In particular, we used POPC, POPS, brain PI4P, brain PI(4,5)P2, liver PI, POPG, and DGS-NTA(Ni). ATTO488-DOPE was purchased from ATTO-TEC (Siegen, Germany). DiD and other chemicals were purchased from Sigma-Aldrich (St. Louis,&#x20;MO).</p>
</sec>
<sec id="s2-9">
<title>Molecular Dynamics Simulations</title>
<p>The atomic coordinates of yeast Sac1 were taken from the crystal structure deposited in the Protein Data Bank (PDB) with the entry code of 3LWT (<xref ref-type="bibr" rid="B25">Manford et&#x20;al., 2010</xref>). Monte Carlo (MC) simulations of Sac1 in contact with a planar membrane were performed within the framework of the Kim&#x2013;Hummer model (<xref ref-type="bibr" rid="B16">Kim and Hummer, 2008</xref>), in which both the membrane and solvent are implicit, whereas crystallized protein domains are treated as rigid bodies. Here, Sac1 was taken as a single rigid body as determined by the crystal structure with the PDB code 3LWT. In the course of the MC simulations, the Sac1 rigid structure was subjected to translations and rotations relative to the membrane surface. Simulated annealing was used to determine the minimum of Sac1&#x2013;membrane interaction energy. The position of the Sac1 rigid structure corresponding to the energy minimum was used as the initial position of Sac1 in MD simulations.</p>
<p>The initial systems for MD simulations were prepared using the Input Generator and the Membrane Builder in the CHARMM-GUI website (<xref ref-type="bibr" rid="B14">Jo et&#x20;al., 2007</xref>; <xref ref-type="bibr" rid="B15">Jo et&#x20;al., 2008</xref>; <xref ref-type="bibr" rid="B40">Wu et&#x20;al., 2014</xref>; <xref ref-type="bibr" rid="B20">Lee et&#x20;al., 2016</xref>). Specifically, two bilayer segments with equal lateral dimensions of 11 by 11&#xa0;nm were formed independently. In one case, the bilayer was composed of 292 POPC and 72 POPS lipids (i.e.,&#x20;with the molar ratio of about 4:1). In the other case, the bilayer was composed of 352 POPC lipids. In both cases, Sac1 was placed on top of the lipid bilayer in accordance with the minimum energy position obtained from the MC simulations. The lipid&#x2013;protein systems were solvated next, and sodium and chloride ions were added to neutralize the systems and to reach a physiological ion concentration of 150&#xa0;mM.</p>
<p>The MD simulations were performed using NAMD 2.10 with CHARMM36 force field and the TIP3P model for water molecules (<xref ref-type="bibr" rid="B23">MacKerell et&#x20;al., 1998</xref>; <xref ref-type="bibr" rid="B17">Klauda et&#x20;al., 2010</xref>; <xref ref-type="bibr" rid="B1">Best et&#x20;al., 2012</xref>; <xref ref-type="bibr" rid="B18">Klauda et&#x20;al., 2012</xref>; <xref ref-type="bibr" rid="B29">Phillips et&#x20;al., 2020</xref>). Temperature was kept at 303&#xa0;K through a Langevin thermostat with a damping coefficient of 1/ps. Pressure was maintained at 1&#xa0;atm using the Langevin piston Nose&#x2013;Hoover method with a damping timescale of 25 fs and an oscillation period of 50 fs. Short-range nonbonded interactions were cutoff smoothly between one and 1.2&#xa0;nm. Long-range electrostatic interactions were computed using the particle mesh Ewald method with a grid spacing of 0.1&#xa0;nm. Simulations were performed with an integration time step of 2&#x20;fs.</p>
<p>The initial systems for MD simulations were energy-minimized using a conjugate gradient method and then equilibrated in a standard procedure using input files provided by the CHARMM-GUI Membrane Builder (<xref ref-type="bibr" rid="B14">Jo et&#x20;al., 2007</xref>; <xref ref-type="bibr" rid="B40">Wu et&#x20;al., 2014</xref>). Next, for each of the two simulation systems (i.e.,&#x20;POPC-POPS-Sac1 and POPC-Sac1), we performed two production runs of 60&#xa0;ns each, amounting to 240&#xa0;ns of MD data for analysis. The trajectories were visualized and analyzed using VMD (<xref ref-type="bibr" rid="B11">Humphrey et&#x20;al., 1996</xref>). To determine contacts between amino acid residues and lipid molecules, we used a simple distance criterion, namely, an amino acid residue was taken to be in contact with a lipid molecule if the minimum distance between non-hydrogen atoms of the residue and the lipid was smaller than 0.45&#xa0;nm.</p>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>Results</title>
<p>Lipid transport can be observed as pertinent to each of the potential participants, that is, all lipids that can be potentially transported by Osh6. We have examined the transport of both PS and PI4P that have been identified to bind Osh6. In addition, we also tested PIP2 as a potential cargo of the Osh6 molecule.</p>
<sec id="s3-1">
<title>PS Transport</title>
<p>We examined the transport of PS first because Osh6 was reported to be a PS transporter in several studies (<xref ref-type="bibr" rid="B24">Maeda et&#x20;al., 2013</xref>; <xref ref-type="bibr" rid="B34">von Filseck et&#x20;al., 2015a</xref>). We have utilized a FRET assay to monitor the amount of transported PS in the acceptor membrane by observing the energy transfer between CFP fused to a PS biosensor [C2 domain from Lactadherin C (LactC) ] and ATTO488-DOPE in the acceptor membrane (<xref ref-type="fig" rid="F1">Figures 1A,B</xref>). To validate the utility of this assay, we determined the concentration range of PS at which the assay provides a dynamic response to PS change. These preliminary experiments revealed that the FRET assay can be used between 0 and 2&#xa0;mol% of PS in the membrane of large unilamellar vesicles (LUVs). PS increases above this level do not cause a further increase in the FRET signal (<xref ref-type="sec" rid="s10">Supplementary Figure&#x20;S1</xref>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>PS transport. <bold>(A)</bold> and <bold>(B)</bold> Schematics of the FRET assay with and without Sac1. Temporal evolution of the PS transport to LUVs composed of POPC solely <bold>(C)</bold>, to LUVs with 5&#xa0;mol% of PI4P <bold>(D)</bold>, to LUVs with 5&#xa0;mol% of PIP2&#x20;<bold>(E)</bold>. LUVs are either not decorated by Sac1 (solid curves), or are decorated by Sac1 wt (dotted curves), or by Sac1 C392S (dashed curves). <bold>(C&#x2013;E)</bold> Representative curves out of at least three independent experiments. <bold>(F)</bold> Reconstitution of the PS transport to the membranes of GUVs. Red GUVs are PS-accepting GUVs, and presence of the green PS biosensor on the red GUVs refers to transport.</p>
</caption>
<graphic xlink:href="fmolb-08-747601-g001.tif"/>
</fig>
<p>Having a well-characterized assay in hands, we examined the PS transport to various acceptor membranes differing in the lipid composition. <xref ref-type="fig" rid="F1">Figure&#x20;1C</xref> clearly shows that PS transport along a concentration gradient to an acceptor membrane composed only of POPC (and the fluorescence acceptor molecules, ATTO488-DOPE) is readily detectable. Simply viewed, the Osh6 molecules pick up PS molecules from any membrane of high PS content and deliver them to the membrane of low PS content, which is a process simply driven by entropy. If this were to happen in intact cells, it would eliminate the high PS concentration in the PM. We reasoned that PI4P or PIP2, both present in the PM, would inhibit the undesirable backward PS transport between the PM and the ER. To test this, we applied the same amount of PI4P in the PS-accepting membrane and PS in the donating membrane and found that the PS transport was significantly impaired (<xref ref-type="fig" rid="F1">Figure&#x20;1D</xref>). Moreover, PIP2 can also bind to the Osh6 transport pocket, but, in contrast to PI4P, the presence of PIP2 in the acceptor membrane inhibits PS transport only slightly (<xref ref-type="fig" rid="F1">Figure&#x20;1E</xref>).</p>
<p>We then focused on the role of the Sac1 phosphatase. We tethered Sac1 to the donor PS membrane (<xref ref-type="fig" rid="F1">Figure&#x20;1B</xref>) by His-tag on the C-terminus of the cytosolic part of Sac1 that binds to DGS-NTA(Ni) present in the PS-donating membrane. Upon Sac1 attachment, the PS transported to the PI4P-containing membranes can be fully restored (<xref ref-type="fig" rid="F1">Figure&#x20;1D</xref> dotted curve). The specificity of this process was demonstrated as the catalytically dead mutant of Sac1(C392S) was without effect (<xref ref-type="fig" rid="F1">Figure&#x20;1D</xref> dashed curve). This suggests that the presence of PI4P in the acceptor membrane can keep the PS gradient between the plasma membrane and ER. The hydrolysis of PI4P in the ER membrane by Sac1 shifts the PI4P-binding equilibrium and leads to the release of the transported PI4P molecule and freeing up the lipid-binding pocket of Osh6 for further binding of PS to deliver it to the PM. By the cooperation of PI4P and Sac1, transportation of PS can be effectively regulated. <xref ref-type="fig" rid="F1">Figures 1C,E</xref> show a small drop in the PS transport efficiencies upon decoration of the PS-containing membrane with Sac1 (wt as well as C392S mutant). This may be a result of protein crowding on the membrane, that is, lower access of Osh6 to the&#x20;cargo.</p>
<p>Also, we reconstituted the PS transport into the membrane of giant unilamellar vesicles (GUVs) (<xref ref-type="fig" rid="F1">Figure&#x20;1F</xref>). The PS-accepting GUVs were doped with DiD (red circles), and the PS-donating membranes were labeled by the PS biosensor C2-CFP (green circles). The transport was observed approximately 30&#xa0;min upon Osh6 addition. The appearance of green color on the originally red membranes refers on the transport process in progress.</p>
</sec>
<sec id="s3-2">
<title>PI4P Transport</title>
<p>PI4P transport was monitored by fluorescence cross-correlation spectroscopy (FCCS). Donor LUVs containing PI4P were mixed with PI4P-accepting LUVs of varied lipid composition also containing the fluorescence marker DiD. The PI4P biosensor, SidC-CFP, was localized only at the donor membrane at time zero. Upon Osh6 addition, the biosensor started to move to the acceptor LUVs following the PI4P transport (<xref ref-type="fig" rid="F2">Figure&#x20;2A</xref>). As a result of this transport process, we could observe the cross-correlation of the fluorescence signal of CFP and DiD. Importantly, also for this assay, we established its dynamic range which was found to be between 0 and 1&#xa0;mol% (<xref ref-type="sec" rid="s10">Supplementary Figure S1</xref>) in agreement with the nanomolar-binding affinity of SidC and PI4P (<xref ref-type="bibr" rid="B22">Luo et&#x20;al., 2015</xref>).</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>PI4P transport observed by the FCCS assay. <bold>(A)</bold> Schematics of the FCCS assay. <bold>(B)</bold> Temporal evolution of the transport-monitoring parameter <italic>G</italic>
<sub>cc</sub>/<italic>G</italic>
<sub>R</sub>. PI4P was transported to the LUVs with an increasing amount of PS: POPC only (full squares), 10&#xa0;mol% of PS (half-filled squares), and 20&#xa0;mol% of PS (empty squares). The PI4P-donating membrane (5&#xa0;mol% PI4P) was also enriched in PIP2 (5&#xa0;mol%): wo/w PIP2&#x2014;black/red squares. The error bars stand for the standard deviation. <bold>(C)</bold> Reconstitution of the PI4P transport to the membranes of GUVs. Red GUVs are PI4P-accepting GUVs, and the presence of the green PI4P biosensor on the red GUVs refers on the transport. The presented data represent the trends that were observed in minimum three independent experiments.</p>
</caption>
<graphic xlink:href="fmolb-08-747601-g002.tif"/>
</fig>
<p>We already established that PI4P can efficiently block the PS transport along the gradient, that is, Osh6 significantly favors PI4P binding in comparison to PS. Next, we focused on the PI4P transport and its dependence on the PS concentration in the PI4P-accepting membrane.</p>
<p>Our FCCS analysis revealed that PI4P is also transported along its gradient by Osh6 although the speed of PI4P transport was slower than that of PS (<xref ref-type="fig" rid="F2">Figure&#x20;2B</xref>). In the case of the PS transport (<xref ref-type="fig" rid="F1">Figure&#x20;1</xref>), saturation of the FRET response happens within few minutes, whereas in the case of SidC that senses PI4P, only small amount of PI4P with a continuous increase appears in the target membrane even after 30&#xa0;min upon Osh6 addition.</p>
<p>PIP2 is the major phosphoinositide of the plasma membrane (<xref ref-type="bibr" rid="B33">van Meer et&#x20;al., 2008</xref>); we have enriched the PI4P-containing LUVs by PIP2 in the 1&#x2013;1&#xa0;M ratio (5&#xa0;mol% each). However, this alteration did not have any effect on the PI4P transport, suggesting lower affinity of PIP2 to Osh6 than PI4P (<xref ref-type="fig" rid="F2">Figure&#x20;2B</xref> red curves).</p>
<p>PI4P transport was also reconstructed in the membranes of GUVs (<xref ref-type="fig" rid="F2">Figure&#x20;2C</xref>), PI4P donor GUVs were labeled by the SidC-CFP, and PI4P acceptor GUVs contained the fluorescence marker DiD. The appearance of the double-labeled vesicles visualized approximately 30&#xa0;min upon Osh6 addition reports on the transport. The visual inspection of the GUVs shows a higher CFP signal in the acceptor GUVs that contained 20&#xa0;mol% of&#x20;PS.</p>
<p>We also used the FCCS assay to explore the impact of Sac1 on the PI4P transport. In this setup, Sac1 was again tethered to the PI4P-accepting membrane. We expected to see the result of two tightly linked processes, that is, the transfer of PI4P followed by its hydrolysis in the acceptor membrane. Because FCCS counts the double-labeled LUVs, PI4P accumulation in the acceptor vesicles as a result of its transfer from donor membranes would appear as an increase in the cross-correlation amplitude, whereas PI4P hydrolysis would result in a drop of the cross-correlation amplitude. Altogether, the transfer coupled to consequent hydrolysis shows a smaller increase of the cross-correlation amplitude than the situation when the hydrolysis does not&#x20;occur.</p>
<p>Remarkably, the effect of Sac1 showed a dependence on the lipid composition of the PI4P acceptor membrane. In a bilayer composed of POPC only, we did not observe any effect of Sac1 (<xref ref-type="fig" rid="F3">Figure&#x20;3A</xref> dark gray curve). However, a negatively charged membrane containing 20&#xa0;mol% PS supported much enhanced Sac1-mediated PI4P hydrolysis (<xref ref-type="fig" rid="F3">Figure&#x20;3A</xref> orange curve). In the case of the Sac1 catalytically dead mutant, the amount of transported PI4P was almost identical to that of the control without Sac1 (<xref ref-type="fig" rid="F3">Figure&#x20;3A</xref> gray curve) for PC and dark yellow curve for 20&#xa0;mol% of PS. These results suggest that the presence of PS in the acceptor membrane changes the activity of Sac1 and consequently the elimination of&#x20;PI4P.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>PI4P transport followed by Sac1 dephosphorylation viewed by FCCS. <bold>(A)</bold> The temporal evolution of the transport-monitoring parameter <italic>G</italic>
<sub>cc</sub>/<italic>G</italic>
<sub>R</sub> and <bold>(B)</bold> the temporal evolution of the amplitude of the CFP autocorrelation function, and the parameter reporting on the biosensor binding <italic>G</italic>
<sub>G</sub>. The PI4P-accepting membrane was composed of POPC and the following: i) not decorated by Sac1 (black squares), ii) decorated by Sac1 wt (dark gray squares), and iii) decorated by Sac1 C392S (light gray squares). The PI4P-accepting membrane was composed of POPC/PS (20&#xa0;mol% of PS) and the following: i) not decorated by Sac1 (red squares), ii) decorated by Sac1 wt (orange squares), and iii) decorated by Sac1 C392S (dark yellow squares). The error bars stand for the standard deviation. The presented data represent the trends that were observed in minimum three independent experiments.</p>
</caption>
<graphic xlink:href="fmolb-08-747601-g003.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F3">Figure&#x20;3B</xref> displays the temporal evolution of the autocorrelation amplitude of the CFP signal. The autocorrelation amplitude is inversely proportional to the number of independently moving fluorescent particles (LUV with several dye molecules is seen as one particle) present in the solution. During the PI4P transport, acceptor LUVs were getting enriched in PI4P, leading to an increase of the number of particles recognized by the PI4P biosensor SidC-CFP; as a result, the autocorrelation amplitude lowers (<xref ref-type="fig" rid="F3">Figure&#x20;3B</xref>, curves: black, light dark, light gray, red, and dark yellow). After the transported PI4P is hydrolyzed in the acceptor membrane, there are fewer vesicles newly decorated by SidC-CFP, and at the same time, the PI4P level in the donating vesicles lowers, which leads to their lower fluorescence signal, and in the solution, the free PI4P biosensor SidC-CFP eventually appears. This leads to an initial increase in the autocorrelation amplitude (<xref ref-type="fig" rid="F3">Figure&#x20;3B</xref> orange curve). The explanation of the amplitude increase is not trivial as also the brightness (number of dye molecules per LUV) of the individually moving fluorescent molecules come into play; however, as explained in SI <xref ref-type="fig" rid="F2">Figure&#x20;2</xref>, the rise of the CFP-SidC autocorrelation amplitude fits well to the expected process.</p>
</sec>
<sec id="s3-3">
<title>Regulation of Sac1 Activity</title>
<p>The transport experiments indicate that the catalytic activity of Sac1 is significantly affected by the presence of PS in the membrane where Sac1 is located and PI4P hydrolysis occurs. This observation suggests that the PS in the ER membrane may stimulate Sac1. Therefore, we decided to investigate this process in greater detail. We have examined the hydrolysis of PI4P by Sac1 tethered to membranes of various lipid compositions. We aimed to determine whether the activity is regulated by PS specifically or other negatively charged lipids would act similarly.</p>
<p>For this, we prepared LUVs that contained 0.1&#xa0;mol% PI4P and 20&#xa0;mol% of a selected negatively charged lipid (PS, PI, and PG). The amount of PI4P was chosen to be low to match the sensitivity range of our biosensor SidC. The LUVs also contained a fluorescence marker DiD. Upon addition of the recombinant protein Sac1, we monitored teh decrease of the cross-correlation amplitude between the SidC-CFP and the DiD signal.</p>
<p>Our experiment clearly showed that many of the negatively charged lipids probe (PS, PI, and PG) had a similar and significant impact on the activity of the Sac1 enzyme (<xref ref-type="fig" rid="F4">Figure&#x20;4B</xref>, solid red, green, orange, and black squares for PS, PI, PG, and PC, respectively). LUVs that were not decorated by Sac1 (Sac1 was just in solution, not tethered to LUVs, and the LUVs did not contain DGC-NTA(Ni) for the His-tag attachment) did not display any significant change of the PI4P level in the membrane (measured by the cross-correlation, <xref ref-type="fig" rid="F4">Figure&#x20;4B</xref>, empty squares). However, all the negatively charged lipids significantly supported the reaction. This refers to electrostatic forces that change PI4P hydrolysis by Sac1 rather than a specific PS-related interaction. As the His-tag attachment to the bilayer seems to be crucial for the hydrolysis, the enzyme in our hands has to work mainly in the <italic>cis</italic>-regime (on the same membrane, which it is attached to) (<xref ref-type="bibr" rid="B21">Li and Xu, 2019</xref>). The <italic>trans</italic>-regime (acting on the opposing membrane) would be hardly visible under our conditions since even the unattached enzyme does not significantly hydrolyze PI4P. Even though not excluding the <italic>trans</italic>-regime totally&#x2014;the membranes in cells are more adjacent to each other compared to freely diffusing LUVs&#x2014;our data support the <italic>cis-</italic>activity of Sac1 in agreement with the study of <xref ref-type="bibr" rid="B41">Zewe et&#x20;al. (2018)</xref>.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>PI4P dephosphorylation by Sac1 followed by FCCS. <bold>(A)</bold> Schematics of the FCCS experiment. <bold>(B)</bold> Temporal evolution of the normalized <italic>G</italic>
<sub>cc</sub> amplitude referring on the SidC-CFP release (PI4P dephosphorylation) from the membranes of LUVs. The membrane where the reaction was followed was composed of POPC (black squares), POPC/PS (20&#xa0;mol% of PS) (red squares), POPC/PI (20&#xa0;mol% of PI) (green squares), and POPC/PG (20&#xa0;mol% of PG) (orange squares). The membrane either contained DGS-NTA(Ni) for Sac1 attachment through His-tag (solid squares), or DGS-NTA(Ni) was omitted and Sac1 remained in the solution (empty squares). The amount of PI4P in all the examined membranes was 0.1&#xa0;mol%. The error bars stand for the standard deviation. The presented data represent the trends that were observed in minimum three independent experiments.</p>
</caption>
<graphic xlink:href="fmolb-08-747601-g004.tif"/>
</fig>
<p>To gain molecular insights into the impact of electrostatics on Sac1 interactions with lipid membranes, we performed two series of MD simulations: in the first one, the membrane was composed of POPC (80%) and POPS (20%) lipids, and in the second one, the membrane contained only POPC molecules. <xref ref-type="fig" rid="F5">Figure&#x20;5</xref> illustrates the differences in the positioning of Sac1 at the membrane as observed in the two series of simulations. In the MD simulations with the POPC-POPS bilayer (<xref ref-type="fig" rid="F5">Figures 5A,C</xref>), Sac1 was found to be anchored to the membrane by three loops consisting of amino acid residues 269-273, 306-311, and 336-344. Interestingly, Arg 398 was also found to make transient contacts with the lipids. In contrast, in the MD simulations with the POPC bilayer (<xref ref-type="fig" rid="F5">Figures 5B,D</xref>), Sac1 was observed to vary its orientation relative to the membrane surface, which resulted in more transient contacts of Sac1 with the lipids. The Sac1 regions that were found to make transient contacts with the POPC bilayer were mainly within helix 10 and the loop 336-344. Importantly, Arg 398 did not make any contacts with POPC molecules and was observed to point away from the membrane surface. These differences in the Sac1 position at the POPC-POPS (<xref ref-type="fig" rid="F5">Figures 5A,C</xref>) and POPC (<xref ref-type="fig" rid="F5">Figures 5B,D</xref>) bilayers provide molecular insights into the possible impact of negatively charged lipids on the enzyme activity.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>
<bold>(A,B)</bold> Snapshots from the MD simulations of Sac1 at (A) POPC-POPS and (B) POPC bilayer. The POPC and POPS molecules are shown as thin sticks in gray and black, respectively. The Sac1 molecule is shown in blue except for the regions that make strong (shown in red) or transient (shown in pink) contacts with lipids. The catalytic P-loop is shown in green. Helices 2 and 10 are indicated to show the differences in how Sac1 is oriented relative to <bold>(A)</bold> POPC-POPS and <bold>(B)</bold> POPC bilayers. Arg 398 is shown in stick representation. <bold>(C,D)</bold> Histograms showing which amino acid residues and how frequently make contacts with <bold>(C)</bold> POPC-POPC and <bold>(D)</bold> POPS bilayers.</p>
</caption>
<graphic xlink:href="fmolb-08-747601-g005.tif"/>
</fig>
</sec>
<sec id="s3-4">
<title>PIP2 Transport</title>
<p>PIP2 is an important plasma membrane phosphoinositide that can potentially also affect the PS transport process. It has been shown that it indeed can bind the ORD domain of ORP5/8 (<xref ref-type="bibr" rid="B9">Ghai et&#x20;al., 2017</xref>) belonging to the same subfamily. We have examined the ability of Osh6 to transport PIP2. From the variety of described PIP2-specific biosensors (<xref ref-type="bibr" rid="B31">Saad et&#x20;al., 2006</xref>; <xref ref-type="bibr" rid="B37">Wan et&#x20;al., 2015</xref>; <xref ref-type="bibr" rid="B10">Hertel et&#x20;al., 2020</xref>), we chose a C2 domain from granuphilin fused to CFP (C2<sub>granuphilin</sub>-CFP) (<xref ref-type="bibr" rid="B37">Wan et&#x20;al., 2015</xref>) that had a linear response to PIP2 concentration in the range 0&#x2013;5&#xa0;mol% of PIP2 (<xref ref-type="sec" rid="s10">Supplementary Figure&#x20;S1</xref>).</p>
<p>Our assay revealed that PIP2 transport indeed occurs similarly to the transport of PI4P, and it is sensitized by the presence of PS in the PIP2-accepting membrane. PIP2 transport is slower when PI4P also is present in the same membrane, illustrating that PIP2 and PI4P compete for the same ligand-binding pocket of Osh6, or it binds to a different binding site (<xref ref-type="bibr" rid="B38">Wang et&#x20;al., 2019</xref>) where it might regulate the entry of PI4P to the binging pocket. The coexistence of both PPIns in the plasma membrane suggests that both of them could potentially act as a cargo for the PS exchange; however, PI4P would be used preferentially due to its higher affinity for Osh6. Also, as seen in our previous experiments (<xref ref-type="fig" rid="F1">Figure&#x20;1</xref>), unlike PI4P, PIP2 cannot significantly prevent PS from being transported along its gradient.</p>
<p>The transport experiment consists of two steps: i) the extraction of the lipid from the donor membrane and ii) the insertion of the cargo lipid in the accepting membrane. These two processes cannot be fully separated when the label is only in the accepting membrane. However, this approach can be modified and the label can be also placed in the donor membrane in a separate experiment, and the extraction of PIP2 from the donor membrane can be monitored upon addition of the transporter (<xref ref-type="fig" rid="F6">Figure&#x20;6B</xref>). The combination of the two outcomes allows us to distinguish whether PIP2 was transported to the acceptor membrane or whether it stayed bound to the transporting protein. If Osh6 is added solely to the PIP2-containing vesicles, the PIP2 level in the donor membranes significantly drops (<xref ref-type="fig" rid="F6">Figure&#x20;6B</xref>, dotted curve). The drop approximately corresponds to the amount of added protein, that is, Osh6 binding capacity was around 80% of the available PIP2). When PC LUVs were added, all PIP2 was removed from the donor membrane (<xref ref-type="fig" rid="F6">Figure&#x20;6B</xref>, black curve). However, there is only a little of PIP2 at the accepting membrane (<xref ref-type="fig" rid="F6">Figure&#x20;6A</xref>). Altogether, our data show that PIP2 remains bound to Osh6. If PIP2 LUVs are combined with 20&#xa0;mol% PS-containing LUVs, the extraction significantly slows down as PS competes with PIP2 for the lipid-binding pocket, and in the end, PS is replaced by PIP2 (red curve). Eventually, we combined PIP2-containing LUVs with PI4P-containing LUVs (blue curve); even though it does not correspond to the real biological situation, it helped us see that PIP2 extraction is blocked even more effectively with PI4P as it is a stronger Osh6 binder. PI4P however cannot be replaced by PIP2 and thus stays bound in the pocket.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>PIP2 transport and extraction. <bold>(A)</bold> Temporal evolution of the transport readout. The PIP2-accepting membrane contained increasing amount of PS: 0&#xa0;mol% (black squares), 5&#xa0;mol% (red squares), 10&#xa0;mol% (green squares), and 20&#xa0;mol% (blue squares). The PIP2-donating membrane contained 5&#xa0;mol% of PIP2 (solid squares), or both PIP2 and PI4P, 5&#xa0;mol% each (empty squares). <bold>(B)</bold> Temporal evolution of the normalized cross-correlation amplitude upon Osh6 addition corresponding to the PIP2 biosensor release from the PIP2-donating LUVs. The PIP2-containing LUVs (5&#xa0;mol%) were mixed with LUVs composed of the following: i) POPC: black squares, ii) POPC/PS (20&#xa0;mol% of PS): red squares, and iii) POPC/PI4P (5&#xa0;mol% of PI4P): blue squares. The presented data represent the trends that were observed in minimum three independent experiments.</p>
</caption>
<graphic xlink:href="fmolb-08-747601-g006.tif"/>
</fig>
<p>These experiments revealed the relative tendency to bind the three ligands probed and the consequences for the Osh6-mediated lipid transport.</p>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>Discussion</title>
<p>In our experiments, we have focused on understanding the biophysical determinants of the Osh6-mediated transport of PS. We reconstructed this process in the minimalistic system of artificial lipid bilayers. In our approach, we explored the process as it relates to the three individual participants of the transport: PS, PI4P, and PIP2. For detection of lipid movements, we have employed several lipid biosensors and characterized the range of their dynamic responses in each of the experimental conditions. This is an important consideration as it allows for reasonable planning and interpretation of the experiments.</p>
<p>We first focused on the relative binding affinities of the identified cargo molecules, that is, on the competition of these ligands for the Osh6 binding pocket. The strongest binder is PI4P, which can only be replaced by great excess of PS. Our experiments show that PI4P significantly inhibits the PS transport. This is a critical finding explaining how high PS levels in the PM can be maintained. Since PI4P competitively inserts into the Osh6 binding pocket, it prevents PS transport to be directed back to ER. Binding of PIP2 and PS seems to be of a similar scale in affinity but perhaps of a different mechanism. While PIP2 can affect the PS transport only slightly (<xref ref-type="fig" rid="F1">Figure&#x20;1</xref>), PS excess is necessary for PIP2 replacement in the binding pocket of Osh6 (<xref ref-type="fig" rid="F6">Figure&#x20;6A</xref>), that is, both the ligands compete for Osh6 with approximately equal chances; alternatively, PIP2 uses a different binding mode to the transporter. The PM contains both PIP2 and PI4P (<xref ref-type="bibr" rid="B33">van Meer et&#x20;al., 2008</xref>); therefore, it is more likely that Osh6 will be occupied by PI4P, and thus PIP2 probably does not significantly contribute to the PS transport in the cell, or it modulates the transport by having impact on the lid opening for instance.</p>
<p>Second, we addressed the impact of the Sac1 phosphatase on Osh6 transport function. The presence of Sac1 on the surface of donor membranes restores the PS transport impeded by the presence of PI4P in the acceptor membrane. Apparently, PI4P must be removed at the PS-donating LUVs to continue the transport. We observed that i) PI4P reaches its target membrane prior to its hydrolysis (<xref ref-type="fig" rid="F3">Figure&#x20;3</xref>) and ii) that PI4P hydrolysis by Sac1 was subjected to the presence of negatively charged lipids in the accepting membrane, specifically PS or PI (<xref ref-type="fig" rid="F4">Figure&#x20;4</xref>). These findings together with the orientation of the Sac1 active site with respect to the negatively charged membrane (<xref ref-type="fig" rid="F5">Figure&#x20;5</xref>) show that PI4P is transported and hydrolyzed in two independent steps, that is, it is not hydrolyzed while residing in the Osh6 binding pocket.</p>
<p>Third, we confirm that Sac1 requires PS or PI as an allosteric activator. This was observed previously, and our <italic>in&#x20;vitro</italic> reconstitution system was able to recapitulate this observation (<xref ref-type="bibr" rid="B42">Zhong et&#x20;al., 2012</xref>) and evaluate the PS or PI presence on the ER surface as a prerequisite of the PS transport process.</p>
<p>To sum up (<xref ref-type="fig" rid="F7">Figure&#x20;7</xref>), the most important biophysical determinants of the transport revealed in our study are i) strong PI4P&#x2013;Osh6 binding preference, which disables the PS transport along its concentration gradient; ii) Sac1 activity in the PS donor/PI4P acceptor membrane that removes transported PI4P; and iii) PS and PI activation of Sac1 in the donor membrane, that is, in the&#x20;ER.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Schematic summary of the most important biophysical determinants for the Osh6 transport. <bold>(A)</bold> PS is transported along the gradient spontaneously. <bold>(B)</bold> If PI4P is present in the PS-accepting membrane, Osh6 binding pocket is occupied by PI4P and PS transport is inhibited. <bold>(C)</bold> Sac1 fails to dephosphorylate PI4P if no negatively charged lipid is present in the PI4P-acceptor membrane. <bold>(D)</bold> Both the PS and PI4P transports are occurring if the PI4P release to the acceptor membrane is synchronized with the Sac1 action, which is subjected to the PS (or PI, PG) presence in the acceptor membrane.</p>
</caption>
<graphic xlink:href="fmolb-08-747601-g007.tif"/>
</fig>
<p>In the study of <xref ref-type="bibr" rid="B34">von Filseck et&#x20;al. (2015a)</xref>, the PS/PI4P exchange mechanism is proposed. Our study in agreement with that shows that all the investigated lipids are to a small extent transported spontaneously along the concentration gradient without a need for exchange. However, our experiments proved the strong tendency of Osh6 to transport PS along its concentration gradient with no need of exchange. In addition, we showed that PI4P instead of enhancing of the PS transport along its concentration gradient, it suppresses it. This suppression can be compensated by Sac1 action on the accepting membrane. Altogether, our data explain the directionality of the PS transport and provide rationale for the absence of the backward PS shuffling in the yeast membranes.</p>
<p>Recently, a study emphasizing a strong impact of the aliphatic chains on the transport scenario under non-exchange and exchange conditions (<xref ref-type="bibr" rid="B13">Ikhlef et&#x20;al., 2021</xref>) appeared. The study shows that PS transport can be both accelerated and decelerated with PI4P depending on the saturation of the aliphatic chains and thus on the ligand-binding affinity. It is however worth noticing that the lipid translocation from the membrane environment to the Osh6 binding pocket is not only subjected to the lipid&#x2013;protein interaction but also to the lipid&#x2013;lipid interactions that occur in the membrane from which it is pulled out. Therefore, it is not surprising that our experiments done in significantly less fluid bilayer compared to the previous study (POPC versus DOPC) in some aspects differ as the composition of the lipid phase also has an impact on the lipid propensity for extraction by the transporter.</p>
<p>Our study addresses biophysical determinants of the transport of a single domain protein. Most proteins of the family, however, also have other domains that target the proteins in-between the two membranes, mostly to the membrane contact sites. Usually, the PH domain recognizes the PM, and the FFAT motif binds the VAP proteins in ER, or alternatively, ORP proteins have a transmembrane helix that anchors them directly to ER (<xref ref-type="bibr" rid="B30">Raychaudhuri and Prinz, 2010</xref>). <italic>In vitro</italic>, Osh6 transports the cargo molecules just as a result of their binding preferences and of a cooperation between the transporter and Sac1. <italic>In vivo</italic>, it requires Ist2, a membrane tether that localizes it, similarly as the other domains in more complex ORPs, to the membrane contact sites where it acts. Ist2 does not only localize the transporter but this interaction has also been shown to be critical for the PS transport as such and for the lipid metabolism (<xref ref-type="bibr" rid="B6">D&#x27;Ambrosio et&#x20;al., 2020</xref>; <xref ref-type="bibr" rid="B39">Wong et&#x20;al., 2021</xref>).</p>
</sec>
</body>
<back>
<sec id="s5">
<title>Data Availability Statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="sec" rid="s9">Supplementary Material</xref>; further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s6">
<title>Author Contributions</title>
<p>AE, EB, and JH designed the experiments. BR performed the simulations. JH and EB wrote the&#x20;paper.</p>
</sec>
<sec id="s7">
<title>Funding</title>
<p>JH and AE would like to acknowledge the support of the Czech Science Foundation by grant number 19-18917S. Further, the research was supported jointly by the Czech Science Foundation, grant number 21-27735K (to EB), and by the National Science Centre, Poland, grant number 2020/02/Y/NZ1/00020 (to BR), within the international CEUS-UNISONO program.</p>
</sec>
<sec sec-type="COI-statement" id="s8">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s9">
<title>Publisher&#x2019;s Note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s10">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fmolb.2021.747601/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fmolb.2021.747601/full&#x23;supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="datasheet1.pdf" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Best</surname>
<given-names>R. B.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Shim</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lopes</surname>
<given-names>P. E. M.</given-names>
</name>
<name>
<surname>Mittal</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Feig</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2012</year>). <article-title>Optimization of the Additive CHARMM All-Atom Protein Force Field Targeting Improved Sampling of the Backbone &#x3d5;, &#x3c8; and Side-Chain &#x3c7;1 and &#x3c7;2 Dihedral Angles</article-title>. <source>J.&#x20;Chem. Theor. Comput.</source> <volume>8</volume> (<issue>9</issue>), <fpage>3257</fpage>&#x2013;<lpage>3273</lpage>. <pub-id pub-id-type="doi">10.1021/ct300400x</pub-id> </citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Boura</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Hurley</surname>
<given-names>J.&#x20;H.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Structural Basis for Membrane Targeting by the MVB12-Associated -prism Domain of the Human ESCRT-I MVB12 Subunit</article-title>. <source>Proc. Natl. Acad. Sci.</source> <volume>109</volume> (<issue>6</issue>), <fpage>1901</fpage>&#x2013;<lpage>1906</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1117597109</pub-id> </citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cai</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Deng</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Horenkamp</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Reinisch</surname>
<given-names>K. M.</given-names>
</name>
<name>
<surname>Burd</surname>
<given-names>C. G.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Sac1-Vps74 Structure Reveals a Mechanism to Terminate Phosphoinositide Signaling in the Golgi Apparatus</article-title>. <source>J.&#x20;Cel Biol.</source> <volume>206</volume> (<issue>4</issue>), <fpage>485</fpage>&#x2013;<lpage>491</lpage>. <pub-id pub-id-type="doi">10.1083/jcb.201404041</pub-id> </citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chalupska</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Eisenreichova</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>R&#xf3;&#x17c;ycki</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Rezabkova</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Humpolickova</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Klima</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>Structural Analysis of Phosphatidylinositol 4-kinase III&#x3b2; (PI4KB) - 14-3-3 Protein Complex Reveals Internal Flexibility and Explains 14-3-3 Mediated Protection from Degradation <italic>In Vitro</italic>
</article-title>. <source>J.&#x20;Struct. Biol.</source> <volume>200</volume> (<issue>1</issue>), <fpage>36</fpage>&#x2013;<lpage>44</lpage>. <pub-id pub-id-type="doi">10.1016/j.jsb.2017.08.006</pub-id> </citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chung</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Torta</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Masai</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Lucast</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Czapla</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Tanner</surname>
<given-names>L. B.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>PI4P/phosphatidylserine Countertransport at ORP5- and ORP8-Mediated ER-Plasma Membrane Contacts</article-title>. <source>Science</source> <volume>349</volume> (<issue>6246</issue>), <fpage>428</fpage>&#x2013;<lpage>432</lpage>. <pub-id pub-id-type="doi">10.1126/science.aab1370</pub-id> </citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>D&#x27;Ambrosio</surname>
<given-names>J.&#x20;M.</given-names>
</name>
<name>
<surname>Alban&#xe8;se</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Lipp</surname>
<given-names>N. F.</given-names>
</name>
<name>
<surname>Fleuriot</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Debayle</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Drin</surname>
<given-names>G.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Osh6 Requires Ist2 for Localization to ER-PM Contacts and Efficient Phosphatidylserine Transport in Budding Yeast</article-title>. <source>J.&#x20;Cel Sci</source> <volume>133</volume> (<issue>11</issue>), <fpage>jcs243733</fpage>. <pub-id pub-id-type="doi">10.1242/jcs.243733</pub-id> </citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>de Saint-Jean</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Delfosse</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Douguet</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Chicanne</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Payrastre</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Bourguet</surname>
<given-names>W.</given-names>
</name>
<etal/>
</person-group> (<year>2011</year>). <article-title>Osh4p Exchanges Sterols for Phosphatidylinositol 4-phosphate between Lipid Bilayers</article-title>. <source>J.&#x20;Cel Biol.</source> <volume>195</volume> (<issue>6</issue>), <fpage>965</fpage>&#x2013;<lpage>978</lpage>. <pub-id pub-id-type="doi">10.1083/jcb.201104062</pub-id> </citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dubankova</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Humpolickova</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Klima</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Boura</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Negative Charge and Membrane-Tethered Viral 3B Cooperate to Recruit Viral RNA Dependent RNA Polymerase 3D Pol</article-title>. <source>Sci. Rep.</source> <volume>7</volume>, <fpage>17309</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-017-17621-6</pub-id> </citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ghai</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ferguson</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Brown</surname>
<given-names>A. J.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>ORP5 and ORP8 Bind Phosphatidylinositol-4, 5-biphosphate (PtdIns(4,5)P 2) and Regulate its Level at the Plasma Membrane</article-title>. <source>Nat. Commun.</source> <volume>8</volume>, <fpage>757</fpage>. <pub-id pub-id-type="doi">10.1038/s41467-017-00861-5</pub-id> </citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hertel</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Pott</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Mehta</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Fluorescent Biosensors for Multiplexed Imaging of Phosphoinositide Dynamics</article-title>. <source>ACS Chem. Biol.</source> <volume>15</volume> (<issue>1</issue>), <fpage>33</fpage>&#x2013;<lpage>38</lpage>. <pub-id pub-id-type="doi">10.1021/acschembio.9b00691</pub-id> </citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Humphrey</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Dalke</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Schulten</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>1996</year>). <article-title>VMD: Visual Molecular Dynamics</article-title>. <source>J.&#x20;Mol. Graphics</source> <volume>14</volume> (<issue>1</issue>), <fpage>33</fpage>&#x2013;<lpage>38</lpage>. <pub-id pub-id-type="doi">10.1016/0263-7855(96)00018-5</pub-id> </citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Humpolickova</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Mejdrov&#xe1;</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Matousova</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Nencka</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Boura</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Fluorescent Inhibitors as Tools to Characterize Enzymes: Case Study of the Lipid Kinase Phosphatidylinositol 4-Kinase III&#x3b2; (PI4KB)</article-title>. <source>J.&#x20;Med. Chem.</source> <volume>60</volume> (<issue>1</issue>), <fpage>119</fpage>&#x2013;<lpage>127</lpage>. <pub-id pub-id-type="doi">10.1021/acs.jmedchem.6b01466</pub-id> </citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ikhlef</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Lipp</surname>
<given-names>N.-F.</given-names>
</name>
<name>
<surname>Delfosse</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Fuggetta</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Bourguet</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Magdeleine</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Functional Analyses of phosphatidylserine/PI(4)P Exchangers with Diverse Lipid Species and Membrane Contexts Set Unanticipated Rules on Lipid Transfer</article-title>. <source>Preprint</source>. <pub-id pub-id-type="doi">10.1101/2021.07.16.452025</pub-id> </citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jo</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Im</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Automated Builder and Database of Protein/Membrane Complexes for Molecular Dynamics Simulations</article-title>. <source>Plos One</source> <volume>2</volume> (<issue>9</issue>), <fpage>e880</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0000880</pub-id> </citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jo</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Iyer</surname>
<given-names>V. G.</given-names>
</name>
<name>
<surname>Im</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>CHARMM-GUI: A Web-Based Graphical User Interface for CHARMM</article-title>. <source>J.&#x20;Comput. Chem.</source> <volume>29</volume> (<issue>11</issue>), <fpage>1859</fpage>&#x2013;<lpage>1865</lpage>. <pub-id pub-id-type="doi">10.1002/jcc.20945</pub-id> </citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname>
<given-names>Y. C.</given-names>
</name>
<name>
<surname>Hummer</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Coarse-grained Models for Simulations of Multiprotein Complexes: Application to Ubiquitin Binding</article-title>. <source>J.&#x20;Mol. Biol.</source> <volume>375</volume> (<issue>5</issue>), <fpage>1416</fpage>&#x2013;<lpage>1433</lpage>. <pub-id pub-id-type="doi">10.1016/j.jmb.2007.11.063</pub-id> </citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Klauda</surname>
<given-names>J.&#x20;B.</given-names>
</name>
<name>
<surname>Venable</surname>
<given-names>R. M.</given-names>
</name>
<name>
<surname>Freites</surname>
<given-names>J.&#x20;A.</given-names>
</name>
<name>
<surname>O&#x2019;Connor</surname>
<given-names>J.&#x20;W.</given-names>
</name>
<name>
<surname>Tobias</surname>
<given-names>D. J.</given-names>
</name>
<name>
<surname>Mondragon-Ramirez</surname>
<given-names>C.</given-names>
</name>
<etal/>
</person-group> (<year>2010</year>). <article-title>Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid Types</article-title>. <source>J.&#x20;Phys. Chem. B</source> <volume>114</volume> (<issue>23</issue>), <fpage>7830</fpage>&#x2013;<lpage>7843</lpage>. <pub-id pub-id-type="doi">10.1021/jp101759q</pub-id> </citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Klauda</surname>
<given-names>J.&#x20;B.</given-names>
</name>
<name>
<surname>Monje</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Im</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Improving the CHARMM Force Field for Polyunsaturated Fatty Acid Chains</article-title>. <source>J.&#x20;Phys. Chem. B</source> <volume>116</volume> (<issue>31</issue>), <fpage>9424</fpage>&#x2013;<lpage>9431</lpage>. <pub-id pub-id-type="doi">10.1021/jp304056p</pub-id> </citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Klima</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Chalupska</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>R&#xf3;&#x17c;ycki</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Humpolickova</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Rezabkova</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Silhan</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>Kobuviral Non-Structural 3A Proteins Act as Molecular Harnesses to Hijack the Host ACBD3 Protein</article-title>. <source>Structure</source> <volume>25</volume> (<issue>2</issue>), <fpage>219</fpage>&#x2013;<lpage>230</lpage>. <pub-id pub-id-type="doi">10.1016/j.str.2016.11.021</pub-id> </citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lee</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Swails</surname>
<given-names>J.&#x20;M.</given-names>
</name>
<name>
<surname>Yeom</surname>
<given-names>M. S.</given-names>
</name>
<name>
<surname>Eastman</surname>
<given-names>P. K.</given-names>
</name>
<name>
<surname>Lemkul</surname>
<given-names>J.&#x20;A.</given-names>
</name>
<etal/>
</person-group> (<year>2016</year>). <article-title>CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field</article-title>. <source>J.&#x20;Chem. Theor. Comput.</source> <volume>12</volume> (<issue>1</issue>), <fpage>405</fpage>&#x2013;<lpage>413</lpage>. <pub-id pub-id-type="doi">10.1021/acs.jctc.5b00935</pub-id> </citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>S. H.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>G. K.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Size-Dependent Mechanics of the Adherens Junction Mediated by Cooperative Trans and Cis Bindings</article-title>. <source>J.&#x20;Appl. Mech.-Trans. Asme</source> <volume>86</volume> (<issue>7</issue>), <fpage>071011</fpage>. <pub-id pub-id-type="doi">10.1115/1.4043287</pub-id> </citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Luo</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wasilko</surname>
<given-names>D. J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>Z. Q.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>Structure of the Legionella Virulence Factor, SidC Reveals a Unique PI(4)P-Specific Binding Domain Essential for its Targeting to the Bacterial Phagosome</article-title>. <source>Plos Pathog.</source> <volume>11</volume> (<issue>6</issue>), <fpage>e1004965</fpage>. <pub-id pub-id-type="doi">10.1371/journal.ppat.1004965</pub-id> </citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>MacKerell</surname>
<given-names>A. D.</given-names>
</name>
<name>
<surname>Bashford</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Bellott</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Dunbrack</surname>
<given-names>R. L.</given-names>
</name>
<name>
<surname>Evanseck</surname>
<given-names>J.&#x20;D.</given-names>
</name>
<name>
<surname>Field</surname>
<given-names>M. J.</given-names>
</name>
<etal/>
</person-group> (<year>1998</year>). <article-title>All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins&#x2020;</article-title>. <source>J.&#x20;Phys. Chem. B</source> <volume>102</volume> (<issue>18</issue>), <fpage>3586</fpage>&#x2013;<lpage>3616</lpage>. <pub-id pub-id-type="doi">10.1021/jp973084f</pub-id> </citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Maeda</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Anand</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Chiapparino</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Kumar</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Poletto</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Kaksonen</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2013</year>). <article-title>Interactome Map Uncovers Phosphatidylserine Transport by Oxysterol-Binding Proteins</article-title>. <source>Nature</source> <volume>501</volume> (<issue>7466</issue>), <fpage>257</fpage>&#x2013;<lpage>261</lpage>. <pub-id pub-id-type="doi">10.1038/nature12430</pub-id> </citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Manford</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Xia</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Saxena</surname>
<given-names>A. K.</given-names>
</name>
<name>
<surname>Stefan</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Emr</surname>
<given-names>S. D.</given-names>
</name>
<etal/>
</person-group> (<year>2010</year>). <article-title>Crystal Structure of the Yeast Sac1: Implications for its Phosphoinositide Phosphatase Function</article-title>. <source>Embo J.</source> <volume>29</volume> (<issue>9</issue>), <fpage>1489</fpage>&#x2013;<lpage>1498</lpage>. <pub-id pub-id-type="doi">10.1038/emboj.2010.57</pub-id> </citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mesmin</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Bigay</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Polidori</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Jamecna</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Lacas&#x2010;Gervais</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Antonny</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Sterol Transfer, PI 4P Consumption, and Control of Membrane Lipid Order by Endogenous OSBP</article-title>. <source>Embo J.</source> <volume>36</volume> (<issue>21</issue>), <fpage>3156</fpage>&#x2013;<lpage>3174</lpage>. <pub-id pub-id-type="doi">10.15252/embj.201796687</pub-id> </citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Moser von Filseck</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>&#xa0;opi&#xa0;</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Delfosse</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Vanni</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Jackson</surname>
<given-names>C. L.</given-names>
</name>
<name>
<surname>Bourguet</surname>
<given-names>W.</given-names>
</name>
<etal/>
</person-group> (<year>2015a</year>). <article-title>Phosphatidylserine Transport by ORP/Osh Proteins Is Driven by Phosphatidylinositol 4-phosphate</article-title>. <source>Science</source> <volume>349</volume> (<issue>6246</issue>), <fpage>432</fpage>&#x2013;<lpage>436</lpage>. <pub-id pub-id-type="doi">10.1126/science.aab1346</pub-id> </citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Moser von Filseck</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Vanni</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Mesmin</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Antonny</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Drin</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2015b</year>). <article-title>A Phosphatidylinositol-4-Phosphate Powered Exchange Mechanism to Create a Lipid Gradient between Membranes</article-title>. <source>Nat. Commun.</source> <volume>6</volume>, <fpage>6671</fpage>. <pub-id pub-id-type="doi">10.1038/ncomms7671</pub-id> </citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Olkkonen</surname>
<given-names>V. M.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Oxysterol-Binding Proteins: Sterol and Phosphoinositide Sensors Coordinating Transport, Signaling and Metabolism</article-title>. <source>Prog. Lipid Res.</source> <volume>52</volume> (<issue>4</issue>), <fpage>529</fpage>&#x2013;<lpage>538</lpage>. <pub-id pub-id-type="doi">10.1016/j.plipres.2013.06.004</pub-id> </citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Olkkonen</surname>
<given-names>V.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>OSBP-Related Proteins: Liganding by Glycerophospholipids Opens New Insight into Their Function</article-title>. <source>Molecules</source> <volume>18</volume> (<issue>11</issue>), <fpage>13666</fpage>&#x2013;<lpage>13679</lpage>. <pub-id pub-id-type="doi">10.3390/molecules181113666</pub-id> </citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Phillips</surname>
<given-names>J.&#x20;C.</given-names>
</name>
<name>
<surname>Hardy</surname>
<given-names>D. J.</given-names>
</name>
<name>
<surname>Maia</surname>
<given-names>J.&#x20;D. C.</given-names>
</name>
<name>
<surname>Stone</surname>
<given-names>J.&#x20;E.</given-names>
</name>
<name>
<surname>Ribeiro</surname>
<given-names>J.&#x20;V.</given-names>
</name>
<name>
<surname>Bernardi</surname>
<given-names>R. C.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Scalable Molecular Dynamics on CPU and GPU Architectures with NAMD</article-title>. <source>J.&#x20;Chem. Phys.</source> <volume>153</volume> (<issue>4</issue>), <fpage>044130</fpage>. <pub-id pub-id-type="doi">10.1063/5.0014475</pub-id> </citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Raychaudhuri</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Prinz</surname>
<given-names>W. A.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>The Diverse Functions of Oxysterol-Binding Proteins</article-title>. <source>Annu. Rev. Cel Dev. Biol.</source> <volume>26</volume>, <fpage>157</fpage>&#x2013;<lpage>177</lpage>. <pub-id pub-id-type="doi">10.1146/annurev.cellbio.042308.113334</pub-id> </citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Saad</surname>
<given-names>J.&#x20;S.</given-names>
</name>
<name>
<surname>Miller</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Tai</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Ghanam</surname>
<given-names>R. H.</given-names>
</name>
<name>
<surname>Summers</surname>
<given-names>M. F.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Structural Basis for Targeting HIV-1 Gag Proteins to the Plasma Membrane for Virus Assembly</article-title>. <source>Proc. Natl. Acad. Sci.</source> <volume>103</volume> (<issue>30</issue>), <fpage>11364</fpage>&#x2013;<lpage>11369</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.0602818103</pub-id> </citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Skerle</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Humpolickova</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Johnson</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Rampirova</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Polachova</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Fliegl</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Membrane Protein Dimerization in Cell-Derived Lipid Membranes Measured by FRET with MC Simulations</article-title>. <source>Biophysical J.</source> <volume>118</volume> (<issue>8</issue>), <fpage>1861</fpage>&#x2013;<lpage>1875</lpage>. <pub-id pub-id-type="doi">10.1016/j.bpj.2020.03.011</pub-id> </citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>van Meer</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Voelker</surname>
<given-names>D. R.</given-names>
</name>
<name>
<surname>Feigenson</surname>
<given-names>G. W.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Membrane Lipids: Where They Are and How They Behave</article-title>. <source>Nat. Rev. Mol. Cel Biol</source> <volume>9</volume> (<issue>2</issue>), <fpage>112</fpage>&#x2013;<lpage>124</lpage>. <pub-id pub-id-type="doi">10.1038/nrm2330</pub-id> </citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wahl</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Gregor</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Patting</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Enderlein</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Fast Calculation of Fluorescence Correlation Data with Asynchronous Time-Correlated Single-Photon Counting</article-title>. <source>Opt. Express</source> <volume>11</volume> (<issue>26</issue>), <fpage>3583</fpage>&#x2013;<lpage>3591</lpage>. <pub-id pub-id-type="doi">10.1364/oe.11.003583</pub-id> </citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wan</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Chu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>Insights into the Molecular Recognition of the Granuphilin C2A Domain with PI(4,5)P2</article-title>. <source>Chem. Phys. Lipids</source> <volume>186</volume>, <fpage>61</fpage>&#x2013;<lpage>67</lpage>. <pub-id pub-id-type="doi">10.1016/j.chemphyslip.2015.01.003</pub-id> </citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Qi</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Rae</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>ORP2 Delivers Cholesterol to the Plasma Membrane in Exchange for Phosphatidylinositol 4, 5-Bisphosphate (PI(4,5)P2)</article-title>. <source>Mol. Cel</source> <volume>73</volume> (<issue>3</issue>), <fpage>458</fpage>&#x2013;<lpage>473</lpage>. <pub-id pub-id-type="doi">10.1016/j.molcel.2018.11.014</pub-id> </citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wong</surname>
<given-names>A. K. O.</given-names>
</name>
<name>
<surname>Young</surname>
<given-names>B. P.</given-names>
</name>
<name>
<surname>Loewen</surname>
<given-names>C. J.&#x20;R.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Ist2 Recruits the Lipid Transporters Osh6/7 to ER-PM Contacts to Maintain Phospholipid Metabolism</article-title>. <source>J.&#x20;Cel Biol.</source> <volume>220</volume> (<issue>9</issue>), <fpage>e201910161</fpage>. <pub-id pub-id-type="doi">10.1083/jcb.201910161</pub-id> </citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>E. L.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Jo</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Rui</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>K. C.</given-names>
</name>
<name>
<surname>D&#xe1;vila-Contreras</surname>
<given-names>E. M.</given-names>
</name>
<etal/>
</person-group> (<year>2014</year>). <article-title>CHARMM-GUIMembrane Buildertoward Realistic Biological Membrane Simulations</article-title>. <source>J.&#x20;Comput. Chem.</source> <volume>35</volume> (<issue>27</issue>), <fpage>1997</fpage>&#x2013;<lpage>2004</lpage>. <pub-id pub-id-type="doi">10.1002/jcc.23702</pub-id> </citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zewe</surname>
<given-names>J.&#x20;P.</given-names>
</name>
<name>
<surname>Wills</surname>
<given-names>R. C.</given-names>
</name>
<name>
<surname>Sangappa</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Goulden</surname>
<given-names>B. D.</given-names>
</name>
<name>
<surname>Hammond</surname>
<given-names>G. R.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>SAC1 Degrades its Lipid Substrate PtdIns4P in the Endoplasmic Reticulum to Maintain a Steep Chemical Gradient with Donor Membranes</article-title>. <source>Elife</source> <volume>7</volume>, <fpage>e35588</fpage>. <pub-id pub-id-type="doi">10.7554/eLife.35588</pub-id> </citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhong</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Hsu</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Stefan</surname>
<given-names>C. J.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Patel</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Cosgrove</surname>
<given-names>M. S.</given-names>
</name>
<etal/>
</person-group> (<year>2012</year>). <article-title>Allosteric Activation of the Phosphoinositide Phosphatase Sac1 by Anionic Phospholipids</article-title>. <source>Biochemistry</source> <volume>51</volume> (<issue>15</issue>), <fpage>3170</fpage>&#x2013;<lpage>3177</lpage>. <pub-id pub-id-type="doi">10.1021/bi300086c</pub-id> </citation>
</ref>
</ref-list>
</back>
</article>