Identification and assessment of cardiolipin interactions with E. coli inner membrane proteins

Effect of cardiolipin on the lamellarity and elongation of liposomes hydrated in PBS

Lamellarity and shape are important factors in the formation of vesicles and determine their role in biological systems and pharmaceutical applications. Cardiolipin (CL) is a major lipid in many biological membranes and exerts a great influence on their structural organization due to its particular structure and physico-chemical properties. Here, we used small-angle X-ray and neutron scattering to study the effects of CL with different acyl chain lengths and saturations (CL14:0, CL18:1, CL18:2) on vesicle morphology and lamellarity in membrane models containing mixtures of phosphatidylcholine and phosphatidylethanolamine with different acyl chain lengths and saturations (C14:0 and C 18:1). Measurements were performed in the presence of Phosphate Buffer Saline (PBS), at 37°C, to better reflect physiological conditions, which resulted in strong effects on vesicle morphology, depending on the type and amount of CL used. The presence of small quantities of CL (from 2.5%) reduced inter-membrane correlations and increased perturbation of the membrane, an effect which is enhanced in the presence of matched shorter saturated acyl chains, and mainly unilamellar vesicles (ULV) are formed. In extruded vesicles, employed for SANS experiments, flattened vesicles are observed partly due to the hypertonic effect of PBS, but also influenced by the type of CL added. Our experimental data from SAXS and SANS revealed a strong dependence on CL content in shaping the membrane microstructure, with an apparent optimum in the PC:CL mixture in terms of promoting reduced correlations, preferred curvature and elongation. However, the use of PBS caused distinct differences from previously published studies in water in terms of vesicle shape, and highlights the need to investigate vesicle formation under physiological conditions in order to be able to draw conclusions about membrane formation in biological systems.

Assessing the Martini 3 protein model: A review of its path and potential

Coarse-grained (CG) protein models have become indispensable tools for studying many biological protein details, from conformational dynamics to the organization of protein macro-complexes, and even the interaction of proteins with other molecules. The Martini force field is one of the most widely used CG models for bio-molecular simulations, partly because of the enormous success of its protein model. With the recent release of a new and improved version of the Martini force field - Martini 3 - a new iteration of its protein model was also made available. The Martini 3 protein force field is an evolution of its Martini 2 counterpart, aimed at improving many of the shortcomings that had been previously identified. In this mini-review, we first provide a general overview of the model and then focus on the successful advances made in the short time since its release, many of which would not have been possible before. Furthermore, we discuss reported limitations, potential directions for model improvement and comment on what the likely future development and application avenues are.

A computational study to assess the pathogenicity of single or combinations of missense variants on respiratory complex I

Variants found in the respiratory complex I (CI) subunit genes encoded by mitochondrial DNA can cause severe genetic diseases. However, it is difficult to establish a priori whether a single or a combination of CI variants may impact oxidative phosphorylation. Here we propose a computational approach based on coarse-grained molecular dynamics simulations aimed at investigating new CI variants. One of the primary CI variants associated with the Leber hereditary optic neuropathy (m.14484T>C/MT-ND6) was used as a test case and was investigated alone or in combination with two additional rare CI variants whose role remains uncertain. We found that the primary variant positioned in the E-channel region, which is fundamental for CI function, stiffens the enzyme dynamics. Moreover, a new mechanism for the transition between π- and α-conformation in the helix carrying the primary variant is proposed. This may have implications for the E-channel opening/closing mechanism. Finally, our findings show that one of the rare variants, located next to the primary one, further worsens the stiffening, while the other rare variant does not affect CI function. This approach may be extended to other variants candidate to exert a pathogenic impact on CI dynamics, or to investigate the interaction of multiple variants.

Structure and mechanism of the K+/H+ exchanger KefC

Intracellular potassium (K+) homeostasis is fundamental to cell viability. In addition to channels, K+ levels are maintained by various ion transporters. One major family is the proton-driven K+ efflux transporters, which in gram-negative bacteria is important for detoxification and in plants is critical for efficient photosynthesis and growth. Despite their importance, the structure and molecular basis for K+-selectivity is poorly understood. Here, we report ~3.1 Å resolution cryo-EM structures of the Escherichia coli glutathione (GSH)-gated K+ efflux transporter KefC in complex with AMP, AMP/GSH and an ion-binding variant. KefC forms a homodimer similar to the inward-facing conformation of Na+/H+ antiporter NapA. By structural assignment of a coordinated K+ ion, MD simulations, and SSM-based electrophysiology, we demonstrate how ion-binding in KefC is adapted for binding a dehydrated K+ ion. KefC harbors C-terminal regulator of K+ conductance (RCK) domains, as present in some bacterial K+-ion channels. The domain-swapped helices in the RCK domains bind AMP and GSH and they inhibit transport by directly interacting with the ion-transporter module. Taken together, we propose that KefC is activated by detachment of the RCK domains and that ion selectivity exploits the biophysical properties likewise adapted by K+-ion-channels.

Multifaceted Activity of Fabimycin: Insights from Molecular Dynamics Studies on Bacterial Membrane Models

Membranes─cells' essential scaffolds─are valid molecular targets for substances with an antimicrobial effect. While certain substances, such as octenidine, have been developed to target membranes for antimicrobial purposes, the recently reported molecule, fabimycin (F2B)─a novel agent targeting drug-resistant Gram-negative bacteria─has not received adequate attention regarding its activity on membranes in the literature. The following study aims to investigate the effects of F2B on different bacterial membrane models, including simple planar bilayers and more complex bilayer systems that mimic the Escherichia coli shell equipped with double inner and outer bilayers. Our results show that F2B exhibited more pronounced interactions with bacterial membrane systems compared to the control PC system. Furthermore, we observed significant changes in local membrane property homeostasis in both the inner and outer membrane models, specifically in the case of lateral diffusion, membrane thickness, and membrane resilience (compressibility, tilt). Finally, our results showed that the effect of F2B differed in a complex system and a single membrane system. Our study provides new insights into the multifaceted activity of F2B, demonstrating its potential to disrupt bacterial membrane homeostasis, indicating that its activity extends the currently known mechanism of FabI enzyme inhibition. This disruption, coupled with the ability of F2B to penetrate the outer membrane layers, sheds new light on the behavior of this antimicrobial molecule. This highlights the importance of the interaction with the membrane, crucial in combating bacterial infections, particularly those caused by drug-resistant strains.

Martini-3 Coarse-Grained Models for the Bacterial Lipopolysaccharide Outer Membrane of Escherichia coli

The outer lipopolysaccharide (LPS) membrane of Gram-negative bacteria forms the main barrier for transport of antimicrobial molecules into the bacterial cell. In this study we develop coarse-grained models for the outer membrane of Escherichia coli in the Martini-3 framework. The coarse-grained model force field was parametrized and validated using all-atom simulations of symmetric membranes of lipid A and rough LPS as well as a complete asymmetric membrane of LPS with the O-antigen. The bonded parameters were obtained using an iterative refinement procedure with target bonded distributions obtained from all-atom simulations. The membrane thickness, area of the LPS, and density distributions for the different regions as well as the water and ion densities in Martini-3 simulations show excellent agreement with the all-atom data. Additionally the solvent accessible surface area for individual molecules in water was found to be in good agreement. The binding of calcium ions with phosphate and carboxylate moieties of LPS is accurately captured in the Martini-3 model, indicative of the integrity of the highly negatively charged LPS molecules in the outer membranes of Gram-negative bacteria. The melting transition of the coarse-grained lipid A membrane model was found to occur between 300 and 310 K, and the model captured variations in area per LPS, order parameter, and membrane thickness across the melting transition. Our study reveals that the proposed Martini-3 models for LPS are able to capture the physicochemical balance of the complex sugar architecture of the outer membrane of Escherichia coli. The coarse-grained models developed in this study would be useful for determining membrane protein interactions and permeation of potential antimicrobials through bacterial membranes at mesoscopic spatial and temporal scales.

Improved Protein Dynamics and Hydration in the Martini3 Coarse-Grain Model

The Martini coarse-grain force-field has emerged as an important framework to probe cellular processes at experimentally relevant time- and length-scales. However, the recently developed version, the Martini3 force-field with the implemented Go̅ model (Martini3Go̅), as well as previous variants of the Martini model have not been benchmarked and rigorously tested for globular proteins. In this study, we consider three globular proteins, ubiquitin, lysozyme, and cofilin, and compare protein dynamics and hydration with observables from experiments and all-atom simulations. We show that the Martini3Go̅ model is able to accurately model the structural and dynamic features of small globular proteins. Overall, the structural integrity of the proteins is maintained, as validated by contact maps, radii of gyration (Rg), and SAXS profiles. The chemical shifts predicted from the ensemble sampled in the simulations are consistent with the experimental data. Further, a good match is observed in the protein-water interaction energetics, and the hydration levels of the residues are similar to atomistic simulations. However, the protein-water interaction dynamics is not accurately represented and appears to depend on the protein structural complexity, residue specificity, and water dynamics. Our work is a step toward testing and assessing the Martini3Go̅ model and provides insights into future efforts to refine Martini models with improved solvation effects and better correspondence to the underlying all-atom systems.

Disruption of CFAP418 interaction with lipids causes widespread abnormal membrane-associated cellular processes in retinal degenerations

Syndromic ciliopathies and retinal degenerations are large heterogeneous groups of genetic diseases. Pathogenic variants in the CFAP418 gene may cause both disorders, and its protein sequence is evolutionarily conserved. However, the disease mechanism underlying CFAP418 mutations has not been explored. Here, we apply quantitative lipidomic, proteomic, and phosphoproteomic profiling and affinity purification coupled with mass spectrometry to address the molecular function of CFAP418 in the retina. We show that CFAP418 protein binds to the lipid metabolism precursor phosphatidic acid (PA) and mitochondrion-specific lipid cardiolipin but does not form a tight and static complex with proteins. Loss of Cfap418 in mice disturbs membrane lipid homeostasis and membrane-protein associations, which subsequently causes mitochondrial defects and membrane-remodeling abnormalities across multiple vesicular trafficking pathways in photoreceptors, especially the endosomal sorting complexes required for transport (ESCRT) pathway. Ablation of Cfap418 also increases the activity of PA-binding protein kinase Cα in the retina. Overall, our results indicate that membrane lipid imbalance is a pathological mechanism underlying syndromic ciliopathies and retinal degenerations which is associated with other known causative genes of these diseases.

LipIDens: simulation assisted interpretation of lipid densities in cryo-EM structures of membrane proteins

Cryo-electron microscopy (cryo-EM) enables the determination of membrane protein structures in native-like environments. Characterising how membrane proteins interact with the surrounding membrane lipid environment is assisted by resolution of lipid-like densities visible in cryo-EM maps. Nevertheless, establishing the molecular identity of putative lipid and/or detergent densities remains challenging. Here we present LipIDens, a pipeline for molecular dynamics (MD) simulation-assisted interpretation of lipid and lipid-like densities in cryo-EM structures. The pipeline integrates the implementation and analysis of multi-scale MD simulations for identification, ranking and refinement of lipid binding poses which superpose onto cryo-EM map densities. Thus, LipIDens enables direct integration of experimental and computational structural approaches to facilitate the interpretation of lipid-like cryo-EM densities and to reveal the molecular identities of protein-lipid interactions within a bilayer environment. We demonstrate this by application of our open-source LipIDens code to ten diverse membrane protein structures which exhibit lipid-like densities.

Rotation of the c-Ring Promotes the Curvature Sorting of Monomeric ATP Synthases

ATP synthases are proteins that catalyse the formation of ATP through the rotatory movement of their membrane-spanning subunit. In mitochondria, ATP synthases are found to arrange as dimers at the high-curved edges of cristae. Here, a direct link is explored between the rotatory movement of ATP synthases and their preference for curved membranes. An active curvature sorting of ATP synthases in lipid nanotubes pulled from giant vesicles is found. Coarse-grained simulations confirm the curvature-seeking behaviour of rotating ATP synthases, promoting reversible and frequent protein-protein contacts. The formation of transient protein dimers relies on the membrane-mediated attractive interaction of the order of 1.5 kB T produced by a hydrophobic mismatch upon protein rotation. Transient dimers are sustained by a conic-like arrangement characterized by a wedge angle of θ ≈ 50°, producing a dynamic coupling between protein shape and membrane curvature. The results suggest a new role of the rotational movement of ATP synthases for their dynamic self-assembly in biological membranes.

OPA1 helical structures give perspective to mitochondrial dysfunction

Dominant optic atrophy is one of the leading causes of childhood blindness. Around 60-80% of cases1 are caused by mutations of the gene that encodes optic atrophy protein 1 (OPA1), a protein that has a key role in inner mitochondrial membrane fusion and remodelling of cristae and is crucial for the dynamic organization and regulation of mitochondria2. Mutations in OPA1 result in the dysregulation of the GTPase-mediated fusion process of the mitochondrial inner and outer membranes3. Here we used cryo-electron microscopy methods to solve helical structures of OPA1 assembled on lipid membrane tubes, in the presence and absence of nucleotide. These helical assemblies organize into densely packed protein rungs with minimal inter-rung connectivity, and exhibit nucleotide-dependent dimerization of the GTPase domains-a hallmark of the dynamin superfamily of proteins4. OPA1 also contains several unique secondary structures in the paddle domain that strengthen its membrane association, including membrane-inserting helices. The structural features identified in this study shed light on the effects of pathogenic point mutations on protein folding, inter-protein assembly and membrane interactions. Furthermore, mutations that disrupt the assembly interfaces and membrane binding of OPA1 cause mitochondrial fragmentation in cell-based assays, providing evidence of the biological relevance of these interactions.

Regulation of Bacterial Two-Component Systems by Cardiolipin

The regulation of membrane protein activity for cellular functions is critically dependent on the composition of phospholipid membranes. Cardiolipin, a unique phospholipid found in bacterial membranes and mitochondrial membranes of eukaryotes, plays a crucial role in stabilizing membrane proteins and maintaining their function. In the human pathogen Staphylococcus aureus, the SaeRS two-component system (TCS) controls the expression of key virulence factors essential for the bacterium's virulence. The SaeS sensor kinase activates the SaeR response regulator via phosphoryl transfer to bind its gene target promoters. In this study, we report that cardiolipin is critical for sustaining the full activity of SaeRS and other TCSs in S. aureus. The sensor kinase protein SaeS binds directly to cardiolipin and phosphatidylglycerol, enabling SaeS activity. Elimination of cardiolipin from the membrane reduces SaeS kinase activity, indicating that bacterial cardiolipin is necessary for modulating the kinase activities of SaeS and other sensor kinases during infection. Moreover, the deletion of cardiolipin synthase genes cls1 and cls2 leads to reduced cytotoxicity to human neutrophils and lower virulence in a mouse model of infection. These findings suggest a model where cardiolipin modulates the kinase activity of SaeS and other sensor kinases after infection to adapt to the hostile environment of the host and expand our knowledge of how phospholipids contribute to membrane protein function.

Mutational analysis of a conserved positive charge in the c-ring of E. coli ATP synthase

F1Fo ATP synthase is a ubiquitous molecular motor that utilizes a rotary mechanism to synthesize adenosine triphosphate (ATP), the fundamental energy currency of life. The membrane-embedded Fo motor converts the electrochemical gradient of protons into rotation, which is then used to drive the conformational changes in the soluble F1 motor that catalyze ATP synthesis. In E. coli, the Fo motor is composed of a c10 ring (rotor) alongside subunit a (stator), which together provide two aqueous half channels that facilitate proton translocation. Previous work has suggested that Arg50 and Thr51 on the cytoplasmic side of each subunit c are involved in the proton translocation process, and positive charge is conserved in this region of subunit c. To further investigate the role of these residues and the chemical requirements for activity at these positions, we generated 13 substitution mutants and assayed their in vitro ATP synthesis, H+ pumping, and passive H+ permeability activities, as well as the ability of mutants to carry out oxidative phosphorylation in vivo. While polar and hydrophobic mutations were generally tolerated in either position, introduction of negative charge or removal of polarity caused a substantial defect. We discuss the possible effects of altered electrostatics on the interaction between the rotor and stator, water structure in the aqueous channel, and interaction of the rotor with cardiolipin.

Renovating a double fence with or without notifying the next door and across the street neighbors: why the biogenic cytoplasmic membrane of Gram-negative bacteria display asymmetry?

The complex two-membrane organization of the envelope of Gram-negative bacteria imposes an unique biosynthetic and topological constraints that can affect translocation of lipids and proteins synthesized on the cytoplasm facing leaflet of the cytoplasmic (inner) membrane (IM), across the IM and between the IM and outer membrane (OM). Balanced growth of two membranes and continuous loss of phospholipids in the periplasmic leaflet of the IM as metabolic precursors for envelope components and for translocation to the OM requires a constant supply of phospholipids in the IM cytosolic leaflet. At present we have no explanation as to why the biogenic E. coli IM displays asymmetry. Lipid asymmetry is largely related to highly entropically disfavored, unequal headgroup and acyl group asymmetries which are usually actively maintained by active mechanisms. However, these mechanisms are largely unknown for bacteria. Alternatively, lipid asymmetry in biogenic IM could be metabolically controlled in order to maintain uniform bilayer growth and asymmetric transmembrane arrangement by balancing temporally the net rates of synthesis and flip-flop, inter IM and OM bidirectional flows and bilayer chemical and physical properties as spontaneous response. Does such flippase-less or 'lipid only", 'passive' mechanism of generation and maintenance of lipid asymmetry exists in the IM? The driving force for IM asymmetry can arise from the packing requirements imposed upon the bilayer system during cell division through disproportional distribution of two negatively curved phospholipids, phosphatidylethanolamine and cardiolipin, with consistent reciprocal tendency to increase and decrease lipid order in each membrane leaflet respectively.

Regulation of Bacterial Two-Component Systems by Cardiolipin

The composition of phospholipid membranes is critical to regulating the activity of membrane proteins for cellular functions. Cardiolipin is a unique phospholipid present within the bacterial membrane and mitochondria of eukaryotes and plays a role in maintaining the function and stabilization of membrane proteins. Here, we report that, in the human pathogen Staphylococcus aureus, cardiolipin is required for full activity of the SaeRS two-component system (TCS). Deletion of the cardiolipin synthase genes, cls1 , and cls2 , reduces the basal activity of SaeRS and other TCSs. Cardiolipin is an indispensable requisite for Sae activation mediated by human neutrophil peptides (HNPs) in the stationary growth phase but not mandatory for Sae induction in the exponential growth phase. Ectopic expression with cls2 , but not with cls1 , in the cls1 cls2 double mutant fully restores Sae activity. Elimination of cardiolipin from the membranes results in decreased kinase activity of the sensor protein SaeS. Purified SaeS protein directly binds to cardiolipin as well as phosphatidylglycerol. A strain lacking cls2 or cls1cls2 renders S. aureus less cytotoxic to human neutrophils and less virulent in a mouse model of infection. Our findings suggest that cardiolipin enables a pathogen to confer virulence by modulating the kinase activity of SaeS and other sensor kinases upon infection.

Supramolecular organization and dynamics of mannosylated phosphatidylinositol lipids in the mycobacterial plasma membrane

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis (TB), a disease that claims ~1.6 million lives annually. The current treatment regime is long and expensive, and missed doses contribute to drug resistance. Therefore, development of new anti-TB drugs remains one of the highest public health priorities. Mtb has evolved a complex cell envelope that represents a formidable barrier to antibiotics. The Mtb cell envelop consists of four distinct layers enriched for Mtb specific lipids and glycans. Although the outer membrane, comprised of mycolic acid esters, has been extensively studied, less is known about the plasma membrane, which also plays a critical role in impacting antibiotic efficacy. The Mtb plasma membrane has a unique lipid composition, with mannosylated phosphatidylinositol lipids (phosphatidyl-myoinositol mannosides, PIMs) comprising more than 50% of the lipids. However, the role of PIMs in the structure and function of the membrane remains elusive. Here, we used multiscale molecular dynamics (MD) simulations to understand the structure-function relationship of the PIM lipid family and decipher how they self-organize to shape the biophysical properties of mycobacterial plasma membranes. We assess both symmetric and asymmetric assemblies of the Mtb plasma membrane and compare this with residue distributions of Mtb integral membrane protein structures. To further validate the model, we tested known anti-TB drugs and demonstrated that our models agree with experimental results. Thus, our work sheds new light on the organization of the mycobacterial plasma membrane. This paves the way for future studies on antibiotic development and understanding Mtb membrane protein function.

Cardiolipin, and not monolysocardiolipin, preferentially binds to the interface of complexes III and IV

The mitochondrial electron transport chain comprises a series of protein complexes embedded in the inner mitochondrial membrane that generate a proton motive force via oxidative phosphorylation, ultimately generating ATP. These protein complexes can oligomerize to form larger structures called supercomplexes. Cardiolipin (CL), a conical lipid, unique within eukaryotes to the inner mitochondrial membrane, has proven essential in maintaining the stability and function of supercomplexes. Monolysocardiolipin (MLCL) is a CL variant that accumulates in people with Barth syndrome (BTHS). BTHS is caused by defects in CL biosynthesis and characterised by abnormal mitochondrial bioenergetics and destabilised supercomplexes. However, the mechanisms by which MLCL causes pathogenesis remain unclear. Here, multiscale molecular dynamics characterise the interactions of CL and MLCL with yeast and mammalian mitochondrial supercomplexes containing complex III (CIII) and complex IV (CIV). Coarse-grained simulations reveal that both CL and MLCL bind to sites at the interface between CIII and CIV of the supercomplex. Free energy perturbation calculations show that MLCL interaction is weaker than that of CL and suggest that interaction with CIV drives this difference. Atomistic contact analyses show that, although interaction with CIII is similar for CL and MLCL, CIV makes more contacts with CL than MLCL, demonstrating that CL is a more successful "glue" between the two complexes. Simulations of the human CIII₂CIV supercomplex show that this interface site is maintained between species. Our study suggests that MLCL accumulation in people with BTHS disrupts supercomplex stability by formation of relatively weak interactions at the interface lipid binding site.

Mitochondrial membrane models built from native lipid extracts: Interfacial and transport properties

The mitochondrion is an essential organelle enclosed by two membranes whose functionalities depend on their very specific protein and lipid compositions. Proteins from the outer mitochondrial membrane (OMM) are specialized in mitochondrial dynamics and mitophagy, whereas proteins of the inner mitochondrial membrane (IMM) have dedicated functions in cellular respiration and apoptosis. As for lipids, the OMM is enriched in glycerophosphatidyl choline but cardiolipin is exclusively found within the IMM. Though the lipid topology and distribution of the OMM and IMM are known since more than four decades, little is known about the interfacial and dynamic properties of the IMM and OMM lipid extracts. Here we build monolayers, supported bilayers and giant unilamellar vesicles (GUVs) of native OMM and IMM lipids extracts from porcine heart. Additionally, we perform a comparative analysis on the interfacial, phase immiscibility and mechanical properties of both types of extract. Our results show that IMM lipids form more expanded and softer membranes than OMM lipids, allowing a better understanding of the physicochemical and biophysical properties of mitochondrial membranes.

New Roles for HAMP Domains: the Tri-HAMP Region of Pseudomonas aeruginosa Aer2 Controls Receptor Signaling and Cellular Localization

The Aer2 chemoreceptor from Pseudomonas aeruginosa is an O₂ sensor involved in stress responses, virulence, and tuning the behavior of the chemotaxis (Che) system. Aer2 is the sole receptor of the Che2 system. It is soluble, but membrane associated, and forms complexes at the cell pole during stationary phase. The domain arrangement of Aer2 is unusual, with a PAS sensing domain sandwiched between five HAMP domains, followed by a C-terminal kinase-control output domain. The first three HAMP domains form a poly-HAMP chain N-terminal to the PAS sensing domain. HAMP domains are often located between signal input and output domains, where they transduce signals. Given that HAMP1 to 3 reside N-terminal to the input-output pathway, we undertook a systematic examination of their function in Aer2. We found that HAMP1 to 3 influence PAS signaling over a considerable distance, as the majority of HAMP1, 2 and 3 mutations, and deletions of helical phase stutters, led to nonresponsive signal-off or off-biased receptors. PAS signal-on lesions that mimic activated Aer2 also failed to override N-terminal HAMP signal-off replacements. This indicates that HAMP1 to 3 are critical coupling partners for PAS signaling and likely function as a cohesive unit and moveable scaffold to correctly orient and poise PAS dimers for O₂-mediated signaling in Aer2. HAMP1 additionally controlled the clustering and polar localization of Aer2 in P. aeruginosa. Localization was not driven by HAMP1 charge, and HAMP1 signal-off mutants still localized. Employing HAMP as a clustering and localization determinant, as well as a facilitator of PAS signaling, are newly recognized roles for HAMP domains. IMPORTANCE P. aeruginosa is an opportunistic pathogen that interprets environmental stimuli via 26 chemoreceptors that signal through 4 distinct chemosensory systems. The second chemosensory system, Che2, contains a receptor named Aer2 that senses O₂ and mediates stress responses and virulence and tunes chemotactic behavior. Aer2 is membrane associated, but soluble, and has three N-terminal HAMP domains (HAMP1 to 3) that reside outside the signal input-output pathway of Aer2. In this study, we determined that HAMP1 to 3 facilitate O₂-dependent signaling from the PAS sensing domain and that HAMP1 controls the formation of Aer2-containing polar foci in P. aeruginosa. Both of these are newly recognized roles for HAMP domains that may be applicable to other non-signal-transducing HAMP domains and poly-HAMP chains.

The AAA+ ATPase RavA and its binding partner ViaA modulate E. coli aminoglycoside sensitivity through interaction with the inner membrane

Enteric bacteria have to adapt to environmental stresses in the human gastrointestinal tract such as acid and nutrient stress, oxygen limitation and exposure to antibiotics. Membrane lipid composition has recently emerged as a key factor for stress adaptation. The E. coli ravA-viaA operon is essential for aminoglycoside bactericidal activity under anaerobiosis but its mechanism of action is unclear. Here we characterise the VWA domain-protein ViaA and its interaction with the AAA+ ATPase RavA, and find that both proteins localise at the inner cell membrane. We demonstrate that RavA and ViaA target specific phospholipids and subsequently identify their lipid-binding sites. We further show that mutations abolishing interaction with lipids restore induced changes in cell membrane morphology and lipid composition. Finally we reveal that these mutations render E. coli gentamicin-resistant under fumarate respiration conditions. Our work thus uncovers a ravA-viaA-based pathway which is mobilised in response to aminoglycosides under anaerobiosis and engaged in cell membrane regulation.

The Hidden Intricacies of Aquaporins: Remarkable Details in a Common Structural Scaffold

Evolution turned aquaporins (AQPs) into the most efficient facilitators of passive water flow through cell membranes at no expense of solute discrimination. In spite of a plethora of solved AQP structures, many structural details remain hidden. Here, by combining extensive sequence- and structural-based analysis of a unique set of 20 non-redundant high-resolution structures and molecular dynamics simulations of four representatives, key aspects of AQP stability, gating, selectivity, pore geometry, and oligomerization, with a potential impact on channel functionality, are identified. The general view of AQPs possessing a continuous open water pore is challenged and it is depicted that AQPs' selectivity is not exclusively shaped by pore-lining residues but also by the relative arrangement of transmembrane helices. Moreover, this analysis reveals that hydrophobic interactions constitute the main determinant of protein thermal stability. Finally, a numbering scheme of the conserved AQP scaffold is established, facilitating direct comparison of, for example, disease-causing mutations and prediction of potential structural consequences. Additionally, the results pave the way for the design of optimized AQP water channels to be utilized in biotechnological applications.

Systematic simulation of the interactions of pleckstrin homology domains with membranes

Pleckstrin homology (PH) domains can recruit proteins to membranes by recognition of phosphatidylinositol phosphate (PIP) lipids. Several family members are linked to diseases including cancer. We report the systematic simulation of the interactions of 100 mammalian PH domains with PIP-containing membranes. The observed PIP interaction hotspots recapitulate crystallographic binding sites and reveal a number of insights: (i) The β1 and β2 strands and their connecting loop constitute the primary PIP interaction site but are typically supplemented by interactions at the β3-β4 and β5-β6 loops; (ii) we reveal exceptional cases such as the Exoc8 PH domain; (iii) PH domains adopt different membrane-bound orientations and induce clustering of anionic lipids; and (iv) beyond family-level insights, our dataset sheds new light on individual PH domains, e.g., by providing molecular detail of secondary PIP binding sites. This work provides a global view of PH domain/membrane association involving multivalent association with anionic lipids.

Cardiolipin Biosynthesis Genes Are Not Required for Salmonella enterica Serovar Typhimurium Pathogenesis in C57BL/6J Mice

Salmonella enterica serovar Typhimurium is an intracellular pathogen that parasitizes macrophages from within a vacuole. The vacuolar environment prompts the bacterium to regulate the lipid composition of the outer membrane (OM), and this influences host inflammation. S. Typhimurium regulates the levels of acidic glycerophospholipids known as cardiolipins (CL) within the OM, and mitochondrial CL molecules can prime and activate host inflammasomes. However, the contribution of S. Typhimurium’s CL biosynthesis genes to intracellular survival, inflammasome activation, and pathogenesis had not been examined. S. Typhimurium genes encode three CL synthases. Single, double, and triple mutants were constructed. Similar to other Enterobacteriaceae, ClsA is the primary CL synthase for S. Typhimurium during logarithmic growth, while ClsB and ClsC contribute CL production in stationary phase. It was necessary to delete all three genes to diminish the CL content of the envelope. Despite being devoid of CL molecules, ΔclsABC mutants were highly virulent during oral and systemic infection for C57BL/6J mice. In macrophages, ΔclsA, ΔclsB, ΔclsC, and ΔclsAC mutants behaved like the wild type, whereas ΔclsAB, ΔclsBC, and ΔclsABC mutants were attenuated and elicited reduced amounts of secreted interleukin-1 beta (IL-1β), IL-18, and lactate dehydrogenase. Hence, when clsA and clsC are deleted, clsB is necessary and sufficient to promote intracellular survival and inflammasome activation. Similarly, when clsB is deleted, clsA and clsC are necessary and sufficient. Therefore, the three CL synthase genes cooperatively and redundantly influence S. Typhimurium inflammasome activation and intracellular survival in C57BL/6J mouse macrophages but are dispensable for virulence in mice.

IMPORTANCESalmonella enterica serovar Typhimurium is a pathogenic Gram-negative bacterium that regulates the cardiolipin (CL) and lipopolysaccharide (LPS) composition of the outer membrane (OM) during infection. Mitochondrial CL molecules activate the inflammasome and its effector caspase-1, which initiates an inflammatory process called pyroptosis. Purified bacterial CL molecules also influence LPS activation of Toll-like receptor 4 (Tlr4). S. Typhimurium resides within macrophage vacuoles and activates Tlr4 and the inflammasome during infection. However, the contribution of the three bacterial CL synthase genes (cls) to microbial pathogenesis and inflammation had not been tested. This study supports that the genes encoding the CL synthases work coordinately to promote intracellular survival in macrophages and to activate the inflammasome but do not influence inflammatory cytokine production downstream of Tlr4 or virulence in C57BL/6J mice. The macrophage phenotypes are not directly attributable to CL production but are caused by deleting specific combinations of cls gene products.

A Cholesterol Dimer Stabilizes the Inactivated State of an Inward-Rectifier Potassium Channel

Cholesterol oligomers reside in multiple membrane protein X-ray crystal structures. Yet, there is no direct link between these oligomers and a biological function. Here we present the structural and functional details of a cholesterol dimer that stabilizes the inactivated state of an inward-rectifier potassium channel KirBac1.1. K⁺ efflux assays confirm that high cholesterol concentration reduces K⁺ conductance. We then determine the structure of the cholesterol-KirBac1.1 complex using Xplor-NIH simulated annealing calculations driven by solid-state NMR distance measurements. These calculations identified an α-α cholesterol dimer docked to a cleft formed by adjacent subunits of the homotetrameric protein. We compare these results to coarse grain molecular dynamics simulations. This is one of the first examples of a cholesterol oligomer performing a distinct biological function and structural characterization of a conserved promiscuous lipid binding region.

What have molecular simulations contributed to understanding of Gram-negative bacterial cell envelopes?

Bacterial cell envelopes are compositionally complex and crowded and while highly dynamic in some areas, their molecular motion is very limited, to the point of being almost static in others. Therefore, it is no real surprise that studying them at high resolution across a range of temporal and spatial scales requires a number of different techniques. Details at atomistic to molecular scales for up to tens of microseconds are now within range for molecular dynamics simulations. Here we review how such simulations have contributed to our current understanding of the cell envelopes of Gram-negative bacteria.

PyLipID: A Python Package for Analysis of Protein-Lipid Interactions from Molecular Dynamics Simulations

Lipids play important modulatory and structural roles for membrane proteins. Molecular dynamics simulations are frequently used to provide insights into the nature of these protein-lipid interactions. Systematic comparative analysis requires tools that provide algorithms for objective assessment of such interactions. We introduce PyLipID, a Python package for the identification and characterization of specific lipid interactions and binding sites on membrane proteins from molecular dynamics simulations. PyLipID uses a community analysis approach for binding site detection, calculating lipid residence times for both the individual protein residues and the detected binding sites. To assist structural analysis, PyLipID produces representative bound lipid poses from simulation data, using a density-based scoring function. To estimate residue contacts robustly, PyLipID uses a dual-cutoff scheme to differentiate between lipid conformational rearrangements while bound from full dissociation events. In addition to the characterization of protein-lipid interactions, PyLipID is applicable to analysis of the interactions of membrane proteins with other ligands. By combining automated analysis, efficient algorithms, and open-source distribution, PyLipID facilitates the systematic analysis of lipid interactions from large simulation data sets of multiple species of membrane proteins.

Leaflet Asymmetry Modeling in the Lipid Composition of Escherichia coli Cytoplasmic Membranes

Lipid composition asymmetry between leaflets is important to cell function and plays a key role in the "positive inside" rule in transmembrane proteins. In this work, Escherichia coli inner plasma membrane models reflecting this asymmetry have been investigated at the early-log and stationary stages during the bacterial lifecycle using all-atom molecular dynamics simulations. The CHARMM36 lipid force field is used, and selected membrane properties are tested for variations between two leaflets and whole membranes. Our models include bacterial lipids with a cyclopropane moiety on the sn-2 acyl chain in the stationary membrane model. The PE/PG ratio for two leaflets reflects the "positive inside" rule of membrane proteins, set to 6.8 and 2.8 for the inner and outer leaflets of the two models, respectively. We are the first to model leaflet asymmetry in the lipid composition of E. coli cytoplasmic membranes and observe the effect on membrane properties in leaflets and whole membranes. Specifically, our results show that for the stationary phase bilayer, the surface area per lipid (SA/lipid) is larger, the thickness (2D_C and D_B) is smaller, the tilt angle is larger, the tilt modulus is smaller, and the deuterium order parameters (S_CD) of sn-1 and sn-2 tails are lower, compared to the early-log stage. Moreover, the stationary stage bilayer has a positive spontaneous curvature, while the early-log stage has a near flat spontaneous curvature. For leaflet asymmetry, the inner leaflet has a larger SA/lipid, a smaller thickness, a smaller elastic tilt modulus (a larger tilt angle), and lower S_CD, compared to the outer leaflet in both stages. Moreover, an asymmetric membrane involves a lipid tilt and a lateral extension, varying from a reference state of a pre-equilibrium membrane. This work encourages a more profound exploration of leaflet asymmetry in various other membrane models and how this might affect the structure and function of membrane-associated peptides and proteins.

Relative Affinities of Protein-Cholesterol Interactions from Equilibrium Molecular Dynamics Simulations

Specific interactions of lipids with membrane proteins contribute to protein stability and function. Multiple lipid interactions surrounding a membrane protein are often identified in molecular dynamics (MD) simulations and are, increasingly, resolved in cryo-electron microscopy (cryo-EM) densities. Determining the relative importance of specific interaction sites is aided by determination of lipid binding affinities using experimental or simulation methods. Here, we develop a method for determining protein-lipid binding affinities from equilibrium coarse-grained MD simulations using binding saturation curves, designed to mimic experimental protocols. We apply this method to directly obtain affinities for cholesterol binding to multiple sites on a range of membrane proteins and compare our results with free energies obtained from density-based equilibrium methods and with potential of mean force calculations, getting good agreement with respect to the ranking of affinities for different sites. Thus, our binding saturation method provides a robust, high-throughput alternative for determining the relative consequence of individual sites seen in, e.g., cryo-EM derived membrane protein structures surrounded by an array of ancillary lipid densities.

Deciphering ion transport and ATPase coupling in the intersubunit tunnel of KdpFABC

KdpFABC, a high-affinity K⁺ pump, combines the ion channel KdpA and the P-type ATPase KdpB to secure survival at K⁺ limitation. Here, we apply a combination of cryo-EM, biochemical assays, and MD simulations to illuminate the mechanisms underlying transport and the coupling to ATP hydrolysis. We show that ions are transported via an intersubunit tunnel through KdpA and KdpB. At the subunit interface, the tunnel is constricted by a phenylalanine, which, by polarized cation-π stacking, controls K⁺ entry into the canonical substrate binding site (CBS) of KdpB. Within the CBS, ATPase coupling is mediated by the charge distribution between an aspartate and a lysine. Interestingly, individual elements of the ion translocation mechanism of KdpFABC identified here are conserved among a wide variety of P-type ATPases from different families. This leads us to the hypothesis that KdpB might represent an early descendant of a common ancestor of cation pumps.