1
|
Oliveira AC, Filipe HAL, Geraldes CF, Voth GA, Moreno MJ, Loura LMS. Interaction of MRI Contrast Agent [Gd(DOTA)] - with Lipid Membranes: A Molecular Dynamics Study. Inorg Chem 2024; 63:10897-10914. [PMID: 38795015 PMCID: PMC11186012 DOI: 10.1021/acs.inorgchem.4c00972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2024] [Revised: 04/16/2024] [Accepted: 04/30/2024] [Indexed: 05/27/2024]
Abstract
Contrast agents are important imaging probes in clinical MRI, allowing the identification of anatomic changes that otherwise would not be possible. Intensive research on the development of new contrast agents is being made to image specific pathological markers or sense local biochemical changes. The most widely used MRI contrast agents are based on gadolinium(III) complexes. Due to their very high charge density, they have low permeability through tight biological barriers such as the blood-brain barrier, hampering their application in the diagnosis of neurological disorders. In this study, we explore the interaction between the widely used contrast agent [Gd(DOTA)]- (Dotarem) and POPC lipid bilayers by means of molecular dynamics simulations. This metal complex is a standard reference where several chemical modifications have been introduced to improve key properties such as bioavailability and targeting. The simulations unveil detailed insights into the agent's interaction with the lipid bilayer, offering perspectives beyond experimental methods. Various properties, including the impact on global and local bilayer properties, were analyzed. As expected, the results indicate a low partition coefficient (KP) and high permeation barrier for this reference compound. Nevertheless, favorable interactions are established with the membrane leading to moderately long residence times. While coordination of one inner-sphere water molecule is maintained for the membrane-associated chelate, the physical-chemical attributes of [Gd(DOTA)]- as a MRI contrast agent are affected. Namely, increases in the rotational correlation times and in the residence time of the inner-sphere water are observed, with the former expected to significantly increase the water proton relaxivity. This work establishes a reference framework for the use of simulations to guide the rational design of new contrast agents with improved relaxivity and bioavailability and for the development of liposome-based formulations for use as imaging probes or theranostic agents.
Collapse
Affiliation(s)
- Alexandre C. Oliveira
- Coimbra
Chemistry Centre, Institute of Molecular
Sciences (CQC-IMS), 3004-535 Coimbra, Portugal
- Department
of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal
| | - Hugo A. L. Filipe
- Coimbra
Chemistry Centre, Institute of Molecular
Sciences (CQC-IMS), 3004-535 Coimbra, Portugal
- CPIRN-IPG—Center
of Potential and Innovation of Natural Resources, Polytechnic Institute of Guarda, 6300-559 Guarda, Portugal
| | - Carlos F.G.C. Geraldes
- Coimbra
Chemistry Centre, Institute of Molecular
Sciences (CQC-IMS), 3004-535 Coimbra, Portugal
- Department
of Life Sciences, University of Coimbra, Calçada Martim de Freitas, 3000-393 Coimbra, Portugal
- CIBIT/ICNAS
- Instituto de Ciências Nucleares Aplicadas à Saúde, Pólo das Ciências
da Saúde, Azinhaga de Santa Comba, 3000-548 Coimbra, Portugal
| | - Gregory A. Voth
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, James Franck
Institute, and Institute for Biophysical Dynamics, University of Chicago, Chicago, Illinois 60637, United States
| | - Maria João Moreno
- Coimbra
Chemistry Centre, Institute of Molecular
Sciences (CQC-IMS), 3004-535 Coimbra, Portugal
- Department
of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal
- CNC−Center
for Neuroscience and Cell Biology, University
of Coimbra, 3004-517 Coimbra, Portugal
| | - Luís M. S. Loura
- Coimbra
Chemistry Centre, Institute of Molecular
Sciences (CQC-IMS), 3004-535 Coimbra, Portugal
- Faculty
of Pharmacy, University of Coimbra, 3000-548 Coimbra, Portugal
- CNC−Center
for Neuroscience and Cell Biology, University
of Coimbra, 3004-517 Coimbra, Portugal
| |
Collapse
|
2
|
Steshin IS, Vasyankin AV, Shirokova EA, Rozhkov AV, Livshits GD, Panteleev SV, Radchenko EV, Ignatov SK, Palyulin VA. Free Energy Barriers for Passive Drug Transport through the Mycobacterium tuberculosis Outer Membrane: A Molecular Dynamics Study. Int J Mol Sci 2024; 25:1006. [PMID: 38256079 PMCID: PMC10815926 DOI: 10.3390/ijms25021006] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2023] [Revised: 01/04/2024] [Accepted: 01/09/2024] [Indexed: 01/24/2024] Open
Abstract
The emergence of multi-drug-resistant tuberculosis strains poses a significant challenge to modern medicine. The development of new antituberculosis drugs is hindered by the low permeability of many active compounds through the extremely strong bacterial cell wall of mycobacteria. In order to estimate the ability of potential antimycobacterial agents to diffuse through the outer mycolate membrane, the free energy profiles, the corresponding activation barriers, and possible permeability modes of passive transport for a series of known antibiotics, modern antituberculosis drugs, and prospective active drug-like molecules were determined using molecular dynamics simulations with the all-atom force field and potential of mean-force calculations. The membranes of different chemical and conformational compositions, density, thickness, and ionization states were examined. The typical activation barriers for the low-mass molecules penetrating through the most realistic membrane model were 6-13 kcal/mol for isoniazid, pyrazinamide, and etambutol, and 19 and 25 kcal/mol for bedaquilin and rifampicin. The barriers for the ionized molecules are usually in the range of 37-63 kcal/mol. The linear regression models were derived from the obtained data, allowing one to estimate the permeability barriers from simple physicochemical parameters of the diffusing molecules, notably lipophilicity and molecular polarizability.
Collapse
Affiliation(s)
- Ilya S. Steshin
- Department of Chemistry, Lobachevsky State University of Nizhny Novgorod, Nizhny Novgorod 603022, Russia; (I.S.S.); (A.V.V.); (E.A.S.); (A.V.R.); (G.D.L.); (S.V.P.); (E.V.R.)
| | - Alexander V. Vasyankin
- Department of Chemistry, Lobachevsky State University of Nizhny Novgorod, Nizhny Novgorod 603022, Russia; (I.S.S.); (A.V.V.); (E.A.S.); (A.V.R.); (G.D.L.); (S.V.P.); (E.V.R.)
| | - Ekaterina A. Shirokova
- Department of Chemistry, Lobachevsky State University of Nizhny Novgorod, Nizhny Novgorod 603022, Russia; (I.S.S.); (A.V.V.); (E.A.S.); (A.V.R.); (G.D.L.); (S.V.P.); (E.V.R.)
| | - Alexey V. Rozhkov
- Department of Chemistry, Lobachevsky State University of Nizhny Novgorod, Nizhny Novgorod 603022, Russia; (I.S.S.); (A.V.V.); (E.A.S.); (A.V.R.); (G.D.L.); (S.V.P.); (E.V.R.)
| | - Grigory D. Livshits
- Department of Chemistry, Lobachevsky State University of Nizhny Novgorod, Nizhny Novgorod 603022, Russia; (I.S.S.); (A.V.V.); (E.A.S.); (A.V.R.); (G.D.L.); (S.V.P.); (E.V.R.)
| | - Sergey V. Panteleev
- Department of Chemistry, Lobachevsky State University of Nizhny Novgorod, Nizhny Novgorod 603022, Russia; (I.S.S.); (A.V.V.); (E.A.S.); (A.V.R.); (G.D.L.); (S.V.P.); (E.V.R.)
| | - Eugene V. Radchenko
- Department of Chemistry, Lobachevsky State University of Nizhny Novgorod, Nizhny Novgorod 603022, Russia; (I.S.S.); (A.V.V.); (E.A.S.); (A.V.R.); (G.D.L.); (S.V.P.); (E.V.R.)
- Department of Chemistry, Lomonosov Moscow State University, Leninskie Gory 1/3, Moscow 119991, Russia
| | - Stanislav K. Ignatov
- Department of Chemistry, Lobachevsky State University of Nizhny Novgorod, Nizhny Novgorod 603022, Russia; (I.S.S.); (A.V.V.); (E.A.S.); (A.V.R.); (G.D.L.); (S.V.P.); (E.V.R.)
| | - Vladimir A. Palyulin
- Department of Chemistry, Lobachevsky State University of Nizhny Novgorod, Nizhny Novgorod 603022, Russia; (I.S.S.); (A.V.V.); (E.A.S.); (A.V.R.); (G.D.L.); (S.V.P.); (E.V.R.)
- Department of Chemistry, Lomonosov Moscow State University, Leninskie Gory 1/3, Moscow 119991, Russia
| |
Collapse
|
3
|
Filipe HAL, Loura LMS, Moreno MJ. Permeation of a Homologous Series of NBD-Labeled Fatty Amines through Lipid Bilayers: A Molecular Dynamics Study. MEMBRANES 2023; 13:551. [PMID: 37367755 DOI: 10.3390/membranes13060551] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/05/2023] [Revised: 05/18/2023] [Accepted: 05/20/2023] [Indexed: 06/28/2023]
Abstract
Permeation through biomembranes is ubiquitous for drugs to reach their active sites. Asymmetry of the cell plasma membrane (PM) has been described as having an important role in this process. Here we describe the interaction of a homologous series of 7-nitrobenz-2-oxa-1,3-diazol-4-yl (NBD)-labeled amphiphiles (NBD-Cn, n = 4 to 16) with lipid bilayers of different compositions (1-palmitoyl, 2-oleoyl-sn-glycero-3-phosphocholine (POPC):cholesterol (1:1) and palmitoylated sphingomyelin (SpM):cholesterol (6:4)), including an asymmetric bilayer. Both unrestrained and umbrella sampling (US) simulations (at varying distances to the bilayer center) were carried out. The free energy profile of NBD-Cn at different depths in the membrane was obtained from the US simulations. The behavior of the amphiphiles during the permeation process was described regarding their orientation, chain elongation, and H-bonding to lipid and water molecules. Permeability coefficients were also calculated for the different amphiphiles of the series, using the inhomogeneous solubility-diffusion model (ISDM). Quantitative agreement with values obtained from kinetic modeling of the permeation process could not be obtained. However, for the longer, and more hydrophobic amphiphiles, the variation trend along the homologous series was qualitatively better matched by the ISDM when the equilibrium location of each amphiphile was taken as reference (ΔG = 0), compared to the usual choice of bulk water.
Collapse
Affiliation(s)
- Hugo A L Filipe
- Coimbra Chemistry Center, Institute of Molecular Sciences (CQC-IMS), University of Coimbra, 3004-535 Coimbra, Portugal
- CPIRN-IPG-Center of Potential and Innovation of Natural Resources, Polytechnic Institute of Guarda, 6300-559 Guarda, Portugal
| | - Luís M S Loura
- Coimbra Chemistry Center, Institute of Molecular Sciences (CQC-IMS), University of Coimbra, 3004-535 Coimbra, Portugal
- CNC-Center for Neuroscience and Cell Biology, University of Coimbra, 3004-535 Coimbra, Portugal
- Faculty of Pharmacy, University of Coimbra, 3000-548 Coimbra, Portugal
| | - Maria João Moreno
- Coimbra Chemistry Center, Institute of Molecular Sciences (CQC-IMS), University of Coimbra, 3004-535 Coimbra, Portugal
- CNC-Center for Neuroscience and Cell Biology, University of Coimbra, 3004-535 Coimbra, Portugal
- Department of Chemistry, Faculty of Sciences and Technology, University of Coimbra, 3004-535 Coimbra, Portugal
| |
Collapse
|
4
|
Benmameri M, Chantemargue B, Humeau A, Trouillas P, Fabre G. MemCross: Accelerated Weight Histogram method to assess membrane permeability. BIOCHIMICA ET BIOPHYSICA ACTA. BIOMEMBRANES 2023; 1865:184120. [PMID: 36669638 DOI: 10.1016/j.bbamem.2023.184120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/28/2022] [Revised: 01/04/2023] [Accepted: 01/05/2023] [Indexed: 01/19/2023]
Abstract
Passive permeation events across biological membranes are determining steps in the pharmacokinetics of xenobiotics. To reach an accurate and rapid prediction of membrane permeation coefficients of drugs is a complex challenge, which can efficiently support drug discovery. Such predictions are indeed highly valuable as they may guide the selection of potential leads with optimum bioavailabilities prior to synthesis. Theoretical models exist to predict these coefficients. Many of them are based on molecular dynamics (MD) simulations, which allow calculation of permeation coefficients through the evaluation of both the potential of mean force (PMF) and the diffusivity profiles. However, these simulations still require intensive computational efforts, and novel methodologies should be developed and benchmarked. Free energy perturbation (FEP) method was recently shown to estimate PMF with a significantly reduced computational cost compared to the adaptive biasing force method. This benchmarking was achieved with small molecules, namely short-chain alcohols. Here, we show that to estimate the PMF of bulkier, drug-like xenobiotics, conformational sampling is a critical issue. To reach a sufficient sampling with FEP calculations requires a relatively long time-scale, which can lower the benefits related to the computational gain. In the present work, the Accelerated Weight Histogram (AWH) method was employed for the first time in all-atom membrane models. The AWH-based protocol, named MemCross, appears affordable to estimate PMF profiles of a series of drug-like xenobiotics, compared to other enhanced sampling methods. The continuous exploration of the crossing pathway by MemCross also allows modeling subdiffusion by computing fractional diffusivity profiles. The method is also versatile as its input parameters are largely insensitive to the molecule properties. It also ensures a detailed description of the molecule orientations along the permeation pathway, picturing all intermolecular interactions at an atomic resolution. Here, MemCross was applied on a series of 12 xenobiotics, including four weak acids, and a coherent structure-activity relationship was established.
Collapse
Affiliation(s)
| | | | | | - Patrick Trouillas
- INSERM, UMR 1248, F-87000 Limoges, France; CATRIN RCPTM, 779 00 Olomouc, Holice, Czech Republic
| | | |
Collapse
|
5
|
Benetis NP, Paloncýová M, Knippenberg S. Multiscale Modeling Unravels the Influence of Biomembranes on the Photochemical Properties of Embedded Anti-Oxidative Polyphenolic and Phenanthroline Chelating Dyes. J Phys Chem B 2023; 127:212-227. [PMID: 36563093 DOI: 10.1021/acs.jpcb.2c07072] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The embedding of caffeate methyl ester, the flavonoids luteolin and quercetin, and the o-phenanthroline and neocuproine in a liquid disordered lipid bilayer has been studied through extensive atomistic calculations. The location and the orientation of these bio-active antioxidants are explained and analyzed. While the two phenanthrolines strongly associate with the lipid tail region, the other three compounds are rather found among the head groups. The simulations showcase conformational changes of the flavonoids. Through the use of a hybrid quantum mechanics-molecular mechanics scheme and supported by a profound benchmarking of the electronic excited-state method for these compounds, the influence of the anisotropic environment on the compounds' optical properties is analyzed. Influences of surrounding water molecules and of the polar parts of the lipids on the transition dipole moments and excited-state dipole moments are weighted with respect to a change in conformation. The current study highlights the importance of the mapping of molecular interactions in model membranes and pinpoints properties, which can be biomedically used to discriminate and detect different lipid environments.
Collapse
Affiliation(s)
| | - Markéta Paloncýová
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technologies and Research Institute, Palacký University Olomouc, Křížkovského 8, Olomouc779 00, Czech Republic
| | - Stefan Knippenberg
- Hasselt University, Theory Lab, Agoralaan Building D, 3590Diepenbeek, Belgium.,Spectroscopy, Quantum Chemistry and Atmospheric Remote Sensing (SQUARES), Université Libre de Bruxelles, 50 Avenue F. Roosevelt, C.P. 160/09, B-1050Brussels, Belgium
| |
Collapse
|
6
|
Paloncýová M, Pykal M, Kührová P, Banáš P, Šponer J, Otyepka M. Computer Aided Development of Nucleic Acid Applications in Nanotechnologies. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2022; 18:e2204408. [PMID: 36216589 DOI: 10.1002/smll.202204408] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/23/2022] [Revised: 09/12/2022] [Indexed: 06/16/2023]
Abstract
Utilization of nucleic acids (NAs) in nanotechnologies and nanotechnology-related applications is a growing field with broad application potential, ranging from biosensing up to targeted cell delivery. Computer simulations are useful techniques that can aid design and speed up development in this field. This review focuses on computer simulations of hybrid nanomaterials composed of NAs and other components. Current state-of-the-art molecular dynamics simulations, empirical force fields (FFs), and coarse-grained approaches for the description of deoxyribonucleic acid and ribonucleic acid are critically discussed. Challenges in combining biomacromolecular and nanomaterial FFs are emphasized. Recent applications of simulations for modeling NAs and their interactions with nano- and biomaterials are overviewed in the fields of sensing applications, targeted delivery, and NA templated materials. Future perspectives of development are also highlighted.
Collapse
Affiliation(s)
- Markéta Paloncýová
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Martin Pykal
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Petra Kührová
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Pavel Banáš
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Jiří Šponer
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
- Institute of Biophysics of the Czech Academy of Sciences, v. v. i., Královopolská 135, Brno, 612 65, Czech Republic
| | - Michal Otyepka
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
- IT4Innovations, VŠB - Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba, 708 00, Czech Republic
| |
Collapse
|
7
|
Nencini R, Ollila OHS. Charged Small Molecule Binding to Membranes in MD Simulations Evaluated against NMR Experiments. J Phys Chem B 2022; 126:6955-6963. [PMID: 36063117 PMCID: PMC9483918 DOI: 10.1021/acs.jpcb.2c05024] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Interactions of charged molecules with biomembranes regulate many of their biological activities, but their binding affinities to lipid bilayers are difficult to measure experimentally and model theoretically. Classical molecular dynamics (MD) simulations have the potential to capture the complex interactions determining how charged biomolecules interact with membranes, but systematic overbinding of sodium and calcium cations in standard MD simulations raises the question of how accurately force fields capture the interactions between lipid membranes and charged biomolecules. Here, we evaluate the binding of positively charged small molecules, etidocaine, and tetraphenylphosphonium to a phosphatidylcholine (POPC) lipid bilayer using the changes in lipid head-group order parameters. We observed that these molecules behave oppositely to calcium and sodium ions when binding to membranes: (i) their binding affinities are not overestimated by standard force field parameters, (ii) implicit inclusion of electronic polarizability increases their binding affinity, and (iii) they penetrate into the hydrophobic membrane core. Our results can be explained by distinct binding mechanisms of charged small molecules with hydrophobic moieties and monoatomic ions. The binding of the former is driven by hydrophobic effects, while the latter has direct electrostatic interactions with lipids. In addition to elucidating how different kinds of charged biomolecules bind to membranes, we deliver tools for further development of MD simulation parameters and methodology.
Collapse
Affiliation(s)
- Ricky Nencini
- Institute of Biotechnology, University of Helsinki, 00014 Helsinki, Finland
| | - O H Samuli Ollila
- Institute of Biotechnology, University of Helsinki, 00014 Helsinki, Finland
| |
Collapse
|
8
|
Antila HS, Kav B, Miettinen MS, Martinez-Seara H, Jungwirth P, Ollila OHS. Emerging Era of Biomolecular Membrane Simulations: Automated Physically-Justified Force Field Development and Quality-Evaluated Databanks. J Phys Chem B 2022. [DOI: 10.1021/acs.jpcb.2c01954] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Hanne S. Antila
- Department of Biomaterials, Max Planck Institute of Colloids and Interfaces, 14424 Potsdam, Germany
| | - Batuhan Kav
- Institute of Biological Information Processing, Structural Biochemistry (IBI-7), Forschungszentrum
Jülich, Wilhelm-Johnen-Str., 52425 Jülich, Germany
| | - Markus S. Miettinen
- Computational Biology Unit, Department of Informatics, University of Bergen, 5008 Bergen, Norway
- Department of Chemistry, University of Bergen, 5020 Bergen, Norway
| | - Hector Martinez-Seara
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nam. 2, 16000 Prague 6, Czech Republic
| | - Pavel Jungwirth
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nam. 2, 16000 Prague 6, Czech Republic
| | - O. H. Samuli Ollila
- Institute of Biotechonology, University of Helsinki, Helsinki 00014, Finland
| |
Collapse
|
9
|
Zhou M, Yang H, Li H, Gu L, Zhou Y, Li M. The effects of molecular weight and orientation on the membrane permeation and partitioning of polycyclic aromatic hydrocarbons: a computational study. Phys Chem Chem Phys 2022; 24:2158-2166. [PMID: 35005759 DOI: 10.1039/d1cp04777a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Membrane permeation and the partitioning of polycyclic aromatic hydrocarbons (PAHs) are crucial aspects affecting their carcinogenicity and mutagenicity. However, a clear understanding of these processes is still rare due to the difficulty of determining the details experimentally. Here, the interactions between PAHs and lipid bilayers were studied by molecular simulations, mainly to check the influence of molecular weight and orientation. The liposome-water partition coefficient (KLW), transmembrane time (τ), and permeability coefficient (P) of the PAHs were calculated by integrating free energy profiles from umbrella sampling. For selected PAHs, the membrane adsorption is a spontaneous process. The preferred location is near the CC bond and the orientation is related to the molecular structure. The P values of all the PAHs are basically the same order of magnitude, which means that the molecular weight contributes little to the process. As for KLW and τ, they show obvious increases with different molecular weights. Unconstrained simulations showed that a flat orientation on the membrane surface would prevent PAHs from being transported through the membrane. Highly hydrophobic driving forces are not always good for the absorption of PAHs, especially the formation of aggregates. In addition, the orientations and energetic barriers of PAHs near the midplane of the lipid bilayer explain the different transitions of high- and low-weight PAHs. This work provides molecular level details relating to the interactions of PAHs with lipid membranes, with significance for understanding the health effects of PAHs.
Collapse
Affiliation(s)
- Mi Zhou
- School of Materials Science & Engineering, Beijing Institute of Technology, Beijing 100081, China.,Institute of Chemical Materials, Chinese Academy of Engineering and Physics, 621900 Mianyang, China.
| | - Hong Yang
- Institute of Chemical Materials, Chinese Academy of Engineering and Physics, 621900 Mianyang, China.
| | - Huarong Li
- Institute of Chemical Materials, Chinese Academy of Engineering and Physics, 621900 Mianyang, China.
| | - Lingzhi Gu
- Institute of Chemical Materials, Chinese Academy of Engineering and Physics, 621900 Mianyang, China.
| | - Yang Zhou
- Institute of Chemical Materials, Chinese Academy of Engineering and Physics, 621900 Mianyang, China.
| | - Ming Li
- Institute of Chemical Materials, Chinese Academy of Engineering and Physics, 621900 Mianyang, China.
| |
Collapse
|
10
|
Róg T, Girych M, Bunker A. Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design. Pharmaceuticals (Basel) 2021; 14:1062. [PMID: 34681286 PMCID: PMC8537670 DOI: 10.3390/ph14101062] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2021] [Revised: 10/14/2021] [Accepted: 10/15/2021] [Indexed: 11/17/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Tomasz Róg
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland;
| | - Mykhailo Girych
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland;
| | - Alex Bunker
- Drug Research Program, Division of Pharmaceutical Biosciences, Faculty of Pharmacy, University of Helsinki, 00014 Helsinki, Finland;
| |
Collapse
|
11
|
Szlenk CT, Gc JB, Natesan S. Membrane-facilitated receptor access and binding mechanisms of long-acting β2-adrenergic receptor (β2-AR) agonists. Mol Pharmacol 2021; 100:406-427. [PMID: 34334369 DOI: 10.1124/molpharm.121.000285] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2021] [Accepted: 07/20/2021] [Indexed: 11/22/2022] Open
Abstract
The drugs salmeterol, formoterol, and salbutamol constitute the frontline treatment for asthma and other chronic pulmonary diseases. These drugs activate the β2-adrenergic receptors (β2-AR), a class A G-protein-coupled receptor (GPCR) and differ significantly in their clinical onset and duration of actions. According to the "microkinetic model," the long duration of action of salmeterol and formoterol compared to salbutamol were attributed, at least in part, to their high lipophilicity and increased local concentrations in the membrane near the receptor. However, the structural and molecular bases of how the lipophilic drugs reach the binding site of the receptor from the surrounding membrane remain unknown. Using a variety of classical and enhanced molecular dynamics simulation techniques, we investigated the membrane partitioning characteristics, binding, and unbinding mechanisms of the ligands. The obtained results offer remarkable insight into the functional role of membrane lipids in the ligand association process. Strikingly, salmeterol entered the binding site from the bilayer through transmembrane helices 1 and 7. The entry was preceded by membrane-facilitated rearrangement and presentation of its phenyl-alkoxy-alkyl tail as a passkey to an access route gated by F193, a residue known critical for salmeterol's affinity. Formoterol's access is through the aqueous path shared by other β2-AR agents. We observed a novel secondary path for salbutamol that is distinct from its primary route. Our study offers a mechanistic description for the membrane-facilitated access and binding of ligands to β2-AR and establishes a groundwork for recognizing membrane lipids as an integral component in the molecular recognition process. Significance Statement The cell membrane's functional role behind the duration of action of long-acting β2-adrenergic receptor (β2-AR) agonists such as salmeterol has been a subject of debate for a long time. We investigated the binding and unbinding mechanisms of the three commonly used β2-AR agonists, salmeterol, formoterol, and salbutamol, using advanced simulation techniques. The obtained results offer unprecedented insights into the active role of membrane lipids in facilitating access and binding of the ligands, affecting the molecular recognition process and their pharmacology.
Collapse
Affiliation(s)
| | | | - Senthil Natesan
- Pharmaceutical Sciences, Washington State University, United States
| |
Collapse
|
12
|
Molecular simulations of lipid membrane partitioning and translocation by bacterial quorum sensing modulators. PLoS One 2021; 16:e0246187. [PMID: 33561158 PMCID: PMC7872223 DOI: 10.1371/journal.pone.0246187] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2020] [Accepted: 01/15/2021] [Indexed: 12/18/2022] Open
Abstract
Quorum sensing (QS) is a bacterial communication process mediated by both native and non-native small-molecule quorum sensing modulators (QSMs), many of which have been synthesized to disrupt QS pathways. While structure-activity relationships have been developed to relate QSM structure to the activation or inhibition of QS receptors, less is known about the transport mechanisms that enable QSMs to cross the lipid membrane and access intracellular receptors. In this study, we used atomistic MD simulations and an implicit solvent model, called COSMOmic, to analyze the partitioning and translocation of QSMs across lipid bilayers. We performed umbrella sampling at atomistic resolution to calculate partitioning and translocation free energies for a set of naturally occurring QSMs, then used COSMOmic to screen the water-membrane partition and translocation free energies for 50 native and non-native QSMs that target LasR, one of the LuxR family of quorum-sensing receptors. This screening procedure revealed the influence of systematic changes to head and tail group structures on membrane partitioning and translocation free energies at a significantly reduced computational cost compared to atomistic MD simulations. Comparisons with previously determined QSM activities suggest that QSMs that are least likely to partition into the bilayer are also less active. This work thus demonstrates the ability of the computational protocol to interrogate QSM-bilayer interactions which may help guide the design of new QSMs with engineered membrane interactions.
Collapse
|
13
|
Coimbra JTS, Feghali R, Ribeiro RP, Ramos MJ, Fernandes PA. The importance of intramolecular hydrogen bonds on the translocation of the small drug piracetam through a lipid bilayer. RSC Adv 2020; 11:899-908. [PMID: 35423709 PMCID: PMC8693363 DOI: 10.1039/d0ra09995c] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Accepted: 12/02/2020] [Indexed: 12/15/2022] Open
Abstract
The number of hydrogen bond donors and acceptors is a fundamental molecular descriptor to predict the oral bioavailability of small drug candidates. In fact, the most widely used oral bioavailability rules (such as the Lipinsky's rule-of-five and the Veber rules) make use of this molecular descriptor. It is generally assumed that hydrogen bond donors and acceptors impact on passive diffusion across cell membranes, a fundamental event during drug absorption and distribution. Although the relationship between the number of these motifs and the probability of having good oral bioavailability has been studied and described for more than 20 years, little attention has been given to their spatial distribution in the molecule. In this paper, we used molecular dynamics to describe the effect of intramolecular hydrogen bonding on the passive diffusion of a small drug (piracetam) through a lipid membrane. The results indicated that the formation of an intramolecular hydrogen bond decreases the barrier for translocation by ca. 4 kcal mol-1 and increases the permeability of the tested molecule, partially compensating the desolvation penalty arising from the penetration of the drug into the biological membrane core. This effect was apparent in simulations where the formation of this interaction was prevented with the help of modified potentials, and in simulations with a similar compound to piracetam that was not able to form this intramolecular hydrogen bond due to a larger distance between the hydrogen bond donor and acceptor groups. These results were also supported by coarse-grained methods, which are becoming an important resource for sampling a larger chemical space of molecules, with reduced computational effort. Furthermore, entropy and enthalpy derived profiles were also obtained as the compounds translocated across the membrane, suggesting that, even though the process of formation of internal hydrogen bonds is entropically unfavorable, the enthalpic gain is such that the formation of these interactions is beneficial for the passive diffusion across cell membranes.
Collapse
Affiliation(s)
- João T S Coimbra
- LAQV, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto Rua do Campo Alegre, s/n 4169-007 Porto Portugal
| | - Ralph Feghali
- LAQV, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto Rua do Campo Alegre, s/n 4169-007 Porto Portugal
| | - Rui P Ribeiro
- LAQV, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto Rua do Campo Alegre, s/n 4169-007 Porto Portugal
| | - Maria J Ramos
- LAQV, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto Rua do Campo Alegre, s/n 4169-007 Porto Portugal
| | - Pedro A Fernandes
- LAQV, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto Rua do Campo Alegre, s/n 4169-007 Porto Portugal
| |
Collapse
|
14
|
Novel Positive Allosteric Modulators of µ Opioid Receptor-Insight from In Silico and In Vivo Studies. Int J Mol Sci 2020; 21:ijms21228463. [PMID: 33187107 PMCID: PMC7697543 DOI: 10.3390/ijms21228463] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2020] [Revised: 11/02/2020] [Accepted: 11/06/2020] [Indexed: 02/06/2023] Open
Abstract
Opioids are the drugs of choice in severe pain management. Unfortunately, their use involves serious, potentially lethal side effects. Therefore, efforts in opioid drug design turn toward safer and more effective mechanisms, including allosteric modulation. In this study, molecular dynamics simulations in silico and ‘writhing’ tests in vivo were used to characterize potential allosteric mechanism of two previously reported compounds. The results suggest that investigated compounds bind to μ opioid receptor in an allosteric site, augmenting action of morphine at subeffective doses, and exerting antinociceptive effect alone at higher doses. Detailed analysis of in silico calculations suggests that first of the compounds behaves more like allosteric agonist, while the second compound acts mainly as a positive allosteric modulator.
Collapse
|
15
|
Grote F, Lyubartsev AP. Optimization of Slipids Force Field Parameters Describing Headgroups of Phospholipids. J Phys Chem B 2020; 124:8784-8793. [PMID: 32926619 PMCID: PMC7584372 DOI: 10.1021/acs.jpcb.0c06386] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The molecular mechanics force field Slipids developed in a series of works by Jämbeck and Lyubartsev (J. Phys. Chem. B 2012, 116, 3164-3179; J. Chem. Theory Comput. 2012, 8, 2938-2948) generally provides a good description of various lipid bilayer systems. However, it was also found that order parameters of C-H bonds in the glycerol moiety of the phosphatidylcholine headgroup deviate significantly from NMR results. In this work, the dihedral force field parameters have been reparameterized in order to improve the agreement with experiment. For this purpose, we have computed energies for a large amount of lipid headgroup conformations using density functional theory on the B3P86/cc-pvqz level and optimized dihedral angle parameters simultaneously to provide the best fit to the quantum chemical energies. The new parameter set was validated for three lipid bilayer systems against a number of experimental properties including order parameters, area per lipid, scattering form factors, bilayer thickness, area compressibility and lateral diffusion coefficients. In addition, the order parameter dependence on cholesterol content in the POPC bilayer was investigated. It is shown that the new force field significantly improves agreement with the experimental order parameters for the lipid headgroup while keeping good agreement with other experimentally measured properties.
Collapse
Affiliation(s)
- Fredrik Grote
- Department of Materials and Environmental Chemistry, Stockholm University, SE 106 91, Stockholm, Sweden
| | - Alexander P Lyubartsev
- Department of Materials and Environmental Chemistry, Stockholm University, SE 106 91, Stockholm, Sweden
| |
Collapse
|
16
|
Rogers JR, Geissler PL. Breakage of Hydrophobic Contacts Limits the Rate of Passive Lipid Exchange between Membranes. J Phys Chem B 2020; 124:5884-5898. [DOI: 10.1021/acs.jpcb.0c04139] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Julia R. Rogers
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Phillip L. Geissler
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
17
|
Loschwitz J, Olubiyi OO, Hub JS, Strodel B, Poojari CS. Computer simulations of protein-membrane systems. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:273-403. [PMID: 32145948 PMCID: PMC7109768 DOI: 10.1016/bs.pmbts.2020.01.001] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.
Collapse
Affiliation(s)
- Jennifer Loschwitz
- Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany; Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany
| | - Olujide O Olubiyi
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany; Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Obafemi Awolowo University, Ile-Ife, Nigeria
| | - Jochen S Hub
- Theoretical Physics and Center for Biophysics, Saarland University, Saarbrücken, Germany
| | - Birgit Strodel
- Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany; Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany
| | - Chetan S Poojari
- Theoretical Physics and Center for Biophysics, Saarland University, Saarbrücken, Germany.
| |
Collapse
|
18
|
Ermilova I, Lyubartsev AP. Modelling of interactions between Aβ(25-35) peptide and phospholipid bilayers: effects of cholesterol and lipid saturation. RSC Adv 2020; 10:3902-3915. [PMID: 35492630 PMCID: PMC9048594 DOI: 10.1039/c9ra06424a] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2019] [Accepted: 12/14/2019] [Indexed: 11/25/2022] Open
Abstract
Aggregation of amyloid beta (Aβ) peptides in neuronal membranes is a known promoter of Alzheimer’s disease. To gain insight into the molecular details of Aβ peptide aggregation and its effect on model neuronal membranes, we carried out molecular dynamics simulations of the Aβ(25–35) fragment of the amyloid precursor protein in phospholipid bilayers composed of either fully saturated or highly unsaturated lipids, in the presence or absence of cholesterol. It was found that the peptide does not penetrate through any of the considered membranes, but can reside in the headgroup region and upper part of the lipid tails showing a clear preference to a polyunsaturated cholesterol-free membrane. Due to the ordering and condensing effect upon addition of cholesterol, membranes become more rigid facilitating peptide aggregation on the surface. Except for the case of the cholesterol-free saturated lipid bilayer, the peptides have a small effect on the membrane structure and ordering. It was also found that the most “active” amino-acid for peptide–lipid and peptide–cholesterol interaction is methionine-35, followed by asparagine-27 and serine-26, which form hydrogen bonds between peptides and polar atoms of lipid headgroups. These amino acids are also primarily responsible for peptide aggregation. This work will be relevant for designing strategies to develop drugs to combat Alzheimer’s disease. Molecular dynamics simulations of Aβ(25–35) peptides in phospholipid bilayers are carried out to investigate the effect of polyunsaturated lipids and cholesterol on aggregation of the peptides. ![]()
Collapse
Affiliation(s)
- Inna Ermilova
- Department of Materials and Environmental Chemistry, Stockholm University Stockholm Sweden +46 8161193
| | - Alexander P Lyubartsev
- Department of Materials and Environmental Chemistry, Stockholm University Stockholm Sweden +46 8161193
| |
Collapse
|
19
|
Juračka J, Šrejber M, Melíková M, Bazgier V, Berka K. MolMeDB: Molecules on Membranes Database. DATABASE-THE JOURNAL OF BIOLOGICAL DATABASES AND CURATION 2020; 2019:5523873. [PMID: 31250015 PMCID: PMC6597476 DOI: 10.1093/database/baz078] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/16/2018] [Revised: 05/22/2019] [Accepted: 05/23/2019] [Indexed: 12/16/2022]
Abstract
Biological membranes act as barriers or reservoirs for many compounds within the human body. As such, they play an important role in pharmacokinetics and pharmacodynamics of drugs and other molecular species. Until now, most membrane/drug interactions have been inferred from simple partitioning between octanol and water phases. However, the observed variability in membrane composition and among compounds themselves stretches beyond such simplification as there are multiple drug–membrane interactions. Numerous experimental and theoretical approaches are used to determine the molecule–membrane interactions with variable accuracy, but there is no open resource for their critical comparison. For this reason, we have built Molecules on Membranes Database (MolMeDB), which gathers data about over 3600 compound–membrane interactions including partitioning, penetration and positioning. The data have been collected from scientific articles published in peer-reviewed journals and complemented by in-house calculations from high-throughput COSMOmic approach to set up a baseline for further comparison. The data in MolMeDB are fully searchable and browsable by means of name, SMILES, membrane, method or dataset and we offer the collected data openly for further reuse and we are open to further additions. MolMeDB can be a powerful tool that could help researchers better understand the role of membranes and to compare individual approaches used for the study of molecule/membrane interactions.
Collapse
Affiliation(s)
- Jakub Juračka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, Tř. 17, listopadu 12, 771 46 Olomouc, Czech Republic
| | - Martin Šrejber
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, Tř. 17, listopadu 12, 771 46 Olomouc, Czech Republic
| | - Michaela Melíková
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, Tř. 17, listopadu 12, 771 46 Olomouc, Czech Republic
| | - Václav Bazgier
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, Tř. 17, listopadu 12, 771 46 Olomouc, Czech Republic
| | - Karel Berka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, Tř. 17, listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
20
|
Paloncýová M, Aniander G, Larsson E, Knippenberg S. Cyanine dyes with tail length asymmetry enhance photoselection: A multiscale study on DiD probes in a liquid disordered membrane. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2020; 224:117329. [PMID: 31326855 DOI: 10.1016/j.saa.2019.117329] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/23/2019] [Revised: 06/25/2019] [Accepted: 06/28/2019] [Indexed: 06/10/2023]
Abstract
Visualization of membrane domains like lipid rafts in natural or artificial membranes is a crucial task for cell biology. For this purpose, fluorescence microscopy is often used. Since fluorescing probes in lipid membranes partition specifically in e.g. local liquid disordered or liquid ordered environments, the consequent changes in their orientation and location are both theoretically and experimentally of interest. Here we focused on a liquid disordered membrane phase and performed molecular dynamics (MD) simulations of the indocarbocyanine DiD probes by varying the length of the attached alkyl tails and also the length of the cyanine backbone. From the probed compounds in a DOPC lipid bilayer at ambient temperature, a varying orientation of the transition dipole moment was observed, which is crucial for fluorescence microscopy and which, through photoselection, was found to be surprisingly more effective for asymmetric probes than for the symmetric ones. Furthermore, we observed that the orientation of the probes was dependent on the tail length; with the methyls or propyls attached, DiD oriented with its tails facing the water, contrary to the ones with longer tails. With advanced hybrid QM/MM calculations we show that the different local environment for differently oriented probes affected the one-photon absorption spectra, that was blue-shifted for the short-tailed DiD with respect to the DiDs with longer tails. We show here that the presented probes can be successfully used for fluorescence microscopy and we believe that the described properties bring further insight for the experimental use of these probes.
Collapse
Affiliation(s)
- Markéta Paloncýová
- Department of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, Stockholm, Sweden; Biomedical Research Institute, Hasselt University, Agoralaan Building C, 3590 Diepenbeek, Belgium
| | - Gustav Aniander
- Department of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, Stockholm, Sweden
| | - Emma Larsson
- Department of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, Stockholm, Sweden
| | - Stefan Knippenberg
- Department of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, Stockholm, Sweden; Biomedical Research Institute, Hasselt University, Agoralaan Building C, 3590 Diepenbeek, Belgium; RCPTM, Department of Physical Chemistry, Fac. Sciences, Palacký University, Olomouc, Czech Republic; Research group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, 2610 Wilrijk, Antwerp, Belgium.
| |
Collapse
|
21
|
Ntarakas N, Ermilova I, Lyubartsev AP. Effect of lipid saturation on amyloid-beta peptide partitioning and aggregation in neuronal membranes: molecular dynamics simulations. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2019; 48:813-824. [PMID: 31655893 PMCID: PMC6853862 DOI: 10.1007/s00249-019-01407-x] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/27/2019] [Revised: 10/02/2019] [Accepted: 10/13/2019] [Indexed: 12/05/2022]
Abstract
Aggregation of amyloid-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta $$\end{document}β (Aβ) peptides, cleaved from the amyloid precursor protein, is known as a precursor of the Alzheimer’s disease (AD). It is also known that Alzheimer’s disease is characterized by a substantial decrease of the amount of polyunsaturated lipids in the neuronal membranes of the frontal gray matter. To get insight into possible interconnection of these phenomena, we have carried out molecular dynamics simulations of two fragments of A\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta $$\end{document}β peptide, A\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta $$\end{document}β\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$_{1-28}$$\end{document}1-28 and A\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta $$\end{document}β\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$_{26-40}$$\end{document}26-40, in four different lipid bilayers: two monocomponent ones (14:0-14:0 PC, 18:0-22:6 PC), and two bilayers containing mixtures of 18:0-18:0 PE, 22:6-22:6 PE, 16:0-16:0 PC and 18:1-18:1 PC lipids of composition mimicking neuronal membranes in a “healthy” and “AD” brain. The simulations showed that the presence of lipids with highly unsaturated 22:6cis fatty acids chains strongly affects the interaction of amyloid-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta $$\end{document}β peptides with lipid membranes. The polyunsaturated lipids cause stronger adsorption of A\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta $$\end{document}β-peptides by the membrane and lead to weaker binding between peptides when the latter form aggregates. This difference in the behaviour observed in monocomponent bilayers is propagated in a similar fashion to the mixed membranes mimicking composition of neuronal membranes in “healthy” and “AD” brains, with “healthy” membrane having higher fraction of polyunsaturated lipids. Our simulations give strong indication that it can be physical–chemical background of the interconnection between amyloid fibrillization causing Alzheimer’s disease, and content of polyunsaturated lipids in the neuronal membranes.
Collapse
Affiliation(s)
- Nikolaos Ntarakas
- Department of Materials and Environmental Chemistry, Stockholm's University, 10691, Stockholm, Sweden
| | - Inna Ermilova
- Department of Materials and Environmental Chemistry, Stockholm's University, 10691, Stockholm, Sweden
| | - Alexander P Lyubartsev
- Department of Materials and Environmental Chemistry, Stockholm's University, 10691, Stockholm, Sweden.
| |
Collapse
|
22
|
Osella S, Smisdom N, Ameloot M, Knippenberg S. Conformational Changes as Driving Force for Phase Recognition: The Case of Laurdan. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:11471-11481. [PMID: 31403301 DOI: 10.1021/acs.langmuir.9b01840] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The development of a universal probe to assess the phase of a lipid membrane is one of the most ambitious goals for fluorescence spectroscopy. The ability of a well-known molecule as Laurdan to reach this aim is here exploited as the behavior of the probe is fully characterized in a dipalmitoylphosphatidylcholine (DPPC) solid gel (So) phase by means of molecular dynamics simulations. Laurdan can take two conformations, depending on whether the carbonyl oxygen points toward the β-position of the naphthalene core (Conf-I) or to the α-position (Conf-II). We observe that Conf-I has an elongated form in this environment, whereas Conf-II takes an L-shape. Interestingly, our theoretical calculations show that these two conformations behave in an opposite way from what is reported in the literature for a DPPC membrane in a liquid disordered (Ld) phase, where Conf-I assumes an L-shape and Conf-II is elongated. Moreover, our results show that in DPPC (So) no intermixing between the conformations is present, whereas it has been seen in a fluid environment such as DOPC (Ld). Through a careful analysis of angle distributions and by means of the rotational autocorrelation function, we predict that the two conformers of Laurdan behave differently in different membrane environments.
Collapse
Affiliation(s)
- Silvio Osella
- Chemical and Biological Systems Simulation Lab, Centre of New Technologies , University of Warsaw , Banacha 2C , 02-097 Warsaw , Poland
| | - Nick Smisdom
- Biomedical Research Institute , Hasselt University , Agoralaan Building C , 3590 Diepenbeek , Belgium
| | - Marcel Ameloot
- Biomedical Research Institute , Hasselt University , Agoralaan Building C , 3590 Diepenbeek , Belgium
| | - Stefan Knippenberg
- Biomedical Research Institute , Hasselt University , Agoralaan Building C , 3590 Diepenbeek , Belgium
- Department of Theoretical Chemistry and Biology, School of Biotechnology , Royal Institute of Technology , SE-10691 Stockholm , Sweden
- RCPTM, Department of Physical Chemistry, Fac. Sciences , Palacký University , 771 46 Olomouc , Czech Republic
| |
Collapse
|
23
|
Yang H, Li H, Liu L, Zhou Y, Long X. Molecular Simulation Studies on the Interactions of 2,4,6-Trinitrotoluene and Its Metabolites with Lipid Membranes. J Phys Chem B 2019; 123:6481-6491. [DOI: 10.1021/acs.jpcb.9b03033] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Hong Yang
- School of Material Science and Engineering, Tsinghua University, Beijing 100084, China
- Institute of Chemical Materials, China Academy of Engineering and Physics, Mianyang 621900, China
| | - Huarong Li
- Institute of Chemical Materials, China Academy of Engineering and Physics, Mianyang 621900, China
| | - Liu Liu
- Institute of Chemical Materials, China Academy of Engineering and Physics, Mianyang 621900, China
| | - Yang Zhou
- Institute of Chemical Materials, China Academy of Engineering and Physics, Mianyang 621900, China
| | - Xinping Long
- Institute of Chemical Materials, China Academy of Engineering and Physics, Mianyang 621900, China
| |
Collapse
|
24
|
Subramaniam R, Lynch S, Cen Y, Balaz S. Polarity of Hydrated Phosphatidylcholine Headgroups. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:8460-8471. [PMID: 31244216 PMCID: PMC6853183 DOI: 10.1021/acs.langmuir.8b03992] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The headgroup (H) stratum (sometimes called the polar region) of membrane bilayers is a relevant yet poorly understood solvation phase for small molecules and macromolecules interacting with the membranes. Solvation of compounds in bilayer strata is characterized experimentally by wide- and small-angle X-ray scattering, neutron diffraction, and various NMR techniques. The quantification is tedious and only available for a limited set of small molecules. Our recently published model of liposome partitioning of small molecules shows that solvation of compounds in the H-stratum of fluid phosphatidylcholine (PC) bilayers correlates well with their solvation in hydrated diacetyl phosphatidylcholine (DAcPC), and solvation in the core (C) depends in a similar way on that in n-hexadecane. These two correlations became a basis for a model describing the location of compounds in the H- and C-strata and at the connecting interface as a nonlinear function of the fragment solvation characteristics of the compounds. In this study, refractivity of hydrated DAcPC phases with varying water contents was measured and polarity was determined using the steady-state fluorescence of indole and Nile Red. The results were compared with the published data obtained by other techniques for PC bilayers in liposomes or on solid supports. The demonstrated qualitative agreement, as well as the polarity and refractivity dependencies on the DAcPC concentration, supports the suitability of hydrated DAcPC as the H-stratum surrogate. Interestingly, depending on hydrations typical for the H-strata of fluid PC bilayers, the dielectric constant could decrease significantly from 31.0 to 7.3 for 16 and 8 water molecules per headgroup, respectively. Although additional experiments are needed for confirmation, this observation could help set proper dielectric constant magnitudes in continuum-based computational models of accumulation and crossing of the PC bilayers with varying hydration levels thanks to the temperature or the structure of fatty acid chains.
Collapse
Affiliation(s)
| | | | | | - Stefan Balaz
- Corresponding author: Stefan Balaz, Albany College of Pharmacy and Health Sciences, Vermont Campus, Department of Pharmaceutical Sciences, 261 Mountain View Road, Colchester, VT 05446, United States, phone 802-735-2615,
| |
Collapse
|
25
|
Hexagonal boron nitride nanosheet as novel drug delivery system for anticancer drugs: Insights from DFT calculations and molecular dynamics simulations. J Mol Graph Model 2019; 89:50-59. [DOI: 10.1016/j.jmgm.2019.02.012] [Citation(s) in RCA: 74] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2018] [Revised: 12/28/2018] [Accepted: 02/25/2019] [Indexed: 11/19/2022]
|
26
|
Rivel T, Ramseyer C, Yesylevskyy S. The asymmetry of plasma membranes and their cholesterol content influence the uptake of cisplatin. Sci Rep 2019; 9:5627. [PMID: 30948733 PMCID: PMC6449338 DOI: 10.1038/s41598-019-41903-w] [Citation(s) in RCA: 81] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2018] [Accepted: 03/18/2019] [Indexed: 12/20/2022] Open
Abstract
The composition of the plasma membrane of malignant cells is thought to influence the cellular uptake of cisplatin and to take part in developing resistance to this widespread anti-cancer drug. In this work we study the permeation of cisplatin through the model membranes of normal and cancer cells using molecular dynamics simulations. A special attention is paid to lipid asymmetry and cholesterol content of the membranes. The loss of lipid asymmetry, which is common for cancer cells, leads to a decrease in their permeability to cisplatin by one order of magnitude in comparison to the membranes of normal cells. The change in the cholesterol molar ratio from 0% to 33% also decreases the permeability of the membrane by approximately one order of magnitude. The permeability of pure DOPC membrane is 5-6 orders of magnitude higher than one of the membranes with realistic lipid composition, which makes it as an inadequate model for the studies of drug permeability.
Collapse
Affiliation(s)
- Timothée Rivel
- Laboratoire Chrono Environnement UMR CNRS 6249, Université de Bourgogne Franche-Comté, 16 route de Gray, 25030, Besançon, Cedex, France.
| | - Christophe Ramseyer
- Laboratoire Chrono Environnement UMR CNRS 6249, Université de Bourgogne Franche-Comté, 16 route de Gray, 25030, Besançon, Cedex, France
| | - Semen Yesylevskyy
- Laboratoire Chrono Environnement UMR CNRS 6249, Université de Bourgogne Franche-Comté, 16 route de Gray, 25030, Besançon, Cedex, France.,Department of Physics of Biological Systems, Institute of Physics of the National Academy of Sciences of Ukraine, Prospect Nauky 46, 03028, Kyiv, Ukraine
| |
Collapse
|
27
|
Paloncýová M, Ameloot M, Knippenberg S. Orientational distribution of DPH in lipid membranes: a comparison of molecular dynamics calculations and experimental time-resolved anisotropy experiments. Phys Chem Chem Phys 2019; 21:7594-7604. [DOI: 10.1039/c8cp07754a] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
The behavior of the fluorescent probe diphenylhexatriene (DPH) in different lipid phases is investigated. The rotational autocorrelation functions are calculated in order to model the time-resolved fluorescence anisotropy decay. The role of the order parameters is discussed.
Collapse
Affiliation(s)
- Markéta Paloncýová
- Department of Theoretical Chemistry and Biology
- School of Engineering Sciences in Chemistry
- Biotechnology and Health
- Royal Institute of Technology
- Stockholm
| | - Marcel Ameloot
- Biomedical Research Institute, Hasselt University
- 3590 Diepenbeek
- Belgium
| | - Stefan Knippenberg
- Department of Theoretical Chemistry and Biology
- School of Engineering Sciences in Chemistry
- Biotechnology and Health
- Royal Institute of Technology
- Stockholm
| |
Collapse
|
28
|
Osella S, Di Meo F, Murugan NA, Fabre G, Ameloot M, Trouillas P, Knippenberg S. Combining (Non)linear Optical and Fluorescence Analysis of DiD To Enhance Lipid Phase Recognition. J Chem Theory Comput 2018; 14:5350-5359. [DOI: 10.1021/acs.jctc.8b00553] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Silvio Osella
- Centre of New Technologies, University of Warsaw, Banacha 2C, 02-097 Warsaw, Poland
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, Royal Institute of Technology, SE-10691 Stockholm, Sweden
| | - Florent Di Meo
- Faculty of Pharmacy, UMR 1248 INSERM, Limoges University, 2 rue du Docteur Marcland, 87025 Limoges Cedex, France
| | - N. Arul Murugan
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, Royal Institute of Technology, SE-10691 Stockholm, Sweden
| | - Gabin Fabre
- LCSN-EA1069, Faculty of Pharmacy, Limoges University, 2, rue du Dr. Marcland, 87025 Limoges Cedex, France
| | - Marcel Ameloot
- Biomedical Research Institute, Hasselt University, B-3590, Diepenbeek, Belgium
| | - Patrick Trouillas
- Faculty of Pharmacy, UMR 1248 INSERM, Limoges University, 2 rue du Docteur Marcland, 87025 Limoges Cedex, France
- Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Stefan Knippenberg
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, Royal Institute of Technology, SE-10691 Stockholm, Sweden
- Biomedical Research Institute, Hasselt University, B-3590, Diepenbeek, Belgium
| |
Collapse
|
29
|
Knippenberg S, Fabre G, Osella S, Di Meo F, Paloncýová M, Ameloot M, Trouillas P. Atomistic Picture of Fluorescent Probes with Hydrocarbon Tails in Lipid Bilayer Membranes: An Investigation of Selective Affinities and Fluorescent Anisotropies in Different Environmental Phases. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2018; 34:9072-9084. [PMID: 29983063 DOI: 10.1021/acs.langmuir.8b01164] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
By reverting to spectroscopy, changes in the biological environment of a fluorescent probe can be monitored and the presence of various phases of the surrounding lipid bilayer membranes can be detected. However, it is currently not always clear in which phase the probe resides. The well-known orange 1,1'-dioctadecyl-3,3,3',3'-tetramethylindodicarbo-cyanine perchlorate (DiI-C18(5)) fluorophore, for instance, and the new, blue BODIPY (4,4-difluoro-4-bora-3 a,4 a-diaza- s-indacene) derivative were experimentally seen to target and highlight identical parts of giant unilamellar vesicles of various compositions, comprising mixtures of dipalmitoylphosphatidylcholine (DPPC), dioleoylphosphatidylcholine (DOPC), sphingomyelin (SM), and cholesterol (Chol). However, it was not clear which of the coexisting membrane phases were visualized (Bacalum et al., Langmuir. 2016, 32, 3495). The present study addresses this issue by utilizing large-scale molecular dynamics simulations and the z-constraint method, which allows evaluating Gibbs free-energy profiles. The current calculations give an indication why, at room temperature, both BODIPY and DiI-C18(5) probes prefer the gel (So) phase in DOPC/DPPC (2:3 molar ratio) and the liquid-ordered (Lo) phase in DOPC/SM/Chol (1:2:1 molar ratio) mixtures. This study highlights the important differences in orientation and location and therefore in efficiency between the probes when they are used in fluorescence microscopy to screen various lipid bilayer membrane phases. Dependent on the lipid composition, the angle between the transition-state dipole moments of both probes and the normal to the membrane is found to deviate clearly from 90°. It is seen that the DiI-C18(5) probe is located in the headgroup region of the SM/Chol mixture, in close contact with water molecules. A fluorescence anisotropy study also indicates that DiI-C18(5) gives rise to a distinctive behavior in the SM/Chol membrane compared to the other considered membranes. The latter behavior has not been seen for the studied BODIPY probe, which is located deeper in the membrane.
Collapse
Affiliation(s)
- S Knippenberg
- Department of Theoretical Chemistry and Biology , KTH Royal Institute of Technology , Roslagstullsbacken 15 , S-106 91 Stockholm , Sweden
- Biomedical Research Institute , Hasselt University , Agoralaan Building C , 3590 Diepenbeek , Belgium
| | - G Fabre
- LCSN-EA1069, Faculty of Pharmacy , Limoges University , 2 rue du Dr. Marcland , 87025 Limoges Cedex , France
| | - S Osella
- Centre of New Technologies , University of Warsaw , Banacha 2C , 02-097 Warsaw , Poland
| | - F Di Meo
- Faculty of Pharmacy , INSERM UMR 1248, Limoges University , 2 rue du Docteur Marcland , 87025 Limoges Cedex , France
| | - M Paloncýová
- Department of Theoretical Chemistry and Biology , KTH Royal Institute of Technology , Roslagstullsbacken 15 , S-106 91 Stockholm , Sweden
| | - M Ameloot
- Biomedical Research Institute , Hasselt University , Agoralaan Building C , 3590 Diepenbeek , Belgium
| | - P Trouillas
- Faculty of Pharmacy , INSERM UMR 1248, Limoges University , 2 rue du Docteur Marcland , 87025 Limoges Cedex , France
- Centre of Advanced Technologies and Materials, Faculty of Science , Palacký University , tř. 17 listopadu 12 , 771 46 Olomouc , Czech Republic
| |
Collapse
|
30
|
Location of fluorescent probes (2′-hydroxy derivatives of 2,5-diaryl-1,3-oxazole) in lipid membrane studied by fluorescence spectroscopy and molecular dynamics simulation. Biophys Chem 2018; 235:9-18. [DOI: 10.1016/j.bpc.2018.01.005] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2018] [Revised: 01/26/2018] [Accepted: 01/28/2018] [Indexed: 11/20/2022]
|
31
|
Lu Y, Shi XF, Salsbury FR, Derreumaux P. Influence of electric field on the amyloid-β(29-42) peptides embedded in a membrane bilayer. J Chem Phys 2018; 148:045105. [DOI: 10.1063/1.5018459] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yan Lu
- School of Physics and Optoelectronic Engineering, Xidian University, Xi’an 710071, China
| | - Xiao-Feng Shi
- School of Physics and Optoelectronic Engineering, Xidian University, Xi’an 710071, China
| | - Freddie R. Salsbury
- Department of Physics, Wake Forest University, Winston-Salem, North Carolina 27106, USA
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, Institut de Biologie Physico-Chimique (IBPC), UPR9080 CNRS, Université Paris Diderot, Sorbonne Paris Cité, 13 Rue Pierre et Marie Curie, 75005 Paris, France
| |
Collapse
|
32
|
Coimbra JTS, Brás NF, Fernandes PA, Rangel M, Ramos MJ. Membrane partition of bis-(3-hydroxy-4-pyridinonato) zinc(ii) complexes revealed by molecular dynamics simulations. RSC Adv 2018; 8:27081-27090. [PMID: 35539964 PMCID: PMC9083369 DOI: 10.1039/c8ra03602k] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2018] [Accepted: 07/23/2018] [Indexed: 12/24/2022] Open
Abstract
The class of 3-hydroxy-4-pyridinone ligands is widely known and valuable for biomedical and pharmaceutical purposes. Their chelating properties towards biologically-relevant transition metal ions highlight their potential biomedical utility. A set of 3-hydroxy-4-pyridinone Zn(ii) complexes at different concentrations was studied for their ability to interact with lipid phases. We employed umbrella sampling simulations to attain the potential-of-mean force for a set of ligands and one Zn(ii) complex, as these permeated a 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) hydrated bilayer system. In addition, we used conventional molecular dynamics simulations to study the behavior of various Zn(ii) complexes in hydrated bilayer systems. This work discusses: (i) the partition of 3-hydroxy-4-pyridinone ligands to bilayer phases; (ii) self-aggregation in crowded environments of Zn(ii) complexes; and (iii) possible mechanisms for the membrane translocation of Zn(ii) complexes. We observed distinct interactions for the studied complexes, and distinct membrane partition coefficients (Kmem) depending on the considered ligand. The more hydrophobic ligand, 1-hexyl-3-hydroxy-2-methyl-4(1H)-pyridinone, partitioned more favorably to lipid phases (at least two orders of magnitude higher Kmem when compared to the other ligands), and the corresponding Zn(ii) complex was also prone to self-aggregation when an increased concentration of the complex was employed. We also observed that the inclusion of a coordinated water molecule in the parameterization of the Zn(ii) coordination sphere, as proposed in the available crystallographic structure of the complex, decreased the partition coefficient and membrane permeability for the tested complex. The membrane partition of hydroxypyridinones and of zinc complexes explored by molecular dynamics.![]()
Collapse
Affiliation(s)
- João T. S. Coimbra
- UCIBIO
- REQUIMTE
- Departamento de Química e Bioquímica
- Faculdade de Ciências
- Universidade do Porto
| | - Natércia F. Brás
- UCIBIO
- REQUIMTE
- Departamento de Química e Bioquímica
- Faculdade de Ciências
- Universidade do Porto
| | - Pedro A. Fernandes
- UCIBIO
- REQUIMTE
- Departamento de Química e Bioquímica
- Faculdade de Ciências
- Universidade do Porto
| | - Maria Rangel
- LAQV
- REQUIMTE
- Instituto de Ciências Biomédicas de Abel Salazar
- Universidade do Porto
- 4050-313 Porto
| | - Maria J. Ramos
- UCIBIO
- REQUIMTE
- Departamento de Química e Bioquímica
- Faculdade de Ciências
- Universidade do Porto
| |
Collapse
|
33
|
Abstract
From the pharmacological point of view, allosteric modulators may present numerous advantages over orthosteric ligands. Growing availability of novel tools and experimental data provides a tempting opportunity to apply computational methods to improve known modulators and design novel ones. However, recent progress in understanding of complexity of allostery increases awareness of problems involved in design of modulators with desired properties. Deeper insight into phenomena such as probe dependence, altering signaling bias with minor changes in ligand structure, as well as influence of subtle endogenous allosteric factors turns out to be fundamental. These effects make the design of a modulator with precise pharmacological outcome a very challenging task, and need to be taken into consideration throughout the design process. In this chapter, we focus on nuances of targeting GPCR allosteric sites in computational drug design efforts, in particular with application of docking, virtual screening, and molecular dynamics.
Collapse
|
34
|
Bartuzi D, Kaczor AA, Matosiuk D. Molecular mechanisms of allosteric probe dependence in μ opioid receptor. J Biomol Struct Dyn 2017; 37:36-47. [PMID: 29241414 DOI: 10.1080/07391102.2017.1417914] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
Allostery is one of the most important features of proteins. It greatly contributes to the complexity of life, since it enables possibility of precise tuning of protein function, as well as performing more than one function per protein. Probe dependence is one of the unique features of allostery. It allows a protein to respond differently to the same allosteric modulator when different drugs or transmitters are bound. Unfortunately, allosteric mechanisms are difficult to investigate experimentally. Instead, they can be reproduced artificially in simulations. We simulated in silico a native-like cell membrane fragment with an active-state human μ opioid receptor (MOR) in order to investigate diverse effects of a receptor's positive allosteric modulator on various agonists. Particular emphasis on native-likeness of the environment was put. We managed to reproduce the experimentally observed effects, which allowed us to take deeper insight into their underlying mechanisms. We found an allosteric pathway in the receptor, leading from the ligand binding site to the intracellular, effector site. We observed that the modulator affected the pathway, inducing different resultant responses for full and partial agonists.
Collapse
Affiliation(s)
- Damian Bartuzi
- a Department of Synthesis and Chemical Technology of Pharmaceutical Substances with Computer Modelling Lab , Medical University of Lublin , Lublin , Poland
| | - Agnieszka A Kaczor
- a Department of Synthesis and Chemical Technology of Pharmaceutical Substances with Computer Modelling Lab , Medical University of Lublin , Lublin , Poland.,b School of Pharmacy , University of Eastern Finland , Kuopio , Finland
| | - Dariusz Matosiuk
- a Department of Synthesis and Chemical Technology of Pharmaceutical Substances with Computer Modelling Lab , Medical University of Lublin , Lublin , Poland
| |
Collapse
|
35
|
Navrátilová V, Paloncýová M, Berka K, Mise S, Haga Y, Matsumura C, Sakaki T, Inui H, Otyepka M. Molecular insights into the role of a distal F240A mutation that alters CYP1A1 activity towards persistent organic pollutants. Biochim Biophys Acta Gen Subj 2017; 1861:2852-2860. [DOI: 10.1016/j.bbagen.2017.08.002] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2017] [Revised: 06/30/2017] [Accepted: 08/01/2017] [Indexed: 01/12/2023]
|
36
|
Ermilova I, Stenberg S, Lyubartsev AP. Quantum chemical and molecular dynamics modelling of hydroxylated polybrominated diphenyl ethers. Phys Chem Chem Phys 2017; 19:28263-28274. [PMID: 29028067 DOI: 10.1039/c7cp03471g] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A series of 19 hydroxylated polybrominated diphenyl ethers (OH-PBDEs) have been studied using density functional theory (DFT) and molecular dynamics simulations with the purpose of investigating eventual correlations between their physicochemical properties and toxic action. Dissociation constants (pKa), solvation free energies and octanol-water partition coefficients (log P) have been computed. Additionally, metadynamics simulations of OH-PBDEs passing through a lipid bilayer have been carried out for four OH-PBDE species. No correlations between computed pKa values and toxicity data have been found. Medium correlations were found between partition coefficients and the ability of OH-PBDEs to alter membrane potential in cell cultures, which is attributed to higher uptake of molecules with larger log P parameters. It was also demonstrated that in lipid bilayers, OH-PBDE molecules differ in their orientational distributions and can adopt different conformations which can affect the uptake of these molecules and influence the pathways of their toxic action.
Collapse
Affiliation(s)
- Inna Ermilova
- Department of Materials and Environmental Chemistry, Stockholm University, SE 106 91, Stockholm, Sweden.
| | | | | |
Collapse
|
37
|
Predicting Critical Micelle Concentrations with Molecular Dynamics Simulations and COSMOmic. CHEM-ING-TECH 2017. [DOI: 10.1002/cite.201700061] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
38
|
Van Oosten B, Marquardt D, Harroun TA. Testing High Concentrations of Membrane Active Antibiotic Chlorhexidine Via Computational Titration and Calorimetry. J Phys Chem B 2017; 121:4657-4668. [DOI: 10.1021/acs.jpcb.6b12510] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Brad Van Oosten
- Department
of Physics, Brock University, St. Catharines, Ontario L2S 3A1, Canada
| | - Drew Marquardt
- Institute
of Molecular Biosciences, University of Graz, 8010 Graz, Austria
| | - Thad A. Harroun
- Department
of Physics, Brock University, St. Catharines, Ontario L2S 3A1, Canada
| |
Collapse
|
39
|
Coimbra JTS, Fernandes PA, Ramos MJ. Revisiting Partition in Hydrated Bilayer Systems. J Chem Theory Comput 2017; 13:2290-2299. [PMID: 28388088 DOI: 10.1021/acs.jctc.6b01258] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
We discuss potential-of-mean force convergence issues involving the methods for generating starting structures of umbrella sampling simulations for the study of the water-bilayer permeation of drug-like compounds. We provide a proof of concept of the benefits of potential-energy modification methods, in particular the flooding technique, to enhance the sampling and improve the convergence of the potential-of-mean force. In addition, we also discuss how these potential modifiers could introduce the necessary bias to induce the correct description of other relevant internal degrees of freedom for the compounds, thus providing a more accurate description of the permeation process. Effective methodological approximations are taken into consideration, more specifically the inclusion of atomic polarization, and experimental data are used in the validation of computed values whenever possible. The parametrization of the atomic point charges is always a relevant issue in biomolecular simulations. We discuss the influence of a molecule's conformation-dependent atomic point charges and their impact on the predicted potential-of-mean force.
Collapse
Affiliation(s)
- João T S Coimbra
- UCIBIO, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto , Rua do Campo Alegre, s/n, 4169-007 Porto, Portugal
| | - Pedro A Fernandes
- UCIBIO, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto , Rua do Campo Alegre, s/n, 4169-007 Porto, Portugal
| | - Maria J Ramos
- UCIBIO, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto , Rua do Campo Alegre, s/n, 4169-007 Porto, Portugal
| |
Collapse
|
40
|
Diffusion of the small, very polar, drug piracetam through a lipid bilayer: an MD simulation study. Theor Chem Acc 2017. [DOI: 10.1007/s00214-017-2073-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
41
|
Fosso-Tande J, Black C, G. Aller S, Lu L, D. Hills Jr R. Simulation of lipid-protein interactions with the CgProt force field. AIMS MOLECULAR SCIENCE 2017. [DOI: 10.3934/molsci.2017.3.352] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
|
42
|
Lopes D, Jakobtorweihen S, Nunes C, Sarmento B, Reis S. Shedding light on the puzzle of drug-membrane interactions: Experimental techniques and molecular dynamics simulations. Prog Lipid Res 2017; 65:24-44. [DOI: 10.1016/j.plipres.2016.12.001] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2016] [Revised: 11/30/2016] [Accepted: 12/03/2016] [Indexed: 12/20/2022]
|
43
|
Dickson CJ, Hornak V, Pearlstein RA, Duca JS. Structure–Kinetic Relationships of Passive Membrane Permeation from Multiscale Modeling. J Am Chem Soc 2016; 139:442-452. [DOI: 10.1021/jacs.6b11215] [Citation(s) in RCA: 80] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Callum J. Dickson
- Computer-Aided Drug Discovery,
Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Viktor Hornak
- Computer-Aided Drug Discovery,
Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Robert A. Pearlstein
- Computer-Aided Drug Discovery,
Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Jose S. Duca
- Computer-Aided Drug Discovery,
Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
44
|
Ermilova I, Lyubartsev AP. Extension of the Slipids Force Field to Polyunsaturated Lipids. J Phys Chem B 2016; 120:12826-12842. [DOI: 10.1021/acs.jpcb.6b05422] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Inna Ermilova
- Department of Materials and
Environmental Chemistry, Stockholm University, SE 106 91 Stockholm, Sweden
| | - Alexander P. Lyubartsev
- Department of Materials and
Environmental Chemistry, Stockholm University, SE 106 91 Stockholm, Sweden
| |
Collapse
|
45
|
Holmboe M, Larsson P, Anwar J, Bergström CAS. Partitioning into Colloidal Structures of Fasted State Intestinal Fluid Studied by Molecular Dynamics Simulations. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2016; 32:12732-12740. [PMID: 27934534 PMCID: PMC5271571 DOI: 10.1021/acs.langmuir.6b03008] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/11/2016] [Revised: 11/04/2016] [Indexed: 06/06/2023]
Abstract
We performed molecular dynamics (MD) simulations to obtain insights into the structure and molecular interactions of colloidal structures present in fasted state intestinal fluid. Drug partitioning and interaction were studied with a mixed system of the bile salt taurocholate (TCH) and 1,2-dilinoleoyl-sn-glycero-3-phosphocholine (DLiPC). Spontaneous aggregation of TCH and DLiPC from unconstrained MD simulations at the united-atom level using the Berger/Gromos54A7 force fields demonstrated that intermolecular hydrogen bonding between TCH molecules was an important factor in determining the overall TCH and DLiPC configuration. In bilayered systems, these intermolecular hydrogen bonds resulted in embedded transmembrane TCH clusters. Free energy simulations using the umbrella sampling technique revealed that the stability of these transmembrane TCH clusters was superior when they consisted of 3 or 4 TCH per bilayer leaflet. All-atom simulations using the Slipids/GAFF force fields showed that the TCH embedded in the bilayer decreased the energy barrier to penetrate the bilayer (ΔGpen) for water, ethanol, and carbamazepine, but not for the more lipophilic felodipine and danazol. This suggests that diffusion of hydrophilic to moderately lipophilic molecules through the bilayer is facilitated by the embedded TCH molecules. However, the effect of embedded TCH on the overall lipid/water partitioning was significant for danazol, indicating that the incorporation of TCH plays a crucial role for the partitioning of lipophilic solutes into e.g. lipidic vesicles existing in fasted state intestinal fluids. To conclude, the MD simulations revealed important intermolecular interactions in lipidic bilayers, both between the bile components themselves and with the drug molecules.
Collapse
Affiliation(s)
- Michael Holmboe
- Department
of Pharmacy, Uppsala University, Uppsala
Biomedical Centre, P.O. Box 580, SE-751 23 Uppsala, Sweden
- Department
of Chemistry, Umeå University, Umeå, Sweden
| | - Per Larsson
- Department
of Pharmacy, Uppsala University, Uppsala
Biomedical Centre, P.O. Box 580, SE-751 23 Uppsala, Sweden
| | - Jamshed Anwar
- Chemical
Theory & Computation, Department of Chemistry, University of Lancaster, Lancaster LA1 4YB, U.K.
| | - Christel A. S. Bergström
- Department
of Pharmacy, Uppsala University, Uppsala
Biomedical Centre, P.O. Box 580, SE-751 23 Uppsala, Sweden
| |
Collapse
|
46
|
McKiernan KA, Wang LP, Pande VS. Training and Validation of a Liquid-Crystalline Phospholipid Bilayer Force Field. J Chem Theory Comput 2016; 12:5960-5967. [DOI: 10.1021/acs.jctc.6b00801] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Lee-Ping Wang
- Department
of Chemistry, University of California—Davis, Davis, California 95618, United States
| | | |
Collapse
|
47
|
Pearlstein RA, Dickson CJ, Hornak V. Contributions of the membrane dipole potential to the function of voltage-gated cation channels and modulation by small molecule potentiators. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1859:177-194. [PMID: 27836643 DOI: 10.1016/j.bbamem.2016.11.005] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/26/2016] [Revised: 10/28/2016] [Accepted: 11/06/2016] [Indexed: 01/27/2023]
Abstract
The membrane dipole potential (Ψd) constitutes one of three electrical potentials generated by cell membranes. Ψd arises from the unfavorable parallel alignment of phospholipid and water dipoles, and varies in magnitude both longitudinally and laterally across the bilayer according to membrane composition and phospholipid packing density. In this work, we propose that dynamic counter-balancing between Ψd and the transmembrane potential (ΔΨm) governs the conformational state transitions of voltage-gated ion channels. Ψd consists of 1) static outer, and dynamic inner leaflet components (Ψd(extra) and Ψd(intra), respectively); and 2) a transmembrane component (ΔΨd(inner-outer)), ariing from differences in intra- and extracellular leaflet composition. Ψd(intra), which transitions between high and low energy states (Ψd(intra, high) and Ψd(intra, low)) as a function of channel conformation, is transduced by the pore domain. ΔΨd(inner-outer) is transduced by the voltage-sensing (VS) domain in summation with ΔΨm. Potentiation of voltage-gated ion channels is of interest for the treatment of cardiac, neuronal, and other disorders arising from inherited/acquired ion channel dysfunction. Potentiators are widely believed to alter the rates and voltage-dependencies of channel gating transitions by binding to pockets in the membrane-facing and other regions of ion channel targets. Here, we propose that potentiators alter Ψd(intra) and/or Ψd(extra), thereby increasing or decreasing the energy barriers governing channel gating transitions. We used quantum mechanical and molecular dynamics (MD) simulations to predict the overall Ψd-modulating effects of a series of published positive hERG potentiators partitioned into model DOPC bilayers. Our findings suggest a strong correlation between the magnitude of Ψd-lowering and positive hERG potentiation across the series.
Collapse
Affiliation(s)
- Robert A Pearlstein
- Global Discovery Chemistry, Computer-Aided Drug Discovery, Novartis Institutes for BioMedical Research, 181 Mass Ave., Cambridge, MA 02139, USA.
| | - Callum J Dickson
- Global Discovery Chemistry, Computer-Aided Drug Discovery, Novartis Institutes for BioMedical Research, 181 Mass Ave., Cambridge, MA 02139, USA
| | - Viktor Hornak
- Global Discovery Chemistry, Computer-Aided Drug Discovery, Novartis Institutes for BioMedical Research, 181 Mass Ave., Cambridge, MA 02139, USA
| |
Collapse
|
48
|
Mayne CG, Arcario MJ, Mahinthichaichan P, Baylon JL, Vermaas JV, Navidpour L, Wen PC, Thangapandian S, Tajkhorshid E. The cellular membrane as a mediator for small molecule interaction with membrane proteins. BIOCHIMICA ET BIOPHYSICA ACTA 2016; 1858:2290-2304. [PMID: 27163493 PMCID: PMC4983535 DOI: 10.1016/j.bbamem.2016.04.016] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/18/2016] [Revised: 04/26/2016] [Accepted: 04/27/2016] [Indexed: 01/05/2023]
Abstract
The cellular membrane constitutes the first element that encounters a wide variety of molecular species to which a cell might be exposed. Hosting a large number of structurally and functionally diverse proteins associated with this key metabolic compartment, the membrane not only directly controls the traffic of various molecules in and out of the cell, it also participates in such diverse and important processes as signal transduction and chemical processing of incoming molecular species. In this article, we present a number of cases where details of interaction of small molecular species such as drugs with the membrane, which are often experimentally inaccessible, have been studied using advanced molecular simulation techniques. We have selected systems in which partitioning of the small molecule with the membrane constitutes a key step for its final biological function, often binding to and interacting with a protein associated with the membrane. These examples demonstrate that membrane partitioning is not only important for the overall distribution of drugs and other small molecules into different compartments of the body, it may also play a key role in determining the efficiency and the mode of interaction of the drug with its target protein. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.
Collapse
Affiliation(s)
- Christopher G Mayne
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States.
| | - Mark J Arcario
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States; Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, United States; College of Medicine, University of Illinois at Urbana-Champaign, United States.
| | - Paween Mahinthichaichan
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States; Department of Biochemistry, University of Illinois at Urbana-Champaign, United States.
| | - Javier L Baylon
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States; Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, United States.
| | - Josh V Vermaas
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States; Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, United States.
| | - Latifeh Navidpour
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States.
| | - Po-Chao Wen
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States.
| | - Sundarapandian Thangapandian
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States; Department of Biochemistry, University of Illinois at Urbana-Champaign, United States.
| | - Emad Tajkhorshid
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, United States; Department of Biochemistry, University of Illinois at Urbana-Champaign, United States; Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, United States; College of Medicine, University of Illinois at Urbana-Champaign, United States.
| |
Collapse
|
49
|
Genheden S, Eriksson LA. Estimation of Liposome Penetration Barriers of Drug Molecules with All-Atom and Coarse-Grained Models. J Chem Theory Comput 2016; 12:4651-61. [PMID: 27541708 DOI: 10.1021/acs.jctc.6b00557] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Liposomes are common carriers of drug molecules, providing enhanced delivery and accumulation of hydrophilic agents or larger biomolecules. Molecular simulations can be used to estimate key features of the drug molecules upon interaction with the liposomes, such as penetration barriers and localization. Herein, we investigate several aspects of the computational estimation of penetration barriers, viz. the potential of mean force (PMFs) along a vector spanning the membrane. First, we provide an evaluation of the all-atom (AA) and coarse-grained (CG) parametrization of 5-aminolevulinic acid (5-ALA) and two of its alkyl esters by computing n-octanol/water partition coefficients. We find that the CG parametrization of the esters performs significantly better than the CG model of 5-ALA, highlighting the difficulty to coarse-grain small, polar molecules. However, the expected trend in partition coefficients is reproduced also with the CG models. Second, we compare PMFs in a small membrane slab described with either the AA or CG models. Here, we are able to reproduce the all-atom PMF of 5-ALA with CG. However, for the alkyl esters it is unfortunately not possible to correctly reproduce both the depth and the penetration barrier of the PMF seen in the AA simulations with any of the tested CG models. We argue that it is more important to choose a CG parametrization that reproduces the depth of the PMF. Third, we compare, using the CG model, PMFs in the membrane slab with PMFs in a large, realistic liposome. We find similar depths but slightly different penetration barriers most likely due to differences in the lipid density along the membrane axis. Finally, we compute PMFs in liposomes with three different lipid compositions. Unfortunately, differences in the PMFs could not be quantified, and it remains to be investigated to what extent liposome simulations can fully reproduce experimental findings.
Collapse
Affiliation(s)
- Samuel Genheden
- Department of Chemistry and Molecular Biology, University of Gothenburg , Box 462, SE-405 30 Göteborg, Sweden
| | - Leif A Eriksson
- Department of Chemistry and Molecular Biology, University of Gothenburg , Box 462, SE-405 30 Göteborg, Sweden
| |
Collapse
|
50
|
Lee CT, Comer J, Herndon C, Leung N, Pavlova A, Swift RV, Tung C, Rowley CN, Amaro RE, Chipot C, Wang Y, Gumbart JC. Simulation-Based Approaches for Determining Membrane Permeability of Small Compounds. J Chem Inf Model 2016; 56:721-33. [PMID: 27043429 DOI: 10.1021/acs.jcim.6b00022] [Citation(s) in RCA: 150] [Impact Index Per Article: 18.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Predicting the rate of nonfacilitated permeation of solutes across lipid bilayers is important to drug design, toxicology, and signaling. These rates can be estimated using molecular dynamics simulations combined with the inhomogeneous solubility-diffusion model, which requires calculation of the potential of mean force and position-dependent diffusivity of the solute along the transmembrane axis. In this paper, we assess the efficiency and accuracy of several methods for the calculation of the permeability of a model DMPC bilayer to urea, benzoic acid, and codeine. We compare umbrella sampling, replica exchange umbrella sampling, adaptive biasing force, and multiple-walker adaptive biasing force for the calculation of the transmembrane PMF. No definitive advantage for any of these methods in their ability to predict the membrane permeability coefficient Pm was found, provided that a sufficiently long equilibration is performed. For diffusivities, a Bayesian inference method was compared to a generalized Langevin method, both being sensitive to chosen parameters and the slow relaxation of membrane defects. Agreement within 1.5 log units of the computed Pm with experiment is found for all permeants and methods. Remaining discrepancies can likely be attributed to limitations of the force field as well as slowly relaxing collective movements within the lipid environment. Numerical calculations based on model profiles show that Pm can be reliably estimated from only a few data points, leading to recommendations for calculating Pm from simulations.
Collapse
Affiliation(s)
- Christopher T Lee
- Department of Chemistry and Biochemistry, University of California, San Diego , 9500 Gilman Drive, La Jolla, California 92093-0340, United States
| | - Jeffrey Comer
- Nanotechnology Innovation Center of Kansas State, Institute of Computational Comparative Medicine, Department of Anatomy and Physiology, Kansas State University , P-213 Mosier Hall, Manhattan, Kansas 66506, United States
| | - Conner Herndon
- School of Physics, Georgia Institute of Technology , 837 State Street, Atlanta, Georgia 30332, United States
| | - Nelson Leung
- Department of Physics, The Chinese University of Hong Kong , Shatin, Hong Kong SAR, China
| | - Anna Pavlova
- School of Physics, Georgia Institute of Technology , 837 State Street, Atlanta, Georgia 30332, United States
| | - Robert V Swift
- Department of Chemistry and Biochemistry, University of California, San Diego , 9500 Gilman Drive, La Jolla, California 92093-0340, United States
| | - Chris Tung
- Department of Physics, The Chinese University of Hong Kong , Shatin, Hong Kong SAR, China
| | - Christopher N Rowley
- Department of Chemistry, Memorial University of Newfoundland , St. John's, NL A1B 3X7 Canada
| | - Rommie E Amaro
- Department of Chemistry and Biochemistry, University of California, San Diego , 9500 Gilman Drive, La Jolla, California 92093-0340, United States
| | - Christophe Chipot
- Laboratoire International Associé Centre National de la Recherche Scientifique and University of Illinois at Urbana-Champaign, UMR n° 7565, Université de Lorraine , B.P. 70239, 54506 Vandœuvre-lès-Nancy, France.,Beckman Institute for Advanced Science and Technology and Department of Physics, University of Illinois at Urbana-Champaign , 405 North Mathews, Urbana, Illinois 61801, United States
| | - Yi Wang
- Department of Physics, The Chinese University of Hong Kong , Shatin, Hong Kong SAR, China.,Shenzhen Research Institute, The Chinese University of Hong Kong , Shatin, Hong Kong SAR, China
| | - James C Gumbart
- School of Physics, Georgia Institute of Technology , 837 State Street, Atlanta, Georgia 30332, United States.,School of Chemistry and Biochemistry, Georgia Institute of Technology , 901 Atlantic Drive NW, Atlanta, Georgia 30332, United States
| |
Collapse
|