The intracellular environment is overcrowded with a range of molecules (small and large), all of which influence protein conformation. As a result, understanding how proteins fold and stay functional in such crowded conditions is essential. Several in vitro experiments have looked into the effects of macromolecular crowding on different proteins. However, there are hardly any reports regarding small molecular crowders used alone and in mixtures to observe their effects on the structure and stability of the proteins, which mimics of the cellular conditions. Here we investigate the effect of different mixtures of crowders, ethylene glycol (EG) and its polymer polyethylene glycol (PEG 400 Da) on the structural and thermal stability of myoglobin (Mb). Our results show that monomer (EG) has no significant effect on the structure of Mb, while the polymer disrupts its structure and decreases its stability. Conversely, the additive effect of crowders showed structural refolding of the protein to some extent. Moreover, the calorimetric binding studies of the protein showed very weak interactions with the mixture of crowders. Usually, we can assume that soft interactions induce structural perturbations while exclusion volume effects stabilize the protein structure; therefore, we hypothesize that under in vivo crowded conditions, both phenomena occur and maintain the stability and function of proteins.

Studying drug-protein interactions has gained significant attention lately, and this is because the majority of drugs interact with proteins, thereby altering their structure and, moreover, their functionality. Rivastigmine tartrate (RT) is a drug that is in use for mild to moderate Alzheimer therapy. This study was targeted to characterize the interaction between human transferrin (hTf) and RT by employing spectroscopy, isothermal titration calorimetry (ITC), and molecular docking studies. Experimental results of fluorescence quenching of hTf induced by RT implied the formation of a static complex between hTf and RT. Further elucidation of the observed fluorescence data retorting Stern-Volmer and modified Stern-Volmer resulted in binding constants for hTf-RT complex of the order 104 M-1 over the studied temperatures. Thermodynamic parameters of hTf-RT interaction were elucidated further by employing these obtained binding constant values. It was quite evident from obtained thermodynamic attributes that RT spontaneously binds to hTf with a postulated existence of hydrogen bonding or Van der Waals forces. Further, Circular dichroism spectroscopy (CD) also confirmed RT-hTf complex formation owing to upward movement of CD spectra in the presence of RT. ITC profiles advocated the existence of reaction to be spontaneous. Moreover, molecular docking further revealed that the important residues play a pivotal role in RT-hTf interaction. The findings of this study can be of a significant benefit to the drug-designing industry in this disease-prone era.

One of the significant proteins that have attracted research groups due to virtue of being a potent selective anticancer drug target and property of triggering apoptosis upon release in cytoplasm is cytochrome c (cyt c). The mechanical transformations due to the macromolecular crowding in membrane in the mammalian cell are proposed to be useful inductors of changes in volume. It is very interesting to know that mitochondrial function were observed to be improved by polyethylene glycol (PEG) interaction, which in turn inhibits the cyt c (a pro-apoptotic cell death factor). In this work, the effect of polyethylene glycol of molecular weight 4 kilo Dalton (PEG 4 kDa) was investigated to highlight the structural transformations (tertiary and secondary structure) in cyt c using a choice of spectroscopic techniques (including UV-Vis absorption, near-UV, far-UV and Soret circular dichroism and fluorescence spectroscopy), which shows noteworthy shifts in the secondary and tertiary structures at higher concentrations of PEG 4 kDa with small changes in the heme-globular interactions. The size distribution changes of native protein treated with various concentrations of the crowder were observed and analyzed by dynamic light scattering (DLS). The interaction studies of the crowder with the protein was observed and analyzed by FTIR, isothermal titration calorimetry, time resolved fluorescence and molecular docking. The investigations suggested that the structural changes in the protein occurred due to soft interactions of PEG 4 kDa, which usually destabilizes proteins. The experimental evidence in this study proposed that crowding could be another approach to mechanical super-competition and free of certain markers that could aid in the identification and control of various diseases. This study suggests that crowders at specific concentrations, which softly interact with proteins, can be exploited as remedy for various diseases.

Due to the lack of efficient therapeutic options and clinical trial limitations, the FDA-approved drugs can be a good choice to handle Coronavirus disease (COVID-19). Many reports have enough evidence for the use of FDA-approved drugs which have inhibitory potential against target proteins of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Here, we utilized a structure-based drug design approach to find possible drug candidates from the existing pool of FDA-approved drugs and checked their effectiveness against the SARS-CoV-2. We performed virtual screening of the FDA-approved drugs against the main protease (Mpro) of SARS-CoV-2, an essential enzyme, and a potential drug target. Using well-defined computational methods, we identified Glecaprevir and Maraviroc (MVC) as the best inhibitors of SARS-CoV-2 Mpro. Both drugs bind to the substrate-binding pocket of SARS-CoV-2 Mpro and form a significant number of non-covalent interactions. Glecaprevir and MVC bind to the conserved residues of substrate-binding pocket of SARS-CoV-2 Mpro. This work provides sufficient evidence for the use of Glecaprevir and MVC for the therapeutic management of COVID-19 after experimental validation and clinical manifestations.

Integrin-linked kinase (ILK) is a member of Ser/Thr kinase which interacts to the cytoplasmic domain of β-integrins, and thereby induces apoptosis. ILK is considered as potential drug target because it's direct involvement in the tumor progression. Here, we have performed molecular docking followed by 100 ns MD simulation to understand the mechanism of interaction of ILK with the ellagic acid (EA). EA is well known for its antiproliferative and antioxidant properties in cancer cell lines and animal models. We have observed that EA binds to the active site cavity of ILK and causes conformational changes in the ILK structure. The orientation of EA in the active pocket of ILK showed to have least RMSD values and stable. The average binding energy ILK-EA complex calculated during MMPBSA was -191.267 kJ/mol, indicating a relatively strong binding affinity. The actual binding affinity of EA to ILK was measured by fluorescence spectroscopy and Kb and n values were 9.28 μM and 1.9264 (~2), respectively. The IC50 values for EA were 26.22 ± 0.12 μM for MCF-7 and 38.45 ± 2.42 μM for HepG2 cells, estimated by MTT assay. Our findings are helpful to design EA-based novel inhibitors of ILK which have potential to attenuate tumor progression.

Investigation of changes in thermal stabilities and structures of proteins in the presence of different co-solutes (ligands) is an integral part in the basic research, discovery, and development of drugs. Ethylene glycol (EG) is known to be toxic and causes teratogenic, inducing primarily skeletal and external malformations and other diseases. The effect of EG on the structure and thermal stability of myoglobin (Mb) was studied using various spectroscopic techniques at pH 7.0 and two different temperatures. As revealed by circular dichroism, Trp fluorescence, nano-DSF, and absorption (UV and visible) measurements, EG (i) has no significant effect on secondary and tertiary structures of Mb at 25 °C, and (ii) it decreases the thermal stability of the protein, which increases with increasing concentration of EG. As revealed by ANS (8-anilino-1-naphthalene sulfonic acid) fluorescence measurements, heat-induced denatured protein has newly exposed hydrophobic patches that bind to ANS. Isothermal titration calorimetry revealed that the interaction between EG and Mb is temperature dependent; the preferential interaction of EG is entropy driven at low temperature, 298 K (25 °C), and it is enthalpy driven at higher temperature, 343 K (70 °C). Molecular docking study showed that EG interacts with side chains of amino acid residues of Mb through van der Waals interactions and hydrogen bonding.

Sphingosine kinase 1 (SphK1) has recently gained attention as a potential drug target for its association with cancer and other inflammatory diseases. Here, we have investigated the binding affinity of dietary phytochemicals viz., ursolic acid, capsaicin, DL-α tocopherol acetate, quercetin, vanillin, citral, limonin and simvastatin with the SphK1. Docking studies revealed that all these compounds bind to the SphK1 with varying affinities. Fluorescence binding and isothermal titration calorimetric measurements suggested that quercetin and capsaicin bind to SphK1 with an excellent affinity, and significantly inhibits its activity with an admirable IC50 values. The binding mechanism of quercetin was assessed by docking and molecular dynamics simulation studies for 100 ns in detail. We found that quercetin acts as a lipid substrate competitive inhibitor, and it interacts with important residues of active-site pocket through hydrogen bonds and other non-covalent interactions. Quercetin forms a stable complex with SphK1 without inducing any significant conformational changes in the protein structure. In conclusion, we infer that quercetin and capsaicin provide a chemical scaffold to develop potent and selective inhibitors of SphK1 after required modifications for the clinical management of cancer.

Vitamin D deficiency significantly correlates with the severity of SARS-CoV-2 infection. Molecular docking-based virtual screening studies predict that novel vitamin D and related lumisterol hydroxymetabolites are able to bind to the active sites of two SARS-CoV-2 transcription machinery enzymes with high affinity. These enzymes are the main protease (Mpro) and RNA-dependent RNA polymerase (RdRP), which play important roles in viral replication and establishing infection. Based on predicted binding affinities and specific interactions, we identified 10 vitamin D3 (D3) and lumisterol (L3) analogs as likely binding partners of SARS-CoV-2 Mpro and RdRP and, therefore, tested their ability to inhibit these enzymes. Activity measurements demonstrated that 25(OH)L3, 24(OH)L3, and 20(OH)7DHC are the most effective of the hydroxymetabolites tested at inhibiting the activity of SARS-CoV-2 Mpro causing 10%-19% inhibition. These same derivatives as well as other hydroxylumisterols and hydroxyvitamin D3 metabolites inhibited RdRP by 50%-60%. Thus, inhibition of these enzymes by vitamin D and lumisterol metabolites may provide a novel approach to hindering the SARS-CoV-2 infection.NEW & NOTEWORTHY Active forms of vitamin D and lumisterol can inhibit SARS-CoV-2 replication machinery enzymes, which indicates that novel vitamin D and lumisterol metabolites are candidates for antiviral drug research.

Sphingosine kinase 1 (SphK1) is one of the well-studied drug targets for cancer and inflammatory diseases. Recently discovered small-molecule inhibitors of SphK1 have been recommended in cancer therapeutics; however, selectivity and potency of first-generation inhibitors are great challenge. In search of effective SphK1 inhibitors, a set of small molecules have been designed and synthesized bearing urea, sulfonylurea, sulfonamide, and sulfonyltriurea groups. The binding affinity of these inhibitors was measured by fluorescence-binding assay and isothermal titration calorimetry. Compounds 1, 5, 6, and 7 showed an admirable binding affinity to the SphK1 in the sub-micromolar range and significantly inhibited SphK1 activity with admirable IC50 values. Molecular docking studies revealed that these compounds fit well into the sphingosine binding pocket of SphK1 and formed significant number of hydrogen bonds and van der Waals interactions. These molecules may be exploited as potent and selective inhibitors of SphK1 that could be implicated in cancer therapeutics after the required in vivo validation.

To combat emerging antimicrobial-resistant microbes, there is an urgent need to develop new antimicrobials with better therapeutic profiles. For this, a series of 13 new spiropyrrolidine derivatives were designed, synthesized, characterized and evaluated for their in vitro antimicrobial, antioxidant and antidiabetic potential. Antimicrobial results revealed that the designed compounds displayed good activity against clinical isolated strains, with 5d being the most potent (MIC 3.95 mM against Staphylococcus aureus ATCC 25923) compared to tetracycline (MIC 576.01 mM). The antioxidant activity was assessed by trapping DPPH, ABTS and FRAP assays. The results suggest remarkable antioxidant potential of all synthesized compounds, particularly 5c, exhibiting the strongest activity with IC50 of 3.26 ± 0.32 mM (DPPH), 7.03 ± 0.07 mM (ABTS) and 3.69 ± 0.72 mM (FRAP). Tested for their α-amylase inhibitory effect, the examined analogues display a variable degree of α-amylase activity with IC50 ranging between 0.55 ± 0.38 mM and 2.19 ± 0.23 mM compared to acarbose (IC50 1.19 ± 0.02 mM), with the most active compounds being 5d, followed by 5c and 5j, affording IC50 of 0.55 ± 0.38 mM, 0.92 ± 0.10 mM, and 0.95 ± 0.14 mM, respectively. Preliminary structure-activity relationships revealed the importance of such substituents in enhancing the activity. Furthermore, the ADME screening test was applied to optimize the physicochemical properties and determine their drug-like characteristics. Binding interactions and stability between ligands and active residues of the investigated enzymes were confirmed through molecular docking and dynamic simulation study. These findings provided guidance for further developing leading new spiropyrrolidine scaffolds with improved dual antimicrobial and antidiabetic activities.

Calcium/calmodulin-dependent protein kinase IV (CaMKIV) is an upstream regulator of CaMKK-CaMKIV signaling cascade that activates various transcription factors, thereby regulating several cellular activities including, neuronal communication and immune response. Owing to the abnormal expression in cancer and neurodegenerative diseases, the CaMKIV has been considered a potential drug target. In the present study, we checked the binding affinity of plant-derived natural compounds viz., quercetin, ellagic acid (EA), simvastatin, capsaicin, ursolic acid, DL-α-tocopherol acetate, and limonin towards CaMKIV. Molecular docking and fluorescence binding studies showed that EA and quercetin bind to the CaMKIV with a considerable affinity in comparison to other compounds. Enzyme inhibition assay revealed that both EA and quercetin inhibit CaMKIV activity with their IC50 values in the micromolar range. To get atomistic insights into the mode of interactions, inhibition mechanism, and the stability of the CaMKIV-ligand complex, a 100 ns MD simulation analysis was performed. Both EA and quercetin bind to the catalytically important residues of active site pocket of CaMKIV forming enough stabilizing interactions presumably inhibiting enzyme activity. Moreover, no significant structural change in the CaMKIV was observed upon binding of EA and quercetin. In conclusion, this study illustrates the application of phytoconstituents in the development of therapeutic molecules targeting CaMKIV having implications in cancer and neurodegenerative diseases after in vivo validation.

Metabolic reprogramming is a key attribute of cancer progression. An altered expression of pyruvate kinase M2 (PKM2), a phosphotyrosine-binding protein is observed in many human cancers. PKM2 plays a vital role in metabolic reprogramming, transcription and cell cycle progression and thus is deliberated as an attractive target in anticancer drug development. The expression of PKM2 is essential for aerobic glycolysis and cell proliferation, especially in cancer cells, facilitating selective targeting of PKM2 in cell metabolism for cancer therapeutics. We have screened a virtual library of phytochemicals from the IMPPAT (Indian Medicinal Plants, Phytochemistry and Therapeutics) database of Indian medicinal plants to identify potential activators of PKM2. The initial screening was carried out for the physicochemical properties of the compounds, and then structure-based molecular docking was performed to select compounds based on their binding affinity towards PKM2. Subsequently, the ADMET (absorption, distribution, metabolism, excretion and toxicity) properties, PAINS (Pan-assay interference compounds) patterns, and PASS evaluation were carried out to find more potent hits against PKM2. Here, Tuberosin was identified from the screening process bearing appreciable binding affinity toward the PKM2-binding pocket and showed a worthy set of drug-like properties. Finally, molecular dynamics simulation for 100 ns was performed, which showed decent stability of the protein-ligand complex and relatival conformational dynamics throughout the trajectory. The study suggests that modulating PKM2 with natural compounds is an attractive approach in treating human malignancy after required validation.

Thiazolidinedione derivatives (TZDs) have attracted attention because of their pharmacological effects. For example, certain TZDs have been reported to ameliorate type II diabetes by binding and activating PPARs (peroxisome proliferator-activated receptors). Nonetheless, no information is available on the interaction between the heterocyclic 2, 4-thiazolidinedione (2,4-TZD) moiety and serum albumin, which could affect the pharmacokinetics and pharmacodynamics of TZDs. In this study, we investigated the binding of 2,4-TZD to human serum albumin (HSA). Intrinsic fluorescence spectroscopy revealed a 1:1 binding stoichiometry between 2,4-TZD and HSA with a binding constant (Kb) of 1.69 ± 0.15 × 103 M-1 at 298 K. Isothermal titration calorimetry studies showed that 2,4-TZD/HSA binding was an exothermic and spontaneous reaction. Molecular docking analysis revealed that 2,4-TZD binds to HSA subdomain IB and that the complex formed is stabilized by van der Waal's interactions and hydrogen bonds. Molecular dynamics simulation confirmed the stability of the HSA-TZD complex. Further, circular dichroism and 3D fluorescence studies showed that the global conformation of HSA was slightly altered by 2,4-TZD binding, enhancing its stability. The results obtained herein further help in understanding the pharmacokinetic properties of thiazolidinedione.

Microtubule affinity-regulating kinase (MARK4) plays a key role in Alzheimer's disease (AD) development as its overexpression is directly linked to increased tau phosphorylation. MARK4 is a potential drug target of AD and is thus its structural features are employed in the development of new therapeutic molecules. Donepezil (DP) and rivastigmine tartrate (RT) are acetylcholinesterase (AChE) inhibitors and are used to treat symptomatic patients of mild to moderate AD. In keeping with the therapeutic implications of DP and RT in AD, we performed binding studies of these drugs with the MARK4. Both DP and RT bound to MARK4 with a binding constant (K) of 107 M-1. The temperature dependency of binding parameters revealed MARK-DP complex to be guided by static mode while MARK-RT complex to be guided by both static and dynamic quenching. Both drugs inhibited MARK4 with IC50 values of 5.3 μM (DP) and 6.74 μM (RT). The evaluation of associated enthalpy change (ΔH) and entropy change (ΔS) implied the complex formation to be driven by hydrogen bonding making it seemingly strong and specific. Isothermal titration calorimetry further advocated a spontaneous binding. In vitro observations were further complemented by the calculation of binding free energy by molecular docking and interactions with the functionally-important residues of the active site pocket of MARK4. This study signifies the implications of AChE inhibitors, RT, and DP in Alzheimer's therapy targeting MARK4.

Protein stability is a requisite in the field of biotechnology, cell biology and drug design. To understand effects of amino acid substitutions, computational models are preferred to save time and expenses. As a systemically important, highly abundant, stable protein, the knowledge of Cu/Zn Superoxide dismutase1 (SOD1) is important, making it a suitable test case for genotype-phenotype correlation in understanding ALS. Here, we report performance of eight protein stability calculators (PoPMuSiC 3.1, I-Mutant 2.0, I-Mutant 3.0, CUPSAT, FoldX, mCSM, BeatMusic and ENCoM) against 54 experimental stability changes due to mutations of SOD1. Four different high-resolution structures were used to test structure sensitivity that may affect protein calculations. Bland-Altman plot was also used to assess agreement between stability analyses. Overall, PoPMuSiC and FoldX emerge as the best methods in this benchmark. The relative performance of all the eight methods was very much structure independent, and also displayed less structural sensitivity. We also analyzed patient's data in relation to experimental and computed protein stabilities for mutations of human SOD1. Correlation between disease phenotypes and stability changes suggest that the changes in SOD1 stability correlate with ALS patient survival times. Thus, the results clearly demonstrate the importance of protein stability in SOD1 pathogenicity.

Sphingosine kinase 1 (SphK1) is an oncogenic lipid kinase that catalyzes the formation of sphingosine-1-phosphate via phosphorylation of sphingosine and known to play a crucial role in angiogenesis, lymphocyte trafficking, signal transduction pathways, and response to apoptotic stimuli. SphK1 has received attention because of its involvement in varying types of cancer and inflammatory diseases such as rheumatoid arthritis, diabetes, renal fibrosis, pulmonary fibrosis, asthma, and neurodegenerative disorders. In the malignancies of breast, lung, uterus, ovary, kidney, and leukemia, overexpression of SphK1 has been reported and thus considered as a potential drug target. In this study, we have performed virtual high-throughput screening of ∼90,000 natural products from the ZINC database to find potential SphK1-inhibitors. Initially, the hits were selected by applying absorption, distribution, metabolism, excretion, and toxicity properties, Lipinski's rule, and PAINS filters. Further, docking analysis was performed to estimate the binding affinities and specificity to find safe and effective preclinical leads against SphK1. Two compounds, ZINC05434006 and ZINC04260971, bearing appreciable binding affinity and SphK1 selectivity were selected for 100 ns molecular dynamics (MD) simulations under explicit water conditions. The all-atom MD simulation results suggested that the ZINC05434006 and ZINC04260971 binding induces a slight structural change and stabilizes the SphK1 structure. In conclusion, we propose natural compounds, ZINC05434006 and ZINC04260971, as potential inhibitors of SphK1, which may be further exploited as potential leads to develop effective therapeutics against SphK1-associated diseases including cancer after in vitro and in vivo validations.

Microtubule affinity regulating kinase (MARK4) is a potential drug target for different types of cancer as it controls the early step of cell division. In this study, we have screened a series of natural compounds and finally identified rosmarinic acid (RA) as a potential inhibitor of MARK4. Molecular docking and 500 ns all-atom simulation studies suggested that RA binds to the active site pocket of MARK4, forming enough number of non-covalent interactions with critical residues and MARK4-RA complex is stable throughout the simulation trajectory. RA shows an excellent binding affinity to the MARK4 with a binding constant (K) of 107 M-1. Furthermore, RA significantly inhibits MARK4 activity (IC50 = 6.204 µM). The evaluation of enthalpy change (∆H) and entropy change (∆S) suggested that the MARK4-RA complex formation is driven by hydrogen bonding and thus complexation process is seemingly specific. The consequence of MARK4 inhibition by RA was further evaluated by cell-based tau-phosphorylation studies, which suggested that RA inhibited the phosphorylation of tau. The treatment of cancer cells with RA significantly controls cell growth and subsequently induces apoptosis. Our study provides a rationale for the therapeutic evaluation of RA and RA-based inhibitors in MARK4 associated cancers and other diseases.

Cyclin-Dependent Kinase 6 (CDK6) plays an important role in cancer progression, and thus, it is considered as an attractive drug target in anticancer therapeutics. This study presents an evaluation of dietary phytochemicals, capsaicin, tocopherol, rosmarinic acid, ursolic acid, ellagic acid (EA), limonene, caffeic acid, and ferulic acid for their potential to inhibit the activity of CDK6. Molecular docking and fluorescence binding studies revealed appreciable binding affinities of these compounds to the CDK6. Among them, EA shows the highest binding affinity for CDK6, and thus a molecular dynamics simulation study of 200 ns was performed to get deeper insights into the binding mechanism and stability of the CDK6-EA complex. Fluorescence binding studies revealed that EA binds to the CDK6 with a binding constant of K = 107 M-1 and subsequently inhibits its enzyme activity with an IC50 value of 3.053 µM. Analysis of thermodynamic parameters of CDK6-EA complex formation suggested a hydrophobic interaction driven process. The treatment of EA decreases the colonization of cancer cells and induces apoptosis. Moreover, the expression of CDK6 has been downregulated in EA-treated human breast cancer cell lines. In conclusion, this study establishes EA as a potent CDK6 inhibitor that can be further evaluated in CDK6 directed anticancer therapies.

Serum and glucocorticoid-regulated kinase 1 (SGK1) is a serine/threonine kinase that works under acute transcriptional control by several stimuli, including serum and glucocorticoids. It plays a significant role in the cancer progression and metastasis, as it regulates inflammation, apoptosis, hormone release, neuro-excitability, and cell proliferation. SGK1 has recently been considered as a potential drug target for cancer, diabetes, and neurodegenerative diseases. In the present study, we have performed structure-based virtual high-throughput screening of natural compounds from the ZINC database to find potential inhibitors of SGK1. Initially, hits were selected based on their physicochemical, absorption, distribution, metabolism, excretion, and toxicity (ADMET), and other drug-like properties. Afterwards, PAINS filter, binding affinities estimation, and interaction analysis were performed to find safe and effective hits. We found four compounds bearing appreciable binding affinity and specificity towards the binding pocket of SGK1. The docking results were complemented by all-atom molecular dynamics simulation for 100 ns, followed by MM/PBSA, and principal component analysis to investigate the conformational changes, stability, and interaction mechanism of SGK1 in-complex with the selected compound ZINC00319000. Molecular dynamics simulation results suggested that the binding of ZINC00319000 stabilizes the SGK1 structure, and it leads to fewer conformational changes. In conclusion, the identified compound ZINC00319000 might be further exploited as a scaffold to develop promising inhibitors of SGK1 for the therapeutic management of associated diseases, including cancer.

Helicobacter pylori is a Gram-negative bacteria that is responsible for gastritis in human. Its spiral flagellated body helps in locomotion and colonization in the host environment. It is capable of living in the highly acidic environment of the stomach with the help of acid adaptive genes. The genome of H. pylori 26695 strain contains 1,555 coding genes that encode 1,445 proteins. Out of these, 340 proteins are characterized as hypothetical proteins (HP). This study involves extensive analysis of the HPs using an established pipeline which comprises various bioinformatics tools and databases to find out probable functions of the HPs and identification of virulence factors. After extensive analysis of all the 340 HPs, we found that 104 HPs are showing characteristic similarities with the proteins with known functions. Thus, on the basis of such similarities, we assigned probable functions to 104 HPs with high confidence and precision. All the predicted HPs contain representative members of diverse functional classes of proteins such as enzymes, transporters, binding proteins, regulatory proteins, proteins involved in cellular processes and other proteins with miscellaneous functions. Therefore, we classified 104 HPs into aforementioned functional groups. During the virulence factors analysis of the HPs, we found 11 HPs are showing significant virulence. The identification of virulence proteins with the help their predicted functions may pave the way for drug target estimation and development of effective drug to counter the activity of that protein.