bioinformatics analysis of the structure and function of NADPH‑cytochrome p450 reductase of Plasmodium vivax

  • Authors:
    • Zhigang Fan
    • Kaijie Li
    • Lingmin Zhang
    • Fan Chen
    • Qiang Wu
    • Na Li
    • Saifeng Zhong
    • Guifen Lin
    • Guogang Yan
  • View Affiliations

  • Published online on: March 4, 2013     https://doi.org/10.3892/br.2013.71
  • Pages: 425-427
Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )
Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )


Abstract

The structure of NADPH‑cytochrome p450 reductase (CPR) of Plasmodium falciparum (P. falciparum or Pf) has been determined using bioinformatics analysis. However, that of Plasmodium vivax (P. vivax or Pv) has not yet been determined. This study aimed to analyze the structure and function of PvCPR using bioinformatics analysis. The results demonstrated that PvCPR was an unstable and alkaline enzyme located in the cytoplasm of parasites with a signal peptide. It possessed seven types of signal sites and eight protein‑protein binding sites, and had a tertiary structure resembling a forceps with a single wing, which differed from that of PfCPR. It also had nine linear B‑cell epitopes and 10 antigenicity sites, which were not homologous with the amino acid sequence of Homo sapiens (H. sapiens or Hs) CPR and six fragments that were similar to fragments of immune‑related protein sequences from H. sapiens. Therefore, the function of PvCPR may be different from that of PfCPR, and PvCPR may participate in the immune escape of P. vivax.

Introduction

Plasmodium vivax (P. vivax or Pv) malaria epidemics occur mainly on the American and Asian continents (1). P. vivax is responsible for 74% of cases on the American continent (2) and for >90% of the cases in China (3). The function of NADPH-cytochrome p450 reductase (CPR) in mammals involves electron transport. The structure and function of Plasmodium falciparum (P. falciparum or Pf) CPR have been analyzed using bioinformatics (4), while those of PvCPR have not. Therefore, the aim of this study was to investigate the structure and function of PvCPR using bioinformatics analysis.

Materials and methods

The methods used to study the amino acid (aa) sequence, physicochemical properties, subcellular localization, secondary protein structure and 3D structure have been previously described (4). The PvCPR signal peptide was predicted using the Signal P program (http://www.cbs.dtu.dk/services/SignalP/). Identification and comparison of similar sequences were performed using BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi). The location of the linear B-cell epitopes and antigenicity sites were predicted with http://www.cbs.dtu.dk/services/BepiPred/ and http://tools.immuneepitope.org/tools/bcell/iedb_input, respectively. Identification and comparison of similar sequences were performed with http://blast.ncbi.nlm.nih.gov/Blast.cgi.

Results

Physical and chemical characteristics of PvCPR

The aa sequence of PvCPR (NCBI reference sequence: XP_001614764.1) contained 746 aa with a molecular weight of 85.175 kDa and a predicted isoelectric point of 8.72. The estimated half-lives of the protein in mammals, yeast and E. coli were 30, >20 and >10 h, respectively, and its instability coefficient was 42.89, which was characteristic of an unstable protein. The grand average of hydrophobicity of PvCPR was −0.375. Therefore, PvCPR was an unstable and alkaline enzyme.

Subcellular localization of PvCPR

The probability of localization of PvCPR was predicted to be 43.5% in the mitochondria, 17.4% in the cytoplasm, 13.0% in the vacuole, 4.3% in the Golgi apparatus, 4.3% in the endoplasmic reticulum, 4.3% in the nucleus and 13.0% outside the cell. PvCPR was predicted to localize in the mitochondria of P. vivax by the k-NN prediction; however, the reliability of PvCPR location in the cytoplasm was 76.7 according to Reinhardt’s method for cytoplasmic/nuclear discrimination. PvCPR was predicted to consist of a signal peptide (probability of 0.719). The likely cleavage site was predicted to be located between aa 23 and 24 or between aa 24 and 25.

Molecular evolutionary conserved domains of PvCPR

PvCPR was most closely related to P. knowlesi (Pk) CPR (NCBI reference sequence: XP_002258503.1) with 82% identity and 89% similarity, but not to P. falciparum(4). PvCPR contained one complete FMN_red superfamily and one incomplete FNR_like superfamily, similar to PfCPR.

Structure of PvCPR

The secondary structure of PvCPR comprised 31.77% α-helices, 15.15% β-strands and 53.08% loops. In addition, this protein possessed no transmembrane region. Approximately 50% of the aa were exposed on the surface of the protein. Five N-glycosylation sites, three cAMP- and cGMP-dependent protein kinase phosphorylation sites, 12 protein kinase C phosphorylation sites, eight casein kinase II phosphorylation sites, four tyrosine kinase phosphorylation sites, seven N-myristoylation sites, one amidation site and eight protein-protein binding sites were identified in PvCPR. Fifteen functional sites and six protein-protein binding sites were located in the FNR_like superfamily region and one protein-protein binding site was located in the nuclear organizer region (NOR), which was predicted to span between aa 203–320.

The tertiary structure of the protein (aa 45–746) resembled a forceps with a single wing (Fig. 1). The FMN_red super-family and FNR_like superfamily formed a large and deep indentation. A section of the NOR (aa 253–290) formed a single wing. Of note, most of the functional sites were scattered on the surface of the protein.

Linear B-cell epitopes and antigenicity sites of PvCPR

PvCPR was located in the cytoplasm of parasites and belonged to a secretory protein. It was not known whether PvCPR was an ideal target for diagnosis, vaccines or drugs. Therefore, the linear B-cell epitopes, antigenicity sites and unique sites of PvCPR were analyzed. There were 26 linear B-cell epitopes and 32 antigenicity sites in the aa sequence of PvCPR (Table I). Nine linear B-cell epitopes were not homologous with the aa sequence of Homo sapiens (H. sapiens or Hs) CPR, as follows: aa 1–3, aa 250–251, aa 253–290, aa 314–322, aa 347–349, aa 593, aa 630, aa 725–730 and aa 736–737. Epitopes aa 137–142, aa 376–416, aa 571–579 and aa 654–662 were less homologous with those of HsCPR compared to the remaining linear B-cell epitopes. The following 10 antigenicity sites were not homologous with the aa sequence of HsCPR: aa 6–23, aa 28–38, aa 40–46, aa 204–234, aa 236–244, aa 246–254, aa 281–287, aa 323–335, aa 391–397 and aa 563–571. In the tertiary protein structure, eight linear B-cell epitopes (between aa 1–746) were located on the surface and one antigenicity site (aa 563–571), was located on the inner side of the surface, whereas six antigenicity sites (located between aa 45–746) were located on the surface.

Table I

Linear B-cell epitopes and antigenicity sites of PvCPR.

Table I

Linear B-cell epitopes and antigenicity sites of PvCPR.

Linear B-cell epitopes (aa sequences)
Antigenicity sites (aa sequences)
1–349–5867–776–2328–3840–4660–66
117–125137–142183–19388–99106–117125–136145–152
250–251253–290314–322155–175204–234236–244246–254
347–349360–367376–416281–287323–335337–345369–375
498–515525541–549391–397401–443451–456463–468
571–577593605–608474–480487–508517–536544–549
630654–662673–675551–558563–571579–596607–617
708–712714–715725–730621–628634–641663–672675–706
736–737740–741

[i] PvCPR, Plasmodium vivax NADPH-cytochrome p450 reductase.

Unique fragments among the aa sequence of PvCPR

The aa fragments of PvCPR which had functional sites and no homologous fragments with the aa sequence of HsCPR inclu ded aa 1–51, aa 135–138, aa 140–147, aa 205–357, aa 366–370, aa 373–407, aa 555–577 and aa 721–730. There was no unique aa fragment that was non-homologous with other protein sequences from H. sapiens. Six fragments were similar to fragments of immune-related protein sequences from H. sapiens. These immune-related proteins included immunoglobulin heavy constant γ 1, Toll/interleukin-1 receptor-like protein 3, toll-like receptor 5 (TLR5), interleukin-4 receptor (IL-4R) and T-cell receptor β variable 3.

Discussion

The physical and chemical characteristics of PvCPR were similar to those of PfCPR (4). PvCPR was located in the cytoplasm of parasites and contained a signal peptide, which differed from that of PfCPR. The tertiary structure of PvCPR was different from that of PfCPR. The changes in the protein structure during evolution may have led to loss or gain of certain functions (5). PvCPR, as opposed to PfCPR, possessed three cAMP- and cGMP-dependent protein kinase phosphorylation sites and one amidation site. However, it did not exhibit vacuole target and bileucine motifs, in contrast to PfCPR. Therefore, the function of PvCPR may differ from that of PfCPR, which may be responsible for the differences in physiology and pathogenicity between P. vivax and P. falciparum.

PvCPR had nine linear B-cell epitopes and 10 antigenicity sites, which were not homologous with the aa sequence of HsCPR. Assuming that PvCPR is the key protein of P. vivax, these sites may be ideal targets for vaccines and drugs. However, if this is not the case, these sites may be ideal diagnostic markers of P. vivax. In addition, six fragments of PvCPR were similar to protein fragments from H. sapiens, which participate in the innate and adaptive immune responses, such as IL-4R and TLR5 (6,7). Therefore, PvCPR may be involved in the immune escape of P. vivax.

Abbreviations:

P. vivax or Pv

Plasmodium vivax

P. falciparum or Pf

Plasmodium falciparum

H. sapiens or Hs

Homo sapiens

CPR

NADPH-cytochrome p450 reductase

NOR

nuclear organizer region

TLR

toll-like receptor

Acknowledgements

The authors thank colleagues at the School of Tropical and Laboratory Medicine and Hainan Provincial Key Laboratory of Tropical Medicine and the teachers from the Department of Immunology of Tongji Medical College in Huazhong University of Science and Technology. This study was supported by Hainan Natural Science Fund (no. 812199) and Scientific Research Funds of Hainan Medical University in 2011 (no. 2010-014).

References

1 

Mueller I, Galinski MR, Baird JK, et al: Key gaps in the knowledge of Plasmodium vivax, a neglected human malaria parasite. Lancet Infect Dis. 9:555–666. 2009.PubMed/NCBI

2 

Valencia SH, Rodríguez DC, Acero DL, et al: Platform for Plasmodium vivax vaccine discovery and development. Mem Inst Oswaldo Cruz. 106(Suppl 1): 179–192. 2011.

3 

Zhou SS, Wang Y and Xia ZG: Malaria situation in the People’s Republic Of China in 2009. Chin J Parasitol Parasit Dis. 29:1–3. 2011.

4 

Fan Z, Zhang L, Yan G, et al: Bioinformatics analysis for structure and function of CPR of Plasmodium falciparum. Asian Pac J Trop Med. 4:85–87. 2011. View Article : Google Scholar : PubMed/NCBI

5 

Bridgham JT, Ortlund EA and Thomton JW: An epistatic ratchet constrains the direction of glucocorticoid receptor evolution. Nature. 461:515–519. 2009. View Article : Google Scholar : PubMed/NCBI

6 

Burt BM, Bader A, Winter D, et al: Expression of interleukin-4 receptor alpha in human pleural mesothelioma is associated with poor survival and promotion of tumor inflammation. Clin Cancer Res. 18:1568–1577. 2012. View Article : Google Scholar : PubMed/NCBI

7 

Descamps D, Le Gars M, Balloy V, et al: Toll-like receptor 5 (TLR5), IL-1β secretion, and asparagine endopeptidase are critical factors for alveolar macrophage phagocytosis and bacterial killing. Proc Natl Acad Sci USA. 109:1619–1624. 2012.

Related Articles

Journal Cover

May-June 2013
Volume 1 Issue 3

Print ISSN: 2049-9434
Online ISSN:2049-9442

Sign up for eToc alerts

Recommend to Library

Copy and paste a formatted citation
x
Spandidos Publications style
Fan Z, Li K, Zhang L, Chen F, Wu Q, Li N, Zhong S, Lin G and Yan G: bioinformatics analysis of the structure and function of NADPH‑cytochrome p450 reductase of Plasmodium vivax. Biomed Rep 1: 425-427, 2013.
APA
Fan, Z., Li, K., Zhang, L., Chen, F., Wu, Q., Li, N. ... Yan, G. (2013). bioinformatics analysis of the structure and function of NADPH‑cytochrome p450 reductase of Plasmodium vivax. Biomedical Reports, 1, 425-427. https://doi.org/10.3892/br.2013.71
MLA
Fan, Z., Li, K., Zhang, L., Chen, F., Wu, Q., Li, N., Zhong, S., Lin, G., Yan, G."bioinformatics analysis of the structure and function of NADPH‑cytochrome p450 reductase of Plasmodium vivax". Biomedical Reports 1.3 (2013): 425-427.
Chicago
Fan, Z., Li, K., Zhang, L., Chen, F., Wu, Q., Li, N., Zhong, S., Lin, G., Yan, G."bioinformatics analysis of the structure and function of NADPH‑cytochrome p450 reductase of Plasmodium vivax". Biomedical Reports 1, no. 3 (2013): 425-427. https://doi.org/10.3892/br.2013.71