1
|
Chompunud Na Ayudhya C, Graidist P, Tipmanee V. Role of CSF1R 550th-tryptophan in kusunokinin and CSF1R inhibitor binding and ligand-induced structural effect. Sci Rep 2024; 14:12531. [PMID: 38822100 PMCID: PMC11143223 DOI: 10.1038/s41598-024-63505-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2024] [Accepted: 05/29/2024] [Indexed: 06/02/2024] Open
Abstract
Binding affinity is an important factor in drug design to improve drug-target selectivity and specificity. In this study, in silico techniques based on molecular docking followed by molecular dynamics (MD) simulations were utilized to identify the key residue(s) for CSF1R binding affinity among 14 pan-tyrosine kinase inhibitors and 15 CSF1R-specific inhibitors. We found tryptophan at position 550 (W550) on the CSF1R binding site interacted with the inhibitors' aromatic ring in a π-π way that made the ligands better at binding. Upon W550-Alanine substitution (W550A), the binding affinity of trans-(-)-kusunokinin and imatinib to CSF1R was significantly decreased. However, in terms of structural features, W550 did not significantly affect overall CSF1R structure, but provided destabilizing effect upon mutation. The W550A also did not either cause ligand to change its binding site or conformational changes due to ligand binding. As a result of our findings, the π-π interaction with W550's aromatic ring could be still the choice for increasing binding affinity to CSF1R. Nevertheless, our study showed that the increasing binding to W550 of the design ligand may not ensure CSF1R specificity and inhibition since W550-ligand bound state did not induce significantly conformational change into inactive state.
Collapse
Affiliation(s)
- Chompunud Chompunud Na Ayudhya
- Department of Biomedical Sciences and Biomedical Engineering, Faculty of Medicine, Prince of Songkla University, Hat Yai, 90100, Songkhla, Thailand
| | - Potchanapond Graidist
- Department of Biomedical Sciences and Biomedical Engineering, Faculty of Medicine, Prince of Songkla University, Hat Yai, 90100, Songkhla, Thailand
- Bioactivity Testing Center, Faculty of Medicine, Prince of Songkla University, Hat Yai, 90100, Songkhla, Thailand
| | - Varomyalin Tipmanee
- Department of Biomedical Sciences and Biomedical Engineering, Faculty of Medicine, Prince of Songkla University, Hat Yai, 90100, Songkhla, Thailand.
- Bioactivity Testing Center, Faculty of Medicine, Prince of Songkla University, Hat Yai, 90100, Songkhla, Thailand.
| |
Collapse
|
2
|
Balasundaram A, Ramireddy S, S UK, D TK, Tayubi IA, Zayed H, C GPD. A new horizon in the phosphorylated sites of AGA: the structural impact of C163S mutation in aspartylglucosaminuria through molecular dynamics simulation. J Biomol Struct Dyn 2024; 42:4313-4324. [PMID: 37334725 DOI: 10.1080/07391102.2023.2220798] [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: 03/18/2023] [Accepted: 05/28/2023] [Indexed: 06/20/2023]
Abstract
Aspartylglucosaminuria (AGU) is a lysosomal storage disorder caused by insufficient aspartylglucosaminidase (AGA) activity leading to chronic neurodegeneration. We utilized the PhosphoSitePlus tool to identify the AGA protein's phosphorylation sites. The phosphorylation was induced on the specific residue of the three-dimensional AGA protein, and the structural changes upon phosphorylation were studied via molecular dynamics simulation. Furthermore, the structural behaviour of C163S mutation and C163S mutation with adjacent phosphorylation was investigated. We have examined the structural impact of phosphorylated forms and C163S mutation in AGA. Molecular dynamics simulations (200 ns) exposed patterns of deviation, fluctuation, and change in compactness of Y178 phosphorylated AGA protein (Y178-p), T215 phosphorylated AGA protein (T215-p), T324 phosphorylated AGA protein (T324-p), C163S mutant AGA protein (C163S), and C163S mutation with Y178 phosphorylated AGA protein (C163S-Y178-p). Y178-p, T215-p, and C163S mutation demonstrated an increase in intramolecular hydrogen bonds, leading to greater compactness of the AGA forms. Principle component analysis (PCA) and Gibbs free energy of the phosphorylated/C163S mutation structures exhibit transition in motion/orientation than Wild type (WT). T215-p may be more dominant among these than the other studied phosphorylated forms. It might contribute to hydrolyzing L-asparagine functioning as an asparaginase, thereby regulating neurotransmitter activity. This study revealed structural insights into the phosphorylation of Y178, T215, and T324 in AGA protein. Additionally, it exposed the structural changes of the C163S mutation and C163S-Y178-p of AGA protein. This research will shed light on a better understanding of AGA's phosphorylated mechanism.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Ambritha Balasundaram
- Laboratory of Integrative Genomics, Department of Integrative Biology, School of Bio Sciences and Technology, Vellore Institute of Technology (VIT), Vellore, Tamil Nadu, India
| | - Sriroopreddy Ramireddy
- Laboratory of Integrative Genomics, Department of Integrative Biology, School of Bio Sciences and Technology, Vellore Institute of Technology (VIT), Vellore, Tamil Nadu, India
| | - Udhaya Kumar S
- Laboratory of Integrative Genomics, Department of Integrative Biology, School of Bio Sciences and Technology, Vellore Institute of Technology (VIT), Vellore, Tamil Nadu, India
| | - Thirumal Kumar D
- Faculty of Allied Health Sciences, Meenakshi Academy of Higher Education and Research, Chennai, Tamil Nadu, India
| | - Iftikhar Aslam Tayubi
- Department of Computer Science, Faculty of Computing and Information Technology, King Abdulaziz University, Jeddah, Kingdom of Saudi Arabia
| | - Hatem Zayed
- Department of Biomedical Sciences College of Health Sciences, QU Health, Qatar University, Doha, Qatar
| | - George Priya Doss C
- Laboratory of Integrative Genomics, Department of Integrative Biology, School of Bio Sciences and Technology, Vellore Institute of Technology (VIT), Vellore, Tamil Nadu, India
| |
Collapse
|
3
|
Gomatam A, Coutinho E. A chirality-sensitive approach to predict chemical transfer across the human placental barrier. Toxicol Lett 2024; 394:66-75. [PMID: 38423482 DOI: 10.1016/j.toxlet.2024.02.012] [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: 10/10/2023] [Revised: 02/24/2024] [Accepted: 02/26/2024] [Indexed: 03/02/2024]
Abstract
The placenta is a membrane that separates the fetus from the maternal circulation, and in addition to protecting the fetus, plays a key role in fetal growth and development. With increasing drug use in pregnancy, it is imperative that reliable models of estimating placental permeability and safety be established. In vitro methods and animal models such as rodent placenta are limited in application since the species-specific nature of the placenta prevents meaningful extrapolations to humans. In this regard, in silico approaches such as quantitative structure-property relationships (QSPRs) are useful alternatives. However, despite evidence that drug transport across the placenta is stereoselective (i.e., governed by the spatial arrangement of the atoms in a molecule), many QSPR models for placental transfer have been built using 2D descriptors that do not account for chirality and stereochemistry. In this study, we apply a chirality-sensitive and proven QSPR methodology titled "EigenValue ANalySis" (EVANS) to build QSPR models for placental transfer. We deploy EVANS along with robust machine learning algorithms to build (i) regression models on a dataset of environmental chemicals (dataset PD I) followed by (ii) classification models on a set of drug-like compounds (dataset PD II). The best models were found to achieve state-of-the-art performance, with the support vector machine algorithm returning rtrain2=0.85,rtest2=0.75 for PD I, and the logistic regression algorithm giving accuracy 0.88 and F1 score 0.93 for PD II. The best models were interpreted with the Shapley Additive Explanations paradigm, and it was found that autocorrelation descriptors are crucial for modelling placental permeability. In conclusion, we demonstrate the need of a chirality-sensitive approach for modelling placental transfer of chemicals, and present two predictive QSPR models that may reliably be used for prediction of placental transfer.
Collapse
Affiliation(s)
- Anish Gomatam
- Department of Pharmaceutical Chemistry, Bombay College of Pharmacy, Kalina, Mumbai, India.
| | - Evans Coutinho
- Department of Pharmaceutical Chemistry, Bombay College of Pharmacy, Kalina, Mumbai, India; St. John Institute of Pharmacy and Research, Palghar (E), Palghar, India
| |
Collapse
|
4
|
Mora Lagares L, Pérez-Castillo Y, Novič M. Exploring the dynamics of the ABCB1 membrane transporter P-glycoprotein in the presence of ATP and active/non-active compounds through molecular dynamics simulations. Toxicology 2024; 502:153732. [PMID: 38272384 DOI: 10.1016/j.tox.2024.153732] [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: 12/20/2023] [Revised: 01/15/2024] [Accepted: 01/17/2024] [Indexed: 01/27/2024]
Abstract
P-glycoprotein (Pgp) is a member of the ATP-binding cassette family of transporters that confers multidrug resistance to cancer cells and is actively involved in the pharmacokinetics and toxicokinetics of a big variety of drugs. Extensive studies have provided insights into the binding of many compounds, but the precise mechanism of translocation across the membrane remains unknown; in this context, the major challenge has been to understand the basis for its polyspecificity. In this study, molecular dynamics (MD) simulations of human P-gp (hP-gp) in an explicit membrane-and-water environment were performed to investigate the dynamic behavior of the transporter in the presence of different compounds (active and inactive) in the binding pocket and ATP molecules within the nucleotide binding domains (NBDs). The complexes studied involve four compounds: cyclosporin A (CSA), amiodarone (AMI), pamidronate (APD), and valproic acid (VPA). While CSA and AMI are known to interact with P-gp, APD and VPA do not. The results highlighted how the presence of ATP notably contributed to increased flexibility of key residues in NBD1 of active systems, indicating potential conformational changes activating the translocation mechanism. MD simulations reveal how these domains adapt and respond to the presence of different substrates, as well as the influence of ATP binding on their flexibility. Furthermore, distinctive behavior was observed in the presence of active and inactive compounds, particularly in the arrangement of ATP between NBDs, supporting the proposed nucleotide sandwich dimer mechanism for ATP binding. This study provides comprehensive insights into P-gp behavior with various ligands and ATP, offering implications for drug development, toxicity assessment and demonstrating the validity of the results derived from the MD simulations.
Collapse
Affiliation(s)
- Liadys Mora Lagares
- Theory Department, Laboratory for Cheminformatics, National Institute of Chemistry, 1000 Ljubljana, Slovenia.
| | - Yunierkis Pérez-Castillo
- Bio-Cheminformatics Research Group and Escuela de Ciencias Físicas y Matemáticas, Universidad de Las, Américas, Quito 170513, Ecuador
| | - Marjana Novič
- Theory Department, Laboratory for Cheminformatics, National Institute of Chemistry, 1000 Ljubljana, Slovenia.
| |
Collapse
|
5
|
Dixit B, Vranken W, Ghysels A. Conformational dynamics of α-1 acid glycoprotein (AGP) in cancer: A comparative study of glycosylated and unglycosylated AGP. Proteins 2024; 92:246-264. [PMID: 37837263 DOI: 10.1002/prot.26607] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2023] [Revised: 09/01/2023] [Accepted: 09/25/2023] [Indexed: 10/15/2023]
Abstract
α-1 acid glycoprotein (AGP) is one of the most abundant plasma proteins. It fulfills two important functions: immunomodulation, and binding to various drugs and receptors. These different functions are closely associated and modulated via changes in glycosylation and cancer missense mutations. From a structural point of view, glycans alter the local biophysical properties of the protein leading to a diverse ligand-binding spectrum. However, glycans can typically not be observed in the resolved X-ray crystallography structure of AGP due to their high flexibility and microheterogeneity, so limiting our understanding of AGP's conformational dynamics 70 years after its discovery. We here investigate how mutations and glycosylation interfere with AGP's conformational dynamics changing its biophysical behavior, by using molecular dynamics (MD) simulations and sequence-based dynamics predictions. The MD trajectories show that glycosylation decreases the local backbone flexibility of AGP and increases the flexibility of distant regions through allosteric effects. We observe that mutations near the glycosylation site affect glycan's conformational preferences. Thus, we conclude that mutations control glycan dynamics which modulates the protein's backbone flexibility directly affecting its accessibility. These findings may assist in the drug design targeting AGP's glycosylation and mutations in cancer.
Collapse
Affiliation(s)
- Bhawna Dixit
- IBiTech-BioMMeda Group, Ghent University, Ghent, Belgium
- Interuniversity Institute of Bioinformatics in Brussels, ULB-VUB, Brussels, Belgium
- Structural Biology Brussels, Vrije Universiteit Brussel, Brussels, Belgium
| | - Wim Vranken
- Interuniversity Institute of Bioinformatics in Brussels, ULB-VUB, Brussels, Belgium
- Structural Biology Brussels, Vrije Universiteit Brussel, Brussels, Belgium
| | - An Ghysels
- IBiTech-BioMMeda Group, Ghent University, Ghent, Belgium
| |
Collapse
|
6
|
Grandits M, Ecker GF. Ligand- and Structure-based Approaches for Transmembrane Transporter Modeling. Curr Drug Res Rev 2024; 16:81-93. [PMID: 37157206 DOI: 10.2174/2589977515666230508123041] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2022] [Revised: 03/15/2023] [Accepted: 03/28/2023] [Indexed: 05/10/2023]
Abstract
The study of transporter proteins is key to understanding the mechanism behind multidrug resistance and drug-drug interactions causing severe side effects. While ATP-binding transporters are well-studied, solute carriers illustrate an understudied family with a high number of orphan proteins. To study these transporters, in silico methods can be used to shed light on the basic molecular machinery by studying protein-ligand interactions. Nowadays, computational methods are an integral part of the drug discovery and development process. In this short review, computational approaches, such as machine learning, are discussed, which try to tackle interactions between transport proteins and certain compounds to locate target proteins. Furthermore, a few cases of selected members of the ATP binding transporter and solute carrier family are covered, which are of high interest in clinical drug interaction studies, especially for regulatory agencies. The strengths and limitations of ligand-based and structure-based methods are discussed to highlight their applicability for different studies. Furthermore, the combination of multiple approaches can improve the information obtained to find crucial amino acids that explain important interactions of protein-ligand complexes in more detail. This allows the design of drug candidates with increased activity towards a target protein, which further helps to support future synthetic efforts.
Collapse
Affiliation(s)
- Melanie Grandits
- Department of Pharmaceutical Sciences, University of Vienna, Vienna, Austria
| | - Gerhard F Ecker
- Department of Pharmaceutical Sciences, University of Vienna, Vienna, Austria
| |
Collapse
|
7
|
Jorgensen C, Ulmschneider MB, Searson PC. Modeling Substrate Entry into the P-Glycoprotein Efflux Pump at the Blood-Brain Barrier. J Med Chem 2023; 66:16615-16627. [PMID: 38097510 DOI: 10.1021/acs.jmedchem.3c01069] [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: 12/29/2023]
Abstract
We report molecular dynamics simulations of rhodamine entry into the central binding cavity of P-gp in the inward open conformation. Rhodamine can enter the inner volume via passive transport across the luminal membrane or lateral diffusion in the lipid bilayer. Entry into the inner volume is determined by the aperture angle at the apex of the protein, with a critical angle of 27° for rhodamine. The central binding cavity has an aqueous phase with a few lipids, which significantly reduces substrate diffusion. Within the central binding cavity, we identified regions with relatively weak binding, suggesting that the combination of reduced mobility and weak substrate binding confines rhodamine to enable the completion of the efflux cycle. Tariquidar, a P-gp inhibitor, aggregates at the lower arms of the P-gp, suggesting that inhibition involves steric hindrance of entry into the inner volume and/or steric hindrance of access of ATP to the nucleotide-binding domains.
Collapse
Affiliation(s)
- Christian Jorgensen
- Institute for NanoBioTechnology, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | | | - Peter C Searson
- Institute for NanoBioTechnology, Johns Hopkins University, Baltimore, Maryland 21218, United States
- Department of Materials Science and Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
8
|
Zuccoli JR, Pagnotta PA, Melito VA, Lavandera JV, Parera VE, Buzaleh AM. An Improved Technique for Genotyping the ABCB1 Gene Variant of Exon 21. Methods Protoc 2023; 6:53. [PMID: 37367997 DOI: 10.3390/mps6030053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2023] [Revised: 04/21/2023] [Accepted: 05/02/2023] [Indexed: 06/28/2023] Open
Abstract
The Multidrug Resistance protein (ABCB1, MDR1) is involved in the transport of xenobiotics and antiretroviral drugs. Some variants of the ABCB1 gene are of clinical importance; among them, exon 12 (c.1236C>T, rs1128503), 21 (c.2677G>T/A, rs2032582), and 26 (c.3435C>T, rs1045642) have a high incidence in Caucasians. Several protocols have been used for genotyping the exon 21 variants, such as allele-specific PCR-RFLP using adapted primer to generate a digestion site for several enzymes and automatic sequencing to detect the SNVs, TaqMan Allele Discrimination assay and High-Resolution Melter analysis (HRMA). The aim was to describe a new approach to genotype the three variants c.2677G>T/A for the exon 21 doing only one PCR with the corresponding primers and the digestion of the PCR product with two restriction enzymes: BrsI to identify A allele and BseYI to differentiate between G or T. An improvement of this methodology was also described. The proposal technique here described is demonstrated to be very efficient, easy, fast, reproducible, and cost-effective.
Collapse
Affiliation(s)
- Johanna Romina Zuccoli
- Centro de Investigaciones sobre Porfirinas y Porfirias (CIPYP), CONICET, Hospital de Clínicas José de San Martín, Universidad de Buenos Aires, Buenos Aires 1120, Argentina
| | - Priscila Ayelén Pagnotta
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires 2610, Argentina
- Institute of Biology and Experimental Medicine (IBYME), CONICET, Buenos Aires 2490, Argentina
| | - Viviana Alicia Melito
- Centro de Investigaciones sobre Porfirinas y Porfirias (CIPYP), CONICET, Hospital de Clínicas José de San Martín, Universidad de Buenos Aires, Buenos Aires 1120, Argentina
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires 2610, Argentina
| | - Jimena Verónica Lavandera
- Cátedra de Bromatología y Nutrición, Facultad de Bioquímica y Ciencias Biológicas, Universidad Nacional del Litoral, Santa Fe 3000, Argentina
| | - Victoria Estela Parera
- Centro de Investigaciones sobre Porfirinas y Porfirias (CIPYP), CONICET, Hospital de Clínicas José de San Martín, Universidad de Buenos Aires, Buenos Aires 1120, Argentina
| | - Ana María Buzaleh
- Centro de Investigaciones sobre Porfirinas y Porfirias (CIPYP), CONICET, Hospital de Clínicas José de San Martín, Universidad de Buenos Aires, Buenos Aires 1120, Argentina
- Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires 2610, Argentina
| |
Collapse
|
9
|
Process Modelling of Protein Crystallisation: A Case Study of Lysozyme. Chem Eng Res Des 2023. [DOI: 10.1016/j.cherd.2023.02.016] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
|
10
|
Interaction of a Homologous Series of Amphiphiles with P-glycoprotein in a Membrane Environment-Contributions of Polar and Non-Polar Interactions. Pharmaceutics 2023; 15:pharmaceutics15010174. [PMID: 36678803 PMCID: PMC9862096 DOI: 10.3390/pharmaceutics15010174] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Revised: 12/06/2022] [Accepted: 12/26/2022] [Indexed: 01/05/2023] Open
Abstract
The transport of drugs by efflux transporters in biomembranes limits their bioavailability and is a major determinant of drug resistance development by cancer cells and pathogens. A large number of chemically dissimilar drugs are transported, and despite extensive studies, the molecular determinants of substrate specificity are still not well understood. In this work, we explore the role of polar and non-polar interactions on the interaction of a homologous series of fluorescent amphiphiles with the efflux transporter P-glycoprotein. The interaction of the amphiphiles with P-glycoprotein is evaluated through effects on ATPase activity, efficiency in inhibition of [125I]-IAAP binding, and partition to the whole native membranes containing the transporter. The results were complemented with partition to model membranes with a representative lipid composition, and details on the interactions established were obtained from MD simulations. We show that when the total concentration of amphiphile is considered, the binding parameters obtained are apparent and do not reflect the affinity for P-gp. A new formalism is proposed that includes sequestration of the amphiphiles in the lipid bilayer and the possible binding of several molecules in P-gp's substrate-binding pocket. The intrinsic binding affinity thus obtained is essentially independent of amphiphile hydrophobicity, highlighting the importance of polar interactions. An increase in the lipophilicity and amphiphilicity led to a more efficient association with the lipid bilayer, which maintains the non-polar groups of the amphiphiles in the bilayer, while the polar groups interact with P-gp's binding pocket. The presence of several amphiphiles in this orientation is proposed as a mechanism for inhibition of P-pg function.
Collapse
|
11
|
Williams A, Zhan CG. Fast Prediction of Binding Affinities of SARS-CoV-2 Spike Protein and Its Mutants with Antibodies through Intermolecular Interaction Modeling-Based Machine Learning. J Phys Chem B 2022; 126:5194-5206. [PMID: 35817617 PMCID: PMC9301770 DOI: 10.1021/acs.jpcb.2c02123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Revised: 06/28/2022] [Indexed: 11/30/2022]
Abstract
Since the introduction of the novel SARS-CoV-2 virus (COVID-19) in late 2019, various new variants have appeared with mutations that confer resistance to the vaccines and monoclonal antibodies that were developed in response to the wild-type virus. As we continue through the pandemic, an accurate and efficient methodology is needed to help predict the effects certain mutations will have on both our currently produced therapeutics and those that are in development. Using published cryo-electron microscopy and X-ray crystallography structures of the spike receptor binding domain region with currently known antibodies, in the present study, we created and cross-validated an intermolecular interaction modeling-based multi-layer perceptron machine learning approach that can accurately predict the mutation-caused shifts in the binding affinity between the spike protein (wild-type or mutant) and various antibodies. This validated artificial intelligence (AI) model was used to predict the binding affinity (Kd) of reported SARS-CoV-2 antibodies with various variants of concern, including the most recently identified "Deltamicron" (or "Deltacron") variant. This AI model may be employed in the future to predict the Kd of developed novel antibody therapeutics to overcome the challenging antibody resistance issue and develop structural bases for the effects of both current and new mutants of the spike protein. In addition, the similar AI strategy and approach based on modeling of the intermolecular interactions may be useful in development of machine learning models predicting binding affinities for other protein-protein binding systems, including other antibodies binding with their antigens.
Collapse
Affiliation(s)
- Alexander
H. Williams
- Molecular
Modeling and Biopharmaceutical Center, University
of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
- Department
of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
| | - Chang-Guo Zhan
- Molecular
Modeling and Biopharmaceutical Center, University
of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
- Department
of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
| |
Collapse
|
12
|
Williams AH, Zhan CG. Generalized Methodology for the Quick Prediction of Variant SARS-CoV-2 Spike Protein Binding Affinities with Human Angiotensin-Converting Enzyme II. J Phys Chem B 2022; 126:2353-2360. [PMID: 35315274 PMCID: PMC8982491 DOI: 10.1021/acs.jpcb.1c10718] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Revised: 02/24/2022] [Indexed: 12/25/2022]
Abstract
Variants of the SARS-CoV-2 virus continue to remain a threat 2 years from the beginning of the pandemic. As more variants arise, and the B.1.1.529 (Omicron) variant threatens to create another wave of infections, a method is needed to predict the binding affinity of the spike protein quickly and accurately with human angiotensin-converting enzyme II (ACE2). We present an accurate and convenient energy minimization/molecular mechanics Poisson-Boltzmann surface area methodology previously used with engineered ACE2 therapeutics to predict the binding affinity of the Omicron variant. Without any additional data from the variants discovered after the publication of our first model, the methodology can accurately predict the binding of the spike/ACE2 variant complexes. From this methodology, we predicted that the Omicron variant spike has a Kd of ∼22.69 nM (which is very close to the experimental Kd of 20.63 nM published during the review process of the current report) and that spike protein of the new "Stealth" Omicron variant (BA.2) will display a Kd of ∼12.9 nM with the wild-type ACE2 protein. This methodology can be used with as-yet discovered variants, allowing for quick determinations regarding the variant's infectivity versus either the wild-type virus or its variants.
Collapse
Affiliation(s)
- Alexander H. Williams
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| | - Chang-Guo Zhan
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| |
Collapse
|
13
|
|
14
|
Santos SP, Lisboa AB, Silva FS, Tiwari S, Azevedo V, Cruz ÁA, Silva ES, Pinheiro CS, Alcantara-Neves NM, Pacheco LG. Rationally designed hypoallergenic mutant variants of the house dust mite allergen Der p 21. Biochim Biophys Acta Gen Subj 2022; 1866:130096. [DOI: 10.1016/j.bbagen.2022.130096] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2021] [Revised: 12/22/2021] [Accepted: 01/17/2022] [Indexed: 11/29/2022]
|
15
|
Stephan S, Dyga M, Alabd Alhafez I, Lenhard J, Urbassek HM, Hasse H. Reproducibility of atomistic friction computer experiments: a molecular dynamics simulation study. MOLECULAR SIMULATION 2021. [DOI: 10.1080/08927022.2021.1987430] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Maximilian Dyga
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Iyad Alabd Alhafez
- Physics Department and Research Center (OPTIMAS), TU Kaiserslautern, Kaiserslautern, Germany
| | - Johannes Lenhard
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Herbert M. Urbassek
- Physics Department and Research Center (OPTIMAS), TU Kaiserslautern, Kaiserslautern, Germany
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| |
Collapse
|
16
|
Wang L, O'Mara ML. Effect of the Force Field on Molecular Dynamics Simulations of the Multidrug Efflux Protein P-Glycoprotein. J Chem Theory Comput 2021; 17:6491-6508. [PMID: 34506133 DOI: 10.1021/acs.jctc.1c00414] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Molecular dynamics (MD) simulations have been used extensively to study P-glycoprotein (P-gp), a flexible multidrug transporter that is a key player in the development of multidrug resistance to chemotherapeutics. A substantial body of literature has grown from simulation studies that have employed various simulation conditions and parameters, including AMBER, CHARMM, OPLS, GROMOS, and coarse-grained force fields, drawing conclusions from simulations spanning hundreds of nanoseconds. Each force field is typically parametrized and validated on different data and observables, usually of small molecules and peptides; there have been few comparisons of force field performance on large protein-membrane systems. Here we compare the conformational ensembles of P-gp embedded in a POPC/cholesterol bilayer generated over 500 ns of replicate simulation with five force fields from popular biomolecular families: AMBER 99SB-ILDN, CHARMM 36, OPLS-AA/L, GROMOS 54A7, and MARTINI. We find considerable differences among the ensembles with little conformational overlap, although they correspond to similar extents to structural data obtained from electron paramagnetic resonance and cross-linking studies. Moreover, each trajectory was still sampling new conformations at a high rate after 500 ns of simulation, suggesting the need for more sampling. This work highlights the need to consider known limitations of the force field used (e.g., biases toward certain secondary structures) and the simulation itself (e.g., whether sufficient sampling has been achieved) when interpreting accumulated results of simulation studies of P-gp and other transport proteins.
Collapse
Affiliation(s)
- Lily Wang
- Research School of Chemistry, College of Science, Australian National University, Canberra, ACT 2601, Australia
| | - Megan L O'Mara
- Research School of Chemistry, College of Science, Australian National University, Canberra, ACT 2601, Australia
| |
Collapse
|
17
|
Balli OI, Uversky VN, Durdagi S, Coskuner-Weber O. Challenges and limitations in the studies of glycoproteins: A computational chemist's perspective. Proteins 2021; 90:322-339. [PMID: 34549826 DOI: 10.1002/prot.26242] [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: 06/23/2021] [Revised: 08/24/2021] [Accepted: 09/07/2021] [Indexed: 11/08/2022]
Abstract
Experimenters face challenges and limitations while analyzing glycoproteins due to their high flexibility, stereochemistry, anisotropic effects, and hydration phenomena. Computational studies complement experiments and have been used in characterization of the structural properties of glycoproteins. However, recent investigations revealed that computational studies face significant challenges as well. Here, we introduce and discuss some of these challenges and weaknesses in the investigations of glycoproteins. We also present requirements of future developments in computational biochemistry and computational biology areas that could be necessary for providing more accurate structural property analyses of glycoproteins using computational tools. Further theoretical strategies that need to be and can be developed are discussed herein.
Collapse
Affiliation(s)
- Oyku Irem Balli
- Molecular Biotechnology, Turkish-German University, Istanbul, Turkey
| | - Vladimir N Uversky
- Department of Molecular Medicine and USF Health Byrd Alzheimer's Institute, Morsani College of Medicine, University of South Florida, Tampa, Florida, USA
| | - Serdar Durdagi
- Computational Biology and Molecular Simulations Laboratory, Department of Biophysics, School of Medicine, Bahcesehir University, Istanbul, Turkey
| | | |
Collapse
|
18
|
The milk-derived lactoferrin inhibits V-ATPase activity by targeting its V1 domain. Int J Biol Macromol 2021; 186:54-70. [PMID: 34237360 DOI: 10.1016/j.ijbiomac.2021.06.200] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2021] [Revised: 06/20/2021] [Accepted: 06/29/2021] [Indexed: 11/20/2022]
Abstract
Lactoferrin (Lf), a bioactive milk protein, exhibits strong anticancer and antifungal activities. The search for Lf targets and mechanisms of action is of utmost importance to enhance its effective applications. A common feature among Lf-treated cancer and fungal cells is the inhibition of a proton pump called V-ATPase. Lf-driven V-ATPase inhibition leads to cytosolic acidification, ultimately causing cell death of cancer and fungal cells. Given that a detailed elucidation of how Lf and V-ATPase interact is still missing, herein we aimed to fill this gap by employing a five-stage computational approach. Molecular dynamics simulations of both proteins were performed to obtain a robust sampling of their conformational landscape, followed by clustering, which allowed retrieving representative structures, to then perform protein-protein docking. Subsequently, molecular dynamics simulations of the docked complexes and free binding energy calculations were carried out to evaluate the dynamic binding process and build a final ranking based on the binding affinities. Detailed atomist analysis of the top ranked complexes clearly indicates that Lf binds to the V1 cytosolic domain of V-ATPase. Particularly, our data suggest that Lf binds to the interfaces between A/B subunits, where the ATP hydrolysis occurs, thus inhibiting this process. The free energy decomposition analysis further identified key binding residues that will certainly aid in the rational design of follow-up experimental studies, hence bridging computational and experimental biochemistry.
Collapse
|
19
|
Mechanistic insights into TNFR1/MADD death domains in Alzheimer's disease through conformational molecular dynamic analysis. Sci Rep 2021; 11:12256. [PMID: 34112868 PMCID: PMC8192743 DOI: 10.1038/s41598-021-91606-4] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2021] [Accepted: 05/27/2021] [Indexed: 01/22/2023] Open
Abstract
Proteins are tiny players involved in the activation and deactivation of multiple signaling cascades through interactions in cells. The TNFR1 and MADD interact with each other and mediate downstream protein signaling pathways which cause neuronal cell death and Alzheimer’s disease. In the current study, a molecular docking approach was employed to explore the interactive behavior of TNFR1 and MADD proteins and their role in the activation of downstream signaling pathways. The computational sequential and structural conformational results revealed that Asp400, Arg58, Arg59 were common residues of TNFR1 and MADD which are involved in the activation of downstream signaling pathways. Aspartic acid in negatively charged residues is involved in the biosynthesis of protein. However, arginine is a positively charged residue with the potential to interact with oppositely charged amino acids. Furthermore, our molecular dynamic simulation results also ensured the stability of the backbone of TNFR1 and MADD death domains (DDs) in binding interactions. This DDs interaction mediates some conformational changes in TNFR1 which leads to the activation of mediators proteins in the cellular signaling pathways. Taken together, a better understanding of TNFR1 and MADD receptors and their activated signaling cascade may help treat Alzheimer’s disease. The death domains of TNFR1 and MADD could be used as a novel pharmacological target for the treatment of Alzheimer’s disease by inhibiting the MAPK pathway.
Collapse
|
20
|
Williams AH, Zhan CG. Fast Prediction of Binding Affinities of the SARS-CoV-2 Spike Protein Mutant N501Y (UK Variant) with ACE2 and Miniprotein Drug Candidates. J Phys Chem B 2021; 125:4330-4336. [PMID: 33881861 PMCID: PMC8084269 DOI: 10.1021/acs.jpcb.1c00869] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2021] [Revised: 03/30/2021] [Indexed: 12/28/2022]
Abstract
A recently identified variant of SARS-CoV-2 virus, known as the United Kingdom (UK) variant (lineage B.1.1.7), has an N501Y mutation on its spike protein. SARS-CoV-2 spike protein binds with angiotensin-converting enzyme 2 (ACE2), a key protein for the viral entry into the host cells. Here, we report an efficient computational approach, including the simple energy minimizations and binding free energy calculations, starting from an experimental structure of the binding complex along with experimental calibration of the calculated binding free energies, to rapidly and reliably predict the binding affinities of the N501Y mutant with human ACE2 (hACE2) and recently reported miniprotein and hACE2 decoy (CTC-445.2) drug candidates. It has been demonstrated that the N501Y mutation markedly increases the ACE2-spike protein binding affinity (Kd) from 22 to 0.44 nM, which could partially explain why the UK variant is more infectious. The miniproteins are predicted to have ∼10,000- to 100,000-fold diminished binding affinities with the N501Y mutant, creating a need for design of novel therapeutic candidates to overcome the N501Y mutation-induced drug resistance. The N501Y mutation is also predicted to decrease the binding affinity of a hACE2 decoy (CTC-445.2) binding with the spike protein by ∼200-fold. This convenient computational approach along with experimental calibration may be similarly used in the future to predict the binding affinities of potential new variants of the spike protein.
Collapse
Affiliation(s)
- Alexander H. Williams
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| | - Chang-Guo Zhan
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| |
Collapse
|
21
|
Kapla J, Rodríguez-Espigares I, Ballante F, Selent J, Carlsson J. Can molecular dynamics simulations improve the structural accuracy and virtual screening performance of GPCR models? PLoS Comput Biol 2021; 17:e1008936. [PMID: 33983933 PMCID: PMC8186765 DOI: 10.1371/journal.pcbi.1008936] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2020] [Revised: 06/08/2021] [Accepted: 04/02/2021] [Indexed: 01/14/2023] Open
Abstract
The determination of G protein-coupled receptor (GPCR) structures at atomic resolution has improved understanding of cellular signaling and will accelerate the development of new drug candidates. However, experimental structures still remain unavailable for a majority of the GPCR family. GPCR structures and their interactions with ligands can also be modelled computationally, but such predictions have limited accuracy. In this work, we explored if molecular dynamics (MD) simulations could be used to refine the accuracy of in silico models of receptor-ligand complexes that were submitted to a community-wide assessment of GPCR structure prediction (GPCR Dock). Two simulation protocols were used to refine 30 models of the D3 dopamine receptor (D3R) in complex with an antagonist. Close to 60 μs of simulation time was generated and the resulting MD refined models were compared to a D3R crystal structure. In the MD simulations, the receptor models generally drifted further away from the crystal structure conformation. However, MD refinement was able to improve the accuracy of the ligand binding mode. The best refinement protocol improved agreement with the experimentally observed ligand binding mode for a majority of the models. Receptor structures with improved virtual screening performance, which was assessed by molecular docking of ligands and decoys, could also be identified among the MD refined models. Application of weak restraints to the transmembrane helixes in the MD simulations further improved predictions of the ligand binding mode and second extracellular loop. These results provide guidelines for application of MD refinement in prediction of GPCR-ligand complexes and directions for further method development.
Collapse
Affiliation(s)
- Jon Kapla
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden
| | - Ismael Rodríguez-Espigares
- Research Programme on Biomedical Informatics (GRIB), Department of Experimental and Health Sciences of Pompeu Fabra University (UPF), Hospital del Mar Medical Research Institute (IMIM), Barcelona, Spain
| | - Flavio Ballante
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden
| | - Jana Selent
- Research Programme on Biomedical Informatics (GRIB), Department of Experimental and Health Sciences of Pompeu Fabra University (UPF), Hospital del Mar Medical Research Institute (IMIM), Barcelona, Spain
| | - Jens Carlsson
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden
| |
Collapse
|
22
|
Marques SM, Šupolíková L, Molčanová L, Šmejkal K, Bednar D, Slaninová I. Screening of Natural Compounds as P-Glycoprotein Inhibitors against Multidrug Resistance. Biomedicines 2021; 9:biomedicines9040357. [PMID: 33808505 PMCID: PMC8066904 DOI: 10.3390/biomedicines9040357] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2021] [Revised: 03/25/2021] [Accepted: 03/26/2021] [Indexed: 12/25/2022] Open
Abstract
Multidrug resistance (MDR) is a common problem when fighting cancer with chemotherapy. P-glycoprotein (P-gp, or MDR1) is an active pump responsible for the efflux of xenobiotics out of the cell, including anti-cancer drugs. It is a validated target against MDR. No crystal structure of the human P-gp is available to date, and only recently several cryo-EM structures have been solved. In this paper, we present a comprehensive computational approach that includes constructing the full-length three-dimensional structure of the human P-gp and its refinement using molecular dynamics. We assessed its flexibility and conformational diversity, compiling a dynamical ensemble that was used to dock a set of lignan compounds, previously reported as active P-gp inhibitors, and disclose their binding modes. Based on the statistical analysis of the docking results, we selected a system for performing the structure-based virtual screening of new potential P-gp inhibitors. We tested the method on a library of 87 natural flavonoids described in the literature, and 10 of those were experimentally assayed. The results reproduced the theoretical predictions only partially due to various possible factors. However, at least two of the predicted natural flavonoids were demonstrated to be effective P-gp inhibitors. They were able to increase the accumulation of doxorubicin inside the human promyelocytic leukemia HL60/MDR cells overexpressing P-gp and potentiate the antiproliferative activity of this anti-cancer drug.
Collapse
Affiliation(s)
- Sérgio M. Marques
- Loschmidt Laboratories, Department of Experimental Biology and Research Centre for Toxic Compounds in the Environment RECETOX, Faculty of Science, Masaryk University, Kamenice 5/C13, 625-00 Brno, Czech Republic;
- International Clinical Research Center, St. Anne’s University Hospital Brno, Pekarska 53, 656-91 Brno, Czech Republic
| | - Lucie Šupolíková
- Department of Biology, Faculty of Medicine, Masaryk University, Kamenice 5/A6, 625-00 Brno, Czech Republic;
| | - Lenka Molčanová
- Department of Natural Drugs, Faculty of Pharmacy, Masaryk University, Palackého 1946/1, 612-00 Brno, Czech Republic; (L.M.); (K.Š.)
| | - Karel Šmejkal
- Department of Natural Drugs, Faculty of Pharmacy, Masaryk University, Palackého 1946/1, 612-00 Brno, Czech Republic; (L.M.); (K.Š.)
| | - David Bednar
- Loschmidt Laboratories, Department of Experimental Biology and Research Centre for Toxic Compounds in the Environment RECETOX, Faculty of Science, Masaryk University, Kamenice 5/C13, 625-00 Brno, Czech Republic;
- International Clinical Research Center, St. Anne’s University Hospital Brno, Pekarska 53, 656-91 Brno, Czech Republic
- Correspondence: (D.B.); (I.S.); Tel.: +420-549492616 (D.B.); +420-549496985 (I.S.)
| | - Iva Slaninová
- Department of Biology, Faculty of Medicine, Masaryk University, Kamenice 5/A6, 625-00 Brno, Czech Republic;
- Correspondence: (D.B.); (I.S.); Tel.: +420-549492616 (D.B.); +420-549496985 (I.S.)
| |
Collapse
|
23
|
Callaghan R, Gelissen IC, George AM, Hartz AMS. Mamma Mia, P-glycoprotein binds again. FEBS Lett 2020; 594:4076-4084. [PMID: 33022784 PMCID: PMC8731231 DOI: 10.1002/1873-3468.13951] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2020] [Revised: 09/09/2020] [Accepted: 09/29/2020] [Indexed: 01/03/2023]
Abstract
The levels of amyloid peptides in the brain are regulated by a clearance pathway from neurons to the blood-brain barrier. The first step is thought to involve diffusion from the plasma membrane to the interstitium. However, amyloid peptides are hydrophobic and avidly intercalate within membranes. The ABC transporter P-glycoprotein is implicated in the clearance of amyloid peptides across the blood-brain, but its role at neurons is undetermined. We here propose that P-glycoprotein mediates 'exit' of amyloid peptides from neurons. Indeed, amyloid peptides have physicochemical similarities to substrates of P-glycoprotein, but their larger size represents a conundrum. This review probes the plausibility of a mechanism for amyloid peptide transport by P-glycoprotein exploiting evolving biochemical and structural models.
Collapse
Affiliation(s)
- Richard Callaghan
- Research School of Biology, and the Medical School, Australian National University, Canberra, ACT, Australia
| | - Ingrid C Gelissen
- School of Pharmacy, Faculty of Medicine and Health, University of Sydney, NSW, Australia
| | - Anthony M George
- School of Life Sciences, University of Technology Sydney, Broadway, NSW, Australia
| | - Anika M S Hartz
- Sanders-Brown Center on Aging, University of Kentucky, Lexington, KY, USA
| |
Collapse
|
24
|
Göddeke H, Schäfer LV. Capturing Substrate Translocation in an ABC Exporter at the Atomic Level. J Am Chem Soc 2020; 142:12791-12801. [PMID: 32578427 DOI: 10.1021/jacs.0c05502] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
ATP-binding cassette (ABC) transporters chemomechanically couple ATP binding and hydrolysis to large-scale conformational changes, ultimately leading to substrate translocation across biological membranes. Despite recent progress in the structure determination of substrate-bound ABC exporters, the inherently dynamic mechanism of substrate transport remains unclear at the atomic level. In this work, we capture substrate translocation in the heterodimeric ABC exporter TM287/288 from the hyperthermophilic bacterium Thermotoga maritima using all-atom molecular dynamics (MD) simulations. Unguided multimicrosecond simulations at 375 K show how the drugs daunorubicin and verapamil, which were initially docked into the ABC transporter, get translocated through the exporter by following its large-scale alternating access conformational transitions between an inward-facing (IF) and an outward-facing (OF) conformation. Triggered by the affinity difference due to differential solvation of the binding cavity in the IF and OF conformations, the substrates unbind from the OF transporter and partition into the lipid bilayer. While daunorubicin is stably inserted into the outer leaflet of the bilayer, verapamil dynamically flip flops between the bilayer leaflets, possibly rendering its net transport futile.
Collapse
Affiliation(s)
- Hendrik Göddeke
- Theoretical Chemistry, 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
|
25
|
Morsy MA, El-Sheikh AAK, Ibrahim ARN, Venugopala KN, Kandeel M. In silico and in vitro identification of secoisolariciresinol as a re-sensitizer of P-glycoprotein-dependent doxorubicin-resistance NCI/ADR-RES cancer cells. PeerJ 2020; 8:e9163. [PMID: 32566390 PMCID: PMC7293189 DOI: 10.7717/peerj.9163] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2020] [Accepted: 04/18/2020] [Indexed: 12/27/2022] Open
Abstract
P-glycoprotein (P-gp) is one of the highly expressed cancer cell efflux transporters that cause the failure of chemotherapy. To reverse P-gp induced multidrug resistance, we employed a flaxseed-derived lignan; secoisolariciresinol (SECO) that acts as an inhibitor of breast cancer resistance protein; another efflux transporter that shares some substrate/inhibitor specificity with P-gp. Molecular dynamics (MD) simulation identified SECO as a possible P-gp inhibitor. Comparing root mean square deviation (RMSD) of P-gp bound with SECO with that bound to its standard inhibitor verapamil showed that fluctuations in RMSD were lower in P-gp bound to SECO demonstrating higher stability of the complex of P-gp with SECO. In addition, the superimposition of P-gp structures after MD simulation showed that the nucleotide-binding domains of P-gp bound to SECO undertook a more central closer position compared with that bound to verapamil. Using rhodamine efflux assay on NCI/ADR-RES cancer cells, SECO was confirmed as a P-gp inhibitor, where cells treated with 25 or 50 µM of SECO showed significantly higher fluorescence intensity compared to control. Using MTT assay, SECO alone showed dose-dependent cytotoxicity, where 25 or 50 µM of SECO caused significantly less NCI/ADR-RES cellular viability compared to control. Furthermore, when 50 µM of SECO was added to doxorubicin (DOX), an anticancer drug, SECO significantly enhanced DOX-induced cytotoxicity compared to DOX alone. The combination index calculated by CompuSyn software indicated synergism between DOX and SECO. Our results suggest SECO as a novel P-gp inhibitor that can re-sensitize cancer cells during DOX chemotherapy.
Collapse
Affiliation(s)
- Mohamed A Morsy
- Department of Pharmaceutical Sciences/College of Clinical Pharmacy, King Faisal University, Al-Ahsa, Eastern Region, Saudi Arabia.,Department of Pharmacology/Faculty of Medicine, Minia University, El-Minia, Egypt
| | - Azza A K El-Sheikh
- Department of Pharmacology/Faculty of Medicine, Minia University, El-Minia, Egypt.,Basic Health Sciences Department/Faculty of Medicine, Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia
| | - Ahmed R N Ibrahim
- Department of Clinical Pharmacy/College of Pharmacy, King Khalid University, Abha, Saudi Arabia.,Department of Biochemistry/Faculty of Pharmacy, Minia University, El-Minia, Egypt
| | - Katharigatta N Venugopala
- Department of Pharmaceutical Sciences/College of Clinical Pharmacy, King Faisal University, Al-Ahsa, Eastern Region, Saudi Arabia.,Department of Biotechnology and Food Technology, Durban University of Technology, Durban, South Africa
| | - Mahmoud Kandeel
- Department of Biomedical Sciences/College of Veterinary Medicine, King Faisal University, Al-Ahsa, Eastern Region, Saudi Arabia.,Department of Pharmacology/Faculty of Veterinary Medicine, Kafrelsheikh University, Kafr El-Sheikh, Egypt
| |
Collapse
|
26
|
Carey Hulyer AR, Briggs DA, O'Mara ML, Kerr ID, Harmer JR, Callaghan R. Cross-linking, DEER-spectroscopy and molecular dynamics confirm the inward facing state of P-glycoprotein in a lipid membrane. J Struct Biol 2020; 211:107513. [PMID: 32339763 DOI: 10.1016/j.jsb.2020.107513] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2019] [Revised: 04/09/2020] [Accepted: 04/23/2020] [Indexed: 02/07/2023]
Abstract
The drug efflux pump P-glycoprotein (P-gp) displays a complex transport mechanism involving multiple drug binding sites and two centres for nucleotide hydrolysis. Elucidating the molecular mechanism of transport remains elusive and the availability of P-gp structures in distinct natural and ligand trapped conformations will accelerate our understanding. The present investigation sought to provide biochemical data to validate specific features of these structures; with particular focus on the transmembrane domain that provides the transport conduit. Hence our focus was on transmembrane helices six and twelve (TM6/TM12), which are believed to participate in drug binding, as they line the central transport conduit and provide a direct link to the catalytic centres. A series of P-gp mutants were generated with a single cysteine in both TM6 and TM12 to facilitate measurement of inter-helical distances using cross-linking and DEER strategies. Experimental results were compared to published structures per se and those refined by MD simulations. This analysis revealed that the refined inward-facing murine structure (4M1M) of P-gp provides a good representation of the proximity, topography and relative motions of TM6 and TM12 in reconstituted human P-gp.
Collapse
Affiliation(s)
- Alex R Carey Hulyer
- Research School of Biology, and the Medical School, Australian National University, Canberra, ACT 2601, Australia
| | - Deborah A Briggs
- Centre for Biochemistry and Cell Biology, School of Biomedical Sciences, University of Nottingham, Queen's Medical Centre, Nottingham, UK
| | - Megan L O'Mara
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| | - Ian D Kerr
- Centre for Biochemistry and Cell Biology, School of Biomedical Sciences, University of Nottingham, Queen's Medical Centre, Nottingham, UK
| | - Jeffrey R Harmer
- The Centre for Advanced Imaging, The University of Queensland, St Lucia, QLD 4072, Australia
| | - Richard Callaghan
- Research School of Biology, and the Medical School, Australian National University, Canberra, ACT 2601, Australia.
| |
Collapse
|
27
|
Dallavalle S, Dobričić V, Lazzarato L, Gazzano E, Machuqueiro M, Pajeva I, Tsakovska I, Zidar N, Fruttero R. Improvement of conventional anti-cancer drugs as new tools against multidrug resistant tumors. Drug Resist Updat 2020; 50:100682. [PMID: 32087558 DOI: 10.1016/j.drup.2020.100682] [Citation(s) in RCA: 145] [Impact Index Per Article: 36.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2019] [Revised: 01/28/2020] [Accepted: 02/05/2020] [Indexed: 02/07/2023]
Abstract
Multidrug resistance (MDR) is the dominant cause of the failure of cancer chemotherapy. The design of antitumor drugs that are able to evade MDR is rapidly evolving, showing that this area of biomedical research attracts great interest in the scientific community. The current review explores promising recent approaches that have been developed with the aim of circumventing or overcoming MDR. Encouraging results have been obtained in the investigation of the MDR-modulating properties of various classes of natural compounds and their analogues. Inhibition of P-gp or downregulation of its expression have proven to be the main mechanisms by which MDR can be surmounted. The use of hybrid molecules that are able to simultaneously interact with two or more cancer cell targets is currently being explored as a means to circumvent drug resistance. This strategy is based on the design of hybrid compounds that are obtained either by merging the structural features of separate drugs, or by conjugating two drugs or pharmacophores via cleavable/non-cleavable linkers. The approach is highly promising due to the pharmacokinetic and pharmacodynamic advantages that can be achieved over the independent administration of the two individual components. However, it should be stressed that the task of obtaining successful multivalent drugs is a very challenging one. The conjugation of anticancer agents with nitric oxide (NO) donors has recently been developed, creating a particular class of hybrid that can combat tumor drug resistance. Appropriate NO donors have been shown to reverse drug resistance via nitration of ABC transporters and by interfering with a number of metabolic enzymes and signaling pathways. In fact, hybrid compounds that are produced by covalently attaching NO-donors and antitumor drugs have been shown to elicit a synergistic cytotoxic effect in a variety of drug resistant cancer cell lines. Another strategy to circumvent MDR is based on nanocarrier-mediated transport and the controlled release of chemotherapeutic drugs and P-gp inhibitors. Their pharmacokinetics are governed by the nanoparticle or polymer carrier and make use of the enhanced permeation and retention (EPR) effect, which can increase selective delivery to cancer cells. These systems are usually internalized by cancer cells via endocytosis and accumulate in endosomes and lysosomes, thus preventing rapid efflux. Other modalities to combat MDR are described in this review, including the pharmaco-modulation of acridine, which is a well-known scaffold in the development of bioactive compounds, the use of natural compounds as means to reverse MDR, and the conjugation of anticancer drugs with carriers that target specific tumor-cell components. Finally, the outstanding potential of in silico structure-based methods as a means to evaluate the ability of antitumor drugs to interact with drug transporters is also highlighted in this review. Structure-based design methods, which utilize 3D structural data of proteins and their complexes with ligands, are the most effective of the in silico methods available, as they provide a prediction regarding the interaction between transport proteins and their substrates and inhibitors. The recently resolved X-ray structure of human P-gp can help predict the interaction sites of designed compounds, providing insight into their binding mode and directing possible rational modifications to prevent them from becoming P-gp drug substrates. In summary, although major efforts were invested in the search for new tools to combat drug resistant tumors, they all require further implementation and methodological development. Further investigation and progress in the abovementioned strategies will provide significant advances in the rational combat against cancer MDR.
Collapse
Affiliation(s)
- Sabrina Dallavalle
- Department of Food, Environmental and Nutritional Sciences, Università degli Studi di Milano, Via Celoria 2, 20133 Milan, Italy
| | - Vladimir Dobričić
- Department of Pharmaceutical Chemistry, University of Belgrade, Faculty of Pharmacy, Vojvode Stepe 450, 11000 Belgrade, Serbia
| | - Loretta Lazzarato
- Department of Drug Science and Technology, Università degli Studi di Torino, Via Pietro Giuria 9, 10125 Turin, Italy
| | - Elena Gazzano
- Department of Oncology, Università degli Studi di Torino, Via Santena 5/bis, 10126 Turin, Italy
| | - Miguel Machuqueiro
- BioISI-Biosystems & Integrative Sciences Institute, Faculty of Sciences, University of Lisbon, C8 Building, Campo Grande, 1749-016, Lisbon, Portugal; Departamento de Química e Bioquímica, Faculdade de Ciências da Universidade de Lisboa, Campo Grande, 1749-016, Lisbon, Portugal
| | - Ilza Pajeva
- QSAR and Molecular Modelling Department, Institute of Biophysics and Biomedical Engineering, Bulgarian Academy of Sciences, Acad. G. Bonchev Str., Block 105, 1113 Sofia, Bulgaria
| | - Ivanka Tsakovska
- QSAR and Molecular Modelling Department, Institute of Biophysics and Biomedical Engineering, Bulgarian Academy of Sciences, Acad. G. Bonchev Str., Block 105, 1113 Sofia, Bulgaria
| | - Nace Zidar
- University of Ljubljana, Faculty of Pharmacy, Aškerčeva cesta 7, 1000 Ljubljana, Slovenia
| | - Roberta Fruttero
- Department of Drug Science and Technology, Università degli Studi di Torino, Via Pietro Giuria 9, 10125 Turin, Italy.
| |
Collapse
|
28
|
Wang L, Zhang L, Liu F, Sun Y. Molecular Energetics of Doxorubicin Pumping by Human P-Glycoprotein. J Chem Inf Model 2019; 59:3889-3898. [DOI: 10.1021/acs.jcim.9b00429] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Lijie Wang
- Department of Biochemical Engineering and Key Laboratory of Systems Bioengineering of the Ministry of Education, School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China
| | - Lin Zhang
- Department of Biochemical Engineering and Key Laboratory of Systems Bioengineering of the Ministry of Education, School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China
| | - Fufeng Liu
- Department of Biochemical Engineering and Key Laboratory of Systems Bioengineering of the Ministry of Education, School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China
- Key Laboratory of Industrial Fermentation Microbiology of Ministry of Education, Tianjin Key Laboratory of Industrial Microbiology, National Engineering Laboratory for Industrial Enzymes, National and Local United Engineering Lab of Metabolic Control Fermentation Technology, College of Biotechnology, Tianjin University of Science & Technology, Tianjin 300457, China
| | - Yan Sun
- Department of Biochemical Engineering and Key Laboratory of Systems Bioengineering of the Ministry of Education, School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China
| |
Collapse
|
29
|
Sakamoto M, Suzuki H, Yura K. Relationship between conformation shift and disease related variation sites in ATP-binding cassette transporter proteins. Biophys Physicobiol 2019; 16:68-79. [PMID: 30923664 PMCID: PMC6435017 DOI: 10.2142/biophysico.16.0_68] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/25/2018] [Accepted: 01/07/2019] [Indexed: 12/13/2022] Open
Abstract
Transport of small molecules across the cell membrane is a crucial biological mechanism for the maintenance of the cell activity. ABC transporter family is a huge group in the transporter membrane proteins and actively transports the substrates using the energy derived from ATP hydrolysis. In humans, there are 48 distinct genes for ABC transporters. A variation of a single amino acid in the amino acid sequence of ABC transporter has been known to be linked with certain disease. The mechanism of the onset of the disease by the variation is, however, still unclear. Recent progress in the method to measure the structures of huge membrane proteins has enabled determination of the 3D structures of ABC transporters and the accumulation of coordinate data of ABC transporter has enabled us to obtain clues for the onset of the disease caused by a single variation of amino acid residue. We compared the structures of ABC transporter in apo and ATP-binding forms and found a possible conformation shift around pivot-like residues in the transmembrane domains. When this conformation change in ABC transporter and the location of pathogenic variation were compared, we found a reasonable match between the two, explaining the onset of the disease by the variation. They likely cause impairment of the pivot-like movement, weakening of ATP binding and weakening of membrane surface interactions. These findings will give a new interpretation of the variations on ABC transporter genes and pave a way to analyse the effect of variation on protein structure and function.
Collapse
Affiliation(s)
- Mika Sakamoto
- Graduate School of Humanities and Sciences, Ochanomizu University, Bunkyo-ku, Tokyo 112-8610, Japan
| | - Hirofumi Suzuki
- School of Advanced Science and Engineering, Waseda University, Shinjuku-ku, Tokyo 169-0072, Japan
| | - Kei Yura
- Graduate School of Humanities and Sciences, Ochanomizu University, Bunkyo-ku, Tokyo 112-8610, Japan.,School of Advanced Science and Engineering, Waseda University, Shinjuku-ku, Tokyo 169-0072, Japan.,Center for Simulation Science and Informational Biology, Ochanomizu University, Bunkyo-ku, Tokyo 112-8610, Japan
| |
Collapse
|
30
|
Subramanian N, Schumann-Gillett A, Mark AE, O’Mara ML. Probing the Pharmacological Binding Sites of P-Glycoprotein Using Umbrella Sampling Simulations. J Chem Inf Model 2018; 59:2287-2298. [DOI: 10.1021/acs.jcim.8b00624] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Nandhitha Subramanian
- School of Chemistry and Molecular Biosciences (SCMB), University of Queensland, Brisbane, QLD 4072, Australia
- Research School of Chemistry (RSC), Australian National University, Canberra, ACT 2601, Australia
| | | | - Alan E. Mark
- School of Chemistry and Molecular Biosciences (SCMB), University of Queensland, Brisbane, QLD 4072, Australia
- The Institute for Molecular Biosciences (IMB), University of Queensland, Brisbane, QLD 4072, Australia
| | - Megan L. O’Mara
- School of Chemistry and Molecular Biosciences (SCMB), University of Queensland, Brisbane, QLD 4072, Australia
| |
Collapse
|
31
|
Mollazadeh S, Sahebkar A, Hadizadeh F, Behravan J, Arabzadeh S. Structural and functional aspects of P-glycoprotein and its inhibitors. Life Sci 2018; 214:118-123. [PMID: 30449449 DOI: 10.1016/j.lfs.2018.10.048] [Citation(s) in RCA: 103] [Impact Index Per Article: 17.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2018] [Revised: 10/12/2018] [Accepted: 10/23/2018] [Indexed: 12/20/2022]
Abstract
P-glycoprotein (P-gp) is a member of ATP-binding cassette (ABC) superfamily which extrudes chemotherapeutic agents out of the cell. Suppression of this efflux activity has been the subject of numerous attempts to develop P-gp inhibitors. The aim of this review is to present up-to-date information on the structural and functional aspects of P-gp and its known inhibitors. The data presented also provide some information on drug discovery approaches for candidate P-gp inhibitors. Nucleotide-binding domains (NBDs) and drug-binding domains (DBDs) have been extensively studied to gain more information about P-gp inhibition and it looks that the ATPase activity of this pump has been the most attractive target for designing inhibitors. Hydrophobic and π-π (aromatic) interactions between P-gp binding domains and inhibitors are dominant intermolecular forces that have been reported in many studies using different methods. Many synthetic and natural products have been found to possess inhibitory or modulatory effects on drug transporter proteins. Log P value is an important factor in studying these inhibitors and has a crucial role on absorption, distribution, metabolism, and excretion (ADME) properties of candidate P-gp inhibitors.
Collapse
Affiliation(s)
- Shirin Mollazadeh
- Department of Medicinal Chemistry, School of Pharmacy, Mashhad University of Medical Sciences, Mashhad, Iran; Biotechnology Research Center, Pharmaceutical Technology Institute, Mashhad University of Medical Sciences, Mashhad, Iran
| | - Amirhossein Sahebkar
- Biotechnology Research Center, Pharmaceutical Technology Institute, Mashhad University of Medical Sciences, Mashhad, Iran; Neurogenic Inflammation Research Center, Mashhad University of Medical Sciences, Mashhad, Iran; School of Pharmacy, Mashhad University of Medical Sciences, Mashhad, Iran
| | - Farzin Hadizadeh
- Department of Medicinal Chemistry, School of Pharmacy, Mashhad University of Medical Sciences, Mashhad, Iran; Biotechnology Research Center, Pharmaceutical Technology Institute, Mashhad University of Medical Sciences, Mashhad, Iran.
| | - Javad Behravan
- Biotechnology Research Center, Pharmaceutical Technology Institute, Mashhad University of Medical Sciences, Mashhad, Iran; School of Pharmacy, Mashhad University of Medical Sciences, Mashhad, Iran.
| | - Sepideh Arabzadeh
- Biotechnology Research Center, Pharmaceutical Technology Institute, Mashhad University of Medical Sciences, Mashhad, Iran
| |
Collapse
|
32
|
Csizmadia G, Farkas B, Spagina Z, Tordai H, Hegedűs T. Quantitative comparison of ABC membrane protein type I exporter structures in a standardized way. Comput Struct Biotechnol J 2018; 16:396-403. [PMID: 30425800 PMCID: PMC6222291 DOI: 10.1016/j.csbj.2018.10.008] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2018] [Revised: 10/12/2018] [Accepted: 10/13/2018] [Indexed: 12/24/2022] Open
Abstract
An increasing number of ABC membrane protein structures are determined by cryo-electron microscopy and X-ray crystallography, consequently identifying differences between their conformations has become an arising issue. Therefore, we propose to define standardized measures for ABC Type I exporter structure characterization. We set conformational vectors, conftors, which describe the relative orientation of domains and can highlight structural differences. In addition, continuum electrostatics calculations were performed to characterize the energetics of membrane insertion illuminating functionally crucial regions. In summary, the proposed metrics contribute to deeper understanding of ABC membrane proteins' structural features, structure validation, and analysis of movements observed in a molecular dynamics trajectory. Moreover, the concept of standardized metrics can be applied not only to ABC membrane protein structures (http://conftors.hegelab.org).
Collapse
Key Words
- ABC proteins
- ABC, ATP binding cassette
- CFTR, cystic fibrosis transmembrane conductance regulator
- CG, coarse grained
- CH, coupling helix
- COG, center of geometry
- ICD, intracellular domain
- Membrane proteins
- NBD, nucleotide binding domain
- Quantitative structural properties
- Structure comparison
- Structure validation
- TH, transmembrane helix
- TM, transmembrane
- TMD, transmembrane domain
Collapse
Affiliation(s)
- Georgina Csizmadia
- MTA-SE Molecular Biophysics Research Group, Hungarian Academy of Sciences, Budapest, Hungary.,Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
| | - Bianka Farkas
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary.,Faculty of Information Technology, Pázmány Péter Catholic University, Budapest, Hungary
| | - Zoltán Spagina
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary.,Faculty of Information Technology, Pázmány Péter Catholic University, Budapest, Hungary
| | - Hedvig Tordai
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
| | - Tamás Hegedűs
- MTA-SE Molecular Biophysics Research Group, Hungarian Academy of Sciences, Budapest, Hungary.,Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, Hungary
| |
Collapse
|
33
|
Sigdel KP, Wilt LA, Marsh BP, Roberts AG, King GM. The conformation and dynamics of P-glycoprotein in a lipid bilayer investigated by atomic force microscopy. Biochem Pharmacol 2018; 156:302-311. [PMID: 30121251 DOI: 10.1016/j.bcp.2018.08.017] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2018] [Accepted: 08/14/2018] [Indexed: 12/17/2022]
Abstract
The membrane-bound P-glycoprotein (Pgp) transporter plays a major role in human disease and drug disposition because of its ability to efflux a chemically diverse range of drugs through ATP hydrolysis and ligand-induced conformational changes. Deciphering these structural changes is key to understanding the molecular basis of transport and to developing molecules that can modulate efflux. Here, atomic force microscopy (AFM) is used to directly image individual Pgp transporter molecules in a lipid bilayer under physiological pH and ambient temperature. Analysis of the Pgp AFM images revealed "small" and "large" protrusions from the lipid bilayer with significant differences in protrusion height and volume. The geometry of these "small" and "large" protrusions correlated to the predicted extracellular (EC) and cytosolic (C) domains of the Pgp X-ray crystal structure, respectively. To assign these protrusions, simulated AFM images were produced from the Pgp X-ray crystal structures with membrane planes defined by three computational approaches, and a simulated 80 Å AFM cantilever tip. The theoretical AFM images of the EC and C domains had similar heights and volumes to the "small" and "large" protrusions in the experimental AFM images, respectively. The assignment of the protrusions in the AFM images to the EC and C domains was confirmed by changes in protrusion volume by Pgp-specific antibodies. The Pgp domains showed a considerable degree of conformational dynamics in time resolved AFM images. With this information, a model of Pgp conformational dynamics in a lipid bilayer is proposed within the context of the known Pgp X-ray crystal structures.
Collapse
Affiliation(s)
- K P Sigdel
- Department of Physics and Astronomy, University of Missouri, Columbia, MO 65211, United States
| | - L A Wilt
- Department of Pharmaceutical and Biomedical Sciences, University of Georgia, Athens, GA 30602, United States
| | - B P Marsh
- Department of Physics and Astronomy, University of Missouri, Columbia, MO 65211, United States
| | - A G Roberts
- Department of Pharmaceutical and Biomedical Sciences, University of Georgia, Athens, GA 30602, United States.
| | - G M King
- Department of Physics and Astronomy, University of Missouri, Columbia, MO 65211, United States; Joint with Department of Biochemistry, University of Missouri, Columbia, MO 65211, United States.
| |
Collapse
|