
 the house dust mites ( hdm ) are major sources of indoor allergens for humans , which induce asthma , rhinitis , dermatitis , and other allergic diseases .
their major allergens ( dermatophagoides pteronyssinus [ der p ] and dermatophagoides farinae [ der f ] ) coexist in most geographical regions with a high proportion ( up to 85% ) of asthmatics being typically hdm allergic ; hence , sensitization is attributed as a risk factor for developing asthma .
recently , a birth cohort study showed that sensitization to hdm at age of 2 years was associated with current wheeze at age of 12 years in both monosensitized and polysensitized hdm - sensitized children . in previous studies , fourteen d. farinae allergens
( der f 13 , 6 , 7 , 10 , 11 , 1318 , and 22 ) were reported before other seventeen allergens belonging to twelve different groups were identified by a procedure of proteomics combined with two - dimensional immunoblotting from d. farinae extracts [ 3 , 4 ] . among the novel identified d. farinae allergens , der f 25 is a triosephosphate isomerase ( tpi ) with a molecular weight of 34  kda , showing 75.6% by immunoblotting and 60% by skin prick positive reaction to dust mite allergic patients , respectively .
tpi is an enzyme ( ec 5.3.1.1 ) that catalyzes the reversible interconversion of the triose phosphate isomers dihydroxyacetone phosphate and d - glyceraldehyde 3-phosphate .
it has been found in nearly every organism searched for the enzyme , including animals such as mammals and insects as well as in fungi , plants , and bacteria .
moreover , some tpis have been identified as an allergen in fish , midges , crustaceans , and various plants [ 612 ] . currently , specific immunotherapy is the only allergen - specific approach for its treatment of mite allergy .
the administration of increasing doses of allergen extracts to patients is the method most commonly applied .
. it could induce severe anaphylactic side reactions or lead to sensitization towards new allergens present in the mixture [ 13 , 14 ] .
different strategies have been designed to try to overcome these negative effects , as the use of allergen - derived b cell peptides , allergen - derived t cell epitope containing peptides , or vaccination with allergen - encoding dna . known
epitopes for some of these mite allergens are described in detail in cui 's review . however , there is no report about the epitope of der f 25 allergen . in the present study ,
we firstly identified the b and t cell epitopes of der f 25 allergen by in silico approach .
it implied their potential utility in a peptide - based vaccine design for mite allergy .
the complete amino acid sequence of der f 25 was acquired from the nucleotide database of ncbi ( http://www.ncbi.nlm.nih.gov/ ) with the accession number of kc305500.1 .
the amino acid sequence was also used as query to search for homologous sequences through the swiss - prot / trembl ( uniprot ) ( http://www.uniprot.org/ ) and tblastn in ncbi ( http://blast.ncbi.nlm.nih.gov/blast.cgi ) .
phylogenetic tree was obtained by using ml ( maximum - likelihood ) method on the basis of the jtt amino acid sequence distance implemented in mega 5.1 ; the reliability was evaluated by the bootstrap method with 1000 replications .
the possible domains and characteristic motifs and patterns contained in der f 25 were investigated by pfam v27.0 ( http://pfam.sanger.ac.uk/ ) , prosite ( http://prosite.expasy.org/scanprosite/ ) , interpro v46.0 ( http://www.ebi.ac.uk/interpro/ ) , and superfamily v1.75 ( http://supfam.cs.bris.ac.uk/superfamily/index.html ) .
physiochemical analysis including molecular weight , theoretical pi , amino acid composition , instability index , aliphatic index , and grand average of hydropathicity ( gravy ) of der f 25 was performed by using protparam tool ( http://web.expasy.org/protparam/ ) .
der f 25 characteristic pattern was checked for original sequence and further analysis was performed to highlight the presence of functional motifs by using the prosite database ( http://prosite.expasy.org/ ) . biologically meaningful motifs and susceptibility to posttranslational modifications
phosphorylation motifs with more than 80% of probability of occurrence were analyzed by using netphos v2.0 ( http://www.cbs.dtu.dk/services/netphos/ ) and netphosk v1.0 ( http://www.cbs.dtu.dk/services/netphosk/ ) .
der f 25 secondary structural elements were predicted by psipred ( http://bioinf.cs.ucl.ac.uk/psipred/ ) , which threads sequence segments through protein data bank ( pdb ) library ( http://www.rcsb.org/ ) to identify conserved substructures .
furthermore , the secondary structure elements were also identified with the result obtained with netsurfp ver .
as well , the homologous templates suitable for der f 25 were selected by psi - blast server ( http://blast.ncbi.nlm.nih.gov/blast.cgi ) and swiss - model server ( http://swissmodel.expasy.org/ ) [ 25 , 26 ] . the best template was retrieved from the results of previous methods and used for homology modeling .
der f 25 modeled protein structure was built through alignment mode in swiss - model using the complete amino acid sequence .
an initial structural model was generated and checked for recognition of errors in 3d structure by procheck , errat , and verify 3d   programs in structural analysis and verification server ( saves ) ( http://nihserver.mbi.ucla.edu/saves/ ) .
the final model structure quality of der f 25 was assessed by qmean , by checking protein stereology with prosa program   and the protein energy with anolea ( http://protein.bio.puc.cl/cardex/servers/anolea/ ) .
the ramachandran plot for all the models was generated , showing the majority of the protein residues in the favored regions .
der f 25 model was submitted to consurf server ( http://consurf.tau.ac.il/ ) in order to generate evolutionary related conservation scores helping to identify functional regions in the proteins .
functional and structural key residues in der f 25 sequence were confirmed by conseq server .
apbs molecular modeling software implemented in pymol 0.99 was used to investigate the electrostatic poisson - boltzmann ( pb ) potentials of der f 25 model structure .
amber99 in pdb2pqr server ( http://nbcr-222.ucsd.edu/pdb2pqr_1.8/ ) was used to assign the charges and radii to all of the atoms ( including hydrogens ) .
fine grid spaces of 0.35   were used to solve the linearized pb equation in sequential focusing multigrid calculations in a mesh of 130 points per dimension at 310.00  k. the dielectric constants were 2.0 and 80.0 for the protein and water .
the output mesh was processed in the scalar opendx format to render isocontours and maps onto the surfaces with pymol 0.99 .
potential values are given in units of kt per unit charge ( k boltzmann 's constant ; t temperature ) .
three immunoinformatics tools including dnastar protean system , bioinformatics predicted antigenic peptides ( bpap ) system ( http://imed.med.ucm.es/tools/antigenic.pl ) , and bepipred 1.0 server ( http://www.cbs.dtu.dk/services/bepipred/ ) were used to predicate the b cell epitopes of der f 25 .
the ultimate consensus epitope results were obtained by combining the results of the three tools together with the method published earlier . in the dnastar protean system , four properties ( hydrophilicity , flexibility , accessibility , and antigenicity ) of the amino acid sequence were chosen as parameters for epitopes prediction .
the bpap system and the bepipred 1.0 server only need the amino acid sequence and provide more straightforward results which are combined with physicochemical properties of amino acids such as hydrophilicity , flexibility , accessibility , turns , and exposed surface .
t cell epitopes are principally predicted indirectly by identifying the binding of peptide fragments to the mhc complexes .
the binding significance of each peptide to the given mhc molecule is based on the estimated strength of binding exhibited by a predicted nested core peptide at a set threshold level . for hla - dr - based t cell epitope prediction , the artificial neural network - based alignment ( nn - align ) method netmhciipan-3.0 ( http://www.cbs.dtu.dk/services/netmhciipan/ )   was applied . for hla - dq alleles , netmhcii-2.2 ( http://www.cbs.dtu.dk/services/netmhcii/ )  
hla - dr 101 , hla - dr 301 , hla - dr 401 , and hla - dr 501 were used to predict hla - dr - based t cell epitope prediction .
the ultimate hla - dr - based t cell epitope results were obtained by combining those four results together that if three of them showed epitope , then the consensus result was epitope .
this method was also used in hla - dq - based t cell epitope prediction .
hla - dqa10101-dqb10501 , hla - dqa10501-dqb10201 , hla - dqa10501-dqb10301 , and hla - dqa10102-dqb10602 were used to predict hla - dq - based t cell epitope prediction . as a result
, the ultimate consensus epitope results were obtained by combining the results of the hla - dr - based t cell epitope and hla - dq - based t cell epitope .
b cell and t cell epitopes identified by computational tools were mapped onto linear sequence and on the three - dimensional model of der f 25 to determine their position and secondary structure elements involved .
the amino acid sequence of der f 25 was obtained from the nucleotide database of ncbi .
uniprot and tblastn were used to search the homologous sequences of der f 25 . as a result ,
thirty - six sequences were obtained and in order to determine the relationships between der f 25 and its homologous sequences , phylogenetic analysis was performed and the evolutionary tree inferred by the ml method was showed in figure 1 .
phylogenetic analysis result showed that there are proteins including der f 25 clustered into the same group , belonging to tips .
moreover , domain analysis results showed that der f 25 belongs to the tim phosphate binding superfamily ( superfamily number ssf51351 and interpro number ipr016040 ) and tpi family ( superfamily number ssf51352 and interpro number ipr000652 ) . after searching for characteristic motifs or patterns
, we found that der f 25 exhibited a tpi pattern , ps00171 ( 162172 , ayepvwaigtg ) ( figure 2 ) .
phosphorylation sites including two ser ( 95 and 221 ) and two thr ( 146 and 171 ) residues were predicted and showed in figure 2 .
two types of kinases ( pkc for 95 , 146 , and 171 and dnapk for 221 ) were predicted to be phosphorylated for der f 25 complete sequence .
the primary structure of der f 25 contained 247 amino acids and the molecular weight is 27134.1 .
the theoretical pi is 6.24 and the aliphatic index is 95.06 . the gravy is 0.103 meaning that der f 25 exhibited hydrophilic character .
the instability index is 30.57 meaning that the sequence of der f 25 is stable .
searching for the proteins with known tertiary structure in the pdb yielded tenebrio molitor tpi ( pdb accession number : 2i9e ) showing the highest sequence identity ( 74% ) with der f 25 .
the swiss - model server was also used to identify the best possible template and found a high score of 365 and very low e - value of e-101 for 2i9e template .
hence , the 2i9e template was used for homology modeling . as indicated by the ramachandran plot ( figure 3(b ) ) , 93% residues in der f 25 model were within the most favored regions , 7% residues in the additional allowed region , 0% residues in the generously allowed regions , and 0% residues in the disallowed region ; 93.4% residues in 2i9e template were within the most favored regions , 6.6% residues in the additional allowed region , 0% residues in the generously allowed regions , and 0% residues in the disallowed region . the goodness factor ( g - factor ) based on the observed distribution of stereochemical parameters ( main chain bond angles , bond length , and phi - psi torsion angles ) returned accurate values for a reliable model ( table 1 ) , in comparison with the template 2i9e ( table 1 ) .
as indicated by the errat program , the result showed that the overall quality factor is 97.034 for der f 25 and 97.789 for 2i9e meaning that both of the two structures have good high resolution .
as indicated by the verify 3d program , the result showed that 99.6% and 99.8% of der f 25 and 2i9e template residues had an average 3d ( atomic model)-1d ( amino acid sequence ) score > 0.2 also meaning that those two structures were good .
prosa analysis returned z - scores of 10.1 and 10.24 for der f 25 and 2i9e ( table 1 ) , respectively .
q values for der f 25 and 2i9e structures are 0.772 and 0.798 , respectively .
root mean square deviations ( rmsd ) between der f 25 structural model and 2i9e template c backbone are 0.062   ( table 1 ) . based on these validations
secondary structure prediction of der f 25 with psipred identified ten -helices and seven -sheets ( figure 2 ) in der f 25 .
the best template 2i9e was used for homology modeling ; the overall 3d structure of der f 25 was shown in figure 3(a ) .
sequence polymorphism was responsible for the changes in the spatial distribution of the skeleton alpha carbons , which is reflected in differences between the structures of der f 25 and 2i9e .
a superposition of the der f 25 with the 2i9e template is shown in figure 3(b ) and the values for superimposed c are 0.062   . as a tpi protein , der f 25 has two active sites ; his in the 94th position is an electrophile while glu in the 164th position is the proton acceptor .
it has two substrate binding sites , the asn in the 10th position and the lys in the 12th position .
moreover , the characteristic pattern predicted by scanprosite tool is shown in figure 3(c ) .
consurf conservational analysis of structural and functional key amino acids showed that the der f 25 protein surfaces were not well conserved , with almost forty high variability residues in different superficial areas .
surface electrostatic potential analysis reveals several prominent charged residues , with half of the side exhibiting large positive values ( blue regions ) and the other half showing predominantly negative values ( red regions ) ( figure 3(d ) ) .
in addition , the existence of regions with high hydrophobicity also provides strong evidence for epitope identification .
based on these sequence properties , the final predicting regions of dnastar were 1119 , 2835 , 6877 , 95111 , 130139 , 173187 , 193197 , and 211225 . also , the predicted results of the bpap system were 2026 , 3252 , 5470 , 8693 , 138168 , 175209 , and 218242 and bepipred 1.0 server were 1118 , 3034 , 7178 , 99107 , 132138 , 168183 , 194199 , and 210224 .
furthermore , the final potential b cell epitopes of der f 25 were selected on the basis of the results of these three tools .
the ultimate results of the three immunoinformatics tools finally predicted eight peptides ( 1118 , 3035 , 7177 , 99107 , 132138 , 173187 , 193197 , and 211224 ) ( table 2 ) and these peptides were also shown in figures 2 and 4 .
netmhciipan-3.0 and netmhcii-2.2 were used to identify the t cell epitope of der f 25 . for hla - dr - based t cell epitope prediction ,
the final predicting regions of hla - dr 101 , hla - dr 301 , hla - dr 401 , and hla - dr 501 were shown in table 1 and the ultimate results of hla - dr - based t cell epitope prediction finally predicted five peptides ( 2634 , 3854 , 6674 , 142151 , and 239247 ) .
for hla - dq alleles , the final results of hla - dqa10101-dqb10501 , hla - dqa10501-dqb10201 , hla - dqa10501-dqb10301 , and hla - dqa10102-dqb10602 were also shown in table 3 and the ultimate results of these four methods finally predicted one peptide , 3948 . as a result , der f 25 was predicted to have five t cell epitope sequences including 2634 , 3854 , 6674 , 142151 , and 239247 ( table 2 ) as shown in figures 2 and 4 .
the prevalence of human atopic disorders including allergic rhinitis , asthma , and atopic dermatitis is increasing during the past several decades .
house dust mite allergies constitute more than 50% of allergic patients and often have severe forms of respiratory allergy , such as asthma .
characterization of mite allergens will be beneficial in the diagnosis and treatment of mite - induced atopic illnesses . among the identified allergens ,
der f 25 is a new protein with a molecular weight of 34  kda , representing the major allergen in d. farinae .
the objective of this study was to predict the b and t cell epitopes of der f 25 .
firstly , in order to better understand the structure and function of der f 25 , we analyzed the basic sequence properties and studied the 2d and 3d structures of der f 25 .
phylogenetic analysis result showed that der f 25 protein clustered into the tips group ; domain analysis also proved a strong evidence illustrating that der f 25 belongs to the tpi family . in 2d structure analysis
, it is clearly shown that der f 25 composed of ten -helices and seven -sheets .
the 3d structure of der f 25 was performed by homology modelling which was widely used in many areas of structure - based analysis and study .
pdb server was used to search templates of der f 25 and found that the structure of tenebrio molitor tpi ( 2i9e ) was the best template with the highest identity . also , small dissimilates in rmsd were observed between der f 25 and 2i9e template .
the built model structure is feasible by the ramachandran plot analysis , errat program , verify 3d program prosa analysis , q values , and rmsd .
the similar methods for the structures modelling were also successfully conducted in other allergens of api si and api sii , der f 5 , ole e 2   ole e 11 , and ole e 12 .
based on the conservational analysis of the primary sequence of der f 25 , it is found that almost forty high variability residues sit in different superficial areas of protein surface .
the active sites his in 94th position and glu in 164th position as well as the substrate binding sites asn in 10th position and lys in 12th position are the completely conserved sites . 
 in silico prediction has already become a familiar and useful tool for selecting epitopes from immunological relevant proteins , which can save the expense of synthetic peptides and the working time .
recently , many algorithms have been developed to predict b cell epitopes on a protein sequence based on propensity values of amino acid properties of hydrophilicity , antigenicity , segmental mobility , flexibility , and accessibility . in the present study , we used three algorithms ( dnastar protean system , bpap , and
the previous study showed that the use of bioinformatics approach to predict b cell epitopes correlated well with the experimental approach .
earlier study showed that allergen epitopes were comprised of a high proportion of hydrophobic amino acids .
the amino acids ala , ser , asn , gly , and particularly lys play a key role in the ige binding allergenic epitopes . in our results ,
nearly half of the total residues lying in b cell epitopes were hydrophobic ( table 2 ) .
moreover , each predicted b cell epitope has one or more special five amino acids and the common residues in all b cell epitopes were gly and lys ( table 2 ) .
electrostatic interactions are known to determine the orientation of the molecules and stabilize antigen - antibody complexes .
surface electrostatic potential analysis result showed that a great part of der f 25 side exhibits large positive values ( blue regions ) .
most parts of b cell epitopes are distributed in the blue regions and showed a strong negative potential . as a result ,
eight peptides ( 1118 , 3035 , 7177 , 99107 , 132138 , 173187 , 193197 , and 211224 ) were predicted as the b cell epitopes .
however , these b cell epitopes need further investigation in clinical samples . in the last several years
, some methods have substantially improved their accuracy to predict t cell epitopes such as netmhcpan-3.0 and netmhcii-2.2 .
netmhcpan-3.0 is based on artificial neural networks and is trained on 52,062 quantitative peptide binding data covering all hla as well as two mouse molecules . in this study , it was used to predict the hla - dr - based t cell epitopes . for hla - dq - based t cell epitopes prediction ,
although limited binding - affinity data are available for hla - dq , it was recently reported to provide the best performance in predicting this locus . as a result , netmhciipan-3.0 and netmhcii-2.2
were used to predict the t cell epitopes in der f 25 allergens and predicted 7 potential t cell epitope sequences including 2634 , 3854 , 6674 , 142151 , and 239247 . despite the high accuracy of these predictions , this approach has not yet been applied to peptide - based vaccine development for allergic diseases .
allergen - specific immunotherapy ( sit ) represents the only allergen - specific and disease - modifying approach with long lasting effects for the treatment of allergic patients .
however , sit can induce side effects , ranging from mild and local to severe and life - threatening symptoms , such as anaphylactic shock .
severe side effects are frequently observed in patients with house dust mite ( hdm ) allergy . the continuous exposure to hdm
allergens further complicates the treatment of patients with hdm allergy . additionally , the quality of natural hdm allergen extracts and vaccines based on these extracts is often poor .
attenuated allergenic molecules , that is , hypoallergens or synthetic peptide fragments , have been used as high dose and safer alternatives to conventional extract - based sit .
vaccination with a combination of small peptides that together extend across the entire native allergenic protein theoretically could preserve t cell activation while avoiding ige - based immune responses .
ige recognizes conformational epitopes of larger peptides ( b cell epitopes ) and proteins while t cell receptors recognize small linear peptides of 8 to 10 amino acids ( t cell epitope ) . by immunizing with small peptides ,
t cell activation could occur while ige binding would be lost [ 55 , 56 ] .
then , we predicted b and t cell epitopes of der f 25 allergen , the major allergen in hdm , using in silico method which can be used to benefit allergen immunotherapies and reduce the frequency of allergic reactions .
in this study , we have a better understanding of the 2d and 3d structures of der f 25 and have predicted eight b cell epitopes ( 1118 , 3035 , 7177 , 99107 , 132138 , 173187 , 193197 , and 211224 ) and five t cell epitope including 2634 , 3854 , 6674 , 142151 , and 239247 of this tpi .
all these results can be used to benefit allergen immunotherapies and reduce the frequency of mite allergic reactions .
 the house dust mites are major sources of indoor allergens for humans , which induce asthma , rhinitis , dermatitis , and other allergic diseases . 
 der f 25 is a triosephosphate isomerase , representing the major allergen identified in dermatophagoides farinae . 
 the objective of this study was to predict the b and t cell epitopes of der f 25 . in the present study 
 , we analyzed the physiochemical properties , function motifs and domains , and structural - based detailed features of der f 25 and predicted the b cell linear epitopes of der f 25 by dnastar protean system , bpap , and bepipred 1.0 server and the t cell epitopes by netmhciipan-3.0 and netmhcii-2.2 . as a result , the sequence and structure analysis identified that der f 25 belongs to the triosephosphate isomerase family and exhibited a triosephosphate isomerase pattern ( ps001371 ) . 
 eight b cell epitopes ( 1118 , 3035 , 7177 , 99107 , 132138 , 173187 , 193197 , and 211224 ) and five t cell epitopes including 2634 , 3854 , 6674 , 142151 , and 239247 were predicted in this study . these results can be used to benefit allergen immunotherapies and reduce the frequency of mite allergic reactions . 
1. Introduction
2. Methodology
3. Results
4. Discussion
5. Conclusion