1
|
Bannon MS, Ellena JF, Gourishankar AS, Marsh SR, Trevisan-Silva D, Sherman NE, Jourdan LJ, Gourdie RG, Letteri RA. Multi-site esterification: a tunable, reversible strategy to tailor therapeutic peptides for delivery. MOLECULAR SYSTEMS DESIGN & ENGINEERING 2024:d4me00072b. [PMID: 39281343 PMCID: PMC11395315 DOI: 10.1039/d4me00072b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/22/2024] [Accepted: 08/20/2024] [Indexed: 09/18/2024]
Abstract
Peptides are naturally potent and selective therapeutics with massive potential; however, low cell membrane permeability limits their clinical implementation, particularly for hydrophilic, anionic peptides with intracellular targets. To overcome this limitation, esterification of anionic carboxylic acids on therapeutic peptides can simultaneously increase hydrophobicity and net charge to facilitate cell internalization, whereafter installed esters can be cleaved hydrolytically to restore activity. To date, however, most esterified therapeutics contain either a single esterification site or multiple esters randomly incorporated on multiple sites. This investigation provides molecular engineering insight into how the number and position of esters installed onto the therapeutic peptide α carboxyl terminus 11 (αCT11, RPRPDDLEI) with 4 esterification sites affect hydrophobicity and the hydrolysis process that reverts the peptide to its original form. After installing methyl esters onto αCT11 using Fischer esterification, we isolated 5 distinct products and used 2D nuclear magnetic resonance spectroscopy, reverse-phase high performance liquid chromatography, and mass spectrometry to determine which residues were esterified in each and the resulting increase in hydrophobicity. We found esterifying the C-terminal isoleucine to impart the largest increase in hydrophobicity. Monitoring ester hydrolysis showed the C-terminal isoleucine ester to be the most hydrolytically stable, followed by the glutamic acid, whereas esters on aspartic acids hydrolyze rapidly. LC-MS revealed the formation of transient intramolecular aspartimides prior to hydrolysis to carboxylic acids. In vitro proof-of-concept experiments showed esterifying αCT11 to increase cell migration into a scratch, highlighting the potential of multi-site esterification as a tunable, reversible strategy to enable the delivery of therapeutic peptides.
Collapse
Affiliation(s)
- Mark S Bannon
- Department of Chemical Engineering, University of Virginia Charlottesville VA 22903 USA +1 434 243 3628
| | - Jeffrey F Ellena
- Biomolecular Magnetic Resonance Facility, School of Medicine, University of Virginia Charlottesville VA 22903 USA
| | - Aditi S Gourishankar
- Department of Chemical Engineering, University of Virginia Charlottesville VA 22903 USA +1 434 243 3628
| | - Spencer R Marsh
- Fralin Biomedical Institute, Virginia Tech Carillion School of Medicine Roanoke VA 24016 USA
| | - Dilza Trevisan-Silva
- Biomolecular Analysis Facility, School of Medicine, University of Virginia Charlottesville VA 22903 USA
| | - Nicholas E Sherman
- Biomolecular Analysis Facility, School of Medicine, University of Virginia Charlottesville VA 22903 USA
| | - L Jane Jourdan
- Fralin Biomedical Institute, Virginia Tech Carillion School of Medicine Roanoke VA 24016 USA
| | - Robert G Gourdie
- Fralin Biomedical Institute, Virginia Tech Carillion School of Medicine Roanoke VA 24016 USA
| | - Rachel A Letteri
- Department of Chemical Engineering, University of Virginia Charlottesville VA 22903 USA +1 434 243 3628
| |
Collapse
|
2
|
Pires CL, Moreno MJ. Improving the Accuracy of Permeability Data to Gain Predictive Power: Assessing Sources of Variability in Assays Using Cell Monolayers. MEMBRANES 2024; 14:157. [PMID: 39057665 PMCID: PMC11278619 DOI: 10.3390/membranes14070157] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/04/2024] [Revised: 07/05/2024] [Accepted: 07/11/2024] [Indexed: 07/28/2024]
Abstract
The ability to predict the rate of permeation of new compounds across biological membranes is of high importance for their success as drugs, as it determines their efficacy, pharmacokinetics, and safety profile. In vitro permeability assays using Caco-2 monolayers are commonly employed to assess permeability across the intestinal epithelium, with an extensive number of apparent permeability coefficient (Papp) values available in the literature and a significant fraction collected in databases. The compilation of these Papp values for large datasets allows for the application of artificial intelligence tools for establishing quantitative structure-permeability relationships (QSPRs) to predict the permeability of new compounds from their structural properties. One of the main challenges that hinders the development of accurate predictions is the existence of multiple Papp values for the same compound, mostly caused by differences in the experimental protocols employed. This review addresses the magnitude of the variability within and between laboratories to interpret its impact on QSPR modelling, systematically and quantitatively assessing the most common sources of variability. This review emphasizes the importance of compiling consistent Papp data and suggests strategies that may be used to obtain such data, contributing to the establishment of robust QSPRs with enhanced predictive power.
Collapse
Affiliation(s)
- Cristiana L. Pires
- Coimbra Chemistry Center—Institute of Molecular Sciences (CQC-IMS), University of Coimbra, 3004-535 Coimbra, Portugal
- Chemistry Department, Faculty of Science and Technology, University of Coimbra, 3004-535 Coimbra, Portugal
| | - Maria João Moreno
- Coimbra Chemistry Center—Institute of Molecular Sciences (CQC-IMS), University of Coimbra, 3004-535 Coimbra, Portugal
- Chemistry Department, Faculty of Science and Technology, University of Coimbra, 3004-535 Coimbra, Portugal
| |
Collapse
|
3
|
Muscat S, Errico S, Danani A, Chiti F, Grasso G. Leveraging Machine Learning-Guided Molecular Simulations Coupled with Experimental Data to Decipher Membrane Binding Mechanisms of Aminosterols. J Chem Theory Comput 2024. [PMID: 38979909 PMCID: PMC11447954 DOI: 10.1021/acs.jctc.4c00127] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
Understanding the molecular mechanisms of the interactions between specific compounds and cellular membranes is essential for numerous biotechnological applications, including targeted drug delivery, elucidation of the drug mechanism of action, pathogen identification, and novel antibiotic development. However, estimation of the free energy landscape associated with solute binding to realistic biological systems is still a challenging task. In this work, we leverage the Time-lagged Independent Component Analysis (TICA) in combination with neural networks (NN) through the Deep-TICA approach for determining the free energy associated with the membrane insertion processes of two natural aminosterol compounds, trodusquemine (TRO), and squalamine (SQ). These compounds are particularly noteworthy because they interact with the outer layer of neuron membranes, protecting them from the toxic action of misfolded proteins involved in neurodegenerative disorders, in both their monomeric and oligomeric forms. We demonstrate how this strategy could be used to generate an effective collective variable for describing solute absorption in the membrane and for estimating free energy landscape of translocation via on-the-fly probability enhanced sampling (OPES) method. In this context, the computational protocol allowed an exhaustive characterization of the aminosterol entry pathway into a neuron-like lipid bilayer. Furthermore, it provided accurate prediction of membrane binding affinities, in close agreement with the experimental binding data obtained by using fluorescently labeled aminosterols and large unilamellar vesicles (LUVs). The findings contribute significantly to our understanding of aminosterol entry pathways and aminosterol-lipid membrane interactions. Finally, the computational methods deployed in this study further demonstrate considerable potential for investigating membrane binding processes.
Collapse
Affiliation(s)
- Stefano Muscat
- Dalle Molle Institute for Artificial Intelligence IDSIA USI-SUPSI, Via la Santa 1 ,Lugano-Viganello 6962, Switzerland
| | - Silvia Errico
- Department of Experimental and Clinical Biomedical Sciences, Section of Biochemistry, University of Florence, Florence 50134, Italy
| | - Andrea Danani
- Dalle Molle Institute for Artificial Intelligence IDSIA USI-SUPSI, Via la Santa 1 ,Lugano-Viganello 6962, Switzerland
| | - Fabrizio Chiti
- Department of Experimental and Clinical Biomedical Sciences, Section of Biochemistry, University of Florence, Florence 50134, Italy
| | - Gianvito Grasso
- Dalle Molle Institute for Artificial Intelligence IDSIA USI-SUPSI, Via la Santa 1 ,Lugano-Viganello 6962, Switzerland
| |
Collapse
|
4
|
Zhu J, Graziotto ME, Cottam V, Hawtrey T, Adair LD, Trist BG, Pham NTH, Rouaen JRC, Ohno C, Heisler M, Vittorio O, Double KL, New EJ. Near-Infrared Ratiometric Fluorescent Probe for Detecting Endogenous Cu 2+ in the Brain. ACS Sens 2024; 9:2858-2868. [PMID: 38787339 DOI: 10.1021/acssensors.3c02549] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2024]
Abstract
Copper participates in a range of critical functions in the nervous system and human brain. Disturbances in brain copper content is strongly associated with neurological diseases. For example, changes in the level and distribution of copper are reported in neuroblastoma, Alzheimer's disease, and Lewy body disorders, such as Parkinson disease and dementia with Lewy bodies (DLB). There is a need for more sensitive techniques to measure intracellular copper levels to have a better understanding of the role of copper homeostasis in neuronal disorders. Here, we report a reaction-based near-infrared (NIR) ratiometric fluorescent probe CyCu1 for imaging Cu2+ in biological samples. High stability and selectivity of CyCu1 enabled the probe to be deployed as a sensor in a range of systems, including SH-SY5Y cells and neuroblastoma tumors. Furthermore, it can be used in plant cells, reporting on copper added to Arabidopsis roots. We also used CyCu1 to explore Cu2+ levels and distribution in post-mortem brain tissues from patients with DLB. We found significant decreases in Cu2+ content in the cytoplasm, neurons, and extraneuronal space in the degenerating substantia nigra in DLB compared with healthy age-matched control tissues. These findings enhance our understanding of Cu2+ dysregulation in Lewy body disorders. Our probe also shows promise as a photoacoustic imaging agent, with potential for applications in bimodal imaging.
Collapse
Affiliation(s)
- Jianping Zhu
- School of Chemistry, The University of Sydney, Sydney, NSW 2006, Australia
| | - Marcus E Graziotto
- School of Chemistry, The University of Sydney, Sydney, NSW 2006, Australia
| | - Veronica Cottam
- Brain and Mind Centre and School of Medical Sciences (Neuroscience), The University of Sydney, Sydney, NSW 2006, Australia
| | - Tom Hawtrey
- School of Chemistry, The University of Sydney, Sydney, NSW 2006, Australia
- Sydney Nano Institute, The University of Sydney, Sydney, NSW 2006, Australia
- Australian Research Council Centre of Excellence for Innovations in Peptide and Protein Science, The University of Sydney, Sydney, NSW 2006, Australia
| | - Liam D Adair
- School of Chemistry, The University of Sydney, Sydney, NSW 2006, Australia
- Sydney Nano Institute, The University of Sydney, Sydney, NSW 2006, Australia
- Australian Research Council Centre of Excellence for Innovations in Peptide and Protein Science, The University of Sydney, Sydney, NSW 2006, Australia
| | - Benjamin G Trist
- Brain and Mind Centre and School of Medical Sciences (Neuroscience), The University of Sydney, Sydney, NSW 2006, Australia
| | - Nguyen T H Pham
- Sydney Imaging, Core Research Facility, The University of Sydney, Sydney, NSW 2006, Australia
| | - Jourdin R C Rouaen
- Children's Cancer Institute, Lowy Cancer Research Centre, University of New South Wales Sydney, Randwick, NSW 2052, Australia
| | - Carolyn Ohno
- School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW 2006, Australia
| | - Marcus Heisler
- School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW 2006, Australia
| | - Orazio Vittorio
- Children's Cancer Institute, Lowy Cancer Research Centre, University of New South Wales Sydney, Randwick, NSW 2052, Australia
- School of Biomedical Sciences, University of New South Wales, Kensington, NSW 2031, Australia
| | - Kay L Double
- Brain and Mind Centre and School of Medical Sciences (Neuroscience), The University of Sydney, Sydney, NSW 2006, Australia
| | - Elizabeth J New
- School of Chemistry, The University of Sydney, Sydney, NSW 2006, Australia
- Sydney Nano Institute, The University of Sydney, Sydney, NSW 2006, Australia
- Australian Research Council Centre of Excellence for Innovations in Peptide and Protein Science, The University of Sydney, Sydney, NSW 2006, Australia
| |
Collapse
|
5
|
Dutta A, Sepehri A, Lazaridis T. Putative Pore Structures of Amyloid β 25-35 in Lipid Bilayers. Biochemistry 2023; 62:2549-2558. [PMID: 37582191 DOI: 10.1021/acs.biochem.3c00323] [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] [Indexed: 08/17/2023]
Abstract
The amyloid β peptide aggregates to form extracellular plaques in the brains of Alzheimer's disease patients. Certain of its fragments have been found to have similar properties to those of the full-length peptide. The best-studied of these is 25-35, which aggregates into fibrils, is toxic to neurons, and forms ion channels in synthetic lipid bilayers. Here, we investigate possible pore-forming structures of oligomers of this peptide in a POPC/POPG membrane. We consider octameric and decameric β-barrels of different topology, strand orientation, and shear, evaluate their stability in an implicit membrane model, and subject the best models to multimicrosecond all-atom molecular dynamics simulations. We find two decameric structures that are kinetically stable in membranes on this time scale: an imperfectly closed antiparallel β-barrel with K28 in the pore lumen and a short parallel β-barrel with K28 toward the membrane interface. Both structures exhibit dehydrated gaps in the pore lumen, which are larger for the antiparallel barrel. Based on these results, the experimental cation selectivity, the dependence of ion channel activity on voltage direction, and certain mutation data, the parallel model seems more compatible with experimental data.
Collapse
Affiliation(s)
- Ankita Dutta
- Department of Chemistry, City College of New York/CUNY, 160 Convent Avenue, New York, New York 10031, United States
- Graduate Program in Biochemistry, The Graduate Center, City University of New York, 365 Fifth Avenue, New York, New York 10016, United States
| | - Aliasghar Sepehri
- Department of Chemistry, City College of New York/CUNY, 160 Convent Avenue, New York, New York 10031, United States
| | - Themis Lazaridis
- Department of Chemistry, City College of New York/CUNY, 160 Convent Avenue, New York, New York 10031, United States
- Graduate Programs in Chemistry, Biochemistry, and Physics The Graduate Center, City University of New York, 365 Fifth Avenue, New York, New York 10016, United States
| |
Collapse
|
6
|
Šlachtová V, Bellová S, La-Venia A, Galeta J, Dračínský M, Chalupský K, Dvořáková A, Mertlíková-Kaiserová H, Rukovanský P, Dzijak R, Vrabel M. Triazinium Ligation: Bioorthogonal Reaction of N1-Alkyl 1,2,4-Triazinium Salts. Angew Chem Int Ed Engl 2023; 62:e202306828. [PMID: 37436086 DOI: 10.1002/anie.202306828] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2023] [Revised: 07/12/2023] [Accepted: 07/12/2023] [Indexed: 07/13/2023]
Abstract
The development of reagents that can selectively react in complex biological media is an important challenge. Here we show that N1-alkylation of 1,2,4-triazines yields the corresponding triazinium salts, which are three orders of magnitude more reactive in reactions with strained alkynes than the parent 1,2,4-triazines. This powerful bioorthogonal ligation enables efficient modification of peptides and proteins. The positively charged N1-alkyl triazinium salts exhibit favorable cell permeability, which makes them superior for intracellular fluorescent labeling applications when compared to analogous 1,2,4,5-tetrazines. Due to their high reactivity, stability, synthetic accessibility and improved water solubility, the new ionic heterodienes represent a valuable addition to the repertoire of existing modern bioorthogonal reagents.
Collapse
Affiliation(s)
- Veronika Šlachtová
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám. 2, 16000, Prague, Czech Republic
| | - Simona Bellová
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám. 2, 16000, Prague, Czech Republic
| | - Agustina La-Venia
- Current address: Instituto de Química Rosario, Facultad de Ciencias Bioquímicas y Farmacéuticas, Universidad Nacional de Rosario-CONICET, Suipacha 531, S2002LRK, Rosario, Argentina
| | - Juraj Galeta
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám. 2, 16000, Prague, Czech Republic
| | - Martin Dračínský
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám. 2, 16000, Prague, Czech Republic
| | - Karel Chalupský
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám. 2, 16000, Prague, Czech Republic
| | - Alexandra Dvořáková
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám. 2, 16000, Prague, Czech Republic
| | - Helena Mertlíková-Kaiserová
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám. 2, 16000, Prague, Czech Republic
| | - Peter Rukovanský
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám. 2, 16000, Prague, Czech Republic
| | - Rastislav Dzijak
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám. 2, 16000, Prague, Czech Republic
| | - Milan Vrabel
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám. 2, 16000, Prague, Czech Republic
| |
Collapse
|
7
|
Aguilar-Toalá JE, Vidal-Limon A, Liceaga AM, Zambrano-Zaragoza ML, Quintanar-Guerrero D. Application of Molecular Dynamics Simulations to Determine Interactions between Canary Seed ( Phalaris canariensis L.) Bioactive Peptides and Skin-Aging Enzymes. Int J Mol Sci 2023; 24:13420. [PMID: 37686226 PMCID: PMC10487734 DOI: 10.3390/ijms241713420] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2023] [Revised: 08/22/2023] [Accepted: 08/28/2023] [Indexed: 09/10/2023] Open
Abstract
Food bioactive peptides are well recognized for their health benefits such as antimicrobial, antioxidant, and antihypertensive benefits, among others. Their drug-like behavior has led to their potential use in targeting skin-related aging factors like the inhibition of enzymes related with the skin-aging process. In this study, canary seed peptides (CSP) after simulated gastrointestinal digestion (<3 kDa) were fractioned by RP-HPLC and their enzyme-inhibition activity towards elastase and tyrosinase was evaluated in vitro. CSP inhibited elastase (IC50 = 6.2 mg/mL) and tyrosinase (IC50 = 6.1 mg/mL), while the hydrophobic fraction-VI (0.2 mg/mL) showed the highest inhibition towards elastase (93%) and tyrosinase (67%). The peptide fraction with the highest inhibition was further characterized by a multilevel in silico workflow, including physicochemical descriptor calculations, antioxidant activity predictions, and molecular dynamics-ensemble docking towards elastase and tyrosinase. To gain insights into the skin permeation process during molecular dynamics simulations, based on their docking scores, five peptides (GGWH, VPPH, EGLEPNHRVE, FLPH, and RPVNKYTPPQ) were identified to have favorable intermolecular interactions, such as hydrogen bonding of polar residues (W, H, and K) to lipid polar groups and 2-3 Å van der Waals close contact of hydrophobic aliphatic residues (P, V, and L). These interactions can play a critical role for the passive insertion of peptides into stratum corneum model skin-membranes, suggesting a promising application of CSP for skin-aging treatments.
Collapse
Affiliation(s)
- José E. Aguilar-Toalá
- Departamento de Ciencias de la Alimentación, División de Ciencias Biológicas y de la Salud, Universidad Autónoma Metropolitana, Unidad Lerma. Av. de las Garzas 10. Col. El Panteón, Lerma de Villada 52005, Estado de México, Mexico;
| | - Abraham Vidal-Limon
- Red de Estudios Moleculares Avanzados, Instituto de Ecología A.C. (INECOL), Carretera Antigua a Coatepec 351, Xalapa 91073, Veracruz, Mexico
| | - Andrea M. Liceaga
- Protein Chemistry and Bioactive Peptides Laboratory, Purdue University, 745 Agriculture Mall, West Lafayette, IN 47907, USA
| | - Maria L. Zambrano-Zaragoza
- Laboratorio de Procesos de Transformación y Tecnologías Emergentes de Alimentos-UIM, FES-Cuautitlán, Universidad Nacional Autónoma de México, Cuautitlán Izcalli 54714, Estado de México, Mexico;
| | - David Quintanar-Guerrero
- Laboratorio de Posgrado en Tecnología Farmacéutica, FES-Cuautitlán, Universidad Nacional Autónoma de México, Av. 1o de Mayo s/n, Cuautitlán Izcalli 54714, Estado de México, Mexico;
| |
Collapse
|
8
|
Harada R, Morita R, Shigeta Y. Free-Energy Profiles for Membrane Permeation of Compounds Calculated Using Rare-Event Sampling Methods. J Chem Inf Model 2023; 63:259-269. [PMID: 36574612 DOI: 10.1021/acs.jcim.2c01097] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
The free-energy profile of a compound is an essential measurement in evaluating the membrane permeation process by means of theoretical methods. Computationally, molecular dynamics (MD) simulation allows the free-energy profile calculation. However, MD simulations frequently fail to sample membrane permeation because they are rare events induced in longer timescales than the accessible timescale of MD, leading to an insufficient conformational search to calculate an incorrect free-energy profile. To achieve a sufficient conformational search, several enhanced sampling methods have been developed and elucidated the membrane permeation process. In addition to these enhanced sampling methods, we proposed a simple yet powerful free-energy calculation of a compound for the membrane permeation process based on originally rare-event sampling methods developed by us. Our methods have a weak dependency on external biases and their optimizations to promote the membrane permeation process. Based on distributed computing, our methods only require the selection of initial structures and their conformational resampling, whereas the enhanced sampling methods may be required to adjust external biases. Furthermore, our methods efficiently search membrane permeation processes with simple scripts without modifying any MD program. As demonstrations, we calculated the free-energy profiles of seven linear compounds for their membrane permeation based on a hybrid conformational search using two rare-event sampling methods, that is, (1) parallel cascade selection MD (PaCS-MD) and (2) outlier flooding method (OFLOOD), combined with a Markov state model (MSM) construction. In the first step, PaCS-MD generated initial membrane permeation paths of a compound. In the second step, OFLOOD expanded the unsearched conformational area around the initial paths, allowing for a broad conformational search. Finally, the trajectories were employed to construct reliable MSMs, enabling correct free-energy profile calculations. Furthermore, we estimated the membrane permeability coefficients of all compounds by constructing the reliable MSMs for their membrane permeation. In conclusion, the calculated coefficients were qualitatively correlated with the experimental measurements (correlation coefficient (R2) = 0.8689), indicating that the hybrid conformational search successfully calculated the free-energy profiles and membrane permeability coefficients of the seven compounds.
Collapse
Affiliation(s)
- Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki305-8577, Japan
| | - Rikuri Morita
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki305-8577, Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki305-8577, Japan
| |
Collapse
|
9
|
Abbas G, Cardenas AE, Elber R. The Structures of Heterogeneous Membranes and Their Interactions with an Anticancer Peptide: A Molecular Dynamics Study. Life (Basel) 2022; 12:1473. [PMID: 36294908 PMCID: PMC9604715 DOI: 10.3390/life12101473] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2022] [Revised: 09/16/2022] [Accepted: 09/20/2022] [Indexed: 12/02/2022] Open
Abstract
We conduct molecular dynamics simulations of model heterogeneous membranes and their interactions with a 24-amino acid peptide-NAF-144-67. NAF-144-67 is an anticancer peptide that selectively permeates and kills malignant cells; it does not permeate normal cells. We examine three membranes with different binary mixtures of lipids, DOPC-DOPA, DOPC-DOPS, and DOPC-DOPE, with a single peptide embedded in each as models for the diversity of biological membranes. We illustrate that the peptide organization in the membrane depends on the types of nearby phospholipids and is influenced by the charge and size of the head groups. The present study sheds light on early events of permeation and the mechanisms by which an amphiphilic peptide crosses from an aqueous solution to a hydrophobic membrane. Understanding the translocation mechanism is likely to help the design of new permeants.
Collapse
Affiliation(s)
- Ghulam Abbas
- National Center for Bioinformatics, Quaid-i-Azam University, Islamabad 45320, Pakistan or
| | - Alfredo E. Cardenas
- Oden Institute for Computational and Engineering Sciences, University of Texas at Austin, Austin, TX 78712, USA
| | - Ron Elber
- Oden Institute for Computational and Engineering Sciences, University of Texas at Austin, Austin, TX 78712, USA
- Department of Chemistry, University of Texas at Austin, Austin, TX 78712, USA
| |
Collapse
|
10
|
Povilaitis SC, Fathizadeh A, Kogan M, Elber R, Webb LJ. Design of Peptides for Membrane Insertion: The Critical Role of Charge Separation. J Phys Chem B 2022; 126:6454-6463. [PMID: 35997537 PMCID: PMC9541189 DOI: 10.1021/acs.jpcb.2c04615] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A physical understanding of membrane permeation and translocation by small, positively charged molecules can illuminate cell penetrating peptide mechanisms of entry and inform drug design. We have previously investigated the permeation of the doubly charged peptide WKW and proposed a defect-assisted permeation mechanism where a small molecule with +2 charge can achieve a metastable state spanning the bilayer by forming a membrane defect with charges stabilized by phospholipid phosphate groups. Here, we investigate the membrane permeation of two doubly charged peptides, WWK and WWWK, with charges separated by different lengths. Through complementary experiments and molecular dynamics simulations, we show that membrane permeation was an order of magnitude more favorable when charges were separated by an ∼2-3 Å greater distance on WWWK compared to WWK. These results agree with the previously proposed defect-assisted permeation mechanism, where a greater distance between positive charges would require a less extreme membrane defect to stabilize the membrane-spanning metastable state. We discuss the implications of these results in understanding the membrane permeation of cell-penetrating peptides and other small, positively charged membrane permeants.
Collapse
Affiliation(s)
- Sydney C. Povilaitis
- Department of Chemistry, University of Texas at Austin, Austin, TX 78712, United States
| | - Arman Fathizadeh
- Oden Institute for Computational Engineering and Science, University of Texas at Austin, Austin, Texas 78712, United States
| | - Molly Kogan
- Department of Chemistry, University of Texas at Austin, Austin, TX 78712, United States
| | - Ron Elber
- Department of Chemistry, University of Texas at Austin, Austin, TX 78712, United States
- Oden Institute for Computational Engineering and Science, University of Texas at Austin, Austin, Texas 78712, United States
| | - Lauren J. Webb
- Department of Chemistry, University of Texas at Austin, Austin, TX 78712, United States
| |
Collapse
|
11
|
Molecular dynamics simulations of a central nervous system-penetrant drug AZD3759 with lipid bilayer. J Mol Model 2022; 28:261. [DOI: 10.1007/s00894-022-05266-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2022] [Accepted: 08/11/2022] [Indexed: 10/15/2022]
|
12
|
Cardenas AE, Drexler CI, Nechushtai R, Mittler R, Friedler A, Webb LJ, Elber R. Peptide Permeation across a Phosphocholine Membrane: An Atomically Detailed Mechanism Determined through Simulations and Supported by Experimentation. J Phys Chem B 2022; 126:2834-2849. [PMID: 35388695 PMCID: PMC9074375 DOI: 10.1021/acs.jpcb.1c10966] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Cell-penetrating peptides (CPPs) facilitate translocation across biological membranes and are of significant biological and medical interest. Several CPPs can permeate into specific cells and organelles. We examine the incorporation and translocation of a novel anticancer CPP in a dioleoylphosphatidylcholine (DOPC) lipid bilayer membrane. The peptide, NAF-144-67, is a short fragment of a transmembrane protein, consisting of hydrophobic N-terminal and charged C-terminal segments. Experiments using fluorescently labeled NAF-144-67 in ∼100 nm DOPC vesicles and atomically detailed simulations conducted with Milestoning support a model in which a significant barrier for peptide-membrane entry is found at the interface between the aqueous solution and membrane. The initial step is the insertion of the N-terminal segment and the hydrophobic helix into the membrane, passing the hydrophilic head groups. Both experiments and simulations suggest that the free energy difference in the first step of the permeation mechanism in which the hydrophobic helix crosses the phospholipid head groups is -0.4 kcal mol-1 slightly favoring motion into the membrane. Milestoning calculations of the mean first passage time and the committor function underscore the existence of an early polar barrier followed by a diffusive barrierless motion in the lipid tail region. Permeation events are coupled to membrane fluctuations that are examined in detail. Our study opens the way to investigate in atomistic resolution the molecular mechanism, kinetics, and thermodynamics of CPP permeation to diverse membranes.
Collapse
Affiliation(s)
- Alfredo E. Cardenas
- Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX 78712, USA
| | - Chad I. Drexler
- Department of Chemistry, The University of Texas at Austin, Austin, TX 78712, USA
| | - Rachel Nechushtai
- The Alexander Silberman Institute of Life Sciences, The Hebrew University of Jerusalem, Edmond J. Safra Campus, Givat-Ram, Jerusalem 91904, Israel
| | - Ron Mittler
- The Department of Surgery, University of Missouri School of Medicine. Christopher S. Bond Life Sciences Center, University of Missouri. 1201 Rollins St, Columbia, MO 65201, USA
| | - Assaf Friedler
- The Institute of Chemistry, The Hebrew University of Jerusalem, Edmond J. Safra Campus, Givat-Ram, Jerusalem 91904, Israel
| | - Lauren J. Webb
- Department of Chemistry, The University of Texas at Austin, Austin, TX 78712, USA
| | - Ron Elber
- Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX 78712, USA
- Department of Chemistry, The University of Texas at Austin, Austin, TX 78712, USA
| |
Collapse
|
13
|
Chandra A, Datta A. A Peptide-Based Fluorescent Sensor for Anionic Phospholipids. ACS OMEGA 2022; 7:10347-10354. [PMID: 35382295 PMCID: PMC8973094 DOI: 10.1021/acsomega.1c06981] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Accepted: 02/25/2022] [Indexed: 06/14/2023]
Abstract
Anionic phospholipids are key cell signal mediators. The distribution of these lipids on the cell membrane and intracellular organelle membranes guides the recruitment of signaling proteins leading to the regulation of cellular processes. Hence, fluorescent sensors that can detect anionic phospholipids within living cells can provide a handle into revealing molecular mechanisms underlying lipid-mediated signal regulation. A major challenge in the detection of anionic phospholipids is related to the presence of these phospholipids mostly in the inner leaflet of the plasma membrane and in the membranes of intracellular organelles. Hence, cell-permeable sensors would provide an advantage by enabling the rapid detection and tracking of intracellular pools of anionic phospholipids. We have developed a peptide-based, cell-permeable, water-soluble, and ratiometric fluorescent sensor that entered cells within 15 min of incubation via the endosomal machinery and showed punctate labeling in the cytoplasm. The probe could also be introduced into living cells via lipofection, which allows bypassing of endosomal uptake, to image anionic phospholipids in the cell membrane. We validated the ability of the sensor toward detection of intracellular anionic phospholipids by colocalization studies with a fluorescently tagged lipid and a protein-based anionic phospholipid sensor. Further, the sensor could image the externalization of anionic phospholipids during programmed cell death, indicating the ability of the probe toward detection of both intra- and extracellular anionic phospholipids based on the biological context.
Collapse
|
14
|
Saur JS, Wirtz SN, Schilling NA, Krismer B, Peschel A, Grond S. Distinct Lugdunins from a New Efficient Synthesis and Broad Exploitation of Its MRSA-Antimicrobial Structure. J Med Chem 2021; 64:4034-4058. [PMID: 33779184 DOI: 10.1021/acs.jmedchem.0c02170] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
A new solid-phase peptide synthesis and bioprofiling of the antimicrobial activity of lugdunin, a fibupeptide, enable a comprehensive structure-activity relationship (SAR) study (MRSA Staphylococcus aureus). Distinct lugdunin analogues with variation of the three important amino acids Val2, Trp3, and Leu4 are readily available based on the established high-output synthesis. This efficient synthesis concept takes advantage of the presynthesized thiazolidine building block. To gain further knowledge of SAR, d-Val2, and d-Leu4 were replaced with aliphatic amino acids. For l-Trp3 derivatization, a set of non-natural aromatic amino acids with manifold substitution and annulation patterns precisely shows structural imperatives, starting from the exchange of d-Val6 → d-Trp6 with a 2-fold improved biological activity. d-Trp6-lugdunin analogues with additional variation of d-Val2 and d-Leu4 residues were designed and synthesized followed by antimicrobial profiling. For the first time, these SAR studies deliver valuable information on the tolerance of other amino acids to d-Val2, l-Trp3, and d-Leu4 in the sequence of lugdunin.
Collapse
Affiliation(s)
- Julian S Saur
- Institute of Organic Chemistry, Eberhard Karls University Tuebingen, Auf der Morgenstelle 18, 72076 Tuebingen, Germany
| | - Sebastian N Wirtz
- Institute of Organic Chemistry, Eberhard Karls University Tuebingen, Auf der Morgenstelle 18, 72076 Tuebingen, Germany
| | - Nadine A Schilling
- Institute of Organic Chemistry, Eberhard Karls University Tuebingen, Auf der Morgenstelle 18, 72076 Tuebingen, Germany
| | - Bernhard Krismer
- Cluster of Excellence EXC 2124 Controlling Microbes to Fight Infections, Eberhard Karls University Tuebingen, Auf der Morgenstelle 28, 72076 Tuebingen, Germany.,Interfaculty Institute of Microbiology and Infection Medicine, German Center for Infection Research (DZIF), Eberhard Karls University Tuebingen, Auf der Morgenstelle 28, 72076 Tuebingen, Germany.,German Center for Infection Research (DZIF), Eberhard Karls University Tuebingen, Auf der Morgenstelle 28, 72076 Tuebingen, Germany
| | - Andreas Peschel
- Cluster of Excellence EXC 2124 Controlling Microbes to Fight Infections, Eberhard Karls University Tuebingen, Auf der Morgenstelle 28, 72076 Tuebingen, Germany.,Interfaculty Institute of Microbiology and Infection Medicine, German Center for Infection Research (DZIF), Eberhard Karls University Tuebingen, Auf der Morgenstelle 28, 72076 Tuebingen, Germany.,German Center for Infection Research (DZIF), Eberhard Karls University Tuebingen, Auf der Morgenstelle 28, 72076 Tuebingen, Germany
| | - Stephanie Grond
- Institute of Organic Chemistry, Eberhard Karls University Tuebingen, Auf der Morgenstelle 18, 72076 Tuebingen, Germany.,Cluster of Excellence EXC 2124 Controlling Microbes to Fight Infections, Eberhard Karls University Tuebingen, Auf der Morgenstelle 28, 72076 Tuebingen, Germany
| |
Collapse
|
15
|
Elber R, Fathizadeh A, Ma P, Wang H. Modeling molecular kinetics with Milestoning. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1512] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Ron Elber
- Department of Chemistry, The Oden Institute for Computational Engineering and Sciences University of Texas at Austin Austin Texas USA
| | - Arman Fathizadeh
- The Oden Institute for Computational Engineering and Sciences University of Texas at Austin Austin Texas USA
| | - Piao Ma
- Department of Chemistry University of Texas at Austin Austin Texas USA
| | - Hao Wang
- The Oden Institute for Computational Engineering and Sciences University of Texas at Austin Austin Texas USA
| |
Collapse
|
16
|
Kornspan JD, Kosower NS, Vaisid T, Katzhandler J, Rottem S. Novel synthetic lipopeptides derived from Mycoplasma hyorhinis upregulate calpastatin in SH-SY5Y neuroblastoma cells and induce a neuroprotective effect against amyloid-β-peptide toxicity. FEMS Microbiol Lett 2020; 367:5824629. [PMID: 32329786 DOI: 10.1093/femsle/fnaa073] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2019] [Accepted: 04/22/2020] [Indexed: 11/14/2022] Open
Abstract
Previously, we showed that contamination of SH-SY5Y neuroblastoma cells by Mycoplasma hyorhinis strains NDMh and MCLD leads to increased levels of calpastatin (the endogenous, specific inhibitor of the Ca2+-dependent protease calpain), resulting in inhibition of calpain activation. We have found that the increased calpastatin level is promoted by the lipoprotein fraction (MhLpp) of the mycoplasmal membrane. Here, we present MhLpp-based novel synthetic lipopeptides that induce upregulation of calpastatin in SH-SY5Y neuroblastoma cells, leading to protection of the treated cells against Ca2+/amyloid-β-peptide toxicity. These lipopeptides present a new class of promising agents against calpain-induced cell toxicity.
Collapse
Affiliation(s)
- Jonathan D Kornspan
- Department of Microbiology and Molecular Genetics, The Hebrew University-Hadassah Medical School, Jerusalem 91120, Israel
| | - Nechama S Kosower
- Department of Human Molecular Genetics & Biochemistry, Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv 69978, Israel
| | - Tali Vaisid
- Department of Human Molecular Genetics & Biochemistry, Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv 69978, Israel
| | | | - Shlomo Rottem
- Department of Microbiology and Molecular Genetics, The Hebrew University-Hadassah Medical School, Jerusalem 91120, Israel
| |
Collapse
|
17
|
Dahlgren D, Sjöblom M, Hedeland M, Lennernäs H. The In Vivo Effect of Transcellular Permeation Enhancers on the Intestinal Permeability of Two Peptide Drugs Enalaprilat and Hexarelin. Pharmaceutics 2020; 12:pharmaceutics12020099. [PMID: 31991924 PMCID: PMC7076382 DOI: 10.3390/pharmaceutics12020099] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2019] [Revised: 01/20/2020] [Accepted: 01/22/2020] [Indexed: 12/14/2022] Open
Abstract
Permeation enhancers like sodium dodecyl sulfate (SDS) and caprate increase the intestinal permeability of small model peptide compounds, such as enalaprilat (349 Da). However, their effects remain to be investigated for larger low-permeability peptide drugs, such as hexarelin (887 Da). The objective of this single-pass perfusion study in rat was to investigate the effect of SDS at 5 mg/mL and of caprate administered at different luminal concentrations (5, 10, and 20 mg/mL) and pH (6.5 and 7.4). The small intestinal permeability of enalaprilat increased by 8- and 9-fold with SDS at 5 mg/mL and with caprate at 10 and 20 mg/mL but only at pH 7.4, where the free dissolved caprate concentration is higher than at pH 6.5 (5 vs. 2 mg/mL). Neither SDS nor caprate at any of the investigated luminal concentrations enhanced absorption of the larger peptide hexarelin. These results show that caprate requires doses above its saturation concentration (a reservoir suspension) to enhance absorption, most likely because dissolved caprate itself is rapidly absorbed. The absent effect on hexarelin may partly explain why the use of permeation enhancers for enabling oral peptide delivery has largely failed to evolve from in vitro evaluations into approved oral products. It is obvious that more innovative and effective drug delivery strategies are needed for this class of drugs.
Collapse
Affiliation(s)
- David Dahlgren
- Department of Pharmacy, Uppsala University, 751 23 Uppsala, Sweden;
| | - Markus Sjöblom
- Department of Neuroscience, Uppsala University, 751 23 Uppsala, Sweden;
| | - Mikael Hedeland
- Department of Medicinal Chemistry, Uppsala University, 751 23 Uppsala, Sweden;
- National Veterinary Institute (SVA), 751 89 Uppsala, Sweden
| | - Hans Lennernäs
- Department of Pharmacy, Uppsala University, 751 23 Uppsala, Sweden;
- Correspondence: ; Tel.: +46-18-471-4317
| |
Collapse
|
18
|
Guan X, Wei DQ, Hu D. Free Energy Calculations on the Water-Chain-Assisted and the Dehydration Mechanisms of Transmembrane Ion Permeation. J Chem Theory Comput 2019; 16:700-710. [DOI: 10.1021/acs.jctc.9b00671] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
|
19
|
Fathizadeh A, Kogan M, Anderson CM, Webb LJ, Elber R. Defect-Assisted Permeation Through a Phospholipid Membrane: Experimental and Computational Study of the Peptide WKW. J Phys Chem B 2019; 123:6792-6798. [PMID: 31304755 PMCID: PMC6687544 DOI: 10.1021/acs.jpcb.9b05414] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
We investigate membrane permeation by the peptide WKW that is amidated at its C-terminus and therefore carries a positive charge of +2. To facilitate an efficient calculation, we introduce a novel set of simple coarse variables that measure permeation depth and membrane distortion. The phospholipid head groups shift toward the center of the membrane, following the permeating peptide, and create a defect that assists permeation. The Milestoning algorithm was used in the new coarse space to compute the free-energy profile and the mean first passage time. The barrier was lower than expected from a simple continuum estimate. This behavior is consistent with the known behavior of positively charged cell-penetrating peptides, and is explained by a detailed mechanism of defect formation and propagation revealed by the simulations.
Collapse
Affiliation(s)
- Arman Fathizadeh
- Oden Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin. TX, 78712
| | - Molly Kogan
- Department of Chemistry, University of Texas at Austin, Austin TX, 78712
| | - Cari M. Anderson
- Department of Chemistry, University of Texas at Austin, Austin TX, 78712
| | - Lauren J. Webb
- Department of Chemistry, University of Texas at Austin, Austin TX, 78712
| | - Ron Elber
- Oden Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin. TX, 78712
- Department of Chemistry, University of Texas at Austin, Austin TX, 78712
| |
Collapse
|
20
|
Lomize AL, Pogozheva ID. Physics-Based Method for Modeling Passive Membrane Permeability and Translocation Pathways of Bioactive Molecules. J Chem Inf Model 2019; 59:3198-3213. [PMID: 31259555 DOI: 10.1021/acs.jcim.9b00224] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Assessment of permeability is a critical step in the drug development process for selection of drug candidates with favorable ADME properties. We have developed a novel physics-based method for fast computational modeling of passive permeation of diverse classes of molecules across lipid membranes. The method is based on heterogeneous solubility-diffusion theory and operates with all-atom 3D structures of solutes and the anisotropic solvent model of the lipid bilayer characterized by transbilayer profiles of dielectric and hydrogen bonding capacity parameters. The optimal translocation pathway of a solute is determined by moving an ensemble of representative conformations of the molecule through the dioleoyl-phosphatidylcholine (DOPC) bilayer and optimizing their rotational orientations in every point of the transmembrane trajectory. The method calculates (1) the membrane-bound state of the solute molecule; (2) free energy profile of the solute along the permeation pathway; and (3) the permeability coefficient obtained by integration over the transbilayer energy profile and assuming a constant size-dependent diffusivity along the membrane normal. The accuracy of the predictions was evaluated against experimental permeability coefficients measured in pure lipid membranes (for 78 compounds, R2 was 0.88 and rmse was 1.15 log units), PAMPA-DS (for 280 compounds, R2 was 0.75 and rmse was 1.59 log units), BBB (for 182 compounds, R2 was 0.69 and rmse was 0.87 log units), and Caco-2/MDCK assays (for 165 compounds, R2 was 0.52 and rmse was 0.89 log units).
Collapse
Affiliation(s)
- Andrei L Lomize
- Department of Medicinal Chemistry, College of Pharmacy , University of Michigan , 428 Church Street , Ann Arbor , Michigan 48109-1065 , United States
| | - Irina D Pogozheva
- Department of Medicinal Chemistry, College of Pharmacy , University of Michigan , 428 Church Street , Ann Arbor , Michigan 48109-1065 , United States
| |
Collapse
|
21
|
Pires F, Geraldo VPN, Rodrigues B, Granada-Flor AD, de Almeida RFM, Oliveira ON, Victor BL, Machuqueiro M, Raposo M. Evaluation of EGCG Loading Capacity in DMPC Membranes. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:6771-6781. [PMID: 31006246 DOI: 10.1021/acs.langmuir.9b00372] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Catechins are molecules with potential use in different pathologies such as diabetes and cancer, but their pharmaceutical applications are often hindered by their instability in the bloodstream. This issue can be circumvented using liposomes as their nanocarriers for in vivo delivery. In this work, we studied the molecular details of (-)-epigallocatechin-3-gallate (EGCG) interacting with 1,2-dimyristoyl- sn-glycero-3-phosphocholine (DMPC) monolayer/bilayer systems to understand the catechin loading ability and liposome stability, using experimental and computational techniques. The molecular dynamics simulations show the EGCG molecules deep inside the lipid bilayer, positioned below the lipid ester groups, generating a concentration-dependent lipid condensation. This effect was also inferred from the surface pressure isotherms of DMPC monolayers. In the polarization-modulated infrared reflection absorption spectra assays, the predominant effect at higher concentrations of EGCG (e.g., 20 mol %) was an increase in lipid tail disorder. The steady-state fluorescence data confirmed this disordered state, indicating that the catechin-induced liposome aggregation outweighs the condensation effects. Therefore, by adding more than 10 mol % EGCG to the liposomes, a destabilization of the vesicles occurs with the ensuing release of entrapped catechins. The loading capacity for DMPC seems to be limited by its disordered lipid arrangements, typical of a fluid phase. To further increase the clinical usefulness of liposomes, lipid bilayers with more stable and organized assemblies should be employed to avoid aggregation at large concentrations of catechin.
Collapse
Affiliation(s)
- Filipa Pires
- Departamento de Física, CEFITEC, Faculdade de Ciências e Tecnologia , Universidade Nova de Lisboa , 2829-516 Caparica , Portugal
| | - Vananélia P N Geraldo
- Instituto de Física de São Carlos , Universidade de São Paulo , 13560-970 Sao Carlos , Brazil
| | - Bárbara Rodrigues
- Departamento de Física, CEFITEC, Faculdade de Ciências e Tecnologia , Universidade Nova de Lisboa , 2829-516 Caparica , Portugal
| | - António de Granada-Flor
- Departamento de Química e Bioquímica, Centro de Química e Bioquímica , Faculdade de Ciências da Universidade de Lisboa , Campo Grande, 1749-016 Lisboa , Portugal
| | - Rodrigo F M de Almeida
- Departamento de Química e Bioquímica, Centro de Química e Bioquímica , Faculdade de Ciências da Universidade de Lisboa , Campo Grande, 1749-016 Lisboa , Portugal
| | - Osvaldo N Oliveira
- Instituto de Física de São Carlos , Universidade de São Paulo , 13560-970 Sao Carlos , Brazil
| | - Bruno L Victor
- Departamento de Química e Bioquímica, Centro de Química e Bioquímica , Faculdade de Ciências da Universidade de Lisboa , Campo Grande, 1749-016 Lisboa , Portugal
| | - Miguel Machuqueiro
- Departamento de Química e Bioquímica, Centro de Química e Bioquímica , Faculdade de Ciências da Universidade de Lisboa , Campo Grande, 1749-016 Lisboa , Portugal
| | - Maria Raposo
- Departamento de Física, CEFITEC, Faculdade de Ciências e Tecnologia , Universidade Nova de Lisboa , 2829-516 Caparica , Portugal
| |
Collapse
|
22
|
Cardenas AE, Anderson CM, Elber R, Webb LJ. Partition of Positively and Negatively Charged Tryptophan Ions in Membranes with Inverted Phospholipid Heads: Simulations and Experiments. J Phys Chem B 2019; 123:3272-3281. [PMID: 30912653 DOI: 10.1021/acs.jpcb.9b00754] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A joint experimental and computational study illustrates that the partitioning of positively and negatively charged tryptophan in a phospholipid bilayer is significantly altered by a reversal in the head group dipole arrangement. Experiments were conducted using tryptophan as a fluorescent reporter of its local environment. Based on the experimental design in a recent publication ( Anderson , C. M. ; Cardenas , A. ; Elber , R. ; Webb , L. J. J. Phys. Chem. B 2018 , 123 , 170 - 179 ), we were able to determine that the arrangement of the head group dipole altered the degree of partitioning of charged tryptophan in the lipid bilayer. In parallel, atomically detailed simulations were performed for the two membrane systems. The simulation results are in accord with the experimental findings and support a simple molecular partition mechanism of electrostatic interactions with the head groups, glycerol linkers, and interfacial water dipoles.
Collapse
|
23
|
Anderson CM, Cardenas A, Elber R, Webb LJ. Preferential Equilibrium Partitioning of Positively Charged Tryptophan into Phosphatidylcholine Bilayer Membranes. J Phys Chem B 2019; 123:170-179. [PMID: 30481465 PMCID: PMC6331081 DOI: 10.1021/acs.jpcb.8b09872] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2018] [Revised: 11/19/2018] [Indexed: 11/29/2022]
Abstract
The interactions between small molecules and lipid bilayers play a critical role in the function of cellular membranes. Understanding how a small molecule interacts with the lipid bilayer differently based on its charge reveals primordial mechanisms of transport across membranes and assists in the design of drug molecules that can penetrate cells. We have previously reported that tryptophan permeated through a phosphatidylcholine lipid bilayer membrane at a faster rate when it was positively charged (Trp+) than when negatively charged (Trp-), which corresponded to a lower potential of mean force (PMF) barrier determined through simulations. In this report, we demonstrate that Trp+ partitions into the lipid bilayer membrane to a greater degree than Trp- by interacting with the ester linkage of a phosphatidylcholine lipid, where it is stabilized by the electron withdrawing glycerol functional group. These results are in agreement with tryptophan's known role as an anchor for transmembrane proteins, though the tendency for binding of a positively charged tryptophan is surprising. We discuss the implications of our results on the mechanisms of unassisted permeation and penetration of small molecules within and across lipid bilayer membranes based on molecular charge, shape, and molecular interactions within the bilayer structure.
Collapse
Affiliation(s)
- Cari M. Anderson
- Department
of Chemistry, Institute for Computational Engineering and Sciences, Institute for Cellular
and Molecular Biology, Texas Materials Institute, The University of Texas at Austin, 2506 Speedway STOP A5300, Austin, Texas 78712, United States
| | - Alfredo Cardenas
- Department
of Chemistry, Institute for Computational Engineering and Sciences, Institute for Cellular
and Molecular Biology, Texas Materials Institute, The University of Texas at Austin, 2506 Speedway STOP A5300, Austin, Texas 78712, United States
| | - Ron Elber
- Department
of Chemistry, Institute for Computational Engineering and Sciences, Institute for Cellular
and Molecular Biology, Texas Materials Institute, The University of Texas at Austin, 2506 Speedway STOP A5300, Austin, Texas 78712, United States
| | - Lauren J. Webb
- Department
of Chemistry, Institute for Computational Engineering and Sciences, Institute for Cellular
and Molecular Biology, Texas Materials Institute, The University of Texas at Austin, 2506 Speedway STOP A5300, Austin, Texas 78712, United States
| |
Collapse
|
24
|
Fathizadeh A, Elber R. Ion Permeation through a Phospholipid Membrane: Transition State, Path Splitting, and Calculation of Permeability. J Chem Theory Comput 2019; 15:720-730. [PMID: 30474968 PMCID: PMC6467798 DOI: 10.1021/acs.jctc.8b00882] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We investigate the thermodynamics and kinetics of the permeation of a potassium ion through a phospholipid membrane. We illustrate that the conventional reaction coordinate (the position of the ion along the normal to the membrane plane) is insufficient to capture essential elements of the process. It is necessary to add coarse variables that measure membrane distortion. New coarse variables are suggested, and a two-dimensional coarse-space is proposed to describe the permeation. We illustrate path splitting and two transition states of comparable barrier heights. The alternative pathways differ by the extent of water solvation of the ion-phosphate pairs. The permeation process cannot be described by a local one-dimensional reaction coordinate, and a network formulation is more appropriate. We use Milestoning with Voronoi tessellation in two dimensions to quantify the equilibrium and rate of the permeation of the positively charged ion. The permeation coefficient is computed and compared favorably to experiment.
Collapse
Affiliation(s)
- Arman Fathizadeh
- Institute for Computational Engineering and Sciences , University of Texas at Austin , Austin , Texas 78712 , United States
| | - Ron Elber
- Institute for Computational Engineering and Sciences , University of Texas at Austin , Austin , Texas 78712 , United States
- Department of Chemistry , University of Texas at Austin , Austin , Texas 78712 , United States
| |
Collapse
|
25
|
Jagger BR, Lee CT, Amaro RE. Quantitative Ranking of Ligand Binding Kinetics with a Multiscale Milestoning Simulation Approach. J Phys Chem Lett 2018; 9:4941-4948. [PMID: 30070844 PMCID: PMC6443090 DOI: 10.1021/acs.jpclett.8b02047] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/07/2023]
Abstract
Efficient prediction and ranking of small molecule binders by their kinetic ( kon and koff) and thermodynamic ( Δ G) properties can be a valuable metric for drug lead optimization, as these quantities are often indicators of in vivo efficacy. We have previously described a hybrid molecular dynamics, Brownian dynamics, and milestoning model, Simulation Enabled Estimation of Kinetic Rates (SEEKR), that can predict kon's, koff's, and Δ G's. Here we demonstrate the effectiveness of this approach for ranking a series of seven small molecule compounds for the model system, β-cyclodextrin, based on predicted kon's and koff's. We compare our results using SEEKR to experimentally determined rates as well as rates calculated using long time scale molecular dynamics simulations and show that SEEKR can effectively rank the compounds by koff and Δ G with reduced computational cost. We also provide a discussion of convergence properties and sensitivities of calculations with SEEKR to establish "best practices" for its future use.
Collapse
Affiliation(s)
- Benjamin R Jagger
- Department of Chemistry and Biochemistry , University of California, San Diego , 9500 Gilman Drive , La Jolla , California 92093-0340 , United States
| | - Christopher T Lee
- Department of Chemistry and Biochemistry , University of California, San Diego , 9500 Gilman Drive , La Jolla , California 92093-0340 , United States
| | - Rommie E Amaro
- Department of Chemistry and Biochemistry , University of California, San Diego , 9500 Gilman Drive , La Jolla , California 92093-0340 , United States
| |
Collapse
|
26
|
Abstract
The kinetics of biochemical and biophysical events determined the course of life processes and attracted considerable interest and research. For example, modeling of biological networks and cellular responses relies on the availability of information on rate coefficients. Atomically detailed simulations hold the promise of supplementing experimental data to obtain a more complete kinetic picture. However, simulations at biological time scales are challenging. Typical computer resources are insufficient to provide the ensemble of trajectories at the correct length that is required for straightforward calculations of time scales. In the last years, new technologies emerged that make atomically detailed simulations of rate coefficients possible. Instead of computing complete trajectories from reactants to products, these approaches launch a large number of short trajectories at different positions. Since the trajectories are short, they are computed trivially in parallel on modern computer architecture. The starting and termination positions of the short trajectories are chosen, following statistical mechanics theory, to enhance efficiency. These trajectories are analyzed. The analysis produces accurate estimates of time scales as long as hours. The theory of Milestoning that exploits the use of short trajectories is discussed, and several applications are described.
Collapse
|
27
|
Cao Z, Bian Y, Hu G, Zhao L, Kong Z, Yang Y, Wang J, Zhou Y. Bias-Exchange Metadynamics Simulation of Membrane Permeation of 20 Amino Acids. Int J Mol Sci 2018; 19:E885. [PMID: 29547563 PMCID: PMC5877746 DOI: 10.3390/ijms19030885] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2018] [Revised: 03/11/2018] [Accepted: 03/12/2018] [Indexed: 11/16/2022] Open
Abstract
Thermodynamics of the permeation of amino acids from water to lipid bilayers is an important first step for understanding the mechanism of cell-permeating peptides and the thermodynamics of membrane protein structure and stability. In this work, we employed bias-exchange metadynamics simulations to simulate the membrane permeation of all 20 amino acids from water to the center of a dipalmitoylphosphatidylcholine (DPPC) membrane (consists of 256 lipids) by using both directional and torsion angles for conformational sampling. The overall accuracy for the free energy profiles obtained is supported by significant correlation coefficients (correlation coefficient at 0.5-0.6) between our results and previous experimental or computational studies. The free energy profiles indicated that (1) polar amino acids have larger free energy barriers than nonpolar amino acids; (2) negatively charged amino acids are the most difficult to enter into the membrane; and (3) conformational transitions for many amino acids during membrane crossing is the key for reduced free energy barriers. These results represent the first set of simulated free energy profiles of membrane crossing for all 20 amino acids.
Collapse
Affiliation(s)
- Zanxia Cao
- Shandong Provincial Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China.
- College of Physics and Electronic Information, Dezhou University, Dezhou 253023, China.
| | - Yunqiang Bian
- Shandong Provincial Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China.
| | - Guodong Hu
- Shandong Provincial Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China.
- College of Physics and Electronic Information, Dezhou University, Dezhou 253023, China.
| | - Liling Zhao
- Shandong Provincial Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China.
- College of Physics and Electronic Information, Dezhou University, Dezhou 253023, China.
| | - Zhenzhen Kong
- Shandong Provincial Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China.
- College of Life Science, Shandong Normal University, Jinan 250014, China.
| | - Yuedong Yang
- Institute for Glycomics and School of Information and Communication Technology, Griffith University, Parklands Dr, Southport, QLD 4222, Australia.
- School of Data and Computer Science, Sun Yat-sen University, Guangzhou 510275, China.
| | - Jihua Wang
- Shandong Provincial Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China.
- College of Physics and Electronic Information, Dezhou University, Dezhou 253023, China.
| | - Yaoqi Zhou
- Shandong Provincial Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China.
- Institute for Glycomics and School of Information and Communication Technology, Griffith University, Parklands Dr, Southport, QLD 4222, Australia.
| |
Collapse
|
28
|
Pokhrel N, Maibaum L. Free Energy Calculations of Membrane Permeation: Challenges Due to Strong Headgroup-Solute Interactions. J Chem Theory Comput 2018; 14:1762-1771. [PMID: 29406707 DOI: 10.1021/acs.jctc.7b01159] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Understanding how different classes of molecules move across biological membranes is a prerequisite to predicting a solute's permeation rate, which is a critical factor in the fields of drug design and pharmacology. We use biased molecular dynamics computer simulations to calculate and compare the free energy profiles of translocation of several small molecules across 1,2-dioleoyl- sn-glycero-3-phosphocholine (DOPC) lipid bilayers as a first step toward determining the most efficient method for free energy calculations. We study the translocation of arginine, a sodium ion, alanine, and a single water molecule using the metadynamics, umbrella sampling, and replica exchange umbrella sampling techniques. Within the fixed lengths of our simulations, we find that all methods produce similar results for charge-neutral permeants, but not for polar or positively charged molecules. We identify the long relaxation time scale of electrostatic interactions between lipid headgroups and the solute to be the principal cause of this difference and show that this slow process can lead to an erroneous dependence of computed free energy profiles on the initial system configuration. We demonstrate the use of committor analysis to validate the proper sampling of the presumed transition state, which in our simulations is achieved only in replica exchange calculations. On the basis of these results we provide some useful guidance to perform and evaluate free energy calculations of membrane permeation.
Collapse
Affiliation(s)
- Nihit Pokhrel
- Department of Chemistry , University of Washington , Seattle , Washington 98195 , United States
| | - Lutz Maibaum
- Department of Chemistry , University of Washington , Seattle , Washington 98195 , United States
| |
Collapse
|
29
|
Lee BL, Kuczera K. Simulating the free energy of passive membrane permeation for small molecules. MOLECULAR SIMULATION 2017. [DOI: 10.1080/08927022.2017.1407029] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Brent L. Lee
- Department of Chemistry, The University of Kansas , Lawrence, KS, USA
| | - Krzysztof Kuczera
- Department of Chemistry, The University of Kansas , Lawrence, KS, USA
- Department of Molecular Biosciences, The University of Kansas , Lawrence, KS, USA
| |
Collapse
|
30
|
Wei C, Pohorille A. Sequence-Dependent Interfacial Adsorption and Permeation of Dipeptides across Phospholipid Membranes. J Phys Chem B 2017; 121:9859-9867. [PMID: 28982244 DOI: 10.1021/acs.jpcb.7b08238] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
We investigate permeation of three blocked dipeptides with different side chain polarity across a phospholipid membrane and their behavior at the water-membrane interface by way of molecular dynamics simulations. Hydrophilic serine-serine dipeptide is found to desorb from the interface to aqueous phase, whereas hydrophobic phenylalanine-leucine and amphiphilic serine-leucine tend to accumulate at the interface with a free energy minimum of -3 kcal/mol. All three dipeptides exhibit free energy barriers to permeation across the membrane located at the center of the bilayer. The height of the barrier is strongly sequence dependent and increases with the dipeptide polarity. It is equal to 3.5, 6.4, and 10.0 kcal/mol for phenylalanine-leucine, serine-leucine, and serine-serine, respectively. The corresponding permeability coefficients are equal to 4.6 × 10-3, 4.5 × 10-5, and 8.7 × 10-8 cm/s. The apparent insensitivity of membrane permeability to hydrophobicity of dipeptides, found in some experiments, is attributed to neglecting corrections for unstirred water layers near membrane surface, which are significant for hydrophobic species. Different hydrophobicity of the dipeptides also influences their conformations and orientations, both at the interface and inside the membrane. In particular, penetration of hydrophilic serine-serine dipeptide causes the formation of water-filled defects in the bilayer. These results are relevant to the delivery of peptide-based therapeutic agents.
Collapse
Affiliation(s)
- Chenyu Wei
- NASA Ames Research Center, Mail Stop 239-4, Moffett Field, California 94035, United States.,Department of Pharmaceutical Chemistry, University of California, San Francisco , San Francisco, California 94143, United States
| | - Andrew Pohorille
- NASA Ames Research Center, Mail Stop 239-4, Moffett Field, California 94035, United States.,Department of Pharmaceutical Chemistry, University of California, San Francisco , San Francisco, California 94143, United States
| |
Collapse
|
31
|
Shrestha R, Anderson CM, Cardenas AE, Elber R, Webb LJ. Direct Measurement of the Effect of Cholesterol and 6-Ketocholestanol on the Membrane Dipole Electric Field Using Vibrational Stark Effect Spectroscopy Coupled with Molecular Dynamics Simulations. J Phys Chem B 2017; 121:3424-3436. [PMID: 28071910 PMCID: PMC5398937 DOI: 10.1021/acs.jpcb.6b09007] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Biological membranes are heterogeneous structures with complex electrostatic profiles arising from lipids, sterols, membrane proteins, and water molecules. We investigated the effect of cholesterol and its derivative 6-ketocholestanol (6-kc) on membrane electrostatics by directly measuring the dipole electric field (F⃗d) within lipid bilayers containing cholesterol or 6-kc at concentrations of 0-40 mol% through the vibrational Stark effect (VSE). We found that adding low concentrations of cholesterol, up to ∼10 mol %, increases F⃗d, while adding more cholesterol up to 40 mol% lowers F⃗d. In contrast, we measured a monotonic increase in F⃗d as 6-kc concentration increased. We propose that this membrane electric field is affected by multiple factors: the polarity of the sterol molecules, the reorientation of the phospholipid dipole due to sterol, and the impact of the sterol on hydrogen bonding with surface water. We used molecular dynamics simulations to examine the distribution of phospholipids, sterol, and helix in bilayers containing these sterols. At low concentrations, we observed clustering of sterols near the vibrational probe whereas at high concentrations, we observed spatial correlation between the positions of the sterol molecules. This work demonstrates how a one-atom difference in a sterol changes the physicochemical and electric field properties of the bilayer.
Collapse
Affiliation(s)
- Rebika Shrestha
- Department of Chemistry, ‡Institute for Cell and Molecular Biology, §Center for Nano- and Molecular Science and Technology, and ∥Institute for Computational Engineering and Sciences, The University of Texas at Austin , Austin, Texas 78712, United States
| | - Cari M Anderson
- Department of Chemistry, ‡Institute for Cell and Molecular Biology, §Center for Nano- and Molecular Science and Technology, and ∥Institute for Computational Engineering and Sciences, The University of Texas at Austin , Austin, Texas 78712, United States
| | - Alfredo E Cardenas
- Department of Chemistry, ‡Institute for Cell and Molecular Biology, §Center for Nano- and Molecular Science and Technology, and ∥Institute for Computational Engineering and Sciences, The University of Texas at Austin , Austin, Texas 78712, United States
| | - Ron Elber
- Department of Chemistry, ‡Institute for Cell and Molecular Biology, §Center for Nano- and Molecular Science and Technology, and ∥Institute for Computational Engineering and Sciences, The University of Texas at Austin , Austin, Texas 78712, United States
| | - Lauren J Webb
- Department of Chemistry, ‡Institute for Cell and Molecular Biology, §Center for Nano- and Molecular Science and Technology, and ∥Institute for Computational Engineering and Sciences, The University of Texas at Austin , Austin, Texas 78712, United States
| |
Collapse
|
32
|
Hub JS, Awasthi N. Probing a Continuous Polar Defect: A Reaction Coordinate for Pore Formation in Lipid Membranes. J Chem Theory Comput 2017; 13:2352-2366. [PMID: 28376619 DOI: 10.1021/acs.jctc.7b00106] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Various biophysical processes involve the formation of aqueous pores over lipid membranes, including processes of membrane fusion, antimicrobial peptide activity, lipid flip-flop, and membrane permeation. Reliable and efficient free-energy calculations of pore formation using molecular dynamics simulations remained challenging due to the lack of good reaction coordinates (RCs) for pore formation. We present a new RC for pore formation that probes the formation and rupture of a continuous polar defect over the membrane. Potential of mean force (PMF) calculations along the new RC rapidly converge and exhibit no hysteresis between pore-opening and pore-closing pathways, in contrast to calculations based on previous RCs. We show that restraints along the new RC may restrain the system tightly to the transition state of pore formation, rationalizing the absence of hysteresis. We observe that the PMF of pore formation in a tension-free membrane of dimyristoylphosphatidylcholine (DMPC) reveals a free-energy barrier for pore nucleation, confirming a long-hypothesized metastable prepore state. We test the influence of the lipid force field, the cutoff distance used for Lennard-Jones interactions, and the lateral membrane size on the free energies of pore formation. In contrast to PMF calculations based on previous RCs, we find that such parameters have a relatively small influence on the free energies of pore nucleation. However, the metastability of the open pore in DMPC may depend on such parameters. The RC has been implemented into an extension of the GROMACS simulation software. The new RC allows for reliable and computationally efficient free-energy calculations of pore formation in lipid membranes.
Collapse
Affiliation(s)
- Jochen S Hub
- Institute for Microbiology and Genetics, University of Göttingen , Justus-von-Liebig Weg 11, 37077 Goettingen, Germany
| | - Neha Awasthi
- Institute for Microbiology and Genetics, University of Göttingen , Justus-von-Liebig Weg 11, 37077 Goettingen, Germany
| |
Collapse
|
33
|
Casas J, Ibarguren M, Álvarez R, Terés S, Lladó V, Piotto SP, Concilio S, Busquets X, López DJ, Escribá PV. G protein-membrane interactions II: Effect of G protein-linked lipids on membrane structure and G protein-membrane interactions. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2017; 1859:1526-1535. [PMID: 28411171 DOI: 10.1016/j.bbamem.2017.04.005] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/11/2017] [Revised: 03/21/2017] [Accepted: 04/09/2017] [Indexed: 10/19/2022]
Abstract
G proteins often bear myristoyl, palmitoyl and isoprenyl moieties, which favor their association with the membrane and their accumulation in G Protein Coupled Receptor-rich microdomains. These lipids influence the biophysical properties of membranes and thereby modulate G protein binding to bilayers. In this context, we showed here that geranylgeraniol, but neither myristate nor palmitate, increased the inverted hexagonal (HII) phase propensity of phosphatidylethanolamine-containing membranes. While myristate and palmitate preferentially associated with phosphatidylcholine membranes, geranylgeraniol favored nonlamellar-prone membranes. In addition, Gαi1 monomers had a higher affinity for lamellar phases, while Gβγ and Gαβγ showed a marked preference for nonlamellar prone membranes. Moreover, geranylgeraniol enhanced the binding of G protein dimers and trimers to phosphatidylethanolamine-containing membranes, yet it decreased that of monomers. By contrast, both myristate and palmitate increased the Gαi1 preference for lamellar membranes. Palmitoylation reinforced the binding of the monomer to PC membranes and myristoylation decreased its binding to PE-enriched bilayer. Finally, binding of dimers and trimers to lamellar-prone membranes was decreased by palmitate and myristate, but it was increased in nonlamellar-prone bilayers. These results demonstrate that co/post-translational G protein lipid modifications regulate the membrane lipid structure and that they influence the physico-chemical properties of membranes, which in part explains why G protein subunits sort to different plasma membrane domains. This article is part of a Special Issue entitled: Membrane Lipid Therapy: Drugs Targeting Biomembranes edited by Pablo V. Escribá.
Collapse
Affiliation(s)
- Jesús Casas
- Laboratory of Molecular Cell Biomedicine, Department of Biology, University of the Balearic Islands, E-07122 Palma de Mallorca, Spain
| | - Maitane Ibarguren
- Laboratory of Molecular Cell Biomedicine, Department of Biology, University of the Balearic Islands, E-07122 Palma de Mallorca, Spain; Lipopharma Therapeutics, S.L., ParcBit, 07121 Palma de Mallorca, Spain
| | - Rafael Álvarez
- Laboratory of Molecular Cell Biomedicine, Department of Biology, University of the Balearic Islands, E-07122 Palma de Mallorca, Spain
| | - Silvia Terés
- Laboratory of Molecular Cell Biomedicine, Department of Biology, University of the Balearic Islands, E-07122 Palma de Mallorca, Spain; Lipopharma Therapeutics, S.L., ParcBit, 07121 Palma de Mallorca, Spain
| | - Victoria Lladó
- Laboratory of Molecular Cell Biomedicine, Department of Biology, University of the Balearic Islands, E-07122 Palma de Mallorca, Spain
| | - Stefano P Piotto
- Department of Pharmacy, University of Salerno, Via Ponte don Melillo, 84084 Fisciano, SA, Italy
| | - Simona Concilio
- Department of Industrial Engineering, University of Salerno, Via Ponte don Melillo, 84084 Fisciano, SA, Italy
| | - Xavier Busquets
- Laboratory of Molecular Cell Biomedicine, Department of Biology, University of the Balearic Islands, E-07122 Palma de Mallorca, Spain
| | - David J López
- Laboratory of Molecular Cell Biomedicine, Department of Biology, University of the Balearic Islands, E-07122 Palma de Mallorca, Spain; Lipopharma Therapeutics, S.L., ParcBit, 07121 Palma de Mallorca, Spain.
| | - Pablo V Escribá
- Laboratory of Molecular Cell Biomedicine, Department of Biology, University of the Balearic Islands, E-07122 Palma de Mallorca, Spain
| |
Collapse
|
34
|
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
|
35
|
McHugh SM, Rogers JR, Solomon SA, Yu H, Lin YS. Computational methods to design cyclic peptides. Curr Opin Chem Biol 2016; 34:95-102. [PMID: 27592259 DOI: 10.1016/j.cbpa.2016.08.004] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2016] [Revised: 08/08/2016] [Accepted: 08/10/2016] [Indexed: 10/21/2022]
Abstract
Cyclic peptides (CPs) are promising modulators of protein-protein interactions (PPIs), but their application remains challenging. It is currently difficult to predict the structures and bioavailability of CPs. The ability to design CPs using computer modeling would greatly facilitate the development of CPs as potent PPI modulators for fundamental studies and as potential therapeutics. Herein, we describe computational methods to generate CP libraries for virtual screening, as well as current efforts to accurately predict the conformations adopted by CPs. These advances are making it possible to envision robust computational design of active CPs. However, unique properties of CPs pose significant challenges associated with sampling CP conformational space and accurately describing CP energetics. These major obstacles to structure prediction likely must be solved before robust design of active CPs can be reliably achieved.
Collapse
Affiliation(s)
- Sean M McHugh
- Department of Chemistry, Tufts University, Medford, MA 02155, United States
| | - Julia R Rogers
- Department of Chemistry, Tufts University, Medford, MA 02155, United States
| | - Sarah A Solomon
- Department of Chemistry, Tufts University, Medford, MA 02155, United States
| | - Hongtao Yu
- Department of Chemistry, Tufts University, Medford, MA 02155, United States
| | - Yu-Shan Lin
- Department of Chemistry, Tufts University, Medford, MA 02155, United States.
| |
Collapse
|
36
|
Lee BL, Kuczera K, Middaugh CR, Jas GS. Permeation of the three aromatic dipeptides through lipid bilayers: Experimental and computational study. J Chem Phys 2016; 144:245103. [DOI: 10.1063/1.4954241] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Affiliation(s)
- Brent L. Lee
- Department of Chemistry, The University of Kansas, Lawrence, Kansas 66045, USA
| | - Krzysztof Kuczera
- Department of Chemistry, The University of Kansas, Lawrence, Kansas 66045, USA
- Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66045, USA
| | - C. Russell Middaugh
- Department of Pharmaceutical Chemistry, The University of Kansas, Lawrence, Kansas 66047, USA
| | - Gouri S. Jas
- Department of Pharmaceutical Chemistry, The University of Kansas, Lawrence, Kansas 66047, USA
| |
Collapse
|
37
|
Awasthi N, Hub JS. Simulations of Pore Formation in Lipid Membranes: Reaction Coordinates, Convergence, Hysteresis, and Finite-Size Effects. J Chem Theory Comput 2016; 12:3261-9. [DOI: 10.1021/acs.jctc.6b00369] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Neha Awasthi
- Institute for Microbiology
and Genetics, Georg-August-Universität, Justus-von-Liebig Weg 11, 37077 Göttingen, Germany
| | - Jochen S. Hub
- Institute for Microbiology
and Genetics, Georg-August-Universität, Justus-von-Liebig Weg 11, 37077 Göttingen, Germany
| |
Collapse
|
38
|
Permeability across lipid membranes. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1858:2254-2265. [PMID: 27085977 DOI: 10.1016/j.bbamem.2016.03.032] [Citation(s) in RCA: 188] [Impact Index Per Article: 23.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/30/2016] [Revised: 03/28/2016] [Accepted: 03/29/2016] [Indexed: 11/22/2022]
Abstract
Molecular permeation through lipid membranes is a fundamental biological process that is important for small neutral molecules and drug molecules. Precise characterization of free energy surface and diffusion coefficients along the permeation pathway is required in order to predict molecular permeability and elucidate the molecular mechanisms of permeation. Several recent technical developments, including improved molecular models and efficient sampling schemes, are illustrated in this review. For larger penetrants, explicit consideration of multiple collective variables, including orientational, conformational degrees of freedom, are required to be considered in addition to the distance from the membrane center along the membrane normal. Although computationally demanding, this method can provide significant insights into the molecular mechanisms of permeation for molecules of medical and pharmaceutical importance. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.
Collapse
|
39
|
Jas GS, Rentchler EC, Słowicka AM, Hermansen JR, Johnson CK, Middaugh CR, Kuczera K. Reorientation Motion and Preferential Interactions of a Peptide in Denaturants and Osmolyte. J Phys Chem B 2016; 120:3089-99. [DOI: 10.1021/acs.jpcb.6b00028] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Gouri S. Jas
- Department
of Pharmaceutical Chemistry, The University of Kansas, Lawrence, Kansas 66047, United States
| | - Eric C. Rentchler
- Department
of Chemistry, The University of Kansas, Lawrence, Kansas 66045, United States
| | - Agnieszka M. Słowicka
- Institute
of Fundamental Technological Research, Polish Academy of Sciences, ul. Pawińskiego 5B, 02-106 Warsaw, Poland
| | - John R. Hermansen
- School
of Medicine, Central University of the Caribbean, Bayamon, PR 00956, United States
| | - Carey K. Johnson
- Department
of Chemistry, The University of Kansas, Lawrence, Kansas 66045, United States
| | - C. Russell Middaugh
- Department
of Pharmaceutical Chemistry, The University of Kansas, Lawrence, Kansas 66047, United States
| | - Krzysztof Kuczera
- Department
of Chemistry, The University of Kansas, Lawrence, Kansas 66045, United States
- Department
of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66045, United States
| |
Collapse
|
40
|
Aristoff D, Bello-Rivas JM, Elber R. A MATHEMATICAL FRAMEWORK FOR EXACT MILESTONING. MULTISCALE MODELING & SIMULATION : A SIAM INTERDISCIPLINARY JOURNAL 2016; 14:301-322. [PMID: 27239166 PMCID: PMC4879838 DOI: 10.1137/15m102157x] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
We give a mathematical framework for Exact Milestoning, a recently introduced algorithm for mapping a continuous time stochastic process into a Markov chain or semi-Markov process that can be efficiently simulated and analyzed. We generalize the setting of Exact Milestoning and give explicit error bounds for the error in the Milestoning equation for mean first passage times.
Collapse
Affiliation(s)
- David Aristoff
- Department of Mathematics, Colorado State University, Fort Collins, CO
| | - Juan M Bello-Rivas
- Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin, TX
| | - Ron Elber
- Institute for Computational Engineering and Sciences, Department of Chemistry, University of Texas at Austin, Austin, TX
| |
Collapse
|
41
|
Teixeira VH, Vila-Viçosa D, Reis PBPS, Machuqueiro M. pKa Values of Titrable Amino Acids at the Water/Membrane Interface. J Chem Theory Comput 2016; 12:930-4. [DOI: 10.1021/acs.jctc.5b01114] [Citation(s) in RCA: 51] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Vitor H. Teixeira
- Centro de Química
e Bioquímica, Departamento de Química e Bioquímica,
Faculdade de Ciências, Universidade de Lisboa, 1749-016 Lisboa, Portugal
| | - Diogo Vila-Viçosa
- Centro de Química
e Bioquímica, Departamento de Química e Bioquímica,
Faculdade de Ciências, Universidade de Lisboa, 1749-016 Lisboa, Portugal
| | - Pedro B. P. S. Reis
- Centro de Química
e Bioquímica, Departamento de Química e Bioquímica,
Faculdade de Ciências, Universidade de Lisboa, 1749-016 Lisboa, Portugal
| | - Miguel Machuqueiro
- Centro de Química
e Bioquímica, Departamento de Química e Bioquímica,
Faculdade de Ciências, Universidade de Lisboa, 1749-016 Lisboa, Portugal
| |
Collapse
|
42
|
Oruç T, Küçük SE, Sezer D. Lipid bilayer permeation of aliphatic amine and carboxylic acid drugs: rates of insertion, translocation and dissociation from MD simulations. Phys Chem Chem Phys 2016; 18:24511-25. [DOI: 10.1039/c6cp05278a] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The overabundance of drugs containing aliphatic amine and carboxylic acid groups is rationalized in terms of their membrane permeability.
Collapse
Affiliation(s)
- Tuğçe Oruç
- Faculty of Engineering and Natural Sciences
- Sabanc University
- 34956 Istanbul
- Turkey
| | - Sami Emre Küçük
- Faculty of Engineering and Natural Sciences
- Sabanc University
- 34956 Istanbul
- Turkey
| | - Deniz Sezer
- Faculty of Engineering and Natural Sciences
- Sabanc University
- 34956 Istanbul
- Turkey
| |
Collapse
|
43
|
Allolio C, Baxova K, Vazdar M, Jungwirth P. Guanidinium Pairing Facilitates Membrane Translocation. J Phys Chem B 2015; 120:143-53. [DOI: 10.1021/acs.jpcb.5b10404] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Christoph Allolio
- Institute
of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, CZ-166 10 Prague 6, Czech Republic
- Institut
für Physikalische and Theoretische Chemie, Universität Regensburg, 93040 Regensburg, Germany
| | - Katarina Baxova
- Institute
of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, CZ-166 10 Prague 6, Czech Republic
| | - Mario Vazdar
- Institut
Rudjer
Bošković, Bijenička
cesta 54, 10000 Zagreb, Croatia
| | - Pavel Jungwirth
- Institute
of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, CZ-166 10 Prague 6, Czech Republic
- Department
of Physics, Tampere University of Technology, P.O. Box 692, FI-33101 Tampere, Finland
| |
Collapse
|
44
|
Awoonor-Williams E, Rowley CN. Molecular simulation of nonfacilitated membrane permeation. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2015; 1858:1672-87. [PMID: 26706099 DOI: 10.1016/j.bbamem.2015.12.014] [Citation(s) in RCA: 83] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/13/2015] [Revised: 12/05/2015] [Accepted: 12/09/2015] [Indexed: 12/29/2022]
Abstract
This is a review. Non-electrolytic compounds typically cross cell membranes by passive diffusion. The rate of permeation is dependent on the chemical properties of the solute and the composition of the lipid bilayer membrane. Predicting the permeability coefficient of a solute is important in pharmaceutical chemistry and toxicology. Molecular simulation has proven to be a valuable tool for modeling permeation of solutes through a lipid bilayer. In particular, the solubility-diffusion model has allowed for the quantitative calculation of permeability coefficients. The underlying theory and computational methods used to calculate membrane permeability are reviewed. We also discuss applications of these methods to examine the permeability of solutes and the effect of membrane composition on permeability. The application of coarse grain and polarizable models is discussed. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.
Collapse
Affiliation(s)
- Ernest Awoonor-Williams
- Department of Chemistry, Memorial University of Newfoundland, St. John's, NL, A1B 3X7 Canada
| | - Christopher N Rowley
- Department of Chemistry, Memorial University of Newfoundland, St. John's, NL, A1B 3X7 Canada.
| |
Collapse
|