Dynamic Water Hydrogen-Bond Networks at the Interface of a Lipid Membrane Containing Palmitoyl-Oleoyl Phosphatidylglycerol

"Head-to-Toe" Lipid Properties Govern the Binding and Cargo Transfer of High-Density Lipoprotein

The viscoelastic properties of biological membranes are crucial in controlling cellular functions and are determined primarily by the lipids' composition and structure. This work studies these properties by varying the structure of the constituting lipids in order to influence their interaction with high-density lipoprotein (HDL) particles. Various fluorescence-based techniques were applied to study lipid domains, membrane order, and the overall lateral as well as the molecule-internal glycerol region mobility in HDL-membrane interactions (i.e., binding and/or cargo transfer). The analysis of interactions with HDL particles and various lipid phases revealed that both fully fluid and some gel-phase lipids preferentially interact with HDL particles, although differences were observed in protein binding and cargo exchange. Both interactions were reduced with ordered lipid mixtures containing cholesterol. To investigate the mechanism, membranes were prepared from single-lipid components, enabling step-by-step modification of the lipid building blocks. On a biophysical level, the different mixtures displayed varying stiffness, fluidity, and hydrogen bond network changes. Increased glycerol mobility and a strengthened hydrogen bond network enhanced anchoring interactions, while fluid membranes with a reduced water network facilitated cargo transfer. In summary, the data indicate that different lipid classes are involved depending on the type of interaction, whether anchoring or cargo transfer.

Nanoscale lipid domains determine the dynamic molecular portraits of mixed DOPC/DOPS bilayers in a fluid phase: A computational insight

Lateral heterogeneity, or mosaicity, is a fundamental property inherent to cell membranes that is crucial for their functioning. While microscopic inhomogeneities (e.g. rafts) are easily detected experimentally, lipid domains with nanoscale dimensions (nanoclusters of nanodomains, NDs) resist reliable characterization by instrumental methods. In such a case, important insight can be gained via computer modeling. Here, NDs composed of lipid's head groups in the mixed zwitterionic dioleoylphosphatidylcholine (DOPC) and negatively charged dioleoylphosphatidylserine (DOPS) bilayers were studied by molecular dynamics. A new algorithm has been developed to identify NDs. Unlike most similar methods, it implicitly considers the heterogeneous distribution of lipid head atomic density and does not require subjectively chosen parameters. In DOPS-rich membranes, lipids form more compact and stable NDs due to strong interlipid interactions. In DOPC-rich systems, NDs arise due to the "packing" effect of weakly bound lipid heads. The clustering picture is related to the physical properties of the bilayer surface: DOPS-rich systems show more pronounced surface heterogeneity of hydrophilic/hydrophobic regions compared to DOPC-rich ones. The results obtained are important for the effective quantitative characterization of the "dynamic molecular portrait" of a membrane surface - its "fingerprint" characterizing dynamical distribution of its physicochemical properties.

Graph-based analysis of H-bond networks and unsupervised learning reveal conformational coupling in prion peptide segments

In this study, we employed a comprehensive computational approach to investigate the physical chemistry of the water networks surrounding hydrated peptide segments, as derived from molecular dynamics simulations. Our analysis uncovers a complex interplay of direct and water-mediated hydrogen bonds that intricately weave through the peptides. We demonstrate that these hydrogen bond networks encode critical information about the peptides' conformational behavior, with the dimensionality of these networks showing sensitivity to the peptides' conformations. Additionally, we estimated the free-energy landscape of the peptides across various conformations, revealing that their structures are predominantly characterized by unfolded, partially folded, and folded configurations, resulting in broad and rugged free-energy surfaces due to the numerous degrees of freedom contributed by the surrounding solvent. Importantly, the structured nature of this free-energy landscape becomes obscured when conventional collective variables, such as the number of hydrogen bonds, are used. Our findings provide new insights into the molecular mechanisms that couple protein and solvent degrees of freedom, highlighting their significance in the functioning of biological systems.

A New Theory about Interfacial Proton Diffusion Revisited: The Commonly Accepted Laws of Electrostatics and Diffusion Prevail

The high propensity of protons to stay at interfaces has attracted much attention over the decades. It enables long-range interfacial proton diffusion without relying on titratable residues or electrostatic attraction. As a result, various phenomena manifest themselves, ranging from spillover in material sciences to local proton circuits between proton pumps and ATP synthases in bioenergetics. In an attempt to replace all existing theoretical and experimental insight into the origin of protons' preference for interfaces, TELP, the "Transmembrane Electrostatically-Localized Protons" hypothesis, has been proposed. The TELP hypothesis envisions static H+ and OH- layers on opposite sides of interfaces that are up to 75 µm thick. Yet, the separation at which the electrostatic interaction between two elementary charges is comparable in magnitude to the thermal energy is more than two orders of magnitude smaller and, as a result, the H+ and OH- layers cannot mutually stabilize each other, rendering proton accumulation at the interface energetically unfavorable. We show that (i) the law of electroneutrality, (ii) Fick's law of diffusion, and (iii) Coulomb's law prevail. Using them does not hinder but helps to interpret previously published experimental results, and also helps us understand the high entropy release barrier enabling long-range proton diffusion along the membrane surface.

Graph-Based Analyses of Dynamic Water-Mediated Hydrogen-Bond Networks in Phosphatidylserine: Cholesterol Membranes

Phosphatidylserine lipids are anionic molecules present in eukaryotic plasma membranes, where they have essential physiological roles. The altered distribution of phosphatidylserine in cells such as apoptotic cancer cells, which, unlike healthy cells, expose phosphatidylserine, is of direct interest for the development of biomarkers. We present here applications of a recently implemented Depth-First-Search graph algorithm to dissect the dynamics of transient water-mediated lipid clusters at the interface of a model bilayer composed of 1-palmytoyl-2-oleoyl-sn-glycero-2-phosphatidylserine (POPS) and cholesterol. Relative to a reference POPS bilayer without cholesterol, in the POPS:cholesterol bilayer there is a somewhat less frequent sampling of relatively complex and extended water-mediated hydrogen-bond networks of POPS headgroups. The analysis protocol used here is more generally applicable to other lipid:cholesterol bilayers.

Real Space and Time Imaging of Collective Headgroup Dipole Motions in Zwitterionic Lipid Bilayers

Lipid bilayers are supramolecular structures responsible for a range of processes, such as transmembrane transport of ions and solutes, and sorting and replication of genetic materials, to name just a few. Some of these processes are transient and currently, cannot be visualized in real space and time. Here, we developed an approach using 1D, 2D, and 3D Van Hove correlation functions to image collective headgroup dipole motions in zwitterionic phospholipid bilayers. We show that both 2D and 3D spatiotemporal images of headgroup dipoles are consistent with commonly understood dynamic features of fluids. However, analysis of the 1D Van Hove function reveals lateral transient and re-emergent collective dynamics of the headgroup dipoles-occurring at picosecond time scales-that transmit and dissipate heat at longer times, due to relaxation processes. At the same time, the headgroup dipoles also generate membrane surface undulations due a collective tilting of the headgroup dipoles. A continuous intensity band of headgroup dipole spatiotemporal correlations-at nanometer length and nanosecond time scales-indicates that dipoles undergo stretching and squeezing elastic deformations. Importantly, the above mentioned intrinsic headgroup dipole motions can be externally stimulated at GHz-frequency scale, enhancing their flexoelectric and piezoelectric capabilities (i.e., increased conversion efficiency of mechanical energy into electric energy). In conclusion, we discuss how lipid membranes can provide molecular-level insights about biological learning and memory, and as platforms for the development of the next generation of neuromorphic computers.

Microfluidic synthesis of multilayered lipid-polymer hybrid nanoparticles for the formulation of low solubility drugs

Hybrid phospholipid/block copolymer membranes where polymers and lipids are molecularly mixed or phase-separated into polymer-rich and lipid-rich domains are promising drug delivery materials. Harnessing the chemical diversity of polymers and the biocompatability of lipids is a compelling approach to design the next generation of drug carriers. Here, we report on the development of a microfluidics-based strategy analogous to produce lipid nanoparticles (LNPs) for the nanomanufacturing of multilayered hybrid nanoparticles (HNPs). Using X-ray scattering, Cryo-electron, and polarized microscopy we show that phosphatidylcholine (PC) and PBD-b-PEO (poly(butadiene-block-ethylene oxide)) hybrid membranes can be nanomanufactured by microfluidics into HNPs with dense and multilayered cores which are ideal carriers of low-solubility drugs of the Biopharmaceutical Classification System (BCS) II and IV such as antimalarial DSM265 and Paclitaxel, respectively.

Integrating In Vitro and In Silico Analysis of a Cationic Antimicrobial Peptide Interaction with Model Membranes of Colistin-Resistant Pseudomonas aeruginosa Strains

Bacterial antibiotic resistance is a serious global public health concern. Infections caused by colistin-resistant Pseudomonas aeruginosa (CRPa) strains represent a serious threat due to their considerable morbidity and mortality rates, since most of the current empirical antibiotic therapies are ineffective against these strains. Accordingly, cationic antimicrobial peptides (CAMPs) have emerged as promising alternatives to control resistant bacteria. In this study, the interaction of a CAMP derived from cecropin D-like (∆M2) with model membranes mimicking bacterial biomembranes of wild-type (WTPa) strains of P. aeruginosa and CRPa was evaluated through in vitro and in silico approaches. In vitro interaction was determined by infrared spectroscopy, whereas in silico molecular dynamics was performed to predict specific interactions between amino acids of ∆M2 and lipids of model membrane systems. Experimental analysis showed this peptide interacted with the lipids of bacterial-like model membranes of WTPa and CRPa. In both cases, an increase in the concentration of peptides induced an increase in the phase transition temperature of the lipid systems. On the other hand, the peptides in solution underwent a transition from a random to a helical secondary structure after interacting with the membranes mostly favored in the CRPa system. The α-helix structure percentage for ΔM2 interacting with WTPa and CRPa lipid systems was 6.4 and 33.2%, respectively. Finally, molecular dynamics showed ∆M2 to have the most affinities toward the phospholipids palmitoyl-oleyl-phosphatidylglycerol (POPG) and palmitoyl-oleoyl-phosphatidylethanolamine (POPE) that mimic membranes of WTPa and CRPa, respectively. This work provides clues for elucidating the membrane-associated mechanism of action of ∆M2 against colistin-susceptible and -resistant strains of Pseudomonas aeruginosa.

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.

Examining the Effect of Charged Lipids on Mitochondrial Outer Membrane Dynamics Using Atomistic Simulations

The outer mitochondrial membrane (OMM) is involved in multiple cellular functions such as apoptosis, inflammation and signaling via its membrane-associated and -embedded proteins. Despite the central role of the OMM in these vital phenomena, the structure and dynamics of the membrane have regularly been investigated in silico using simple two-component models. Accordingly, the aim was to generate the realistic multi-component model of the OMM and inspect its properties using atomistic molecular dynamics (MD) simulations. All major lipid components, phosphatidylinositol (PI), phosphatidylcholine (PC), phosphatidylethanolamine (PE), and phosphatidylserine (PS), were included in the probed OMM models. Because increased levels of anionic PS lipids have potential effects on schizophrenia and, more specifically, on monoamine oxidase B enzyme activity, the effect of varying the PS concentration was explored. The MD simulations indicate that the complex membrane lipid composition (MLC) behavior is notably different from the two-component PC-PE model. The MLC changes caused relatively minor effects on the membrane structural properties such as membrane thickness or area per lipid; however, notable effects could be seen with the dynamical parameters at the water-membrane interface. Increase of PS levels appears to slow down lateral diffusion of all lipids and, in general, the presence of anionic lipids reduced hydration and slowed down the PE headgroup rotation. In addition, sodium ions could neutralize the membrane surface, when PI was the main anionic component; however, a similar effect was not seen for high PS levels. Based on these results, it is advisable for future studies on the OMM and its protein or ligand partners, especially when wanting to replicate the correct properties on the water-membrane interface, to use models that are sufficiently complex, containing anionic lipid types, PI in particular.

Translocation of silica nanospheres through giant unilamellar vesicles (GUVs) induced by a high frequency electromagnetic field

Membrane model systems capable of mimicking live cell membranes were used for the first time in studying the effects arising from electromagnetic fields (EMFs) of 18 GHz where membrane permeability was observed following exposure. A present lack of understanding of the mechanisms that drive such a rapid change in membrane permeabilization as well as any structural or dynamic changes imparted on biomolecules affected by high-frequency electromagnetic irradiation limits the use of 18 GHz EMFs in biomedical applications. A phospholipid, 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) labelled with a fluorescent marker 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine-N-(lissamine rhodamine B sulfonyl) (rhodamine-DOPE) was used in constructing the giant unilamellar vesicles (GUVs). After three cycles of exposure, enhanced membrane permeability was observed by the internalisation of hydrophilic silica nanospheres of 23.5 nm and their clusters. All-atom molecular dynamics simulations of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membranes exposed to high frequency electric fields of different field strengths showed that within the simulation timeframe only extremely high strength fields were able to cause an increase in the interfacial water dynamics characterized by water dipole realignments. However, a lower strength, high frequency EMF induced changes of the water hydrogen bond network, which may contribute to the mechanisms that facilitate membrane permeabilization in a longer timeframe.

A nexus of intrinsic dynamics underlies translocase priming

The cytoplasmic ATPase SecA and the membrane-embedded SecYEG channel assemble to form the Sec translocase. How this interaction primes and catalytically activates the translocase remains unclear. We show that priming exploits a nexus of intrinsic dynamics in SecA. Using atomistic simulations, smFRET, and HDX-MS, we reveal multiple dynamic islands that cross-talk with domain and quaternary motions. These dynamic elements are functionally important and conserved. Central to the nexus is a slender stem through which rotation of the preprotein clamp of SecA is biased by ATPase domain motions between open and closed clamping states. An H-bonded framework covering most of SecA enables multi-tier dynamics and conformational alterations with minimal energy input. As a result, cognate ligands select preexisting conformations and alter local dynamics to regulate catalytic activity and clamp motions. These events prime the translocase for high-affinity reception of non-folded preprotein clients. Dynamics nexuses are likely universal and essential in multi-liganded proteins.

Using Graphs of Dynamic Hydrogen-Bond Networks To Dissect Conformational Coupling in a Protein Motor

DExD/H-box proteins are soluble enzymes that couple binding and hydrolysis of adenosine triphosphate (ATP) with reactions involving RNA metabolism or bind and push newly synthesized proteins across bacterial cell membranes. Knowledge of the reaction mechanism of these enzymes could help the development of new therapeutics. In order to explore the mechanism of long-distance conformational coupling in SecA, the DEAD-box motor of the Sec protein secretion in bacteria, we implemented algorithms that provide simplified graph representations of the protein's dynamic hydrogen-bond networks. We find that mutations near the nucleotide-binding site or changes of the nucleotide-binding state of SecA associate with altered dynamics at the preprotein binding domain and identify extended networks of hydrogen bonds that connect the active site of SecA to the region where SecA binds newly synthesized secretory proteins. Water molecules participate in hydrogen-bonded water chains that bridge functional domains of SecA and could contribute to long-distance conformational coupling.

Reactions at Biomembrane Interfaces

Membranes surrounding the biological cell and its internal compartments host proteins that catalyze chemical reactions essential for the functioning of the cell. Rather than being a passive structural matrix that holds membrane-embedded proteins in place, the membrane can largely shape the conformational energy landscape of membrane proteins and impact the energetics of their chemical reaction. Here, we highlight the challenges in understanding how lipids impact the conformational energy landscape of macromolecular membrane complexes whose functioning involves chemical reactions including proton transfer. We review here advances in our understanding of how chemical reactions occur at membrane interfaces gleaned with both theoretical and experimental advances using simple protein systems as guides. Our perspective is that of bridging experiments with theory to understand general physicochemical principles of membrane reactions, with a long term goal of furthering our understanding of the role of the lipids on the functioning of complex macromolecular assemblies at the membrane interface.

Intrinsic Membrane Permeability to Small Molecules

Spontaneous solute and solvent permeation through membranes is of vital importance to human life, be it gas exchange in red blood cells, metabolite excretion, drug/toxin uptake, or water homeostasis. Knowledge of the underlying molecular mechanisms is the sine qua non of every functional assignment to membrane transporters. The basis of our current solubility diffusion model was laid by Meyer and Overton. It correlates the solubility of a substance in an organic phase with its membrane permeability. Since then, a wide range of studies challenging this rule have appeared. Commonly, the discrepancies have their origin in ill-used measurement approaches, as we demonstrate on the example of membrane CO₂ transport. On the basis of the insight that scanning electrochemical microscopy offered into solute concentration distributions in immediate membrane vicinity of planar membranes, we analyzed the interplay between chemical reactions and diffusion for solvent transport, weak acid permeation, and enzymatic reactions adjacent to membranes. We conclude that buffer reactions must also be considered in spectroscopic investigations of weak acid transport in vesicular suspensions. The evaluation of energetic contributions to membrane translocation of charged species demonstrates the compatibility of the resulting membrane current with the solubility diffusion model. A local partition coefficient that depends on membrane penetration depth governs spontaneous membrane translocation of both charged and uncharged molecules. It is determined not only by the solubility in an organic phase but also by other factors like cholesterol concentration and intrinsic electric membrane potentials.

Dynamic water bridging and proton transfer at a surface carboxylate cluster of photosystem II

A hydrogen-bond cluster at a negatively-charged protein interface with a bound protein and long-lived waters might be a proton storage site.

Physiological Calcium Concentrations Slow Dynamics at the Lipid-Water Interface

Phospholipids can interact strongly with ions at physiological concentrations, and these interactions can alter membrane properties. Here, we describe the effects of calcium ions on the dynamics in phospholipid membranes. We used a combination of time-resolved ultrafast two-dimensional infrared spectroscopy and molecular dynamics simulations. We found that millimolar Ca2+ concentrations lead to slower fluctuations in the local environment at the lipid-water interface of membranes with phosphatidylserine. The effect was only observed in bilayers containing anionic phosphatidylserine; membranes composed of only zwitterionic phosphatidylcholine did not experience a slowdown. Local water dynamics were measured using the ester groups as label-free probes and were found to be up to 50% slower with 2.5 mM Ca2+. Molecular dynamics simulations show that Ca2+ primarily binds to the carboxylate group of phosphatidylserines. These findings have implications for apoptotic and diseased cells in which phosphatidylserine is exposed to extracellular calcium and for the biophysical effects of divalent cations on lipid bilayers.