Antibody informatics for drug discovery.

More and more antibody therapeutics are being approved every year, mainly due to their high efficacy and antigen selectivity. However, it is still difficult to identify the antigen, and thereby the function, of an antibody if no other information is available. There are obstacles inherent to the antibody science in every project in antibody drug discovery. Recent experimental technologies allow for the rapid generation of large-scale data on antibody sequences, affinity, potency, structures, and biological functions; this should accelerate drug discovery research. Therefore, a robust bioinformatic infrastructure for these large data sets has become necessary. In this article, we first identify and discuss the typical obstacles faced during the antibody drug discovery process. We then summarize the current status of three sub-fields of antibody informatics as follows: (i) recent progress in technologies for antibody rational design using computational approaches to affinity and stability improvement, as well as ab-initio and homology-based antibody modeling; (ii) resources for antibody sequences, structures, and immune epitopes and open drug discovery resources for development of antibody drugs; and (iii) antibody numbering and IMGT. Here, we review "antibody informatics," which may integrate the above three fields so that bridging the gaps between industrial needs and academic solutions can be accelerated. This article is part of a Special Issue entitled: Recent advances in molecular engineering of antibody.

Kotai Antibody Builder: automated high-resolution structural modeling of antibodies.

MOTIVATION: Kotai Antibody Builder is a Web service for tertiary structural modeling of antibody variable regions. It consists of three main steps: hybrid template selection by sequence alignment and canonical rules, 3D rendering of alignments and CDR-H3 loop modeling. For the last step, in addition to rule-based heuristics used to build the initial model, a refinement option is available that uses fragment assembly followed by knowledge-based scoring. Using targets from the Second Antibody Modeling Assessment, we demonstrate that Kotai Antibody Builder generates models with an overall accuracy equal to that of the best-performing semi-automated predictors using expert knowledge. AVAILABILITY AND IMPLEMENTATION: Kotai Antibody Builder is available at http://kotaiab.org CONTACT: standley@ifrec.osaka-u.ac.jp.

BACTERIAL IDENTIFICATION USING SSRA ENCODING TRANSFER-MESSENGER RNA.

Abstract. Ribosomal DNA (rDNA) sequences are widely used for phylogenetic and bacterial identification. However, rDNA of different species often reveals similar or identical same sequences. This study employed the bacterial stable small RNA (ssrA) gene encoding transfer-messenger RNA (tmRNA) as a tool for identification of Staphylococcus aureus, Enterococcus spp, Pseudomonas spp and Enterobacteriaceae from clinical isolates as representative groups using PCR and species specific primers. The method correctly identified 11 standard strains and 99 clinical isolates. Quantitative PCR revealed a limit of detection of 10(-5) µg of DNA for S. aureus and Enterococcus spp, and 10(-6) µg for Pseudomonas spp and Enterobacteriaceae. Further studies with a greater number of bacteria especially from clinical samples will need to be undertaken before this bacterial molecular marker can be applied in a clinical setting.

High-resolution modeling of antibody structures by a combination of bioinformatics, expert knowledge, and molecular simulations.

In the second antibody modeling assessment, we used a semiautomated template-based structure modeling approach for 11 blinded antibody variable region (Fv) targets. The structural modeling method involved several steps, including template selection for framework and canonical structures of complementary determining regions (CDRs), homology modeling, energy minimization, and expert inspection. The submitted models for Fv modeling in Stage 1 had the lowest average backbone root mean square deviation (RMSD) (1.06 Å). Comparison to crystal structures showed the most accurate Fv models were generated for 4 out of 11 targets. We found that the successful modeling in Stage 1 mainly was due to expert-guided template selection for CDRs, especially for CDR-H3, based on our previously proposed empirical method (H3-rules) and the use of position specific scoring matrix-based scoring. Loop refinement using fragment assembly and multicanonical molecular dynamics (McMD) was applied to CDR-H3 loop modeling in Stage 2. Fragment assembly and McMD produced putative structural ensembles with low free energy values that were scored based on the OSCAR all-atom force field and conformation density in principal component analysis space, respectively, as well as the degree of consensus between the two sampling methods. The quality of 8 out of 10 targets improved as compared with Stage 1. For 4 out of 10 Stage-2 targets, our method generated top-scoring models with RMSD values of less than 1 Å. In this article, we discuss the strengths and weaknesses of our approach as well as possible directions for improvement to generate better predictions in the future.

Modulation of O-antigen chain length by the wzz gene in Escherichia coli O157 influences its sensitivities to serum complement.

Escherichia coli O157 strains belonging to a distinct lineage and expressing different O-antigen (Oag) lengths were isolated. Although the function of wzz in E. coli has not been adequately investigated, this gene is known to be associated with regulation of Oag length. Using E. coli O157:H7 ATCC43888 (wild-type), several wzz mutants of E. coli O157, including a wzz deletion mutant, were generated and the relationship between the length of Oag modulated by the wzz gene and sensitivities to serum complement investigated. SDS-PAGE, immunoblot analyses and sensitivity tests to human serum complement were performed on these strains. The lengths of the O157-antigen could be modulated by the wzz gene mutations and were classified into long, intermediate and short groups. The short chain mutant was more serum sensitive than the wild-type strain and the other wzz mutants (P < 0.001). In conclusion, Oag chain length modulated by the wzz gene in E. coli O157 influences its sensitivities to serum complement. The present findings suggest that E. coli O157 strains with intermediate or long length Oag chains might show greater resistance to serum complement than those with short chains.

Increasing the melting temperature of VHH with the in silico free energy score.

VHH, the antigen-binding fragment of a heavy chain-only antibody, is a useful component of antibody-based therapeutics. Thermal stability, represented by the melting temperature (Tm), is one of the key factors affecting the developability of antibody-based therapeutics. In this study, we examined whether the in silico free energy score dStability can be used to design mutants with improved Tm compared to the anti-lysozyme VHH, D3-L11. After verifying that exhaustive mutagenesis was inefficient for improving Tm, we performed a two-round rational approach that combined dStability calculations with a small number of experiments. This method improved the Tm by more than 5 °C in several single mutants including A79I. It reduced the affinity for the antigen by less than 1.6-fold. We speculate that stabilization of A79I required exquisite compatibility among neighboring residues to fill in the internal cavity in the protein. Given that we identified only one mutation that could simultaneously improve Tm and almost maintain affinity, we concluded that achieving both is extremely difficult, even with single mutations that are not located in the paratope. Therefore, we recommend using a variety of approaches when trying to achieve such a feat. Our method will be a useful complementary approach to other existing methods.

Computational design, construction, and characterization of a set of specificity determining residues in protein-protein interactions.

Proteins interact with different partners to perform different functions and it is important to elucidate the determinants of partner specificity in protein complex formation. Although methods for detecting specificity determining positions have been developed previously, direct experimental evidence for these amino acid residues is scarce, and the lack of information has prevented further computational studies. In this article, we constructed a dataset that is likely to exhibit specificity in protein complex formation, based on available crystal structures and several intuitive ideas about interaction profiles and functional subclasses. We then defined a "structure-based specificity determining position (sbSDP)" as a set of equivalent residues in a protein family showing a large variation in their interaction energy with different partners. We investigated sequence and structural features of sbSDPs and demonstrated that their amino acid propensities significantly differed from those of other interacting residues and that the importance of many of these residues for determining specificity had been verified experimentally.

Systematic classification of CDR-L3 in antibodies: implications of the light chain subtypes and the VL-VH interface.

Antibody modeling is widely used for the analysis of antibody-antigen interactions and for the design of potent antibody drugs. The antibody combining site is composed of six complementarity determining regions (CDRs). The CDRs, except for CDR-H3, which is the most diverse CDR, form limited numbers of canonical structures, which can be identified from the amino acid sequences. A method to classify the CDR-H3 structure from its amino acid sequence was previously proposed. However, since those CDR structures were classified, many more antibody crystal structures have been determined. We performed systematic analyses of the CDR-L3 structures and found novel canonical structures, and we also classified a previously identified canonical structure into two subtypes. In addition, two differently defined canonical structures in the kappa and lambda subtypes were classified into the same canonical structure. We also identified a key residue in CDR-L3, which determines the conformation of CDR-H3. Several analyses of CDR-L3 loops longer than nine residues were performed. These new findings should be useful for structural modeling and are eventually expected to accelerate the design of antibody drugs.

Exploring designability of electrostatic complementarity at an antigen-antibody interface directed by mutagenesis, biophysical analysis, and molecular dynamics simulations.

Antibodies protect organisms from a huge variety of foreign antigens. Antibody diversity originates from both genetic and structural levels. Antigen recognition relies on complementarity between antigen-antibody interfaces. Recent methodological advances in structural biology and the accompanying rapid increase of the number of crystal structures of proteins have enabled atomic-level manipulation of protein structures to effect alterations in function. In this study, we explored the designability of electrostatic complementarity at an antigen-antibody interface on the basis of a crystal structure of the complex. We designed several variants with altered charged residues at the interface and characterized the designed variants by surface plasmon resonance, circular dichroism, differential scanning calorimetry, and molecular dynamics simulations. Both successes and failures of the structure-based design are discussed. The variants that compensate electrostatic interactions can restore the interface complementarity, enabling the cognate antigen-antibody binding. Retrospectively, we also show that these mutational effects could be predicted by the simulations. Our study demonstrates the importance of charged residues on the physical properties of this antigen-antibody interaction and suggests that computational approaches can facilitate design of antibodies that recognize a weakly immunogenic antigen.

Structural classification of CDR-H3 revisited: a lesson in antibody modeling.

Among the six complementarity-determining regions (CDRs) in the variable domains of an antibody, the third CDR of the heavy chain (CDR-H3), which lies in the center of the antigen-binding site, plays a particularly important role in antigen recognition. CDR-H3 shows significant variability in its length, sequence, and structure. Although difficult, model building of this segment is the most critical step in antibody modeling. Since our first proposal of the "H3-rules," which classify CDR-H3 structure based on amino acid sequence, the number of experimentally determined antibody structures has increased. Here, we revise these H3-rules and propose an improved classification scheme for CDR-H3 structure modeling. In addition, we determine the common features of CDR-H3 in antibody drugs as well as discuss the concept of "antibody druggability," which can be applied as an indicator of antibody evaluation during drug discovery.

[Interface between bioinformatics and docking study].

We describe the prospects of bioinformatics for drug discovery and discuss the current status, problems, and future direction of the interface between bioinformatics and docking studies. We also describe our recent work on sequence and structure analysis using the guanidino-modifying enzymes superfamily as a good example.

The guanidino-group modifying enzymes: structural basis for their diversity and commonality.

The guanidino-group modifying enzyme (GME) superfamily contains many drug targets, including metabolic enzymes from pathogenic microorganisms as well as key regulatory proteins from higher eukaryotes. These enzymes, despite their diverse sequences, adopt the common alpha/beta propeller fold and catalyze the modification of (methylated) guanidino groups. Our structural superposition and structure-based alignment for the GMEs have identified key residues that are involved in the catalysis and substrate binding. We have shown that conserved guanidino-carboxyl interactions are utilized in two different ways; the acidic residues in the catalytic site form hydrogen bonds to the substrate guanidino group, and the enzyme Arg residues at several key positions recognize the carboxyl group of the substrate and fix its orientation. Based on this observation, we have proposed rules for classifying the GME sequences and predicting their molecular function from the conservation of the key acidic and Arg residues. Other novel motifs have been identified, which involve residues that are not in direct contact with the substrate but are likely to stabilize the active-site conformation through hydrogen-bonding networks. In addition, we have examined the domain architecture of the GMEs. Although most members consist of a single catalytic domain, fold recognition analysis has identified a likely bifunctional enzyme from a cyanobacterium. It has also revealed common immunoglobulin-like beta-sandwich domains found in the enzymes that recognize protein substrates. These findings will be useful for predicting the precise mechanism of action for potential novel targets and designing therapeutic compounds against them.

Comparative study on the structure and cytopathogenic activity of HIV Vpr/Vpx proteins.

The three-dimensional (3-D) structure of human immunodeficiency virus type 2 (HIV-2) Vpr/Vpx was predicted by homology modeling based on the NMR structure of human immunodeficiency virus type 1 (HIV-1) Vpr. The three proteins similarly have three major amphipathic alpha-helices. In contrast to HIV-1 Vpr, Vpr/Vpx of HIV-2 have a long N-terminal loop and clustered prolines in the second half of the C-terminal loop. HIV-2 Vpx uniquely contains a long region between the second and third major helices, and bears several glycines in the first half of the C-terminal loop. Instead of the glycines, there is a group of hydrophilic amino acids and arginines in the corresponding regions of the two Vprs. To compare the cytopathogenic potentials of HIV-1 Vpr and HIV-2 Vpr/Vpx, we examined the production of luciferase as a marker of cell damage. We further analyzed the characteristics of cells transduced with vpr/vpx genes driven by an inducible promoter. The results obtained clearly show that structurally similar, but distinct, HIV Vpr/Vpx proteins are detrimental to target cells.

High titers of IgG antibodies against Plasmodium falciparum serine repeat antigen 5 (SERA5) are associated with protection against severe malaria in Ugandan children.

Plasmodium falciparum serine repeat antigen (SERA5) is a promising asexual blood stage malaria candidate vaccine. However, there is a paucity of information about natural immune responses to SERA5 in children from malaria-endemic regions. We undertook a hospital-based case-control study of severe malaria in Apac District, Northern Uganda, in children 6-59 months of age. The commonest symptoms observed in children with severe malaria (SM) were respiratory distress (53.4%) and prostration (40.4%) followed by circulatory collapse (7.4%), severe anemia (Hb < 5 g/dL, 7.0%), and seizures (2.6%). None of the SM children had impaired consciousness, coma, or cerebral malaria. We measured serum IgG antibodies using a recombinant construct of SERA5 (SE36) in enzyme-linked immunosorbent assays. High titers of IgG anti-SE36 were associated with protection against severe malaria in children under 5 years old.

Antibody titres and boosting after natural malaria infection in BK-SE36 vaccine responders during a follow-up study in Uganda.

The malaria vaccine BK-SE36 is a recombinant protein (SE36) based on the Honduras 1 serine repeat antigen-5 of Plasmodium falciparum, adsorbed to aluminium hydroxide gel. The phase Ib trial in Uganda demonstrated the safety and immunogenicity of BK-SE36. Ancillary analysis in the follow-up study of 6-20 year-old volunteers suggest significant differences in time to first episodes of clinical malaria in vaccinees compared to placebo/control group. Here, we aimed to get further insights into the association of anti-SE36 antibody titres and natural P. falciparum infection. Children who received BK-SE36 and whose antibody titres against SE36 increased by ≥1.92-fold after vaccination were categorised as responders. Most responders did not have or only had a single episode of natural P. falciparum infection. Notably, responders who did not experience infection had relatively high anti-SE36 antibody titres post-second vaccination compared to those who were infected. The anti-SE36 antibody titres of the responders who experienced malaria were boosted after infection and they had lower risk of reinfection. These findings show that anti-SE36 antibody titres induced by BK-SE36 vaccination offered protection against malaria. The vaccine is now being evaluated in a phase Ib trial in children less than 5 years old.

Prediction of the structure and function of AstA and AstB, the first two enzymes of the arginine succinyltransferase pathway of arginine catabolism.

Arginine succinyltransferase and succinylarginine dihydrolase catalyze the first two steps of arginine catabolism by the arginine succinyltransferase pathway. This route is the only major arginine catabolic pathway in Escherichia coli including its pathogenic strains O157 and CFT073. We have used fold recognition tools and identified novel homologies between each of these two enzymes and proteins of known three-dimensional structure: arginine succinyltransferase belongs to the acyl-CoA N-acyltransferase superfamily and succinylarginine dihydrolase belongs to the amidinotransferase superfamily. These findings shed light on the structures, catalytic mechanisms and evolution of diverse enzymes involved in arginine catabolism.

Affinity improvement of a therapeutic antibody by structure-based computational design: generation of electrostatic interactions in the transition state stabilizes the antibody-antigen complex.

The optimization of antibodies is a desirable goal towards the development of better therapeutic strategies. The antibody 11K2 was previously developed as a therapeutic tool for inflammatory diseases, and displays very high affinity (4.6 pM) for its antigen the chemokine MCP-1 (monocyte chemo-attractant protein-1). We have employed a virtual library of mutations of 11K2 to identify antibody variants of potentially higher affinity, and to establish benchmarks in the engineering of a mature therapeutic antibody. The most promising candidates identified in the virtual screening were examined by surface plasmon resonance to validate the computational predictions, and to characterize their binding affinity and key thermodynamic properties in detail. Only mutations in the light-chain of the antibody are effective at enhancing its affinity for the antigen in vitro, suggesting that the interaction surface of the heavy-chain (dominated by the hot-spot residue Phe101) is not amenable to optimization. The single-mutation with the highest affinity is L-N31R (4.6-fold higher affinity than wild-type antibody). Importantly, all the single-mutations showing increase affinity incorporate a charged residue (Arg, Asp, or Glu). The characterization of the relevant thermodynamic parameters clarifies the energetic mechanism. Essentially, the formation of new electrostatic interactions early in the binding reaction coordinate (transition state or earlier) benefits the durability of the antibody-antigen complex. The combination of in silico calculations and thermodynamic analysis is an effective strategy to improve the affinity of a matured therapeutic antibody.

Hematological and biochemical data obtained in rural northern Uganda.

Reference intervals for common hematological and clinical chemistry parameters constitute an important basis for health care. Moreover, with increasing priority in drug and vaccine development for infectious diseases in Africa, the first priority is the safety evaluation and tolerability of the candidate interventions in healthy populations. To accurately assess health status and address adverse events, clinical reference intervals in the target population are necessary. We report on hematological and biochemical indices from healthy volunteers who participated in a clinical trial in Lira, northern Uganda. Median and nonparametric 95% percentiles on five hematology and 15 biochemistry analytes are shown. Although most hematological analytes conformed to reported reference intervals and trends in Africa, literature review from different African countries highlight the need for a region-specific children reference interval that can be appropriate for the population.

Association of naturally acquired IgG antibodies against Plasmodium falciparum serine repeat antigen-5 with reduced placental parasitemia and normal birth weight in pregnant Ugandan women: a pilot study.

Plasmodium falciparum infection during pregnancy contributes substantially to malaria burden in both mothers and offspring. Analysis of naturally acquired immune responses that confer protection against parasitemia and clinical disease is important to guide vaccine evaluation as well as identify immune correlates. Unfortunately, few studies have addressed the relationship between immune responses to malaria vaccine candidate antigens and protection against adverse effects on pregnant women and newborn birth weight. This study examines the relationship of maternal antibody responses to serine repeat antigen-5 (SE36) and merozoite surface protein-1 (MSP119 and MSP142) with placental parasitemia and birth weight. In a peri-urban setting in Uganda, pregnant women without placental parasites have high median ODs for antibodies against SE36 (P<0.001). Naturally acquired anti-SE36 IgG was most prevalent in women without placental parasitemia (P<0.001). Furthermore, pregnant women with significantly high levels of anti-SE36 IgG delivered babies with normal birth weights (P<0.001). That antibody to SE36 was associated with both a reduced risk of placental parasitemia and resulting normal birth weight in newborns suggests some protective role. In contrast, although antibody to MSP142 was also associated with reduced placental parasitemia and immune responses to both MSP119 and MSP142 may be of importance, there was no association between anti-MSP119 antibodies and infant birth weight outcomes. This study highlights the need for conducting further studies to investigate the association of antibodies against SE36 and outcomes of malaria infection in pregnant women.

Phase 1b randomized trial and follow-up study in Uganda of the blood-stage malaria vaccine candidate BK-SE36.

BACKGROUND: Up to now a malaria vaccine remains elusive. The Plasmodium falciparum serine repeat antigen-5 formulated with aluminum hydroxyl gel (BK-SE36) is a blood-stage malaria vaccine candidate that has undergone phase 1a trial in malaria-naive Japanese adults. We have now assessed the safety and immunogenicity of BK-SE36 in a malaria endemic area in Northern Uganda. METHODS: We performed a two-stage, randomized, single-blinded, placebo-controlled phase 1b trial (Current Controlled trials ISRCTN71619711). A computer-generated sequence randomized healthy subjects for 2 subcutaneous injections at 21-day intervals in Stage1 (21-40 year-olds) to 1-mL BK-SE36 (BKSE1.0) (n = 36) or saline (n = 20) and in Stage2 (6-20 year-olds) to BKSE1.0 (n = 33), 0.5-mL BK-SE36 (BKSE0.5) (n = 33), or saline (n = 18). Subjects and laboratory personnel were blinded. Safety and antibody responses 21-days post-second vaccination (Day42) were assessed. Post-trial, to compare the risk of malaria episodes 130-365 days post-second vaccination, Stage2 subjects were age-matched to 50 control individuals. RESULTS: Nearly all subjects who received BK-SE36 had induration (Stage1, n = 33, 92%; Stage2, n = 63, 96%) as a local adverse event. No serious adverse event related to BK-SE36 was reported. Pre-existing anti-SE36 antibody titers negatively correlated with vaccination-induced antibody response. At Day42, change in antibody titers was significant for seronegative adults (1.95-fold higher than baseline [95% CI, 1.56-2.43], p = 0.004) and 6-10 year-olds (5.71-fold [95% CI, 2.38-13.72], p = 0.002) vaccinated with BKSE1.0. Immunogenicity response to BKSE0.5 was low and not significant (1.55-fold [95% CI, 1.24-1.94], p = 0.75). In the ancillary analysis, cumulative incidence of first malaria episodes with ≥5000 parasites/µL was 7 cases/33 subjects in BKSE1.0 and 10 cases/33 subjects in BKSE0.5 vs. 29 cases/66 subjects in the control group. Risk ratio for BKSE1.0 was 0.48 (95% CI, 0.24-0.98; p = 0.04). CONCLUSION: BK-SE36 is safe and immunogenic. The promising potential of BK-SE36, observed in the follow-up study, warrants a double-blind phase 1/2b trial in children under 5 years. TRIAL REGISTRATION: Controlled-Trials.com ISRCTN71619711.