1
|
Abdel-Naim AB, Kumar P, Bazuhair MA, Rizg WY, Niyazi HA, Alkuwaity K, Niyazi HA, Alharthy SA, Harakeh S, Haque S, Prakash A, Kumar V. Computational insights into dynamics and conformational stability of N-acetylmannosamine kinase mutations. J Biomol Struct Dyn 2024:1-11. [PMID: 38502682 DOI: 10.1080/07391102.2024.2323702] [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: 05/02/2023] [Accepted: 02/21/2024] [Indexed: 03/21/2024]
Abstract
The activity of UDP-GlcNAc 2-epimerase/ManNAc kinase (GNE) is essential for the biosynthesis of sialic acid, which is involved in cellular processes in health and diseases. GNE contains an N-terminal epimerase domain and a C-terminal kinase domain (N-acetylmannosamine kinase, MNK). Mutations of the GNE protein led to hypoactivity of the enzyme and cause sialurea or autosomal recessive inclusion body myopathy/Nonaka myopathy. Here, we used all-atom molecular dynamics (MD) simulations to comprehend the folding, dynamics and conformational stability of MNK variants, including the wild type (WT) and three mutants (H677R, V696M and H677R/V696M). The deleterious and destabilizing nature of MNK mutants were predicted using different prediction tools. Results predicted that mutations modulate the stability, flexibility and function of MNK. The effect of mutations on the conformational stability and dynamics of MNK was next studied through the free-energy landscape (FEL), hydrogen-bonds and secondary structure changes. The FEL results show that the mutations interfere with various conformational transitions in both WT and mutants, exposing the structural underpinnings of protein destabilization and unfolding brought on by mutation. We discover that, when compared to the other two mutations, V696M and H677R/V696M, H677R has the most harmful effects. These findings have a strong correlation with published experimental studies that demonstrate how these mutations disrupt MNK activity. Hence, this computational study describes the structural details to unravel the mutant effects at the atomistic resolution and has implications for understanding the GNE's physiological and pathological role.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Ashraf B Abdel-Naim
- Department of Pharmacology and Toxicology, Faculty of Pharmacy, King Abdulaziz University, Jeddah, Saudi Arabia
- Mohamed Saeed Tamer Chair for Pharmaceutical Industries, Faculty of Pharmacy, King Abdulaziz University, Jeddah, Saudi Arabia
- Center of Excellence for Drug Research and Pharmaceutical Industries, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Pawan Kumar
- School of Computational & Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Mohammed A Bazuhair
- Center of Excellence for Drug Research and Pharmaceutical Industries, King Abdulaziz University, Jeddah, Saudi Arabia
- Department of Clinical Pharmacology, Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Waleed Y Rizg
- Mohamed Saeed Tamer Chair for Pharmaceutical Industries, Faculty of Pharmacy, King Abdulaziz University, Jeddah, Saudi Arabia
- Center of Excellence for Drug Research and Pharmaceutical Industries, King Abdulaziz University, Jeddah, Saudi Arabia
- Department of Pharmaceutics, Faculty of Pharmacy, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Hatoon A Niyazi
- Department of Pharmacology and Toxicology, Faculty of Pharmacy, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Khalil Alkuwaity
- Department of Pharmacology and Toxicology, Faculty of Pharmacy, King Abdulaziz University, Jeddah, Saudi Arabia
- Vaccines and Immunotherapy Unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Hanouf A Niyazi
- Department of Pharmacology and Toxicology, Faculty of Pharmacy, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Saif A Alharthy
- Department of Medical Laboratory Sciences, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah, Saudi Arabia
- Toxicology and Forensic Sciences Unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Steve Harakeh
- King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia
- Yousef Abdul Latif Jameel Scientific Chair of Prophetic Medicine Application, Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Shafiul Haque
- Research and Scientific Studies Unit, College of Nursing and Health Sciences, Jazan University, Jazan, Saudi Arabia
- Gilbert and Rose-Marie Chagoury School of Medicine, Lebanese American University, Beirut, Lebanon
- Centre of Medical and Bio-Allied Health Sciences Research, Ajman University, Ajman, United Arab Emirates
| | - Amresh Prakash
- Amity Institute of Integrative Sciences and Health, Amity University Haryana, Gurgaon, India
| | - Vijay Kumar
- Amity Institute of Neuropsychology & Neurosciences, Amity University, Noida, India
| |
Collapse
|
2
|
Alamri SH, Haque S, Alghamdi BS, Tayeb HO, Azhari S, Farsi RM, Elmokadem A, Alamri TA, Harakeh S, Prakash A, Kumar V. Comprehensive mapping of mutations in TDP-43 and α-Synuclein that affect stability and binding. J Biomol Struct Dyn 2023:1-13. [PMID: 38126188 DOI: 10.1080/07391102.2023.2293258] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2023] [Accepted: 11/11/2023] [Indexed: 12/23/2023]
Abstract
Abnormal aggregation and amyloid inclusions of TAR DNA-binding protein 43 (TDP-43) and α-Synuclein (α-Syn) are frequently co-observed in amyotrophic lateral sclerosis, Parkinson's disease, and Alzheimer's disease. Several reports showed TDP-43 C-terminal domain (CTD) and α-Syn interact with each other and the aggregates of these two proteins colocalized together in different cellular and animal models. Molecular dynamics simulation was conducted to elucidate the stability of the TDP-43 and Syn complex structure. The interfacial mutations in protein complexes changes the stability and binding affinity of the protein that may cause diseases. Here, we have utilized the computational saturation mutagenesis approach including structure-based stability and binding energy calculations to compute the systemic effects of missense mutations of TDP-43 CTD and α-Syn on protein stability and binding affinity. Most of the interfacial mutations of CTD and α-Syn were found to destabilize the protein and reduced the protein binding affinity. The results thus shed light on the functional consequences of missense mutations observed in TDP-43 associated proteinopathies and may provide the mechanisms of co-morbidities involving these two proteins.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Sultan H Alamri
- Department of Family Medicine, Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Shafiul Haque
- Research and Scientific Studies Unit, College of Nursing and Allied Health Sciences, Jazan University, Jazan, Saudi Arabia
- Gilbert and Rose-Marie Chagoury School of Medicine, Lebanese American University, Beirut, Lebanon
- Centre of Medical and Bio-Allied Health Sciences Research, Ajman University, Ajman, United Arab Emirates
| | - Badra S Alghamdi
- Department of Physiology, Neuroscience Unit, Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
- Pre-Clinical Research Unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Haythum O Tayeb
- The Mind and Brain Studies Initiative, Neuroscience Research Unit, Department of Neurology, Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Shereen Azhari
- Department of Biological Sciences, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Reem M Farsi
- Department of Biological Sciences, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Abear Elmokadem
- Department of Hematology/Pediatric Oncology, King Abdulaziz University Hospital, Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Turki A Alamri
- Family and Community Medicine Department, Faculty of Medicine in Rabigh, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Steve Harakeh
- King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Jeddah, Saudi Arabia
- Yousef Abdul Latif Jameel Scientific Chair of Prophetic Medicine Application, Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Amresh Prakash
- Amity Institute of Integrative Sciences and Health (AIISH), Amity University Haryana, Gurgaon, India
| | - Vijay Kumar
- Amity Institute of Neuropsychology & Neurosciences, Amity University, Noida, India
| |
Collapse
|
3
|
Haque S, Kumar P, Mathkor DM, Bantun F, Jalal NA, Mufti AH, Prakash A, Kumar V. In silico evaluation of the inhibitory potential of nucleocapsid inhibitors of SARS-CoV-2: a binding and energetic perspective. J Biomol Struct Dyn 2023; 41:9797-9807. [PMID: 36379684 DOI: 10.1080/07391102.2022.2146752] [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/17/2022] [Accepted: 11/07/2022] [Indexed: 11/17/2022]
Abstract
The COVID-19 outbreak brought on by the SARS-CoV-2 virus continued to infect a sizable population worldwide. The SARS-CoV-2 nucleocapsid (N) protein is the most conserved RNA-binding structural protein and is a desirable target because of its involvement in viral transcription and replication. Based on this aspect, this study focused to repurpose antiviral compounds approved or in development for treating COVID-19. The inhibitors chosen are either FDA-approved or are currently being studied in clinical trials against COVID-19. Initially, they were designed to target stress granules and other RNA biology. We have utilized structure-based molecular docking and all-atom molecular dynamics (MD) simulation approach to investigate in detail the binding energy and binding modes of the different anti-N inhibitors to N protein. The result showed that five drugs including Silmitasterib, Ninetanidinb, Ternatin, Luteolin, Fedratinib, PJ34, and Zotatafin were found interacting with RNA binding sites as well as to predicted protein interface with higher binding energy. Overall, drug binding increases the stability of the complex with maximum stability found in the order, Silmitasertib > PJ34 > Zotatatafin. In addition, the frustration changes due to drug binding brings a decrease in local frustration and this decrease is mainly observed in α-helix, β3, β5, and β6 strands and are important for drug binding. Our in-silico data suggest that an effective interaction occurs for some of the tested drugs and prompt their further validation to reduce the rapid outspreading of SARS-CoV-2.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Shafiul Haque
- Research and Scientific Studies Unit, College of Nursing and Allied Health Sciences, Jazan University, Jazan, Saudi Arabia
| | - Pawan Kumar
- School of Computational & Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Darin Mansor Mathkor
- Research and Scientific Studies Unit, College of Nursing and Allied Health Sciences, Jazan University, Jazan, Saudi Arabia
| | - Farkad Bantun
- Department of Microbiology, Faculty of Medicine, Umm Al-Qura University, Makkah, Saudi Arabia
| | - Naif A Jalal
- Department of Microbiology, Faculty of Medicine, Umm Al-Qura University, Makkah, Saudi Arabia
| | - Ahmad Hasan Mufti
- Medical Genetics Department, Faculty of Medicine, Umm Al-Qura University, Makkah, Saudi Arabia
| | - Amresh Prakash
- Amity Institute of Integrative Sciences and Health, Amity University Haryana, Gurgaon, India
| | - Vijay Kumar
- Amity Institute of Neuropsychology & Neurosciences, Amity University, Noida, Uttar Pradesh, India
| |
Collapse
|
4
|
Tan Y, Chen Y, Liu X, Tang Y, Lao Z, Wei G. Dissecting how ALS-associated D290V mutation enhances pathogenic aggregation of hnRNPA2 286-291 peptides: Dynamics and conformational ensembles. Int J Biol Macromol 2023; 241:124659. [PMID: 37119915 DOI: 10.1016/j.ijbiomac.2023.124659] [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: 03/20/2023] [Revised: 04/22/2023] [Accepted: 04/24/2023] [Indexed: 05/01/2023]
Abstract
The aggregation of RNA binding proteins, including hnRNPA1/2, TDP-43 and FUS, is heavily implicated in causing or increasing disease risk for a series of neurodegenerative diseases such as amyotrophic lateral sclerosis (ALS). A recent experimental study demonstrated that an ALS-related D290V mutation in the low complexity domain (LCD) of hnRNPA2 can enhance the aggregation propensity of wild type (WT) hnRNPA2286-291 peptide. However, the underlying molecular mechanisms remain elusive. Herein, we investigated effects of D290V mutation on aggregation dynamics of hnRNPA2286-291 peptide and the conformational ensemble of hnRNPA2286-291 oligomers by performing all-atom molecular dynamic and replica-exchange molecular dynamic simulations. Our simulations demonstrate that D290V mutation greatly reduces the dynamics of hnRNPA2286-291 peptide and that D290V oligomers possess higher compactness and β-sheet content than WT, indicative of mutation-enhanced aggregation capability. Specifically, D290V mutation strengthens inter-peptide hydrophobic, main-chain hydrogen bonding and side-chain aromatic stacking interactions. Those interactions collectively lead to the enhancement of aggregation capability of hnRNPA2286-291 peptides. Overall, our study provides insights into the dynamics and thermodynamic mechanisms underlying D290V-induced disease-causing aggregation of hnRNPA2286-291, which could contribute to better understanding of the transitions from reversible condensates to irreversible pathogenic aggregates of hnRNPA2 LCD in ALS-related diseases.
Collapse
Affiliation(s)
- Yuan Tan
- Department of Physics, Fudan University, Shanghai 200438, People's Republic of China; State Key Laboratory of Surface Physics, Fudan University, Shanghai 200438, People's Republic of China; Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200438, People's Republic of China
| | - Yujie Chen
- Department of Physics, Fudan University, Shanghai 200438, People's Republic of China; State Key Laboratory of Surface Physics, Fudan University, Shanghai 200438, People's Republic of China; Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200438, People's Republic of China
| | - Xianshi Liu
- Department of Physics, Fudan University, Shanghai 200438, People's Republic of China; State Key Laboratory of Surface Physics, Fudan University, Shanghai 200438, People's Republic of China; Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200438, People's Republic of China
| | - Yiming Tang
- Department of Physics, Fudan University, Shanghai 200438, People's Republic of China; State Key Laboratory of Surface Physics, Fudan University, Shanghai 200438, People's Republic of China; Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200438, People's Republic of China
| | - Zenghui Lao
- Department of Physics, Fudan University, Shanghai 200438, People's Republic of China; State Key Laboratory of Surface Physics, Fudan University, Shanghai 200438, People's Republic of China; Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200438, People's Republic of China
| | - Guanghong Wei
- Department of Physics, Fudan University, Shanghai 200438, People's Republic of China; State Key Laboratory of Surface Physics, Fudan University, Shanghai 200438, People's Republic of China; Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200438, People's Republic of China.
| |
Collapse
|
5
|
Liu X, Li X, Qiao Q, Li F, Wei G. ALS-Linked A315T and A315E Mutations Enhance β-Barrel Formation of the TDP-43 307-319 Hexamer: A REST2 Simulation Study. ACS Chem Neurosci 2023; 14:1310-1320. [PMID: 36888995 DOI: 10.1021/acschemneuro.3c00012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/10/2023] Open
Abstract
Pathogenic mutations of transactivation response element DNA-binding protein 43 (TDP-43) are closely linked with amyotrophic lateral sclerosis (ALS). It was recently reported that two ALS-linked familial mutants A315T and A315E of TDP-43307-319 peptides can self-assemble into oligomers including tetramers, hexamers, and octamers, among which hexamers were suggested to form the β-barrel structure. However, due to the transient nature of oligomers, their conformational properties and the atomic mechanisms underlying the β-barrel formation remain largely elusive. Herein, we investigated the hexameric conformational distributions of the wild-type (WT) TDP-43307-319 fragment and its A315T and A315E mutants by performing all-atom explicit-solvent replica exchange with solute tempering 2 simulations. Our simulations reveal that each peptide can self-assemble into diverse conformations including ordered β-barrels, bilayer β-sheets and/or monolayer β-sheets, and disordered complexes. A315T and A315E mutants display higher propensity to form β-barrel structures than the WT, which provides atomic explanation for their enhanced neurotoxicity reported previously. Detailed interaction analysis shows that A315T and A315E mutations increase inter-molecular interactions. Also, the β-barrel structures formed by the three different peptides are stabilized by distinct inter-peptide side-chain hydrogen bonding, hydrophobic, and aromatic stacking interactions. This study demonstrates the enhanced β-barrel formation of the TDP-43307-319 hexamer by the pathogenic A315T and A315E mutations and reveals the underlying molecular determinants, which may be helpful for in-depth understanding of the ALS-mutation-induced neurotoxicity of TDP-43 protein.
Collapse
Affiliation(s)
- Xianshi Liu
- Department of Physics, State Key Laboratory of Surface Physics, and Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200438, People's Republic of China
| | - Xuhua Li
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, People's Republic of China
| | - Qin Qiao
- Digital Medical Research Center, School of Basic Medical Sciences, Fudan University, Shanghai 200032, China.,Shanghai Key Laboratory of Medical Image Computing and Computer Assisted Intervention, Shanghai 200032, China
| | - Fangying Li
- Department of Physics, State Key Laboratory of Surface Physics, and Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200438, People's Republic of China
| | - Guanghong Wei
- Department of Physics, State Key Laboratory of Surface Physics, and Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200438, People's Republic of China
| |
Collapse
|
6
|
Liu X, Lao Z, Li X, Dong X, Wei G. ALS-associated A315E and A315pT variants exhibit distinct mechanisms in inducing irreversible aggregation of TDP-43 312-317 peptides. Phys Chem Chem Phys 2022; 24:16263-16273. [PMID: 35758309 DOI: 10.1039/d2cp01625g] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Amyotrophic lateral sclerosis (ALS) is intensively associated with insoluble aggregates formed by transactivation response element DNA-binding protein 43 (TDP-43) in the cytoplasm of neuron cells. A recent experimental study reported that two ALS-linked familial variants, A315E and A315pT (pT, phosphorylated threonine), can induce irreversible aggregation of the TDP-43 312NFGAFS317 segment (TDP-43312-317). However, the underlying molecular mechanism remains largely elusive. Here, we investigated the early aggregation process of the wild type (WT) 312NFGAFS317 segment and its A315E and A315pT variants by performing multiple microsecond all-atom molecular dynamics simulations. Our simulations show that the two variants display lower fluidity than WT, consistent with their decreased labilities observed in previous denaturation assay experiments. Despite each of the two variants carrying one negative charge, unexpectedly, we find that both A315E mutation and A315pT phosphorylation enhance intermolecular interactions and result in the formation of more compact oligomers. Compared to WT, A315E oligomers possess low β-sheet content but a compact hydrophobic core, while A315pT oligomers have high β-sheet content and large β-sheets. Side chain hydrogen-bonding and hydrophobic interactions as well as N312-E315 salt bridges contribute most to the increased aggregation propensity of the A315E mutant. By contrast, main chain and side chain hydrogen-bonding interactions, side chain hydrophobic and aromatic interactions, are crucial to the enhanced aggregation capability of the A315pT variant. These results indicate that glutamate mutation and phosphorylation at position 315 induce the irreversible aggregation of TDP-43312-317 peptides through differential mechanisms, which remind us that we should be careful in the investigation of the phosphorylation effect on protein aggregation by using phosphomimetic substitutions. This study provides mechanistic insights into the A315E/A315pT-induced irreversible aggregation of TDP-43312-317, which may be helpful for the in-depth understanding of ALS-mutation/phosphorylation-associated liquid-to-solid phase transition of TDP-43 protein aggregates.
Collapse
Affiliation(s)
- Xianshi Liu
- Department of Physics, State Key Laboratory of Surface Physics, and Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200438, P. R. China.
| | - Zenghui Lao
- Department of Physics, State Key Laboratory of Surface Physics, and Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200438, P. R. China.
| | - Xuhua Li
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, P. R. China
| | - Xuewei Dong
- Department of Physics, State Key Laboratory of Surface Physics, and Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200438, P. R. China.
| | - Guanghong Wei
- Department of Physics, State Key Laboratory of Surface Physics, and Key Laboratory for Computational Physical Sciences (Ministry of Education), Fudan University, Shanghai 200438, P. R. China.
| |
Collapse
|
7
|
Sharma A, Dey P. Novel insights into the structural changes induced by disease-associated mutations in TDP-43: a computational approach. J Biomol Struct Dyn 2022:1-11. [PMID: 35751132 DOI: 10.1080/07391102.2022.2092551] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Abstract
Over the last two decades, the pathogenic aggregation of TAR DNA-binding protein 43 (TDP-43) is found to be strongly associated with several fatal neurodegenerative diseases such as amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTD), etc. While the mutations and truncation in TDP-43 protein have been suggested to be responsible for TDP-43 pathogenesis by accelerating the aggregation process, the effects of these mutations on the bio-mechanism of pathological TDP-43 protein remained poorly understood. Investigating this at the molecular level, we formulized an integrated workflow of molecular dynamic simulation and machine learning models (MD-ML). By performing an extensive structural analysis of three disease-related mutations (i.e., I168A, D169G, and I168A-D169G) in the conserved RNA recognition motifs (RRM1) of TDP-43, we observed that the I168A-D169G double mutant delineates the highest packing of the protein inner core as compared to the other mutations, which may indicate more stability and higher chances of pathogenesis. Moreover, through our MD-ML workflow, we identified the biological descriptors of TDP-43 which includes the interacting residue pairs and individual protein residues that influence the stability of the protein and could be experimentally evaluated to develop potential therapeutic strategies.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Abhibhav Sharma
- School of Computer and System Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Pinki Dey
- School of Biotechnology and Biomolecular Sciences, University of New South Wales, Sydney, Australia
| |
Collapse
|
8
|
Li R, Singh R, Kashav T, Yang C, Sharma RD, Lynn AM, Prasad R, Prakash A, Kumar V. Computational Insights of Unfolding of N-Terminal Domain of TDP-43 Reveal the Conformational Heterogeneity in the Unfolding Pathway. Front Mol Neurosci 2022; 15:822863. [PMID: 35548668 PMCID: PMC9083116 DOI: 10.3389/fnmol.2022.822863] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2021] [Accepted: 03/22/2022] [Indexed: 02/05/2023] Open
Abstract
TDP-43 proteinopathies is a disease hallmark that characterizes amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD). The N-terminal domain of TDP-43 (NTD) is important to both TDP-43 physiology and TDP-43 proteinopathy. However, its folding and dimerization process is still poorly characterized. In the present study, we have investigated the folding/unfolding of NTD employing all-atom molecular dynamics (MD) simulations in 8 M dimethylsulfoxide (DMSO) at high temperatures. The MD results showed that the unfolding of the NTD at high temperature evolves through the formation of a number of conformational states differing in their stability and free energy. The presence of structurally heterogeneous population of intermediate ensembles was further characterized by the different extents of solvent exposure of Trp80 during unfolding. We suggest that these non-natives unfolded intermediate ensembles may facilitate NTD oligomerization and subsequently TDP-43 oligomerization, which might lead to the formation of irreversible pathological aggregates, characteristics of disease pathogenesis.
Collapse
Affiliation(s)
- Ruiting Li
- School of Engineering, Guangzhou College of Technology and Business, Guangzhou, China
| | - Ruhar Singh
- School of Computational & Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Tara Kashav
- Department of Life Science, Central University of South Bihar, Gaya, India
| | - ChunMin Yang
- School of Engineering, Guangzhou College of Technology and Business, Guangzhou, China
| | - Ravi Datta Sharma
- Amity Institute of Biotechnology, Amity University Haryana, Gurgaon, India
| | - Andrew M. Lynn
- School of Computational & Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Rajendra Prasad
- Amity Institute of Biotechnology, Amity University Haryana, Gurgaon, India
| | - Amresh Prakash
- Amity Institute of Integrative Sciences and Health (AIISH), Amity University Haryana, Gurgaon, India
- *Correspondence: Vijay Kumar Amresh Prakash
| | - Vijay Kumar
- Amity Institute of Neuropsychology & Neurosciences (AINN), Amity University, Noida, India
- *Correspondence: Vijay Kumar Amresh Prakash
| |
Collapse
|
9
|
Ali N, Prasad K, AlAsmari AF, Alharbi M, Rashid S, Kumar V. Genomics-guided targeting of stress granule proteins G3BP1/2 to inhibit SARS-CoV-2 propagation. Int J Biol Macromol 2021; 190:636-648. [PMID: 34517025 PMCID: PMC8431879 DOI: 10.1016/j.ijbiomac.2021.09.018] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2021] [Revised: 08/07/2021] [Accepted: 09/03/2021] [Indexed: 12/27/2022]
Abstract
SARS-CoV-2 nucleocapsid (N) protein undergoes RNA-induced phase separation (LLPS) and sequesters the host key stress granule (SG) proteins, Ras-GTPase-activating protein SH3-domain-binding protein 1 and 2 (G3BP1 and G3BP2) to inhibit SG formation. This will allow viral packaging and propagation in host cells. Based on a genomic-guided meta-analysis, here we identify upstream regulatory elements modulating the expression of G3BP1 and G3BP2 (collectively called G3BP1/2). Using this strategy, we have identified FOXA1, YY1, SYK, E2F-1, and TGFBR2 as activators and SIN3A, SRF, and AKT-1 as repressors of G3BP1/2 genes. Panels of the activators and repressors were then used to identify drugs that change their gene expression signatures. Two drugs, imatinib, and decitabine have been identified as putative modulators of G3BP1/2 genes and their regulators, suggesting their role as COVID-19 mitigation agents. Molecular docking analysis suggests that both drugs bind to G3BP1/2 with a much higher affinity than the SARS-CoV-2 N protein. This study reports imatinib and decitabine as candidate drugs against N protein and G3BP1/2 protein.
Collapse
Affiliation(s)
- Nemat Ali
- Department of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia
| | - Kartikay Prasad
- Amity Institute of Neuropsychology & Neurosciences, Amity University, Noida, UP 201303, India
| | - Abdullah F AlAsmari
- Department of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia
| | - Metab Alharbi
- Department of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia
| | - Summya Rashid
- Department of Pharmacology & Toxicology, College of Pharmacy Girls Section, Prince Sattam Bin AbdulAziz University, Al kharj 11942, Saudi Arabia
| | - Vijay Kumar
- Amity Institute of Neuropsychology & Neurosciences, Amity University, Noida, UP 201303, India.
| |
Collapse
|
10
|
Deng X, Wang S, Han Z, Gong W, Liu Y, Li C. Dynamics of binding interactions of TDP-43 and RNA: An equally weighted multiscale elastic network model study. Proteins 2021; 90:589-600. [PMID: 34599611 DOI: 10.1002/prot.26255] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2021] [Revised: 09/15/2021] [Accepted: 09/21/2021] [Indexed: 01/03/2023]
Abstract
Transactive response DNA binding protein 43 (TDP-43), an alternative-splicing regulator, can specifically bind long UG-rich RNAs, associated with a range of neurodegenerative diseases. Upon binding RNA, TDP-43 undergoes a large conformational change with two RNA recognition motifs (RRMs) connected by a long linker rearranged, strengthening the binding affinity of TDP-43 with RNA. We extend the equally weighted multiscale elastic network model (ewmENM), including its Gaussian network model (ewmGNM) and Anisotropic network model (ewmANM), with the multiscale effect of interactions considered, to the characterization of the dynamics of binding interactions of TDP-43 and RNA. The results reveal upon RNA binding a loss of flexibility occurs to TDP-43's loop3 segments rich in positively charged residues and C-terminal of high flexibility, suggesting their anchoring RNA, induced fit and conformational adjustment roles in recognizing RNA. Additionally, based on movement coupling analyses, it is found that RNA binding strengthens the interactions among intra-RRM β-sheets and between RRMs partially through the linker's mediating role, which stabilizes RNA binding interface, facilitating RNA binding efficiency. In addition, utilizing our proposed thermodynamic cycle method combined with ewmGNM, we identify the key residues for RNA binding whose perturbations induce a large change in binding free energy. We identify not only the residues important for specific binding, but also the ones critical for the conformational rearrangement between RRMs. Furthermore, molecular dynamics simulations are also performed to validate and further interpret the ENM-based results. The study demonstrates a useful avenue to utilize ewmENM to investigate the protein-RNA interaction dynamics characteristics.
Collapse
Affiliation(s)
- Xueqing Deng
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing, China
| | - Shihao Wang
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing, China
| | - Zhongjie Han
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing, China
| | - Weikang Gong
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing, China
| | - Yang Liu
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing, China
| | - Chunhua Li
- Faculty of Environmental and Life Sciences, Beijing University of Technology, Beijing, China
| |
Collapse
|
11
|
Liu W, Li C, Shan J, Wang Y, Chen G. Insights into the aggregation mechanism of RNA recognition motif domains in TDP-43: a theoretical exploration. ROYAL SOCIETY OPEN SCIENCE 2021; 8:210160. [PMID: 34457335 PMCID: PMC8371369 DOI: 10.1098/rsos.210160] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/04/2021] [Accepted: 07/19/2021] [Indexed: 05/10/2023]
Abstract
The transactive response DNA-binding protein 43 (TDP-43) is associated with several diseases such as amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD) due to pathogenic aggregations. In this work, we examined the dimer, tetramer and hexamer models built from the RRM domains of TDP-43 using molecular dynamics simulations in combination with the protein-protein docking. Our results showed that the formations of the dimer models are mainly achieved by the interactions of the RRM1 domains. The parallel β-sheet layers between the RRM1 domains provide most of the binding sites in these oligomer models, and thus play an important role in the aggregation process. The approaching of the parallel β-sheet layers from small oligomer models gradually expand to large ones through the allosteric communication between the α1/α2 helices of the RRM1 domains, which maintains the binding affinities and interactions in the larger oligomer models. Using the repeatable-superimposing method based on the tetramer models, we proposed a new aggregation mechanism of RRM domains in TDP-43, which could well characterize the formation of the large aggregation models with the repeated, helical and rope-like structures. These new insights help to understand the amyloid-like aggregation phenomena of TDP-43 protein in ALS and FTLD diseases.
Collapse
Affiliation(s)
- Wei Liu
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Chaoqun Li
- Hebei Key Laboratory of Heterocyclic Compounds, College of Chemistry, Chemical Engineering and Materials, Handan University, Handan 056005, Hebei Province, People's Republic of China
| | - Jiankai Shan
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Yan Wang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| | - Guangju Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People's Republic of China
| |
Collapse
|
12
|
Jairajpuri DS, Hussain A, Nasreen K, Mohammad T, Anjum F, Tabish Rehman M, Mustafa Hasan G, Alajmi MF, Imtaiyaz Hassan M. Identification of natural compounds as potent inhibitors of SARS-CoV-2 main protease using combined docking and molecular dynamics simulations. Saudi J Biol Sci 2021; 28:2423-2431. [PMID: 33526965 PMCID: PMC7839507 DOI: 10.1016/j.sjbs.2021.01.040] [Citation(s) in RCA: 36] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2020] [Revised: 01/17/2021] [Accepted: 01/19/2021] [Indexed: 12/12/2022] Open
Abstract
Coronavirus disease 2019 (COVID-19) has emerged from China and globally affected the entire population through the human-to-human transmission of a newly emerged virus called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The genome of SARS-CoV-2 encodes several proteins that are essential for multiplication and pathogenesis. The main protease (Mpro or 3CLpro) of SARS-CoV-2 plays a central role in its pathogenesis and thus is considered as an attractive drug target for the drug design and development of small-molecule inhibitors. We have employed an extensive structure-based high-throughput virtual screening to discover potential natural compounds from the ZINC database which could inhibit the Mpro of SARS-CoV-2. Initially, the hits were selected on the basis of their physicochemical and drug-like properties. Subsequently, the PAINS filter, estimation of binding affinities using molecular docking, and interaction analyses were performed to find safe and potential inhibitors of SARS-CoV-2 Mpro. We have identified ZINC02123811 (1-(3-(2,5,9-trimethyl-7-oxo-3-phenyl-7H-furo[3,2-g]chromen-6-yl)propanoyl)piperidine-4-carboxamide), a natural compound bearing appreciable affinity, efficiency, and specificity towards the binding pocket of SARS-CoV-2 Mpro. The identified compound showed a set of drug-like properties and preferentially binds to the active site of SARS-CoV-2 Mpro. All-atom molecular dynamics (MD) simulations were performed to evaluate the conformational dynamics, stability and interaction mechanism of Mpro with ZINC02123811. MD simulation results indicated that Mpro with ZINC02123811 forms a stable complex throughout the trajectory of 100 ns. These findings suggest that ZINC02123811 may be further exploited as a promising scaffold for the development of potential inhibitors of SARS-CoV-2 Mpro to address COVID-19.
Collapse
Affiliation(s)
- Deeba Shamim Jairajpuri
- Department of Medical Biochemistry, College of Medicine and Medical Sciences, Arabian Gulf University, P.O. Box 22971, Manama, Bahrain
| | - Afzal Hussain
- Department of Pharmacognosy, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia
| | - Khalida Nasreen
- Centre for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia, Jamia Nagar, New Delhi 110025, India
| | - Taj Mohammad
- Centre for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia, Jamia Nagar, New Delhi 110025, India
| | - Farah Anjum
- Department of Clinical Laboratory Sciences, College of Applied Medical Sciences, Taif University, Taif, Saudi Arabia
| | - Md. Tabish Rehman
- Department of Pharmacognosy, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia
| | - Gulam Mustafa Hasan
- Department of Biochemistry, College of Medicine, Prince Sattam Bin Abdulaziz University, P.O. Box 173 Al-Kharj, 11942, Saudi Arabia
| | - Mohamed F. Alajmi
- Department of Pharmacognosy, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia
| | - Md. Imtaiyaz Hassan
- Centre for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia, Jamia Nagar, New Delhi 110025, India
| |
Collapse
|
13
|
Prasad K, Ahamad S, Kanipakam H, Gupta D, Kumar V. Simultaneous Inhibition of SARS-CoV-2 Entry Pathways by Cyclosporine. ACS Chem Neurosci 2021; 12:930-944. [PMID: 33606519 DOI: 10.1021/acschemneuro.1c00019] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
The COVID-19 pandemic caused by SARS-CoV-2 represents a global public health emergency. The entry of SARS-CoV-2 into host cells requires the activation of its spike protein by host cell proteases. The serine protease, TMPRSS2, and cysteine proteases, Cathepsins B/L, activate spike protein and enable SARS-CoV-2 entry to the host cell through two completely different and independent pathways. Therefore, inhibiting either TMPRSS2 or cathepsin B/L may not sufficiently block the virus entry. We here hypothesized that simultaneous targeting of both the entry pathways would be more efficient to block the virus entry rather than targeting the entry pathways individually. To this end, we utilized the network-based drug repurposing analyses to identify the possible common drugs that can target both the entry pathways. This study, for the first time, reports the molecules like cyclosporine, calcitriol, and estradiol as candidate drugs with the binding ability to the host proteases, TMPRSS2, and cathepsin B/L. Next, we analyzed drug-gene and gene-gene interaction networks using 332 human targets of SARS-CoV-2 proteins. The network results indicate that, out of 332 human proteins, cyclosporine interacts with 216 (65%) proteins. Furthermore, we performed molecular docking and all-atom molecular dynamics (MD) simulations to explore the binding of drug with TMPRSS2 and cathepsin L. The molecular docking and MD simulation results showed strong and stable binding of cyclosporine A (CsA) with TMPRSS2 and CTSL genes. The above results indicate cyclosporine as a potential drug molecule, as apart from interacting with SARS-CoV-2 entry receptors, it also interacts with most of SARS-CoV-2 target host genes; thus it could potentially interfere with functions of SARS-CoV-2 proteins in human cells. We here also suggest that these antiviral drugs alone or in combination can simultaneously target both the entry pathways and thus can be considered as a potential treatment option for COVID-19.
Collapse
Affiliation(s)
- Kartikay Prasad
- Amity Institute of Neuropsychology & Neurosciences, Amity University, Noida, UP 201303, India
| | - Shahzaib Ahamad
- Translational Bioinformatics Group, International Centre for Genetic Engineering and Biotechnology (ICGEB), Aruna Asaf Ali Marg, New Delhi 110067, India
| | - Hema Kanipakam
- Translational Bioinformatics Group, International Centre for Genetic Engineering and Biotechnology (ICGEB), Aruna Asaf Ali Marg, New Delhi 110067, India
| | - Dinesh Gupta
- Translational Bioinformatics Group, International Centre for Genetic Engineering and Biotechnology (ICGEB), Aruna Asaf Ali Marg, New Delhi 110067, India
| | - Vijay Kumar
- Amity Institute of Neuropsychology & Neurosciences, Amity University, Noida, UP 201303, India
| |
Collapse
|
14
|
Tethering-induced destabilization and ATP-binding for tandem RRM domains of ALS-causing TDP-43 and hnRNPA1. Sci Rep 2021; 11:1034. [PMID: 33441818 PMCID: PMC7806782 DOI: 10.1038/s41598-020-80524-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2020] [Accepted: 12/21/2020] [Indexed: 12/12/2022] Open
Abstract
TDP-43 and hnRNPA1 contain tandemly-tethered RNA-recognition-motif (RRM) domains, which not only functionally bind an array of nucleic acids, but also participate in aggregation/fibrillation, a pathological hallmark of various human diseases including amyotrophic lateral sclerosis (ALS), frontotemporal dementia (FTD), alzheimer's disease (AD) and Multisystem proteinopathy (MSP). Here, by DSF, NMR and MD simulations we systematically characterized stability, ATP-binding and conformational dynamics of TDP-43 and hnRNPA1 RRM domains in both tethered and isolated forms. The results reveal three key findings: (1) upon tethering TDP-43 RRM domains become dramatically coupled and destabilized with Tm reduced to only 49 °C. (2) ATP specifically binds TDP-43 and hnRNPA1 RRM domains, in which ATP occupies the similar pockets within the conserved nucleic-acid-binding surfaces, with the affinity slightly higher to the tethered than isolated forms. (3) MD simulations indicate that the tethered RRM domains of TDP-43 and hnRNPA1 have higher conformational dynamics than the isolated forms. Two RRM domains become coupled as shown by NMR characterization and analysis of inter-domain correlation motions. The study explains the long-standing puzzle that the tethered TDP-43 RRM1–RRM2 is particularly prone to aggregation/fibrillation, and underscores the general role of ATP in inhibiting aggregation/fibrillation of RRM-containing proteins. The results also rationalize the observation that the risk of aggregation-causing diseases increases with aging.
Collapse
|
15
|
Dahiya R, Mohammad T, Alajmi MF, Rehman MT, Hasan GM, Hussain A, Hassan MI. Insights into the Conserved Regulatory Mechanisms of Human and Yeast Aging. Biomolecules 2020; 10:E882. [PMID: 32526825 PMCID: PMC7355435 DOI: 10.3390/biom10060882] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2020] [Revised: 05/22/2020] [Accepted: 05/27/2020] [Indexed: 12/12/2022] Open
Abstract
Aging represents a significant biological process having strong associations with cancer, diabetes, and neurodegenerative and cardiovascular disorders, which leads to progressive loss of cellular functions and viability. Astonishingly, age-related disorders share several genetic and molecular mechanisms with the normal aging process. Over the last three decades, budding yeast Saccharomyces cerevisiae has emerged as a powerful yet simple model organism for aging research. Genetic approaches using yeast RLS have led to the identification of hundreds of genes impacting lifespan in higher eukaryotes. Numerous interventions to extend yeast lifespan showed an analogous outcome in multi-cellular eukaryotes like fruit flies, nematodes, rodents, and humans. We collected and analyzed a multitude of observations from published literature and provide the contribution of yeast in the understanding of aging hallmarks most applicable to humans. Here, we discuss key pathways and molecular mechanisms that underpin the evolutionarily conserved aging process and summarize the current understanding and clinical applicability of its trajectories. Gathering critical information on aging biology would pave the way for future investigation targeted at the discovery of aging interventions.
Collapse
Affiliation(s)
- Rashmi Dahiya
- Centre for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia, New Delhi 110025, India;
| | - Taj Mohammad
- Centre for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia, New Delhi 110025, India;
| | - Mohamed F. Alajmi
- Department of Pharmacognosy, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia; (M.F.A.); (M.T.R.); (A.H.)
| | - Md. Tabish Rehman
- Department of Pharmacognosy, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia; (M.F.A.); (M.T.R.); (A.H.)
| | - Gulam Mustafa Hasan
- Department of Biochemistry, College of Medicine, Prince Sattam Bin Abdulaziz University, P.O. Box 173, Al-Kharj 11942, Saudi Arabia;
| | - Afzal Hussain
- Department of Pharmacognosy, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia; (M.F.A.); (M.T.R.); (A.H.)
| | - Md. Imtaiyaz Hassan
- Centre for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia, New Delhi 110025, India;
| |
Collapse
|
16
|
Prakash A, Kumar V, Banerjee A, Lynn AM, Prasad R. Structural heterogeneity in RNA recognition motif 2 (RRM2) of TAR DNA-binding protein 43 (TDP-43): clue to amyotrophic lateral sclerosis. J Biomol Struct Dyn 2020; 39:357-367. [DOI: 10.1080/07391102.2020.1714481] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Amresh Prakash
- Amity Institute of Integrative Sciences and Health, Amity University, Gurgaon, India
| | - Vijay Kumar
- Amity Institute of Neuropsychology & Neurosciences, Amity University, Noida, India
| | - Atanu Banerjee
- Amity Institute of Integrative Sciences and Health, Amity University, Gurgaon, India
| | - Andrew M. Lynn
- School of Computational & Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Rajendra Prasad
- Amity Institute of Biotechnology, Amity University, Gurgaon, India
| |
Collapse
|
17
|
Singh R, Meena NK, Das T, Sharma RD, Prakash A, Lynn AM. Delineating the conformational dynamics of intermediate structures on the unfolding pathway of β-lactoglobulin in aqueous urea and dimethyl sulfoxide. J Biomol Struct Dyn 2019; 38:5027-5036. [PMID: 31744390 DOI: 10.1080/07391102.2019.1695669] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
The funnel shaped energy landscape model of the protein folding suggests that progression of folding proceeds through multiple pathways, having the multiple intermediates which leads to multidimensional free-energy surface. Herein, we applied all-atom MD simulation to conduct a comparative study on the structure of β-lactoglobulin (β-LgA) in aqueous mixture of 8 M urea and 8 M dimethyl sulfoxide (DMSO), at different temperatures. The cumulative results of multiple simulations suggest a common unfolding pathway of β-LgA, occurred through the stable and meta-stable intermediates (I), in both urea and DMSO. However, the free-energy landscape (FEL) analyses show that the structural transitions of I-states are energetically different. In urea, FEL shows distinct ensemble of intermediates, I1 and I2, separated by the energy barrier of ∼3.0 kcal mol-1. Similarly, we find the population of two distinct I1 and I2 states in DMSO, however, the I1 appeared transiently around ∼30-35 ns and is short-lived. But, the I2 ensemble is observed structurally compact and long-lived (∼50-150 ns) as compared to unfolding in urea. Furthermore, the I1 and I2 are separated through a high energy barrier of ∼6.0 kcal mol-1. Thus, our results provide the structural insights of intermediates which essentially bear the signature of a different unfolding pathway of β-LgA in urea and DMSO.Abbreviationsβ-LgAβ-lactoglobulinDMSOdimethyl sulfoxideFELfree-energy landscapeGdmClguanidinium chlorideIintermediate stateMGmolten globule statePMEparticle mesh EwaldQfraction of native contactsRMSDroot mean square deviationRMSFroot mean square fluctuationRgradius of gyrationSASAsolvent Accessible Surface AreascSASAthe side chain SASATrptryptophanCommunicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Ruhar Singh
- School of Computational & Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Naveen Kumar Meena
- School of Computational & Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Trishala Das
- School of Computational & Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Ravi Datta Sharma
- Amity Institute of Integrative Sciences and Health, Amity University, Haryana, India
| | - Amresh Prakash
- Amity Institute of Integrative Sciences and Health, Amity University, Haryana, India
| | - Andrew M Lynn
- School of Computational & Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| |
Collapse
|
18
|
Sharma M, Bulusu G, Mitra A. Unfolding Transitions of Peripheral Subunit Binding Domains Show Cooperative Behavior. J Phys Chem B 2019; 123:3441-3451. [PMID: 30958002 DOI: 10.1021/acs.jpcb.9b01114] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Characterization of native, intermediate, and denatured states is crucial for understanding the factors influencing the stability of proteins. We have carried out molecular dynamics simulations to study the unfolding of three peripheral subunit binding domains (PSBDs): E. coli BBL, Bacillus stearothermophilus E3BD, and human hbSBD, at three different temperatures: 300, 330, and 400 K, and in the presence of two solvents: water and 5 M guanidinium hydrochloride (GndCl) solution. These proteins share similar folds, with two parallel helices, maintained via a hydrophobic core comprising residues from their interconnecting loop. BBL is more sensitive to thermal and chemical denaturation in comparison to hbSBD, and E3BD is the most stable of all of the three proteins. The effect of temperature on the stability of these proteins is more pronounced in "water-only" simulations compared to that in the presence of guanidium hydrochloride in high concentrations. Our results show cooperative unfolding transitions of these proteins, which are triggered by an initial melting of the C-terminal helix H2. The consequent loss of interhelical interactions or native contacts, as observed, leads to the subsequent melting of the N-terminal helix H1.
Collapse
Affiliation(s)
- Monika Sharma
- Department of Chemical Sciences , Indian Institute of Science Education and Research (IISER) , Sector 81, Knowledge City , SAS Nagar 140306 , Punjab , India
| | - Gopalakrishnan Bulusu
- TCS Innovation Labs - Hyderabad (Life Sciences Division), Tata Consultancy Services Limited , Hyderabad 500081 , India.,Center of Computational Natural Sciences and Bioinformatics (CCNSB) , International Institute of Information Technology (IIIT) , Hyderabad 500032 , India
| | - Abhijit Mitra
- Center of Computational Natural Sciences and Bioinformatics (CCNSB) , International Institute of Information Technology (IIIT) , Hyderabad 500032 , India
| |
Collapse
|
19
|
Babinchak WM, Li Z, Buck M. Keys to Amyloid City: Computation and NMR Reveal Potential TDP-43 ALS Intermediates. Biophys J 2018; 115:1625-1627. [PMID: 30290957 PMCID: PMC6225042 DOI: 10.1016/j.bpj.2018.09.012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2018] [Accepted: 09/13/2018] [Indexed: 11/26/2022] Open
Affiliation(s)
- W Michael Babinchak
- Department of Physiology and Biophysics, Case Western Reserve University, Cleveland, Ohio
| | - Zhenlu Li
- Department of Physiology and Biophysics, Case Western Reserve University, Cleveland, Ohio
| | - Matthias Buck
- Department of Physiology and Biophysics, Case Western Reserve University, Cleveland, Ohio.
| |
Collapse
|
20
|
Prakash A, Kumar V, Meena NK, Lynn AM. Elucidation of the structural stability and dynamics of heterogeneous intermediate ensembles in unfolding pathway of the N-terminal domain of TDP-43. RSC Adv 2018; 8:19835-19845. [PMID: 35548664 PMCID: PMC9088055 DOI: 10.1039/c8ra03368d] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2018] [Accepted: 05/14/2018] [Indexed: 11/21/2022] Open
Abstract
The N-terminal domain of the RNA binding protein TDP-43 (NTD) is essential to both physiology and proteinopathy; however, elucidation of its folding/unfolding still remains a major quest. In this study, we have investigated the biophysical behavior of intermediate ensembles employing all-atom molecular dynamics simulations in 8 M urea accelerated with high temperatures to achieve unfolded states in a confined computation time. The cumulative results of the 2.75 μs simulations show that unfolding of the NTD at 350 K evolves through different stable and meta-stable intermediate states. The free-energy landscape reveals two meta-stable intermediates (IN and IU) stabilized by non-native interactions, which are largely hydrophilic and highly energetically frustrated. A single buried tryptophan residue, W80, undergoes solvent exposure to different extents during unfolding; this suggests a structurally heterogeneous population of intermediate ensembles. Furthermore, the structure properties of the IN state show a resemblance to the molten globule (MG) state with most of the secondary structures intact. The unfolding of the NTD is initiated by the loss of β-strands, and the unfolded (U) states exhibit a population of non-native α-helices. These non-native unfolded intermediate ensembles may mediate protein oligomerization, leading to the formation of pathological, irreversible aggregates, characteristics of disease pathogenesis.
Collapse
Affiliation(s)
- Amresh Prakash
- School of Computational & Integrative Sciences, Jawaharlal Nehru University New Delhi-110067 India
| | - Vijay Kumar
- Centre for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia Jamia Nagar New Delhi-110025 India
| | - Naveen Kumar Meena
- School of Computational & Integrative Sciences, Jawaharlal Nehru University New Delhi-110067 India
| | - Andrew M Lynn
- School of Computational & Integrative Sciences, Jawaharlal Nehru University New Delhi-110067 India
| |
Collapse
|