Emerging Diversity in Lipid-Protein Interactions

Biomimetic nanoplatforms constructed from dialkylaminostyryl hetarene dyes and phospholipids exhibiting selective fluorescent response to specific proteins

The present work explores the specificity of supramolecular assemblies comprising dialkylaminostyrylhetarene dye molecules incorporated into phosphatidylcholine (PC) or phosphatidylserine (PS) aggregates. In PS-based assemblies, the dyes demonstrate a concentration-dependent fluorescent response, distinguishing anionic proteins such as bovine serum albumin (BSA) and pepsin from lysozyme (LYZ) in aqueous solutions. Conversely, no significant response is observed when the dyes are incorporated into the well-organized bilayers of neutral PC. The fluorescent response arises from the binding of dyes to proteins, leading to the detachment of dye molecules from the assemblies, rather than from the binding of proteins to the assemblies, although the latter process is facilitated by electrostatic attraction. Thus, both the poor ordering of PS molecules and the interfacial arrangement of the dyes are prerequisites for the fluorescent response of dye-PS aggregates. The structure of the dyes significantly impacts the spectral features of dye-PS and dye-protein assemblies. An optimal dye structure has been identified for the recognition of BSA, with a limit of detection (LOD) of 10.8 nM.

Methods for monitoring protein-membrane binding. Comparison based on the interactions between amyloidogenic protein human cystatin C and phospholipid liposomes

A cell membrane is an essential cellular component providing protection against the outer environment. It is also a host for proteins and carbohydrates responsible for, e.g. transporter, receptor, or enzymatic functions. In parallel, the membrane may also be implicated in pathological processes leading, e.g. to the oligomerization of amyloid-forming proteins, a hallmark of i.a. Alzheimer's disease. The increasing need for detailed information on mechanisms driving the amyloid formation and the potential role of cell membranes in the process proves the research on protein-membrane interactions biologically relevant. Considering the potential and limitations of the relatively well established and newly developed methods, this study focused on selecting methods that allow a broad and comprehensive description of interactions between amyloidogenic protein human cystatin C and lipid bilayers. In the first step, dot-blot and ELISA tests were selected as techniques allowing fast screening for protein-ligand interactions. Next, surface plasmon resonance, spectral shift, biolayer interferometry, and switchSENSE® technology were used to determine kinetic parameters and binding constants for interactions between human cystatin C and the selected lipid bilayers. Based on the obtained results we have proposed the most promising candidates for monitoring of interactions and determining affinity between amyloidogenic proteins and membrane mimetics.

Bilayer lipids modulate ligand binding to atypical chemokine receptor 3

Chemokine receptors belong to the large class of G protein-coupled receptors (GPCRs) and are involved in a number of (patho)physiological processes. Previous studies highlighted the importance of membrane lipids for modulating GPCR structure and function. However, the underlying mechanisms of how lipids regulate GPCRs are often poorly understood. Here, we report that anionic lipid bilayers increase the binding affinity of the chemokine CXCL12 for the atypical chemokine receptor 3 (ACKR3) by modulating the CXCL12 binding kinetics. Notably, the anionic bilayer favors CXCL12 over the more positively charged chemokine CXCL11, which we explained by bilayer interactions orienting CXCL12 but not CXCL11 for productive ACKR3 binding. Furthermore, our data suggest a stabilization of active ACKR3 conformations in anionic bilayers. Taken together, the described regulation of chemokine selectivity of ACKR3 by the lipid bilayer proposes an extended version of the classical model of chemokine binding including the lipid environment of the receptor.

Glycosphingolipid synthesis mediates immune evasion in KRAS-driven cancer

Cancer cells frequently alter their lipids to grow and adapt to their environment1-3. Despite the critical functions of lipid metabolism in membrane physiology, signalling and energy production, how specific lipids contribute to tumorigenesis remains incompletely understood. Here, using functional genomics and lipidomic approaches, we identified de novo sphingolipid synthesis as an essential pathway for cancer immune evasion. Synthesis of sphingolipids is surprisingly dispensable for cancer cell proliferation in culture or in immunodeficient mice but required for tumour growth in multiple syngeneic models. Blocking sphingolipid production in cancer cells enhances the anti-proliferative effects of natural killer and CD8+ T cells partly via interferon-γ (IFNγ) signalling. Mechanistically, depletion of glycosphingolipids increases surface levels of IFNγ receptor subunit 1 (IFNGR1), which mediates IFNγ-induced growth arrest and pro-inflammatory signalling. Finally, pharmacological inhibition of glycosphingolipid synthesis synergizes with checkpoint blockade therapy to enhance anti-tumour immune response. Altogether, our work identifies glycosphingolipids as necessary and limiting metabolites for cancer immune evasion.

Depletion of membrane cholesterol modifies structure, dynamic and activation of Nav1.7

Cholesterol is a major component of plasma membranes and plays a significant role in actively regulating the functioning of several membrane proteins in humans. In this study, we focus on the role of cholesterol depletion on the voltage-gated sodium channel Nav1.7, which is primarily expressed in the peripheral sensory neurons and linked to various chronic inherited pain syndromes. Coarse-grained molecular dynamics simulations revealed key dynamic changes of Nav1.7 upon membrane cholesterol depletion: A loss of rigidity in the structural motifs linked to activation and fast-inactivation is observed, suggesting an easier transition of the channel between different gating states. In-vitro whole-cell patch clamp experiments on HEK293t cells expressing Nav1.7 validated these predictions at the functional level: Hyperpolarizing shifts in the voltage-dependence of activation and fast-inactivation were observed along with an acceleration of the time to peak and onset kinetics of fast inactivation. These results underline the critical role of membrane composition, and of cholesterol in particular, in influencing Nav1.7 gating characteristics. Furthermore, our results also point to cholesterol-driven changes of the geometry of drug-binding regions, hinting to a key role of the membrane environment in the regulation of drug effects.

Identification of lipid-specific proteins with high-density lipid-immobilized beads

In biological membranes, lipids often interact with membrane proteins (MPs), regulating the localization and activity of MPs in cells. Although elucidating lipid-MP interactions is critical to comprehend the physiological roles of lipids, a systematic and comprehensive identification of lipid-binding proteins has not been adequately established. Therefore, we report the development of lipid-immobilized beads where lipid molecules were covalently immobilized. Owing to the detergent tolerance, these beads enable screening of water-soluble proteins and MPs, the latter of which typically necessitate surfactants for solubilization. Herein, two sphingolipid species-ceramide and sphingomyelin-which are major constituents of lipid rafts, were immobilized on the beads. We first showed that the density of immobilized lipid molecules on the beads was as high as that of biological lipid membranes. Subsequently, we confirmed that these beads enabled the selective pulldown of known sphingomyelin- or ceramide-binding proteins (lysenin, p24, and CERT) from protein mixtures, including cell lysates. In contrast, commercial sphingomyelin beads, on which lipid molecules are sparsely immobilized through biotin-streptavidin linkage, failed to capture lysenin, a well-known protein that recognizes clustered sphingomyelin molecules. This clearly demonstrates the applicability of our beads for obtaining proteins that recognize not only a single lipid molecule but also lipid clusters or lipid membranes. Finally, we demonstrated the screening of lipid-binding proteins from Neuro2a cell lysates using these beads. This method is expected to significantly contribute to the understanding of interactions between lipids and proteins and to unravel the complexities of lipid diversity.

Emerging Designs and Applications for Biomembrane Biosensors

Nature has inspired the development of biomimetic membrane sensors in which the functionalities of biological molecules, such as proteins and lipids, are harnessed for sensing applications. This review provides an overview of the recent developments for biomembrane sensors compatible with either bulk or planar sensing applications, namely using lipid vesicles or supported lipid bilayers, respectively. We first describe the individual components required for these sensing platforms and the design principles that are considered when constructing them, and we segue into recent applications being implemented across multiple fields. Our goal for this review is to illustrate the versatility of nature's biomembrane toolbox and simultaneously highlight how biosensor platforms can be enhanced by harnessing it.

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.

Curvature Footprints of Transmembrane Proteins in Simulations with the Martini Force Field

Membranes play essential roles in biological systems and are tremendously diverse in the topologies and chemical and elastic properties that define their functions. In many cases, a given membrane may display considerable heterogeneity, with localized clusters of lipids and proteins exhibiting distinct characteristics compared to adjoining regions. These lipid-protein assemblies can span nanometers to micrometers and are associated with cellular processes such as transport and signaling. While lipid-protein assemblages are dynamic, they can be stabilized by coupling between local membrane composition and shape. Due to the inherent difficulty in resolving atomistic details of membrane proteins in their native lipid environments, these complexes are notoriously challenging to study experimentally; however, molecular dynamics (MD) simulations might be a viable alternative. Here, we aim to assess the utility of coarse-grained (CG) MD simulations with the Martini force field for studying membrane curvature induced by transmembrane (TM) proteins that are reported to generate local curvature. The direction and magnitude of curvature induced by five different TM proteins, as well as certain lipid-protein and protein-protein interactions, were found to be in good agreement with available reference data.

Simultaneous Native Mass Spectrometry Analysis of Single and Double Mutants To Probe Lipid Binding to Membrane Proteins

Lipids are critical modulators of membrane protein structure and function. However, it is challenging to investigate the thermodynamics of protein-lipid interactions because lipids can simultaneously bind membrane proteins at different sites with different specificities. Here, we developed a native mass spectrometry (MS) approach using single and double mutants to measure the relative energetic contributions of specific residues on Aquaporin Z (AqpZ) toward cardiolipin (CL) binding. We first mutated potential lipid-binding residues on AqpZ, and mixed mutant and wild-type proteins together with CL. By using native MS to simultaneously resolve lipid binding to the mutant and wild-type proteins in a single spectrum, we directly determined the relative affinities of CL binding, thereby revealing the relative Gibbs free energy change for lipid binding caused by the mutation. Comparing different mutants revealed that W14 contributes to the tightest CL binding site, with R224 contributing to a lower affinity site. Using double mutant cycling, we investigated the synergy between W14 and R224 sites on CL binding. Overall, this novel native MS approach provides unique insights into the binding of lipids to specific sites on membrane proteins.

Improved Highly Mobile Membrane Mimetic Model for Investigating Protein-Cholesterol Interactions

Cholesterol (CHL) plays an integral role in modulating the function and activity of various mammalian membrane proteins. Due to the slow dynamics of lipids, conventional computational studies of protein-CHL interactions rely on either long-time scale atomistic simulations or coarse-grained approximations to sample the process. A highly mobile membrane mimetic (HMMM) has been developed to enhance lipid diffusion and thus used to facilitate the investigation of lipid interactions with peripheral membrane proteins and, with customized in silico solvents to replace phospholipid tails, with integral membrane proteins. Here, we report an updated HMMM model that is able to include CHL, a nonphospholipid component of the membrane, henceforth called HMMM-CHL. To this end, we had to optimize the effect of the customized solvents on CHL behavior in the membrane. Furthermore, the new solvent is compatible with simulations using force-based switching protocols. In the HMMM-CHL, both improved CHL dynamics and accelerated lipid diffusion are integrated. To test the updated model, we have applied it to the characterization of protein-CHL interactions in two membrane protein systems, the human β2-adrenergic receptor (β2AR) and the mitochondrial voltage-dependent anion channel 1 (VDAC-1). Our HMMM-CHL simulations successfully identified CHL binding sites and captured detailed CHL interactions in excellent consistency with experimental data as well as other simulation results, indicating the utility of the improved model in applications where an enhanced sampling of protein-CHL interactions is desired.

Molecular dynamics simulations of lipid-protein interactions in SLC4 proteins

The SLC4 family of secondary bicarbonate transporters is responsible for the transport of HCO3-, CO32-, Cl-, Na+, K+, NH3, and H+, which are necessary for regulation of pH and ion homeostasis. They are widely expressed in numerous tissues throughout the body and function in different cell types with different membrane properties. Potential lipid roles in SLC4 function have been reported in experimental studies, focusing mostly on two members of the family: AE1 (Cl-/HCO3- exchanger) and NBCe1 (Na+-CO32-cotransporter). Previous computational studies of the outward-facing state of AE1 with model lipid membranes revealed enhanced protein-lipid interactions between cholesterol (CHOL) and phosphatidylinositol bisphosphate (PIP2). However, the protein-lipid interactions in other members of the family and other conformation states are still poorly understood and this precludes the detailed studies of a potential regulatory role for lipids in the SLC4 family. In this work, we performed coarse-grained and atomistic molecular dynamics simulations on three members of the SLC4 family with different transport modes: AE1, NBCe1, and NDCBE (an Na+-CO32-/Cl- exchanger), in model HEK293 membranes consisting of CHOL, PIP2, phosphatidylcholine, phosphatidylethanolamine, phosphatidylserine, and sphingomyelin. The recently resolved inward-facing state of AE1 was also included in the simulations. Lipid-protein contact analysis of the simulated trajectories was performed with the ProLint server, which provides a multitude of visualization tools for illustration of areas of enhanced lipid-protein contact and identification of putative lipid binding sites within the protein matrix. We observed enrichment of CHOL and PIP2 around all proteins with subtle differences in their distribution depending on the protein type and conformation state. Putative binding sites were identified for CHOL, PIP2, phosphatidylcholine, and sphingomyelin in the three studied proteins, and their potential roles in the SLC4 transport function, conformational transition, and protein dimerization are discussed.

Three-Dimensional Interaction Homology: Deconstructing Residue-Residue and Residue-Lipid Interactions in Membrane Proteins

A method is described to deconstruct the network of hydropathic interactions within and between a protein's sidechain and its environment into residue-based three-dimensional maps. These maps encode favorable and unfavorable hydrophobic and polar interactions, in terms of spatial positions for optimal interactions, relative interaction strength, as well as character. In addition, these maps are backbone angle-dependent. After map calculation and clustering, a finite number of unique residue sidechain interaction maps exist for each backbone conformation, with the number related to the residue's size and interaction complexity. Structures for soluble proteins (~749,000 residues) and membrane proteins (~387,000 residues) were analyzed, with the latter group being subdivided into three subsets related to the residue's position in the membrane protein: soluble domain, core-facing transmembrane domain, and lipid-facing transmembrane domain. This work suggests that maps representing residue types and their backbone conformation can be reassembled to optimize the medium-to-high resolution details of a protein structure. In particular, the information encoded in maps constructed from the lipid-facing transmembrane residues appears to paint a clear picture of the protein-lipid interactions that are difficult to obtain experimentally.

The Role of Lipid Intrinsic Curvature in the Droplet Interface Bilayer

Model bilayers are constructed from lipids having different intrinsic curvatures using the droplet interface bilayer (DIB) method, and their static physicochemical properties are determined. Geometrical and tensiometric measurements are used to derive the free energy of formation (ΔF) of a two-droplet DIB relative to a pair of isolated aqueous droplets, each decorated with a phospholipid monolayer. The lipid molecules employed have different headgroup sizes but identical hydrophobic tail structure, and each is characterized by an intrinsic curvature value (c0) that increases in absolute value with decreasing size of headgroup. Mixtures of lipids at different ratios were also investigated. The role of curvature stress on the values of ΔF of the respective lipid bilayers in these model membranes is discussed and is illuminated by the observation of a decrement in ΔF that scales as a near linear function of c02. Overall, the results reveal an association that should prove useful in studies of ion channels and other membrane proteins embedded in model droplet bilayer systems that will impact the understanding of protein function in cellular membranes composed of lipids of high and low curvature.

Advances in Membrane Mimetic Systems for Manipulation and Analysis of Membrane Proteins: Detergents, Polymers, Lipids and Scaffolds

Extracting membrane proteins from the hydrophobic environment of the biological membrane, in a physiologically relevant and stable state, suitable for downstream analysis remains a challenge. The traditional route to membrane protein extraction has been to use detergents and the last 15 years or so have seen a veritable explosion in the development of novel detergents with improved properties, making them more suitable for individual proteins and specific applications. There have also been significant advances in the development of encapsulation of membrane proteins in lipid based nanodiscs, either directly from the native membrane using polymers allowing effective capture of the protein and protein-associated membrane lipids, or via reconstitution of detergent extracted and purified protein into nanodiscs of defined lipid composition. All of these advances have been successfully applied to the study of membrane proteins via a range of techniques and there have been some spectacular membrane protein structures solved. In addition, the first detailed structural and biophysical analyses of membrane proteins retained within a biological membrane have been reported. Here we summarise and review the recent advances with respect to these new agents and systems for membrane protein extraction, reconstitution and analysis.

Allosteric Cholesterol Site in Glycine Receptors Characterized through Molecular Simulations

Glycine receptors are pentameric ligand-gated ion channels that conduct chloride ions across postsynaptic membranes to facilitate fast inhibitory neurotransmission. In addition to gating by the glycine agonist, interactions with lipids and other compounds in the surrounding membrane environment modulate their function, but molecular details of these interactions remain unclear, in particular, for cholesterol. Here, we report coarse-grained simulations in a model neuronal membrane for three zebrafish glycine receptor structures representing apparent resting, open, and desensitized states. We then converted the systems to all-atom models to examine detailed lipid interactions. Cholesterol bound to the receptor at an outer-leaflet intersubunit site, with a preference for the open and desensitized versus resting states, indicating that it can bias receptor function. Finally, we used short atomistic simulations and iterative amino acid perturbations to identify residues that may mediate allosteric gating transitions. Frequent cholesterol contacts in atomistic simulations clustered with residues identified by perturbation analysis and overlapped with mutations influencing channel function and pathology. Cholesterol binding at this site was also observed in a recently reported pig heteromeric glycine receptor. These results indicate state-dependent lipid interactions relevant to allosteric transitions of glycine receptors, including specific amino acid contacts applicable to biophysical modeling and pharmaceutical design.

Concentration-Dependent Effects of Curcumin on Membrane Permeability and Structure

Growing evidence suggests that many bioactive molecules can nonspecifically modulate the physicochemical properties of membranes and influence the action of embedded membrane proteins. This study investigates the interactions of curcumin with protein-free model membranes consisting of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and DOPC with cholesterol (4/1 mol ratio). The focus is on the capability of curcumin to modify membrane barrier properties such as water permeability assayed through the droplet interface bilayer (DIB) model membrane. For pure DOPC, our findings show a concentration-dependent biphasic effect: a reduction in water permeability is observed at low concentrations (up to 2 mol %), whereas at high concentrations of curcumin, water permeability increases. In the presence of cholesterol, we observed an overall reduction in water permeability. A combination of complementary experimental methods, including phase transition parameters studied by differential scanning calorimetry (DSC) and structural properties measured by attenuated total reflectance (ATR)-FTIR, provides a deeper understanding of concentration-dependent interactions of curcumin with DOPC bilayers in the absence and presence of cholesterol. Our experimental findings align with a molecular mechanism of curcumin's interaction with model membranes, wherein its effect is contingent on its concentration. At low concentrations, curcumin binds to the lipid-water interface through hydrogen bonding with the phosphate headgroup, thereby obstructing the transport of water molecules. Conversely, at high concentrations, curcumin permeates the acyl chain region, inducing packing disorders and demonstrating evidence of phase separation. Enhanced knowledge of the impact of curcumin on membranes, which, in turn, can affect protein function, is likely to be beneficial for the successful translation of curcumin into effective medicine.

Sources of biases in the in vitro testing of nanomaterials: the role of the biomolecular corona

The biological fate of nanomaterials (NMs) is driven by specific interactions through which biomolecules, naturally adhering onto their surface, engage with cell membrane receptors and intracellular organelles. The molecular composition of this layer, called the biomolecular corona (BMC), depends on both the physical-chemical features of the NMs and the biological media in which the NMs are dispersed and cells grow. In this work, we demonstrate that the widespread use of 10% fetal bovine serum in an in vitro assay cannot recapitulate the complexity of in vivo systemic administration, with NMs being transported by the blood. For this purpose, we undertook a comparative journey involving proteomics, lipidomics, high throughput multiparametric in vitro screening, and single molecular feature analysis to investigate the molecular details behind this in vivo/in vitro bias. Our work indirectly highlights the need to introduce novel, more physiological-like media closer in composition to human plasma to produce realistic in vitro screening data for NMs. We also aim to set the basis to reduce this in vitro-in vivo mismatch, which currently limits the formulation of NMs for clinical settings.

Protein-Mediated Changes in Membrane Fluidity and Ordering: Insights into the Molecular Mechanism and Implications for Cellular Function

Probing protein-membrane interactions is vital for understanding biological functionality for various applications such as drug development, targeted drug delivery, and creation of functional biomaterials for medical and industrial purposes. In this study, we have investigated interaction of Human Serum Albumin (HSA) with two different lipids, dipalmitoylphosphatidylglycerol (dDPPG) and dipalmitoylphosphatidylcholine (dDPPC), using Vibrational Sum Frequency Generation spectroscopy at different membrane fluidity values. In the liquid-expanded (LE) state of the lipid, HSA (at pH 3.5) deeply intercalated lipid chains through a combination of electrostatic and hydrophobic interactions, which resulted in more ordering of the lipid chains. However, in the liquid-condensed (LC) state, protein intercalation is decreased due to tighter lipid packing. Moreover, our findings revealed distinct differences in HSA's interaction with dDPPG and dDPPC lipids. The interaction with dDPPC remained relatively weak compared to dDPPG. These results shed light on the significance of protein mediated changes in lipid characteristics, which hold considerable implications for understanding membrane protein behavior, lipid-mediated cellular processes, and lipid-based biomaterial design.

Cholesterol Dependence of the Conformational Changes in Metabotropic Glutamate Receptor 1

Metabotropic glutamate receptors (mGluRs) are class C G protein-coupled receptors that function as obligate dimers in regulating neurotransmission and synaptic plasticity in the central nervous system. The mGluR1 subtype has been shown to be modulated by the membrane lipid environment, particularly cholesterol, though the molecular mechanisms remain elusive. In this study, we employed all-atom molecular dynamics simulations to investigate the effects of cholesterol on the conformational dynamics of the mGluR1 seven-transmembrane (7TM) domain in an inactive state model. Simulations were performed with three different cholesterol concentrations (0%, 10%, and 25%) in a palmitoyl-oleoyl phosphatidylcholine (POPC) lipid bilayer system. Our results demonstrate that cholesterol induces conformational changes in the mGluR1 dimer more significantly than in the individual protomers. Notably, cholesterol modulates the dynamics and conformations of the TM1 and TM2 helices at the dimer interface. Interestingly, an intermediate cholesterol concentration of 10% elicits more pronounced conformational changes compared to both cholesterol-depleted (0%) and cholesterol-enriched (25%) systems. Specific electrostatic interaction unique to the 10% cholesterol system further corroborate these conformational differences. Given the high sequence conservation of the 7TM domains across mGluR subtypes, the cholesterol-dependent effects observed in mGluR1 are likely applicable to other members of this receptor family. Our findings provide atomistic insights into how cholesterol modulates the conformational landscape of mGluRs, which could impact their function and signaling mechanisms.

Configurational Sampling of All-Atom Solvated Membranes Using Hybrid Nonequilibrium Molecular Dynamics Monte Carlo Simulations

All-atom simulations are a powerful approach to study the structure and dynamics of biological membranes. However, sampling the atomic configurations of inhomogeneous membranes can be challenging due to the slow lateral diffusion of the various constituents. To address this issue, a hybrid nonequilibrium molecular dynamics Monte Carlo (neMD/MC) simulation method is proposed in which randomly chosen lipid molecules are swapped to generate configurations that are subsequently accepted or rejected according to a Metropolis criterion based on the alchemical work for the attempted swap calculated via a short trajectory. A dual-topology framework constraining the common atoms of the exchanging molecules yields a good acceptance probability using switching trajectories as short as 10 ps. The performance of the hybrid neMD/MC algorithm and its ability to sample the distribution of lipids near a transmembrane helix carrying a net charge are illustrated for a binary mixture of charged and zwitterionic lipids.

State-dependent binding of cholesterol and an anionic lipid to the muscle-type Torpedo nicotinic acetylcholine receptor

The ability of the Torpedo nicotinic acetylcholine receptor (nAChR) to undergo agonist-induced conformational transitions requires the presence of cholesterol and/or anionic lipids. Here we use recently solved structures along with multiscale molecular dynamics simulations to examine lipid binding to the nAChR in bilayers that have defined effects on nAChR function. We examine how phosphatidic acid and cholesterol, lipids that support conformational transitions, individually compete for binding with phosphatidylcholine, a lipid that does not. We also examine how the two lipids work synergistically to stabilize an agonist-responsive nAChR. We identify rapidly exchanging lipid binding sites, including both phospholipid sites with a high affinity for phosphatidic acid and promiscuous cholesterol binding sites in the grooves between adjacent transmembrane α-helices. A high affinity cholesterol site is confirmed in the inner leaflet framed by a key tryptophan residue on the MX α-helix. Our data provide insight into the dynamic nature of lipid-nAChR interactions and set the stage for a detailed understanding of the mechanisms by which lipids facilitate nAChR function at the neuromuscular junction.

Analyzing lipid distributions and curvature in molecular dynamics simulations of complex membranes

We describe methods to analyze lipid distributions and curvature in membranes with complex lipid mixtures and embedded membrane proteins. We discuss issues involved in these analyses, available tools to calculate curvature preferences of lipids and proteins, and focus on tools developed in our group for visual analysis of lipid-protein interactions and the analysis of membrane curvature.

The density-threshold affinity: Calculating lipid binding affinities from unbiased coarse-grained molecular dynamics simulations

Many membrane proteins are sensitive to their local lipid environment. As structural methods for membrane proteins have improved, there is growing evidence of direct, specific binding of lipids to protein surfaces. Unfortunately the workhorse of understanding protein-small molecule interactions, the binding affinity for a given site, is experimentally inaccessible for these systems. Coarse-grained molecular dynamics simulations can be used to bridge this gap, and are relatively straightforward to learn. Such simulations allow users to observe spontaneous binding of lipids to membrane proteins and quantify localized densities of individual lipids or lipid fragments. In this chapter we outline a protocol for extracting binding affinities from these localized distributions, known as the "density threshold affinity." The density threshold affinity uses an adaptive and flexible definition of site occupancy that alleviates the need to distinguish between "bound'' lipids and bulk lipids that are simply diffusing through the site. Furthermore, the method allows "bead-level" resolution that is suitable for the case where lipids share binding sites, and circumvents ambiguities about a relevant reference state. This approach provides a convenient and straightforward method for comparing affinities of a single lipid species for multiple sites, multiple lipids for a single site, and/or a single lipid species modeled using multiple forcefields.

Gold nanoparticle-powered screening of membrane protein-specific lipids from complex lipid mixtures

Membrane proteins (MPs) are affected by binding of specific lipids. We previously developed a methodology for systematically analyzing MP-lipid interactions leveraging surface plasmon resonance (SPR). In this method, the gold sensor chip surface was modified with a self-assembled monolayer (SAM), which allowed for a larger amount of MP-immobilization. However, the laborious lipid purification step remained a bottleneck. To address this issue, a new strategy has been developed utilizing gold nanoparticles (AuNPs) instead of the gold sensor chip. AuNPs were coated with SAM, on which MP was covalently anchored. The MP-immobilized AuNPs were mixed with a lipid mixture, and the recovered lipids were quantified by LC-MS. Bacteriorhodopsin (bR) was used as an MP to demonstrate this concept. We optimized immobilization conditions and confirmed the efficient immobilization of bR by dynamic light scattering and electron micrographs. Washing conditions for pulldown experiments were optimized to efficiently remove non-specific lipids. A new binding index was introduced to qualitatively reproduce the known affinity of lipids for bR. Consequently, the low-abundant and least-studied lipid S-TeGD was identified as a candidate for bR-specific lipids. This technique can skip the laborious lipid purification process, accelerating the screening of MP-specific lipids from complex lipid mixtures.

GPCRs in the round: SMA-like copolymers and SMALPs as a platform for investigating GPCRs

G-protein-coupled receptors (GPCRs) are the largest family of membrane proteins, regulate a plethora of physiological responses and are the therapeutic target for 30-40% of clinically-prescribed drugs. They are integral membrane proteins deeply embedded in the plasma membrane where they activate intracellular signalling via coupling to G-proteins and β-arrestin. GPCRs are in intimate association with the bilayer lipids and that lipid environment regulates the signalling functions of GPCRs. This complex lipid 'landscape' is both heterogeneous and dynamic. GPCR function is modulated by bulk membrane properties including membrane fluidity, microdomains, curvature, thickness and asymmetry but GPCRs are also regulated by specific lipid:GPCR binding, including cholesterol and anionic lipids. Understanding the molecular mechanisms whereby GPCR signalling is regulated by lipids is a very active area of research currently. A major advance in membrane protein research in recent years was the application of poly(styrene-co-maleic acid) (SMA) copolymers. These spontaneously generate SMA lipid particles (SMALPs) encapsulating membrane protein in a nano-scale disc of cell membrane, thereby removing the historical need for detergent and preserving lipid:GPCR interaction. The focus of this review is how GPCR-SMALPs are increasing our understanding of GPCR structure and function at the molecular level. Furthermore, an increasing number of 'second generation' SMA-like copolymers have been reported recently. These are reviewed from the context of increasing our understanding of GPCR molecular mechanisms. Moreover, their potential as a novel platform for downstream biophysical and structural analyses is assessed and looking ahead, the translational application of SMA-like copolymers to GPCR drug discovery programmes in the future is considered.

The reproduction process of Gram-positive protocells

Protocells are believed to have existed on early Earth prior to the emergence of prokaryotes. Due to their rudimentary nature, it is widely accepted that these protocells lacked intracellular mechanisms to regulate their reproduction, thereby relying heavily on environmental conditions. To understand protocell reproduction, we adopted a top-down approach of transforming a Gram-positive bacterium into a lipid-vesicle-like state. In this state, cells lacked intrinsic mechanisms to regulate their morphology or reproduction, resembling theoretical propositions on protocells. Subsequently, we grew these proxy-protocells under the environmental conditions of early Earth to understand their impact on protocell reproduction. Despite the lack of molecular biological coordination, cells in our study underwent reproduction in an organized manner. The method and the efficiency of their reproduction can be explained by an interplay between the physicochemical properties of cell constituents and environmental conditions. While the overall reproductive efficiency in these top-down modified cells was lower than their counterparts with a cell wall, the process always resulted in viable daughter cells. Given the simplicity and suitability of this reproduction method to early Earth environmental conditions, we propose that primitive protocells likely reproduced by a process like the one we described below.

Cooperative Membrane Binding of HIV-1 Matrix Proteins

The HIV-1 assembly process begins with a newly synthesized Gag polyprotein being targeted to the inner leaflet of the plasma membrane of the infected cells to form immature viral particles. Gag-membrane interactions are mediated through the myristoylated (Myr) N-terminal matrix (MA) domain of Gag, which eventually multimerize on the membrane to form trimers and higher order oligomers. The study of the structure and dynamics of peripheral membrane proteins like MA has been challenging for both experimental and computational studies due to the complex transient dynamics of protein-membrane interactions. Although the roles of anionic phospholipids (PIP2, PS) and the Myr group in the membrane targeting and stable membrane binding of MA are now well-established, the cooperative interactions between the MA monomers and MA-membrane remain elusive in the context of viral assembly and release. Our present study focuses on the membrane binding dynamics of a higher order oligomeric structure of MA protein (a dimer of trimers), which has not been explored before. Employing time-lagged independent component analysis (tICA) to our microsecond-long trajectories, we investigate conformational changes of the matrix protein induced by membrane binding. Interestingly, the Myr switch of an MA monomer correlates with the conformational switch of adjacent monomers in the same trimer. Together, our findings suggest complex protein dynamics during the formation of the immature HIV-1 lattice; while MA trimerization facilitates Myr insertion, MA trimer-trimer interactions in the immature lattice can hinder the same.

Cholesterol modulates the structural dynamics of the paddle motif loop of KvAP voltage sensor

KvAP is a prokaryotic Kv channel, which has been widely used as a model system to understand voltage- and lipid-dependent gating mechanisms. In phospholipid membranes, the KvAP-VSD adopts the activated/'Up' conformation, whereas the presence of non-phospholipids in membranes favours the structural transition to resting/'Down' state. The S3b-S4 paddle motif loop of KvAP-VSD is functionally important as this participates in protein-protein interactions and is the target for animal toxins. In this study, we have monitored the modulatory role of cholesterol - the physiologically-relevant non-phospholipid - on the organization and dynamics of the S3b-S4 loop of the isolated KvAP-VSD in membranes by site-directed fluorescence approaches using the environmental sensitivity of 7-nitrobenz-2-oxa-1,3-diazol-4-yl-ethylenediamine (NBD) fluorescence. Our results show that cholesterol alters the dynamic nature (rotational and hydration dynamics) of S3b-S4 loop in a segmental fashion, i.e., the residues 110 to 114 and 115 to 117 behave differently in the presence of cholesterol, which is accompanied by considerable change in conformational heterogeneity. Further, quantitative depth measurements using the parallax quenching method reveal that the sensor loop is located at the shallow interfacial region of cholesterol-containing membranes, suggesting that the sensor loop organization is not directly correlated with S4 helix movement. Our results clearly show that cholesterol-induced changes in bilayer properties may not be the predominant factor for the sensor loop's altered structural dynamics, but can be attributed to the conformational change of the KvAP-VSD in cholesterol-containing membranes. Overall, these results are relevant for gating mechanisms, particularly the lipid-dependent gating, of Kv channels in membranes.

Functional regulation of aquaporin dynamics by lipid bilayer composition

With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics are examined. We demonstrate that SoPIP2;1's structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.

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.

Comparative hemolymph proteomic analyses of the freezing and resistance-freezing Ostrinia furnacalis (Guenée)

The Asian corn borer, Ostrinia furnacalis (Guenée) (Lepidoptera: Crambidae), is one of the most harmful pests of maize in Asia. It poses a significant threat to maize production, causing economic losses due to its strong ecological adaptation. In this study, we compared and analyzed the hemolymph proteome between freezing and resistance-freezing O. furnacalis strains using two-dimensional gel electrophoresis to gain insights into the mechanisms of cold resistance. The results revealed that 300-400 hemolymph protein spots were common, with 24 spots showing differences between the two strains. Spectrometry analysis revealed 21 protein spots, including 17 upregulated spots and 4 downregulated ones. The expression of upregulation/downregulation proteins plays a crucial role in the metabolism, energy supply, and defense reaction of insects. Proteomics research not only provides a method for investigating protein expression patterns but also identifies numerous attractive candidates for further exploration.

Membrane-Bound Redox Enzyme Cytochrome bd-I Promotes Carbon Monoxide-Resistant Escherichia coli Growth and Respiration

The terminal oxidases of bacterial aerobic respiratory chains are redox-active electrogenic enzymes that catalyze the four-electron reduction of O2 to 2H2O taking out electrons from quinol or cytochrome c. Living bacteria often deal with carbon monoxide (CO) which can act as both a signaling molecule and a poison. Bacterial terminal oxidases contain hemes; therefore, they are potential targets for CO. However, our knowledge of this issue is limited and contradictory. Here, we investigated the effect of CO on the cell growth and aerobic respiration of three different Escherichia coli mutants, each expressing only one terminal quinol oxidase: cytochrome bd-I, cytochrome bd-II, or cytochrome bo3. We found that following the addition of CO to bd-I-only cells, a minimal effect on growth was observed, whereas the growth of both bd-II-only and bo3-only strains was severely impaired. Consistently, the degree of resistance of aerobic respiration of bd-I-only cells to CO is high, as opposed to high CO sensitivity displayed by bd-II-only and bo3-only cells consuming O2. Such a difference between the oxidases in sensitivity to CO was also observed with isolated membranes of the mutants. Accordingly, O2 consumption of wild-type cells showed relatively low CO sensitivity under conditions favoring the expression of a bd-type oxidase.

Direct Binding Methods to Measure Receptor-Ligand Interactions

G-protein-coupled receptors (GPCRs) contribute to numerous physiological processes via complex network mechanisms. While indirect signaling assays (Ca2+ mobilization, cAMP production, and GTPγS binding) have been useful in identifying and characterizing downstream signaling mechanisms of GPCRs, these methods lack measurements of direct binding affinities, kinetics, binding specificity, and selectivity that are important parameters in GPCR drug discovery. In comparison to existing direct methods that use radio- or fluorescent labels, label-free techniques can closely emulate the native interactions around binding partners. Surface plasmon resonance (SPR) is a label-free technique that utilizes the refractive index (RI) property and is applied widely in quantitative GPCR-ligand binding kinetics measurement including small molecules screening. However, purified GPCRs are further embedded in a synthetic lipid environment which is immobilized through different tags to the SPR sensor surface, resulting in a non-native environment. Here, we introduced a methodology that also uses the RI property to measure binding interactions in a label-free, immobilization-free arrangement. The free-solution technique is successfully applied in quantifying the interaction of bioactive lipids to cognate lipid GPCRs, which is not purified but rather present in near-native conditions, i.e., in milieu of other cytoplasmic lipids and proteins. To further consider the wide applicability of these free-solution approaches in biomolecular interaction research, additional applications on a variety of receptor-ligand pairs are imperative.

Lipid-Centric Design of Plasma Membrane-Mimicking Nanocarriers for Targeted Chemotherapeutic Delivery

The plasma membranes (PM) of mammalian cells contain diverse lipids, proteins, and carbohydrates that are important for systemic recognition and communication in health and disease. Cell membrane coating technology that imparts unique properties of natural plasma membranes to the surface of encapsulated nanoparticles is thus becoming a powerful platform for drug delivery, immunomodulation, and vaccination. However, current coating methods fail to take full advantage of the natural systems because they disrupt the complex and functionally essential features of PMs, most notably the chemical diversity and compositional differences of lipids in two leaflets of the PM. Herein, a new lipid coating approach is reported in which the lipid composition is optimized through a combination of biomimetic and systematic variation approaches for the custom design of nanocarrier systems for precision drug delivery. Nanocarriers coated with the optimized lipids offer unique advantages in terms of bioavailability and efficiency in tumor targeting, tumor penetration, cellular uptake, and drug release. This pilot study provides new insight into the rational design and optimization of nanocarriers for cancer chemotherapeutic drugs and lays the foundation for further customization of cell membrane-mimicking nanocarriers through systematic incorporation of other components.

The bacterial magnesium transporter MgtA reveals highly selective interaction with specific cardiolipin species

The bacterial magnesium transporter A (MgtA) is a specialized P-type ATPase important for Mg2+ import into the cytoplasm; disrupted magnesium homeostasis is linked to intrinsic ribosome instability and antibacterial resistance in Salmonella strains. Here, we show that MgtA has functional specificity for cardiolipin 18:1. Still, it reaches maximum activity only in combination with cardiolipin 16:0, equivalent to the major components of native cardiolipin found in E. coli membranes. Native mass spectrometry indicates the presence of two binding sites for cardiolipin, agreeing with the kinetic studies revealing that a cooperative relationship likely exists between the two cardiolipin variants. This is the first experimental evidence of cooperative effects between lipids of the same class, with only minor variations in their acyl chain composition, acting on a membrane protein. In summary, our results reveal that MgtA exhibits a highly complex interaction with one cardiolipin 18:1 and one cardiolipin 16:0, affecting protein activity and stability, contributing to our understanding of the particular interactions between lipid environment and membrane proteins. Further, a better understanding of Mg2+ homeostasis in bacteria, due to its role as a virulence regulator, will provide further insights into the regulation and mechanism of bacterial infections.

Palmitoylation of the Glucagon-like Peptide-1 Receptor Modulates Cholesterol Interactions at the Receptor-Lipid Microenvironment

The G protein-coupled receptor (GPCR) superfamily of cell surface receptors has been shown to be functionally modulated by post-translational modifications. The glucagon-like peptide receptor-1 (GLP-1R), which is a drug target in diabetes and obesity, undergoes agonist-dependent palmitoyl tail conjugation. The palmitoylation in the C-terminal domain of GLP-1R has been suggested to modulate the receptor-lipid microenvironment. In this work, we have performed coarse-grain molecular dynamics simulations of palmitoylated and nonpalmitoylated GLP-1R to analyze the differential receptor-lipid interactions. Interestingly, the placement and dynamics of the C-terminal domain of GLP-1R are found to be directly dependent on the palmitoyl tail. We observe that both cholesterol and phospholipids interact with the receptor but display differential interactions in the presence and absence of the palmitoyl tail. We characterize important cholesterol-binding sites and validate sites that have been previously reported in experimentally resolved structures of the receptor. We show that the receptor acts like a conduit for cholesterol flip-flop by stabilizing cholesterol in the membrane core. Taken together, our work represents an important step in understanding the molecular effects of lipid modifications in GPCRs.

Cardiolipin binding enhances KcsA channel gating via both its specific and dianion-monoanion interchangeable sites

KcsA is a potassium channel with a plethora of structural and functional information, but its activity in the KcsA-producing actinomycete membranes remains elusive. To determine lipid species involved in channel-modulation, a surface plasmon resonance (SPR)-based methodology, characterized by immobilization of membrane proteins under a membrane environment, was applied. Dianionic cardiolipin (CL) showed extremely higher affinity for KcsA than monoanionic lipids. The SPR experiments further demonstrated that CL bound not only to the N-terminal M0 helix, a lipid-sensor domain, but to the M0 helix-deleted mutant. In contrast, monoanionic lipids interacted primarily with the M0 helix. This indicates the presence of an alternative CL-binding site, plausibly in the transmembrane domain. Single-channel recordings demonstrated that CL enhanced channel opening in an M0-independent manner. Taken together, the action of monoanionic lipids is exclusively mediated by the M0 helix, while CL binds both the M0 helix and its specific site, further enhancing the channel activity.

SIMULTANEOUS NATIVE MASS SPECTROMETRY ANALYSIS OF SINGLE AND DOUBLE MUTANTS TO PROBE LIPID BINDING TO MEMBRANE PROTEINS

Lipids are critical modulators of membrane protein structure and function. However, it is challenging to investigate the thermodynamics of protein-lipid interactions because lipids can simultaneously bind membrane proteins at different sites with different specificities. Here, we developed a native mass spectrometry (MS) approach using single and double mutants to measure the relative energetic contributions of specific residues on Aquaporin Z (AqpZ) toward cardiolipin (CL) binding. We first mutated potential lipid-binding residues on AqpZ, and mixed mutant and wild-type proteins together with CL. By using native MS to simultaneously resolve lipid binding to the mutant and wild-type proteins in a single spectrum, we directly determined the relative affinities of CL binding, thereby revealing the relative Gibbs free energy change for lipid binding caused by the mutation. Comparing different mutants revealed that the W14 contributes to the tightest CL binding site, with R224 contributing to a lower affinity site. Using double mutant cycling, we investigated the synergy between W14 and R224 sites on CL binding. Overall, this novel native MS approach provides unique insights into lipid binding to specific sites on membrane proteins.

Backmapping with Mapping and Isomeric Information

I present a powerful and flexible backmapping tool named Multiscale Simulation Tool (mstool) that converts a coarse-grained (CG) system into all-atom (AA) resolution and only requires AA to CG mapping and isomeric information (cis/trans/dihedral/chiral). The backmapping procedure includes two simple steps: (a) AA atoms are randomly placed near the corresponding CG beads according to the provided mapping scheme. (b) Energy minimization is performed with two modifications in the AA force field (FF). First, nonbonded interactions are replaced with cosine functions to ensure the numerical stability. Second, additional torsions are imposed to maintain the molecules' isomeric properties. To test the simplicity and robustness of the tool, I backmapped multiple membrane and protein CG structures into AA resolution, including a four-bead CG lipid model (resolution increased by a factor of 34) without using intermediate resolution. The tool is freely available at github.com/ksy141/mstool.

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.

Phospholipid Preparations to Characterize Protein-Lipid Interactions In Vitro

The lipid bilayers of the cell are composed of various lipid classes and species. These engage in cell signaling and regulation by recruiting cytosolic proteins to the membrane and interacting with membrane-embedded proteins to alternate their activity and stability. Like lipids, membrane proteins are amphipathic and are stabilized by the hydrophobic forces of the lipid bilayer. Membrane protein-lipid interactions are difficult to investigate since membrane proteins need to be reconstituted in a lipid-mimicking environment. A common and well-established approach is the detergent-based solubilization of the membrane proteins in detergent micelles. Nowadays, nanodiscs and liposomes are used to mimic the lipid bilayer and enable the work with membrane proteins in a more natural environment. However, these protocols need optimization and are labor intensive. The present protocol describes straightforward instructions on how the preparation of lipids is performed and how the lipid detergent mixture is integrated with the membrane protein MARCH5. The lipidation protocol was performed prior to an activity assay specific to membrane-bound E3 ubiquitin ligases and a stability assay that could be used for any membrane protein of choice.

Design and evaluation of tadpole-like conformational antimicrobial peptides

Antimicrobial peptides are promising alternatives to conventional antibiotics. Herein, we report a class of "tadpole-like" peptides consisting of an amphipathic α-helical head and an aromatic tail. A structure-activity relationship (SAR) study of "tadpole-like" temporin-SHf and its analogs revealed that increasing the number of aromatic residues in the tail, introducing Arg to the α-helical head and rearranging the peptide topology dramatically increased antimicrobial activity. Through progressive structural optimization, we obtained two peptides, HT2 and RI-HT2, which exhibited potent antimicrobial activity, no hemolytic activity and cytotoxicity, and no propensity to induce resistance. NMR and molecular dynamics simulations revealed that both peptides indeed adopted "tadpole-like" conformations. Fluorescence experiments and electron microscopy confirmed the membrane targeting mechanisms of the peptides. Our studies not only lead to the discovery of a series of ultrashort peptides with potent broad-spectrum antimicrobial activities, but also provide a new strategy for rational design of novel "tadpole-like" antimicrobial peptides.

Structural modulation of insulin by hydrophobic and hydrophilic molecules

In the bloodstream, insulin interacts with various kinds of molecules, which can alter its structure and modulate its function. In this work, we have synthesized two molecules having extremely hydrophilic and hydrophobic side chains. The effects of hydrophilic and hydrophobic molecules on the binding with insulin have been investigated through a multi-spectroscopic approach. We found that hydrophilic molecules have a slightly higher binding affinity towards insulin. Insulin can bind with the hydrophilic molecules as it binds glucose. The high insulin binding affinity of a hydrophobic molecule indicates its dual nature. The hydrophobic molecule binds at the hydrophobic pocket of the insulin surface, where hydrophilic molecules interact at the polar surface of the insulin. Such binding with the hydrophobic molecule perturbs strongly the secondary structure of the insulin much more in comparison to hydrophilic molecules. Therefore, the stability of insulin decreases in the presence of hydrophobic molecules.

Molecular insights into the effects of focused ultrasound mechanotherapy on lipid bilayers: Unlocking the keys to design effective treatments

Administration of focused ultrasounds (US) represents an attractive complement to classical therapies for a wide range of maladies, from cancer to neurological pathologies, as they are non-invasive, easily targeted, their dosage is easy to control, and they involve low risks. Different mechanisms have been proposed for their activity but the direct effect of their interaction with cell membranes is not well understood at the molecular level. This is in part due to the difficulty of designing experiments able to probe the required spatio-temporal resolutions. Here we use Molecular Dynamics (MD) simulations at two resolution levels and machine learning (ML) classification tools to shed light on the effects that focused US mechanotherapy methods have over a range of lipid bilayers. Our results indicate that the dynamic-structural response of the membrane models to the mechanical perturbations caused by the sound waves strongly depends on the lipid composition. The analyses performed on the MD trajectories contribute to a better understanding of the behavior of lipid membranes, and to open up a path for the rational design of new therapies for the long list of diseases characterized by specific lipid profiles of pathological membrane cells.

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Illuminating the structural basis of human neurokinin 1 receptor (NK1R) antagonism through classical all-atoms molecular dynamics simulations

Advances in structural biology have bestowed insights into the pleiotropic effects of neurokinin 1 receptors (NK1R) in diverse patho-physiological processes, thereby highlighting the potential therapeutic value of antagonists directed against NK1R. Herein, we investigate the mode of antagonist recognition to discern the obscure atomic facets germane for the function and molecular determinants of NK1R. To commence discernment of potent antagonists and the conformational changes in NK1R, induced upon antagonist binding, state-of-the-art classical all-atoms molecular dynamics (MD) simulations in lipid mimetic bilayers have been utilized. MD simulations of structural ensembles reveals the involvement of TM5 and TM6 in tight anchoring of antagonists through a network of interhelical hydrogen-bonds, while, the extracellular loop 2 (ECL2) governs the overall size and nature of the pocket, thereby modulating NK1R. Consistent comparison between experiments and MD simulation results discerns the predominant role of TM3, TM4, and TM6 in lipid-NK1R interaction. Correlation between hydrophobic index and helicity of TM domains elucidates their importance in maintaining the structural stability in addition to regulating NK1R antagonism. Taken together, we anticipate that our computational study marks a comprehensive structural basis of NK1R antagonism in lipid bilayers, which may facilitate designing of new therapeutics against associated diseases targeting human neurokinin receptors.

Martini 3 Coarse-Grained Force Field for Cholesterol

Cholesterol plays a crucial role in biomembranes by regulating various properties, such as fluidity, rigidity, permeability, and organization of lipid bilayers. The latest version of the Martini model, Martini 3, offers significant improvements in interaction balance, molecular packing, and inclusion of new bead types and sizes. However, the release of the new model resulted in the need to reparameterize many core molecules, including cholesterol. Here, we describe the development and validation of a Martini 3 cholesterol model, addressing issues related to its bonded setup, shape, volume, and hydrophobicity. The proposed model mitigates some limitations of its Martini 2 predecessor while maintaining or improving the overall behavior.

The terminal oxidase cytochrome bd-I confers carbon monoxide resistance to Escherichia coli cells

Carbon monoxide (CO) plays a multifaceted role in the physiology of organisms, from poison to signaling molecule. Heme proteins, including terminal oxidases, are plausible CO targets. Three quinol oxidases terminate the branched aerobic respiratory chain of Escherichia coli. These are the heme‑copper cytochrome bo3 and two copper-lacking bd-type cytochromes, bd-I and bd-II. All three enzymes generate a proton motive force during the four-electron oxygen reduction reaction that is used for ATP production. The bd-type oxidases also contribute to mechanisms of bacterial defense against various types of stresses. Here we report that in E. coli cells, at the enzyme concentrations tested, cytochrome bd-I is much more resistant to inhibition by CO than cytochrome bd-II and cytochrome bo3. The apparent half-maximal inhibitory concentration values, IC50, for inhibition of O2 consumption of the membrane-bound bd-II and bo3 oxidases by CO at ~150 μM O2 were estimated to be 187.1 ± 11.1 and 183.3 ± 13.5 μM CO, respectively. Under the same conditions, the maximum inhibition observed with the membrane-bound cytochrome bd-I was 20 ± 10% at ~200 μM CO.

The Important Role of Transporter Structures in Drug Disposition, Efficacy, and Toxicity

The ATP-binding cassette (ABC) and solute carrier (SLC) transporters are critical determinants of drug disposition, clinical efficacy, and toxicity as they specifically mediate the influx and efflux of various substrates and drugs. ABC transporters can modulate the pharmacokinetics of many drugs via mediating the translocation of drugs across biologic membranes. SLC transporters are important drug targets involved in the uptake of a broad range of compounds across the membrane. However, high-resolution experimental structures have been reported for a very limited number of transporters, which limits the study of their physiologic functions. In this review, we collected structural information on ABC and SLC transporters and described the application of computational methods in structure prediction. Taking P-glycoprotein (ABCB1) and serotonin transporter (SLC6A4) as examples, we assessed the pivotal role of structure in transport mechanisms, details of ligand-receptor interactions, drug selectivity, the molecular mechanisms of drug-drug interactions, and differences caused by genetic polymorphisms. The data collected contributes toward safer and more effective pharmacological treatments. SIGNIFICANCE STATEMENT: The experimental structure of ATP-binding cassette and solute carrier transporters was collected, and the application of computational methods in structure prediction was described. P-glycoprotein and serotonin transporter were used as examples to reveal the pivotal role of structure in transport mechanisms, drug selectivity, the molecular mechanisms of drug-drug interactions, and differences caused by genetic polymorphisms.