Isolation and characterization of high affinity and highly stable anti-Chikungunya virus antibodies using ALTHEA Gold Libraries™.

BACKGROUND: More than 3 million infections were attributed to Chikungunya virus (CHIKV) in the 2014-2016 outbreak in Mexico, Central and South America, with over 500 deaths directly or indirectly related to this viral disease. CHIKV outbreaks are recurrent and no vaccine nor approved therapeutics exist to prevent or treat CHIKV infection. Reliable and robust diagnostic methods are thus critical to control future CHIKV outbreaks. Direct CHIKV detection in serum samples via highly specific and high affinity anti-CHIKV antibodies has shown to be an early and effective clinical diagnosis. METHODS: To isolate highly specific and high affinity anti-CHIKV, Chikungunya virions were isolated from serum of a patient in Veracruz, México. After purification and characterization via electron microscopy, SDS-PAGE and binding to well-characterized anti-CHIKV antibodies, UV-inactivated particles were utilized as selector in a solid-phase panning in combination with ALTHEA Gold Libraries™, as source of antibodies. The screening was based on ELISA and Next-Generation Sequencing. RESULTS: The CHIKV isolate showed the typical morphology of the virus. Protein bands in the SDS-PAGE were consistent with the size of CHIKV capsid proteins. UV-inactivated CHIKV particles bound tightly the control antibodies. The lead antibodies here obtained, on the other hand, showed high expression yield, > 95% monomeric content after a single-step Protein A purification, and importantly, had a thermal stability above 75 °C. Most of the antibodies recognized linear epitopes on E2, including the highest affinity antibody called C7. A sandwich ELISA implemented with C7 and a potent neutralizing antibody isolated elsewhere, also specific for E2 but recognizing a discontinuous epitope, showed a dynamic range of 0.2-40.0 mg/mL of UV-inactivated CHIKV purified preparation. The number of CHIKV particles estimated based on the concentration of E2 in the extract suggested that the assay could detect clinically meaningful amounts of CHIKV in serum. CONCLUSIONS: The newly discovered antibodies offer valuable tools for characterization of CHIKV isolates. Therefore, the strategy here followed using whole viral particles and ALTHEA Gold Libraries™ could expedite the discovery and development of antibodies for detection and control of emergent and quickly spreading viral outbreaks.

VIR.II: a new interface with the antibody sequences in the Kabat database.

The Kabat database is the source of information par excellence on antibody sequences. In 1995, we developed an interface with the Kabat database, called VIR. VIR has been very useful in conducting studies aiming to find structure-function relationships in antibodies. Here we report a new version adapted to the World Wide Web, called VIR.II. VIR.II allows searches by type of chain (V(H) or V(L)), by species, and by specificity. The species are selected using a pulldown menu, whereas the specificities can be selected from a list containing the unique specificities reported in the Kabat database. These facilities avoid mistakes and redundancies in the searches. Another feature, and probably the most important one, is that VIR.II introduces a classification of specificities in terms of the chemical and biochemical nature of the antigen, like anti-protein, anti-peptide, anti-hapten, etc. This classification has been useful in discovering patterns in the antigen-binding site of antibodies that correlate with the type of antigen the antibody interacts with. To illustrate this, while showing the capabilities of VIR.II, we analyze all the murine anti-peptide and anti-protein antibody sequences compiled as of July, 2000 in the Kabat database.

Analysis of antibodies of known structure suggests a lack of correspondence between the residues in contact with the antigen and those modified by somatic hypermutation.

Forty unique murine antibody-antigen complexes determined at 2.5 A or less resolution are analyzed to determine whether the residues in direct contact with the antigen are modified by somatic hypermutation. This was done by taking advantage of the recent characterization of the pool of Vkappa germline genes of the mouse. The average number of residues in contact with the antigen in the V(L) gene, which contains the CDRL-1, CDRL-2, and all but one residue of CDRL-3, was six. The average number of somatic mutations was similar (around five). However, as many as 53% of the antibodies did not show somatic replacements of residues in contact with the antigen. Another 28% had only one. Overall, the frequency of antibodies with increasing number of somatic replacements in residues in contact with the antigen decreased exponentially. A possible explanation of this finding is that mutations in the contacting residues have an adverse effect on the antigen-antibody interaction. This implies that most of the observed mutations are those remaining after negative (purifying) selection. Therefore, efficient strategies of site-directed mutagenesis to improve the affinity of antibodies should be focused on residues other than those directly interacting with the antigen.

Phage displayed biomolecules as preventive and therapeutic agents.

Phage display is a powerful technology for selecting and engineering peptides and proteins expressed on the surface of filamentous bacteriophage. The advantages of phage display technology over other research tools and its' great potential have been demonstrated by successful application of phage display in diverse fields of biomedical/clinical research. In this review we will describe some recent developments in phage display, including new expression vectors, display formats, bioselection strategies and applications in pharmaceutical biotechnology. We highlight some important applications of phage display to identify disease- and pathogen-specific biomolecules, making particular emphasis on development of phage display-derived preventive and therapeutic vaccines.

Discrimination between native and intentionally misfolded conformations of proteins: ES/IS, a new method for calculating conformational free energy that uses both dynamics simulations with an explicit solvent and an implicit solvent continuum model.

A new method for calculating the total conformational free energy of proteins in water solvent is presented. The method consists of a relatively brief simulation by molecular dynamics with explicit solvent (ES) molecules to produce a set of microstates of the macroscopic conformation. Conformational energy and entropy are obtained from the simulation, the latter in the quasi-harmonic approximation by analysis of the covariance matrix. The implicit solvent (IS) dielectric continuum model is used to calculate the average solvation free energy as the sum of the free energies of creating the solute-size hydrophobic cavity, of the van der Waals solute-solvent interactions, and of the polarization of water solvent by the solute's charges. The reliability of the solvation free energy depends on a number of factors: the details of arrangement of the protein's charges, especially those near the surface; the definition of the molecular surface; and the method chosen for solving the Poisson equation. Molecular dynamics simulation in explicit solvent relaxes the protein's conformation and allows polar surface groups to assume conformations compatible with interaction with solvent, while averaging of internal energy and solvation free energy tend to enhance the precision. Two recently developed methods--SIMS, for calculation of a smooth invariant molecular surface, and FAMBE, for solution of the Poisson equation via a fast adaptive multigrid boundary element--have been employed. The SIMS and FAMBE programs scale linearly with the number of atoms. SIMS is superior to Connolly's MS (molecular surface) program: it is faster, more accurate, and more stable, and it smooths singularities of the molecular surface. Solvation free energies calculated with these two programs do not depend on molecular position or orientation and are stable along a molecular dynamics trajectory. We have applied this method to calculate the conformational free energy of native and intentionally misfolded globular conformations of proteins (the EMBL set of deliberately misfolded proteins) and have obtained good discrimination in favor of the native conformations in all instances.

Identification of autoimmune thrombocytopenic purpura-related epitopes using phage-display peptide library.

A random heptapeptide phage-displayed library was screened with two serum samples from autoimmune thrombocytopenic purpura (AITP) patients to address the repertoire of autoantigenic epitopes involved in platelet destruction. We obtained a panel of affinity-selected phage clones that have been shown to react in enzyme-linked immunosorbent assay with autoantibodies from other AITP patients. None of the peptides obtained has been described previously as possibly being an epitope for antiplatelet antibodies, and the majority of them did not show any homology with known platelet glycoproteins. We conclude that peptides identified in this study could represent discontinuous epitopes or mimotopes of natural autoantigens. Also, they could be present in still-unknown proteins involved in AITP pathogenesis.

Structural differences between the repertoires of mouse and human germline genes and their evolutionary implications.

Although human and mouse antibodies are similar when one considers their diversification strategies, they differ in the extent to which kappa and lambda light chains are present in their respective variable light chain repertoires. While the Igk-V germline genes are preponderant in mice (95% or more), they comprise only 60% in humans. This may account for differences in the structural repertoire encoded in the Igk-V germline genes of these species. However, this subject has not been properly investigated, partially because a systematic structural characterization of the mouse Igk-V germline genes has not been undertaken. In the present study we compiled all available information on mouse Igk-V germline genes to characterize their structural repertoire. As expected, comparison with the structural repertoire of human Igk-V germline genes indicates differences. The most interesting is that the mouse Igk-V germline gene repertoire is more diverse in structural terms than its human counterpart: the mouse encodes seven canonical structure classes (combination of canonical structures in L1 and L3). In contrast, the human encodes only four. Analysis of the evolutionary relationships of human and mouse Igk-V germline genes led us to propose that the difference reflects a strategy of mice to compensate for the small lambda chain contribution to the repertoire of their variable light chains.

Evolution of the structural repertoire of the human V(H) and Vkappa germline genes.

Variable genes of human Ig are classified in families and clans which reflect the early events of gene duplication in the evolution of the locus. This organization in multiple copies of variable genes plus the somatic processes of recombination and hypermutation allows the immune system to generate an antibody repertoire of great diversity. At present the role that somatic processes play in the generation of that diversity is understood with some detail. It is a matter of hard controversy, however, which selective pressures have shaped the evolution of the germline genes of Ig and, consequently, what the role of this germline component in the generation of the antibody diversity actually is. Previous studies of our group have showed that the structural repertoire of Ig-determined by the canonical structures-is an important factor to determine the recognition properties of the antibodies. Complete knowledge of the sequences of the human V(H) and Vkappa loci is available to analyze the evolution of the structural repertoire of these loci. Two phylogenetic gene trees were built from the functional germline genes and the evolution of the structural repertoire was studied. We report that for both loci the canonical structures are not randomly distributed within the tree. Conversely, it is shown that the evolution of the structural repertoire follows a gradual process of diversification. This indicates a correlation between the evolution of genes and the structural repertoire, although important differences are found in the patterns of evolution of the structural repertoire between V(H) and Vkappa. Based on those results we propose a primordial structural repertoire for V(H) and Vkappa. The general properties and an outline of the three-dimensional structure of this primordial repertoire are given.

The differences between the structural repertoires of VH germ-line gene segments of mice and humans: implication for the molecular mechanism of the immune response.

Although human and murine antibodies are similar when considering their diversification strategies, they differ in the proportion by which kappa and lambda type chains are present in their receptive V, repertoires. It has been shown that this difference implies a divergence in the structural repertoire of the kappa and lambda genes of these species. Nonetheless, the differences in VH have not been systematically studied. In this paper a systematic characterization of the VH structural repertoire of mice is made, so that a comparison with the VH structural repertoire of humans, described in detail elsewhere, could be properly accomplished. Our study shows the structural repertoire of mice to be dominated by canonical structure class 1-2 (approximately 60%), while in humans the dominant one is class 1-3 (approximately 40%). Analysis of the evolutionary relationships between human and mice suggest that this divergence may have a functional meaning. The implications of such findings are discussed.

Antibody-antigen recognition: a canonical structure paradigm.

The antibodies of known three-dimensional structure exhibit a definite number of conformations (canonical structures) for five of six hypervariable loops. In the present study it was found that approximately 85% of the immunoglobulin sequences analyzed fall into a small number of canonical structure combinations, representing only 3% of the total possible. These structures were classified into six distinct groups, depending on the type of antigen with which they interact. Within each loop, the positions responsible for maintaining these canonical structures show a use frequency of amino acids that fits an inverse power law, whereas the use frequency of the amino acids responsible for the detailed antigenic specificity follows an exponential distribution. We propose an evolutionary interpretation that connects these data, using the fact that the inverse power law is generated by statistical processes of the type that yield a wealth curve and the fact that exponential distribution is generated by processes that are not biased by past history.

Immunodominant synthetic peptides of Taenia crassiceps in murine and human cysticercosis.

The screening of a cDNA library of Taenia crassiceps revealed a clone designated KETc7 that induced high levels of protection against murine cysticercosis in previous experiments. The molecular structure of the deduced 100-amino acid sequence of the corresponding proline-rich polypeptide was studied to detect potentially immunologically active epitopes. Several candidate epitopes were identified, three of which were synthesized by solid-phase peptide synthesis and used as antigens in enzyme-linked immunosorbent assay (ELISA) for detection of specific antibodies in a selected panel of sera from mice infected with Taenia crassiceps and pigs infected with Taenia solium, as well as in the serum and cerebrospinal fluid of human patients with neurocysticercosis. The three peptides detected antibodies in serum from all infected mice. Seven of nine sera from patients with neurocysticercosis reacted strongly with peptide GK-3, and four of them with peptides GK-1 and GK-2. A lower reactivity was observed in sera from experimentally infected pigs. Peptide GK-3 reacted also with 45 out of 77 cerebrospinal fluids (CSF) from patients with confirmed neurocysticercosis and with 14 out of 68 CSF from control patients with other neurological disorders. This is the first report on synthetic peptides that are prominent in the humoral response of murine, porcine and human cysticercosis. Their identification implies finer molecular tools in the exploration of this form of host-parasite relationship, as well as hints to their application in immunodiagnosis and in vaccine design.

Frequency analysis of amino acids in the recognition regions of T-cell receptors.

In immunoglobulins (Igs), key amino acids in the Complementarity Determining Regions (CDR) are responsible for maintaining specific conformations called canonical structures. In T-cell receptors (TCRs), protein members of the Ig superfamily, the corresponding residues for maintaining these canonical structures have not been found. In previous studies we have found in Igs that the frequency of use of amino acids in some positions of the CDRs follows an inverse power law distribution, while the frequency of amino acids in the rest of the positions of the CDRs follows an exponential law distribution. The positions that follow the inverse power law distribution are precisely those involved in maintaining the canonical structures, while those positions for which the distribution fits the exponential distribution are those that should be properly involved in the recognition mechanism. In this paper, when the same analysis is applied to the use frequency of amino acids on the CDRs of TCRs, it is found that some positions that have been previously identified as having a structural role are those fitting the inverse power law. That finding combined with the cooperative or long-range interaction properties of systems that follow the inverse power law leads us to propose that the lack of determined key residues in certain positions is compensated by "equivalent' residues in other positions within the CDRs in order to maintain the canonical structures. Other positions that follow the exponential distribution are those which can be involved in the recognition process. These results coincide with a computer-generated model of TCR/peptide/MHC interaction previously published by the authors.

Canonical structure repertoire of the antigen-binding site of immunoglobulins suggests strong geometrical restrictions associated to the mechanism of immune recognition.

Is the structural repertoire of immunoglobulins free to adopt an almost infinite number of conformations to build the diversity of the immune response or does it take advantage of only a few conformations? In this paper we study this question by applying the canonical structure model to characterize the structural repertoire of immunoglobulins. The results found, indicate that only ten combinations out of the 300 possible different canonical structure classes (combinations of canonical structures), make up 87% of 381 sequences analyzed. This suggests that the structural repertoire of immunoglobulins is restricted to the preferential use of a small number of canonical structure classes. The possible functional significance of these results was studied by analyzing the correspondence between the observed canonical structural repertoire implicit in Ig sequences and the types of antigens recognized. Two different sets of canonical structure classes were distinguished: one with preference for some specific types of antigens like proteins, polysaccharides or haptens, and the other with multi-specific binding capabilities. Analysis of antibodies of known three-dimensional structure shows that for two specific classes, the canonical conformations of H2 and L1 determine the geometrical characteristics of the antigen-binding site, while at least in one multi-specific class, the changes in the general geometry of the antigen-binding site are produced by different conformations of H3. Implications of these results for the molecular recognition process mediated by immunoglobulins are discussed.

Molecular modeling of a T-cell receptor bound to a major histocompatibility complex molecule: implications for T-cell recognition.

The main functions of the T-cell receptor (TCR) involve its specific interaction with short and linear antigenic peptides bound to the major histocompatibility complex (MHC) molecules. In the absence of a 3D structure for TCR and for the TCR/peptide/MHC complex, several attempts to characterize the structural components of the TCR/peptide/MHC interaction have been made. However, this subject is still troublesome. In this paper a computer-based 3D model for a TCR/peptide/MHC complex (5C.C7/moth cytochrome c [MCC] peptide 93-103/I-Ek) was obtained. The complex surface shows a high complementarity between the 5C.C7 structure and the peptide/I-Ek molecule. The mapping of residues involved in the TCR/peptide/MHC interaction shows close agreement with mutational experiments (Jorgensen JL, Reay PA, Ehrich EW, Davis MM, 1992b, Annu Rev Immunol 10:835-873). Moreover, the results are consistent with a recent variability analysis of TCR sequences using three variability indexes (Almagro JC, Zenteno-Cuevas R, Vargas-Madrazo E, Lara-Ochoa F, 1995b, Int J Pept Protein Res 45:180-186). Accordingly, the 3D model of the 5C.C7/MCC peptide 93-103/I-Ek complex provides a framework to generate testable hypotheses about TCR recognition. Thus, starting from this model, the role played by each loop that forms the peptide/MHC binding site of the TCR is discussed.

Distributions of the use frequencies of amino acids in the hypervariable regions of immunoglobulins.

Two types of distributions for the frequencies of occurrence of amino acids in each position of hypervariable regions CDR-1 and CDR-2 were obtained for 2,000 immunoglobulins. The results show that some positions fit an inverse power-law distribution, while others fit an exponential-type distribution. As a result of comparison with structural data in the literature it is proposed that sites in which the frequency distribution fits the inverse power law are critical to maintaining canonical shapes of the recognition regions or are involved in modulating these canonical conformations, while those sites where the distribution fits the exponential law are those which should be exclusively involved in the recognition mechanism.

Structural repertoire in VH pseudogenes of immunoglobulins: comparison with human germline genes and human amino acid sequences.

In the pool of human immunoglobulin VH gene segments, pseudogenes amount to roughly 30% of the total number of genes. Some of them are highly conserved among unrelated individuals. These facts suggest a possible functional role for pseudogenes in the human immune response diversity. This paper intends to provide additional information about the structure of VH pseudogene sequences to evaluate the possible role of pseudogenes in the immune response. Mutations capable of altering framework stability in human VH pseudogenes were analyzed. Results indicate that VH pseudogenes are about 14 times as divergent as human VH functional germline genes on the one hand, and four times as divergent in the case of human VH amino acid sequences on the other. The high number of disruptive mutations in pseudogenes is an expected result because of the lack of functionality of these genes. In the second part of the work we analyze whether or not the same takes place in the positions that determine the existence of canonical structures in the hypervariable loops in VH pseudogenes. An extension of such analysis is applied to all species with reported VH pseudogenes. In contrast with results concerning framework positions, 69% of known human VH pseudogenes have canonical structures in the first hypervariable loop, while 48% do so in the second one. Comparison of these results with those found in human VH functional germline genes and human VH amino acid sequences shows that in the former as many as 100% and in the latter 96% have canonical structures. In VH amino acid sequences the result is similar to pseudogenes for H1. For H2, such value lies between the percentage of germline genes (96%) and the percentage of pseudogenes (48%). The possible significance of the existence of canonical structures in the hypervariable loops of VH pseudogenes is discussed.

Variability analysis of the T-cell receptors using three variability indexes.

In the absence of a three-dimensional structure for TCR molecules, several attempts to identify their hypervariable regions by variability methods have been made; this subjects is still troublesome. In this paper three different variability indexes were used: (i) the Kabat index, which is the classical measure of sequence variability, (ii) the modified Kabat index, successfully used in the beta-chain of T-cell receptors and (iii) an information-theoretical entropy concept, recently proposed as an improved measure of the variability. In order to identify the hypervariable regions in the TCR sequences, a Fourier filtering was applied on each variability profile. Results show that the three variability indexes have distinct resolutions for different levels of variability. Thus, the simultaneous use of these indexes compensates for the deficiency of any one of them in estimating variability. Applying the Fourier filtering, it is found that the hypervariable regions here identified, roughly coincide with the defined CDR-2 and CDR-3 in TCR by analogy with Ig. However, no hypervariable in the CDR-1 of alpha- and beta-chains was found. The study on the influence of sample size in variability analysis, indicates that results are independent of the sample size. Considering current structural models of TCR-peptide-MHC interaction, one can suggest that the low-variability characteristics of these regions is inherently related to the interaction with relatively conserved region on the alpha-helices of MHC.

VIR: a computational tool for analysis of immunoglobulin sequences.

In this paper a microcomputer software named VIR (Variable domains of the Immune Receptors) is reported. This package can be used in sequence studies of immunoglobulin variable domains. The main features of the VIR software in the sequences management are: (1) ease of information recovery/extraction from amino acid sequences; and (2) its capability to obtain multiple sequence alignments with predefined characteristics (i.e. specie and/or specificity). As an analytical tool, the VIR package employs such multiple sequence alignments to compute: (1) tables showing amino acid frequencies; (2) three variability indexes; (3) identity matrices; (4) random samples; and (5) sequences with possible canonical structures. Thus the software reported here is proposed as a useful tool to carry out detailed studies of immunoglobulin variable domains.

Patterns in the complementary determining regions of immunoglobulins (CDRs).

An analysis of the frequency of use of amino acids on the CDR-1 and CDR-2 of 1500 immunoglobulins showed that the frequencies of amino acids in different positions could be fitted by two types of distribution. For some positions the frequencies were fitted by an inverse power law and for other positions by an exponential distribution. In order to see whether the more frequently used amino acids for specific positions had physicochemical properties or attributes in common, they were clustered using an algorithm normally applied to artificial intelligence problems. It was found that the amino acids in those positions fitted by the inverse power law have similar hydrophobicity and volume, which are commonly attributes of amino acids in structural positions. Thus, if these positions are critical to maintaining the structural features of the CDR domains, the rest of the positions should be either properly involved in the recognition process or irrelevant. The frequencies of amino acids in these recognition positions were fitted by the exponential law, and it was found by the clustering analysis that these amino acids share properties of a more general type, such as capability of forming hydrogen bonds, polarity, etc. This suggests that at least part of the recognition mechanism requires general properties rather than specific amino acids. Amino acids sharing the required attributes for each one of these positions are then used with random frequency.