1
|
Juneja P, Rashid N, Abul Qais F, Tanwar S, Sultan I, Ahmad F, Rehman SU. Alternative splicing generates a novel ferroportin isoform with a shorter C-terminal and an intact iron- and hepcidin-binding property. IUBMB Life 2024; 76:523-533. [PMID: 38348962 DOI: 10.1002/iub.2809] [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: 12/20/2023] [Accepted: 01/09/2024] [Indexed: 08/15/2024]
Abstract
Ferroportin (FPN) is a transmembrane protein and is the only known iron exporter that helps in maintaining iron homeostasis in vertebrates. To maintain stable iron equilibrium in the body, ferroportin works in conjunction with a peptide called hepcidin. In this study, we have identified an alternatively spliced novel isoform of the human SLC40A1 gene, which encodes for the FPN protein and is found to be expressed in different tissues. The novel transcript has an alternate last exon and encodes 31-amino acid long peptide sequence that replaces 104 amino acids at C-terminal in the novel transcript. Molecular modelling and molecular dynamics (MD) simulation studies revealed key structural features of the novel isoform (FPN-N). FPN-N was predicted to have 12 transmembrane domains similar to the reported isoform (FPN), despite being much smaller in size. FPN-N was found to interact with hepcidin, a key regulator of ferroportin activity. Also, the iron-binding sites were retained in the novel isoform as revealed by the MD simulation of FPN-N in bilipid membrane. The novel isoform identified in this study may play important role in iron homeostasis. However, further studies are required to characterize the FPN-N isoform and decipher its role inside the cell.
Collapse
Affiliation(s)
- Pallavi Juneja
- Department of Biochemistry, School of Chemical and Life Sciences, New Delhi, India
| | - Naira Rashid
- Department of Biochemistry, School of Chemical and Life Sciences, New Delhi, India
| | - Faizan Abul Qais
- Department of Agricultural Microbiology, Faculty of Agricultural Sciences, Aligarh Muslim University, Aligarh, India
| | - Supriya Tanwar
- Department of Biochemistry, School of Chemical and Life Sciences, New Delhi, India
| | - Insha Sultan
- Department of Biochemistry, School of Chemical and Life Sciences, New Delhi, India
| | - Faizan Ahmad
- Department of Biochemistry, School of Chemical and Life Sciences, New Delhi, India
| | - Sayeed Ur Rehman
- Department of Biochemistry, School of Chemical and Life Sciences, New Delhi, India
| |
Collapse
|
2
|
Rout M, Dey S, Mishra S, Panda S, Singh MK, Sinha R, Dehury B, Pati S. Machine learning and classical MD simulation to identify inhibitors against the P37 envelope protein of monkeypox virus. J Biomol Struct Dyn 2024; 42:3935-3948. [PMID: 37221882 DOI: 10.1080/07391102.2023.2216290] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2023] [Accepted: 05/16/2023] [Indexed: 05/25/2023]
Abstract
Monkeypox virus (MPXV) outbreak is a serious public health concern that requires international attention. P37 of MPXV plays a pivotal role in DNA replication and acts as one of the promising targets for antiviral drug design. In this study, we intent to screen potential analogs of existing FDA approved drugs of MPXV against P37 using state-of-the-art machine learning and computational biophysical techniques. AlphaFold2 guided all-atoms molecular dynamics simulations optimized P37 structure is used for molecular docking and binding free energy calculations. Similar to members of Phospholipase-D family , the predicted P37 structure also adopts a β-α-β-α-β sandwich fold, harbouring strongly conserved HxKxxxxD motif. The binding pocket comprises of Tyr48, Lys86, His115, Lys117, Ser130, Asn132, Trp280, Asn240, His325, Lys327 and Tyr346 forming strong hydrogen bonds and dense hydrophobic contacts with the screened analogs and is surrounded by positively charged patches. Loops connecting the two domains and C-terminal region exhibit high degree of flexibility. In some structural ensembles, the partial disorderness in the C-terminal region is presumed to be due to its low confidence score, acquired during structure prediction. Transition from loop to β-strands (244-254 aa) in P37-Cidofovir and its analog complexes advocates the need for further investigations. MD simulations support the accuracy of the molecular docking results, indicating the potential of analogs as potent binders of P37. Taken together, our results provide preferable understanding of molecular recognition and dynamics of ligand-bound states of P37, offering opportunities for development of new antivirals against MPXV. However, the need of in vitro and in vivo assays for confirmation of these results still persists.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Madhusmita Rout
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Nalco Square, Bhubaneswar, Odisha, India
| | - Suchanda Dey
- Biomics and Biodiversity Lab, Siksha 'O' Anusandhan (deemed to be) University, Bhubaneswar, Odisha, India
| | - Sarbani Mishra
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Nalco Square, Bhubaneswar, Odisha, India
| | - Sunita Panda
- Mycology Division, ICMR-Regional Medical Research Centre, Nalco Square, Bhubaneswar, Odisha, India
| | - Mahender Kumar Singh
- Data Science Laboratory, National Brain Research Centre, Gurgaon, Haryana, India
| | - Rohan Sinha
- Computer Science, National Institute of Technology Patna, Patna, India
| | - Budheswar Dehury
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Nalco Square, Bhubaneswar, Odisha, India
| | - Sanghamitra Pati
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Nalco Square, Bhubaneswar, Odisha, India
| |
Collapse
|
3
|
Dey S, Rout M, Pati S, Singh MK, Dehury B, Subudhi E. All-atoms molecular dynamics study to screen potent efflux pump inhibitors against KpnE protein of Klebsiella pneumoniae. J Biomol Struct Dyn 2024; 42:3492-3506. [PMID: 37218086 DOI: 10.1080/07391102.2023.2214232] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2023] [Accepted: 05/07/2023] [Indexed: 05/24/2023]
Abstract
The Small Multidrug Resistance efflux pump protein KpnE, plays a pivotal role in multi-drug resistance in Klebsiella pneumoniae. Despite well-documented study of its close homolog, EmrE, from Escherichia coli, the mechanism of drug binding to KpnE remains obscure due to the absence of a high-resolution experimental structure. Herein, we exclusively elucidate its structure-function mechanism and report some of the potent inhibitors through drug repurposing. We used molecular dynamics simulation to develop a dimeric structure of KpnE and explore its dynamics in lipid-mimetic bilayers. Our study identified both semi-open and open conformations of KpnE, highlighting its importance in transport process. Electrostatic surface potential map suggests a considerable degree of similarity between KpnE and EmrE at the binding cleft, mostly occupied by negatively charged residues. We identify key amino acids Glu14, Trp63 and Tyr44, indispensable for ligand recognition. Molecular docking and binding free energy calculations recognizes potential inhibitors like acarbose, rutin and labetalol. Further validations are needed to confirm the therapeutic role of these compounds. Altogether, our membrane dynamics study uncovers the crucial charged patches, lipid-binding sites and flexible loop that could potentiate substrate recognition, transport mechanism and pave the way for development of novel inhibitors against K. pneumoniae.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Suchanda Dey
- Biomics and Biodiversity lab, Siksha 'O' Anusandhan (deemed to be) University, Bhubaneswar, Odisha, India
| | - Madhusmita Rout
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Nalco Square, Bhubaneswar, Odisha, India
| | - Sanghamitra Pati
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Nalco Square, Bhubaneswar, Odisha, India
| | - Mahender Kumar Singh
- Data Science Laboratory, National Brain Research Centre, Gurgaon, Haryana, India
| | - Budheswar Dehury
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Nalco Square, Bhubaneswar, Odisha, India
| | - Enketeswara Subudhi
- Biomics and Biodiversity lab, Siksha 'O' Anusandhan (deemed to be) University, Bhubaneswar, Odisha, India
| |
Collapse
|
4
|
Chakraborty S, Deshmukh A, Kesari P, Bhaumik P. Toxoplasma gondii aspartic protease 5: structural basis of substrate binding and inhibition mechanism. J Biomol Struct Dyn 2024:1-16. [PMID: 38424737 DOI: 10.1080/07391102.2024.2322625] [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: 12/08/2023] [Accepted: 02/19/2024] [Indexed: 03/02/2024]
Abstract
Toxoplasma gondii, a worldwide prevalent parasite is responsible for causing toxoplasmosis in almost all warm-blooded animals, including humans. Golgi-resident T. gondii aspartic protease 5 (TgASP5) plays an essential role in the maturation and export of the effector proteins those modulate the host immune system to establish a successful infection. Hence, inhibiting this enzyme can be a possible therapeutic strategy against toxoplasmosis. This is the first report of the detailed structural investigations of the TgASP5 mature enzyme using molecular modeling and an all-atom simulation approach which provide in-depth knowledge of the active site architecture of TgASP5. The analysis of the binding mode of the TEXEL (Toxoplasma EXport Element) substrate to TgASP5 highlighted the importance of the active site residues. Ser505, Ala776 and Tyr689 in the S2 binding pocket are responsible for the specificity towards Arg at the P2 position of TEXEL substrate. The molecular basis of inhibition by the only known inhibitor RRLStatine has been identified, and our results show that it blocks the active site by forming a hydrogen bond with a catalytic aspartate. Besides that, known aspartic protease inhibitors were screened against TgASP5 using docking, MD simulations and MM-PBSA binding energy calculations. The top-ranked inhibitors (SC6, ZY1, QBH) showed higher binding energy than RRLStatine. Understanding the structural basis of substrate recognition and the binding mode of these inhibitors will help to develop potent mechanistic inhibitors against TgASP5. This study will also provide insights into the structural features of pepsin-like aspartic proteases from other apicomplexan parasites for developing antiparasitic agents.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Satadru Chakraborty
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Mumbai, India
| | - Anuradha Deshmukh
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Mumbai, India
| | - Pooja Kesari
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Mumbai, India
| | - Prasenjit Bhaumik
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Mumbai, India
| |
Collapse
|
5
|
Priyadarshinee M, Dehury B, Mishra S, Jena C, Patra M, Mishra NK, Samanta L, Mallick BC. Spectroscopic insights with molecular docking and molecular dynamic simulation studies of anticancer drug 5-Fluorouracil targeting human pyruvate kinase m2. J Biomol Struct Dyn 2024:1-13. [PMID: 38345048 DOI: 10.1080/07391102.2024.2313158] [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: 10/13/2023] [Accepted: 01/27/2024] [Indexed: 02/22/2024]
Abstract
This study was conducted to test the efficacy of 5-fluorouracil (5-FU) as an anticancer drug against the human pyruvate kinase isozyme M2 (PKM2) using spectroscopic, molecular docking and molecular dynamic simulation studies. PKM2 fluorescence quenching studies in the presence of 5-FU performed at three different temperatures indicates dynamic quenching processes with single-set of binding (n ≈ 1) profile. The biomolecular quenching constants (kq) and the effective binding constants (Kb) obtained are shown to increase with temperature. The calculated enthalpy (ΔH) and entropy changes (ΔS) are estimated to be -118.06 kJ/mol and 146.14 kJ/mol/K respectively, which suggest the possible mode of interaction as electrostatic and hydrogen bonding. Further, these values were used to estimate the free energy changes (ΔG) and that increases with temperature. The negative ΔG values clearly indicates spontaneous binding process that stabilizes the complex formed between 5-FU and PKM2. Far-UV CD spectra of PKM2 in the presence of 5-FU shows decrease in α-helix contents which point towards the destabilization of secondary structure that weakens the biological activity of PKM2. The intrinsic fluorescence study and circular dichroism (CD) spectra showed minor conformational changes of PKM2 in the presence of 5-FU. Additionally, the results obtained from molecular docking and all-atom molecular dynamic simulation study supports the insight of the spectroscopic binding studies, and strengthens the dynamic stability of the complex between 5-FU and PKM2 through H-bonding. This study establishes a paradigm of 5-FU-PKM2 complexation and the efficacy of 5-FU that compromises the biological activity of the targeted PKM2.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
| | - Budheswar Dehury
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Nalco Square, Chandrasekharpur, Bhubaneswar, India
| | - Sarbani Mishra
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Nalco Square, Chandrasekharpur, Bhubaneswar, India
| | | | | | - Neeraj K Mishra
- Department of Biotechnology, GITAM University, Vishakhapatnam, India
| | - Luna Samanta
- Department of Zoology, Ravenshaw University, Cuttack, India
| | - Bairagi C Mallick
- Department of Chemistry, Ravenshaw University, Cuttack, India
- Department of Chemistry, Central University of Jharkhand, Ranchi, India
| |
Collapse
|
6
|
Bhatt S, Choudhary P, Chatterjee S, Akhter Y. Comparative analysis of SilA-laccase mediated degradation of ciprofloxacin, norfloxacin and ofloxacin and interpretation of the possible catalytic mechanism. J Biomol Struct Dyn 2024; 42:425-434. [PMID: 37096761 DOI: 10.1080/07391102.2023.2197074] [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/09/2022] [Accepted: 03/13/2023] [Indexed: 04/26/2023]
Abstract
Fluoroquinolones (FQs) are the most commonly used antimicrobial drugs and regardless of their advantages in the healthcare sector, the pollution of these antimicrobial drugs in the environment has big concerns about human and environmental health. The presence of these antibiotic drugs even at the lowest concentrations in the environment has resulted in the emergence and spread of antibiotic resistance. Hence, it is necessary to remediate these pollutants from the environment. Previously alkaline laccase (SilA) from Streptomyces ipomoeae has been demonstrated to show degrading potentials against two of the FQs, Ciprofloxacin (CIP) and Norfloxacin (NOR); however, the molecular mechanism was not elucidated in detail. In this study, we have analyzed the possible molecular catalytic mechanism of FQ degrading SilA-laccase for the degradation of the FQs, CIP, NOR and Ofloxacin (OFL) using three-dimensional protein structure modeling, molecular docking and molecular dynamic (MD) studies. The comparative protein sequence analysis revealed the presence of tetrapeptide conserved catalytic motif, His102-X-His104-Gly105. After evaluating the active site of the enzyme in depth using CDD, COACH and S-site tools, we have identified the catalytic triad composed of three conserved amino acid residues, His102, Val103 and Tyr108 with which ligands interacted during the catalysis process. By analyzing the MD trajectories, it is revealed that the highest degradation potential of SilA is for CIP followed by NOR and OFL. Ultimately, this study provides the possible comparative catalytic mechanism for the degradation of CIP, NOR and OFL by the SilA enzyme.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Sunidhi Bhatt
- Bioremediation and Metabolomics Research Group, Department of Environmental Sciences, Central University of Himachal Pradesh, Distt-Kangara, Himachal Pradesh, India
| | - Priyanka Choudhary
- Bioremediation and Metabolomics Research Group, Department of Environmental Sciences, Central University of Himachal Pradesh, Distt-Kangara, Himachal Pradesh, India
| | - Subhankar Chatterjee
- Bioremediation and Metabolomics Research Group, Department of Ecology & Environmental Sciences, Pondicherry University, Kalapet, Puducherry, India
| | - Yusuf Akhter
- Department of Biotechnology, Babasaheb Bhimrao Ambedkar University, Lucknow, India
| |
Collapse
|
7
|
Rout M, Mishra S, Panda S, Dehury B, Pati S. Lipid and cholesterols modulate the dynamics of SARS-CoV-2 viral ion channel ORF3a and its pathogenic variants. Int J Biol Macromol 2024; 254:127986. [PMID: 37944718 DOI: 10.1016/j.ijbiomac.2023.127986] [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: 08/24/2023] [Revised: 11/06/2023] [Accepted: 11/07/2023] [Indexed: 11/12/2023]
Abstract
SARS-CoV-2 accessory protein, ORF3a is a putative ion channel which immensely contributes to viral pathogenicity by modulating host immune responses and virus-host interactions. Relatively high expression of ORF3a in diseased individuals and implication with inflammasome activation, apoptosis and autophagy inhibition, ratifies as an effective target for developing vaccines and therapeutics. Herein, we present the elusive dynamics of ORF3a-dimeric state using all-atoms molecular dynamics (MD) simulations at μ-seconds scale in a heterogeneous lipid-mimetic system in multiple replicates. Additionally, we also explore the effect of non-synonymous pathogenic mutations on ORF3a ion channel activity and viral pathogenicity in different SARS-CoV-2 variants using various structure-based protein stability (ΔΔG) tools and computational saturation mutagenesis. Our study ascertains the role of phosphatidylcholines and cholesterol in modulating the structure of ORF3a, which perturbs the size and flexibility of the polar cavity that allows permeation of large cations. Discrete trend in ion channel pore radius and area per lipid arises the premise that presence of lipids might also affect the overall conformation of ORF3a. MD structural-ensembles, in some replicates rationalize the crucial role of TM2 in maintaining the native structure of ORF3a. We also infer that loss of structural stability primarily grounds for pathogenicity in more than half of the pathogenic variants of ORF3a. Overall, the effect of mutation on alteration of ion permeability of ORF3a, proposed in this study brings mechanistic insights into variant consequences on viral membrane proteins of SARS-CoV-2, which can be utilized for the development of novel therapeutics to treat COVID-19 and other coronavirus diseases.
Collapse
Affiliation(s)
- Madhusmita Rout
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Chandrasekharpur, Bhubaneswar 751023, Odisha, India
| | - Sarbani Mishra
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Chandrasekharpur, Bhubaneswar 751023, Odisha, India
| | - Sunita Panda
- Mycology Division, ICMR-Regional Medical Research Centre, Chandrasekharpur, Bhubaneswar 751023, Odisha, India
| | - Budheswar Dehury
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Chandrasekharpur, Bhubaneswar 751023, Odisha, India.
| | - Sanghamitra Pati
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Chandrasekharpur, Bhubaneswar 751023, Odisha, India.
| |
Collapse
|
8
|
Cao Y, Wei H, Jiang S, Lu T, Nie P, Yang C, Liu N, Lee I, Meng X, Wang W, Yuan Z. Effect of AQP4 and its palmitoylation on the permeability of exogenous reactive oxygen species: Insights from computational study. Int J Biol Macromol 2023; 253:127568. [PMID: 37866582 DOI: 10.1016/j.ijbiomac.2023.127568] [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: 08/02/2023] [Revised: 10/16/2023] [Accepted: 10/18/2023] [Indexed: 10/24/2023]
Abstract
Aquaporin 4 (AQP4) facilitates the transport of reactive oxygen species (ROS). Both cancer cells and the ionizing radiation microenvironment can induce posttranslational modifications (PTMs) in AQP4, which may affect its permeability to ROS. Because this ROS diffusion process is rapid, microscopic, and instantaneous within and outside cells, conventional experimental methods are inadequate for elucidating the molecular mechanisms involved. In this study, computational methods were employed to investigate the permeability of exogenous ROS mediated by radiation in AQP4 at a molecular scale. We constructed a simulation system incorporating AQP4 and AQP4-Cysp13 in a complex lipid environment with ROS. Long-timescale molecular dynamics simulations were conducted to assess the structural stability of both AQP4 and AQP4-Cysp13. Free energy calculations were utilized to determine the ROS transport capability of the two AQP4 proteins. Computational electrophysiology and channel structural analysis quantitatively evaluated changes in ROS transport capacity under various radiation-induced transmembrane voltage microenvironments. Our findings demonstrate the distinct transport capabilities of AQP4 channels for water molecules and various types of ROS and reveal a decrease in transport efficiency when AQP4 undergoes palmitoylation modification. In addition, we have simulated the radiation-induced alteration of cell membrane voltage, which significantly affected the ROS transport capacity. We propose that this research will enhance the understanding of the molecular mechanisms governing the transport of exogenous ROS by AQP4 and elucidate the influence of palmitoylation on ROS transport. This study will also help clarify how different structural features of AQP4 affect the transport of exogenous ROS mediated by radiotherapy, thereby providing a theoretical molecular basis for the development of new treatment strategies that combine with radiotherapy.
Collapse
Affiliation(s)
- Yipeng Cao
- Key Laboratory of Cancer Prevention and Therapy, Tianjin's Clinical Research Center for Cancer, National Clinical Research Center for Cancer, Tianjin Medical University Cancer Institute and Hospital, 300060, PR China; National Supercomputer Center in Tianjin, 300457, PR China.
| | - Hui Wei
- Key Laboratory of Cancer Prevention and Therapy, Tianjin's Clinical Research Center for Cancer, National Clinical Research Center for Cancer, Tianjin Medical University Cancer Institute and Hospital, 300060, PR China
| | - Shengpeng Jiang
- Key Laboratory of Cancer Prevention and Therapy, Tianjin's Clinical Research Center for Cancer, National Clinical Research Center for Cancer, Tianjin Medical University Cancer Institute and Hospital, 300060, PR China
| | - Tong Lu
- Key Laboratory of Cancer Prevention and Therapy, Tianjin's Clinical Research Center for Cancer, National Clinical Research Center for Cancer, Tianjin Medical University Cancer Institute and Hospital, 300060, PR China
| | - Pengfei Nie
- National Supercomputer Center in Tianjin, 300457, PR China
| | - Chengwen Yang
- Key Laboratory of Cancer Prevention and Therapy, Tianjin's Clinical Research Center for Cancer, National Clinical Research Center for Cancer, Tianjin Medical University Cancer Institute and Hospital, 300060, PR China
| | - Ningbo Liu
- Key Laboratory of Cancer Prevention and Therapy, Tianjin's Clinical Research Center for Cancer, National Clinical Research Center for Cancer, Tianjin Medical University Cancer Institute and Hospital, 300060, PR China
| | - Imshik Lee
- College of Physics, Nankai University, Tianjin 300071, PR China
| | - Xiangfei Meng
- National Supercomputer Center in Tianjin, 300457, PR China.
| | - Wei Wang
- Key Laboratory of Cancer Prevention and Therapy, Tianjin's Clinical Research Center for Cancer, National Clinical Research Center for Cancer, Tianjin Medical University Cancer Institute and Hospital, 300060, PR China.
| | - Zhiyong Yuan
- Key Laboratory of Cancer Prevention and Therapy, Tianjin's Clinical Research Center for Cancer, National Clinical Research Center for Cancer, Tianjin Medical University Cancer Institute and Hospital, 300060, PR China.
| |
Collapse
|
9
|
Mishra S, Rout M, Singh MK, Dehury B, Pati S. Illuminating the structural basis of human neurokinin 1 receptor (NK1R) antagonism through classical all-atoms molecular dynamics simulations. J Cell Biochem 2023; 124:1848-1869. [PMID: 37942587 DOI: 10.1002/jcb.30493] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2023] [Revised: 09/26/2023] [Accepted: 10/16/2023] [Indexed: 11/10/2023]
Abstract
Advances in structural biology have bestowed insights into the pleiotropic effects of neurokinin 1 receptors (NK1R) in diverse patho-physiological processes, thereby highlighting the potential therapeutic value of antagonists directed against NK1R. Herein, we investigate the mode of antagonist recognition to discern the obscure atomic facets germane for the function and molecular determinants of NK1R. To commence discernment of potent antagonists and the conformational changes in NK1R, induced upon antagonist binding, state-of-the-art classical all-atoms molecular dynamics (MD) simulations in lipid mimetic bilayers have been utilized. MD simulations of structural ensembles reveals the involvement of TM5 and TM6 in tight anchoring of antagonists through a network of interhelical hydrogen-bonds, while, the extracellular loop 2 (ECL2) governs the overall size and nature of the pocket, thereby modulating NK1R. Consistent comparison between experiments and MD simulation results discerns the predominant role of TM3, TM4, and TM6 in lipid-NK1R interaction. Correlation between hydrophobic index and helicity of TM domains elucidates their importance in maintaining the structural stability in addition to regulating NK1R antagonism. Taken together, we anticipate that our computational study marks a comprehensive structural basis of NK1R antagonism in lipid bilayers, which may facilitate designing of new therapeutics against associated diseases targeting human neurokinin receptors.
Collapse
Affiliation(s)
- Sarbani Mishra
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Bhubaneswar, Odisha, India
| | - Madhusmita Rout
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Bhubaneswar, Odisha, India
| | - Mahender Kumar Singh
- Data Science Laboratory, National Brain Research Centre, Gurgaon, Haryana, India
| | - Budheswar Dehury
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Bhubaneswar, Odisha, India
| | - Sanghamitra Pati
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Bhubaneswar, Odisha, India
| |
Collapse
|
10
|
Salvatierra A, Mateluna P, Toro G, Solís S, Pimentel P. Genome-Wide Identification and Gene Expression Analysis of Sweet Cherry Aquaporins ( Prunus avium L.) under Abiotic Stresses. Genes (Basel) 2023; 14:genes14040940. [PMID: 37107698 PMCID: PMC10138167 DOI: 10.3390/genes14040940] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2023] [Revised: 04/13/2023] [Accepted: 04/14/2023] [Indexed: 04/29/2023] Open
Abstract
Aquaporins (AQPs) are integral transmembrane proteins well known as channels involved in the mobilization of water, small uncharged molecules and gases. In this work, the main objective was to carry out a comprehensive study of AQP encoding genes in Prunus avium (cv. Mazzard F12/1) on a genome-wide scale and describe their transcriptional behaviors in organs and in response to different abiotic stresses. A total of 28 non-redundant AQP genes were identified in Prunus spp. Genomes, which were phylogenetically grouped into five subfamilies (seven PIPs, eight NIPs, eight TIPs, three SIPs and two XIPs). Bioinformatic analyses revealed a high synteny and remarkable conservation of structural features among orthologs of different Prunus genomes. Several cis-acting regulatory elements (CREs) related to stress regulation were detected (ARE, WRE3, WUN, STRE, LTR, MBS, DRE, AT-rich and TC-rich). The above could be accounting for the expression variations associated with plant organs and, especially, each abiotic stress analyzed. Gene expressions of different PruavAQPs were shown to be preferentially associated with different stresses. PruavXIP2;1 and PruavXIP1;1 were up-regulated in roots at 6 h and 72 h of hypoxia, and in PruavXIP2;1 a slight induction of expression was also detected in leaves. Drought treatment strongly down-regulated PruavTIP4;1 but only in roots. Salt stress exhibited little or no variation in roots, except for PruavNIP4;1 and PruavNIP7;1, which showed remarkable gene repression and induction, respectively. Interestingly, PruavNIP4;1, the AQP most expressed in cherry roots subjected to cold temperatures, also showed this pattern in roots under high salinity. Similarly, PruavNIP4;2 consistently was up-regulated at 72 h of heat and drought treatments. From our evidence is possible to propose candidate genes for the development of molecular markers for selection processes in breeding programs for rootstocks and/or varieties of cherry.
Collapse
Affiliation(s)
- Ariel Salvatierra
- Centro de Estudios Avanzados en Fruticultura (CEAF), Camino Las Parcelas 882, km 105 Ruta 5 Sur, Sector Los Choapinos, Rengo 2940000, Chile
| | - Patricio Mateluna
- Centro de Estudios Avanzados en Fruticultura (CEAF), Camino Las Parcelas 882, km 105 Ruta 5 Sur, Sector Los Choapinos, Rengo 2940000, Chile
| | - Guillermo Toro
- Centro de Estudios Avanzados en Fruticultura (CEAF), Camino Las Parcelas 882, km 105 Ruta 5 Sur, Sector Los Choapinos, Rengo 2940000, Chile
| | - Simón Solís
- Centro de Estudios Avanzados en Fruticultura (CEAF), Camino Las Parcelas 882, km 105 Ruta 5 Sur, Sector Los Choapinos, Rengo 2940000, Chile
| | - Paula Pimentel
- Centro de Estudios Avanzados en Fruticultura (CEAF), Camino Las Parcelas 882, km 105 Ruta 5 Sur, Sector Los Choapinos, Rengo 2940000, Chile
| |
Collapse
|
11
|
Feng S, Park S, Choi YK, Im W. CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications. J Chem Theory Comput 2023; 19:2161-2185. [PMID: 37014931 PMCID: PMC10174225 DOI: 10.1021/acs.jctc.2c01246] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/06/2023]
Abstract
Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.
Collapse
Affiliation(s)
- Shasha Feng
- Departments of Biological Sciences and Chemistry, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Soohyung Park
- Departments of Biological Sciences and Chemistry, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Yeol Kyo Choi
- Departments of Biological Sciences and Chemistry, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Wonpil Im
- Departments of Biological Sciences and Chemistry, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| |
Collapse
|
12
|
Cetintas VB, Duzgun Z, Akalin T, Ozgiray E, Dogan E, Yildirim Z, Akinturk N, Biceroglu H, Ertan Y, Kosova B. Molecular dynamic simulation and functional analysis of pathogenic PTEN mutations in glioblastoma. J Biomol Struct Dyn 2023; 41:11471-11483. [PMID: 36591942 DOI: 10.1080/07391102.2022.2162582] [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: 08/17/2022] [Accepted: 12/20/2022] [Indexed: 01/03/2023]
Abstract
PTEN, a dual-phosphatase and scaffold protein, is one of the most commonly mutated tumour suppressor gene across various cancer types in human. The aim of this study therefore was to investigate the stability, structural and functional effects, and pathogenicity of 12 missense PTEN mutations (R15S, E18G, G36R, N49I, Y68H, I101T, C105F, D109N, V133I, C136Y, R173C and N276S) found by next generation sequencing of the PTEN gene in tissue samples obtained from glioblastoma patients. Computational tools and molecular dynamic simulation programs were used to identify the deleterious effects of these mutations. Furthermore, PTEN mRNA and protein expression levels were evaluated by qRT-PCR, Western Blot, and immunohistochemistry staining methods. Various computational tools predicted strong deleterious effects for the G36R, C105F, C136Y and N276S mutations. Molecular dynamic simulation revealed a significant decrease in protein stability for the Y68H and N276S mutations when compared with the wild type protein; whereas, C105F, D109N, V133I and R173C showed partial stability reduction. Significant residual fluctuations were observed in the R15S, N49I and C136Y mutations and radius of gyration graphs revealed the most compact structure for D109N and least for C136Y. In summary, our study is the first one to show the presence of PTEN E18G, N49I, D109N and N276S mutations in glioblastoma patients; where, D109N is neutral and N276S is a damaging and disease-associated mutation.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
| | - Zekeriya Duzgun
- Department of Medical Biology, Giresun University Faculty of Medicine, Giresun, Turkey
| | - Taner Akalin
- Department of Pathology, Ege University Faculty of Medicine, Izmir, Turkey
| | - Erkin Ozgiray
- Department of Neurosurgery, Ege University Faculty of Medicine, Izmir, Turkey
| | - Eda Dogan
- Department of Medical Biology, Ege University Faculty of Medicine, Izmir, Turkey
| | - Zafer Yildirim
- Department of Medical Biology, Ege University Faculty of Medicine, Izmir, Turkey
| | - Nevhis Akinturk
- Department of Neurosurgery, Ege University Faculty of Medicine, Izmir, Turkey
| | - Huseyin Biceroglu
- Department of Neurosurgery, Ege University Faculty of Medicine, Izmir, Turkey
| | - Yesim Ertan
- Department of Pathology, Ege University Faculty of Medicine, Izmir, Turkey
| | - Buket Kosova
- Department of Medical Biology, Ege University Faculty of Medicine, Izmir, Turkey
| |
Collapse
|
13
|
Structural insights into the RNA interaction with Yam bean Mosaic virus (coat protein) from Pachyrhizus erosus using bioinformatics approach. PLoS One 2022; 17:e0270534. [PMID: 35867657 PMCID: PMC9307209 DOI: 10.1371/journal.pone.0270534] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2022] [Accepted: 06/11/2022] [Indexed: 11/19/2022] Open
Abstract
Plants are constantly threatened by a virus infection, i.e., Potyviruses, the second largest genus of plant viruses which results in several million-dollar losses in various essential crops globally. Yam bean (Pachyrhizus erosus) is considered to be one of the essential tuberous legume crops holding a great potential source of starch. Yam Bean Mosaic Virus (YBMV) of Potyvirus group belonging to the family potyviridae affects Yam bean and several angiosperms both in the tropical and sub-tropical regions causing large economical losses in crops. In this study, we attempted to understand the sequence-structure relationship and mode of RNA binding mechanism in YBMV CP using in silico integrative modeling and all-atoms molecular dynamics (MD) simulations. The assembly of coat protein (CP) subunits from YBMV and the plausible mode of RNA binding were compared with the experimental structure of CP from Watermelon mosaic virus potyvirus (5ODV). The transmembrane helix region is present in the YBMV CP sequence ranging from 76 to 91 amino acids. Like the close structural-homolog, 24 CPs monomeric sub-units formed YBMV a conserved fold. Our computational study showed that ARG124, ARG155, and TYR151 orient towards the inner side of the virion, while, THR122, GLN125, SER92, ASP94 reside towards the outer side of the virion. Despite sharing very low sequence similarity with CPs from other plant viruses, the strongly conserved residues Ser, Arg, and Asp within the RNA binding pocket of YBMV CP indicate the presence of a highly conserved RNA binding site in CPs from different families. Using several bioinformatics tools and comprehensive analysis from MD simulation, our study has provided novel insights into the RNA binding mechanism in YBMV CP. Thus, we anticipate that our findings from this study will be useful for the development of new therapeutic agents against the pathogen, paving the way for researchers to better control this destructive plant virus.
Collapse
|