On the origins of cells: a hypothesis for the evolutionary transitions from abiotic geochemistry to chemoautotrophic prokaryotes,and from prokaryotes to nucleated cells

All life is organized as cells. Physical compartmentation from the environment and self-organization of self-contained redox reactions are the most conserved attributes of living things, hence inorganic matter with such attributes would be life's most likely forebear. We propose that life evolved in structured iron monosulphide precipitates in a seepage site hydrothermal mound at a redox, pH and temperature gradient between sulphide-rich hydrothermal fluid and iron(II)-containing waters of the Hadean ocean floor. The naturally arising, three-dimensional compartmentation observed within fossilized seepage-site metal sulphide precipitates indicates that these inorganic compartments were the precursors of cell walls and membranes found in free-living prokaryotes. The known capability of FeS and NiS to catalyse the synthesis of the acetyl-methylsulphide from carbon monoxide and methylsulphide, constituents of hydrothermal fluid, indicates that pre-biotic syntheses occurred at the inner surfaces of these metal-sulphide-walled compartments, which furthermore restrained reacted products from diffusion into the ocean, providing sufficient concentrations of reactants to forge the transition from geochemistry to biochemistry. The chemistry of what is known as the RNA-world could have taken place within these naturally forming, catalyticwalled compartments to give rise to replicating systems. Sufficient concentrations of precursors to support replication would have been synthesized in situ geochemically and biogeochemically, with FeS (and NiS) centres playing the central catalytic role. The universal ancestor we infer was not a free-living cell, but rather was confined to the naturally chemiosmotic, FeS compartments within which the synthesis of its constituents occurred. The first free-living cells are suggested to have been eubacterial and archaebacterial chemoautotrophs that emerged more than 3.8 Gyr ago from their inorganic confines. We propose that the emergence of these prokaryotic lineages from inorganic confines occurred independently, facilitated by the independent origins of membrane-lipid biosynthesis: isoprenoid ether membranes in the archaebacterial and fatty acid ester membranes in the eubacterial lineage. The eukaryotes, all of which are ancestrally heterotrophs and possess eubacterial lipids, are suggested to have arisen ca. 2 Gyr ago through symbiosis involving an autotrophic archaebacterial host and a heterotrophic eubacterial symbiont, the common ancestor of mitochondria and hydrogenosomes. The attributes shared by all prokaryotes are viewed as inheritances from their confined universal ancestor. The attributes that distinguish eubacteria and archaebacteria, yet are uniform within the groups, are viewed as relics of their phase of differentiation after divergence from the non-free-living universal ancestor and before the origin of the free-living chemoautotrophic lifestyle. The attributes shared by eukaryotes with eubacteria and archaebacteria, respectively, are viewed as inheritances via symbiosis. The attributes unique to eukaryotes are viewed as inventions specific to their lineage. The origin of the eukaryotic endomembrane system and nuclear membrane are suggested to be the fortuitous result of the expression of genes for eubacterial membrane lipid synthesis by an archaebacterial genetic apparatus in a compartment that was not fully prepared to accommodate such compounds, resulting in vesicles of eubacterial lipids that accumulated in the cytosol around their site of synthesis. Under these premises, the most ancient divide in the living world is that between eubacteria and archaebacteria, yet the steepest evolutionary grade is that between prokaryotes and eukaryotes.     

Unusual pressure dependence of the lateral motion of pyrene-labeled phosphatidylcholine in bipolar lipid vesicles.

The lateral mobility of a pyrene-labeled phosphatidylcholine probe in liposomes containing archaebacterial bipolar lipids has been studied isothermally as a function of pressure. The pressure-dependence of the probe mobility, R, is found to be slightly positive or zero in the temperature range of 17 - 48 degrees C. At temperatures > 48 degrees C, R becomes negative and decreases with temperature. The data indicate that lateral mobility only becomes appreciable at high temperatures. In addition, the R values obtained with other lipid membranes are much lower than that obtained with bipolar liposomes, implying that the membranes of archaebacterial liposomes are laterally immobile, as compared to other lipid membranes.     

Molecular mechanisms of water and solute transport across archaebacterial lipid membranes

Archaebacteria thrive in environments characterized by anaeobiosis, saturated salt, and both high and low extremes of temperature and pH. The bulk of their membrane lipids are polar, characterized by the archaeal structural features typified by ether linkage of the glycerol backbone to isoprenoid chains of constant length, often fully saturated, and with sn-2,3 stereochemistry opposite that of glycerolipids of Bacteria and Eukarya. Also unique to these bacteria are macrocyclic archaeol and membrane spanning caldarchaeol lipids that are found in some extreme thermophiles and methanogens. To define the barrier function of archaebacterial membranes and to examine the effects of these unique structural features on permeabilities, we investigated the water, solute (urea and glycerol), proton, and ammonia permeability of liposomes formed by these lipids. Both the macrocyclic archaeol and caldarchaeol lipids reduced the water, ammonia, urea, and glycerol permeability of liposomes significantly (6-120-fold) compared with diphytanylphosphatidylcholine liposomes. The presence of the ether bond and phytanyl chains did not significantly affect these permeabilities. However, the apparent proton permeability was reduced 3-fold by the presence of an ether bond. The presence of macrocyclic archaeol and caldarchaeol structures further reduced apparent proton permeabilities by 10-17-fold. These results indicate that the limiting mobility of the midplane hydrocarbon region of the membranes formed by macrocyclic archaeol and caldarchaeol lipids play a significant role in reducing the permeability properties of the lipid membrane. In addition, it appears that substituting ether for ester bonds presents an additional barrier to proton flux.     

Ether polar lipids of methanogenic bacteria: structures, comparative aspects, and biosyntheses.

Complete structures of nearly 40 ether polar lipids from seven species of methanogens have been elucidated during the past 10 years. Three kinds of variations of core lipids, macrocyclic archaeol and two hydroxyarchaeols, were identified, in addition to the usual archaeol and caldarchaeol (for the nomenclature of archaeal [archaebacterial] ether lipids, see the text). Polar head groups of methanogen phospholipids include ethanolamine, serine, inositol, N-acetylglucosamine, dimethyl- and trimethylaminopentanetetrol, and glucosaminylinositol. Glucose is the sole hexose moiety of glycolipids in most methanogens, and galactose and mannose have been found in a few species. Methanogen lipids are characterized by their diversity in phosphate-containing polar head groups and core lipids, which in turn can be used for chemotaxonomy of methanogens. This was shown by preliminary simplified analyses of lipid component residues. Core lipid analysis by high-pressure liquid chromatography provides a method of determining the methanogenic biomass in natural samples. There has been significant progress in the biosynthetic studies of methanogen lipids in recent years. In vivo incorporation experiments have led to delineation of the outline of the synthetic route of the diphytanylglycerol ether core. The mechanisms of biosynthesis of tetraether lipids and various polar lipids, and cell-free systems of either lipid synthesis, however, remain to be elucidated. The significance and the origin of archaeal ether lipids is discussed in terms of the lipid composition of bacteria living in a wide variety of environments, the oxygen requirement for biosynthesis of hydrocarbon chains, and the physicochemical properties and functions of lipids as membrane constituents.     

Biology, ecology, and biotechnological applications of anaerobic bacteria adapted to environmental stresses in temperature, pH, salinity, or substrates.

Anaerobic bacteria include diverse species that can grow at environmental extremes of temperature, pH, salinity, substrate toxicity, or available free energy. The first evolved archaebacterial and eubacterial species appear to have been anaerobes adapted to high temperatures. Thermoanaerobes and their stable enzymes have served as model systems for basic and applied studies of microbial cellulose and starch degradation, methanogenesis, ethanologenesis, acetogenesis, autotrophic CO2 fixation, saccharidases, hydrogenases, and alcohol dehydrogenases. Anaerobes, unlike aerobes, appear to have evolved more energy-conserving mechanisms for physiological adaptation to environmental stresses such as novel enzyme activities and stabilities and novel membrane lipid compositions and functions. Anaerobic syntrophs do not have similar aerobic bacterial counterparts. The metabolic end products of syntrophs are potent thermodynamic inhibitors of energy conservation mechanisms, and they require coordinated consumption by a second partner organism for species growth. Anaerobes adapted to environmental stresses and their enzymes have biotechnological applications in organic waste treatment systems and chemical and fuel production systems based on biomass-derived substrates or syngas. These kinds of anaerobes have only recently been examined by biologists, and considerably more study is required before they are fully appreciated by science and technology.     

Morphology, biophysical properties and protein-mediated fusion of archaeosomes

As variance from standard phospholipids of eubacteria and eukaryotes, archaebacterial diether phospholipids contain branched alcohol chains (phytanol) linked to glycerol exclusively with ether bonds. Giant vesicles (GVs) constituted of different species of archaebacterial diether phospholipids and glycolipids (archaeosomes) were prepared by electroformation and observed under a phase contrast and/or fluorescence microscope. Archaebacterial lipids and different mixtures of archaebacterial and standard lipids formed GVs which were analysed for size, yield and ability to adhere to each other due to the mediating effects of certain plasma proteins. GVs constituted of different proportions of archaeal or standard phosphatidylcholine were compared. In nonarchaebacterial GVs (in form of multilamellar lipid vesicles, MLVs) the main transition was detected at T_m = 34. 2°C with an enthalpy of ΔH = 0.68 kcal/mol, whereas in archaebacterial GVs (MLVs) we did not observe the main phase transition in the range between 10 and 70°C. GVs constituted of archaebacterial lipids were subject to attractive interaction mediated by beta 2 glycoprotein I and by heparin. The adhesion constant of beta 2 glycoprotein I – mediated adhesion determined from adhesion angle between adhered GVs was in the range of 10⁻⁸ J/m². In the course of protein mediated adhesion, lateral segregation of the membrane components and presence of thin tubular membranous structures were observed. The ability of archaebacterial diether lipids to combine with standard lipids in bilayers and their compatibility with adhesion-mediating molecules offer further evidence that archaebacterial lipids are appropriate for the design of drug carriers.     

HIV-1, lipid rafts, and antibodies to liposomes: implications for anti-viral-neutralizing antibodies

The human immunodeficiency virus type 1 (HIV-1) is an enveloped virus with a lipid bilayer that contains several glycoproteins that are anchored in, or closely associated with, the membrane surface. The envelope proteins have complex interactions with the lipids both on the host cells and on the target cells. The processes of budding from host cells and entry into target cells occur at sites on the plasma membrane, known as lipid rafts, that represent specialized regions that are rich in cholesterol and sphingolipids. Although the envelope glycoproteins are antigenic molecules that potentially might be used for development of broadly neutralizing antibodies in a vaccine to HIV-1, the development of such antibodies that have broad specificities against primary field isolates of virus has been largely thwarted to date by the ability of the envelope proteins to evade the immune system through various mechanisms. In this review, the interactions of HIV-1 with membrane lipids are summarized. Liposomes are commonly used as models for understanding interactions of proteins with membrane lipids; and liposomes have also been used both as carriers for vaccines, and as antigens for induction of antibodies to liposomal lipids. The possibility is proposed that liposomal lipids, or liposome-protein combinations, could be useful as antigens for inducing broadly neutralizing antibodies to HIV-1.     

Nano-scale resolution of native retinal rod disk membranes reveals differences in lipid composition

Photoreceptors rely on distinct membrane compartments to support their specialized function. Unlike protein localization, identification of critical differences in membrane content has not yet been expanded to lipids, due to the difficulty of isolating domain-specific samples. We have overcome this by using SMA to coimmunopurify membrane proteins and their native lipids from two regions of photoreceptor ROS disks. Each sample''s copurified lipids were subjected to untargeted lipidomic and fatty acid analysis. Extensive differences between center (rhodopsin) and rim (ABCA4 and PRPH2/ROM1) samples included a lower PC to PE ratio and increased LC- and VLC-PUFAs in the center relative to the rim region, which was enriched in shorter, saturated FAs. The comparatively few differences between the two rim samples likely reflect specific protein–lipid interactions. High-resolution profiling of the ROS disk lipid composition gives new insights into how intricate membrane structure and protein activity are balanced within the ROS, and provides a model for future studies of other complex cellular structures.     

A method for the separation and characterization of archaebacterial signature ether lipids

A reproducible high performance liquid chromatography (HPLC) method for the separation of diethers and tetraethers isolated from archaebacterial phospholipids is reported. Fourier transform infrared spectroscopy was used for structural confirmation of these signature lipids. A mixture of tetraethers from a thermoacidophilic archaebacteria was resolved into three major components by the normal phase separation. These components were differentiated by Fourier self-deconvolution of infrared spectra. The application of the HPLC technique to environmental samples may provide an accurate assessment of archaebacterial biomass in various microbial communities.     

HIV-1, lipid rafts,and antibodies to liposomes: implications for anti-viral-neutralizing antibodies (Review)

The human immunodeficiency virus type 1 (HIV-1) is an enveloped virus with a lipid bilayer that contains several glycoproteins that are anchored in, or closely associated with, the membrane surface. The envelope proteins have complex interactions with the lipids both on the host cells and on the target cells. The processes of budding from host cells and entry into target cells occur at sites on the plasma membrane, known as lipid rafts, that represent specialized regions that are rich in cholesterol and sphingolipids. Although the envelope glycoproteins are antigenic molecules that potentially might be used for development of broadly neutralizing antibodies in a vaccine to HIV-1, the development of such antibodies that have broad specificities against primary field isolates of virus has been largely thwarted to date by the ability of the envelope proteins to evade the immune system through various mechanisms. In this review, the interactions of HIV-1 with membrane lipids are summarized. Liposomes are commonly used as models for understanding interactions of proteins with membrane lipids; and liposomes have also been used both as carriers for vaccines, and as antigens for induction of antibodies to liposomal lipids. The possibility is proposed that liposomal lipids, or liposome-protein combinations, could be useful as antigens for inducing broadly neutralizing antibodies to HIV-1.     

Fluid lipid fraction in rod outer segment membrane

Rod outer segment membrane is analyzed using the spin label technique by means of two probes. The solubility of the first label, 2,2,6,6-tetramethylpiperidin-1-oxyl, is correlated with the membrane fluidity which is measured using a stearic acid spin probe. The two values are compared to the solubility-fluidity relationship which characterizes a model system in which all lipids are in a fluid state. The analysis leads to the conclusion that only two thirds of the membrane lipids are fluid. This conclusion is reinforced by the observation that partial lipid removal leaves rigid lipids associated with the rhodopsin molecules.     

Molecular phylogenies based on ribosomal protein L11, L1, L10, and L12 sequences

Summary Available sequences that correspond to the E. coli ribosomal proteins L11, L1, L10, and L12 from eubacteria, archaebacteria, and eukaryotes have been aligned. The alignments were analyzed qualitatively for shared structural features and for conservation of deletions or insertions. The alignments were further subjected to quantitative phylogenetic analysis, and the amino acid identity between selected pairs of sequences was calculated. In general, eubacteria, archaebacteria, and eukaryotes each form coherent and well-resolved nonoverlapping phylogenetic domains. The degree of diversity of the four proteins between the three groups is not uniform. For L11, the eubacterial and archaebacterial proteins are very similar whereas the eukaryotic L11 is clearly less similar. In contrast, in the case of the L12 proteins and to a lesser extent the L10 proteins, the archaebacterial and eukaryotic proteins are similar whereas the eubacterial proteins are different. The eukaryotic L1 equivalent protein has yet to be identified. If the root of the universal tree is near or within the eubacterial domain, our ribosomal protein-based phylogenies indicate that archaebacteria are monophyletic. The eukaryotic lineage appears to originate either near or within the archaebacterial domain. Correspondence to: P. Dennis     

Fluid lipid fraction in rod outer segment membrane.

Rod outer segment membrane is analyzed using the spin resonance label technique by means of two probes. The solubility of the first label,2,2,6,6-tetramethylpiperidin-1-oxyl, is correlated with the membrane fluidity which is measured using a stearic acid spin probe. The two values are compared to the solubility-fluidity relationship which characterizes a model system in which all lipids are in a fluid state. The analysis leads to the conclusion that only two thirds of the membrane lipids are fluid. This conclusion is reinforced by the observation that partial lipid removal leaves rigid lipids associated with the rhodopsin molecules.     

Archaeobacterial ether lipid liposomes (archaeosomes) as novel vaccine and drug delivery systems

Liposomes are artificial, spherical, closed vesicles consisting of one or more lipid bilayer(s). Liposomes made from ester phospholipids have been studied extensively over the last 3 decades as artificial membrane models. Considerable interest has been generated for applications of liposomes in medicine, including their use as diagnostic reagents, as carrier vehicles in vaccine formulations, or as delivery systems for drugs, genes, or cancer imaging agents. The objective of this article is to review the properties and potential applications of novel liposomes made from the membrane lipids of Archaeobacteria (Archaea). These lipids are unique and distinct from those encountered in Eukarya and Bacteria. Polar glycerolipids make up the bulk of the membrane lipids, with the remaining neutral lipids being primarily squalenes and other hydrocarbons. The polar lipids consist of regularly branched, and usually fully saturated, phytanyl chains of 20, 25, or 40 carbon length, with the 20 and 40 being most common. The phytanyl chains are attached via ether bonds to the sn-2,3 carbons of the glycerol backbone(s). It has been shown only recently that total polar lipids of archaeobacteria, and purified lipid fractions therefrom, can form liposomes. We refer to liposomes made with any lipid composition that includes ether lipids characteristic of Archaeobacteria as archaeosomes to distinguish them from vesicles made from the conventional lipids obtained from eukaryotic or eubacterial sources or their synthetic analogs. In general, archaeosomes demonstrate relatively higher stabilities to oxidative stress, high temperature, alkaline pH, action of phospholipases, bile salts, and serum proteins. Some archaeosome formulations can be sterilized by autoclaving, without problems such as fusion or aggregation of the vesicles. The uptake of archaeosomes by phagocytic cells can be up to 50-fold greater than that of conventional liposome formulations. Studies in mice have indicated that systemic administration of several test antigens entrapped within certain archaeosome compositions give humoral immune responses that are comparable to those obtained with the potent but toxic Freund's adjuvant. Archaeosome compositions can be selected to give a prolonged, sustained immune response, and the generation of a memory response. Tissue distribution studies of archaeosomes administered via various systemic and peroral routes indicate potential for targeting to specific organs. All in vitro and in vivo studies performed to date indicate that archaeosomes are safe and do not invoke any noticeable toxicity in mice. The stability, tissue distribution profiles, and adjuvant activity of archaeosome formulations indicate that they may offer a superior alternative to the use of conventional liposomes, at least for some biotechnology applications.     

Particle weights and protein composition of the ribosomal subunits of the extremely thermoacidophilic archaebacterium Caldariella acidophila.

1. The ribosomal subunits of one thermoacidophilic archaebacterium (Caldariella acidophila) and of two reference eubacterial species (Bacillus acidocaldarius, Escherichia coli) were compared with respect to ribosome mass and protein composition by (i) equilibrium-density sedimentation of the particles in CsCl and (ii) gel-electrophoretic estimations of the molecular weights of the protein and the rRNA. 2. By either procedure, it is estimated that synthetically active archaebacterial 30S subunits (52% protein by wt.) are appreciably richer in protein than the corresponding eubacterial particles (31% protein by wt.) 3. The greater protein content of the archaebacterial 30S subunits is accounted for by both a larger number and a greater average molecular weight of the subunit proteins; specifically, C. acidophila 30S subunits yield 28 proteins whose combined mass is 0.6 X 10(6) Da, compared with 20 proteins totalling 0.35 X 10(6) Da mass for eubacterial 30S subunits. 4. No differences in protein number are detected among the large subunits, but C. acidophila 50S subunits exhibit a greater number-average molecular weight of their protein components than do eubacterial 50S particles. 5. Particle weights estimated by either buoyant-density data, or molecular weights of rRNA plus protein, agree to within less than 2%. By either procedure C. acidophila 30S subunits 1.15 X 10(6) Da mass) are estimated to be about 300 000 Da heavier than their eubacterial counterparts (0.87 X 10(6) Da mass); a smaller difference. 0.15 X 10(6) Da, exists between the archaebacterial and the eubacterial 50S subunits (respectively 1.8 X 10(6) and 1.65 X 10(6) Da). It is concluded that the heavier-than-eubacterial mass of the C. acidophila ribosomes resides principally in their smaller subunits.     

Molecular mechanisms of spontaneous curvature and softening in complex lipid bilayer mixtures

Membrane reshaping is an essential biological process. The chemical composition of lipid membranes determines their mechanical properties and thus the energetics of their shape. Hundreds of distinct lipid species make up native bilayers, and this diversity complicates efforts to uncover what compositional factors drive membrane stability in cells. Simplifying assumptions, therefore, are used to generate quantitative predictions of bilayer dynamics based on lipid composition. One assumption commonly used is that “per lipid” mechanical properties are both additive and constant—that they are an intrinsic property of lipids independent of the surrounding composition. Related to this is the assumption that lipid bulkiness, or “shape,” determines its curvature preference, independently of context. In this study, all-atom molecular dynamics simulations on three separate multilipid systems were used to explicitly test these assumptions, applying methodology recently developed to isolate properties of single lipids or nanometer-scale patches of lipids. The curvature preference experienced by populations of lipid conformations were inferred from their redistribution on a dynamically fluctuating bilayer. Representative populations were extracted by both structural similarity and semi-automated hidden Markov model analysis. The curvature preferences of lipid dimers were then determined and compared with an additive model that combines the monomer curvature preference of both the individual lipids. In all three systems, we identified conformational subpopulations of lipid dimers that showed non-additive curvature preference, in each case mediated by a special chemical interaction (e.g., hydrogen bonding). Our study highlights the importance of specific chemical interactions between lipids in multicomponent bilayers and the impact of interactions on bilayer stiffness. We identify two mechanisms of bilayer softening: diffusional softening, driven by the dynamic coupling between lipid distributions and membrane undulations, and conformational softening, driven by the inter-conversion between distinct dimeric conformations.     

Sequence alignment and evolutionary comparison of the L10 equivalent and L12 equivalent ribosomal proteins from archaebacteria,eubacteria, and eucaryotes

Summary The genes corresponding to the L10 and L12 equivalent ribosomal proteins (L10e and L12e) ofEscherichia coli have been cloned and sequenced from two widely divergent species of archaebacteria,Halobacterium cutirubrum andSulfolobus solfataricus. The deduced amino acid sequences of the L10e and L12e proteins have been compared to each other and to available eubacterial and eucaryotic sequences. We have identified the hyman P0 protein as the eucaryotic L10e. The L10e proteins from the three kingdoms were found to be colinear. The eubacterial L10e protein is much shorter than the archaebacterial-eucaryotic proteins because of two large deletions, one internal and one at the carboxy terminus. The archaebacterial and eucaryotic L12e proteins were also colinear; the eubacterial protein is homologous to the archaebacterial and eucaryotic L12e proteins, but has suffered rearrangement through what appear to be gene fusion events. Intraspecies comparisons between L10e and L12e sequences indicate the archaebacterial and eucaryotic L10e proteins contain a partial copy of the L12e protein fused to their carboxy terminus. In the eubacteria most of this fusion has been removed by the carboxy terminal deletion. Within the L12e-derived region, a 26-amino acid-long internal modular sequence reiterated thrice in the archaebacterial L10e, twice in the eucaryotic L10e, and once in the eubacterial L10e was discovered. This modular sequence also appears to be present as a single copy in all L12e proteins and may play a role in L12e dimerization, L10e–L12e complex formation, and the function of L10e–L12e complex in translation. From these sequence comparisons a model depicting the evolutionary progression of the L10e and L12e genes and proteins from the primordial state to the contemporary archaebacterial, eucaryotic, and eubacterial states is presented.     

Microbial signature lipid biomarker analysis – an approach that is still preferred,even amid various method modifications

The lipid composition of microbial communities can indicate their response to changes in the surrounding environment induced by anthropogenic practices, chemical contamination or climatic conditions. A considerable number of analytical techniques exist for the examination of microbial lipids. This article reviews a selection of methods available for environmental samples as applied for lipid extraction, fractionation, derivatization and quantification. The discussion focuses on the origin of the standard methods, the different modified versions developed for investigation of microbial lipids, as well as the advantages and limitations of each. Current modifications to standard methods show a number of improvements for each of the different steps associated with analysis. The advantages and disadvantages of lipid analysis compared to other popular techniques are clarified. Accordingly, the preferential utilization of signature lipid biomarker analysis in current research is considered. It is clear from recent literature that this technique stays relevant – mainly for the variety of microbial properties that can be determined in a single analysis.     

Effects of temperature on ether lipid composition of Caldariella acidophila

The composition of the ether lipids of a strain of Caldariella acidophila, with respect to the different numbers of cyclizations of the biphytanyl components, is shown to differ between the various complex lipid classes, but the degree of cyclization increases systematically with the growth temperature in the range 75–89°. The problem of distinguishing adaptive from phyletic features in archaebacterial lipids is considered.     

Algorithm to catalogue topologies of dynamic lipid hydrogen-bond networks

Lipid membrane interfaces host reactions essential for the functioning of cells. The hydrogen-bonding environment at the membrane interface is particularly important for binding of proteins, drug molecules, and ions. We present here the implementation and applications of a depth-first search algorithm that analyzes dynamic lipid interaction networks. Lipid hydrogen-bond networks sampled transiently during simulations of lipid bilayers are clustered according to main types of topologies that characterize three-dimensional arrangements of lipids connected to each other via short water bridges. We characterize the dynamics of hydrogen-bonded lipid clusters in simulations of model POPE and POPE:POPG membranes that are often used for bacterial membrane proteins, in a model of the Escherichia coli membrane with six different lipid types, and in POPS membranes. We find that all lipids sample dynamic hydrogen-bonded networks with linear, star, or circular arrangements of the lipid headgroups, and larger networks with combinations of these three types of topologies. Overall, linear lipid-water bridges tend to be short. Water-mediated lipid clusters in all membranes with PE lipids tend to be somewhat small, with about four lipids in all membranes studied here. POPS membranes allow circular arrangements of three POPS lipids to be sampled frequently, and complex arrangements of linear, star, and circular paths may also be sampled. These findings suggest a molecular picture of the membrane interface whereby lipid molecules transiently connect in clusters with somewhat small spatial extension.