1
|
Shahpouri P, Mehralitabar H, Kheirabadi M, Kazemi Noureini S. Potential suppression of multidrug-resistance-associated protein 1 by coumarin derivatives: an insight from molecular docking and MD simulation studies. J Biomol Struct Dyn 2023:1-17. [PMID: 37667877 DOI: 10.1080/07391102.2023.2250456] [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: 04/24/2023] [Accepted: 08/15/2023] [Indexed: 09/06/2023]
Abstract
Human MRP1 protein plays a vital role in cancer multidrug resistance. Coumarins show promising pharmacological properties. Virtual screening, ADMET, molecular docking and molecular dynamics (MD) simulations were utilized as pharmacoinformatic tools to identify potential MRP1 inhibitors among coumarin derivatives. Using in silico ADMET, 50 hits were further investigated for their selectivity toward the nucleotide-binding domains (NBDs) of MRP1 using molecular docking. Accordingly, coumarin, its symmetrical ketone derivative Lig. No. 4, and Reversan were candidates for focused docking study with the NBDs domains compared with ATP. The result indicates that Lig. No. 4, with the best binding score, interacts with NBDs via hydrogen bonds with residues: GLN713, LYS684, GLY683, CYS682 in NBD1, and GLY1432, GLY771, SER769 and GLN1374 in NBD2, which mostly overlap with ATP binding residues. Moreover, doxorubicin (Doxo) was docked to the transmembrane domains (TMDs) active site of MRP1. Doxo interaction with TMDs was subjected to MD simulation in the NBDs free and occupied with Lig. No. 4 states. The results showed that Doxo interacts more strongly with TMD residues in inward facing feature of TMDs helices. However, when Lig. No. 4 exists in NBDs, Doxo interactions are different, and TMD helices show more outward-facing conformation. This result may suggest a partial competitive inhibition mechanism for the Lig. No. 4 on MRP1 compared with ATP. So, it may inhibit active complex formation by interfering with ATP entrance to NBDs and locking MRP1 conformation in outward-facing mode. This study suggests a valuable coumarin derivative that can be further investigated for potent MRP1 inhibitors.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Parisa Shahpouri
- Department of Biology, Faculty of Science, Hakim Sabzevari University, Sabzevar, Iran
| | - Havva Mehralitabar
- Department of Biology, Faculty of Science, Hakim Sabzevari University, Sabzevar, Iran
| | - Mitra Kheirabadi
- Department of Biology, Faculty of Science, Hakim Sabzevari University, Sabzevar, Iran
| | | |
Collapse
|
2
|
Deshpande SH, Bagewadi ZK, Khan TMY, Mahnashi MH, Shaikh IA, Alshehery S, Khan AA, Patil VS, Roy S. Exploring the Potential of Phytocompounds for Targeting Epigenetic Mechanisms in Rheumatoid Arthritis: An In Silico Study Using Similarity Indexing. Molecules 2023; 28:molecules28062430. [PMID: 36985402 PMCID: PMC10051859 DOI: 10.3390/molecules28062430] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Revised: 03/02/2023] [Accepted: 03/04/2023] [Indexed: 03/30/2023] Open
Abstract
Finding structurally similar compounds in compound databases is highly efficient and is widely used in present-day drug discovery methodology. The most-trusted and -followed similarity indexing method is Tanimoto similarity indexing. Epigenetic proteins like histone deacetylases (HDACs) inhibitors are traditionally used to target cancer, but have only been investigated very recently for their possible effectiveness against rheumatoid arthritis (RA). The synthetic drugs that have been identified and used for the inhibition of HDACs include SAHA, which is being used to inhibit the activity of HDACs of different classes. SAHA was chosen as a compound of high importance as it is reported to inhibit the activity of many HDAC types. Similarity searching using the UNPD database as a reference identified aglaithioduline from the Aglaia leptantha compound as having a ~70% similarity of molecular fingerprints with SAHA, based on the Tanimoto indexing method using ChemmineR. Aglaithioduline is abundantly present in the shell and fruits of A. leptantha. In silico studies with aglaithioduline were carried out against the HDAC8 protein target and showed a binding affinity of -8.5 kcal mol. The complex was further subjected to molecular dynamics simulation using Gromacs. The RMSD, RMSF, compactness and SASA plots of the target with aglaithioduline, in comparison with the co-crystallized ligand (SAHA) system, showed a very stable configuration. The results of the study are supportive of the usage of A. leptantha and A. edulis in Indian traditional medicine for the treatment of pain-related ailments similar to RA. Our study therefore calls for further investigation of A. leptantha and A. edulis for their potential use against RA by targeting epigenetic changes, using in vivo and in vitro studies.
Collapse
Affiliation(s)
- Sanjay H Deshpande
- Department of Biotechnology, KLE Technological University, Hubballi 580031, India
| | - Zabin K Bagewadi
- Department of Biotechnology, KLE Technological University, Hubballi 580031, India
| | - T M Yunus Khan
- Department of Mechanical Engineering, College of Engineering, King Khalid University, Abha 61421, Saudi Arabia
| | - Mater H Mahnashi
- Department of Pharmaceutical Chemistry, College of Pharmacy, Najran University, Najran 66462, Saudi Arabia
| | - Ibrahim Ahmed Shaikh
- Department of Pharmacology, College of Pharmacy, Najran University, Najran 66462, Saudi Arabia
| | - Sultan Alshehery
- Department of Mechanical Engineering, College of Engineering, King Khalid University, Abha 61421, Saudi Arabia
| | - Aejaz A Khan
- Department of General Science, Ibn Sina National College for Medical Studies, Jeddah 22421, Saudi Arabia
| | - Vishal S Patil
- ICMR-National Institute of Traditional Medicine, Belagavi 590010, India
| | - Subarna Roy
- ICMR-National Institute of Traditional Medicine, Belagavi 590010, India
| |
Collapse
|
3
|
Schoeman D, Cloete R, Fielding BC. The Flexible, Extended Coil of the PDZ-Binding Motif of the Three Deadly Human Coronavirus E Proteins Plays a Role in Pathogenicity. Viruses 2022; 14:v14081707. [PMID: 36016329 PMCID: PMC9416557 DOI: 10.3390/v14081707] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2022] [Revised: 07/22/2022] [Accepted: 07/29/2022] [Indexed: 02/04/2023] Open
Abstract
The less virulent human (h) coronaviruses (CoVs) 229E, NL63, OC43, and HKU1 cause mild, self-limiting respiratory tract infections, while the more virulent SARS-CoV-1, MERS-CoV, and SARS-CoV-2 have caused severe outbreaks. The CoV envelope (E) protein, an important contributor to the pathogenesis of severe hCoV infections, may provide insight into this disparate severity of the disease. We, therefore, generated full-length E protein models for SARS-CoV-1 and -2, MERS-CoV, HCoV-229E, and HCoV-NL63 and docked C-terminal peptides of each model to the PDZ domain of the human PALS1 protein. The PDZ-binding motif (PBM) of the SARS-CoV-1 and -2 and MERS-CoV models adopted a more flexible, extended coil, while the HCoV-229E and HCoV-NL63 models adopted a less flexible alpha helix. All the E peptides docked to PALS1 occupied the same binding site and the more virulent hCoV E peptides generally interacted more stably with PALS1 than the less virulent ones. We hypothesize that the increased flexibility of the PBM in the more virulent hCoVs facilitates more stable binding to various host proteins, thereby contributing to more severe disease. This is the first paper to model full-length 3D structures for both the more virulent and less virulent hCoV E proteins, providing novel insights for possible drug and/or vaccine development.
Collapse
Affiliation(s)
- Dewald Schoeman
- Molecular Biology and Virology Research Laboratory, Department of Medical Biosciences, University of the Western Cape, Private Bag X17, Bellville, Cape Town 7535, South Africa;
| | - Ruben Cloete
- South African Medical Research Council Bioinformatics Unit, South African National Bioinformatics Institute, University of the Western Cape, Private Bag X17, Bellville, Cape Town 7535, South Africa;
| | - Burtram C. Fielding
- Molecular Biology and Virology Research Laboratory, Department of Medical Biosciences, University of the Western Cape, Private Bag X17, Bellville, Cape Town 7535, South Africa;
- Correspondence:
| |
Collapse
|
4
|
Sowlati-Hashjin S, Gandhi A, Garton M. Dawn of a New Era for Membrane Protein Design. BIODESIGN RESEARCH 2022; 2022:9791435. [PMID: 37850134 PMCID: PMC10521746 DOI: 10.34133/2022/9791435] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2022] [Accepted: 03/20/2022] [Indexed: 10/19/2023] Open
Abstract
A major advancement has recently occurred in the ability to predict protein secondary structure from sequence using artificial neural networks. This new accessibility to high-quality predicted structures provides a big opportunity for the protein design community. It is particularly welcome for membrane protein design, where the scarcity of solved structures has been a major limitation of the field for decades. Here, we review the work done to date on the membrane protein design and set out established and emerging tools that can be used to most effectively exploit this new access to structures.
Collapse
Affiliation(s)
- Shahin Sowlati-Hashjin
- Institute of Biomedical Engineering, University of Toronto, 164 College Street, Toronto, ON, Canada, M5S 3E2
| | - Aanshi Gandhi
- Institute of Biomedical Engineering, University of Toronto, 164 College Street, Toronto, ON, Canada, M5S 3E2
| | - Michael Garton
- Institute of Biomedical Engineering, University of Toronto, 164 College Street, Toronto, ON, Canada, M5S 3E2
| |
Collapse
|
5
|
Ardalan A, Sowlati-Hashjin S, Oduwoye H, Uwumarenogie SO, Karttunen M, Smith MD, Jelokhani-Niaraki M. Biphasic Proton Transport Mechanism for Uncoupling Proteins. J Phys Chem B 2021; 125:9130-9144. [PMID: 34365794 DOI: 10.1021/acs.jpcb.1c04766] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
It has been suggested that uncoupling proteins (UCPs) transport protons via interconversion between two conformational states: one in the "cytoplasmic state" and the other in the "matrix state". Matrix and cytoplasmic salt-bridge networks are key controllers of these states. This study proposes a mechanism for proton transport in tetrameric UCP2, with focus on the role of the matrix network. Eleven mutants were prepared to disrupt (K → Q or D → N mutations) or alter (K → D and D → K mutations) the salt-bridges in the matrix network. Proteins were recombinantly expressed in Escherichia coli membrane, reconstituted in model lipid membranes, and their structures and functions were analyzed by gel electrophoresis, circular dichroism spectroscopy, fluorescence assays, as well as molecular dynamics simulations. It is shown that the UCP2 matrix network contains five salt-bridges (rather than the previously reported three), and the matrix network can regulate the proton transport by holding the protein's transmembrane helices in close proximity, limiting the movement of the activator fatty acid(s). A biphasic two-state molecular model is proposed for proton transport in tetrameric (a dimer of stable dimers) UCP2, in which all the monomers are functional, and monomers in each dimer are in the same transport mode. Purine nucleotide (e.g., ATP) can occlude the internal pore of the monomeric units of UCP tetramers via interacting with positive residues at or in the proximity of the matrix network (K38, K141, K239, R88, R185, and R279) and prevent switching between cytoplasmic and matrix states, thus inhibiting the proton transport. This study provides new insights into the mechanism of proton transport and regulation in UCPs.
Collapse
Affiliation(s)
- Afshan Ardalan
- Department of Chemistry and Biochemistry, Wilfrid Laurier University, Waterloo, Ontario N2L 3C5, Canada
| | - Shahin Sowlati-Hashjin
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 3K7, Canada.,The Center for Advanced Materials and Biomaterials Research, The University of Western Ontario, London, Ontario N6K 3K7, Canada
| | - Habib Oduwoye
- Department of Chemistry and Biochemistry, Wilfrid Laurier University, Waterloo, Ontario N2L 3C5, Canada
| | - Stephanie O Uwumarenogie
- Department of Chemistry and Biochemistry, Wilfrid Laurier University, Waterloo, Ontario N2L 3C5, Canada
| | - Mikko Karttunen
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 3K7, Canada.,The Center for Advanced Materials and Biomaterials Research, The University of Western Ontario, London, Ontario N6K 3K7, Canada.,Department of Physics and Astronomy, The University of Western Ontario, London, Ontario N6A 3K7, Canada
| | - Matthew D Smith
- Department of Biology, Wilfrid Laurier University, Waterloo, Ontario N2L 3C5, Canada
| | - Masoud Jelokhani-Niaraki
- Department of Chemistry and Biochemistry, Wilfrid Laurier University, Waterloo, Ontario N2L 3C5, Canada
| |
Collapse
|
6
|
Telles IM, Bombardelli RK, dos Santos AP, Levin Y. Simulations of electroosmotic flow in charged nanopores using Dissipative Particle Dynamics with Ewald summation. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116263] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
7
|
Razavi L, Raissi H, Farzad F. Assessment of the effect of external and internal triggers on adsorption and release of paclitaxel from the PEI functionalized silicene nanosheet: A molecular dynamic simulation. J Mol Graph Model 2021; 106:107930. [PMID: 34022539 DOI: 10.1016/j.jmgm.2021.107930] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2021] [Revised: 04/14/2021] [Accepted: 04/15/2021] [Indexed: 12/24/2022]
Abstract
In order to examine the adsorption mechanisms of paclitaxel (PTX) on silicene nanosheet (SNS) molecular dynamics (MD) simulations are carried out. The MD outcomes show that the adsorption of PTX on the pristine SNS is exothermic and spontaneous. The interaction between the PTX molecule and the pristine SNS is mainly due to the formation of π-π interactions through their aromatic rings, which are supplemented by X-π (X = N-H, C-H, and CO) interactions. Upon functionalization of SNS by Polyethylenimine (PEI), drug molecules prefer to bind to the nanocarrier instead of the polymer. In the functionalized SNS (f-SNS), the binding energy of the drug with the nanocarrier becomes stronger in comparison to the SNS case (Eads: -2468.91 vs -840.95 kJ/mol). At the acidic condition, protonation of drug and PEI cause that the interaction between PTX and the nanocarrier become weaker and drug molecules could release from the nanocarrier surface. Finally, two f-SNS and protonated f-SNS (f-pSNS) systems are induced by the electric field (EF). Evaluation of the dynamic properties of these systems (with strengths 0.5 and 1 V/nm) shows that the electric field could be acted as a stimulus for drug release from nanocarriers. The obtained results from this study provide valuable information about the loading/release mechanisms of PTX on/from the SNS surface.
Collapse
Affiliation(s)
- Leila Razavi
- Department of Chemistry, University of Birjand, Birjand, Iran.
| | - Heidar Raissi
- Department of Chemistry, University of Birjand, Birjand, Iran.
| | - Farzaneh Farzad
- Department of Chemistry, University of Birjand, Birjand, Iran.
| |
Collapse
|
8
|
Ardalan A, Sowlati-Hashjin S, Uwumarenogie SO, Fish M, Mitchell J, Karttunen M, Smith MD, Jelokhani-Niaraki M. Functional Oligomeric Forms of Uncoupling Protein 2: Strong Evidence for Asymmetry in Protein and Lipid Bilayer Systems. J Phys Chem B 2020; 125:169-183. [PMID: 33373220 DOI: 10.1021/acs.jpcb.0c09422] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Stoichiometry of uncoupling proteins (UCPs) and their coexistence as functional monomeric and associated forms in lipid membranes remain intriguing open questions. In this study, tertiary and quaternary structures of UCP2 were analyzed experimentally and through molecular dynamics (MD) simulations. UCP2 was overexpressed in the inner membrane of Escherichia coli, then purified and reconstituted in lipid vesicles. Structure and proton transport function of UCP2 were characterized by circular dichroism (CD) spectroscopy and fluorescence methods. Findings suggest a tetrameric functional form for UCP2. MD simulations conclude that tetrameric UCP2 is a dimer of dimers, is more stable than its monomeric and dimeric forms, is asymmetrical and induces asymmetry in the membrane's lipid structure, and a biphasic on-off switch between the dimeric units is its possible mode of transport. MD simulations also show that the water density inside the UCP2 monomer is asymmetric, with the cytoplasmic side having a higher water density and a wider radius. In contrast, the structurally comparable adenosine 5'-diphosphate (ADP)/adenosine 5'-triphosphate (ATP) carrier (AAC1) did not form tetramers, implying that tetramerization cannot be generalized to all mitochondrial carriers.
Collapse
Affiliation(s)
- Afshan Ardalan
- Department of Chemistry and Biochemistry, Wilfrid Laurier University, Waterloo, Ontario, Canada N2L 3C5
| | - Shahin Sowlati-Hashjin
- Department of Chemistry, The University of Western Ontario, London, Ontario, Canada N6A 3K7.,Center for Advanced Materials and Biomaterials Research, The University of Western Ontario, London, Ontario, Canada N6K 3K7
| | - Stephanie O Uwumarenogie
- Department of Chemistry and Biochemistry, Wilfrid Laurier University, Waterloo, Ontario, Canada N2L 3C5
| | - Michael Fish
- Department of Chemistry and Biochemistry, Wilfrid Laurier University, Waterloo, Ontario, Canada N2L 3C5.,Department of Biology, Wilfrid Laurier University, Waterloo, Ontario, Canada N2L 3C5
| | - Joel Mitchell
- Department of Chemistry and Biochemistry, Wilfrid Laurier University, Waterloo, Ontario, Canada N2L 3C5
| | - Mikko Karttunen
- Department of Chemistry, The University of Western Ontario, London, Ontario, Canada N6A 3K7.,Center for Advanced Materials and Biomaterials Research, The University of Western Ontario, London, Ontario, Canada N6K 3K7.,Department of Applied Mathematics, The University of Western Ontario, London, Ontario, Canada N6A 5B7
| | - Matthew D Smith
- Department of Biology, Wilfrid Laurier University, Waterloo, Ontario, Canada N2L 3C5
| | - Masoud Jelokhani-Niaraki
- Department of Chemistry and Biochemistry, Wilfrid Laurier University, Waterloo, Ontario, Canada N2L 3C5
| |
Collapse
|
9
|
Bolmatov D, Soloviov D, Zhernenkov M, Zav'yalov D, Mamontov E, Suvorov A, Cai YQ, Katsaras J. Molecular Picture of the Transient Nature of Lipid Rafts. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2020; 36:4887-4896. [PMID: 32259453 DOI: 10.1021/acs.langmuir.0c00125] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
In biological membranes, lipid rafts are now thought to be transient and nanoscopic. However, the mechanism responsible for these nanoscopic assemblies remains poorly understood, even in the case of model membranes. As a result, it has proven extremely challenging to probe the physicochemical properties of lipid rafts at the molecular level. Here, we use all-atom molecular dynamics (MD) simulations and inelastic X-ray scattering (IXS), an intrinsically nanoscale technique, to directly probe the energy transfer and collective short-wavelength dynamics (phonons) of biologically relevant model membranes. We show that the nanoscale propagation of stress in lipid rafts takes place in the form of collective motions made up of longitudinal (compression waves) and transverse (shear waves) molecular vibrations. Importantly, we provide a molecular picture for the so-called van der Waals mediated "force from lipid" [Anishkin, A. et al. Proc. Natl. Acad. Sci. U.S.A. 2014, 111, 7898], a key parameter for the ionic channel mechano-transduction and the mechanism for the lipid transfer of molecular level stress [Aponte-Santamarı́a, C. et al. J. Am. Chem. Soc. 2017, 139, 13588]. Specifically, we describe how lipid rafts are formed and maintained through the propagation of molecular stress, lipid raft rattling dynamics, and a relaxation process. Eventually, the rafts dissipate through the self-diffusion of lipids making up the rafts. We also show that the molecular stress and viscoelastic properties of transient lipid rafts can be modulated through the use of hydrophobic biomolecules such as melatonin and tryptophan. Ultimately, the herein proposed mechanism describing the molecular interactions for the formation and dissolution of lipid rafts may offer insights as to how lipid rafts enable biological function.
Collapse
Affiliation(s)
- Dima Bolmatov
- Large Scale Structures Group, Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
- Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Dmytro Soloviov
- Moscow Institute of Physics and Technology, Dolgoprudny 141701, Russia
- Frank Laboratory of Neutron Physics, Joint Institute for Nuclear Research, Dubna 141980, Russia
- Institute for Safety Problems of Nuclear Power Plants, NAS of Ukraine, Kyiv 03680, Ukraine
| | - Mikhail Zhernenkov
- National Synchrotron Light Source II, Brookhaven National Laboratory, Upton, New York 11973, United States
| | | | - Eugene Mamontov
- Spectroscopy Group, Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Alexey Suvorov
- National Synchrotron Light Source II, Brookhaven National Laboratory, Upton, New York 11973, United States
| | - Yong Q Cai
- National Synchrotron Light Source II, Brookhaven National Laboratory, Upton, New York 11973, United States
| | - John Katsaras
- Large Scale Structures Group, Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
- Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996, United States
| |
Collapse
|
10
|
Molecular dynamic simulations to investigate the structural impact of known drug resistance mutations on HIV-1C Integrase-Dolutegravir binding. PLoS One 2020; 15:e0223464. [PMID: 32379830 PMCID: PMC7205217 DOI: 10.1371/journal.pone.0223464] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2019] [Accepted: 04/21/2020] [Indexed: 12/18/2022] Open
Abstract
Resistance associated mutations (RAMs) threaten the long-term success of combination antiretroviral therapy (cART) outcomes for HIV-1 treatment. HIV-1 Integrase (IN) strand transfer inhibitors (INSTIs) have proven to be a viable option for highly specific HIV-1 therapy. The INSTI, Dolutegravir is recommended by the World Health Organization for use as first-line cART. This study aims to understand how RAMs affect the stability of IN, as well as the binding of the drug Dolutegravir to the catalytic pocket of the protein. A homology model of HIV-1 subtype C IN was successfully constructed and validated. The site directed mutator webserver was used to predict destabilizing and/or stabilizing effects of known RAMs while FoldX confirmed any changes in protein energy upon introduction of mutation. Also, interaction analysis was performed between neighbouring residues. Three mutations known to be associated with Raltegravir, Elvitegravir and Dolutegravir resistance were selected; E92Q, G140S and Y143R, for molecular dynamics simulations. The structural quality assessment indicated high reliability of the HIV-1C IN tetrameric structure, with more than 90% confidence in modelled regions. Change in free energy for the three mutants indicated different effects, while simulation analysis showed G140S to have the largest affect on protein stability and flexibility. This was further supported by weaker non-bonded pairwise interaction energy and binding free energy values between the drug DTG and E92Q, Y143R and G140S mutants suggesting reduced binding affinity, as indicated by interaction analysis in comparison to the WT. Our findings suggest the G140S mutant has the strongest effect on the HIV-1C IN protein structure and Dolutegravir binding. To the best of our knowledge, this is the first study that uses the consensus wild type HIV-1C IN sequence to build an accurate 3D model to understand the effect of three known mutations on DTG drug binding in a South Africa context.
Collapse
|
11
|
Khuntawee W, Sutthibutpong T, Phongphanphanee S, Karttunen M, Wong-Ekkabut J. Molecular dynamics study of natural rubber-fullerene composites: connecting microscopic properties to macroscopic behavior. Phys Chem Chem Phys 2019; 21:19403-19413. [PMID: 31455965 DOI: 10.1039/c9cp03155c] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Macroscopic and microscopic properties of fullerene (C60)-cis-polyisoprene (cis-PI) composites at varying fullerene concentrations were investigated using atomistic molecular dynamics (MD) simulations over microsecond time scales. Results show that the introduction of fullerenes into a polymer matrix increases density, bulk modulus and heat capacity while thermal expansivity decreases. The presence of fullerenes slowed the diffusion of both C60 and cis-PI. Moreover, increasing fullerene concentration results in ordering of the cis-PI chains at the cis-PI-fullerene interfaces and shrinking of bulk PI regions. Free energy calculations of fullerene dimerization suggest that fullerenes disperse at low and aggregate at high fullerene concentrations. Our multi-scaled analysis approach demonstrates the role of 'ordered' regions adjacent to the interface between cis-PI and fullerene in controlling the level of order and mobility of the cis-PI chains. The relationship between the microscopic behavior and the changes in mechanical and thermal properties are discussed. Our study is beneficial for further studies and development of advanced rubber technology for novel, cost-effective, material with very high stiffness and thermal endurance with optimizing conditions of filler contents.
Collapse
Affiliation(s)
- Wasinee Khuntawee
- Computational Biomodelling Laboratory for Agricultural Science and Technology (CBLAST), Faculty of Science, Kasetsart University, Bangkok 10900, Thailand. and Department of Physics, Faculty of Science, Kasetsart University, Bangkok, 10900, Thailand and Thailand Center of Excellence in Physics (ThEP Center), Commission on Higher Education, Bangkok 10400, Thailand
| | - Thana Sutthibutpong
- Computational Biomodelling Laboratory for Agricultural Science and Technology (CBLAST), Faculty of Science, Kasetsart University, Bangkok 10900, Thailand. and Thailand Center of Excellence in Physics (ThEP Center), Commission on Higher Education, Bangkok 10400, Thailand and Theoretical and Computational Science Center (TaCS), Science Laboratory Building, Faculty of Science, King Mongkut's University of Technology Thonburi (KMUTT), 126 Pracha-Uthit Road, Bang Mod, Thrung Khru, Bangkok 10140, Thailand
| | - Saree Phongphanphanee
- Computational Biomodelling Laboratory for Agricultural Science and Technology (CBLAST), Faculty of Science, Kasetsart University, Bangkok 10900, Thailand. and Thailand Center of Excellence in Physics (ThEP Center), Commission on Higher Education, Bangkok 10400, Thailand and Department of Material Science, Faculty of Science, Kasetsart University, Bangkok, 10900, Thailand and Specialized Center of Rubber and Polymer Materials for Agriculture and Industry (RPM), Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
| | - Mikko Karttunen
- Department of Chemistry, The University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 5B7, Canada and Department of Applied Mathematics, The University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 5B7, Canada and The Centre of Advanced Materials and Biomaterials Research, The University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 5B7, Canada
| | - Jirasak Wong-Ekkabut
- Computational Biomodelling Laboratory for Agricultural Science and Technology (CBLAST), Faculty of Science, Kasetsart University, Bangkok 10900, Thailand. and Department of Physics, Faculty of Science, Kasetsart University, Bangkok, 10900, Thailand and Thailand Center of Excellence in Physics (ThEP Center), Commission on Higher Education, Bangkok 10400, Thailand and Specialized Center of Rubber and Polymer Materials for Agriculture and Industry (RPM), Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
| |
Collapse
|
12
|
Kongsema M, Wongkhieo S, Khongkow M, Lam EWF, Boonnoy P, Vongsangnak W, Wong-Ekkabut J. Molecular mechanism of Forkhead box M1 inhibition by thiostrepton in breast cancer cells. Oncol Rep 2019; 42:953-962. [PMID: 31322278 PMCID: PMC6667886 DOI: 10.3892/or.2019.7225] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2019] [Accepted: 07/03/2019] [Indexed: 12/16/2022] Open
Abstract
Breast cancer is the most common type of malignancies in women worldwide, and genotoxic chemotherapeutic drugs are effective by causing DNA damage in cancer cells. However, >90% of patients with metastatic cancer are resistant to chemotherapy. The Forkhead box M1 (FOXM1) transcription factor plays a pivotal role in the resistance of breast cancer cells to chemotherapy by promoting DNA damage repair following genotoxic drug treatment. The aim of the present study was to investigate the inhibition of the FOXM1 protein by thiostrepton, a natural antibiotic produced by the Streptomyces species. Experimental studies were designed to examine the effectiveness of thiostrepton in downregulating FOXM1 mRNA expression and activity, leading to senescence and apoptosis of breast cancer cells. The cytotoxicity of thiostrepton in breast cancer was determined using cell viability assay. Additionally, thiostrepton treatment decreased the mRNA expression of cyclin B1 (CCNB1), a downstream target of FOXM1. The present results indicated that thiostrepton inhibited FOXM1 mRNA expression and its effect on CCNB1. Molecular dynamic simulations were performed to study the interactions between FOXM1-DNA and thiostrepton after molecular docking. The results revealed that the possible mechanism underlying the inhibitory effect of thiostrepton on FOXM1 function was by forming a tight complex with the DNA and FOXM1 via its binding domain. Collectively, these results indicated that thiostrepton is a specific and direct inhibitor of the FOXM1 protein in breast cancer. The findings of the present study may lead to the development of novel therapeutic strategies for breast cancer and help overcome resistance to conventional chemotherapeutic drugs.
Collapse
Affiliation(s)
- Mesayamas Kongsema
- Department of Zoology, Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
| | - Sudtirak Wongkhieo
- Department of Zoology, Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
| | - Mattaka Khongkow
- National Nanotechnology Centre (NANOTEC), National Science and Technology Development Agency, Pathum Thani 12120, Thailand
| | - Eric W-F Lam
- Department of Surgery and Cancer, Imperial College London, Hammersmith Hospital Campus, London W12 0NN, UK
| | - Phansiri Boonnoy
- Computational Biomodelling Laboratory for Agricultural Science and Technology (CBLAST), Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
| | - Wanwipa Vongsangnak
- Department of Zoology, Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
| | - Jirasak Wong-Ekkabut
- Computational Biomodelling Laboratory for Agricultural Science and Technology (CBLAST), Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
| |
Collapse
|
13
|
Boonnoy P, Karttunen M, Wong-ekkabut J. Does α-Tocopherol Flip-Flop Help to Protect Membranes Against Oxidation? J Phys Chem B 2018; 122:10362-10370. [DOI: 10.1021/acs.jpcb.8b09064] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
| | | | - Jirasak Wong-ekkabut
- Thailand Center of Excellence in Physics, Commission on Higher Education, Bangkok 10400, Thailand
| |
Collapse
|
14
|
Talarimoghari M, Baaken G, Hanselmann R, Behrends JC. Size-dependent interaction of a 3-arm star poly(ethylene glycol) with two biological nanopores. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2018; 41:77. [PMID: 29926213 DOI: 10.1140/epje/i2018-11687-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/23/2018] [Accepted: 05/31/2018] [Indexed: 06/08/2023]
Abstract
We use two pore-forming proteins, alpha-hemolysin and aerolysin, to compare the polymer size-dependence of ionic current block by two types of ethyleneglycol polymers: 1) linear and 2) 3-arm star poly(ethylene glycol), both applied as a polydisperse mixture of average mass 1kDa under high salt conditions. The results demonstrate that monomer size sensitivity, as known for linear PEGs, is conserved for the star polymers with only subtle differences in the dependence of the residual conductance on monomer number. To explain this absence of a dominant effect of polymer architecture, we propose that PEG adsorbs to the inner pore wall in a collapsed, salted-out state, likely due to the effect of hydrophobic residues in the pore wall on the availability of water for hydration.
Collapse
Affiliation(s)
- Monasadat Talarimoghari
- Laboratory for Membrane Physiology and Technology, Department of Physiology, Faculty of Medicine, University of Freiburg, Hermann-Herder-Str. 7, 79104, Freiburg, Germany
| | - Gerhard Baaken
- Ionera Technologies GmbH, Hermann-Herder-Str. 7, 79104, Freiburg, Germany
| | - Ralf Hanselmann
- Institute for Macromolecular Chemistry, Stefan-Meier-Str. 31, 79104, Freiburg, Germany
- Freiburg Materials Research Centre, University of Freiburg, Stefan-Meier-Str. 19, 79104, Freiburg, Germany
| | - Jan C Behrends
- Laboratory for Membrane Physiology and Technology, Department of Physiology, Faculty of Medicine, University of Freiburg, Hermann-Herder-Str. 7, 79104, Freiburg, Germany.
- Freiburg Materials Research Centre, University of Freiburg, Stefan-Meier-Str. 19, 79104, Freiburg, Germany.
- Freiburg Centre for Interactive Materials and Bioinspired Technologies, University of Freiburg, Georges-Köhler-Allee 105, 79110, Freiburg, Germany.
| |
Collapse
|
15
|
Karttunen M, Choy WY, Cino EA. Prediction of Binding Energy of Keap1 Interaction Motifs in the Nrf2 Antioxidant Pathway and Design of Potential High-Affinity Peptides. J Phys Chem B 2018; 122:5851-5859. [PMID: 29745220 DOI: 10.1021/acs.jpcb.8b03295] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Nuclear factor erythroid 2-related factor 2 (Nrf2) is a transcription factor and principal regulator of the antioxidant pathway. The Kelch domain of Kelch-like ECH-associated protein 1 (Keap1) binds to motifs in the N-terminal region of Nrf2, promoting its degradation. There is interest in developing ligands that can compete with Nrf2 for binding to Kelch, thereby activating its transcriptional activities and increasing antioxidant levels. Using experimental Δ Gbind values of Kelch-binding motifs determined previously, a revised hydrophobicity-based model was developed for estimating Δ Gbind from amino acid sequence and applied to rank potential uncharacterized Kelch-binding motifs identified from interaction databases and BLAST searches. Model predictions and molecular dynamics (MD) simulations suggested that full-length MAD2A binds Kelch more favorably than a high-affinity 20-mer Nrf2 E78P peptide, but that the motif in isolation is not a particularly strong binder. Endeavoring to develop shorter peptides for activating Nrf2, new designs were created based on the E78P peptide, some of which showed considerable propensity to form binding-competent structures in MD, and were predicted to interact with Kelch more favorably than the E78P peptide. The peptides could be promising new ligands for enhancing the oxidative stress response.
Collapse
Affiliation(s)
- Mikko Karttunen
- Department of Chemistry and Department of Applied Mathematics , The University of Western Ontario , London , Ontario , Canada N6A 5B7
| | - Wing-Yiu Choy
- Department of Biochemistry , The University of Western Ontario , London , Ontario , Canada N6A 5C1
| | - Elio A Cino
- Department of Biochemistry and Immunology , Federal University of Minas Gerais , Belo Horizonte 31270-901 , Brazil
| |
Collapse
|
16
|
Kazemi Noureini S, Kheirabadi M, Masoumi F, Khosrogerdi F, Zarei Y, Suárez-Rozas C, Salas-Norambuena J, Kennedy Cassels B. Telomerase Inhibition by a New Synthetic Derivative of the Aporphine Alkaloid Boldine. Int J Mol Sci 2018; 19:ijms19041239. [PMID: 29671783 PMCID: PMC5979471 DOI: 10.3390/ijms19041239] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2018] [Revised: 04/15/2018] [Accepted: 04/16/2018] [Indexed: 12/16/2022] Open
Abstract
Telomerase, the enzyme responsible for cell immortality, is an important target in anti-cancer drug discovery. Boldine, an abundant aporphine alkaloid of Peumus boldus, is known to inhibit telomerase at non-toxic concentrations. Cytotoxicity of N-benzylsecoboldine hydrochloride (BSB), a synthetic derivative of boldine, was determined using the MTT method in MCF7 and MDA-MB231 cells. Aliquots of cell lysates were incubated with various concentrations of BSB in qTRAP (quantitative telomere repeat amplification protocol)-ligand experiments before substrate elongation by telomerase or amplification by hot-start Taq polymerase. The crystal structure of TERT, the catalytic subunit of telomerase from Tribolium castaneum, was used for docking and molecular dynamics analysis. The qTRAP-ligand data gave an IC50 value of about 0.17 ± 0.1 µM for BSB, roughly 400 times stronger than boldine, while the LD50 in the cytotoxicity assays were 12.5 and 21.88 µM, respectively, in cells treated for 48 h. Although both compounds interacted well with the active site, MD analysis suggests a second binding site with which BSB interacts via two hydrogen bonds, much more strongly than boldine. Theoretical analyses also evaluated the IC50 for BSB as submicromolar. BSB, with greater hydrophobicity and flexibility than boldine, represents a promising structure to inhibit telomerase at non-toxic concentrations.
Collapse
Affiliation(s)
- Sakineh Kazemi Noureini
- Department of Biology, Faculty of Sciences, Hakim Sabzevari University, Sabzevar 9617976487, Iran.
| | - Mitra Kheirabadi
- Department of Biology, Faculty of Sciences, Hakim Sabzevari University, Sabzevar 9617976487, Iran.
| | - Fatima Masoumi
- Department of Biology, Faculty of Sciences, Hakim Sabzevari University, Sabzevar 9617976487, Iran.
| | - Farve Khosrogerdi
- Department of Biology, Faculty of Sciences, Hakim Sabzevari University, Sabzevar 9617976487, Iran.
| | - Younes Zarei
- Department of Biology, Faculty of Sciences, Hakim Sabzevari University, Sabzevar 9617976487, Iran.
| | - Cristian Suárez-Rozas
- Department of Chemistry, Faculty of Sciences, University of Chile, Santiago 1058, Chile.
| | - Julio Salas-Norambuena
- Department of Chemistry, Faculty of Sciences, University of Chile, Santiago 1058, Chile.
| | - Bruce Kennedy Cassels
- Department of Chemistry, Faculty of Sciences, University of Chile, Santiago 1058, Chile.
| |
Collapse
|
17
|
Göddeke H, Timachi MH, Hutter CAJ, Galazzo L, Seeger MA, Karttunen M, Bordignon E, Schäfer LV. Atomistic Mechanism of Large-Scale Conformational Transition in a Heterodimeric ABC Exporter. J Am Chem Soc 2018; 140:4543-4551. [DOI: 10.1021/jacs.7b12944] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Hendrik Göddeke
- Theoretical Chemistry, Faculty of Chemistry and Biochemistry, Ruhr-University Bochum, 44780 Bochum, Germany
| | - M. Hadi Timachi
- EPR Spectroscopy, Faculty of Chemistry and Biochemistry, Ruhr-University Bochum, 44780 Bochum, Germany
| | - Cedric A. J. Hutter
- Institute of Medical Microbiology, University of Zürich, 8006 Zürich, Switzerland
| | - Laura Galazzo
- EPR Spectroscopy, Faculty of Chemistry and Biochemistry, Ruhr-University Bochum, 44780 Bochum, Germany
| | - Markus A. Seeger
- Institute of Medical Microbiology, University of Zürich, 8006 Zürich, Switzerland
| | - Mikko Karttunen
- Department of Chemistry and Department of Applied Mathematics, Western University, London, Ontario N6A 3K7, Canada
| | - Enrica Bordignon
- EPR Spectroscopy, Faculty of Chemistry and Biochemistry, Ruhr-University Bochum, 44780 Bochum, Germany
| | - Lars V. Schäfer
- Theoretical Chemistry, Faculty of Chemistry and Biochemistry, Ruhr-University Bochum, 44780 Bochum, Germany
| |
Collapse
|
18
|
Boonnoy P, Karttunen M, Wong-Ekkabut J. Alpha-tocopherol inhibits pore formation in oxidized bilayers. Phys Chem Chem Phys 2018; 19:5699-5704. [PMID: 28138670 DOI: 10.1039/c6cp08051k] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
In biological membranes, alpha-tocopherols (α-toc; vitamin E) protect polyunsaturated lipids from free radicals. Although the interactions of α-toc with non-oxidized lipid bilayers have been studied, their effects on oxidized bilayers remain unknown. In this study, atomistic molecular dynamics (MD) simulations of oxidized lipid bilayers were performed with varying concentrations of α-toc. Bilayers with 1-palmitoyl-2-lauroyl-sn-glycero-3-phosphocholine (PLPC) lipids and their aldehyde derivatives at a 1 : 1 ratio were studied. Our simulations show that oxidized lipids self-assemble into aggregates with a water pore rapidly developing across the bilayer. The free energy of transporting an α-toc molecule in a bilayer suggests that α-tocs can passively adsorb into it. When α-toc molecules were present at low concentrations in bilayers containing oxidized lipids, water pore formation was slowed down. At high α-toc concentrations, no pores were observed. Based on the simulations, we propose that the mechanism of how α-toc inhibits pore formation in bilayers with oxidized lipids is the following: α-tocs trap the polar groups of the oxidized lipids at the membrane-water interface resulting in a decreased probability of the oxidized lipids making contact with the two leaflets and initiating pore formation. This demonstrates that α-toc molecules not only protect the bilayer from oxidation but also help to stabilize the bilayer after lipid peroxidation occurs. These results will help in designing more efficient molecules to protect membranes from oxidative stress.
Collapse
Affiliation(s)
- Phansiri Boonnoy
- Department of Physics, Faculty of Science, Kasetsart University, Bangkok 10900, Thailand. and Computational Biomodelling Laboratory for Agricultural Science and Technology (CBLAST), Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
| | - Mikko Karttunen
- Department of Mathematics and Computer Science & Institute for Complex Molecular Systems, Eindhoven University of Technology, MetaForum, 5600 MB Eindhoven, The Netherlands and Departments of Chemistry and Applied Mathematics, Western University, 1151 Richmond Street, London, Ontario N6A 5B7, Canada
| | - Jirasak Wong-Ekkabut
- Department of Physics, Faculty of Science, Kasetsart University, Bangkok 10900, Thailand. and Computational Biomodelling Laboratory for Agricultural Science and Technology (CBLAST), Faculty of Science, Kasetsart University, Bangkok 10900, Thailand and Thailand Center of Excellence in Physics (ThEP Center), Commission on Higher Education, Bangkok 10400, Thailand
| |
Collapse
|
19
|
Khuntawee W, Karttunen M, Wong-Ekkabut J. A molecular dynamics study of conformations of beta-cyclodextrin and its eight derivatives in four different solvents. Phys Chem Chem Phys 2018; 19:24219-24229. [PMID: 28848954 DOI: 10.1039/c7cp04009a] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Understanding the atomic level interactions and the resulting structural characteristics is required for developing beta-cyclodextrin (βCD) derivatives for pharmaceutical and other applications. The effect of four different solvents on the structures of the native βCD and its hydrophilic (methylated βCD; MEβCD and hydroxypropyl βCD; HPβCD) and hydrophobic derivatives (ethylated βCD; ETβCD) was explored using molecular dynamics (MD) simulations and solvation free energy calculations. The native βCD, 2-MEβCD, 6-MEβCD, 2,6-DMβCD, 2,3,6-TMβCD, 6-HPβCD, 2,6-HPβCD and 2,6-ETβCD in non-polar solvents (cyclohexane; CHX and octane; OCT) were stably formed in a symmetric cyclic cavity shape through their intramolecular hydrogen bonds. In contrast, βCDs in polar solvents (methanol; MeOH and water; WAT) exhibited large structural changes and fluctuations leading to significant deformations of their cavities. Hydrogen bonding with polar solvents was found to be one of the major contributors to this behavior: solvent-βCD hydrogen bonding strongly competes with intramolecular bonding leading to significant changes in the structural stability of βCDs. An exception to this is the hydrophobic 2,6-ETβCD which retained its spherical cavity in all solvents. Based on this, it is proposed that the 2,6-ETβCD can act as a sustained release drug carrier.
Collapse
Affiliation(s)
- Wasinee Khuntawee
- Department of Physics, Faculty of Science, Kasetsart University, Bangkok 10900, Thailand.
| | | | | |
Collapse
|
20
|
Desikan R, Maiti PK, Ayappa KG. Assessing the Structure and Stability of Transmembrane Oligomeric Intermediates of an α-Helical Toxin. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:11496-11510. [PMID: 28930630 DOI: 10.1021/acs.langmuir.7b02277] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Protein membrane interactions play an important role in our understanding of diverse phenomena ranging from membrane-assisted protein aggregation to oligomerization and folding. Pore-forming toxins (PFTs) are the primary vehicle for infection by several strains of bacteria. These proteins which are expressed in a water-soluble form (monomers) bind to the target membrane and conformationally transform (protomers) and self-assemble to form a multimer transmembrane pore complex through a process of oligomerization. On the basis of the structure of the transmembrane domains, PFTs are broadly classified into β or α toxins. In contrast to β-PFTs, the paucity of available crystal structures coupled with the amphipathic nature of the transmembrane domains has hindered our understanding of α-PFT pore formation. In this article, we use molecular dynamics (MD) simulations to examine the process of pore formation of the bacterial α-PFT, cytolysin A from Escherichia coli (ClyA) in lipid bilayer membranes. Using atomistic MD simulations ranging from 50 to 500 ns, we show that transmembrane oligomeric intermediates or "arcs" form stable proteolipidic complexes consisting of protein arcs with toroidal lipids lining the free edges. By creating initial conditions where the lipids are contained within the arcs, we study the dynamics of spontaneous lipid evacuation and toroidal edge formation. This process occurs on the time scale of tens of nanoseconds, suggesting that once protomers oligomerize, transmembrane arcs are rapidly stabilized to form functional water channels capable of leakage. Using umbrella sampling with a coarse-grained molecular model, we obtain the free energy of insertion of a single protomer into the membrane. A single inserted protomer has a stabilization free energy of -52.9 ± 1.2 kJ/mol and forms a stable transmembrane water channel capable of leakage. Our simulations reveal that arcs are stable and viable intermediates that can occur during the pore-formation pathway for ClyA.
Collapse
Affiliation(s)
- Rajat Desikan
- Department of Chemical Engineering, ‡Centre for Condensed Matter Theory, Department of Physics, and §Centre for Biosystems Science and Engineering, Indian Institute of Science , Bengaluru, India 560012
| | - Prabal K Maiti
- Department of Chemical Engineering, ‡Centre for Condensed Matter Theory, Department of Physics, and §Centre for Biosystems Science and Engineering, Indian Institute of Science , Bengaluru, India 560012
| | - K Ganapathy Ayappa
- Department of Chemical Engineering, ‡Centre for Condensed Matter Theory, Department of Physics, and §Centre for Biosystems Science and Engineering, Indian Institute of Science , Bengaluru, India 560012
| |
Collapse
|
21
|
Temperature influence on the structure and dynamics of polymers at the interface: Atomistic molecular dynamics simulation of atactic polystyrene nanoconfined between graphene surfaces. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.08.109] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
22
|
Desikan R, Patra SM, Sarthak K, Maiti PK, Ayappa KG. Comparison of coarse-grained (MARTINI) and atomistic molecular dynamics simulations of $$\alpha $$ α and $$\beta $$ β toxin nanopores in lipid membranes. J CHEM SCI 2017. [DOI: 10.1007/s12039-017-1316-0] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
|
23
|
Sequence- and Structure-Based Functional Annotation and Assessment of Metabolic Transporters in Aspergillus oryzae: A Representative Case Study. BIOMED RESEARCH INTERNATIONAL 2016; 2016:8124636. [PMID: 27274991 PMCID: PMC4870676 DOI: 10.1155/2016/8124636] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/25/2016] [Accepted: 04/06/2016] [Indexed: 11/17/2022]
Abstract
Aspergillus oryzae is widely used for the industrial production of enzymes. In A. oryzae metabolism, transporters appear to play crucial roles in controlling the flux of molecules for energy generation, nutrients delivery, and waste elimination in the cell. While the A. oryzae genome sequence is available, transporter annotation remains limited and thus the connectivity of metabolic networks is incomplete. In this study, we developed a metabolic annotation strategy to understand the relationship between the sequence, structure, and function for annotation of A. oryzae metabolic transporters. Sequence-based analysis with manual curation showed that 58 genes of 12,096 total genes in the A. oryzae genome encoded metabolic transporters. Under consensus integrative databases, 55 unambiguous metabolic transporter genes were distributed into channels and pores (7 genes), electrochemical potential-driven transporters (33 genes), and primary active transporters (15 genes). To reveal the transporter functional role, a combination of homology modeling and molecular dynamics simulation was implemented to assess the relationship between sequence to structure and structure to function. As in the energy metabolism of A. oryzae, the H+-ATPase encoded by the AO090005000842 gene was selected as a representative case study of multilevel linkage annotation. Our developed strategy can be used for enhancing metabolic network reconstruction.
Collapse
|
24
|
The good, the bad and the user in soft matter simulations. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1858:2529-2538. [PMID: 26862882 DOI: 10.1016/j.bbamem.2016.02.004] [Citation(s) in RCA: 74] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/16/2015] [Revised: 02/01/2016] [Accepted: 02/02/2016] [Indexed: 11/21/2022]
Abstract
Molecular dynamics (MD) simulations have become popular in materials science, biochemistry, biophysics and several other fields. Improvements in computational resources, in quality of force field parameters and algorithms have yielded significant improvements in performance and reliability. On the other hand, no method of research is error free. In this review, we discuss a few examples of errors and artifacts due to various sources and discuss how to avoid them. Besides bringing attention to artifacts and proper practices in simulations, we also aim to provide the reader with a starting point to explore these issues further. In particular, we hope that the discussion encourages researchers to check software, parameters, protocols and, most importantly, their own practices in order to minimize the possibility of errors. The focus here is on practical issues. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.
Collapse
|
25
|
Wong-Ekkabut J, Karttunen M. Molecular dynamics simulation of water permeation through the alpha-hemolysin channel. J Biol Phys 2016; 42:133-46. [PMID: 26264478 PMCID: PMC4713412 DOI: 10.1007/s10867-015-9396-x] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2015] [Accepted: 07/21/2015] [Indexed: 02/05/2023] Open
Abstract
The alpha-hemolysin (AHL) nanochannel is a non-selective channel that allows for uncontrolled transport of small molecules across membranes leading to cell death. Although it is a bacterial toxin, it has promising applications, ranging from drug delivery systems to nano-sensing devices. This study focuses on the transport of water molecules through an AHL nanochannel using molecular dynamics (MD) simulations. Our results show that AHL can quickly transport water across membranes. The first-passage time approach was used to estimate the diffusion coefficient and the mean exit time. To study the energetics of transport, the potential of mean force (PMF) of a water molecule along the AHL nanochannel was calculated. The results show that the energy barriers of water permeation across a nanopore are always positive along the channel and the values are close to thermal energy (kBT). These findings suggest that the observed quick permeation of water is due to small energy barriers and a hydrophobic inner channel surface resulting in smaller friction. We speculate that these physical mechanisms are important in how AHL causes cell death.
Collapse
Affiliation(s)
- Jirasak Wong-Ekkabut
- Department of Physics, Faculty of Science, Kasetsart University, 50 Phahon Yothin Rd, Chatuchak, Bangkok, Thailand, 10900.
| | - Mikko Karttunen
- Department of Chemistry and Waterloo Institute for Nanotechnology, University of Waterloo, 200 University Avenue West, Waterloo, Ontario, Canada, N2L 3G1.
- Department of Mathematics and Computer Science & Institute for Complex Molecular Systems, Eindhoven University of Technology, PO Box 513, MetaForum, 5600 MB, Eindhoven, The Netherlands.
| |
Collapse
|
26
|
Boonnoy P, Jarerattanachat V, Karttunen M, Wong-Ekkabut J. Bilayer Deformation, Pores, and Micellation Induced by Oxidized Lipids. J Phys Chem Lett 2015; 6:4884-8. [PMID: 26673194 DOI: 10.1021/acs.jpclett.5b02405] [Citation(s) in RCA: 71] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
The influence of different oxidized lipids on lipid bilayers was investigated with 16 individual 1 μs atomistic molecular dynamics (MD) simulations. Binary mixtures of lipid bilayers of 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine (PLPC) and its peroxide and aldehyde products were performed at different concentrations. In addition, an asymmetrical short chain lipid, 1-palmitoyl-2-decanoyl-sn-glycero-3-phosphatidylcholine (PDPC), was used to compare the effects of polar/apolar groups in the lipid tail on lipid bilayer. Although water defects occurred with both aldehyde and peroxide lipids, full pore formation was observed only for aldehyde lipids. At medium concentrations the pores were stable. At higher concentrations, however, the pores became unstable and micellation occurred. Data analysis shows that aldehyde lipids' propensity for pore formation is due to their shorter and highly mobile tail. The highly polar peroxide lipids are stabilized by strong hydrogen bonds with interfacial water.
Collapse
Affiliation(s)
- Phansiri Boonnoy
- Department of Physics, Faculty of Science, Kasetsart University , Bangkok, 10900, Thailand
| | - Viwan Jarerattanachat
- Department of Physics, Faculty of Science, Kasetsart University , Bangkok, 10900, Thailand
- Clarendon Laboratory, Department of Physics, University of Oxford , Oxford OX1 3PU, United Kingdom
| | - Mikko Karttunen
- Department of Mathematics and Computer Science & Institute for Complex Molecular Systems, Eindhoven University of Technology, MetaForum , 5600 MB Eindhoven, The Netherlands
| | - Jirasak Wong-Ekkabut
- Department of Physics, Faculty of Science, Kasetsart University , Bangkok, 10900, Thailand
| |
Collapse
|
27
|
Titantah JT, Karttunen M. Hydrophobicity: effect of density and order on water's rotational slowing down. SOFT MATTER 2015; 11:7977-7985. [PMID: 26327528 DOI: 10.1039/c5sm00930h] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Ab initio molecular dynamics (AIMD) simulations of over 4.5 ns were performed in the temperature range of T = 260-350 K with van der Waals corrections to investigate the relationship between local water density and tetrahedral order in bulk water and in the presence of a hydrophobe, tetramethylurea (TMU). We demonstrate that in bulk water, defects consisting of 5- and higher coordinated water are a major contributor to dynamics. Close to a hydrophobe, 3-coordinated defects take over. The co-existence of these defects gives rise to very different local densities. We propose that the slowing down of rotational motion close to a hydrophobe is induced by an interplay between density and order with the slowing down decreasing in the following order: (i) low-density ordered-water, (ii) normal-density ordered-water, (iii) high-density ordered-water and (iv) disordered-water. The proportions of these water environments vary with temperature. These local environments also support the idea of water's polymorphism, i.e., the existence of the high- and low-density states in supercooled water.
Collapse
Affiliation(s)
- John Tatini Titantah
- Department of Applied Mathematics, University of Western Ontario, 1151 Richmond Street North, London, Ontario, Canada N6A 5B7
| | | |
Collapse
|
28
|
Khuntawee W, Wolschann P, Rungrotmongkol T, Wong-ekkabut J, Hannongbua S. Molecular Dynamics Simulations of the Interaction of Beta Cyclodextrin with a Lipid Bilayer. J Chem Inf Model 2015; 55:1894-902. [DOI: 10.1021/acs.jcim.5b00152] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
| | - Peter Wolschann
- Department
of Pharmaceutical Technology and Biopharmaceutics, University of Vienna, 14 Althan Straße 14, Vienna 1090, Austria
- Institute
of Theoretical Chemistry, University of Vienna, Währinger
Straße 17, Vienna 1090, Austria
| | | | - Jirasak Wong-ekkabut
- Department
of Physics, Faculty of Science, Kasetsart University, 50 Phahon
Yothin Road, Chatuchak, Bangkok 10900, Thailand
| | | |
Collapse
|
29
|
Nisoh N, Karttunen M, Monticelli L, Wong-ekkabut J. Lipid monolayer disruption caused by aggregated carbon nanoparticles. RSC Adv 2015. [DOI: 10.1039/c4ra17006g] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Carbon nanoparticles (CNP) have significant impact on the Pulmonary Surfactant (PS), the first biological barrier in the respiratory system.
Collapse
Affiliation(s)
- Nililla Nisoh
- Department of Physics
- Faculty of Science
- Kasetsart University
- Bangkok
- Thailand
| | - Mikko Karttunen
- Department of Chemistry and Waterloo Institute for Nanotechnology
- University of Waterloo
- Waterloo
- Canada
| | | | | |
Collapse
|
30
|
Chen W, Shen JK. Effects of system net charge and electrostatic truncation on all-atom constant pH molecular dynamics. J Comput Chem 2014; 35:1986-96. [PMID: 25142416 PMCID: PMC4165709 DOI: 10.1002/jcc.23713] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2014] [Revised: 07/20/2014] [Accepted: 08/03/2014] [Indexed: 12/21/2022]
Abstract
Constant pH molecular dynamics offers a means to rigorously study the effects of solution pH on dynamical processes. Here, we address two critical questions arising from the most recent developments of the all-atom continuous constant pH molecular dynamics (CpHMD) method: (1) What is the effect of spatial electrostatic truncation on the sampling of protonation states? (2) Is the enforcement of electrical neutrality necessary for constant pH simulations? We first examined how the generalized reaction field and force-shifting schemes modify the electrostatic forces on the titration coordinates. Free energy simulations of model compounds were then carried out to delineate the errors in the deprotonation free energy and salt-bridge stability due to electrostatic truncation and system net charge. Finally, CpHMD titration of a mini-protein HP36 was used to understand the manifestation of the two types of errors in the calculated pK(a) values. The major finding is that enforcing charge neutrality under all pH conditions and at all time via cotitrating ions significantly improves the accuracy of protonation-state sampling. We suggest that such finding is also relevant for simulations with particle mesh Ewald, considering the known artifacts due to charge-compensating background plasma.
Collapse
Affiliation(s)
- Wei Chen
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201
| | - Jana K. Shen
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201
| |
Collapse
|
31
|
Liu B, Hoopes MI, Karttunen M. Molecular Dynamics Simulations of DPPC/CTAB Monolayers at the Air/Water Interface. J Phys Chem B 2014; 118:11723-37. [PMID: 25222268 DOI: 10.1021/jp5050892] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Bin Liu
- Department of Chemistry and
Waterloo Institute for Nanotechnology, University of Waterloo, 200 University
Avenue West, Waterloo, Ontario, Canada N2L 3G1
| | - Matthew I. Hoopes
- Department of Chemistry and
Waterloo Institute for Nanotechnology, University of Waterloo, 200 University
Avenue West, Waterloo, Ontario, Canada N2L 3G1
| | - Mikko Karttunen
- Department of Chemistry and
Waterloo Institute for Nanotechnology, University of Waterloo, 200 University
Avenue West, Waterloo, Ontario, Canada N2L 3G1
| |
Collapse
|
32
|
The challenges of understanding glycolipid functions: An open outlook based on molecular simulations. Biochim Biophys Acta Mol Cell Biol Lipids 2014; 1841:1130-45. [DOI: 10.1016/j.bbalip.2013.12.016] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2013] [Revised: 12/29/2013] [Accepted: 12/30/2013] [Indexed: 11/20/2022]
|
33
|
Carter KA, Shao S, Hoopes MI, Luo D, Ahsan B, Grigoryants VM, Song W, Huang H, Zhang G, Pandey RK, Geng J, Pfeifer BA, Scholes CP, Ortega J, Karttunen M, Lovell JF. Porphyrin-phospholipid liposomes permeabilized by near-infrared light. Nat Commun 2014; 5:3546. [PMID: 24699423 PMCID: PMC3988818 DOI: 10.1038/ncomms4546] [Citation(s) in RCA: 246] [Impact Index Per Article: 24.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2013] [Accepted: 03/03/2014] [Indexed: 12/12/2022] Open
Abstract
The delivery of therapeutic compounds to target tissues is a central challenge in treating disease. Externally controlled drug release systems hold potential to selectively enhance localized delivery. Here we describe liposomes doped with porphyrin–phospholipid that are permeabilized directly by near-infrared light. Molecular dynamics simulations identified a novel light-absorbing monomer esterified from clinically approved components predicted and experimentally demonstrated to give rise to a more stable porphyrin bilayer. Light-induced membrane permeabilization is enabled with liposomal inclusion of 10 molar % porphyrin–phospholipid and occurs in the absence of bulk or nanoscale heating. Liposomes reseal following laser exposure and permeability is modulated by varying porphyrin–phospholipid doping, irradiation intensity or irradiation duration. Porphyrin–phospholipid liposomes demonstrate spatial control of release of entrapped gentamicin and temporal control of release of entrapped fluorophores following intratumoral injection. Following systemic administration, laser irradiation enhances deposition of actively loaded doxorubicin in mouse xenografts, enabling an effective single-treatment antitumour therapy. The delivery of therapeutics using an external trigger is an attractive route for the improvement of targeted disease treatment. Here, the authors have discovered a porphyrin–phospholipid liposome for light-controlled membrane permeabilization and use the system to deliver an anticancer drug in vivo.
Collapse
Affiliation(s)
- Kevin A Carter
- 1] Department of Biomedical Engineering, University at Buffalo, State University of New York, Buffalo, New York 14260, USA [2] Department of Chemical and Biological Engineering, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - Shuai Shao
- 1] Department of Biomedical Engineering, University at Buffalo, State University of New York, Buffalo, New York 14260, USA [2] Department of Chemical and Biological Engineering, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - Matthew I Hoopes
- Department of Chemistry and Waterloo Institute for Nanotechnology, University of Waterloo, Waterloo, Ontario, Canada N2L 3G1
| | - Dandan Luo
- Department of Biomedical Engineering, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - Bilal Ahsan
- Department of Biochemistry and Biomedical Sciences and M. G. DeGroote Institute for Infectious Diseases Research, McMaster University, Hamilton, Ontario, Canada L8S4L8
| | - Vladimir M Grigoryants
- Department of Chemistry, University at Albany, State University of New York, Albany, New York 12222, USA
| | - Wentao Song
- Department of Biomedical Engineering, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - Haoyuan Huang
- 1] Department of Biomedical Engineering, University at Buffalo, State University of New York, Buffalo, New York 14260, USA [2] Department of Chemical and Biological Engineering, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - Guojian Zhang
- Department of Chemical and Biological Engineering, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - Ravindra K Pandey
- PDT Center, Roswell Park Cancer Institute, Buffalo, New York 14263, USA
| | - Jumin Geng
- Department of Biomedical Engineering, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - Blaine A Pfeifer
- Department of Chemical and Biological Engineering, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - Charles P Scholes
- Department of Chemistry, University at Albany, State University of New York, Albany, New York 12222, USA
| | - Joaquin Ortega
- Department of Biochemistry and Biomedical Sciences and M. G. DeGroote Institute for Infectious Diseases Research, McMaster University, Hamilton, Ontario, Canada L8S4L8
| | - Mikko Karttunen
- Department of Chemistry and Waterloo Institute for Nanotechnology, University of Waterloo, Waterloo, Ontario, Canada N2L 3G1
| | - Jonathan F Lovell
- 1] Department of Biomedical Engineering, University at Buffalo, State University of New York, Buffalo, New York 14260, USA [2] Department of Chemical and Biological Engineering, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| |
Collapse
|
34
|
Cino EA, Killoran RC, Karttunen M, Choy WY. Binding of disordered proteins to a protein hub. Sci Rep 2014; 3:2305. [PMID: 23892546 PMCID: PMC3725505 DOI: 10.1038/srep02305] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2013] [Accepted: 07/12/2013] [Indexed: 12/30/2022] Open
Abstract
A small number of proteins, called hubs, have high connectivity and are essential for interactome functionality and integrity. Keap1 is a crucial hub in the oxidative stress response and apoptosis. The Kelch domain of Keap1 preferentially binds to disordered regions of its partners, which share similar binding motifs, but have a wide range of binding affinities. Isothermal titration calorimetry (ITC) and multi-microsecond molecular dynamics (MD) simulations were used to determine the factors that govern the affinity of all currently known disordered binding partners to Kelch. Our results show that the affinities to this hub are largely determined by the extent of preformed bound state-like conformation in the free state structures of these disordered targets. Based on our findings, we have designed a high-affinity peptide that can specifically disrupt the Keap1-NRF2 interaction and has the potential for therapeutic applications.
Collapse
Affiliation(s)
- Elio A Cino
- Department of Biochemistry, The University of Western Ontario, London, Ontario, Canada N6A 5C1
| | | | | | | |
Collapse
|
35
|
Cisneros GA, Karttunen M, Ren P, Sagui C. Classical electrostatics for biomolecular simulations. Chem Rev 2014; 114:779-814. [PMID: 23981057 PMCID: PMC3947274 DOI: 10.1021/cr300461d] [Citation(s) in RCA: 192] [Impact Index Per Article: 19.2] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
36
|
Choi Y, Attwood SJ, Hoopes MI, Drolle E, Karttunen M, Leonenko Z. Melatonin directly interacts with cholesterol and alleviates cholesterol effects in dipalmitoylphosphatidylcholine monolayers. SOFT MATTER 2014; 10:206-213. [PMID: 24651707 DOI: 10.1039/c3sm52064a] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
Melatonin is a pineal hormone that has been shown to have protective effects in several diseases that are associated with cholesterol dysregulation, including cardiovascular disease, Alzheimer's disease, and certain types of cancers. Cholesterol is a major membrane constituent with both a structural and functional influence. It is also known that melatonin readily partitions into cellular membranes. We investigated the effects of melatonin and cholesterol on the structure and physical properties of a 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) monolayer as a simple membrane model using the Langmuir-Blodgett (L-B) monolayer technique and molecular dynamics (MD) simulations. We report that melatonin increases the area per lipid and elastic compressibility of the DPPC monolayer in a concentration dependent manner, while cholesterol has the opposite effect. When both melatonin and cholesterol were present in the monolayer, the compression isotherms showed normalization of the area per molecule towards that of the pure DPPC monolayer, thus indicating that melatonin counteracts and alleviates cholesterol's effects. Atomistic MD simulations of melatonin enriched DPPC systems correlate with our experimental findings and illustrate the structural effects of both cholesterol and melatonin. Our results suggest that melatonin is able to lessen the influence of cholesterol through two different mechanisms. Firstly, we have shown that melatonin has a fluidizing effect on monolayers comprising only lipid molecules. Secondly, we also observe that melatonin interacts directly with cholesterol. Our findings suggest a direct nonspecific interaction of melatonin may be a mechanism involved in reducing cholesterol associated membrane effects, thus suggesting the existence of a new mechanism of melatonin's action. This may have important biological relevance in addition to the well-known anti-oxidative and receptor binding effects.
Collapse
Affiliation(s)
- Youngjik Choi
- Department of Biology, University of Waterloo, 200 University Avenue West, Waterloo, ON, CanadaN2L 3G1.
| | | | | | | | | | | |
Collapse
|
37
|
Cino EA, Choy WY, Karttunen M. Conformational Biases of Linear Motifs. J Phys Chem B 2013; 117:15943-57. [DOI: 10.1021/jp407536p] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Elio A. Cino
- Department
of Chemistry and Waterloo Institute for Nanotechnology, University of Waterloo, Waterloo, Ontario, Canada N2L 3G1
| | - Wing-Yiu Choy
- Department
of Biochemistry, The University of Western Ontario, London, Ontario, Canada N6A 5C1
| | - Mikko Karttunen
- Department
of Chemistry and Waterloo Institute for Nanotechnology, University of Waterloo, Waterloo, Ontario, Canada N2L 3G1
| |
Collapse
|
38
|
α-helical structures drive early stages of self-assembly of amyloidogenic amyloid polypeptide aggregate formation in membranes. Sci Rep 2013; 3:2781. [PMID: 24071712 PMCID: PMC3784961 DOI: 10.1038/srep02781] [Citation(s) in RCA: 83] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2013] [Accepted: 09/03/2013] [Indexed: 12/17/2022] Open
Abstract
The human islet amyloid polypeptide (hIAPP) is the primary component in the toxic islet amyloid deposits in type-2 diabetes. hIAPP self-assembles to aggregates that permeabilize membranes and constitutes amyloid plaques. Uncovering the mechanisms of amyloid self-assembly is the key to understanding amyloid toxicity and treatment. Although structurally similar, hIAPP's rat counterpart, the rat islet amyloid polypeptide (rIAPP), is non-toxic. It has been a puzzle why these peptides behave so differently. We combined multiscale modelling and theory to explain the drastically different dynamics of hIAPP and rIAPP: The differences stem from electrostatic dipolar interactions. hIAPP forms pentameric aggregates with the hydrophobic residues facing the membrane core and stabilizing water-conducting pores. We give predictions for pore sizes, the number of hIAPP peptides, and aggregate morphology. We show the importance of curvature-induced stress at the early stages of hIAPP assembly and the α-helical structures over β-sheets. This agrees with recent fluorescence spectroscopy experiments.
Collapse
|
39
|
Drolle E, Kučerka N, Hoopes M, Choi Y, Katsaras J, Karttunen M, Leonenko Z. Effect of melatonin and cholesterol on the structure of DOPC and DPPC membranes. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2013; 1828:2247-54. [DOI: 10.1016/j.bbamem.2013.05.015] [Citation(s) in RCA: 75] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/11/2013] [Revised: 05/17/2013] [Accepted: 05/17/2013] [Indexed: 12/14/2022]
|
40
|
Sokkar P, Choi SM, Rhee YM. Simple Method for Simulating the Mixture of Atomistic and Coarse-Grained Molecular Systems. J Chem Theory Comput 2013; 9:3728-39. [DOI: 10.1021/ct400091a] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Pandian Sokkar
- Center for
Self-assembly and Complexity, Institute for Basic Science (IBS), Pohang 790-784, Korea
| | - Sun Mi Choi
- Center for
Self-assembly and Complexity, Institute for Basic Science (IBS), Pohang 790-784, Korea
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang 790-784,
Korea
| | - Young Min Rhee
- Center for
Self-assembly and Complexity, Institute for Basic Science (IBS), Pohang 790-784, Korea
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang 790-784,
Korea
| |
Collapse
|
41
|
Jarerattanachat V, Karttunen M, Wong-ekkabut J. Molecular Dynamics Study of Oxidized Lipid Bilayers in NaCl Solution. J Phys Chem B 2013; 117:8490-501. [DOI: 10.1021/jp4040612] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Viwan Jarerattanachat
- Department of Physics, Faculty
of Science, Kasetsart University, 50 Phahon
Yothin Rd, Chatuchak, Bangkok, Thailand
| | - Mikko Karttunen
- Department of Chemistry and
Waterloo Institute for Nanotechnology, University of Waterloo, 200 University Avenue West, Waterloo, Ontario,
Canada
| | - Jirasak Wong-ekkabut
- Department of Physics, Faculty
of Science, Kasetsart University, 50 Phahon
Yothin Rd, Chatuchak, Bangkok, Thailand
| |
Collapse
|
42
|
Pourmousa M, Karttunen M. Early stages of interactions of cell-penetrating peptide penetratin with a DPPC bilayer. Chem Phys Lipids 2013; 169:85-94. [PMID: 23499547 DOI: 10.1016/j.chemphyslip.2013.02.011] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2012] [Revised: 02/28/2013] [Accepted: 02/28/2013] [Indexed: 10/27/2022]
Abstract
Cell-penetrating-peptides (CPPs) can deliver themselves together with a macromolecular cargo into cells and, hence, have promising applications in drug delivery. The detailed physical mechanisms that underlie and determine their cellular uptake remain unknown. We used molecular dynamics (MD) simulations to study the interaction of a well-known CPP, namely penetratin, with a zwitterionic di-palmitoyl-phosphatidyl-choline (DPPC) bilayer. Our study shows that the arginine and lysine residues play a crucial role in peptide-membrane binding through charge-pair and hydrogen bond interactions. We also characterize peptide conformation and show that it remains helical near the N-terminus but can fold into a variety of other conformations in the residues close to the C-terminus. The response of membrane to the peptide is also investigated.
Collapse
Affiliation(s)
- Mohsen Pourmousa
- Department of Applied Mathematics, The University of Western Ontario, London, ON, Canada N6A 5B7
| | | |
Collapse
|
43
|
Masella M, Borgis D, Cuniasse P. A multiscale coarse-grained polarizable solvent model for handling long tail bulk electrostatics. J Comput Chem 2013; 34:1112-24. [PMID: 23382002 DOI: 10.1002/jcc.23237] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2012] [Revised: 11/16/2012] [Accepted: 12/21/2012] [Indexed: 12/28/2022]
Abstract
A multiscale coarse-grained approach able to handle efficiently the solvation of microscopic solutes in extended chemical environment is described. That approach is able to compute readily and efficiently very long-range solute/solvent electrostatic microscopic interactions, up to the 1-μm scale, by considering a reduced amount of computational resources. All the required parameters are assigned to reproduce available data concerning the solvation of single ions. Such a strategy makes it possible to reproduce with good accuracy the solvation properties concerning simple ion pairs in solution (in particular, the asymptotic behavior of the ion pair potentials of mean force). This new method represents an extension of the polarizable pseudoparticle solvent model, which has been recently improved to account for the main features of hydrophobic effects in liquid water (Masella et al., J. Comput. Chem. 2011, 32, 2664). This multiscale approach is well suited to be used for computing the impact of charge changes in free energy computations, in terms of both accuracy and efficiency.
Collapse
Affiliation(s)
- Michel Masella
- Laboratoire de Chimie du Vivant, Service d'ingénierie moléculaire des protéines, Institut de biologie et de technologies de Saclay, Commissariat à l'énergie atomique, Centre de Saclay, 91191 Gif-sur-Yvette Cedex, France.
| | | | | |
Collapse
|
44
|
Pannuzzo M, Milardi D, Raudino A, Karttunen M, La Rosa C. Analytical model and multiscale simulations of Aβ peptide aggregation in lipid membranes: towards a unifying description of conformational transitions, oligomerization and membrane damage. Phys Chem Chem Phys 2013; 15:8940-51. [DOI: 10.1039/c3cp44539a] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
|
45
|
Pourmousa M, Wong-ekkabut J, Patra M, Karttunen M. Molecular dynamic studies of transportan interacting with a DPPC lipid bilayer. J Phys Chem B 2012; 117:230-41. [PMID: 23214935 DOI: 10.1021/jp310255r] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Translocation of peptides through cellular membranes is a fundamental problem in developing antimicrobial peptides and in drug delivery. There is a class of peptides, known as cell-penetrating peptides, that are able to penetrate membranes without disrupting them. They can carry pharmacological compounds, thus a promising strategy for drug delivery. The physical mechanisms that facilitate translocation are not known. We have used large-scale molecular dynamics simulations to study the penetration of transportan across a zwitterionic dipalmitoyl-phosphatidyl-choline (DPPC) bilayer. We obtained the free energy profile for one peptide inside the bilayer and discuss the response of the bilayer to the presence of transportan. We also discuss the importance of lysine residues and speculate on the possible penetration mechanism of the peptide and propose a graded-like penetration process.
Collapse
Affiliation(s)
- Mohsen Pourmousa
- Department of Applied Mathematics, The University of Western Ontario, London, Ontario, Canada N6A 5B7
| | | | | | | |
Collapse
|
46
|
Abstract
Inside cells, the concentration of macromolecules can reach up to 400 g/L. In such crowded environments, proteins are expected to behave differently than in vitro. It has been shown that the stability and the folding rate of a globular protein can be altered by the excluded volume effect produced by a high density of macromolecules. However, macromolecular crowding effects on intrinsically disordered proteins (IDPs) are less explored. These proteins can be extremely dynamic and potentially sample a wide ensemble of conformations under non-denaturing conditions. The dynamic properties of IDPs are intimately related to the timescale of conformational exchange within the ensemble, which govern target recognition and how these proteins function. In this work, we investigated the macromolecular crowding effects on the dynamics of several IDPs by measuring the NMR spin relaxation parameters of three disordered proteins (ProTα, TC1, and α-synuclein) with different extents of residual structures. To aid the interpretation of experimental results, we also performed an MD simulation of ProTα. Based on the MD analysis, a simple model to correlate the observed changes in relaxation rates to the alteration in protein motions under crowding conditions was proposed. Our results show that 1) IDPs remain at least partially disordered despite the presence of high concentration of other macromolecules, 2) the crowded environment has differential effects on the conformational propensity of distinct regions of an IDP, which may lead to selective stabilization of certain target-binding motifs, and 3) the segmental motions of IDPs on the nanosecond timescale are retained under crowded conditions. These findings strongly suggest that IDPs function as dynamic structural ensembles in cellular environments.
Collapse
|
47
|
Nguyen TD, Carrillo JMY, Dobrynin AV, Brown WM. A Case Study of Truncated Electrostatics for Simulation of Polyelectrolyte Brushes on GPU Accelerators. J Chem Theory Comput 2012; 9:73-83. [DOI: 10.1021/ct300718x] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Trung Dac Nguyen
- National Center
for Computational
Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831,
United States
| | - Jan-Michael Y. Carrillo
- Polymer Program,
Institute of
Materials Science and Department of Physics, University of Connecticut,
Storrs, Connecticut 06269, United States
| | - Andrey V. Dobrynin
- Polymer Program,
Institute of
Materials Science and Department of Physics, University of Connecticut,
Storrs, Connecticut 06269, United States
| | - W. Michael Brown
- National Center
for Computational
Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831,
United States
| |
Collapse
|
48
|
Jämbeck JPM, Lyubartsev AP. Another Piece of the Membrane Puzzle: Extending Slipids Further. J Chem Theory Comput 2012; 9:774-84. [PMID: 26589070 DOI: 10.1021/ct300777p] [Citation(s) in RCA: 210] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
To be able to model complex biological membranes in a more realistic manner, the force field Slipids (Stockholm lipids) has been extended to include parameters for sphingomyelin (SM), phosphatidylglycerol (PG), phosphatidylserine (PS) lipids, and cholesterol. Since the parametrization scheme was faithful to the scheme used in previous editions of Slipids, all parameters are consistent and fully compatible. The results of careful validation of a number of key structural properties for one and two component lipid bilayers are in excellent agreement with experiments. Potentials of mean force for transferring water across binary mixtures of lipids and cholesterol were also computed in order to compare water permeability rates to experiments. In agreement with experimental and simulation studies, it was found that the permeability and partitioning of water is affected by cholesterol in lipid bilayers made of saturated lipids to the largest extent. With the extensions of Slipids presented here, it is now possible to study complex systems containing many different lipids and proteins in a fully atomistic resolution in the isothermic-isobaric (NPT) ensemble, which is the proper ensemble for membrane simulations.
Collapse
Affiliation(s)
- Joakim P M Jämbeck
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, Stockholm, SE-10691, Sweden
| | - Alexander P Lyubartsev
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, Stockholm, SE-10691, Sweden
| |
Collapse
|