The cellular membrane as a mediator for small molecule interaction with membrane proteins

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.

Conformational State of Fenamates at the Membrane Interface: A MAS NOESY Study

The present work analyzes the ¹H NOESY MAS NMR spectra of three fenamates (mefenamic, tolfenamic, and flufenamic acids) localized in the lipid-water interface of phosphatidyloleoylphosphatidylcholine (POPC) membranes. The observed cross-peaks in the two-dimensional NMR spectra characterized intramolecular proximities between the hydrogen atoms of the fenamates as well as intermolecular interactions between the fenamates and POPC molecules. The peak amplitude normalization for an improved cross-relaxation (PANIC) approach, the isolated spin-pair approximation (ISPA) model, and the two-position exchange model were used to calculate the interproton distances indicative of specific conformations of the fenamates. The results showed that the proportions of the A+C and B+D conformer groups of mefenamic and tolfenamic acids in the presence of POPC were comparable within the experimental error and amounted to 47.8%/52.2% and 47.7%/52.3%, respectively. In contrast, these proportions for the flufenamic acid conformers differed and amounted to 56.6%/43.4%. This allowed us to conclude that when they bind to the POPC model lipid membrane, fenamate molecules change their conformational equilibria.

Differential dynamics and direct interaction of bound ligands with lipids in multidrug transporter ABCG2

ABCG2 is an ATP-binding cassette (ABC) transporter that extrudes a wide range of xenobiotics and drugs from the cell and contributes to multidrug resistance in cancer cells. Following our recent structural characterization of topotecan-bound ABCG2, here, we present cryo-EM structures of ABCG2 under turnover conditions in complex with a special modulator and slow substrate, tariquidar, in nanodiscs. The structures reveal that similar to topotecan, tariquidar induces two distinct ABCG2 conformations under turnover conditions (turnover-1 and turnover-2). μs-scale molecular dynamics simulations of drug-bound and apo ABCG2 in native-like lipid bilayers, in both topotecan- and tariquidar-bound states, characterize the ligand size as a major determinant of its binding stability. The simulations highlight direct lipid-drug interactions for the smaller topotecan, which exhibits a highly dynamic binding mode. In contrast, the larger tariquidar occupies most of the available volume in the binding pocket, thus leaving little space for lipids to enter the cavity and interact with it. Similarly, when simulating ABCG2 in the apo inward-open state, we also observe spontaneous penetration of phospholipids into the binding cavity. The captured phospholipid diffusion pathway into ABCG2 offers a putative general path to recruit any hydrophobic/amphiphilic substrates directly from the membrane. Our simulations also reveal that ABCG2 rejects cholesterol as a substrate, which is omnipresent in plasma membranes that contain ABCG2. At the same time, cholesterol is found to prohibit the penetration of phospholipids into ABCG2. These molecular findings have direct functional ramifications on ABCG2's function as a transporter.

The interaction of steroids with phospholipid bilayers and membranes

Steroids are critical for various physiological processes and used to treat inflammatory conditions. Steroids act by two distinct pathways. The genomic pathway is initiated by the steroid binding to nuclear receptors while the non-genomic pathway involves plasma membrane receptors. It has been proposed that steroids might also act in a more indirect mechanism by altering biophysical properties of membranes. Yet, little is known about the effect of steroids on membranes, and steroid-membrane interactions are complex and challenging to characterise. The focus of this review is to outline what is currently known about the interactions of steroids with phospholipid bilayers and illustrate the complexity of these systems using cortisone and progesterone as the main examples. The combined findings from current work demonstrate that the hydrophobicity and planarity of the steroid core does not provide a consensus for steroid-membrane interactions. Even small differences in the substituents on the steroid core can result in significant changes in steroid-membrane interactions. Furthermore, steroid-induced changes in phospholipid bilayer properties are often dependent on steroid concentration and lipid composition. This complexity means that currently there is insufficient information to establish a reliable structure-activity relationship to describe the effect of steroids on membrane properties. Future work should address the challenge of connecting the findings from studying the effect of steroids on phospholipid bilayers to cell membranes. Insights from steroid-membrane interactions will benefit our understanding of normal physiology and assist drug development.

Sustained-release ketamine-loaded lipid-particulate system: in vivo assessment in mice

Ketamine is used as an analgesic adjuvant in patients with chronic cancer-related pain. However, ketamine's short half-life requires frequent dose administration. Our aim was to develop a sustained release formulation of ketamine with high loading and to evaluate the in vivo pharmacokinetics and biodistribution in mice. Here, ketamine hydrochloride sustained-release lipid particles (KSL) were developed using the thin-film hydration method. The mean (± SD) encapsulation efficiency (EE) and drug loading (DL) of KSL were 65.6 (± 1.7)% and 72.4 (± 0.5)% respectively, and the mean (± SD) size of the lipid particles and the polydispersity index were 738 (± 137) nm and 0.44 (± 0.02) respectively. The release period of KSL in pH 7.4 medium was 100% complete within 8 h in vitro but a sustained-release profile was observed for more than 5 days after intravenous injection in mice. Importantly, the KSL formulation resulted in a 27-fold increase in terminal half-life, a threefold increase in systemic exposure (AUC_0-∞), and a threefold decrease in clearance compared with the corresponding pharmacokinetics for intravenous ketamine itself. Our findings demonstrate high encapsulation efficiency of ketamine in the sustained-release KSL formulation with prolonged release in mice after systemic dose administration despite 100% in vitro release within 8 h that requires future investigation.

Tepotinib Inhibits Several Drug Efflux Transporters and Biotransformation Enzymes: The Role in Drug-Drug Interactions and Targeting Cytostatic Resistance In Vitro and Ex Vivo

Tepotinib is a novel tyrosine kinase inhibitor recently approved for the treatment of non-small cell lung cancer (NSCLC). In this study, we evaluated the tepotinib's potential to perpetrate pharmacokinetic drug interactions and modulate multidrug resistance (MDR). Accumulation studies showed that tepotinib potently inhibits ABCB1 and ABCG2 efflux transporters, which was confirmed by molecular docking. In addition, tepotinib inhibited several recombinant cytochrome P450 (CYP) isoforms with varying potency. In subsequent drug combination experiments, tepotinib synergistically reversed daunorubicin and mitoxantrone resistance in cells with ABCB1 and ABCG2 overexpression, respectively. Remarkably, MDR-modulatory properties were confirmed in ex vivo explants derived from NSCLC patients. Furthermore, we demonstrated that anticancer effect of tepotinib is not influenced by the presence of ABC transporters associated with MDR, although monolayer transport assays designated it as ABCB1 substrate. Finally, tested drug was observed to have negligible effect on the expression of clinically relevant drug efflux transporters and CYP enzymes. In conclusion, our findings provide complex overview on the tepotinib's drug interaction profile and suggest a promising novel therapeutic strategy for future clinical investigations.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

The Whole Is Bigger than the Sum of Its Parts: Drug Transport in the Context of Two Membranes with Active Efflux

Cell envelope plays a dual role in the life of bacteria by simultaneously protecting it from a hostile environment and facilitating access to beneficial molecules. At the heart of this ability lie the restrictive properties of the cellular membrane augmented by efflux transporters, which preclude intracellular penetration of most molecules except with the help of specialized uptake mediators. Recently, kinetic properties of the cell envelope came into focus driven on one hand by the urgent need in new antibiotics and, on the other hand, by experimental and theoretical advances in studies of transmembrane transport. A notable result from these studies is the development of a kinetic formalism that integrates the Michaelis-Menten behavior of individual transporters with transmembrane diffusion and offers a quantitative basis for the analysis of intracellular penetration of bioactive compounds. This review surveys key experimental and computational approaches to the investigation of transport by individual translocators and in whole cells, summarizes key findings from these studies and outlines implications for antibiotic discovery. Special emphasis is placed on Gram-negative bacteria, whose envelope contains two separate membranes. This feature sets these organisms apart from Gram-positive bacteria and eukaryotic cells by providing them with full benefits of the synergy between slow transmembrane diffusion and active efflux.

Reduction in hypoxia-reoxygenation-induced myocardial mitochondrial damage with exogenous methane

Albeit previous experiments suggest potential anti‐inflammatory effect of exogenous methane (CH₄) in various organs, the mechanism of its bioactivity is not entirely understood. We aimed to investigate the potential mitochondrial effects and the underlying mechanisms of CH₄ in rat cardiomyocytes and mitochondria under simulated ischaemia/reperfusion (sI/R) conditions. Three‐day‐old cultured cardiomyocytes were treated with 2.2% CH₄‐artificial air mixture during 2‐hour‐long reoxygenation following 4‐hour‐long anoxia (sI/R and sI/R + CH₄, n = 6‐6), with normoxic groups serving as controls (SH and SH + CH₄; n = 6‐6). Mitochondrial functions were investigated with high‐resolution respirometry, and mitochondrial membrane injury was detected by cytochrome c release and apoptotic characteristics by using TUNEL staining. CH₄ admixture had no effect on complex II (CII)‐linked respiration under normoxia but significantly decreased the complex I (CI)‐linked oxygen consumption. Nevertheless, addition of CH₄ in the sI/R + CH4 group significantly reduced the respiratory activity of CII in contrast to CI and the CH₄ treatment diminished mitochondrial H₂O₂ production. Substrate‐induced changes to membrane potential were partially preserved by CH₄, and additionally, cytochrome c release and apoptosis of cardiomyocytes were reduced in the CH₄‐treated group. In conclusion, the addition of CH₄ decreases mitochondrial ROS generation via blockade of electron transport at CI and reduces anoxia‐reoxygenation‐induced mitochondrial dysfunction and cardiomyocyte injury in vitro.

The ecological roles of microbial lipopeptides: Where are we going?

Lipopeptides (LPs) are secondary metabolites produced by a diversity of bacteria and fungi. Their unique chemical structure comprises both a peptide and a lipid moiety. LPs are of major biotechnological interest owing to their emulsification, antitumor, immunomodulatory, and antimicrobial activities. To date, these versatile compounds have been applied across multiple industries, from pharmaceuticals through to food processing, cosmetics, agriculture, heavy metal, and hydrocarbon bioremediation. The variety of LP structures and the diversity of the environments from which LP-producing microorganisms have been isolated suggest important functions in their natural environment. However, our understanding of the ecological role of LPs is limited. In this review, the mode of action and the role of LPs in motility, antimicrobial activity, heavy metals removal and biofilm formation are addressed. We include discussion on the need to characterise LPs from a diversity of microorganisms, with a focus on taxa inhabiting 'extreme' environments. We introduce the use of computational target fishing and molecular dynamics simulations as powerful tools to investigate the process of interaction between LPs and cell membranes. Together, these advances will provide new understanding of the mechanism of action of novel LPs, providing greater insights into the roles of LPs in the natural environment.

On the Rational Drug Design for Hypertension through NMR Spectroscopy

Antagonists of the AT1receptor (AT1R) are beneficial molecules that can prevent the peptide hormone angiotensin II from binding and activating the specific receptor causing hypertension in pathological states. This review article summarizes the multifaced applications of solid and liquid state high resolution nuclear magnetic resonance (NMR) spectroscopy in antihypertensive commercial drugs that act as AT1R antagonists. The 3D architecture of these compounds is explored through 2D NOESY spectroscopy and their interactions with micelles and lipid bilayers are described using solid state 13CP/MAS, 31P and 2H static solid state NMR spectroscopy. Due to their hydrophobic character, AT1R antagonists do not exert their optimum profile on the AT1R. Therefore, various vehicles are explored so as to effectively deliver these molecules to the site of action and to enhance their pharmaceutical efficacy. Cyclodextrins and polymers comprise successful examples of effective drug delivery vehicles, widely used for the delivery of hydrophobic drugs to the active site of the receptor. High resolution NMR spectroscopy provides valuable information on the physical-chemical forces that govern these drug:vehicle interactions, knowledge required to get a deeper understanding on the stability of the formed complexes and therefore the appropriateness and usefulness of the drug delivery system. In addition, it provides valuable information on the rational design towards the synthesis of more stable and efficient drug formulations.

Mechanistic Understanding From Molecular Dynamics Simulation in Pharmaceutical Research 1: Drug Delivery

In this review, we outline the growing role that molecular dynamics simulation is able to play as a design tool in drug delivery. We cover both the pharmaceutical and computational backgrounds, in a pedagogical fashion, as this review is designed to be equally accessible to pharmaceutical researchers interested in what this new computational tool is capable of and experts in molecular modeling who wish to pursue pharmaceutical applications as a context for their research. The field has become too broad for us to concisely describe all work that has been carried out; many comprehensive reviews on subtopics of this area are cited. We discuss the insight molecular dynamics modeling has provided in dissolution and solubility, however, the majority of the discussion is focused on nanomedicine: the development of nanoscale drug delivery vehicles. Here we focus on three areas where molecular dynamics modeling has had a particularly strong impact: (1) behavior in the bloodstream and protective polymer corona, (2) Drug loading and controlled release, and (3) Nanoparticle interaction with both model and biological membranes. We conclude with some thoughts on the role that molecular dynamics simulation can grow to play in the development of new drug delivery systems.

In silico approach: biological prediction of nordentatin derivatives as anticancer agent inhibitors in the cAMP pathway

A combination of computational techniques has been carried out to predict the binding of nordentatin derivatives based on pyranocoumarin semi-synthesis with the target protein from the expression of the PDE4B gene. The inhibition of the cAMP pathway is the main target of anti-cancer drugs, which is responsible for uncontrolled cell division in cancer. Modeling was done using a combination of semi-empirical methods and the density functional theory (PM3-DFT/6-31G*/B3LYP) to obtain the optimal structure of a small ligand that could be modeled. Studies on the interaction of the ligands and amino acid residues on protein targets were carried out using a combination of molecular docking and molecular dynamic simulation. Molecular docking based on functional grid scores showed a very good native ligand pose with an RMSD of 0.93 Å in determining the initial coordinates of the ligand-receptor interactions. Furthermore, the amino acid residues responsible for interaction through H-bonds were Tyr103, His104, His177, Met217, and Gln313. The binding free energy (kcal mol-1) results of the candidates were PS-1 (-36.84 ± 0.31), PS-2 (-35.34 ± 0.28), PS-3 (-26.65 ± 0.30), PS-5 (-42.66 ± 0.26), PS-7 (-35.33 ± 0.23), and PS-9 (-32.57 ± 0.20), which are smaller than that of the native ligand Z72 (-24.20 ± 0.19), and thus these have good potential as drugs that can inhibit the cAMP pathway. These results provide theoretical information for the efficient inhibition of the cAMP pathway in the future.

Interaction of Cytochrome C Oxidase with Steroid Hormones

Estradiol, testosterone and other steroid hormones inhibit cytochrome c oxidase (CcO) purified from bovine heart. The inhibition is strongly dependent on concentration of dodecyl-maltoside (DM) in the assay. The plots of K_i vs [DM] are linear for both estradiol and testosterone which may indicate an 1:1 stoichiometry competition between the hormones and the detergent. Binding of estradiol, but not of testosterone, brings about spectral shift of the oxidized CcO consistent with an effect on heme a₃³⁺. We presume that the hormones bind to CcO at the bile acid binding site described by Ferguson-Miller and collaborators. Estradiol is shown to inhibit intraprotein electron transfer between hemes a and a₃. Notably, neither estradiol nor testosterone suppresses the peroxidase activity of CcO. Such a specific mode of action indicates that inhibition of CcO activity by the hormones is associated with impairing proton transfer via the K-proton channel.

PLD2-PI(4,5)P2 interactions in fluid phase membranes: Structural modeling and molecular dynamics simulations

Interaction of phospholipase D2 (PLD2) with phosphatidylinositol (4,5)-bisphosphate (PIP2) is regarded as the critical step of numerous physiological processes. Here we build a full-length model of human PLD2 (hPLD2) combining template-based and ab initio modeling techniques and use microsecond all-atom molecular dynamics (MD) simulations of the protein in contact with a complex membrane to determine hPLD2-PIP2 interactions. MD simulations reveal that the intermolecular interactions preferentially occur between specific PIP2 phosphate groups and hPLD2 residues; the most strongly interacting residues are arginine at the pbox consensus sequence (PX) and pleckstrin homology (PH) domain. Interaction networks indicate formation of clusters at the protein-membrane interface consisting of amino acids, PIP2, and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidic acid (POPA); the largest cluster was in the PH domain.

Molecular Dynamics Simulation of Small Molecules Interacting with Biological Membranes

Cell membranes protect and compartmentalise cells and their organelles. The semi-permeable nature of these membranes controls the exchange of solutes across their structure. Characterising the interaction of small molecules with biological membranes is critical to understanding of physiological processes, drug action and permeation, and many biotechnological applications. This review provides an overview of how molecular simulations are used to study the interaction of small molecules with biological membranes, with a particular focus on the interactions of water, organic compounds, drugs and short peptides with models of plasma cell membrane and stratum corneum lipid bilayers. This review will not delve on other types of membranes which might have different composition and arrangement, such as thylakoid or mitochondrial membranes. The application of unbiased molecular dynamics simulations and enhanced sampling methods such as umbrella sampling, metadynamics and replica exchange are described using key examples. This review demonstrates how state-of-the-art molecular simulations have been used successfully to describe the mechanism of binding and permeation of small molecules with biological membranes, as well as associated changes to the structure and dynamics of these membranes. The review concludes with an outlook on future directions in this field.

Lipidomic analysis of human primary hepatocytes following LXR activation with GW3965 identifies AGXT2L1 as a main target associated to changes in phosphatidylethanolamine

Liver X receptor (LXR) agonists have the potential to alleviate obesity related diseases, particularly atherosclerosis. However, LXRs are transcriptional regulators that induce de novo lipogenesis and lipid accumulation in hepatocytes which represents a serious adverse effect. In this work, we sought to characterize the LXR agonist GW3965 effects on fatty acid (FA) and phospholipid (PL) remodelling and the correlation with gene expression in order to better understand the underlying effects leading to hepatic pathology upon LXR activation. Human primary hepatocytes treated for 48 h with GW3965 were analysed for changes in lipid metabolism gene expression by qPCR, variations in the FA profile was evaluated by GC-FID and in PL profiles using thin layer chromatography, ESI-MS and MS/MS analysis. Changes in cell membrane biochemical properties were studied using bilayer models generated with CHARMM-GUI. ELOLV6 and SCD1 mRNA increase was consistent with higher C16:1 and C18:1n9 at the expense of C16:0 and C18:0. The reduction of C18:2n6 and increase in C20:2n6 was in agreement with ELOVL5 upregulation. Phosphatydilethanolamine (PE) levels tended to decrease and phosphatidylinositol to increase; although differences did not reach significance, they correlated with changes in AGXT2L1, CDS1 and LPIN1 mRNA levels that were increased. The overall effect of GW3965 on PEs molecular profiles was an increase of long-chain polyunsaturated FA chains and a decrease of C16/C18 saturated and monounsaturated FAs chains. Additionally, PC (32:1) and PC (34:2) were decreased, and PC (36:1) and PC (34:1) were increased. AGXT2L1 is an enzyme with strict substrate specificity for phosphoethanolamine, which is converted into ammonia in GW3965-treated hepatocytes and could explain the PE reduction. In summary, LXR activation by GW3965 targets PE biosynthesis and FA elongation/desaturation, which tends to decrease PE in relation to total PL levels, and remodelling of PC and PE molecular species. We identified the human AGXT2L1 gene as induced by LXR activation by both synthetic and endogenous agonist treatment. The increase in acetaldehyde-induced oxidative stress, and in the lipid species identified have the potential to enhance the inflammatory process and impair membrane function. Future studies should focus on inhibition of AGXT2L1 activity with the aim of reverting the steatosis induced by LXR activation.

Ensartinib (X-396) Effectively Modulates Pharmacokinetic Resistance Mediated by ABCB1 and ABCG2 Drug Efflux Transporters and CYP3A4 Biotransformation Enzyme

Ensartinib (X-396) is a promising tyrosine kinase inhibitor currently undergoing advanced clinical evaluation for the treatment of non-small cell lung cancer. In this work, we investigate possible interactions of this promising drug candidate with ATP-binding cassette (ABC) drug efflux transporters and cytochrome P450 biotransformation enzymes (CYPs), which play major roles in multidrug resistance (MDR) and pharmacokinetic drug-drug interactions (DDIs). Accumulation studies showed that ensartinib is a potent inhibitor of ABCB1 and ABCG2 transporters. Additionally, incubation experiments with recombinant CYPs showed that ensartinib significantly inhibits CYP3A4 and CYP2C9. Subsequent molecular docking studies confirmed these findings. Drug combination experiments demonstrated that ensartinib synergistically potentiates the antiproliferative effects of daunorubicin, mitoxantrone, and docetaxel in ABCB1, ABCG2, and CYP3A4-overexpressing cellular models, respectively. Advantageously, ensartinib’s antitumor efficiency was not compromised by the presence of MDR-associated ABC transporters, although it acted as a substrate of ABCB1 in Madin-Darby Canine Kidney II (MDCKII) monolayer transport assays. Finally, we demonstrated that ensartinib had no significant effect on the mRNA-level expression of examined transporters and enzymes in physiological and lung tumor cellular models. In conclusion, ensartinib may perpetrate clinically relevant pharmacokinetic DDIs and modulate ABCB1-, ABCG2-, and CYP3A4-mediated MDR. The in vitro findings presented here will provide a valuable foundation for future in vivo investigations.

Menthol Binding to the Human α4β2 Nicotinic Acetylcholine Receptor Facilitated by Its Strong Partitioning in the Membrane

We utilize various computational methodologies to study menthol's interaction with multiple organic phases, a lipid bilayer, and the human α4β2 nicotinic acetylcholine receptor (nAChR), the most abundant nAChR in the brain. First, force field parameters developed for menthol are validated in alchemical free energy perturbation simulations to calculate solvation free energies of menthol in water, dodecane, and octanol and compare the results against experimental data. Next, umbrella sampling is used to construct the free energy profile of menthol permeation across a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. The results from a flooding simulation designed to study the water-membrane partitioning of menthol in a POPC lipid bilayer are used to determine the penetration depth and the preferred orientation of menthol in the bilayer. Finally, employing both docking and flooding simulations, menthol is shown to bind to different sites on the human α4β2 nAChR. The most likely binding mode of menthol to a desensitized membrane-embedded α4β2 nAChR is identified to be via a membrane-mediated pathway in which menthol binds to the sites at the lipid-protein interface after partitioning in the membrane. A rare but distinct binding mode in which menthol binds to the extracellular opening of receptor's ion permeation pore is also reported.

Computational Dissection of Membrane Transport at a Microscopic Level

Membrane transporters are key gatekeeper proteins at cellular membranes that closely control the traffic of materials. Their function relies on structural rearrangements of varying degrees that facilitate substrate translocation across the membrane. Characterizing these functionally important molecular events at a microscopic level is key to our understanding of membrane transport, yet challenging to achieve experimentally. Recent advances in simulation technology and computing power have rendered molecular dynamics (MD) simulation a powerful biophysical tool to investigate a wide range of dynamical events spanning multiple spatial and temporal scales. Here, we review recent studies of diverse membrane transporters using computational methods, with an emphasis on highlighting the technical challenges, key lessons learned, and new opportunities to illuminate transporter structure and function.

Sequential activation of STIM1 links Ca2+ with luminal domain unfolding

The stromal interaction molecule 1 (STIM1) has two important functions, Ca2+ sensing within the endoplasmic reticulum and activation of the store-operated Ca2+ channel Orai1, enabling plasma-membrane Ca2+ influx. We combined molecular dynamics (MD) simulations with live-cell recordings and determined the sequential Ca2+-dependent conformations of the luminal STIM1 domain upon activation. Furthermore, we identified the residues within the canonical and noncanonical EF-hand domains that can bind to multiple Ca2+ ions. In MD simulations, a single Ca2+ ion was sufficient to stabilize the luminal STIM1 complex. Ca2+ store depletion destabilized the two EF hands, triggering disassembly of the hydrophobic cleft that they form together with the stable SAM domain. Point mutations associated with tubular aggregate myopathy or cancer that targeted the canonical EF hand, and the hydrophobic cleft yielded constitutively clustered STIM1, which was associated with activation of Ca2+ entry through Orai1 channels. On the basis of our results, we present a model of STIM1 Ca2+ binding and refine the currently known initial steps of STIM1 activation on a molecular level.

A Perspective: Active Role of Lipids in Neurotransmitter Dynamics

Synaptic neurotransmission is generally considered as a function of membrane-embedded receptors and ion channels in response to the neurotransmitter (NT) release and binding. This perspective aims to widen the protein-centric view by including another vital component—the synaptic membrane—in the discussion. A vast set of atomistic molecular dynamics simulations and biophysical experiments indicate that NTs are divided into membrane-binding and membrane-nonbinding categories. The binary choice takes place at the water-membrane interface and follows closely the positioning of the receptors’ binding sites in relation to the membrane. Accordingly, when a lipophilic NT is on route to a membrane-buried binding site, it adheres on the membrane and, then, travels along its plane towards the receptor. In contrast, lipophobic NTs, which are destined to bind into receptors with extracellular binding sites, prefer the water phase. This membrane-based sorting splits the neurotransmission into membrane-independent and membrane-dependent mechanisms and should make the NT binding into the receptors more efficient than random diffusion would allow. The potential implications and notable exceptions to the mechanisms are discussed here. Importantly, maintaining specific membrane lipid compositions (MLCs) at the synapses, especially regarding anionic lipids, affect the level of NT-membrane association. These effects provide a plausible link between the MLC imbalances and neurological diseases such as depression or Parkinson’s disease. Moreover, the membrane plays a vital role in other phases of the NT life cycle, including storage and release from the synaptic vesicles, transport from the synaptic cleft, as well as their synthesis and degradation.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Flexible lipid nanomaterials studied by NMR spectroscopy

Advances in solid-state nuclear magnetic resonance spectroscopy inform the emergence of material properties from atomistic-level interactions in membrane lipid nanostructures.

Electrostatic lock in the transport cycle of the multidrug resistance transporter EmrE

EmrE is a small, homodimeric membrane transporter that exploits the established electrochemical proton gradient across the Escherichia coli inner membrane to export toxic polyaromatic cations, prototypical of the wider small-multidrug resistance transporter family. While prior studies have established many fundamental aspects of the specificity and rate of substrate transport in EmrE, low resolution of available structures has hampered identification of the transport coupling mechanism. Here we present a complete, refined atomic structure of EmrE optimized against available cryo-electron microscopy (cryo-EM) data to delineate the critical interactions by which EmrE regulates its conformation during the transport process. With the model, we conduct molecular dynamics simulations of the transporter in explicit membranes to probe EmrE dynamics under different substrate loading and conformational states, representing different intermediates in the transport cycle. The refined model is stable under extended simulation. The water dynamics in simulation indicate that the hydrogen-bonding networks around a pair of solvent-exposed glutamate residues (E14) depend on the loading state of EmrE. One specific hydrogen bond from a tyrosine (Y60) on one monomer to a glutamate (E14) on the opposite monomer is especially critical, as it locks the protein conformation when the glutamate is deprotonated. The hydrogen bond provided by Y60 lowers the [Formula: see text] of one glutamate relative to the other, suggesting both glutamates should be protonated for the hydrogen bond to break and a substrate-free transition to take place. These findings establish the molecular mechanism for the coupling between proton transfer reactions and protein conformation in this proton-coupled secondary transporter.

Micro-pharmacokinetics: Quantifying local drug concentration at live cell membranes

Fundamental equations for determining pharmacological parameters, such as the binding affinity of a ligand for its target receptor, assume a homogeneous distribution of ligand, with concentrations in the immediate vicinity of the receptor being the same as those in the bulk aqueous phase. It is, however, known that drugs are able to interact directly with the plasma membrane, potentially increasing local ligand concentrations around the receptor. We have previously reported an influence of ligand-phospholipid interactions on ligand binding kinetics at the β₂-adrenoceptor, which resulted in distinct "micro-pharmacokinetic" ligand profiles. Here, we directly quantified the local concentration of BODIPY630/650-PEG8-S-propranolol (BY-propranolol), a fluorescent derivative of the classical β-blocker propranolol, at various distances above membranes of single living cells using fluorescence correlation spectroscopy. We show for the first time a significantly increased ligand concentration immediately adjacent to the cell membrane compared to the bulk aqueous phase. We further show a clear role of both the cell membrane and the β₂-adrenoceptor in determining high local BY-propranolol concentrations at the cell surface. These data suggest that the true binding affinity of BY-propranolol for the β₂-adrenoceptor is likely far lower than previously reported and highlights the critical importance of understanding the "micro-pharmacokinetic" profiles of ligands for membrane-associated proteins.

Nanodiscs in Membrane Biochemistry and Biophysics

Membrane proteins play a most important part in metabolism, signaling, cell motility, transport, development, and many other biochemical and biophysical processes which constitute fundamentals of life on the molecular level. Detailed understanding of these processes is necessary for the progress of life sciences and biomedical applications. Nanodiscs provide a new and powerful tool for a broad spectrum of biochemical and biophysical studies of membrane proteins and are commonly acknowledged as an optimal membrane mimetic system that provides control over size, composition, and specific functional modifications on the nanometer scale. In this review we attempted to combine a comprehensive list of various applications of nanodisc technology with systematic analysis of the most attractive features of this system and advantages provided by nanodiscs for structural and mechanistic studies of membrane proteins.

On Atomistic Models for Molecular Oxygen

Molecular oxygen (O₂) is key to all life on earth, as it is constantly cycled via photosynthesis and cellular respiration. Substantial scientific effort has been devoted to understanding every part of this cycle. Classical molecular dynamics (MD) simulations have been used to study some of the key processes involved in cellular respiration: O₂ permeation through alveolar monolayers and cellular membranes, its binding to hemoglobin during transport in the bloodstream, as well as its transport along optimal pathways toward its reduction sites in proteins. Moreover, MD simulations can help interpret the results of several imaging techniques in which O₂ is used because of its paramagnetic nature. However, despite the widespread use of computational models for the O₂ molecule, their performances have never been systematically evaluated. In this paper, we assess the performances of 14 different models of O₂ available in the literature by calculating four thermodynamic properties: density, heat of vaporization, free energy of hydration, and free energy of solvation in hexadecane. For each property, reliable experimental data are available. Most models perform reasonably well in predicting the correct trends, but they fail to reproduce the experimental data quantitatively. We then develop new models for O₂, with and without a quadrupole moment, and compare their behavior with the behavior of previously published models. The new models show significant improvement in terms of density, heat of vaporization, and free energy of hydration. However, quantitative agreement with water-oil partitioning is not reached due to discrepancies between the calculated and measured free energies of solvation in hexadecane. We suggest that classical pairwise-additive models may be inadequate to properly describe the thermodynamics of solvation of apolar species, such as O₂, in apolar solvents.

Molecular Dynamics of Photosystem II Embedded in the Thylakoid Membrane

Photosystem II (PSII) is one of the key protein complexes in photosynthesis. We introduce a coarse grained model of PSII and present the analysis of 60 μs molecular dynamics simulations of PSII in both monomeric and dimeric form, embedded in a thylakoid membrane model that reflects its native lipid composition. We describe in detail the setup of the protein complex and the many natural cofactors and characterize their mobility. Overall we find that the protein subunits and cofactors are more flexible toward the periphery of the complex as well as near the PLQ exchange cavity and at the dimer interface. Of all cofactors, β-carotenes show the highest mobility. Some of the β-carotenes diffuse in and out of the protein complex via the thylakoid membrane. In contrast with the PSII dimer, the monomeric form adopts a tilted conformation in the membrane, with strong interactions between the soluble PsbO subunit and the glycolipid headgroups. Interestingly, the tilted conformation causes buckling of the membrane. Together, our results provide an unprecedented view of PSII dynamics on a microsecond time scale. Our data may be used as basis for the interpretation of experimental data as well as for theoretical models describing exciton energy transfer.