Phosphorylation-induced changes in the energetic frustration in human Tank binding kinase 1

Computational insights into dynamics and conformational stability of N-acetylmannosamine kinase mutations

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.

Identifying potential neuroprotective polyphenols targeting endoplasmic reticulum stress through an in silico approach

The endoplasmic reticulum (ER) is essential in many cellular processes, including protein processing, lipid metabolism, and calcium storage. Dysregulation of ER function has been linked with neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, etc. The primary pathological alteration explicated in the diseases is the accumulation of misfolded proteins in the neuronal cells. ER stress-associated activation of PERK-mediated pro-apoptotic cell death leads to neurodegeneration. In this study, we have primarily screened the potential polyphenols evidenced for neuroprotective activity. The 24 polyphenols were selected to explore their binding affinity towards various proteins of ER cascade such as pPERK (phospho-PERK), EIF2 (Eukaryotic Initiation Factor 2), and ATF4 (Activating Transcription Factor 4). On the basis of binding affinity, four phytopolyphenols were further selected for in-silico ADMET and molecular dynamic simulation. Among them curcumin found to be the most promising and serve as a potential hit against all three targets of ER cascade. The selected proteins' active site has demonstrated high stability of curcumin binding according to molecular dynamics findings. Though curcumin exhibited a significant hit in interaction with targets but needs to be further improved in drug-ability criteria. Thus, seventy derivatives of curcumin scaffold (from the published literature) were also screened with improve in druggability criteria, which showed good interaction with unfolded protein response related targets. The new scaffolds serve considerable potential to be developed as novel polyphenolic lead for neurodegenerative disorders.Communicated by Ramaswamy H. Sarma.

In silico evaluation of the inhibitory potential of nucleocapsid inhibitors of SARS-CoV-2: a binding and energetic perspective

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.

Glutaredoxin proteins from E. coli isoforms were compared in terms of energy frustration

Glutaredoxin (GRXs) protein plays a vital role inside the cell, including redox control of transcription to the cell's antioxidant defense, apoptosis, and cellular differentiation regulation. In this study, we have investigated the energy landscape and characterized the pattern of local frustration in different forms and states of the GRX protein ofE. coli.Analysis was done on the conformational alterations, significant changes in the frustration pattern, and different GRXs such as GRX-II, GRX-III, GRX-II-GSH, and GRX-III-GSH complex. We have found the practice of frustration, and structure was quite similar in the same isoform having different states of protein; however, a significant difference was observed between different isoforms. Moreover, oxidation of GRX-I introduced an extra α-helix increasing the destabilizing interactions within the protein. The study of frustrated contacts on oxidized and reduced GRX and with bound and unbound Glutathione indicates its potential application in activating and regulating the behavior of GRXs.

Energetic and frustration analysis of SARS-CoV-2 nucleocapsid protein mutations

The ongoing COVID-19 spreads worldwide with the ability to evolve in diverse human populations. The nucleocapsid (N) protein is one of the mutational hotspots in the SARS-CoV-2 genome. The N protein is an abundant RNA-binding protein critical for viral genome packaging. It comprises two large domains including the N-terminal domain (NTD) and the C-terminal domain (CTD) linked by the centrally located linker region. Mutations in N protein have been reported to increase the severity of disease by modulating viral transmissibility, replication efficiency as well as virulence properties of the virus in different parts of the world. To study the effect of N protein missense mutations on protein stability, function, and pathogenicity, we analyzed 228 mutations from each domain of N protein. Further, we have studied the effect of mutations on local residual frustration changes in N protein. Out of 228 mutations, 11 mutations were predicted to be deleterious and destabilized. Among these mutations, R32C, R191C, and R203 M mutations fall into disordered regions and show significant change in frustration state. Overall, this work reveals that by altering the energetics and residual frustration, N protein mutations might affect the stability, function, and pathogenicity of the SARS-CoV-2.

Frustration analysis of TBK1 missense mutations reported in ALS/FTD and cancer patients

Tank-binding kinase 1 (TBK1) is a multifunctional kinase having essential roles in cellular processes, autophagy/mitophagy, and selective clearance of damaged proteins. More than 90 mutations in the TBK1 gene are linked with multiple cancer types, amyotrophic lateral sclerosis (ALS), and frontotemporal dementia (FTD). Some of these missense mutations disrupt the abilities of TBK1 to dimerize, associate with the mitophagy receptor optineurin (OPTN), autoactivate, or catalyze phosphorylation. Some mutations may cause severe dysregulation of the pathway, while others induce a limited disruption. Here, we have studied those mutations reported in cancer, ALS and FTD, and subsequently investigated the effect of missense mutations on the structure and function of TBK1 for localized residual frustration change. Out of 33 ALS/FTD causing mutations and 28 oncogenic mutations, 10 mutations and 12 oncogenic mutations showed significant change in the residual frustration. The local frustration plays an important role in the conformation of protein structure in active and inactive kinases. Our analysis reports the change in residual frustration state, conformational change and effect on active and inactive TBK1 function due to ALS/FTD causing and oncogenic missense mutations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03240-0.

Structure-guided design and development of cyclin-dependent kinase 4/6 inhibitors: A review on therapeutic implications

Cyclin-dependent kinase 6 (EC 2.7.11.22) play significant roles in numerous biological processes and triggers cell cycle events. CDK6 controlled the transcriptional regulation. A dysregulated function of CDK6 is linked with the development of progression of multiple tumor types. Thus, it is considered as an effective drug target for cancer therapy. Based on the direct roles of CDK4/6 in tumor development, numerous inhibitors developed as promising anti-cancer agents. CDK4/6 inhibitors regulate the G1 to S transition by preventing Rb phosphorylation and E2F liberation, showing potent anti-cancer activity in several tumors, including HR+/HER2- breast cancer. CDK4/6 inhibitors such as abemaciclib, palbociclib, and ribociclib, control cell cycle, provoke cell senescence, and induces tumor cell disturbance in pre-clinical studies. Here, we discuss the roles of CDK6 in cancer along with the present status of CDK4/6 inhibitors in cancer therapy. We further discussed, how structural features of CDK4/6 could be implicated in the design and development of potential anti-cancer agents. In addition, the therapeutic potential and limitations of available CDK4/6 inhibitors are described in detail. Recent pre-clinical and clinical information for CDK4/6 inhibitors are highlighted. In addition, combination of CDK4/6 inhibitors with other drugs for the therapeutic management of cancer are discussed.

Delineating the Aggregation-Prone Hotspot Regions (Peptides) in the Human Cu/Zn Superoxide Dismutase 1

Amyotrophic lateral sclerosis (ALS) is a fatal, incurable neurodegenerative disease described by progressive degeneration of motor neurons. The most common familial form of ALS (fALS) has been associated with mutations in the Cu/Zn superoxide dismutase (SOD1) gene. Mutation-induced misfolding and aggregation of SOD1 is often found in ALS patients. In this work, we probe the aggregation properties of peptides derived from the SOD1. To examine the source of SOD1 aggregation, we have employed a computational algorithm to identify four peptides from the SOD1 protein sequence that aggregates into a fibril. Aided by computational algorithms, we identified four peptides likely involved in SOD1 fibrillization. These four aggregation-prone peptides were 14VQGIINFE21, 30KVWGSIKGL38, 101DSVISLS107, and 147GVIGIAQ153. In addition, the formation of fibril propensities from the identified peptides was investigated through different biophysical techniques. The atomic structures of two fibril-forming peptides from the C-terminal SOD1 showed that the steric zippers formed by 101DSVISLS107 and 147GVIGIAQ153 vary in their arrangement. We also discovered that fALS mutations in the peptide 147GVIGIAQ153 increased the fibril-forming propensity and altered the steric zipper's packing. Thus, our results suggested that the C-terminal peptides of SOD1 have a central role in amyloid formation and might be involved in forming the structural core of SOD1 aggregation observed in vivo.

Molecular interactions between C9ORF72 and SMCR8: A local energetic frustration perspective

The hexanucleotide repeat expansion in C9orf72 represents a major cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). C9orf72, together with SMCR8 and WDR41, can form a stable complex that regulates autophagy and membrane trafficking. Very recently, the cryo-EM structure of C9orf72-SMCR8-WDR41 helps in understanding the structure-function relationship of C9orf72. This protein complex is indispensable to several cellular processes and is strongly linked to familial ALS and FTD. Understanding the molecular basis of the C9orf72-SMCR8 protein-protein interaction is thus important to comprehend their function. To establish a basis for understanding the relationships between sequence, structure, and function of the C9orf72, this study reports a local frustration analysis on the C9orf72-SMCR8 complex structure. An analysis of local frustration profiles indicated that (1) the structural domains in C9orf72 are minimally-frustrated and relatively conserved, (2) high frustration patches on the protein-protein interface (3) increased frustration in the C-terminal helices involved in the dimerization of C9orf72 structures.

A review on the role of TANK-binding kinase 1 signaling in cancer

TANK-binding kinase 1 (TBK1) regulates various biological processes including, NF-κB signaling, immune response, autophagy, cell division, Ras-mediated oncogenesis, and AKT pro-survival signaling. Enhanced TBK1 activity is associated with autoimmune diseases and cancer, suggesting its role in therapeutic targeting of interferonopathies. In addition, dysregulation of TBK1 activity promotes several inflammatory disorders and oncogenesis. Structural and biochemical study reports provide the molecular process of TBK1 activation and recap the substrate selection about TBK1. This review summarizes recent findings on the molecular mechanisms by which TBK1 is involved in cancer signaling. The IKK-ε and TBK1 are together associated with inflammatory diseases by inducing type I IFNs. Furthermore, TBK1 signaling regulates radiation-induced epithelial-mesenchymal transition by controlling phosphorylation of GSK-3β and expression of Zinc finger E-box-binding homeobox 1, suggesting, TBK1 could be targeted for radiotherapy-induced metastasis therapy. Despite a considerable increase in the list of TBK1 inhibitors, only a few has potential to control cancer. Among them, a compound BX795 is considered a potent and selective inhibitor of TBK1. We discussed the therapeutic potential of small-molecule inhibitors of TBK1, particularly those with high selectivity, which will enable further exploration in the therapeutic management of cancer and inflammatory diseases.

Insights into the biased activity of dextromethorphan and haloperidol towards SARS-CoV-2 NSP6: in silico binding mechanistic analysis

Abstract

The outbreak of novel coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) virus continually led to infect a large population worldwide. SARS-CoV-2 utilizes its NSP6 and Orf9c proteins to interact with sigma receptors that are implicated in lipid remodeling and ER stress response, to infect cells. The drugs targeting the sigma receptors, sigma-1 and sigma-2, have emerged as effective candidates to reduce viral infectivity, and some of them are in clinical trials against COVID-19. The antipsychotic drug, haloperidol, exerts remarkable antiviral activity, but, at the same time, the sigma-1 benzomorphan agonist, dextromethorphan, showed pro-viral activity. To explore the potential mechanisms of biased binding and activity of the two drugs, haloperidol and dextromethorphan towards NSP6, we herein utilized molecular docking–based molecular dynamics simulation studies. Our extensive analysis of the protein-drug interactions, structural and conformational dynamics, residual frustrations, and molecular switches of NSP6-drug complexes indicates that dextromethorphan binding leads to structural destabilization and increase in conformational dynamics and energetic frustrations. On the other hand, the strong binding of haloperidol leads to minimal structural and dynamical perturbations to NSP6. Thus, the structural insights of stronger binding affinity and favorable molecular interactions of haloperidol towards viral NSP6 suggests that haloperidol can be potentially explored as a candidate drug against COVID-19.

Key messages

•Inhibitors of sigma receptors are considered as potent drugs against COVID-19.

•Antipsychotic drug, haloperidol, binds strongly to NSP6 and induces the minimal changes in structure and dynamics of NSP6.

•Dextromethorphan, agonist of sigma receptors, binding leads to overall destabilization of NSP6.

•These two drugs bind with NSP6 differently and also induce differences in the structural and conformational changes that explain their different mechanisms of action.

•Haloperidol can be explored as a candidate drug against COVID-19.

Electronic supplementary material

The online version of this article (10.1007/s00109-020-01980-1) contains supplementary material, which is available to authorized users.

Elucidation of the structural stability and dynamics of heterogeneous intermediate ensembles in unfolding pathway of the N-terminal domain of TDP-43

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.