MHCBI: a pipeline for calculating peptide-MHC binding energy using semi-empirical quantum mechanical methods with explicit/implicit solvent models.

Experimentally estimating peptide-major histocompatibility complex (pMHC) binding affinity has been quite challenging due to the many receptors and the many potential ligands implicated in it. We have thus proposed a straightforward computational methodology considering the different mechanisms involved in pMHC binding to facilitate studying such receptor-ligand interactions. We have developed a pipeline using semi-empirical quantum mechanical methods for calculating pMHC class I and II molecules' binding energy (BE). This pipeline can systematize the methodology for calculating pMHC system BE, enabling the rational design of T-cell epitopes to be used as pharmaceuticals and vaccines.

The molecular basis for peptide-based antimalarial vaccine development targeting erythrocyte invasion by P. falciparum.

This work describes a methodology for developing a minimal, subunit-based, multi-epitope, multi-stage, chemically-synthesised, anti-Plasmodium falciparum malaria vaccine. Some modified high activity binding peptides (mHABPs) derived from functionally relevant P. falciparum MSP, RH5 and AMA-1 conserved amino acid regions (cHABPs) for parasite binding to and invasion of red blood cells (RBC) were selected. They were highly immunogenic as assessed by indirect immunofluorescence (IFA) and Western blot (WB) assays and protective immune response-inducers against malarial challenge in the Aotus monkey experimental model. NetMHCIIpan 4.0 was used for predicting peptide-Aotus/human major histocompatibility class II (MHCII) binding affinity in silico due to the similarity between Aotus and human immune system molecules; ∼50% of Aotus MHCII allele molecules have a counterpart in the human immune system, being Aotus-specific, whilst others enabled recognition of their human counterparts. Some peptides' 1H-NMR-assessed structural conformation was determined to explain residue modifications in mHABPs inducing secondary structure changes. These directly influenced immunological behaviour, thereby highlighting the relationship with MHCII antigen presentation. The data obtained in such functional, immunological, structural and predictive approach suggested that some of these peptides could be excellent components of a fully-protective antimalarial vaccine.

Structural Modelling of KCNQ1 and KCNH2 Double Mutant Proteins, Identified in Two Severe Long QT Syndrome Cases, Reveals New Insights into Cardiac Channelopathies.

Congenital long QT syndrome (LQTS) is a cardiac channelopathy characterized by a prolongation of the QT interval and T-wave abnormalities, caused, in most cases, by mutations in KCNQ1, KCNH2, and SCN5A. Although the predominant pattern of LQTS inheritance is autosomal dominant, compound heterozygous mutations in genes encoding potassium channels have been reported, often with early disease onset and more severe phenotypes. Since the molecular mechanisms underlying severe phenotypes in carriers of compound heterozygous mutations are unknown, it is possible that these compound mutations lead to synergistic or additive alterations to channel structure and function. In this study, all-atom molecular dynamic simulations of KCNQ1 and hERG channels were carried out, including wild-type and channels with compound mutations found in two patients with severe LQTS phenotypes and limited family history of the disease. Because channels can likely incorporate different subunit combinations from different alleles, there are multiple possible configurations of ion channels in LQTS patients. This analysis allowed us to establish the structural impact of different configurations of mutant channels in the activated/open state. Our data suggest that channels with these mutations show moderate changes in folding energy (in most cases of stabilizing character) and changes in channel mobility and volume, differentiating them from each other and from WT. This would indicate possible alterations in K+ ion flow. Hetero-tetrameric mutant channels showed intermediate structural and volume alterations vis-à-vis homo-tetrameric channels. These findings support the hypothesis that hetero-tetrameric channels in patients with compound heterozygous mutations do not necessarily lead to synergistic structural alterations.

Malaria: Paving the way to developing peptide-based vaccines against invasion in infectious diseases.

Malaria remains a large-scale public health problem, killing more than 400,000 people and infecting up to 230 million worldwide, every year. Unfortunately, despite numerous efforts and research concerning vaccine development, results to date have been low and/or strain-specific. This work describes a strategy involving Plasmodium falciparum Duffy binding-like (DBL) and reticulocyte-binding protein homologue (RH) family-derived minimum functional peptides, netMHCIIpan3.2 parental and modified peptides' in silico binding prediction and modeling some Aotus major histocompatibility class II (MHCII) molecules based on known human molecules' structure to understand their differences. These are used to explain peptides' immunological behaviour when used as vaccine components in the Aotus model. Despite the great similarity between human and Aotus immune system molecules, around 50% of Aotus allele molecules lack a counterpart in the human immune system which could lead to an Aotus-specific vaccine. It was also confirmed that functional Plasmodium falciparum' conserved proteins are immunologically silent (in both the animal model and in-silico prediction); they must therefore be modified to elicit an appropriate immune response. Some peptides studied here had the desired behaviour and can thus be considered components of a fully-protective antimalarial vaccine.

Assessing Peptide Binding to MHC II: An Accurate Semiempirical Quantum Mechanics Based Proposal.

Estimating peptide-major histocompatibility complex (pMHC) binding using structural computational methods has an impact on understanding overall immune function triggering adaptive immune responses in MHC class II molecules. We developed a strategy for optimizing pMHC structure interacting with water molecules and for calculating the binding energy of receptor + ligand systems, such as HLA-DR1 + HA, HLA-DR1 + CLIP, HLA-DR2 + MBP, and HLA-DR3 + CLIP, as well as a monosubstitution panel. Taking pMHC's structural properties, we assumed that ΔH ≫ -TΔS would generate a linear model for estimating relative free energy change, using three semiempirical quantum methods (PM6, PM7, and FMO-SCC-DFTB3) along with the implicit solvent models, and considering proteins in neutral and charged states. Likewise, we confirmed our approach's effectiveness in calculating binding energies having high correlation with experimental data and low root-mean-square error (<2 kcal/mol). All in all, our pipeline differentiates weak from strong peptide binders as a reliable method for studying pMHC interactions.

Structural analysis of owl monkey MHC-DR shows that fully-protective malaria vaccine components can be readily used in humans.

More than 50 years ago the owl monkey (genus Aotus) was found to be highly susceptible to developing human malaria, making it an excellent experimental model for this disease. Microbes and parasites' (especially malaria) tremendous genetic variability became resolved during our malaria vaccine development, involving conserved peptides having high host cell binding activity (cHABPs); however, cHABPs are immunologically silent and must be specially modified (mHABPs) to induce a perfect fit into major histocompatibility complex (MHC) molecules (HLA in humans). Since malarial immunity is mainly antibody-mediated and controlled by the HLA-DRB genetic region, ∼1000 Aotus have been molecularly characterised for MHC-DRB, revealing striking similarity between human and Aotus MHC-DRB repertories. Such convergence suggested that a large group of immune protection-inducing protein structures (IMPIPS), highly immunogenic and protection inducers against malarial intravenous challenge in Aotus, could easily be used in humans for inducing full protection against malaria. We highlight the value of a logical and rational methodology for developing a vaccine in an appropriate animal model: Aotus monkeys.

TCR-contacting residues orientation and HLA-DRß* binding preference determine long-lasting protective immunity against malaria.

Fully-protective, long-lasting, immunological (FPLLI) memory against Plasmodium falciparum malaria regarding immune protection-inducing protein structures (IMPIPS) vaccinated into monkeys previously challenged and re-challenged 60 days later with a lethal Aotus monkey-adapted P. falciparum strain was found to be associated with preferential high binding capacity to HLA-DRß1* allelic molecules of the major histocompatibility class II (MHC-II), rather than HLA-DRß3*, ß4*, ß5* alleles. Complete PPIIL 3D structure, a longer distance (26.5 Å ± 1.5 Å) between residues perfectly fitting into HLA-DRß1*PBR pockets 1 and 9, a gauche(-) rotamer orientation in p8 TCR-contacting polar residue and a larger volume of polar p2 residues was also found. This data, in association with previously-described p3 and p7 apolar residues having gauche(+) orientation to form a perfect MHC-II-peptide-TCR complex, determines the stereo-electronic and topochemical characteristics associated with FPLLI immunological memory.

Characterising atypical Candida albicans clinical isolates from six third-level hospitals in Bogotá, Colombia.

BACKGROUND: Candida species are the most frequently found fungal pathogens causing nosocomial disease in a hospital setting. Such species must be correctly identified to ensure that appropriate control measures are taken and that suitable treatment is given for each species. Candida albicans is causing most fungal disease burden worldwide; the challenge lies in differentiating it from emerging atypical, minor and related species such as Candida dubliniensis and Candida africana. The purpose of this study was to compare identification based on MALDI-TOF MS to standard identification systems using a set of nosocomial isolates. METHODS: Eleven nosocomial samples were collected from 6 third-level hospitals in Bogotá, Colombia. All the samples were identified by combining MALDI-TOF MS with morphological characters, carbohydrate assimilation and molecular markers (D1/D2 and HWP1). RESULTS: The present work describes the first collection of atypical Colombian Candida clinical isolates; these were identified as Candida albicans/Candida africana by their MALDI-TOF MS profile. Phenotypical characteristics showed that they were unable to produce chlamydospores, assimilate trehalose, glucosamine, N- acetyl-glucosamine and barely grew at 42 °C, as would be expected for Candida africana. The molecular identification of the D1/D2 region of large subunit ribosomal RNA and HWP1 hyphal cell wall protein 1 sequences from these isolates was consistent with those for Candida albicans. The mass spectra obtained by MALDI-TOF MS were analysed by multi-dimensional scaling (MDS) and cluster analysis, differences being revealed between Candida albicans, Candida africana, Candida dubliniensis reference spectra and two clinical isolate groups which clustered according to the clinical setting, one of them being clearly related to C. albicans. CONCLUSION: This study highlights the importance of using MALDI-TOF MS in combination with morphology, substrate assimilation and molecular markers for characterising Candida albicans-related and atypical C. albicans species, thereby overcoming conventional identification methods. This is the first report of hospital-obtained isolates of this type in Colombia; the approach followed might be useful for gathering knowledge regarding local epidemiology which could, in turn, have an impact on clinical management. The findings highlight the complexity of distinguishing between typical and atypical Candida albicans isolates in hospitals.

Identifying major histocompatibility complex class II-DR molecules in bovine and swine peripheral blood monocyte-derived macrophages using mAb-L243.

Major histocompatibility complex class II (MHC-II) molecules are involved in immune responses against pathogens and vaccine candidates' immunogenicity. Immunopeptidomics for identifying cancer and infection-related antigens and epitopes have benefited from advances in immunopurification methods and mass spectrometry analysis. The mouse anti-MHC-II-DR monoclonal antibody L243 (mAb-L243) has been effective in recognising MHC-II-DR in both human and non-human primates. It has also been shown to cross-react with other animal species, although it has not been tested in livestock. This study used mAb-L243 to identify Staphylococcus aureus and Salmonella enterica serovar Typhimurium peptides binding to cattle and swine macrophage MHC-II-DR molecules using flow cytometry, mass spectrometry and two immunopurification techniques. Antibody cross-reactivity led to identifying expressed MHC-II-DR molecules, together with 10 Staphylococcus aureus peptides in cattle and 13 S. enterica serovar Typhimurium peptides in swine. Such data demonstrates that MHC-II-DR expression and immunocapture approaches using L243 mAb represents a viable strategy for flow cytometry and immunopeptidomics analysis of bovine and swine antigen-presenting cells.

SM-COLSARSPROT: Highly Immunogenic Supramutational Synthetic Peptides Covering the World's Population.

Fifty ~20-amino acid (aa)-long peptides were selected from functionally relevant SARS-CoV-2 S, M, and E proteins for trial B-21 and another 53 common ones, plus some new ones derived from the virus' main genetic variants for complementary trial C-21. Peptide selection was based on tremendous SARS-CoV-2 genetic variability for analysing them concerning vast human immunogenetic polymorphism for developing the first supramutational, Colombian SARS-protection (SM-COLSARSPROT), peptide mixture. Specific physicochemical rules were followed, i.e., aa predilection for polyproline type II left-handed (PPIIL) formation, replacing ß-branched, aromatic aa, short-chain backbone H-bond-forming residues, π-π interactions (n→π* and π-CH), aa interaction with π systems, and molecular fragments able to interact with them, disrupting PPIIL propensity formation. All these modified structures had PPIIL formation propensity to enable target peptide interaction with human leukocyte antigen-DRß1* (HLA-DRß1*) molecules to mediate antigen presentation and induce an appropriate immune response. Such modified peptides were designed for human use; however, they induced high antibody titres against S, M, and E parental mutant peptides and neutralising antibodies when suitably modified and chemically synthesised for immunising 61 major histocompatibility complex class II (MHCII) DNA genotyped Aotus monkeys (matched with their corresponding HLA-DRß1* molecules), predicted to cover 77.5% to 83.1% of the world's population. Such chemically synthesised peptide mixture represents an extremely pure, stable, reliable, and cheap vaccine for COVID-19 pandemic control, providing a new approach for a logical, rational, and soundly established methodology for other vaccine development.

The First Chemically-Synthesised, Highly Immunogenic Anti-SARS-CoV-2 Peptides in DNA Genotyped Aotus Monkeys for Human Use.

Thirty-five peptides selected from functionally-relevant SARS-CoV-2 spike (S), membrane (M), and envelope (E) proteins were suitably modified for immunising MHC class II (MHCII) DNA-genotyped Aotus monkeys and matched with HLA-DRß1* molecules for use in humans. This was aimed at producing the first minimal subunit-based, chemically-synthesised, immunogenic molecules (COLSARSPROT) covering several HLA alleles. They were predicted to cover 48.25% of the world's population for 6 weeks (short-term) and 33.65% for 15 weeks (long-lasting) as they induced very high immunofluorescent antibody (IFA) and ELISA titres against S, M and E parental native peptides, SARS-CoV-2 neutralising antibodies and host cell infection. The same immunological methods that led to identifying new peptides for inclusion in the COLSARSPROT mixture were used for antigenicity studies. Peptides were analysed with serum samples from patients suffering mild or severe SARS-CoV-2 infection, thereby increasing chemically-synthesised peptides' potential coverage for the world populations up to 62.9%. These peptides' 3D structural analysis (by 1H-NMR acquired at 600 to 900 MHz) suggested structural-functional immunological association. This first multi-protein, multi-epitope, minimal subunit-based, chemically-synthesised, highly immunogenic peptide mixture highlights such chemical synthesis methodology's potential for rapidly obtaining very pure, highly reproducible, stable, cheap, easily-modifiable peptides for inducing immune protection against COVID-19, covering a substantial percentage of the human population.

Robust, Comprehensive Molecular, and Phenotypical Characterisation of Atypical Candida albicans Clinical Isolates From Bogotá, Colombia.

Candida albicans is commensal in human microbiota and is known to be the commonest opportunistic pathogen, having variable clinical outcomes that can lead to up to 60% mortality. Such wide clinical behaviour can be attributed to its phenotypical plasticity and high genetic diversity. This study characterised 10 Colombian clinical isolates which had already been identified as C. albicans by molecular tests; however, previous bioinformatics analysis of protein mass spectra and phenotypical characteristics has shown that this group of isolates has atypical behaviour, sharing characteristics of both C. africana and C. albicans. This study was aimed at evaluating atypical isolates' pathogenic capability in the Galleria mellonella model; susceptibility profiles were determined and MLST was used for molecular characterisation. Cluster analysis, enabling unbiased bootstrap to classify the isolates and establish their cluster membership and e-BURST, was used for establishing clonal complexes (CC). Both approaches involved using representative MLST data from the 18 traditional C. albicans clades, as well as C. albicans-associated and minor species. Ten atypical isolates were distributed as follows: 6/10 (B71, B41, B60, R6, R41, and R282) were grouped into a statistically well-supported atypical cluster (AC) and constituted a differentiated CC 6; 2/10 of the isolates were clearly grouped in clade 1 and were concurrent in CC 4 (B80, B44). Another 2/10 atypical isolates were grouped in clade 10 and concurred in CC 7 (R425, R111); most atypical isolates were related to geographically distant isolates and some represented new ST. Isolates B41 and R41 in the AC had greater virulence. Isolate B44 was fluconazole-resistant and was grouped in clade 1. The atypical nature of the isolates studied here was demonstrated by the contrast between phenotypical traits (C. africana-like), molecular markers (C. albicans-like), virulence, and antifungal resistance, highlighting the widely described genetic plasticity for this genus. Our results showed that the atypical isolates forming well-differentiated groups belonged to C. albicans. Our findings could contribute towards developing molecular epidemiology approaches for managing hospital-acquired infection.

Plasmodium vivax Cell Traversal Protein for Ookinetes and Sporozoites (CelTOS) Functionally Restricted Regions Are Involved in Specific Host-Pathogen Interactions.

Following the injection of Plasmodium sporozoites by a female Anopheles mosquito into the dermis, they become engaged on a long journey to hepatic tissue where they must migrate through different types of cell to become established in parasitophorous vacuoles in hepatocytes. Studies have shown that proteins such as cell traversal protein for Plasmodium ookinetes and sporozoites (CelTOS) play a crucial role in cell-traversal ability. Although CelTOS has been extensively studied in various species and included in pre-clinical assays it remains unknown which P. vivax CelTOS (PvCelTOS) regions are key in its interaction with traversed or target cells (Kupffer or hepatocytes) and what type of pressure, association and polymorphism these important regions could have to improve their candidacy as important vaccine antigens. This work has described producing a recombinant PvCelTOS which was recognized by ~30% P. vivax-infected individuals, thereby confirming its ability for inducing a natural immune response. PvCelTOS' genetic diversity in Colombia and its ability to interact with HeLa (traversal cell) and/or HepG2 cell (target cell) external membrane have been assessed. One region in the PvCelTOS amino-terminal region and another in its C-terminus were seen to be participating in host-pathogen interactions. These regions had important functional constraint signals (ω < 0.3 and several sites under negative selection) and were able to inhibit specific rPvCelTOS binding to HeLa cells. This led to suggesting that sequences between aa 41-60 (40833) and 141-160 (40838) represent promising candidates for an anti-P. vivax subunit-based vaccine.

Publisher Correction: Mass & secondary structure propensity of amino acids explain their mutability and evolutionary replacements.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

Mass & secondary structure propensity of amino acids explain their mutability and evolutionary replacements.

Why is an amino acid replacement in a protein accepted during evolution? The answer given by bioinformatics relies on the frequency of change of each amino acid by another one and the propensity of each to remain unchanged. We propose that these replacement rules are recoverable from the secondary structural trends of amino acids. A distance measure between high-resolution Ramachandran distributions reveals that structurally similar residues coincide with those found in substitution matrices such as BLOSUM: Asn ↔ Asp, P™he ↔ Tyr, Lys ↔ Arg, Gln ↔ Glu, Ile ↔ Val, Met → Leu; with Ala, Cys, His, Gly, Ser, Pro, and Thr, as structurally idiosyncratic residues. We also found a high average correlation ([Formula: see text] = 0.85) between thirty amino acid mutability scales and the mutational inertia (I X ), which measures the energetic cost weighted by the number of observations at the most probable amino acid conformation. These results indicate that amino acid substitutions follow two optimally-efficient principles: (a) amino acids interchangeability privileges their secondary structural similarity, and (b) the amino acid mutability depends directly on its biosynthetic energy cost, and inversely with its frequency. These two principles are the underlying rules governing the observed amino acid substitutions.

How frequently do clusters occur in hierarchical clustering analysis? A graph theoretical approach to studying ties in proximity.

BACKGROUND: Hierarchical cluster analysis (HCA) is a widely used classificatory technique in many areas of scientific knowledge. Applications usually yield a dendrogram from an HCA run over a given data set, using a grouping algorithm and a similarity measure. However, even when such parameters are fixed, ties in proximity (i.e. two equidistant clusters from a third one) may produce several different dendrograms, having different possible clustering patterns (different classifications). This situation is usually disregarded and conclusions are based on a single result, leading to questions concerning the permanence of clusters in all the resulting dendrograms; this happens, for example, when using HCA for grouping molecular descriptors to select that less similar ones in QSAR studies. RESULTS: Representing dendrograms in graph theoretical terms allowed us to introduce four measures of cluster frequency in a canonical way, and use them to calculate cluster frequencies over the set of all possible dendrograms, taking all ties in proximity into account. A toy example of well separated clusters was used, as well as a set of 1666 molecular descriptors calculated for a group of molecules having hepatotoxic activity to show how our functions may be used for studying the effect of ties in HCA analysis. Such functions were not restricted to the tie case; the possibility of using them to derive cluster stability measurements on arbitrary sets of dendrograms having the same leaves is discussed, e.g. dendrograms from variations of HCA parameters. It was found that ties occurred frequently, some yielding tens of thousands of dendrograms, even for small data sets. CONCLUSIONS: Our approach was able to detect trends in clustering patterns by offering a simple way of measuring their frequency, which is often very low. This would imply, that inferences and models based on descriptor classifications (e.g. QSAR) are likely to be biased, thereby requiring an assessment of their reliability. Moreover, any classification of molecular descriptors is likely to be far from unique. Our results highlight the need for evaluating the effect of ties on clustering patterns before classification results can be used accurately.Graphical abstractFour cluster contrast functions identifying statistically sound clusters within dendrograms considering ties in proximity.

High level of conservation in Plasmodium vivax merozoite surface protein 4 (PvMSP4).

Plasmodium vivax merozoite surface protein 4 (PvMSP4) has been considered to be a promising malarial vaccine candidate. The antigenic diversity displayed by parasite populations is one of the main factors limiting the efficacy of asexual-stage anti-malarial vaccine candidates. The present study is the first characterising PvMSP4 polymorphism. P. vivax isolates were collected from endemic areas in Colombia and diversity and selection pattern studies were carried out. Overall conservation in this protein was remarkable. Changes were only found in exons I and II, the former only having single nucleotide polymorphisms (SNPs) whilst the latter showed variations in tandem repeat number caused by exon II slippage. The remaining regions (EGF-like domain, GPI anchor and intron) were completely conserved. Selection and neutrality tests carried out over variant exons indicated negative selective forces acting on them. No evidence of intragenic recombination was found. The strong conservation displayed in this molecule by isolates from geographically different regions (Colombia, Salvador and Thailand) stresses its potential importance as a candidate for a vaccine against P. vivax malaria.

Characterising a microsatellite for DRB typing in Aotus vociferans and Aotus nancymaae (Platyrrhini).

Non-human primates belonging to the Aotus genus have been shown to be excellent experimental models for evaluating drugs and vaccine candidates against malaria and other human diseases. The immune system of this animal model must be characterised to assess whether the results obtained here can be extrapolated to humans. Class I and II major histocompatibility complex (MHC) proteins are amongst the most important molecules involved in response to pathogens; in spite of this, the techniques available for genotyping these molecules are usually expensive and/or time-consuming. Previous studies have reported MHC-DRB class II gene typing by microsatellite in Old World primates and humans, showing that such technique provides a fast, reliable and effective alternative to the commonly used ones. Based on this information, a microsatellite present in MHC-DRB intron 2 and its evolutionary patterns were identified in two Aotus species (A. vociferans and A. nancymaae), as well as its potential for genotyping class II MHC-DRB in these primates.

Owl monkey MHC-DRB exon 2 reveals high similarity with several HLA-DRB lineages.

One hundred and ten novel MHC-DRB gene exon 2 nucleotide sequences were sequenced in 96 monkeys from three owl monkey species (67 from Aotus nancymaae, 30 from Aotus nigriceps and 13 from Aotus vociferans). Owl monkeys, like humans, have high MHC-DRB allele polymorphism, revealing a striking similarity with several human allele lineages in the peptide binding region and presenting major convergence with DRB lineages from several Catarrhini (humans, apes and Old World monkeys) rather than with others New World monkeys (Platyrrhini). The parallelism between human and Aotus MHC-DRB reveals additional similarities regarding variability pattern, selection pressure and physicochemical constraints in amino acid replacements. These observations concerning previous findings of similarity between the Aotus immune system molecules and their human counterparts affirm this specie's usefulness as an excellent animal model in biomedical research.

MHC class I genes in the owl monkey: mosaic organisation, convergence and loci diversity.

The MHC class I molecule plays an important role in immune response, pathogen recognition and response against vaccines and self- versus non-self-recognition. Studying MHC class I characteristics thus became a priority when dealing with Aotus to ensure its use as an animal model for biomedical research. Isolation, cloning and sequencing of exons 1-8 from 27 MHC class I alleles obtained from 13 individuals classified as belonging to three owl monkey species (A. nancymaae, A. nigriceps and A. vociferans) were carried out to establish similarities between Aotus MHC class I genes and those expressed by other New and Old World primates. Six Aotus MHC class I sequence groups (Ao-g1, Ao-g2, Ao-g3, Ao-g4, Ao-g5 and Ao-g6) weakly related to non-classical Catarrhini MHC were identified. An allelic lineage was also identified in one A. nancymaae and two A. vociferans monkeys, exhibiting a high degree of conservation, negative selection along the molecule and premature termination of the open reading frame at exon 5 (Ao-g5). These sequences' high conservation suggests that they more likely correspond to a soluble form of Aotus MHC class I molecules than to a new group of processed pseudogenes. Another group, named Ao-g6, exhibited a strong relationship with Catarrhini's classical MHC-B-C loci. Sequence evolution and variability analysis indicated that Aotus MHC class I molecules experience inter-locus gene conversion phenomena, contributing towards their high variability.