Functional Annotation and Molecular modeling of Hypothetical Proteins (HPs) from P.aeruginosa plasmid pUM505: An In silico Approach

Functional Annotation and Molecular modeling of Hypothetical Proteins (HPs) from P.aeruginosa plasmid pUM505: An In silico Approach

Srikant Awasthi1-2, Pragya Saxena2, Hillol Chakdar2, Alok K Srivastava2 and Salman Akhtar1*

1Department of Bioengineering, Integral University, Lucknow, INDIA, 226026

2Microbial Genomics Laboratory, National Bureau of Agriculturally Important Microorganisms (NBAIM), Mau, Uttar Pradesh, INDIA, 275101

Dr. Salman Akhtar Associate Professor

Department of Bioengineering Integral University, Lucknow, India 226026

Abstract:- Structural and function annotation of P. aeruginosa plasmid pUM505isessentially required to facilitate the understanding of mechanisms of pathogenesis and biochemical pathways important for selecting novel therapeutic target. In present study, randomly selected twelve hypothetical protein sequence of P. aeruginosa plasmid pUM505 has been annotated using various In-Silico tools and databases to determine domain family, solubility of protein, ligand binding sites etc. Six out of 12 proteins have been putatively annotated, in which four have been annotated with high confidence. Physio-chemical characterization revealed two proteins are stable. The three-dimensional structure of two important annotated proteins were modeled and their ligand binding sites were identified. Domains and families for six proteins have been found. The analysis revealed that these proteins have antitoxin activity, integrase enzyme activity, conjugal DNA transfer activity, etc. Structural prediction of these proteins and detection of binding sites from this study would indicate a potential target aiding docking studies for therapeutic designing against various diseases.

Keywords: Hypothetical Protein, Functional Annotation, Molecular Modeling, Docking

1. INTRODUCTION

   Pseudomonas aeruginosa, a gram-negative bacteriumis well known for its environmental versatility. Diverse growing habitat includes soil,coastal marine, plant and animal tissues (Khan et al., 2007). P. aeruginosais also well known for its multidrug resistant and is a global threat towards many infections disease. P. aeruginosais a major opportunistic pathogen in humans, causing serious complications caused by infections in patients particularly susceptible like people with immune system deficiencies, victims of skin burns, catheterized patients who suffer urinary tract infections and patients with respirators, causing nosocomial pneumonia (Lyczak et al., 2000). It is the major cause of mortality in patients with cystic fibrosis colonizing the lungs (Williams et al., 2010). Role of plasmids in antibiotic resistance are well establish.Plasmids are circular deoxyribonucleic acid molecules that exist in bacteria, usually independent of the chromosome. The study of plasmids is important to medical microbiology because plasmids can encode genes for antibiotic resistance or virulence factors (Wang et al., 1988). The pUM505 plasmid contains a genomic island with sequence similar to islands found in chromosomes of virulent P. aeruginosa clinical isolates. Plasmid pUM505 contains several genes that encode virulence factors, suggesting that the plasmid may contribute directly to bacterial virulence (Rodríguezetal., 2016). The bacterium's virulence depends on a large number of cell-associated and extracellular factors. The virulence factors play an important pathological role in the colonization, the survival of the bacteria and the invasion of tissues(Wang and Gui, 2013).

   Due to cost-efficiency throughput of genome sequencing has increased enormously resulting thousands of bacterial genomes now available and this number is increasing enormously day by day.Functional annotation of proteomes is a demanding problem(Roberts et al., 2004). A large fraction of proteins is still labeled as hypothetical protein,unknownfunction or with similar terms that imply that there is no functional indication for the ORF.Function annotation of putative uncharacterized HPs for their possible biological activity has emerged as an important focus for computational biology (Kumar et al., 2014; Loewenstein et al., 2009; Shahbaaz et al., 2013).The pUM505 sequence contained 138 complete coding regions, the majority of them encoded on the complementary DNA strand (75%), with respect to the predicted origin of replication (oriV). Most of the identified genes (46%) encode hypothetical proteins (HSPs). Proper structural and functional determination of this huge fraction (46%) is very important to reveal complete understanding of virulence mechanism in P. aeruginosa.Therefor an improved functional annotation of its proteome is of particular urgency. In present study, 12 randomly selected HPs from P.aeruginosa plasmid pUM505 have been annotated with the help of various bioinformatics resources. Moreover, two important annotated HPS have been structurally modeled and characterized.

2. MATERIALS AND METHODS

   1. Sequence retrieval and physiochemical characterization

      Randomly selected twelve HPs which contain standard number of amino acids sequences of P.aeruginosa plasmid pUM505 were retrieved from NCBI (http://www.ncbi.nlm.nih.gov/).For Physio-chemical characterization EXPAXY ProtParam (Gasteigeret al.,2005) tool were used.

   2. Functional characterization

      Functional annotation was performed using conserved domain database (CDD) (MarchlerBauer et al. 2011)Pfam (Bateman et al.,2005) and EGGNOG (Powell et al., 2012).CELLO (Yu et al., 2014) and PSORT B (Gardyet. al., 2003)were used for subcellular localization of proteins. Signal Peptide(Yuet al., 2010)and SecretomePwere used for nonclassicalsecretory pathway in proteins.

   3. 3D structure modelling

      Homology models of HPs were determined using SWISSMODEL (Kieferet al.,2009) and protein homology recognition engine Phyre2 (Kelley and Sternberg,2009).The Phyre server uses a library of known protein structures taken from the structural classification of proteins (SCOP) database and the PDB (Kelley and Sternberg 2009). The top ten scoring alignments were used to construct the three-dimensional structure of each HP.

   4. Validation of predicted model

      The stereo-chemical quality of the modeled structure for HPs was validated with verification server PDBSUM (Laskowski et al., 1993) using PROCHECK (Laskowski et al., 1996). PROCHECK validates the stereochemical quality of a protein structure by analyzing the overall structure and residue-by-residue geometry of proteins. ERRAT server (Li et al., 2007) was used for structure validation.

   5. Active site prediction

   MetaPocket 2.0(http://projects.biotec.tu-dresden.de/metapocket/help.php)was used to find out the ligand binding sites.Proteins are primarily scanned for ligands and it uses the interaction energy between the protein and a simple van der Waals probe to locate vigorously favorable binding sites (Zhang et al., 2011).

3. RESULTS AND DISCUSSION

   1. Sequence analysis and functional annotation

      Physicochemical studies of selected proteins revealed that two out of 12 proteins are stable. Most of the proteins are of acidic in nature (Supplementary Table S1).In pursuit of assigning putative function to hypotheticalproteins,available sequences and functional information from various resources has been interrogated. Six out of 12 selected hypothetical proteins have been putatively annotated usingsequence/domaininformationfromPFAM,functionalinformationfromEGGNOG,pathway(Table1).The presence of different domains in varying combinations in different proteins gives rise to th diverse repertoire of proteins found in nature. Identifying the domains present in a protein can provide insights into the function of that protein.Protein YP_004928038.1 has been annotated as PIN3 family protein,YP_004928003.1 asTraU family,YP_004927980.1 as NA-37 family, and YP_004927989.1 as HNH endonuclease (Table 1). Four proteins have been annotated with high confidant. PIN3 family protein plays avery important role and function as nuclease enzymes that cleave single stranded RNA in a sequence dependent manner. PIN domain contain three nearly invariant acidic residues. PIN-domain proteins found in prokaryotes are the toxic components of toxin-antitoxin operons. These loci provide a control mechanism that helps free-living prokaryotes cope with nutritional stress (Arcuset al., 2011). HNH proteins are involved in DNA homing, restriction, repair, or chromosome degradation.The HNH proteins are not only involved in homing but also carry other biological functions, such as DNA degradation, repair, and restriction (Walker et al., 2002). TraU is anessential protein to conjugal DNA transfer (Moore et al., 1990).All these annotated HPs play crucial role.Sub-cellular localization of protein is an important parameter for protein function in various cellular processes. Sub-cellular localization analysis revealed that 9 proteins are cytoplasmic, one outer membrane protein and two unknown (Table 2).Motifs are signatures of protein families and can be preferably used to define the protein function, particularly in enzyme where motifs are associated with the catalytic function (Bork and Koonin, 1996).

   2. 3D Structure prediction and Validation

      Three dimensional structures of two important annotated proteins have been done with the SWISS Model and Phyre2 server. Based on the results, the stereo chemical evaluation of backbone and dihedral angles of the HSPs revealed that for model HP38; 94.6, and5.4 % residues were within the most favored regions, additionally allowed regions respectively (Fig 1b). No residue was found in generouslyallowed or disallowedregions.For model HP80,84.6,12.8, 1.7 and 0.9% residues fallwithin the most favored regions, additionally allowed regions, generously allowed regions and disallowed regions respectively (Fig 2b).Overall 100 % implies beststereo-chemical quality for both the model. Based on the stereo chemical validation model HP38 showed more robust model than HP80.Secondary structure prediction showed 5 helices, 2 strands, 7 beta turns for model HP03 while 7 helices, 4 strands and 7 beta turns for HP38(Fig 3 a-b).

      ERRAT is a program for verifying protein structures determined by crystallography. ERRAT server was also used for model quality estimation which showed that overall quality factor was 77.11 and 81.32 % for model HP38 andHP80 respectively which confirms as a good model (Fig.1a- 2a).

   3. Active sites in predicted models

      Activesites on a protein is of fundamental importance for a range of applications including molecular docking, de novo drug design and structural identification and comparison of functional sites.During analysis, seven and four binding pockets were

      identified for model HP38 and HP80, respectively (Fig 4). During analysis of NABPs in HP38 model, two largest positive patches were detected viz Patch 1:'GHE', 'SFN', 'LCS', 'FPK', 'PAS', 'CON' Patch 2: 'FPK', 'PAS', 'LCS', 'CON', 'GHE'. For model HP38, two largest positive patches were detected viz Patcp 'GHE', 'SFN', 'LCS', 'FPK', 'PAS', 'CON' Patch 2: 'PAS','FPK', 'LCS','SFN', 'GHE' Patch 3: 'PAS', 'LCS', 'GHE'.Results are shown in Table 3. This data of active binding site residues will give insight into identifying binding interactions and docking with specific ligands.

4. CONCLUSIONS

Our primary sequence-based analysis led to the identification of two HPs as biologically significant, which might be involved as enzymes (antitoxins, conjugal DNA transferase, and oxido-reductase etc.). Furthermore, we successfully predicted the structure of two HPs to describe their functions at the molecular level. The outcome of the present study may facilitate better understanding of the mechanism of virulence, drug resistance, pathogenesis, adaptability to host, tolerance for host immune response and drug discovery for treatment of infections caused by P. aeruginosa.

ACKNOWLEDGEMENTS

The authors gratefully acknowledge the financial assistance under project Application of Microorganisms in Agriculture and Allied Sectors (AMAAS) from Indian Council of Agricultural Research (ICAR), India. The authors further acknowledge the Integral University, Lucknow for providing the necessary support for the completion of this study and allotting its manuscript communication ID.

Conflict of interest: The authors declare that they have no conflict of interest.

REFERENCES

1. Arcus, V.L., McKenzie, J.L., Robson, J., Cook, G.M., 2011. The PIN-domain ribonucleasesand the prokaryotic VapBC toxin-antitoxin array.Protein Eng. Des. Sel. 24 (1-2), 3340.

2. Bateman,A., Coin, L., Durbin, R., Finn, R. D., Hollich, V., GriffithsJones,S., &Studholme, D. J., 2004. The Pfam protein families database.Nucleic acids research, 32(suppl 1), D138-D141.

3. Bork, P.,Koonin, E.V.,1996.Protein sequence motifs. Current opinion instructural biology, 6(3), 366-376.

4. Gardy, J. L., Spencer, C., Wang, K., Ester, M., Tusnady, G. E., Simon, I., Brinkman, F. S., 2003. PSORT-B: Improving protein subcellular localization prediction for Gram-negative bacteria. Nucleic acids research,31(13), 3613-3617.

5. Díaz, D. A., Barreto, G. E., Santos, J. G., 2014. Structural and Functional Prediction of the Hypothetical Protein Pa2481 in Pseudomonas Aeruginosa Pao1. Advances in Computational Biology, Springer International Publishing (pp. 47-55).

6. Gasteiger, E., Hoogland, C., Gattiker, A., Duvaud, S. E., Wilkins, M. R., Appel, R. D., Bairoch, A., 2005. Protein identification and analysis tools on the ExPASy server. Humana Press,(pp. 571-607).

7. Kelley, L. A., Sternberg, M. J., 2009. Protein structure prediction on the Web: a case study using the Phyre server. Nature protocols,4(3), 363-371.

8. Kiefer, F., Arnold, K., Künzli, M., Bordoli, L., Schwede, T., 2009.The SWISS-MODEL Repository and associated resources. Nucleic acids research,37(suppl 1), D387-D392.

9. Khan, N. H., Ishii, Y., Kimata-Kino, N., Esaki, H., Nishino, T., Nishimura, M., Kogure, K., 2007. Isolation of Pseudomonas aeruginosa from open ocean and comparison with freshwater, clinical, and animal isolates. Microbial ecology, 53(2), 173-186.

10. Kumar K., Prakash A., Tasleem M., Islam A., Ahmad F., Hassan M.I., 2014. Functional annotation of putative hypothetical proteins from Candida

    dubliniensis. Gene 543,93100.

11. Laskowski, R. A.,MacArthur, M. W., Moss, D. S., Thornton, J. M., 1993.PROCHECK: a program to check the stereochemical quality of proteinstructures. Journal of applied crystallography, 26(2), 283-291.

12. Laskowski, R. A., Rullmann, J. A. C., MacArthur, M. W., Kaptein, R., Thornton, J. M., 1996. AQUA and PROCHECK-NMR: programs for checking the quality of protein structures solved by NMR. Journal of biomolecular NMR, 8(4), 477-486.

13. Loewenstein, Y., Raimondo, D., Redfern, O.C., Watson, J., Frishman D., Linial M., Orengo C., Thornton J., Tramontano A., 2009. Protein function annotation by homology-based inference. Genome Biology, 10,207.

14. Li, M., Wang, B., 2007. Homology modeling and examination of the effect of the D92E mutation on the H5N1 nonstructural protein NS1 effector domain. Journal of molecular modeling, 13(12), 1237-1244.

15. Lyczak. J. B., Cannon, C. L., Pier, G. B., 2000. Establishmentof Pseudomonas aeruginosa infection: lessons from a versatile opportunist. Microbiology Infection, 2,10511060.

16. Marchler-Bauer, A., Lu, S., Anderson, J. B., Chitsaz, F., Derbyshire, M. K., DeWeese-Scott, C.,Gwadz, M., 2001. CDD: a Conserved Domain Database for the functional annotation of proteins. Nucleic acids research,39(suppl 1), D225-D229.

17. Mayer, L. W., 1988. Use of plasmid profiles in epidemiologic surveillance of disease outbreaks and in tracing the transmission of antibiotic resistance.Clinical microbiology reviews, 1(2), 228-243.

18. Moore, D., Maneewannakul, K., Maneewannakul, S., Wu, J. H., Ippen-Ihler, K., Bradley, D. E., 1990. Characterization of the F-plasmid conjugative transfer gene traU, Journal of Bacteriology,172, 4263-4270.

19. Powell, S., Szklarczyk, D., Trachana, K., Roth, A., Kuhn, M., Muller, J., Jensen, L. J.,2012.eggNOG v3. 0: orthologous groups covering 1133 organisms at 41 different taxonomic ranges. Nucleic acids research, 40(D1), D284-D289.

20. Roberts, R. J., 2004. Identifying protein function -a call for community action. PloS 2, E42.

21. Rodríguez-Andrade, E., Hernández-Ramírez, K. C., Díaz-Peréz, S. P., Díaz-Magaña, A., Chávez-Moctezuma, M. P., Meza-Carmen, V., Ramírez-Díaz,

    M. I., 2016. Genes from pUM505 plasmid contribute to Pseudomonas aeruginosa virulence. Antonievan Leeuwenhoek, 109(3), 389-396.

22. Shahbaaz M., Hassan, M. I., Ahmad, F., 2013. Functional annotation of conserved hypothetical proteins from HaemophilusinfluenzaeRd KW20. PLoS ONE 8, e84263.

23. Wang, H., Tu, F.,Gui, Z., 2013. Virulence factors in Pseudomonas aeruginosa: mechanisms and modes of regulation. Indian Journal of Microbiology, 2, 163-167.

24. Williams, B.J., Dehnbostel, J., Blackwell, T.S., 2010. Pseudomonas aeruginosa: host defence in lung diseases. Respirology 15, 1037 1056.

25. Yu, C. S., Cheng, C. W., Su, W. C., Chang, K. C., Huang, S. W., Hwang, J. K., Lu, C. H.,2014. CELLO2GO: a web server for protein subCELlularLOcalization prediction with functional gene ontology annotation. PloS one,9(6), e99368.

26. Zhang, Z., Li, Y., Lin, B., Schroeder, M., Huang, B., 2011. Identification of cavities on protein surface using multiple computational approaches for drug binding siteprediction. Bioinformatics, 27(15), 2083-2088.

    FIGURE LEGENDS

    Fig1. (a-b) (a)Three-dimensional structure of HP38 model and (b) their stereo-chemical property by Ramchandran plot. All the residues are in most favored region

    Fig2 (a-b)(a)Three-dimensional structure of HP80 model and (b) their stereo-chemical propertyby Ramachandran plot. All the residues are in most favored region

    Fig. 3 (a-b)The secondary structure of (a) HP38 and (b) HP80 showing helices, b-Sheets and b-hairpins. Fig. 4 (a-b)Model quality estimation plot obtained by ERRAT server for (a) model HP38 and (b) model HP80. Fig4(c-d).Active sites (shown in balls) identified in protein (c) HP38 and (d) HP80.

    Fig1. (a-b):- (a)Three dimensional structure of HP38 model and (b) their stereo-chemical property by Ramchandran plot. All the residues are in most favored reg

    1. (b)

Fig2 (a-b):-(a)Three dimensional structure of HP80 model and (b) their stereo-chemical propertyby Ramachandran plot.

Fig. 3 (a-b)The secondary structure of (a) HP38 and (b) HP80 showing helices, b-Sheets and b-hairpins.

1. (b)

   Fig. 4 (a-b):-Model quality estimation plot obtained by ERRAT server for (a) model HP38 and (b) model HP80

   1. (b)

Fig4(c-d).Active sites (shown in balls) identified in protein (c) HP38 and (d) HP80

(c) (d)

Tables

Table1. Predicted Functions of HPs in P.aeruginosa plasmid pUM505

Table2.Subcellular localization of HPs predicted by different bioinformatics tool

Table 3. MetaPocket clusters and their unctional residues

Supplementary Tables

Table 1. Supplementary table S1.Physiochemical characterization of HPs protein