Recent Advances in Modeling Membrane ß-Barrel Proteins Using Molecular Dynamics Simulations: From Their Lipid Environments to Their Assemblies.

Spurred by advances in AI-driven modeling and experimental methods, molecular dynamics simulations are now acting as a platform to integrate these different approaches. This combination of methods is especially useful to understand ß-barrel proteins from the molecular level, e.g., identifying specific interactions with lipids or small molecules, up to assemblies comprised of hundreds of proteins and thousands of lipids. In this minireview, we will discuss recent advances, mainly from the last 5 years, in modeling ß-barrel proteins and their assemblies. These approaches require specific kinds of modeling and potentially different model resolutions that we will first describe in Subheading 1. We will then focus on different aspects of ß-barrel protein modeling: how different types of molecules can diffuse through ß-barrel proteins (Subheading 2); how lipids can interact with these proteins (Subheading 3); how ß-barrel proteins can interact with membrane partners (Subheading 4) or periplasmic extensions and partners (Subheading 5) to form large assemblies.

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.

Direct determination of oligomeric organization of integral membrane proteins and lipids from intact customizable bilayer.

Hierarchical organization of integral membrane proteins (IMP) and lipids at the membrane is essential for regulating myriad downstream signaling. A quantitative understanding of these processes requires both detections of oligomeric organization of IMPs and lipids directly from intact membranes and determination of key membrane components and properties that regulate them. Addressing this, we have developed a platform that enables native mass spectrometry (nMS) analysis of IMP-lipid complexes directly from intact and customizable lipid membranes. Both the lipid composition and membrane properties (such as curvature, tension, and fluidity) of these bilayers can be precisely customized to a target membrane. Subsequent direct nMS analysis of these intact proteolipid vesicles can yield the oligomeric states of the embedded IMPs, identify bound lipids, and determine the membrane properties that can regulate the observed IMP-lipid organization. Applying this method, we show how lipid binding regulates neurotransmitter release and how membrane composition regulates the functional oligomeric state of a transporter.

Mesoscale simulations: An indispensable approach to understand biomembranes.

Computer simulation techniques form a versatile tool, a computational microscope, for exploring biological processes. This tool has been particularly effective in exploring different features of biological membranes. In recent years, thanks to elegant multiscale simulation schemes, some fundamental limitations of investigations by distinct simulation techniques have been resolved. As a result, we are now capable of exploring processes spanning multiple scales beyond the capacity of any single technique. In this perspective, we argue that mesoscale simulations require more attention and must be further developed to fill evident gaps in a quest toward simulating and modeling living cell membranes.

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 CIII2CIV 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.

The guidance and adhesion protein FLRT2 dimerizes in cis via dual small-X3-small transmembrane motifs.

Fibronectin Leucine-rich Repeat Transmembrane (FLRT 1-3) proteins are a family of broadly expressed single-spanning transmembrane receptors that play key roles in development. Their extracellular domains mediate homotypic cell-cell adhesion and heterotypic protein interactions with other receptors to regulate cell adhesion and guidance. These in trans FLRT interactions determine the formation of signaling complexes of varying complexity and function. Whether FLRTs also interact at the surface of the same cell, in cis, remains unknown. Here, molecular dynamics simulations reveal two dimerization motifs in the FLRT2 transmembrane helix. Single particle tracking experiments show that these Small-X3-Small motifs synergize with a third dimerization motif encoded in the extracellular domain to permit the cis association and co-diffusion patterns of FLRT2 receptors on cells. These results may point to a competitive switching mechanism between in cis and in trans interactions, which suggests that homotypic FLRT interaction mirrors the functionalities of classic adhesion molecules.

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.

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

Integral membrane proteins are localized and/or regulated by lipids present in the surrounding bilayer. While bacteria have relatively simple membranes, there is ample evidence that many bacterial proteins bind to specific lipids, especially the anionic lipid cardiolipin. Here, we apply molecular dynamics simulations to assess lipid binding to 42 different Escherichia coli inner membrane proteins. Our data reveal an asymmetry between the membrane leaflets, with increased anionic lipid binding to the inner leaflet regions of the proteins, particularly for cardiolipin. From our simulations, we identify >700 independent cardiolipin binding sites, allowing us to identify the molecular basis of a prototypical cardiolipin binding site, which we validate against structures of bacterial proteins bound to cardiolipin. This allows us to construct a set of metrics for defining a high-affinity cardiolipin binding site on bacterial membrane proteins, paving the way for a heuristic approach to defining other protein-lipid interactions.

Modulation of adenosine A2a receptor oligomerization by receptor activation and PIP2 interactions.

GPCRs have been shown to form oligomers, which generate distinctive signaling outcomes. However, the structural nature of the oligomerization process remains uncertain. We have characterized oligomeric configurations of the adenosine A2a receptor (A2aR) by combining large-scale molecular dynamics simulations with Markov state models. These oligomeric structures may also serve as templates for studying oligomerization of other class A GPCRs. Our simulation data revealed that receptor activation results in enhanced oligomerization, more diverse oligomer populations, and a more connected oligomerization network. The active state conformation of the A2aR shifts protein-protein association interfaces to those involving intracellular loop ICL3 and transmembrane helix TM6. Binding of PIP2 to A2aR stabilizes protein-protein interactions via PIP2-mediated association interfaces. These results indicate that A2aR oligomerization is responsive to the local membrane lipid environment. This, in turn, suggests a modulatory effect on A2aR whereby a given oligomerization profile favors the dynamic formation of specific supramolecular signaling complexes.

Origin and diversification of the cardiolipin biosynthetic pathway in the Eukarya domain.

Cardiolipin (CL) and its precursor phosphatidylglycerol (PG) are important anionic phospholipids widely distributed throughout all domains of life. They have key roles in several cellular processes by shaping membranes and modulating the activity of the proteins inserted into those membranes. They are synthesized by two main pathways, the so-called eukaryotic pathway, exclusively found in mitochondria, and the prokaryotic pathway, present in most bacteria and archaea. In the prokaryotic pathway, the first and the third reactions are catalyzed by phosphatidylglycerol phosphate synthase (Pgps) belonging to the transferase family and cardiolipin synthase (Cls) belonging to the hydrolase family, while in the eukaryotic pathway, those same reactions are catalyzed by unrelated homonymous enzymes: Pgps of the hydrolase family and Cls of the transferase family. Because of the enzymatic arrangement found in both pathways, it seems that the eukaryotic pathway evolved by convergence to the prokaryotic pathway. However, since mitochondria evolved from a bacterial endosymbiont, it would suggest that the eukaryotic pathway arose from the prokaryotic pathway. In this review, it is proposed that the eukaryote pathway evolved directly from a prokaryotic pathway by the neofunctionalization of the bacterial enzymes. Moreover, after the eukaryotic radiation, this pathway was reshaped by horizontal gene transfers or subsequent endosymbiotic processes.

Monolysocardiolipin (MLCL) interactions with mitochondrial membrane proteins.

Monolysocardiolipin (MLCL) is a three-tailed variant of cardiolipin (CL), the signature lipid of mitochondria. MLCL is not normally found in healthy tissue but accumulates in mitochondria of people with Barth syndrome (BTHS), with an overall increase in the MLCL:CL ratio. The reason for MLCL accumulation remains to be fully understood. The effect of MLCL build-up and decreased CL content in causing the characteristics of BTHS are also unclear. In both cases, an understanding of the nature of MLCL interaction with mitochondrial proteins will be key. Recent work has shown that MLCL associates less tightly than CL with proteins in the mitochondrial inner membrane, suggesting that MLCL accumulation is a result of CL degradation, and that the lack of MLCL-protein interactions compromises the stability of the protein-dense mitochondrial inner membrane, leading to a decrease in optimal respiration. There is some data on MLCL-protein interactions for proteins involved in the respiratory chain and in apoptosis, but there remains much to be understood regarding the nature of MLCL-protein interactions. Recent developments in structural, analytical and computational approaches mean that these investigations are now possible. Such an understanding will be key to further insights into how MLCL accumulation impacts mitochondrial membranes. In turn, these insights will help to support the development of therapies for people with BTHS and give a broader understanding of other diseases involving defective CL content.

Computational Investigation of Voltage-Gated Sodium Channel ß3 Subunit Dynamics.

Voltage-gated sodium (Na v ) channels form the basis for the initiation of the action potential in excitable cells by allowing sodium ions to pass through the cell membrane. The Na v channel α subunit is known to function both with and without associated ß subunits. There is increasing evidence that these ß subunits have multiple roles that include not only influencing the voltage-dependent gating but also the ability to alter the spatial distribution of the pore-forming α subunit. Recent structural data has shown possible ways in which ß1 subunits may interact with the α subunit. However, the position of the ß1 subunit would not be compatible with a previous trimer structure of the ß3 subunit. Furthermore, little is currently known about the dynamic behavior of the ß subunits both as individual monomers and as higher order oligomers. Here, we use multiscale molecular dynamics simulations to assess the dynamics of the ß3, and the closely related, ß1 subunit. These findings reveal the spatio-temporal dynamics of ß subunits and should provide a useful framework for interpreting future low-resolution experiments such as atomic force microscopy.

Defining how multiple lipid species interact with inward rectifier potassium (Kir2) channels.

Protein-lipid interactions are a key element of the function of many integral membrane proteins. These potential interactions should be considered alongside the complexity and diversity of membrane lipid composition. Inward rectifier potassium channel (Kir) Kir2.2 has multiple interactions with plasma membrane lipids: Phosphatidylinositol (4, 5)-bisphosphate (PIP2) activates the channel; a secondary anionic lipid site has been identified, which augments the activation by PIP2; and cholesterol inhibits the channel. Molecular dynamics simulations are used to characterize in molecular detail the protein-lipid interactions of Kir2.2 in a model of the complex plasma membrane. Kir2.2 has been simulated with multiple, functionally important lipid species. From our simulations we show that PIP2 interacts most tightly at the crystallographic interaction sites, outcompeting other lipid species at this site. Phosphatidylserine (PS) interacts at the previously identified secondary anionic lipid interaction site, in a PIP2 concentration-dependent manner. There is interplay between these anionic lipids: PS interactions are diminished when PIP2 is not present in the membrane, underlining the need to consider multiple lipid species when investigating protein-lipid interactions.

Lipid-Dependent Regulation of Ion Channels and G Protein-Coupled Receptors: Insights from Structures and Simulations.

Ion channels and G protein-coupled receptors (GPCRs) are regulated by lipids in their membrane environment. Structural studies combined with biophysical and molecular simulation investigations reveal interaction sites for specific lipids on membrane protein structures. For K channels, PIP2 plays a key role in regulating Kv and Kir channels. Likewise, several recent cryo-EM structures of TRP channels have revealed bound lipids, including PIP2 and cholesterol. Among the pentameric ligand-gated ion channel family, structural and biophysical studies suggest the M4 TM helix may act as a lipid sensor, e.g., forming part of the binding sites for neurosteroids on the GABAA receptor. Structures of GPCRs have revealed multiple cholesterol sites, which may modulate both receptor dynamics and receptor oligomerization. PIP2 also interacts with GPCRs and may modulate their interactions with G proteins. Overall, it is evident that multiple lipid binding sites exist on channels and receptors that modulate their function allosterically and are potential druggable sites.

More Favorable Palmitic Acid Over Palmitoleic Acid Modification of Wnt3 Ensures Its Localization and Activity in Plasma Membrane Domains.

While the lateral organization of plasma membrane components has been shown to control binding of Wnt ligands to their receptors preferentially in the ordered membrane domains, the role of posttranslational lipid modification of Wnt on this selective binding is unknown. Here, we identify that the canonical Wnt is presumably acylated by palmitic acid, a saturated 16-carbon fatty acid, at a conserved serine residue. Acylation of Wnt3 is dispensable for its secretion and binding to Fz8 while it is essential for Wnt3's proper binding and domain-like diffusion in the ordered membrane domains. We further unravel that non-palmitoylated Wnt3 is unable to activate Wnt/ß-catenin signaling either in zebrafish embryos or in mammalian cells. Based on these results, we propose that the lipidation of canonical Wnt, presumably by a saturated fatty acid, determines its competence in interacting with the receptors in the appropriate domains of the plasma membrane, ultimately keeping the signaling activity under control.

How nanoscale protein interactions determine the mesoscale dynamic organisation of bacterial outer membrane proteins.

The spatiotemporal organisation of membranes is often characterised by the formation of large protein clusters. In Escherichia coli, outer membrane protein (OMP) clustering leads to OMP islands, the formation of which underpins OMP turnover and drives organisation across the cell envelope. Modelling how OMP islands form in order to understand their origin and outer membrane behaviour has been confounded by the inherent difficulties of simulating large numbers of OMPs over meaningful timescales. Here, we overcome these problems by training a mesoscale model incorporating thousands of OMPs on coarse-grained molecular dynamics simulations. We achieve simulations over timescales that allow direct comparison to experimental data of OMP behaviour. We show that specific interaction surfaces between OMPs are key to the formation of OMP clusters, that OMP clusters present a mesh of moving barriers that confine newly inserted proteins within islands, and that mesoscale simulations recapitulate the restricted diffusion characteristics of OMPs.

Interactions of the EphA2 Kinase Domain with PIPs in Membranes: Implications for Receptor Function.

EphA2 is a member of the receptor tyrosine kinase family. Interactions of the cytoplasmic region of EphA2 with the cell membrane are functionally important and yet remain incompletely characterized. Molecular dynamics simulations combined with biochemical studies reveal the interactions of the transmembrane, juxtamembrane (JM), and kinase domains with the membrane. We describe how the kinase domain is oriented relative to the membrane and how the JM region can modulate this interaction. We highlight the role of phosphatidylinositol phosphates (PIPs) in mediating the interaction of the kinase domain with the membrane and, conversely, how positively charged patches at the kinase surface and in the JM region induce the formation of nanoclusters of PIP molecules in the membrane. Integration of these results with those from previous studies enable computational reconstitution of a near complete EphA2 receptor within a membrane, suggesting a role for receptor-lipid interactions in modulation of EphA2.

Cardiolipin dynamics and binding to conserved residues in the mitochondrial ADP/ATP carrier.

Cardiolipin in eukaryotes is found in the mitochondrial inner membrane, where it interacts with membrane proteins and, although not essential, is necessary for the optimal activity of a number of proteins. One of them is the mitochondrial ADP/ATP carrier, which imports ADP into the mitochondrion and exports ATP. In the crystal structures, cardiolipin is bound to three equivalent sites of the ADP/ATP carrier, but its role is unresolved. Conservation of residues at these cardiolipin binding sites across other members of the mitochondrial carrier superfamily indicates cardiolipin binding is likely to be important for the function of all mitochondrial carriers. Multiscale simulations were performed in a cardiolipin-containing membrane to investigate the dynamics of cardiolipin around the yeast and bovine ADP/ATP carriers in a lipid bilayer and the properties of the cardiolipin-binding sites. In coarse-grain simulations, cardiolipin molecules bound to the carriers for longer periods of time than phosphatidylcholine and phosphatidylethanolamine lipids-with timescales in the tens of microseconds. Three long-lived cardiolipin binding sites overlapped with those in the crystal structures of the carriers. Other shorter-lived cardiolipin interaction sites were identified in both membrane leaflets. However, the timescales of the interactions were of the same order as phosphatidylcholine and phosphatidylethanolamine, suggesting that these sites are not specific for cardiolipin binding. The calculation of lipid binding times and the overlap of the cardiolipin binding sites between the structures and simulations demonstrate the potential of multiscale simulations to investigate the dynamics and behavior of lipids interacting with membrane proteins.

Protein crowding and lipid complexity influence the nanoscale dynamic organization of ion channels in cell membranes.

Cell membranes are crowded and complex environments. To investigate the effect of protein-lipid interactions on dynamic organization in mammalian cell membranes, we have performed coarse-grained molecular dynamics simulations containing >100 copies of an inwardly rectifying potassium (Kir) channel which forms specific interactions with the regulatory lipid phosphatidylinositol 4,5-bisphosphate (PIP2). The tendency of protein molecules to cluster has the effect of organizing the membrane into dynamic compartments. At the same time, the diversity of lipids present has a marked effect on the clustering behavior of ion channels. Sub-diffusion of proteins and lipids is observed. Protein crowding alters the sub-diffusive behavior of proteins and lipids such as PIP2 which interact tightly with Kir channels. Protein crowding also affects bilayer properties, such as membrane undulations and bending rigidity, in a PIP2-dependent manner. This interplay between the diffusion and the dynamic organization of Kir channels may have important implications for channel function.