Psychedelics promote plasticity by directly binding to BDNF receptor TrkB

Psychedelics produce fast and persistent antidepressant effects and induce neuroplasticity resembling the effects of clinically approved antidepressants. We recently reported that pharmacologically diverse antidepressants, including fluoxetine and ketamine, act by binding to TrkB, the receptor for BDNF. Here we show that lysergic acid diethylamide (LSD) and psilocin directly bind to TrkB with affinities 1,000-fold higher than those for other antidepressants, and that psychedelics and antidepressants bind to distinct but partially overlapping sites within the transmembrane domain of TrkB dimers. The effects of psychedelics on neurotrophic signaling, plasticity and antidepressant-like behavior in mice depend on TrkB binding and promotion of endogenous BDNF signaling but are independent of serotonin 2A receptor (5-HT2A) activation, whereas LSD-induced head twitching is dependent on 5-HT2A and independent of TrkB binding. Our data confirm TrkB as a common primary target for antidepressants and suggest that high-affinity TrkB positive allosteric modulators lacking 5-HT2A activity may retain the antidepressant potential of psychedelics without hallucinogenic effects.

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

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

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.

Mechanistic Insight into How PEGylation Reduces the Efficacy of pH-Sensitive Liposomes from Molecular Dynamics Simulations

Liposome-based drug delivery systems composed of DOPE stabilized with cholesteryl hemisuccinate (CHMS) have been proposed as a drug delivery mechanism with pH-triggered release as the anionic form (CHSa) is protonated (CHS) at reduced pH; PEGylation is known to decrease this pH sensitivity. In this manuscript, we set out to use molecular dynamics (MD) simulations with a model with all-atom resolution to provide insight into why incorporation of poly(ethyleneglycol) (PEG) into DOPE–CHMS liposomes reduces their pH sensitivity; we also address two additional questions: (1) How CHSa stabilizes DOPE bilayers into a lamellar conformation at a physiological pH of 7.4? and (2) how the change from CHSa to CHS at acidic pH triggers the destabilization of DOPE bilayers? We found that (A) CHSa stabilizes the DOPE lipid membrane by increasing the hydrophilicity of the bilayer surface, (B) when CHSa changes to CHS by pH reduction, DOPE bilayers are destabilized due to a reduction in bilayer hydrophilicity and a reduction in the area per lipid, and (C) PEG stabilizes DOPE bilayers into the lamellar phase, thus reducing the pH sensitivity of the liposomes by increasing the area per lipid through penetration into the bilayer, which is our main focus.

Cholesterol-recognition motifs in the transmembrane domain of the tyrosine kinase receptor family: The case of TRKB

Cholesterol is an essential constituent of cell membranes. The discovery of cholesterol-recognition amino acid consensus (CRAC) motif in proteins indicated a putative direct, non-covalent interaction between cholesterol and proteins. In the present study, we evaluated the presence of a CRAC motif and its inverted version (CARC) in the transmembrane region (TMR) of the tyrosine kinase receptor family (RTK) in several species using in silico methods. CRAC motifs were found across all species analyzed, while CARC was found only in vertebrates. The tropomyosin-related kinase B (TRKB), a member of the RTK family, through interaction with its endogenous ligand brain-derived neurotrophic factor (BDNF) is a core participant in the neuronal plasticity process and exhibits a CARC motif in its TMR. Upon identifying the conserved CARC motif in the TRKB, we performed molecular dynamics simulations of the mouse TRKB.TMR. The simulations indicated that cholesterol interaction with the TRKB CARC motif occurs mainly at the central Y433 residue. Our binding assay suggested a bell-shaped effect of cholesterol on BDNF interaction with TRKB receptors, and our results suggest that CARC/CRAC motifs may play a role in the function of the RTK family TMR.

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.

Rigorous Computational Study Reveals What Docking Overlooks: Double Trouble from Membrane Association in Protein Kinase C Modulators

Increasing protein kinase C (PKC) activity is of potential therapeutic value. Its activation involves an interaction between the C1 domain and diacylglycerol (DAG) at intracellular membrane surfaces; DAG mimetics hold promise as new drugs. We previously developed the isophthalate derivative HMI-1a3, an effective but highly lipophilic (clogP = 6.46) DAG mimetic. Although a less lipophilic pyrimidine analog, PYR-1gP (clogP = 3.30), gave positive results in computational docking, it unexpectedly presented greatly diminished binding to PKC in vitro. Through more rigorous computational molecular modeling, we reveal that, unlike HMI-1a3, PYR-1gP forms an intramolecular hydrogen bond, which both obstructs binding and reorients PYR-1gP in the membrane in a fashion that prevents it from correctly accessing the PKC C1 domain. Our results highlight the great value of molecular dynamics simulations as a key component for the drug design process of ligands targeting weakly membrane-associated proteins, where simulation in the relevant membrane environment is crucial for obtaining biologically applicable results.

The F1 loop of the talin head domain acts as a gatekeeper in integrin activation and clustering

Integrin activation and clustering by talin are early steps of cell adhesion. Membrane-bound talin head domain and kindlin bind to the β integrin cytoplasmic tail, cooperating to activate the heterodimeric integrin, and the talin head domain induces integrin clustering in the presence of Mn2+ Here we show that kindlin-1 can replace Mn2+ to mediate β3 integrin clustering induced by the talin head, but not that induced by the F2-F3 fragment of talin. Integrin clustering mediated by kindlin-1 and the talin head was lost upon deletion of the flexible loop within the talin head F1 subdomain. Further mutagenesis identified hydrophobic and acidic motifs in the F1 loop responsible for β3 integrin clustering. Modeling, computational and cysteine crosslinking studies showed direct and catalytic interactions of the acidic F1 loop motif with the juxtamembrane domains of α- and β3-integrins, in order to activate the β3 integrin heterodimer, further detailing the mechanism by which the talin-kindlin complex activates and clusters integrins. Moreover, the F1 loop interaction with the β3 integrin tail required the newly identified compact FERM fold of the talin head, which positions the F1 loop next to the inner membrane clasp of the talin-bound integrin heterodimer.This article has an associated First Person interview with the first author of the paper.

Cooperative Effects of an Antifungal Moiety and DMSO on Pore Formation over Lipid Membranes Revealed by Free Energy Calculations

Itraconazole is a triazole drug widely used in the treatment of fungal infections, and it is in clinical trials for treatment of several cancers. However, the drug suffers from poor solubility, while experiments have shown that itraconazole delivery in liposome nanocarriers improves both circulation half-life and tissue distribution. The drug release mechanism from the nanocarrier is still unknown, and it depends on several factors including membrane stability against defect formation. In this work, we used molecular dynamics simulations and potential of mean force (PMF) calculations to quantify the influence of itraconazole on pore formation over lipid membranes, and we compared the effect by itraconazole with a pore-stabilizing effect by the organic solvent dimethyl sulfoxide (DMSO). According to the PMFs, both itraconazole and DMSO greatly reduce the free energy of pore formation, by up to ∼20 kJ mol^-1. However, whereas large concentrations of itraconazole of 8 mol % (relative to lipid) were required, only small concentrations of a few mole % DMSO (relative to water) were sufficient to stabilize pores. In addition, itraconazole and DMSO facilitate pore formation by different mechanisms. Whereas itraconazole predominantly aids the formation of a partial defect with a locally thinned membrane, DMSO mainly stabilizes a transmembrane water needle by shielding it from the hydrophobic core. Notably, the two distinct mechanisms act cooperatively upon adding both itraconazole and DMSO to the membrane, as revealed by an additional reduction of the pore free energy. Overall, our simulations reveal molecular mechanisms and free energies of membrane pore formation by small molecules. We suggest that the stabilization of a locally thinned membrane as well as the shielding of a transmembrane water needle from the hydrophobic membrane core may be a general mechanism by which amphiphilic molecules facilitate pore formation over lipid membranes at sufficient concentrations.

Membrane-Dependent Binding and Entry Mechanism of Dopamine into Its Receptor

Synaptic neurotransmission has recently been proposed to function via either a membrane-independent or a membrane-dependent mechanism, depending on the neurotransmitter type. In the membrane-dependent mechanism, amphipathic neurotransmitters first partition to the lipid headgroup region and then diffuse along the membrane plane to their membrane-buried receptors. However, to date, this mechanism has not been demonstrated for any neurotransmitter-receptor complex. Here, we combined isothermal calorimetry measurements with a diverse set of molecular dynamics simulation methods to investigate the partitioning of an amphipathic neurotransmitter (dopamine) and the mechanism of its entry into the ligand-binding site. Our results show that the binding of dopamine to its receptor is consistent with the membrane-dependent binding and entry mechanism. Both experimental and simulation results showed that dopamine favors binding to lipid membranes especially in the headgroup region. Moreover, our simulations revealed a ligand-entry pathway from the membrane to the binding site. This pathway passes through a lateral gate between transmembrane alpha-helices 5 and 6 on the membrane-facing side of the protein. All in all, our results demonstrate that dopamine binds to its receptor by a membrane-dependent mechanism, and this is complemented by the more traditional binding mechanism directly through the aqueous phase. The results suggest that the membrane-dependent mechanism is common in other synaptic receptors, too.

Syndecan-4 tunes cell mechanics by activating the kindlin-integrin-RhoA pathway

Extensive research over the past decades has identified integrins to be the primary transmembrane receptors that enable cells to respond to external mechanical cues. We reveal here a mechanism whereby syndecan-4 tunes cell mechanics in response to localized tension via a coordinated mechanochemical signalling response that involves activation of two other receptors: epidermal growth factor receptor and β1 integrin. Tension on syndecan-4 induces cell-wide activation of the kindlin-2/β1 integrin/RhoA axis in a PI3K-dependent manner. Furthermore, syndecan-4-mediated tension at the cell-extracellular matrix interface is required for yes-associated protein activation. Extracellular tension on syndecan-4 triggers a conformational change in the cytoplasmic domain, the variable region of which is indispensable for the mechanical adaptation to force, facilitating the assembly of a syndecan-4/α-actinin/F-actin molecular scaffold at the bead adhesion. This mechanotransduction pathway for syndecan-4 should have immediate implications for the broader field of mechanobiology.

Cholesterol Reduces Partitioning of Antifungal Drug Itraconazole into Lipid Bilayers

Cholesterol plays a crucial role in modulating the physicochemical properties of biomembranes, both increasing mechanical strength and decreasing permeability. Cholesterol is also a common component of vesicle-based delivery systems, including liposome-based drug delivery systems (LDSs). However, its effect on the partitioning of drug molecules to lipid membranes is very poorly recognized. Herein, we performed a combined experimental/computational study of the potential for the use of the LDS formulation for the delivery of the antifungal drug itraconazole (ITZ). We consider the addition of cholesterol to the lipid membrane. Since ITZ is only weakly soluble in water, its bioavailability is limited. Use of an LDS has thus been proposed. We studied lipid membranes composed of cholesterol, 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphocholine (POPC), and ITZ using a combination of computational molecular dynamics (MD) simulations of lipid bilayers and Brewster angle microscopy (BAM) experiments of monolayers. Both experimental and computational results show separation of cholesterol and ITZ. Cholesterol has a strong preference to orient parallel to the bilayer normal. However, ITZ, a long and relatively rigid molecule with weakly hydrophilic groups along the backbone, predominantly locates below the interface between the hydrocarbon chain region and the polar region of the membrane, with its backbone oriented parallel to the membrane surface; the orthogonal orientation in the membrane could be the cause of the observed separation. In addition, fluorescence measurements demonstrated that the affinity of ITZ for the lipid membrane is decreased by the presence of cholesterol, which is thus probably not a suitable formulation component of an LDS designed for ITZ delivery.

Behavior of the DPH fluorescence probe in membranes perturbed by drugs

1,6-Diphenyl-1,3,5-hexatriene (DPH) is one of the most commonly used fluorescent probes to study dynamical and structural properties of lipid bilayers and cellular membranes via measuring steady-state or time-resolved fluorescence anisotropy. In this study, we present a limitation in the use of DPH to predict the order of lipid acyl chains when the lipid bilayer is doped with itraconazole (ITZ), an antifungal drug. Our steady-state fluorescence anisotropy measurements showed a significant decrease in fluorescence anisotropy of DPH embedded in the ITZ-containing membrane, suggesting a substantial increase in membrane fluidity, which indirectly indicates a decrease in the order of the hydrocarbon chains. This result or its interpretation is in disagreement with the fluorescence recovery after photobleaching measurements and molecular dynamics (MD) simulation data. The results of these experiments and calculations indicate an increase in the hydrocarbon chain order. The MD simulations of the bilayer containing both ITZ and DPH provide explanations for these observations. Apparently, in the presence of the drug, the DPH molecules are pushed deeper into the hydrophobic membrane core below the lipid double bonds, and the probe predominately adopts the orientation of the ITZ molecules that is parallel to the membrane surface, instead of orienting parallel to the lipid acyl chains. For this reason, DPH anisotropy provides information related to the less ordered central region of the membrane rather than reporting the properties of the upper segments of the lipid acyl chains.

Complex Behavior of Phosphatidylcholine-Phosphatidic Acid Bilayers and Monolayers: Effect of Acyl Chain Unsaturation

Phosphatidic acids (PAs) have many biological functions in biomembranes, e.g., they are involved in the proliferation, differentiation, and transformation of cells. Despite decades of research, the molecular understanding of how PAs affect the properties of biomembranes remains elusive. In this study, we explored the properties of lipid bilayers and monolayers composed of PAs and phosphatidylcholines (PCs) with various acyl chains. For this purpose, the Langmuir monolayer technique and atomistic molecular dynamics (MD) simulations were used to study the miscibility of PA and PC lipids and the molecular organization of mixed bilayers. The monolayer experiments demonstrated that the miscibility of membrane components strongly depends on the structure of the hydrocarbon chains and thus on the overall lipid shape. Interactions between PA and PC molecules vary from repulsive, for systems containing lipids with saturated and unsaturated acyl tails (strongly positive values of the excess free energy of mixing), to attractive, for systems in which all lipid tails are saturated (negative values of the excess free energy of mixing). The MD simulations provided atomistic insight into polar interactions (formation of hydrogen bonds and charge pairs) in PC-PA systems. H-bonding between PA monoanions and PCs in mixed bilayers is infrequent, and the lipid molecules interact mainly via electrostatic interactions. However, the number of charge pairs significantly decreases with the number of unsaturated lipid chains in the PA-PC system. The PA dianions weakly interact with the zwitterionic lipids, but their headgroups are more hydrated as compared to the monoanionic form. The acyl chains in all PC-PA bilayers are more ordered compared to single-component PC systems. In addition, depending on the combination of lipids, we observed a deeper location of the PA phosphate groups compared to the PC phosphate groups, which can alter the presentation of PAs for the peripheral membrane proteins, affecting their accessibility for binding.

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.

Membrane bound COMT isoform is an interfacial enzyme: general mechanism and new drug design paradigm

We have determined the substrate differentiation mechanism between the membrane bound and water soluble isoforms of important drug target catechol-O-methyltransferase.

doGlycans-Tools for Preparing Carbohydrate Structures for Atomistic Simulations of Glycoproteins, Glycolipids, and Carbohydrate Polymers for GROMACS

Carbohydrates constitute a structurally and functionally diverse group of biological molecules and macromolecules. In cells they are involved in, e.g., energy storage, signaling, and cell–cell recognition. All of these phenomena take place in atomistic scales, thus atomistic simulation would be the method of choice to explore how carbohydrates function. However, the progress in the field is limited by the lack of appropriate tools for preparing carbohydrate structures and related topology files for the simulation models. Here we present tools that fill this gap. Applications where the tools discussed in this paper are particularly useful include, among others, the preparation of structures for glycolipids, nanocellulose, and glycans linked to glycoproteins. The molecular structures and simulation files generated by the tools are compatible with GROMACS.

Dynamics and energetics of the mammalian phosphatidylinositol transfer protein phospholipid exchange cycle

Phosphatidylinositol-transfer proteins (PITPs) regulate phosphoinositide signaling in eukaryotic cells. The defining feature of PITPs is their ability to exchange phosphatidylinositol (PtdIns) molecules between membranes, and this property is central to PITP-mediated regulation of lipid signaling. However, the details of the PITP-mediated lipid exchange cycle remain entirely obscure. Here, all-atom molecular dynamics simulations of the mammalian StART-like PtdIns/phosphatidylcholine (PtdCho) transfer protein PITPα, both on membrane bilayers and in solvated systems, informed downstream biochemical analyses that tested key aspects of the hypotheses generated by the molecular dynamics simulations. These studies provided five key insights into the PITPα lipid exchange cycle: (i) interaction of PITPα with the membrane is spontaneous and mediated by four specific protein substructures; (ii) the ability of PITPα to initiate closure around the PtdCho ligand is accompanied by loss of flexibility of two helix/loop regions, as well as of the C-terminal helix; (iii) the energy barrier of phospholipid extraction from the membrane is lowered by a network of hydrogen bonds between the lipid molecule and PITPα; (iv) the trajectory of PtdIns or PtdCho into and through the lipid-binding pocket is chaperoned by sets of PITPα residues conserved throughout the StART-like PITP family; and (v) conformational transitions in the C-terminal helix have specific functional involvements in PtdIns transfer activity. Taken together, these findings provide the first mechanistic description of key aspects of the PITPα PtdIns/PtdCho exchange cycle and offer a rationale for the high conservation of particular sets of residues across evolutionarily distant members of the metazoan StART-like PITP family.

Calcium Assists Dopamine Release by Preventing Aggregation on the Inner Leaflet of Presynaptic Vesicles

In this study, the dopamine-lipid bilayer interactions were probed with three physiologically relevant ion compositions using atomistic molecular dynamics simulations and free energy calculations. The in silico results indicate that calcium is able to decrease significantly the binding of dopamine to a neutral (zwitterionic) phosphatidylcholine lipid bilayer model mimicking the inner leaflet of a presynaptic vesicle. We argue that the observed calcium-induced effect is likely in crucial role in the neurotransmitter release from the presynaptic vesicles docked in the active zone of nerve terminals. The inner leaflets of presynaptic vesicles, which are responsible for releasing neurotransmitters into the synaptic cleft, are mainly composed of neutral lipids such as phosphatidylcholine and phosphatidylethanolamine. The neutrality of the lipid head group region, enhanced by a low pH level, should limit membrane aggregation of transmitters. In addition, the simulations suggest that the high calcium levels inside presynaptic vesicles prevent even the most lipophilic transmitters such as dopamine from adhering to the inner leaflet surface, thus rendering unhindered neurotransmitter release feasible.

Phase Partitioning of GM1 and Its Bodipy-Labeled Analog Determine Their Different Binding to Cholera Toxin

Driven by interactions between lipids and proteins, biological membranes display lateral heterogeneity that manifests itself in a mosaic of liquid-ordered (Lo) or raft, and liquid-disordered (Ld) or non-raft domains with a wide range of different properties and compositions. In giant plasma membrane vesicles and giant unilamellar vesicles, specific binding of Cholera Toxin (CTxB) to GM1 glycolipids is a commonly used strategy to label raft domains or Lo membrane environments. However, these studies often use acyl-chain labeled bodipy-GM1 (bdGM1), whose headgroup accessibility and membrane order or phase partitioning may differ from those of GM1, rendering the interpretation of CTxB binding data quite problematic. To unravel the molecular basis of CTxB binding to GM1 and bdGM1, we explored the partitioning and the headgroup presentation of these gangliosides in the Lo and Ld phases using atomistic molecular dynamics simulations complemented by CTxB binding experiments. The conformation of both GM1 and bdGM1 was shown to be largely similar in the Lo and Ld phases. However, bdGM1 showed reduction in receptor availability when reconstituted into synthetic bilayer mixtures, highlighting that membrane phase partitioning of the gangliosides plays a considerable role in CTxB binding. Our results suggest that the CTxB binding is predominately modulated by the partitioning of the receptor to an appropriate membrane phase. Further, given that the Lo and Ld partitioning of bdGM1 differs from those of GM1, usage of bdGM1 for studying GM1 behavior in cells can lead to invalid interpretation of experimental data.

A computational study suggests that replacing PEG with PMOZ may increase exposure of hydrophobic targeting moiety

In a previous study we showed that the cause of failure of a new, proposed, targeting ligand, the AETP moiety, when attached to a PEGylated liposome, was occlusion by the poly(ethylene glycol) (PEG) layer due to its hydrophobic nature, given that PEG is not entirely hydrophilic. At the time we proposed that possible replacement with a more hydrophilic protective polymer could alleviate this problem. In this study we have used computational molecular dynamics modelling, using a model with all atom resolution, to suggest that a specific alternative protective polymer, poly(2-methyloxazoline) (PMOZ), would perform exactly this function. Our results show that when PEG is replaced by PMOZ the relative exposure to the solvent of AETP is increased to a level even greater than that we found in previous simulations for the RGD peptide, a targeting moiety that has previously been used successfully in PEGylated liposome based therapies. While the AETP moiety itself is no longer under consideration, the results of this computational study have broader significance: the use of PMOZ as an alternative polymer coating to PEG could be efficacious in the context of more hydrophobic targeting ligands. In addition to PMOZ we studied another polyoxazoline, poly(2-ethyloxazoline) (PEOZ), that has also been mooted as a possible alternate protective polymer. It was also found that the RDG peptide occlusion was significantly greater for the case of both oxazolines as opposed to PEG and that, unlike PEG, neither oxazoline entered the membrane. As far as we are aware this is the first time that polyoxazolines have been studied using molecular dynamics simulation with all atom resolution.

Effects of Membrane PEGylation on Entry and Location of Antifungal Drug Itraconazole and Their Pharmacological Implications

Itraconazole (ITZ) is an antifungal agent used clinically to treat mycotic infections. However, its therapeutic effects are limited by low solubility in aqueous media. Liposome-based delivery systems (LDS) have been proposed as a delivery mechanism for ITZ to alleviate this problem. Furthermore, PEGylation, the inclusion in the formulation of a protective "stealth sheath" of poly(ethylene glycol) around carrier particles, is widely used to increase circulation time in the bloodstream and hence efficacy. Together, these themes highlight the importance of mechanistic and structural understanding of ITZ incorporation into liposomes both with and without PEGylation because it can provide a potential foundation for the rational design of LDS-based systems for delivery of ITZ, using alternate protective polymers or formulations. Here we have combined atomistic simulations, cryo-TEM, Langmuir film balance, and fluorescence quenching experiments to explore how ITZ interacts with both pristine and PEGylated liposomes. We found that the drug can be incorporated into conventional and PEGylated liposomes for drug concentrations up to 15 mol % without phase separation. We observed that, in addition to its protective properties, PEGylation significantly increases the stability of liposomes that host ITZ. In a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer without PEGylation, ITZ was found to reside inside the lipid bilayer between the glycerol and the double-bond regions of POPC, adopting a largely parallel orientation along the membrane surface. In a PEGylated liposome, ITZ partitions mainly to the PEG layer. The results provide a solid basis for further development of liposome-based delivery systems.

Calcium Directly Regulates Phosphatidylinositol 4,5-Bisphosphate Headgroup Conformation and Recognition

The orchestrated recognition of phosphoinositides and concomitant intracellular release of Ca²⁺ is pivotal to almost every aspect of cellular processes, including membrane homeostasis, cell division and growth, vesicle trafficking, as well as secretion. Although Ca²⁺ is known to directly impact phosphoinositide clustering, little is known about the molecular basis for this or its significance in cellular signaling. Here, we study the direct interaction of Ca²⁺ with phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂), the main lipid marker of the plasma membrane. Electrokinetic potential measurements of PI(4,5)P₂ containing liposomes reveal that Ca²⁺ as well as Mg²⁺ reduce the zeta potential of liposomes to nearly background levels of pure phosphatidylcholine membranes. Strikingly, lipid recognition by the default PI(4,5)P₂ lipid sensor, phospholipase C delta 1 pleckstrin homology domain (PLC δ1-PH), is completely inhibited in the presence of Ca²⁺, while Mg²⁺ has no effect with 100 nm liposomes and modest effect with giant unilamellar vesicles. Consistent with biochemical data, vibrational sum frequency spectroscopy and atomistic molecular dynamics simulations reveal how Ca²⁺ binding to the PI(4,5)P₂ headgroup and carbonyl regions leads to confined lipid headgroup tilting and conformational rearrangements. We rationalize these findings by the ability of calcium to block a highly specific interaction between PLC δ1-PH and PI(4,5)P₂, encoded within the conformational properties of the lipid itself. Our studies demonstrate the possibility that switchable phosphoinositide conformational states can serve as lipid recognition and controlled cell signaling mechanisms.

Atomistic Molecular Dynamics Simulations of Mitochondrial DNA Polymerase γ: Novel Mechanisms of Function and Pathogenesis

DNA polymerase γ (Pol γ) is a key component of the mitochondrial DNA replisome and an important cause of neurological diseases. Despite the availability of its crystal structures, the molecular mechanism of DNA replication, the switch between polymerase and exonuclease activities, the site of replisomal interactions, and functional effects of patient mutations that do not affect direct catalysis have remained elusive. Here we report the first atomistic classical molecular dynamics simulations of the human Pol γ replicative complex. Our simulation data show that DNA binding triggers remarkable changes in the enzyme structure, including (1) completion of the DNA-binding channel via a dynamic subdomain, which in the apo form blocks the catalytic site, (2) stabilization of the structure through the distal accessory β-subunit, and (3) formation of a putative transient replisome-binding platform in the "intrinsic processivity" subdomain of the enzyme. Our data indicate that noncatalytic mutations may disrupt replisomal interactions, thereby causing Pol γ-associated neurodegenerative disorders.

Glycosylation and Lipids Working in Concert Direct CD2 Ectodomain Orientation and Presentation

Proteins embedded in the plasma membrane mediate interactions with the cell environment and play decisive roles in many signaling events. For cell-cell recognition molecules, it is highly likely that their structures and behavior have been optimized in ways that overcome the limitations of membrane tethering. In particular, the ligand binding regions of these proteins likely need to be maximally exposed. Here we show by means of atomistic simulations of membrane-bound CD2, a small cell adhesion receptor expressed by human T-cells and natural killer cells, that the presentation of its ectodomain is highly dependent on membrane lipids and receptor glycosylation acting in apparent unison. Detailed analysis shows that the underlying mechanism is based on electrostatic interactions complemented by steric interactions between glycans in the protein and the membrane surface. The findings are significant for understanding the factors that render membrane receptors accessible for binding and signaling.

Identifying involvement of Lys251/Asp252 pair in electron transfer and associated proton transfer at the quinone reduction site of Rhodobacter capsulatus cytochrome bc1

Describing dynamics of proton transfers in proteins is challenging, but crucial for understanding processes which use them for biological functions. In cytochrome bc1, one of the key enzymes of respiration or photosynthesis, proton transfers engage in oxidation of quinol (QH2) and reduction of quinone (Q) taking place at two distinct catalytic sites. Here we evaluated by site-directed mutagenesis the contribution of Lys251/Asp252 pair (bacterial numbering) in electron transfers and associated with it proton uptake to the quinone reduction site (Qi site). We showed that the absence of protonable group at position 251 or 252 significantly changes the equilibrium levels of electronic reactions including the Qi-site mediated oxidation of heme bH, reverse reduction of heme bH by quinol and heme bH/Qi semiquinone equilibrium. This implicates the role of H-bonding network in binding of quinone/semiquinone and defining thermodynamic properties of Q/SQ/QH2 triad. The Lys251/Asp252 proton path is disabled only when both protonable groups are removed. With just one protonable residue from this pair, the entrance of protons to the catalytic site is sustained, albeit at lower rates, indicating that protons can travel through parallel routes, possibly involving water molecules. This shows that proton paths display engineering tolerance for change as long as all the elements available for functional cooperation secure efficient proton delivery to the catalytic site.

Distribution and dynamics of quinones in the lipid bilayer mimicking the inner membrane of mitochondria

Quinone and its analogues (Q) constitute an important class of compounds that perform key electron transfer reactions in oxidative- and photo-phosphorylation. In the inner membrane of mitochondria, ubiquinone molecules undergo continuous redox transitions enabling electron transfer between the respiratory complexes. In such a dynamic system undergoing continuous turnover for ATP synthesis, an uninterrupted supply of substrate molecules is absolutely necessary. In the current work, we have performed atomistic molecular dynamics simulations and free energy calculations to assess the structure, dynamics, and localization of quinone and its analogues in a lipid bilayer, whose composition mimics the one in the inner mitochondrial membrane. The results show that there is a strong tendency of both quinone and quinol molecules to localize in the vicinity of the lipids' acyl groups, right under the lipid head group region. Additionally, we observe a second location in the middle of the bilayer where quinone molecules tend to stabilize. Translocation of quinone through a lipid bilayer is very fast and occurs in 10-100ns time scale, whereas the translocation of quinol is at least an order of magnitude slower. We suggest that this has important mechanistic implications given that the localization of Q ensures maximal occupancy of the Q-binding sites or Q-entry points in electron transport chain complexes, thereby maintaining an optimal turnover rate for ATP synthesis.

Design of cholesterol arabinogalactan anchored liposomes for asialoglycoprotein receptor mediated targeting to hepatocellular carcinoma: In silico modeling, in vitro and in vivo evaluation

We have developed active targeting liposomes to deliver anticancer agents to ASGPR which will contribute to effective treatment of hepatocellular carcinoma. Active targeting is achieved through polymeric ligands on the liposome surface. The liposomes were prepared using reverse phase evaporation method and doxorubicin hydrocholoride, a model drug, was loaded using the ammonium sulphate gradient method. Liposomes loaded with DOX were found to have a particle size of 200nm with more than 90% entrapment efficiency. Systems were observed to release the drug in a sustained manner in acidic pH in vitro. Liposomes containing targeting ligands possessed greater and selective toxicity to ASGPR positive HepG2 cell lines due to specific ligand receptor interaction. Bio-distribution studies revealed that liposomes were concentrated in the liver even after 3h of administration, thus providing conclusive evidence of targeting potential for formulated nanosystems. Tumor regression studies indicated greater tumor suppression with targeted liposomes thereby establishing superiority of the liposomal system. In this work, we used a novel methodology to guide the determination of the optimal composition of the targeting liposomes: molecular dynamics (MD) simulation that aided our understanding of the behaviour of the ligand within the bilayer. This can be seen as a demonstration of the utility of this methodology as a rational design tool for active targeting liposome formulation.

Data including GROMACS input files for atomistic molecular dynamics simulations of mixed, asymmetric bilayers including molecular topologies, equilibrated structures, and force field for lipids compatible with OPLS-AA parameters

In this Data in Brief article we provide a data package of GROMACS input files for atomistic molecular dynamics simulations of multicomponent, asymmetric lipid bilayers using the OPLS-AA force field. These data include 14 model bilayers composed of 8 different lipid molecules. The lipids present in these models are: cholesterol (CHOL), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylethanolamine (POPE), 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidyl-ethanolamine (SOPE), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylserine (POPS), 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidylserine (SOPS), N-palmitoyl-D-erythro-sphingosyl-phosphatidylcholine (SM16), and N-lignoceroyl-D-erythro-sphingosyl-phosphatidylcholine (SM24). The bilayers׳ compositions are based on lipidomic studies of PC-3 prostate cancer cells and exosomes discussed in Llorente et al. (2013) [1], showing an increase in the section of long-tail lipid species (SOPS, SOPE, and SM24) in the exosomes. Former knowledge about lipid asymmetry in cell membranes was accounted for in the models, meaning that the model of the inner leaflet is composed of a mixture of PC, PS, PE, and cholesterol, while the extracellular leaflet is composed of SM, PC and cholesterol discussed in Van Meer et al. (2008) [2]. The provided data include lipids׳ topologies, equilibrated structures of asymmetric bilayers, all force field parameters, and input files with parameters describing simulation conditions (md.mdp). The data is associated with the research article “Interdigitation of Long-Chain Sphingomyelin Induces Coupling of Membrane Leaflets in a Cholesterol Dependent Manner” (Róg et al., 2016) [3].

Functionalized lipids and surfactants for specific applications

Synthetic lipids and surfactants that do not exist in biological systems have been used for the last few decades in both basic and applied science. The most notable applications for synthetic lipids and surfactants are drug delivery, gene transfection, as reporting molecules, and as support for structural lipid biology. In this review, we describe the potential of the synergistic combination of computational and experimental methodologies to study the behavior of synthetic lipids and surfactants embedded in lipid membranes and liposomes. We focused on select cases in which molecular dynamics simulations were used to complement experimental studies aiming to understand the structure and properties of new compounds at the atomistic level. We also describe cases in which molecular dynamics simulations were used to design new synthetic lipids and surfactants, as well as emerging fields for the application of these compounds. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Selective effect of cell membrane on synaptic neurotransmission

Atomistic molecular dynamics simulations were performed with 13 non-peptidic neurotransmitters (NTs) in three different membrane environments. The results provide compelling evidence that NTs are divided into membrane-binding and membrane-nonbinding molecules. NTs adhere to the postsynaptic membrane surface whenever the ligand-binding sites of their synaptic receptors are buried in the lipid bilayer. In contrast, NTs that have extracellular ligand-binding sites do not have a similar tendency to adhere to the membrane surface. This finding is a seemingly simple yet important addition to the paradigm of neurotransmission, essentially dividing it into membrane-independent and membrane-dependent mechanisms. Moreover, the simulations also indicate that the lipid composition especially in terms of charged lipids can affect the membrane partitioning of NTs. The revised paradigm, highlighting the importance of cell membrane and specific lipids for neurotransmission, should to be of interest to neuroscientists, drug industry and the general public alike.

Stearylated cycloarginine nanosystems for intracellular delivery – simulations, formulation and proof of concept

Novel cationic agent liposomes performed better in silico translating in higher cellular uptake with reduced toxicity.

Toward Atomistic Resolution Structure of Phosphatidylcholine Headgroup and Glycerol Backbone at Different Ambient Conditions

Phospholipids are essential building blocks of biological membranes. Despite a vast amount of very accurate experimental data, the atomistic resolution structures sampled by the glycerol backbone and choline headgroup in phoshatidylcholine bilayers are not known. Atomistic resolution molecular dynamics simulations have the potential to resolve the structures, and to give an arrestingly intuitive interpretation of the experimental data, but only if the simulations reproduce the data within experimental accuracy. In the present work, we simulated phosphatidylcholine (PC) lipid bilayers with 13 different atomistic models, and compared simulations with NMR experiments in terms of the highly structurally sensitive C-H bond vector order parameters. Focusing on the glycerol backbone and choline headgroups, we showed that the order parameter comparison can be used to judge the atomistic resolution structural accuracy of the models. Accurate models, in turn, allow molecular dynamics simulations to be used as an interpretation tool that translates these NMR data into a dynamic three-dimensional representation of biomolecules in biologically relevant conditions. In addition to lipid bilayers in fully hydrated conditions, we reviewed previous experimental data for dehydrated bilayers and cholesterol-containing bilayers, and interpreted them with simulations. Although none of the existing models reached experimental accuracy, by critically comparing them we were able to distill relevant chemical information: (1) increase of choline order parameters indicates the P-N vector tilting more parallel to the membrane, and (2) cholesterol induces only minor changes to the PC (glycerol backbone) structure. This work has been done as a fully open collaboration, using nmrlipids.blogspot.fi as a communication platform; all the scientific contributions were made publicly on this blog. During the open research process, the repository holding our simulation trajectories and files ( https://zenodo.org/collection/user-nmrlipids ) has become the most extensive publicly available collection of molecular dynamics simulation trajectories of lipid bilayers.

PIP2 and Talin Join Forces to Activate Integrin

Integrins are major players in cell adhesion and migration, and malfunctions in controlling their activity are associated with various diseases. Nevertheless, the details of integrin activation are not completely understood, and the role of lipids in the process is largely unknown. Herein, we show using atomistic molecular dynamics simulations that the interplay of phosphatidylinositol 4,5-bisphosphate (PIP2) and talin may directly alter the conformation of integrin αIIbβ3. Our results provide a new perspective on the role of PIP2 in integrin activation and indicate that the charged PIP2 lipid headgroup can perturb a clasp at the cytoplasmic face of the integrin heterodimer.

Building Synthetic Sterols Computationally - Unlocking the Secrets of Evolution?

Cholesterol is vital in regulating the physical properties of animal cell membranes. While it remains unclear what renders cholesterol so unique, it is known that other sterols are less capable in modulating membrane properties, and there are membrane proteins whose function is dependent on cholesterol. Practical applications of cholesterol include its use in liposomes in drug delivery and cosmetics, cholesterol-based detergents in membrane protein crystallography, its fluorescent analogs in studies of cholesterol transport in cells and tissues, etc. Clearly, in spite of their difficult synthesis, producing the synthetic analogs of cholesterol is of great commercial and scientific interest. In this article, we discuss how synthetic sterols non-existent in nature can be used to elucidate the roles of cholesterol’s structural elements. To this end, we discuss recent atomistic molecular dynamics simulation studies that have predicted new synthetic sterols with properties comparable to those of cholesterol. We also discuss more recent experimental studies that have vindicated these predictions. The paper highlights the strength of computational simulations in making predictions for synthetic biology, thereby guiding experiments.

Role of subunit III and its lipids in the molecular mechanism of cytochrome c oxidase

The terminal respiratory enzyme cytochrome c oxidase (CcO) reduces molecular oxygen to water, and pumps protons across the inner mitochondrial membrane, or the plasma membrane of bacteria. A two-subunit CcO harbors all the elements necessary for oxygen reduction and proton pumping. However, it rapidly undergoes turnover-induced irreversible damage, which is effectively prevented by the presence of subunit III and its tightly bound lipids. We have performed classical atomistic molecular dynamics (MD) simulations on a three-subunit CcO, which show the formation of water wires between the polar head groups of lipid molecules bound to subunit III and the proton uptake site Asp91 (Bos taurus enzyme numbering). Continuum electrostatic calculations suggest that these lipids directly influence the proton affinity of Asp91 by 1-2pK units. We surmise that lipids bound to subunit III influence the rate of proton uptake through the D-pathway, and therefore play a key role in preventing turnover-induced inactivation. Atomistic MD simulations show that subunit III is rapidly hydrated in the absence of internally bound lipids, which is likely to affect the rate of O2 diffusion into the active-site. The role of subunit III with its indigenous lipids in the molecular mechanism of CcO is discussed.

Sec14-nodulin proteins and the patterning of phosphoinositide landmarks for developmental control of membrane morphogenesis

A Sec14-nodulin protein model is used to identify the nodulin domain as a novel phosphoinositide effector module with a role in controlling lateral organization of phosphoinositide. The domain organization of Sec14-nodulin proteins suggests a versatile principle for the bit mapping of membrane surfaces into high-definition lipid-signaling screens.

Polarized membrane morphogenesis is a fundamental activity of eukaryotic cells. This process is essential for the biology of cells and tissues, and its execution demands exquisite temporal coordination of functionally diverse membrane signaling reactions with high spatial resolution. Moreover, mechanisms must exist to establish and preserve such organization in the face of randomizing forces that would diffuse it. Here we identify the conserved AtSfh1 Sec14-nodulin protein as a novel effector of phosphoinositide signaling in the extreme polarized membrane growth program exhibited by growing Arabidopsis root hairs. The data are consistent with Sec14-nodulin proteins controlling the lateral organization of phosphatidylinositol 4,5-bisphosphate (PtdIns(4,5)P₂) landmarks for polarized membrane morphogenesis in plants. This patterning activity requires both the PtdIns(4,5)P₂ binding and homo-oligomerization activities of the AtSfh1 nodulin domain and is an essential aspect of the polarity signaling program in root hairs. Finally, the data suggest a general principle for how the phosphoinositide signaling landscape is physically bit mapped so that eukaryotic cells are able to convert a membrane surface into a high-definition lipid-signaling screen.

Co-exposure with fullerene may strengthen health effects of organic industrial chemicals

In vitro toxicological studies together with atomistic molecular dynamics simulations show that occupational co-exposure with C₆₀ fullerene may strengthen the health effects of organic industrial chemicals. The chemicals studied are acetophenone, benzaldehyde, benzyl alcohol, m-cresol, and toluene which can be used with fullerene as reagents or solvents in industrial processes. Potential co-exposure scenarios include a fullerene dust and organic chemical vapor, or a fullerene solution aerosolized in workplace air. Unfiltered and filtered mixtures of C₆₀ and organic chemicals represent different co-exposure scenarios in in vitro studies where acute cytotoxicity and immunotoxicity of C₆₀ and organic chemicals are tested together and alone by using human THP-1-derived macrophages. Statistically significant co-effects are observed for an unfiltered mixture of benzaldehyde and C₆₀ that is more cytotoxic than benzaldehyde alone, and for a filtered mixture of m-cresol and C₆₀ that is slightly less cytotoxic than m-cresol. Hydrophobicity of chemicals correlates with co-effects when secretion of pro-inflammatory cytokines IL-1β and TNF-α is considered. Complementary atomistic molecular dynamics simulations reveal that C₆₀ co-aggregates with all chemicals in aqueous environment. Stable aggregates have a fullerene-rich core and a chemical-rich surface layer, and while essentially all C₆₀ molecules aggregate together, a portion of organic molecules remains in water.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.

Cholesterol level affects surface charge of lipid membranes in saline solution

Cholesterol is an important component of all biological membranes as well as drug delivery liposomes. We show here that increasing the level of cholesterol in a phospholipid membrane decreases surface charge in the physiological environment. Through molecular dynamics simulation we have shown that increasing the level of cholesterol decreases Na+ ion binding. Complementary experimental ζ--potential measurements have shown a decreased ζ--potential with increasing cholesterol content, indicative of reduced surface charge. Both experiments and simulations have been carried out on both saturated 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC) and monounsaturated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membranes. This result is particularly important because membrane surface charge plays an important role in the interactions of biomembranes with peripheral membrane proteins and drug delivery liposomes with the immune system.