Can Duvelisib and Eganelisib work for both cancer and COVID-19? Molecular-level insights from MD simulations and enhanced samplings.

SARS-CoV-2 has caused severe illness and anxiety worldwide, evolving into more dreadful variants capable of evading the host's immunity. Cytokine storms, led by PI3Kγ, are common in cancer and SARS-CoV-2. Naturally, there is a yearning to see whether any drugs could alleviate cytokine storms for both. Upon investigation, we identified two anticancer drugs, Duvelisib and Eganelisib, that could also work against SARS-CoV-2. This report is the first to decipher their synergic therapeutic effectiveness against COVID-19 and cancer with molecular insights from atomistic simulations. In addition to PI3Kγ, these drugs exhibit specificity for the main protease among all SARS-CoV-2 targets, with significant negative binding free energies and small time-dependent conformational changes of the complexes. Complexation makes active sites and secondary structures highly mechanically stiff, with barely any deformation. Replica simulations estimated large pulling forces in enhanced sampling to dissociate the drugs from Mpro's active site. Furthermore, the radial distribution function (RDF) demonstrated that the therapeutic molecules were closest to the His41 and Cys145 catalytic dyad residues. Finally, analyses implied Duvelisib and Eganelisib as promising dual-purposed anti-COVID and anticancer drugs, potentially targeting Mpro and PI3Kγ to stop virus replication and cytokine storms concomitantly. We also distinguished hotspot residues imparting significant interactions.

In silico evaluation of S-adenosyl-L-homocysteine analogs as inhibitors of nsp14-viral cap N7 methyltranferase and PLpro of SARS-CoV-2: synthesis, molecular docking, physicochemical data, ADMET and molecular dynamics simulations studies.

A series of S-adenosyl-L-homosysteine (SAH) analogs, with modification in the base and sugar moiety, have been designed, synthesized and screened as nsp14 and PLpro inhibitors of severe acute respiratory syndrome corona virus (SARS-CoV-2). The outcomes of ADMET (Adsorption, Distribution, Metabolism, Excretion, and Toxicity) studies demonstrated that the physicochemical properties of all analogs were permissible for development of these SAH analogs as antiviral agents. All molecules were screened against different SARS-CoV-2 targets using molecular docking. The docking results revealed that the SAH analogs interacted well in the active site of nsp14 protein having H-bond interactions with the amino acid residues Arg289, Val290, Asn388, Arg400, Phe401 and π-alkyl interactions with Arg289, Val290 and Phe426 of Nsp14-MTase site. These analogs also formed stable H-bonds with Leu163, Asp165, Arg167, Ser246, Gln270, Tyr274 and Asp303 residues of PLpro proteins and found to be quite stable complexes therefore behaved as probable nsp14 and PLpro inhibitors. Interestingly, analog 3 showed significant in silico activity against the nsp14 N7 methyltransferase of SARS-CoV-2. The molecular dynamics (MD) and post-MD results of analog 3 unambiguously established the higher stability of the nsp14 (N7 MTase):3 complex and also indicated its behavior as probable nsp14 inhibitor like the reference sinefungin. The docking and MD simulations studies also suggested that sinefungin did act as SARS-CoV-2 PLpro inhibitor as well. This study's findings not only underscore the efficacy of the designed SAH analogs as potent inhibitors against crucial SARS-CoV-2 proteins but also pinpoint analog 3 as a particularly promising candidate. All the study provides valuable insights, paving the way for potential advancements in antiviral drug development against SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

HighlightsSAH analogs bearing modified bases and sugar moiety have been synthesized as antivirals against SARS-CoV-2.Molecular dynamics simulation established the stability of ligand-protein complex of analog 3 with nsp14 (N7-MTase) protein of SARS-CoV-2.Molecular docking studies of SAH analogs indicated them as nsp14 N7 methyltranferase as well as the PLpro inhibitors of SARS-CoV-2.The in silico antiviral activity of SAH analogs has been found comparable to the reference drug Sinefungin.

Unmasking an Allosteric Binding Site of the Papain-like Protease in SARS-CoV-2: Molecular Dynamics Simulations of Corticosteroids.

To date, mechanistic insights into many clinical drugs against COVID-19 remain unexplored. Dexamethasone, a corticosteroid, is one of them. While treating the entire corticosteroid database, including vitamins D2 and D3, with cutting-edge computational techniques, several intriguing results are unfolded. From the top-notch candidates, dexamethasone is likely to inhibit the viral main protease (Mpro), with vitamin D3 exhibiting multitarget [Mpro, papain-like protease (PLpro), and nucleocapsid protein (N-pro)] roles and ciclesonide's dynamic flipping disinterring a cryptic allosteric site in the PLpro enzyme. The results rationalize why these drugs improve the health of COVID-19 patients. Understanding an enzyme's secret binding site is essential to understanding how the enzyme works and how to inhibit its function. Ciclesonide's allosteric inhibition could not only jeopardize PLpro's catalytic role in polyprotein processing but also make it less vulnerable to the host body's defense machinery. Hotspot residues in the identified allosteric site could be considered for effective therapeutic designs against PLpro.

Exploring antiviral potency of N-1 substituted pyrimidines against HIV-1 and other DNA/RNA viruses: Design, synthesis, characterization, ADMET analysis, docking, molecular dynamics and biological activity.

A novel series of pyrimidine derivatives, bearing modified benzimidazoles at N-1 position, has been designed, synthesized and screened as NNRTIs against HIV and as broad-spectrum antiviral agents. The molecules were screened against different HIV targets using molecular docking experiment. The docking results indicated that the molecules interacted well with the residues Lys101, Tyr181, Tyr188, Trp229, Phe227 and Tyr318 present in NNIBP of HIV-RT protein, formed quite stable complexes and, thus, behaved as probable NNRTIs. Among these compounds, 2b and 4b showed anti-HIV activity with IC50 values as 6.65 µg/mL (SI = 15.50) and 15.82 µg/mL (SI = 14.26), respectively. Similarly, compound 1a showed inhibitory property against coxsackie virus B4 and compound 3b against different viruses. Molecular dynamics simulation results unequivocally demonstrated the higher stability of the complex HIV-RT:2b than the HIV-RT:nevirapine complex. The MM/PBSA-based binding free energy (-) 114.92 kJ/mol of HIV-RT:2b complex in comparison to that of HIV-RT:nevirapine complex (-) 88.33 kJ/mol, further demonstrated the higher binding strength of 2b and thus, established the potential of compound 2b as a lead molecule as an HIV-RT inhibitor.

Cell Surface Fibroblast Activation Protein-2 (Fap2) of Fusobacterium nucleatum as a Vaccine Candidate for Therapeutic Intervention of Human Colorectal Cancer: An Immunoinformatics Approach.

Colorectal cancer (CRC) is one of the most common cancers and is the second-highest in cancer-related deaths worldwide. The changes in gut homeostasis and microbial dysbiosis lead to the initiation of the tumorigenesis process. Several pathogenic gram-negative bacteria including Fusobacterium nucleatum are the principal contributors to the induction and pathogenesis of CRC. Thus, inhibiting the growth and survival of these pathogens can be a useful intervention strategy. Fibroblast activation protein-2 (Fap2) is an essential membrane protein of F. nucleatum that promotes the adherence of the bacterium to the colon cells, recruitment of immune cells, and induction of tumorigenesis. The present study depicts the design of an in silico vaccine candidate comprising the B-cell and T-cell epitopes of Fap2 for improving cell-mediated and humoral immune responses against CRC. Notably, this vaccine participates in significant protein-protein interactions with human Toll-like receptors, especially with TLR6 reveals, which is most likely to be correlated with its efficacy in eliciting potential immune responses. The immunogenic trait of the designed vaccine was verified by immune simulation approach. The cDNA of the vaccine construct was cloned in silico within the expression vector pET30ax for protein expression. Collectively, the proposed vaccine construct may serve as a promising therapeutic in intervening F. nucleatum-induced human CRC.

Identification and Investigation of a Cryptic Binding Pocket of the P37 Envelope Protein of Monkeypox Virus by Molecular Dynamics Simulations.

The spread of the monkeypox virus has surged during the unchecked COVID-19 epidemic. The most crucial target is the viral envelope protein, p37. However, lacking p37's crystal structure is a significant hurdle to rapid therapeutic discovery and mechanism elucidation. Structural modeling and molecular dynamics (MD) of the enzyme with inhibitors reveal a cryptic pocket occluded in the unbound structure. For the first time, the inhibitor's dynamic flip from the active to the cryptic site enlightens p37's allosteric site, which squeezes the active site, impairing its function. A large force is needed for inhibitor dissociation from the allosteric site, ushering in its biological importance. In addition, hot spot residues identified at both locations and discovered drugs more potent than tecovirimat may enable even more robust inhibitor designs against p37 and accelerate the development of monkeypox therapies.

Hydroxyalkynyl uracil derivatives as NNRTIs against HIV-1: in silico predictions, synthesis, docking and molecular dynamics simulation studies.

To improve rationally the efficacy of the non-nucleoside human immunodeficiency virus (HIV-1) inhibitors, it is important to have a precise and detailed understanding of the HIV-1 reverse transcriptase (RT) and inhibitor interactions. For the 1-[(2-hydroxyethoxy) methyl]-6-(phenylthio) thymine (HEPT) type of nucleoside reverse transcriptase inhibitors (NNRTIs), the H-bond between the N-3H of the inhibitor and the backbone carbonyl group of K101 represents the major hydrophilic interaction. This H-bond contributes to the NNRTI binding affinity. The descriptor analyses of different uracil derivatives proved their good cell internalization. The bioactivity score reflected higher drug likeness score and the ligands showed interesting docking results. All molecules were deeply buried and stabilized into the allosteric site of HIV-1 RT. For majority of molecules, residues Lys101, Lys103, Tyr181 and Tyr188 were identified as key protein residues responsible for generation of H-bond and major interactions were similar to all known NNRTIs while very few molecules interacted with residues Phe227 and Tyr318. The TOPKAT protocol available in Discovery Studio 3.0 was used to predict the pharmacokinetics of the designed uracil derivatives in the human body. The molecular dynamics (MD) and post-MD analyses results reflected that the complex HIVRT:5 appeared to be more stable than the complex HIVRT:HEPT, where HEPT was used as reference. Different uracil derivatives have been synthesized by using uracil as starting material and commercially available propargyl bromide. The N-1 derivative of uracil was further reacted with sodamide and different aldehydes/ketones bearing alkyl and phenyl ring to obtain hydroxyalkynyl uracil derivatives as NNRTIs.Communicated by Ramaswamy H. Sarma.

Designing AbhiSCoVac - A single potential vaccine for all 'corona culprits': Immunoinformatics and immune simulation approaches.

The coronaviridae family has generated highly virulent viruses, including the ones responsible for three major pandemics in last two decades with SARS in 2002, MERS outbreak in 2012 and the current nCOVID19 crisis that has turned the world breadthless. Future outbreaks are also a plausible threat to mankind. As computational biologists, we are committed to address the need for a universal vaccine that can deter all these pathogenic viruses in a single shot. Notably, the spike proteins present in all these viruses function as credible PAMPs that are majorly sensed by human TLR4 receptors. Our study aims to recognize the amino acid sequence(s) of the viral spike proteins that are precisely responsible for interaction with human TLR4 and to screen the immunogenic epitopes present in them to develop a multi-epitope multi-target chimeric vaccine against the coronaviruses. Molecular design of the constructed vaccine peptide is qualified in silico; additionally, molecular docking and molecular dynamics simulation studies collectively reveal strong and stable interactions of the vaccine construct with TLRs and MHC receptors. In silico cloning is performed for proficient expression in bacterial systems. In silico immune simulation of the vaccine indicates highly immunogenic nature of the vaccine construct without any allergic response. The present biocomputational study hereby innovates a vaccine candidate - AbhiSCoVac hypothesized as a potent remedy to combat all the virulent forms of coronaviruses.

Repurposing of FDA-approved drugs as potential inhibitors of the SARS-CoV-2 main protease: Molecular insights into improved therapeutic discovery.

With numerous infections and fatalities, COVID-19 has wreaked havoc around the globe. The main protease (Mpro), which cleaves the polyprotein to form non-structural proteins, thereby helping in the replication of SARS-CoV-2, appears as an attractive target for antiviral therapeutics. As FDA-approved drugs have shown effectiveness in targeting Mpro in previous SARS-CoV(s), molecular docking and virtual screening of existing antiviral, antimalarial, and protease inhibitor drugs were carried out against SARS-CoV-2 Mpro. Among 53 shortlisted drugs with binding energies lower than that of the crystal-bound inhibitor α-ketoamide 13 b (-6.7 kcal/mol), velpatasvir, glecaprevir, grazoprevir, baloxavir marboxil, danoprevir, nelfinavir, and indinavir (-9.1 to -7.5 kcal/mol) were the most significant on the list (hereafter referred to as the 53-list). Molecular dynamics (MD) simulations confirmed the stability of their Mpro complexes, with the MMPBSA binding free energy (ΔGbind) ranging between -124 kJ/mol (glecaprevir) and -28.2 kJ/mol (velpatasvir). Despite having the lowest initial binding energy, velpatasvir exhibited the highest ΔGbind value for escaping the catalytic site during the MD simulations, indicating its reduced efficacy, as observed experimentally. Available inhibition assay data adequately substantiated the computational forecast. Glecaprevir and nelfinavir (ΔGbind = -95.4 kJ/mol) appear to be the most effective antiviral drugs against Mpro. Furthermore, the remaining FDA drugs on the 53-list can be worth considering, since some have already demonstrated antiviral activity against SARS-CoV-2. Hence, theoretical pKi (Ki = inhibitor constant) values for all 53 drugs were provided. Notably, ΔGbind directly correlates with the average distance of the drugs from the His41-Cys145 catalytic dyad of Mpro, providing a roadmap for rapid screening and improving the inhibitor design against SARS-CoV-2 Mpro.

Binding mechanism and structural insights into the identified protein target of COVID-19 and importin-α with in-vitro effective drug ivermectin.

While an FDA approved drug Ivermectin was reported to dramatically reduce the cell line of SARS-CoV-2 by ∼5000 folds within 48 h, the precise mechanism of action and the COVID-19 molecular target involved in interaction with this in-vitro effective drug are unknown yet. Among 12 different COVID-19 targets along with Importin-α studied here, the RNA dependent RNA polymerase (RdRp) with RNA and Helicase NCB site show the strongest affinity to Ivermectin amounting -10.4 kcal/mol and -9.6 kcal/mol, respectively, followed by Importin-α with -9.0 kcal/mol. Molecular dynamics of corresponding protein-drug complexes reveals that the drug bound state of RdRp with RNA has better structural stability than the Helicase NCB site and Importin-α, with MM/PBSA free energy of -187.3 kJ/mol, almost twice that of Helicase (-94.6 kJ/mol) and even lower than that of Importin-α (-156.7 kJ/mol). The selectivity of Ivermectin to RdRp is triggered by a cooperative interaction of RNA-RdRp by ternary complex formation. Identification of the target and its interaction profile with Ivermectin can lead to more powerful drug designs for COVID-19 and experimental exploration.

In-silico evidences on filarial cystatin as a putative ligand of human TLR4.

Cystatin is a small molecular weight immunomodulatory protein of filarial parasite that plays a pivotal role in downregulating the host immune response to prolong the survival of the parasite inside the host body. Hitherto, this protein is familiar as an inhibitor of human proteases. However, growing evidences on the role of cystatin in regulating inflammatory homeostasis prompted us to investigate the molecular reasons behind the explicit anti-inflammatory trait of this protein. We have explored molecular docking and molecular dynamics simulation approaches to explore the interaction of cystatin of Wuchereria bancrofti (causative parasite of human filariasis) with human Toll-like receptors (TLRs). TLRs are the most crucial component of frontline host defence against pathogenic infections including filarial infection. Our in-silico data clearly revealed that cystatin strongly interacts with the extracellular domain of TLR4 (binding energy=-93.5 ± 10 kJ/mol) and this biophysical interaction is mediated by hydrogen bonding and hydrophobic interaction. Molecular dynamics simulation analysis revealed excellent stability of the cystatin-TLR4 complex. Taken together, our data indicated that cystatin appears to be a ligand of TLR4 and we hypothesize that cystatin-TLR4 interaction most likely to play a key role in activating the alternative activation pathways to establish an anti-inflammatory milieu. Thus, the study provokes the development of chemotherapeutics and/or vaccines for targeting the cystatin-TLR4 interaction to disrupt the pathological attributes of human lymphatic filariasis. Our findings are expected to provide a novel dimension to the existing knowledge on filarial immunopathogenesis and it will encourage the scientific communities for experimental validation of the present investigation. Communicated by Ramaswamy H. Sarma.

ACE-2-Derived Biomimetic Peptides for the Inhibition of Spike Protein of SARS-CoV-2.

SARS-CoV-2, a novel coronavirus causing overwhelming death and infection worldwide, has emerged as a pandemic. Compared to its predecessor SARS-CoV, SARS-CoV-2 is more infective for being highly contagious and exhibiting tighter binding with host angiotensin-converting enzyme 2 (hACE-2). The entry of the virus into host cells is mediated by the interaction of its spike protein with hACE-2. Thus, a peptide that has a resemblance to hACE-2 but can overpower the spike protein-hACE-2 interaction will be a potential therapeutic to contain this virus. The non-interacting residues in the receptor-binding domain of hACE-2 have been mutated to generate a library of 136 new peptides. Out of this library, docking and virtual screening discover seven peptides that can exert a stronger interaction with the spike protein than hACE-2. A peptide derived from simultaneous mutation of all the non-interacting residues of hACE-2 yields almost three-fold stronger interaction than hACE-2 and thus turns out here to be the best peptide inhibitor of the novel coronavirus. The binding of the best peptide inhibitor with the spike protein is explored further by molecular dynamics, free energy, and principal component analysis, which demonstrate its efficacy compared to hACE-2. The delivery of the screened inhibitors with nanocarriers like metal-organic frameworks will be worthy of further consideration to boost their efficacy.

Screening and molecular characterization of lethal mutations of human homogentisate 1, 2 dioxigenase.

Alkaptonuria (AKU) is an autosomal recessive disorder, which is caused by a site-specific mutation(s) and thus, impaired the function of Homogentisate-1, 2-dioxygenase (HGD), an essential enzyme for the catabolism of phenylalanine and tyrosine. Among frameshift, intronic, splice-site and missense mutations, the latter has been the most common form of genetic variations for the disease. How do the acquired mutations in HGD correlate with the disease? Systematic staged-screening of some sixty-five mutations, which are known to have a relation with the disease, by GVGD, SIFT, SNAP, PANTHER, SDM, PHD-SNP, Meta-SNP, Pmut and Mutpred methods, showed that mutations, W60G, A122D and V300G are potentially related with the severity of AKU. Detailed analyses on molecular docking and molecular dynamics simulation (MDS) of these mutants against the wild-type HGD reveal the loss of structural and molecular dynamic properties of the enzyme. Further, the observed conformational flexibility in mutants at targeted peptide segments seems to have a relation with the impairment of the function of HGD. Taken together, the study involves a designed computational methodology to analyse the disease-associated nsSNPs for AKU, the knowledge of which seems to have potential applications in drug therapies for the disease in particular and other similar systems in general.Communicated by Ramaswamy H. Sarma.

Binding insight of clinically oriented drug famotidine with the identified potential target of SARS-CoV-2.

Communicated by Ramaswamy H. Sarma.

Ivermectin, Famotidine, and Doxycycline: A Suggested Combinatorial Therapeutic for the Treatment of COVID-19.

At times, combination therapy has proven to be very effective. While no cure is available to date, herein we put forward with rationale and supporting evidence that if administrated simultaneously, a combination of FDA-approved drugs comprising ivermectin, famotidine, and doxycycline may provide robust chemoprophylaxis effective against COVID-19.

Alkylated benzimidazoles: Design, synthesis, docking, DFT analysis, ADMET property, molecular dynamics and activity against HIV and YFV.

A series of alkylated benzimidazole derivatives was synthesized and screened for their anti-HIV, anti-YFV, and broad-spectrum antiviral properties. The physicochemical parameters and drug-like properties of the compounds were assessed first, and then docking studies and MD simulations on HIV-RT allosteric sites were conducted to find the possible mode of their action. DFT analysis was also performed to confirm the nature of the hydrogen bonding interaction of active compounds. The in silico studies indicated that the molecules behaved like possible NNRTIs. The nature - polar or non-polar and position of the substituent present at fifth, sixth, and N-1 positions of the benzimidazole moiety played an important role in determining the antiviral properties of the compounds. Among the various compounds, 2-(5,6-dibromo-2-chloro-1H-benzimidazol-1-yl)ethan-1-ol (3a) showed anti-HIV activity with an appreciably low IC50 value as 0.386 × 10-5µM. Similarly, compound 2b, 3-(2-chloro-5-nitro-1H-benzimidazol-1-yl) propan-1-ol, showed excellent inhibitory property against the yellow fever virus (YFV) with EC50 value as 0.7824 × 10-2µM.

AUTOMINv1.0: an automation for minimization of Protein Data Bank files and its usage.

Global minimal structure of protein/enzyme is energetically compromised that maintains an intricate balance between the rigidity and the flexibility. Such a state makes it interactive to its ligand molecules. Although protein data bank files (PDB) may have achieved the state, in many situations minimization has been crucial to overcome unwanted steric clashes, and other conformational strains. It is more so, when orthologous PDB structures that are intended in a given study, show variations in resolution, R-factor, shell-water contents, loop characteristics etc. Here, a fully automated procedure of minimization would be highly useful. AUTOMINv1.0 is such an automation of minimization that runs on any number of structure files with any number of chains in them along with the inclusion of selective/non-selective shell-waters interacting with polar and or non-polar atom-types of protein. Comparison of the mean binaryitems of salt-bridges of minimized and un-minimized structures (chains > 100) of nucleoside diphosphate kinase from mimi virus shows dramatic improvements in the earlier. Again, the mean steric clashes of 2AZ3.pdb are reduced upon minimization. Remarkably, the observed steric clashes between shell-waters and atom-types of protein are seen to be removed upon minimization. Taken together, AUTOMINv1.0 is an automation of minimization that finds applications in structural bioinformatics.

Insight into SNPs and epitopes of E protein of newly emerged genotype-I isolates of JEV from Midnapur, West Bengal, India.

BACKGROUND: Japanese encephalitis virus (JEV) is a mosquito-borne flavivirus that causes Japanese Encephalitis (JE) and Acute Encephalitis Syndrome (AES) in humans. Genotype-I (as co-circulating cases with Genotype-III) was isolated in 2010 (JEV28, JEV21) and then in 2011 (JEV45) from Midnapur district, West Bengal (WB) for the first time from clinical patients who were previously been vaccinated with live attenuated SA14-14-2 strain. We apply bioinformatics and immunoinformatics on sequence and structure of E protein for analysis of crucial substitutions that might cause the genotypic transition, affecting protein-function and altering specificity of epitopes. RESULTS: Although frequency of substitutions in E glycoprotein of JEV28, JEV21 and JEV45 isolates vary, its homologous patterns remain exactly similar as earlier Japan isolate (Ishikawa). Sequence and 3D model-structure based analyses of E protein show that only four of all substitutions are critical for genotype-I specific effect of which N103K is common among all isolates indicating its role in the transition of genotype-III to genotype-I. Predicted B-cell and T-cell epitopes are seen to harbor these critical substitutions that affect overall conformational stability of the protein. These epitopes were subjected to conservation analyses using a large set of the protein from Asian continent. CONCLUSIONS: The study identifies crucial substitutions that contribute to the emergence of genotype-I. Predicted epitopes harboring these substitutions may alter specificity which might be the reason of reported failure of vaccine. Conservation analysis of these epitopes would be useful for design of genotype-I specific vaccine.

Salt-bridge energetics in halophilic proteins.

Halophilic proteins have greater abundance of acidic over basic and very low bulky hydrophobic residues. Classical electrostatic stabilization was suggested as the key determinant for halophilic adaptation of protein. However, contribution of specific electrostatic interactions (i.e. salt-bridges) to overall stability of halophilic proteins is yet to be understood. To understand this, we use Adaptive-Poison-Boltzmann-Solver Methods along with our home-built automation to workout net as well as associated component energy terms such as desolvation energy, bridge energy and background energy for 275 salt-bridges from 20 extremely halophilic proteins. We then perform extensive statistical analysis on general and energetic attributes on these salt-bridges. On average, 8 salt-bridges per 150 residues protein were observed which is almost twice than earlier report. Overall contributions of salt-bridges are -3.0 kcal mol-1. Majority (78%) of salt-bridges in our dataset are stable and conserved in nature. Although, average contributions of component energy terms are equal, their individual details vary greatly from one another indicating their sensitivity to local micro-environment. Notably, 35% of salt-bridges in our database are buried and stable. Greater desolvation penalty of these buried salt-bridges are counteracted by stable network salt-bridges apart from favorable equal contributions of bridge and background terms. Recruitment of extensive network salt-bridges (46%) with a net contribution of -5.0 kcal mol-1 per salt-bridge, seems to be a halophilic design wherein favorable average contribution of background term (-10 kcal mol-1) exceeds than that of bridge term (-7 kcal mol-1). Interiors of proteins from halophiles are seen to possess relatively higher abundance of charge and polar side chains than that of mesophiles which seems to be satisfied by cooperative network salt-bridges. Overall, our theoretical analyses provide insight into halophilic signature in its specific electrostatic interactions which we hope would help in protein engineering and bioinformatics studies.