Structural diversity of class I MHC-like molecules and its implications in binding specificities

Zinc-a2-Glycoprotein Acts as a Component of PNN to Protect Hippocampal Neurons from Apoptosis

In the adult mouse brain, perineuronal net (PNN), a highly structured extracellular matrix, surrounds subsets of neurons. The AZGP1 gene encodes zinc-2-glycoprotein (ZAG) is a lipid-mobilizing factor. However, its expression and distribution in the adult brain have been controversial. Here, for the first time, we demonstrate that the secreted ZAG is localized to Wisteria floribunda agglutinin (WFA)-positive PNNs around parvalbumin (PV)-expressing interneurons in the hippocampus, cortex, and a number of other PNN-bearing neurons and co-localizes with aggrecan, one of the components of PNNs. Few ZAG-positive nets were seen in the area without WFA staining by chondroitinase ABC (ChABC) which degrades glycosaminoglycans (GAGs) from the chondroitin sulfate proteoglycans (CSPGs) in the PNN. Reanalysis of single-cell sequencing data revealed that ZAG mRNA was mainly expressed in oligodendrocyte lineages, specifically in olfactory sheathing cells. The ZAG receptor β3 adrenergic receptor (β3AR) is also selectively co-localized with PV interneurons and CA2 pyramidal neurons in the hippocampus. In addition, molecular docking provides valuable new insights on how GAGs interfere with ZAG and ZAG/β3AR complex. Finally, our results indicated that human recombinant ZAG could significantly inhibit serum derivation-induced cell apoptosis in HT22 cells. Our combined experimental and theoretical approach raises a unique hypothesis namely that ZAG may be a crucial functional attribute of PNNs in the brain to protect neuronal cell from apoptosis.

Insights into the Antibacterial Activity of Prolactin-Inducible Protein against the Standard and Environmental MDR Bacterial Strains

Background: Prolactin inducible protein (PIP) is a small secretary glycoprotein present in most biological fluids and contributes to various cellular functions, including cell growth, fertility, antitumor, and antifungal activities. Objectives: The present study evaluated the antibacterial activities of recombinant PIP against multiple broad-spectrum MDR bacterial strains. Methods: The PIP gene was cloned, expressed and purified using affinity chromatography. Disk diffusion, broth microdilution, and growth kinetic assays were used to determine the antibacterial activities of PIP. Results: Disk diffusion assay showed that PIP has a minimum and maximum zone of inhibition against E. coli and P. aeruginosa, respectively, compared to the reference drug ampicillin. Furthermore, growth kinetics studies also suggested that PIP significantly inhibited the growth of E. coli and P. aeruginosa. The minimum inhibitory concentration of PIP was 32 µg/mL for E. coli (443), a standard bacterial strain, and 64 µg/mL for Bacillus sp. (LG1), an environmental multidrug-resistant (MDR) strain. The synergistic studies of PIP with ampicillin showed better efficacies towards selected bacterial strains having MDR properties. Conclusion: Our findings suggest that PIP has a broad range of antibacterial activities with important implications in alleviating MDR problems.

Mapping the gene network landscape of Alzheimer's disease through integrating genomics and transcriptomics

Integration of multi-omics data with molecular interaction networks enables elucidation of the pathophysiology of Alzheimer's disease (AD). Using the latest genome-wide association studies (GWAS) including proxy cases and the STRING interactome, we identified an AD network of 142 risk genes and 646 network-proximal genes, many of which were linked to synaptic functions annotated by mouse knockout data. The proximal genes were confirmed to be enriched in a replication GWAS of autopsy-documented cases. By integrating the AD gene network with transcriptomic data of AD and healthy temporal cortices, we identified 17 gene clusters of pathways, such as up-regulated complement activation and lipid metabolism, down-regulated cholinergic activity, and dysregulated RNA metabolism and proteostasis. The relationships among these pathways were further organized by a hierarchy of the AD network pinpointing major parent nodes in graph structure including endocytosis and immune reaction. Control analyses were performed using transcriptomics from cerebellum and a brain-specific interactome. Further integration with cell-specific RNA sequencing data demonstrated genes in our clusters of immunoregulation and complement activation were highly expressed in microglia.

In-silico design of peptides for inhibition of HLA-A*03-KLIETYFSK complex as a new drug design for treatment of multiples sclerosis disease

Multiple sclerosis is recognized as a chronic inflammatory disease. Human leukocyte antigen (HLA) plays an important role in initiating adaptive immune responses. HLA class I is present in almost all nucleated cells and presents the cleaved endogenous peptide antigens to cytotoxic T cells. HLA-A*03 is one of the HLA class I alleles, which is reported as substantially related HLA to MS disease. In 2011, the structure of the HLA-A*03 in complex was identified with an immunodominant proteolipid protein (PLP) epitope (KLIETYFSK). This complex has been reported as an important autoantigen-presenting complex in MS pathogenesis. In this study, new peptides were designed to bind to this complex that may prevent specific pathogenic cytotoxic T cell binding to this autoantigen-presenting complex and CNS demyelination. Herein, 14 new helical peptides containing 19 amino acids were designed and their structures were predicted using the PEP-FOLD server. The binding of each designed peptide to the mentioned complex was then performed. A mutation approach was used by the BeAtMuSiC server to improve the binding affinity of the designed peptide. In each position, amino acid substitutions leading to an increase in the binding affinity of the peptide to the mentioned complex were determined. Finally, the resulting complexes were simulated for 40 ns using AMBER18 software. The results revealed that out of 14 designed peptides, "WRYWWKDWAKQFRQFYRWF" peptide exhibited the highest affinity for binding to the mentioned complex. This peptide can be considered as a potential drug to control multiple sclerosis disease in patients carrying the HLA-A*03 allele.

TTP-mediated regulation of mRNA stability in immune cells contributes to adaptive immunity, immune tolerance and clinical applications

Dendritic cells (DCs) form a sentinel network to induce protective immunity against pathogens or self-tolerance. mRNA stability is an important part of the post-transcriptional regulation (PTR) that controls the maturation and function of DCs. In this review, we summarize the effects of TTP-mediated regulation of mRNA stability in DCs, focusing on DC maturation and antigen presentation, T cell activation and differentiation, immune tolerance and inflammation. We also discuss the potential DC-based immune treatment for HIV⁺ patients through regulation of mRNA stability. This review proposes the regulation of mRNA stability as a novel immune therapy for various inflammatory diseases, such as arthritis and dermatitis.

Bovine Coronavirus and the Associated Diseases

Coronaviruses (CoVs) possess the largest and most complex RNA genome (up to 32 kb) that encodes for 16 non-structural proteins regulating RNA synthesis and modification. Coronaviruses are known to infect a wide range of mammalian and avian species causing remarkably diverse disease syndromes. Variable tissue tropism and the ability to easily cross interspecies barriers are the well-known characteristics of certain CoVs. The 21st century epidemics of severe acute respiratory CoV (SARS-CoV), Middle East respiratory CoV and the ongoing SARS-CoV-2 pandemic further highlight these characteristics and emphasize the relevance of CoVs to the global public health. Bovine CoVs (BCoVs) are betacoronaviruses associated with neonatal calf diarrhea, and with winter dysentery and shipping fever in older cattle. Of interest, no distinct genetic or antigenic markers have been identified in BCoVs associated with these distinct clinical syndromes. In contrast, like other CoVs, BCoVs exist as quasispecies. Besides cattle, BCoVs and bovine-like CoVs were identified in various domestic and wild ruminant species (water buffalo, sheep, goat, dromedary camel, llama, alpaca, deer, wild cattle, antelopes, giraffes, and wild goats), dogs and humans. Surprisingly, bovine-like CoVs also cannot be reliably distinguished from BCoVs using comparative genomics. Additionally, there are historical examples of zoonotic transmission of BCoVs. This article will discuss BCoV pathogenesis, epidemiology, interspecies transmission, immune responses, vaccines, and diagnostics.

Canine Respiratory Coronavirus, Bovine Coronavirus, and Human Coronavirus OC43: Receptors and Attachment Factors

Despite high similarity of canine respiratory coronavirus (CRCoV), bovine coronavirus, (BCoV) and human coronavirus OC43 (HCoV-OC43), these viruses differ in species specificity. For years it was believed that they share receptor specificity, utilizing sialic acids for cell surface attachment, internalization, and entry. Interestingly, careful literature analysis shows that viruses indeed bind to the cell surface via sialic acids, but there is no solid data that these moieties mediate virus entry. In our study, using a number of techniques, we showed that all three viruses are indeed able to bind to sialic acids to a different extent, but these molecules render the cells permissive only for the clinical strain of HCoV-OC43, while for others they serve only as attachment receptors. CRCoV and BCoV appear to employ human leukocyte antigen class I (HLA-1) as the entry receptor. Furthermore, we identified heparan sulfate as an alternative attachment factor, but this may be related to the cell culture adaptation, as in ex vivo conditions, it does not seem to play a significant role. Summarizing, we delineated early events during CRCoV, BCoV, and HCoV-OC43 entry and systematically studied the attachment and entry receptor utilized by these viruses.

Expression and Function of Zinc-α2-Glycoprotein

Zinc-α2-glycoprotein (ZAG), encoded by the AZGP1 gene, is a major histocompatibility complex I molecule and a lipid-mobilizing factor. ZAG has been demonstrated to promote lipid metabolism and glucose utilization, and to regulate insulin sensitivity. Apart from adipose tissue, skeletal muscle, liver, and kidney, ZAG also occurs in brain tissue, but its distribution in brain is debatable. Only a few studies have investigated ZAG in the brain. It has been found in the brains of patients with Krabbe disease and epilepsy, and in the cerebrospinal fluid of patients with Alzheimer disease, frontotemporal lobe dementia, and amyotrophic lateral sclerosis. Both ZAG protein and AZGP1 mRNA are decreased in epilepsy patients and animal models, while overexpression of ZAG suppresses seizure and epileptic discharges in animal models of epilepsy, but knowledge of the specific mechanism of ZAG in epilepsy is limited. In this review, we summarize the known roles and molecular mechanisms of ZAG in lipid metabolism and glucose metabolism, and in the regulation of insulin sensitivity, and discuss the possible mechanisms by which it suppresses epilepsy.

Sequence, structure and evolutionary analysis of cold shock domain proteins, a member of OB fold family

The cold shock domain (CSD) belongs to the oligosaccharide/oligonucleotide-binding fold superfamily which is highly conserved from prokaryotes to higher eukaryotes, and appears to function as RNA chaperones. CSD is involved in diverse cellular processes, including adaptation to low temperatures, nutrient stress, cellular growth and developmental processes. Structural Classification of Proteins (SCOP) database broadly classifies OB fold proteins into 18 different superfamilies, including nucleic acid-binding superfamily (NAB). The NAB is further divided into 17 families together with cold shock DNA-binding protein family (CSDB). The CSDB have more than 240 000 sequences in UniProt database consisting of 32 domains including CSD. Among these domains, CSD is the second largest sequence contributor (> 40 398 sequences). Herein, we have systematically analysed the relative abundance and distribution of CSD proteins based on sequences, structures, repeats and gene ontology (GO) molecular functions in all domains of life. Analysis of sequence distribution suggesting that CSDs are largely found in bacteria (83-94%) with single CSD repeat. However, repeat distribution in eukaryota varies from 1 to 5 in combination with other auxiliary domain that makes CSD proteins functionally more diverse compared to the bacterial counterparts. Further, analysis of repeats distributions on evolutionary scale suggest that existence of CSD in multiple repeats is mainly driven through speciation, gene shuffling and gene duplication events.

Employing proteomics in the study of antigen presentation: an update

INTRODUCTION

Our immune system discriminates self from non-self by examining the peptide cargo of human leukocyte antigen (HLA) molecules displayed on the cell surface. Successful recognition of HLA-bound non-self peptides can induce T cell responses leading to, for example, the destruction of infected cells. Today, largely due to advances in technology, we have an unprecedented capability to identify the nature of these presented peptides and unravel the true complexity of antigen presentation. Areas covered: In addition to conventional linear peptides, HLA molecules also present post-translationally modified sequences comprising a wealth of chemical and structural modifications, including a novel class of noncontiguous spliced peptides. This review focuses on these emerging themes in antigen presentation and how mass spectrometry in particular has contributed to a new view of the antigenic landscape that is presented to the immune system. Expert Commentary: Advances in the sensitivity of mass spectrometers and use of hybrid fragmentation technologies will provide more information-rich spectra of HLA bound peptides leading to more definitive identification of T cell epitopes. Coupled with improvements in sample preparation and new informatics workflows, studies will access novel classes of peptide antigen and allow interrogation of rare and clinically relevant samples.

Molecular basis of the structural stability of hemochromatosis factor E: A combined molecular dynamic simulation and GdmCl-induced denaturation study

Hemochromatosis factor E (HFE) is a member of class I MHC family and plays a significant role in the iron homeostasis. Denaturation of HFE induced by guanidinium chloride (GdmCl) was measured by monitoring changes in [θ]222 (mean residue ellipticity at 222 nm), intrinsic fluorescence emission intensity at 346 nm (F346 ) and the difference absorption coefficient at 287 nm (Δε287) at pH 8.0 and 25°C. Coincidence of denaturation curves of these optical properties suggests that GdmCl-induced denaturation (native (N) state ↔ denatured (D) state) is a two-state process. The GdmCl-induced denaturation was found reversible in the entire concentration range of the denaturant. All denaturation curves were analyzed for ΔGD0, Gibbs free energy change associated with the denaturation equilibrium (N state ↔ D state) in the absence of GdmCl, which is a measure of HFE stability. We further performed molecular dynamics simulation for 40 ns to see the effect of GdmCl on the structural stability of HFE. A well defined correlation was established between in vitro and in silico studies.

Binding studies and biological evaluation of β-carotene as a potential inhibitor of human calcium/calmodulin-dependent protein kinase IV

Human calcium/calmodulin-dependent protein kinase IV (CAMKIV), a member of Ser/Thr kinase family, is associated with cancer, cerebral hypoxia and neurodegenerative diseases. β-carotene is a colored organic compound, abundant in plants and fruits and is used in cancer prevention. Here, we report a strong binding affinity of β-carotene with CAMKIV using molecular docking, fluorescence binding and isothermal titration calorimetry methods. Furthermore, β-carotene also reduces the enzyme activity of CAMKIV moderately as observed during ATPase assay. To see the role of β-carotene on cell proliferation and apoptosis, cancerous cells (HeLa, HuH7and MCF-7) and normal (HEK-293-T) cell lines were used. Admirable anticancer activity of β-carotene was observed. We further performed propidium iodide and DAPI (4',6-diamidino-2-phenylindole) assays to understand the mechanism of anticancer activity of β-carotene at molecular level. Our findings provide a newer insight into the use of β-carotene in cancer prevention and protection via inhibition of CAMKIV by regulating the signaling pathways.

Chitinase from Thermomyces lanuginosus SSBP and its biotechnological applications

Chitinases are ubiquitous class of extracellular enzymes, which have gained attention in the past few years due to their wide biotechnological applications. The effectiveness of conventional insecticides is increasingly compromised by the occurrence of resistance; thus, chitinase offers a potential alternative to the use of chemical fungicides. The thermostable enzymes from thermophilic microorganisms have numerous industrial, medical, environmental and biotechnological applications due to their high stability for temperature and pH. Thermomyces lanuginosus produced a large number of chitinases, of which chitinase I and II are successfully cloned and purified recently. Molecular dynamic simulations revealed that the stability of these enzymes are maintained even at higher temperature. In this review article we have focused on chitinases from different sources, mainly fungal chitinase of T. lanuginosus and its industrial application.

Effect of pH on the stability of hemochromatosis factor E: a combined spectroscopic and molecular dynamics simulation-based study

Hereditary hemochromatosis is an iron overburden condition, which is mainly governed by hereditary hemochromatosis factor E (HFE), a member of major histocompatibility complex class I. To understand the effect of pH on the structure and stability of HFE, we have cloned, expressed, and purified the HFE in the bacterial system and performed circular dichroism, fluorescence, and absorbance measurements at a wide pH range (pH 3.0-11.0). We found that HFE remains stable in the pH range 7.5-11.0 and gets completely acid denatured at low pH values. In this work, we also analyzed the contribution of salt bridges to the stability of HFE. We further performed molecular dynamics simulations for 80 ns at different pH values. An excellent agreement was observed between results from biophysical and MD simulation studies. At lower pH, HFE undergoes denaturation and may be driven toward a degradation pathway, such as ubiquitination. Hence, HFE is not available to bind again with transferrin receptor1 to negatively regulate iron homeostasis. Further we postulated that, might be low pH of cancerous cells helps them to meet their high iron requirement.

Calcium/calmodulin-dependent protein kinase IV: A multifunctional enzyme and potential therapeutic target

The calcium/calmodulin-dependent protein kinase IV (CAMKIV) belongs to the serine/threonine protein kinase family, and is primarily involved in transcriptional regulation in lymphocytes, neurons and male germ cells. CAMKIV operates the signaling cascade and regulates activity of several transcription activators by phosphorylation, which in turn plays pivotal roles in immune response, inflammation and memory consolidation. In this review, we tried to focus on different aspects of CAMKIV to understand the significance of this protein in the biological system. This enzyme is associated with varieties of disorders such as cerebral hypoxia, azoospermia, endometrial and ovarian cancer, systemic lupus, etc., and hence it is considered as a potential therapeutic target. Structure of CAMKIV is comprised of five distinct domains in which kinase domain is responsible for enzyme activity. CAMKIV is involved in varieties of cellular functions such as regulation of gene expression, T-cell maturation, regulation of survival phase of dendritic cells, bone growth and metabolism, memory consolidation, sperm motility, regulation of microtubule dynamics, cell-cycle progression and apoptosis. In this review, we performed an extensive analysis on structure, function and regulation of CAMKIV and associated diseases.

Classification and Functional Analyses of Putative Conserved Proteins from Chlamydophila pneumoniae CWL029

Chlamydophila pneumoniae, a Gram-negative bacterium belongs to the family Chlamydiaceae, is known to cause community-acquired pneumonia and bronchitis. There is a need for genomic analyses of C. pneumoniae as its chronic infections result in reactive airway disease, lung cancer and asthma. Recent advancement in the sequencing techniques led to the generation of large genomic data. In order to utilize these data, sequence-based function predictions were used for annotating the uncharacterized genes. The genome of C. pneumoniae encodes 1052 proteins, which include a group of 366 functionally uncharacterized proteins, known as "hypothetical proteins" (HPs). Functions of these HPs were predicted by utilizing an integrated approach that combines varieties of bioinformatics tools. The functions of 142 proteins were successfully predicted and categorized into different classes of enzymes, transport proteins, binding proteins and virulence factors. Among these functionally annotated HPs, we were able to identify 12 virulent HPs. Furthermore, the HP with the highest virulence score was subjected to molecular dynamics (MD) simulations to better understand their dynamical behavior in explicit water conditions. These analyses could be utilized for an in-depth understanding of virulence mechanism. The functional knowledge of these proteins could be useful in drug design and discovery process of infections caused by C. pneumoniae.

Structure-based function analysis of putative conserved proteins with isomerase activity from Haemophilus influenzae

Haemophilus influenzae, a Gram-negative bacterium and a member of the family Pasteurellaceae, causes chronic bronchitis, bacteremia, meningitis, etc. The H. influenzae is the first organism whose genome was completely sequenced and annotated. Here, we have extensively analyzed the genome of H. influenzae using available proteins structure and function analysis tools. The objective of this analysis is to assign a precise function to hypothetical proteins (HPs) whose functions are not determined so far. Function prediction of these proteins is helpful in precise understanding of mechanisms of pathogenesis and biochemical pathways important for selecting novel therapeutic target. After an extensive analysis of H. Influenzae genome we have found 13 HPs showing high level of sequence and structural similarity to the enzyme isomerase. Consequently, the structures of HPs have been modeled and analyzed to determine their precise functions. We found these HPs are alanine racemase, lysine 2, 3-aminomutase, topoisomerase DNA-binding C4 zinc finger, pseudouridine synthase B, C and E (Rlu B, C and E), hydroxypyruvate isomerase, nucleoside-diphosphate-sugar epimerase, amidophosphoribosyltransferase, aldose-1-epimerase, tautomerase/MIF, Xylose isomerase-like, have TIM barrel domain and sedoheptulose-7-phosphate isomerase like activity, signifying their corresponding functions in the H. influenzae. This work provides a better understanding of the role HPs with isomerase activities in the survival and pathogenesis of H. influenzae.

Structure-based functional annotation of hypothetical proteins from Candida dubliniensis: a quest for potential drug targets

Candida dubliniensis is an emerging pathogenic yeast in humans and infections are usually restricted to mucosal parts of the body. However, its presence in specimens of immunocompromised individuals, especially in HIV-positive patients, is of major medical concern. There is a large fraction of genomes of C. dubliniensis in the database which are uncharacterized for their biochemical, biophysical, and/or cellular functions, and are identified as hypothetical proteins (HPs). Function annotation of Candida genome is, therefore, essentially required to facilitate the understanding of mechanisms of pathogenesis and biochemical pathways important for selecting novel therapeutic target. Here, we carried out an extensive analysis to explain the functional properties of genome, using available protein structure and function analysis tools. We successfully modeled the structures of eight HPs for which a template with moderate sequence similarity was available in the protein data bank. All modeled structures were analyzed and we found that these proteins may act as transporter, kinase, transferase, ketosteroid, isomerase, hydrolase, oxidoreductase, and binding targets for DNA and RNA. Since these unique HPs of Candida showed no homologs in humans, these proteins are expected to be a potential target for future antifungal therapy.

Towards New Drug Targets? Function Prediction of Putative Proteins of Neisseria meningitidis MC58 and Their Virulence Characterization

Neisseria meningitidis is a Gram-negative aerobic diplococcus, responsible for a variety of meningococcal diseases. The genome of N. meningitidis MC58 is comprised of 2114 genes that are translated into 1953 proteins. The 698 genes (∼35%) encode hypothetical proteins (HPs), because no experimental evidence of their biological functions are available. Analyses of these proteins are important to understand their functions in the metabolic networks and may lead to the discovery of novel drug targets against the infections caused by N. meningitidis. This study aimed at the identification and categorization of each HP present in the genome of N. meningitidis MC58 using computational tools. Functions of 363 proteins were predicted with high accuracy among the annotated set of HPs investigated. The reliably predicted 363 HPs were further grouped into 41 different classes of proteins, based on their possible roles in cellular processes such as metabolism, transport, and replication. Our studies revealed that 22 HPs may be involved in the pathogenesis caused by this microorganism. The top two HPs with highest virulence scores were subjected to molecular dynamics (MD) simulations to better understand their conformational behavior in a water environment. We also compared the MD simulation results with other virulent proteins present in N. meningitidis. This study broadens our understanding of the mechanistic pathways of pathogenesis, drug resistance, tolerance, and adaptability for host immune responses to N. meningitidis.

Dynamics of free versus complexed β2-microglobulin and the evolution of interfaces in MHC class I molecules

In major histocompatibility complex (MHC) class I molecules, monomorphic β(2)-microglobulin (β(2)m) is non-covalently bound to a heavy chain (HC) exhibiting a variable degree of polymorphism. β(2)M can stabilize a wide variety of complexes ranging from classical peptide binding to nonclassical lipid presenting MHC class I molecules as well as to MHC class I-like molecules that do not bind small ligands. Here we aim to assess the dynamics of individual regions in free as well as complexed β(2)m and to understand the evolution of the interfaces between β(2)m and different HC. Using human β(2)m and the HLA-B*27:09 complex as a model system, a comparison of free and HC-bound β(2)m by nuclear magnetic resonance spectroscopy was initially carried out. Although some regions retain their flexibility also after complex formation, these studies reveal that most parts of β(2)m gain rigidity upon binding to the HC. Sequence analyses demonstrate that some of the residues exhibiting flexibility participate in evolutionarily conserved β(2)m-HC contacts which are detectable in diverse vertebrate species or characterize a particular group of MHC class I complexes such as peptide- or lipid-binding molecules. Therefore, the spectroscopic experiments and the interface analyses demonstrate that β(2)m fulfills its role of interacting with diverse MHC class I HC as well as effector cell receptors not only by engaging in conserved intermolecular contacts but also by falling back upon key interface residues that exhibit a high degree of flexibility.

Prediction and identification of T cell epitopes in the H5N1 influenza virus nucleoprotein in chicken

T cell epitopes can be used for the accurate monitoring of avian influenza virus (AIV) immune responses and the rational design of vaccines. No T cell epitopes have been previously identified in the H5N1 AIV virus nucleoprotein (NP) in chickens. For the first time, this study used homology modelling techniques to construct three-dimensional structures of the peptide-binding domains of chicken MHC class Ι molecules for four commonly encountered unique haplotypes, i.e., B4, B12, B15, and B19. H5N1 AIV NP was computationally parsed into octapeptides or nonapeptides according to the peptide-binding motifs of MHC class I molecules of the B4, B12, B15 and B19 haplotypes. Seventy-five peptide sequences were modelled and their MHC class I molecule-binding abilities were analysed by molecular docking. Twenty-five peptides (Ten for B4, six for B12, two for B15, and seven for B19) were predicted to be potential T cell epitopes in chicken. Nine of these peptides and one unrelated peptide were manually synthesized and their T cell responses were tested in vitro. Spleen lymphocytes were collected from SPF chickens that had been immunised with a NP-expression plasmid, pCAGGS-NP, and they were stimulated using the synthesized peptides. The secretion of chicken IFN-γ and the proliferation of CD8(+) T cells were tested using an ELISA kit and flow cytometry, respectively. The significant secretion of chicken IFN-γ and proliferation of CD8(+) T lymphocytes increased by 13.7% and 11.9% were monitored in cells stimulated with peptides NP(89-97) and NP(198-206), respectively. The results indicate that peptides NP(89-97) (PKKTGGPIY) and NP(198-206) (KRGINDRNF) are NP T cell epitopes in chicken of certain haplotypes. The method used in this investigation is applicable to predicting T cell epitopes for other antigens in chicken, while this study also extends our understanding of the mechanisms of the immune response to AIV in chicken.