NetMHCpan,a method for MHC class I binding prediction beyond humans

Binding of peptides to major histocompatibility complex (MHC) molecules is the single most selective step in the recognition of pathogens by the cellular immune system. The human MHC genomic region (called HLA) is extremely polymorphic comprising several thousand alleles, each encoding a distinct MHC molecule. The potentially unique specificity of the majority of HLA alleles that have been identified to date remains uncharacterized. Likewise, only a limited number of chimpanzee and rhesus macaque MHC class I molecules have been characterized experimentally. Here, we present NetMHCpan-2.0, a method that generates quantitative predictions of the affinity of any peptide–MHC class I interaction. NetMHCpan-2.0 has been trained on the hitherto largest set of quantitative MHC binding data available, covering HLA-A and HLA-B, as well as chimpanzee, rhesus macaque, gorilla, and mouse MHC class I molecules. We show that the NetMHCpan-2.0 method can accurately predict binding to uncharacterized HLA molecules, including HLA-C and HLA-G. Moreover, NetMHCpan-2.0 is demonstrated to accurately predict peptide binding to chimpanzee and macaque MHC class I molecules. The power of NetMHCpan-2.0 to guide immunologists in interpreting cellular immune responses in large out-bred populations is demonstrated. Further, we used NetMHCpan-2.0 to predict potential binding peptides for the pig MHC class I molecule SLA-1*0401. Ninety-three percent of the predicted peptides were demonstrated to bind stronger than 500 nM. The high performance of NetMHCpan-2.0 for non-human primates documents the method’s ability to provide broad allelic coverage also beyond human MHC molecules. The method is available at . Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Diversification of porcine MHC class II genes: evidence for selective advantage

The major histocompatibility complex (MHC) is an immunological gene-dense region of high diversity in mammalian species. Sus scrofa was domesticated by at least six independent events over Eurasia during the Holocene period. It has been hypothesized that the level and distribution of MHC variation in pig populations reflect genetic selection and environmental influences. In an effort to define the complexity of MHC polymorphisms and the role of selection in the generation of class II gene diversity (DQB, DRB1, and pseudogene ΨDRB3), DNA from globally distributed unrelated domestic pigs of European and Asian origins and a Suidae out-group was analyzed. The number of pseudogene alleles identified (ΨDRB3 33) was greater than those found in the expressed genes (DQB 20 and DRB1 23) but the level of observed heterozygosity (ΨDRB3 0.452, DQB 0.732, and DRB1 0.767) and sequence diversity (ΨDRB3 0.029, DQB 0.062, and DRB1 0.074) were significantly lower in the pseudogene, respectively. The substitution ratios reflected an excess of d _N (DQB 1.476, DRB1 1.724, and ΨDRB3 0.508) and the persistence of expressed gene alleles suggesting the influence of balancing selection, while the pseudogene was undergoing purifying selection. The lack of a clear MHC phylogeographic tree, coupled with close genetic distances observed between the European and Asian populations (DQB 0.047 and DRB1 0.063) suggested that unlike observations using mtDNA, the MHC diversity lacks phylogeographic structure and appears to be globally uniform. Taken together, these results suggest that, despite regional differences in selective breeding and environments, no skewing of MHC diversity has occurred. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

MHC class II polymorphism is associated with a canine SLE-related disease complex

Nova Scotia duck tolling retrievers are predisposed to a SLE-related disease complex including immune-mediated rheumatic disease (IMRD) and steroid-responsive meningitis–arteritis (SRMA). IMRD involves symptoms that resemble those seen in systemic autoimmune rheumatic diseases, such as systemic lupus erythematosus, SLE, or SLE-related diseases, in humans. This disease complex involves persistent lameness, stiffness, mainly after resting, and palpable pain from several joints of extremities. The majority of affected dogs display antinuclear autoantibody (ANA)-reactivity. SRMA is manifested in young dogs with high fever and neck stiffness and can be treated with corticosteroids. We have investigated the possible role of MHC class II as a genetic risk factor in IMRD and SRMA etiology. We performed sequence-based typing of the DLA-DRB1, -DQA1, and -DQB1 class II loci in a total of 176 dogs including 51 IMRD (33 ANA-positive), 49 SRMA cases, and 78 healthy controls (two dogs were both IMRD- and SRMA-affected). Homozygosity for the risk haplotype DRB1*00601/DQA1*005011/DQB1*02001 increased the risk for IMRD (OR?=?4.9; ANA-positive IMRD: OR?=?7.2) compared with all other genotypes. There was a general heterozygote advantage, homozygotes had OR?=?4.4 (ANA-positive IMRD: OR?=?8.9) compared with all heterozygotes. The risk haplotype contains the five amino acid epitope RARAA, known as the shared epitope for rheumatoid arthritis. No association was observed for SRMA. We conclude that DLA class II is a highly significant genetic risk factor for ANA-positive IMRD. The results indicate narrow diversity of DLA II haplotypes and identify an IMRD-related risk haplotype, which becomes highly significant in homozygous dogs.     

Feline herpesvirus‐1 down‐regulates MHC class I expression in an homologous cell system

Cytotoxic T lymphocytes (CTLs) are an essential component of the immune defense against many virus infections. CTLs recognize viral peptides in the context of the major histocompatibility complex (MHC) class I molecules on the surface of infected cells. Many viruses have evolved mechanisms to interfere with MHC class I expression as a means of evading the host immune response. In the present research we have studied the effect of in vitro Feline Herpesvirus 1 (FeHV‐1) infection on MHC class I expression. The results of this study demonstrate that FeHV‐1 down regulates surface expression of MHC class I molecules on infected cells, presumably to evade cytotoxic T‐cell recognition and, perhaps, attenuate induction of immunity. Sensitivity to UV irradiation and insensitivity to a viral DNA synthesis inhibitor, like phosphonacetic acid, revealed that immediate early or early viral gene(s) are responsible. Use of the protein translation inhibitor cycloheximide confirmed that an early gene is primarily responsible. J. Cell. Biochem. 106: 179–185, 2009. © 2008 Wiley‐Liss, Inc.     

Parental relatedness and major histocompatibility effects on early embryo survivorship in Atlantic salmon

Salmon have provided key insights into the relative influence of natural and sexual selection on major histocompatibility complex (MHC) variation. Natural selection on salmon MHC genes has been demonstrated in pathogen studies, and there is evidence of MHC-based mate choice (sexual selection). We tested whether parental MHC genes affect survivorship of juvenile Atlantic salmon (Salmo salar) by quantifying the influence of parental genome-wide relatedness and MHC genotype on survivorship to the swim-up stage. Thirteen microsatellite loci were used to estimate the influence of genome-wide relatedness between parents on offspring survivorship and MHC genotypes were determined by sequencing part of the class IIβ gene. Our results revealed no significant relationship between early offspring survivorship and genome-wide relatedness, predicted MHC heterozygosity, or MHC allelic similarity. Overall, our data are consistent with the contention that excess MHC heterozygosity in Atlantic salmon juveniles is due to sexual selection as well as differential survival of offspring due to MHC genotype.     

The Role of the N-Terminus of the Myosin Essential Light Chain in Cardiac Muscle Contraction

To study the regulation of cardiac muscle contraction by the myosin essential light chain (ELC) and the physiological significance of its N-terminal extension, we generated transgenic (Tg) mice by partially replacing the endogenous mouse ventricular ELC with either the human ventricular ELC wild type (Tg-WT) or its 43-amino-acid N-terminal truncation mutant (Tg-Δ43) in the murine hearts. The mutant protein is similar in sequence to the short ELC variant present in skeletal muscle, and the ELC protein distribution in Tg-Δ43 ventricles resembles that of fast skeletal muscle. Cardiac muscle preparations from Tg-Δ43 mice demonstrate reduced force per cross-sectional area of muscle, which is likely caused by a reduced number of force-generating myosin cross-bridges and/or by decreased force per cross-bridge. As the mice grow older, the contractile force per cross-sectional area further decreases in Tg-Δ43 mice and the mutant hearts develop a phenotype of nonpathologic hypertrophy while still maintaining normal cardiac performance. The myocardium of older Tg-Δ43 mice also exhibits reduced myosin content. Our results suggest that the role of the N-terminal ELC extension is to maintain the integrity of myosin and to modulate force generation by decreasing myosin neck region compliance and promoting strong cross-bridge formation and/or by enhancing myosin attachment to actin.     

Alternative Exon 9-Encoded Relay Domains Affect More than One Communication Pathway in the Drosophila Myosin Head

We investigated the biochemical and biophysical properties of one of the four alternative regions within the Drosophila myosin catalytic domain: the relay domain encoded by exon 9. This domain of the myosin head transmits conformational changes in the nucleotide-binding pocket to the converter domain, which is crucial to coupling catalytic activity with mechanical movement of the lever arm. To study the function of this region, we used chimeric myosins (IFI-9b and EMB-9a), which were generated by exchange of the exon 9-encoded domains between the native embryonic body wall (EMB) and indirect flight muscle isoforms (IFI). Kinetic measurements show that exchange of the exon 9-encoded region alters the kinetic properties of the myosin S1 head. This is reflected in reduced values for ATP-induced actomyosin dissociation rate constant (K₁k₊₂) and ADP affinity (K_AD), measured for the chimeric constructs IFI-9b and EMB-9a, compared to wild-type IFI and EMB values. Homology models indicate that, in addition to affecting the communication pathway between the nucleotide-binding pocket and the converter domain, exchange of the relay domains between IFI and EMB affects the communication pathway between the nucleotide-binding pocket and the actin-binding site in the lower 50-kDa domain (loop 2). These results suggest an important role of the relay domain in the regulation of actomyosin cross-bridge kinetics.     

Coupling Molecular Dynamics Simulations with Experiments for the Rational Design of Indolicidin-Analogous Antimicrobial Peptides

Antimicrobial peptides (AMPs) have attracted much interest in recent years because of their potential use as new-generation antibiotics. Indolicidin (IL) is a 13-residue cationic AMP that is effective against a broad spectrum of bacteria, fungi, and even viruses. Unfortunately, its high hemolytic activity retards its clinical applications. In this study, we adopted molecular dynamics (MD) simulations as an aid toward the rational design of IL analogues exhibiting high antimicrobial activity but low hemolysis. We employed long-timescale, multi-trajectory all-atom MD simulations to investigate the interactions of the peptide IL with model membranes. The lipid bilayer formed by the zwitterionic 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) was chosen as the model erythrocyte membrane; lipid bilayers formed from a mixture of POPC and the negatively charged 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol were chosen to model bacterial membranes. MD simulations with a total simulation time of up to 4 μs revealed the mechanisms of the processes of IL adsorption onto and insertion into the membranes. The packing order of these lipid bilayers presumably correlated to the membrane stability upon IL adsorption and insertion. We used the degree of local membrane thinning and the reduction in the order parameter of the acyl chains of the lipids to characterize the membrane stability. The order of the mixed 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol/POPC lipid bilayer reduced significantly upon the adsorption of IL. On the other hand, although the order of the pure-POPC lipid bilayer was perturbed slightly during the adsorption stage, the value was reduced more dramatically upon the insertion of IL into the membrane's hydrophobic region. The results imply that enhancing IL adsorption on the microbial membrane may amplify its antimicrobial activity, while the degree of hemolysis may be reduced through inhibition of IL insertion into the hydrophobic region of the erythrocyte membrane. In addition, through simulations, we identified the amino acids that are most responsible for the adsorption onto or insertion into the two model membranes. Positive charges are critical to the peptide's adsorption, whereas the presence of hydrophobic Trp8 and Trp9 leads to its deeper insertion. Combining the hypothetical relationships between the membrane disordering and the antimicrobial and hemolytical activities with the simulated results, we designed three new IL-analogous peptides: IL-K7 (Pro7 → Lys), IL-F89 (Trp8 and Trp9 → Phe), and IL-K7F89 (Pro7 → Lys; Trp8 and Trp9 → Phe). The hemolytic activity of IL-F89 is considerably lower than that of IL, whereas the antimicrobial activity of IL-K7 is greatly enhanced. In particular, the de novo peptide IL-K7F89 exhibits higher antimicrobial activity against Escherichia coli; its hemolytic activity decreased to only 10% of that of IL. Our simulated and experimental results correlated well. This approach—coupling MD simulations with experimental design—is a useful strategy toward the rational design of AMPs for potential therapeutic use.     

A flexible docking approach for prediction of T cell receptor–peptide–MHC complexes

T cell receptors (TCRs) are immune proteins that specifically bind to antigenic molecules, which are often foreign peptides presented by major histocompatibility complex proteins (pMHCs), playing a key role in the cellular immune response. To advance our understanding and modeling of this dynamic immunological event, we assembled a protein–protein docking benchmark consisting of 20 structures of crystallized TCR/pMHC complexes for which unbound structures exist for both TCR and pMHC. We used our benchmark to compare predictive performance using several flexible and rigid backbone TCR/pMHC docking protocols. Our flexible TCR docking algorithm, TCRFlexDock, improved predictive success over the fixed backbone protocol, leading to near‐native predictions for 80% of the TCR/pMHC cases among the top 10 models, and 100% of the cases in the top 30 models. We then applied TCRFlexDock to predict the two distinct docking modes recently described for a single TCR bound to two different antigens, and tested several protein modeling scoring functions for prediction of TCR/pMHC binding affinities. This algorithm and benchmark should enable future efforts to predict, and design of uncharacterized TCR/pMHC complexes.