Aligning,analyzing, and visualizing sequences for antibody engineering: Automated recognition of immunoglobulin variable region features

The analysis and comparison of large numbers of immunoglobulin (Ig) sequences that arise during an antibody selection campaign can be time‐consuming and tedious. Typically, the identification and annotation of framework as well as complementarity‐determining regions (CDRs) is based on multiple sequence alignments using standardized numbering schemes, which allow identification of equivalent residues among different family members but often necessitate expert knowledge and manual intervention. Moreover, due to the enormous length variability of some CDRs the benefit of conventional Ig numbering schemes is limited and the calculation of correct sequence alignments can become challenging. Whereas, in principle, a well established set of rules permits the assignment of CDRs from the amino acid sequence alone, no currently available sequence alignment editor provides an algorithm to annotate new Ig sequences accordingly. Here we present a unique pattern matching method implemented into our recently developed ANTIC ALIgN editor that automatically identifies all hypervariable and framework regions in experimentally elucidated antibody sequences using so‐called “regular expressions.” By combination of this widely supported software syntax with the unique capabilities of real‐time aligning, editing and analyzing extended sets of amino acid and/or nucleotide sequences simultaneously on a local workstation, ANTIC ALIgN provides a powerful utility for antibody engineering. Proteins 2016; 85:65–71. © 2016 Wiley Periodicals, Inc.     

High‐Affinity Antibody Detection with a Bivalent Circularized Peptide Containing Antibody‐Binding Domains

Direct chemical labeling of antibody produces molecules with poorly defined modifications. Use of a small antibody‐binding protein as an adapter can simplify antibody functionalization by forming a specific antibody‐bound complex and introducing site‐specific modifications. To stabilize a noncovalent antibody complex that may be used without chemical crosslinking, a bivalent antibody‐binding protein is engineered with an improved affinity of interaction by joining two Z domains with a conformationally flexible linker. The linker is essential for the increase in affinity because it allows simultaneous binding of both domains. The molecule is further circularized using a split intein, creating a novel adapter protein (“lasso”), which binds human immunoglobulin G1 (IgG1) with K _D = 0.53 n m and a dissociation rate that is 55‐ to 84‐fold slower than Z. The lasso contains a unique cysteine for conjugation with a reporter and may be engineered to introduce other functional groups, including a biotin tag and protease recognition sequences. When used in enzyme‐linked immunosorbent assay (ELISA), the lasso generates a stronger reporter signal compared to a secondary antibody and lowers the limit of detection by 12‐fold. The small size of the lasso and a long half‐life of dissociation make the peptide a useful tool in antibody detection and immobilization.     

Multispecificity of a recombinant anti‐ras monoclonal antibody

Recombinant monoclonal antibodies (Ab's) have widespread application as research tools, diagnostic reagents and as biotherapeutics. Whilst studying the cellular molecular switch protein m‐ras, a recombinant monoclonal antibody to m‐ras was generated for use as a research tool. Antibody genes from a single rabbit B cell secreting IgG to an m‐ras specific peptide sequence were expressed in mammalian cells, and monoclonal rabbit IgG binding was characterized by ELISA and peptide array blotting. Although the monoclonal Ab was selected for specificity to m‐ras peptide, it also bound to both recombinant full‐length m‐ras and h‐ras proteins. The cross‐reactive binding of the monoclonal Ab to h‐ras was defined by peptide array blot revealing that the Ab showed preference for peptide sequences containing multiple positively charged amino acid residues. These data reinforce the concept of antibody multispecificity through multiple interactions of the Ab paratope with diverse polypeptides. They also emphasize the importance of immunogen and Ab selection processes when generating recombinant monoclonal Ab's.     

Ⅱ型志贺毒素单克隆抗体可变区基因克隆及单链抗体的表达和功能鉴定

目的从分泌抗肠出血性大肠埃希菌Ⅱ型志贺毒素中和单克隆抗体杂交瘤细胞株S2C4中克隆抗体可变区基因,构建单链抗体（ScFv）,进行原核表达,并对其功能进行鉴定。方法从杂交瘤细胞株S2C4中提取总RNA,逆转录成cDNA。在cDNA3’-OH末端添加poly．G。PCR扩增包括5’非翻译区和信号肽序列在内的抗体重、轻链可变区基因VH和VL,PCR产物装入T—A载体测序。根据测序结果,设计引物分别扩增VH和VL编码区,再通过重叠PCR,在VH和VL．编码区基因之间引入连接链,构建Scn基因,并克隆到表达载体pComb3xSS中。重组载体导入E．coliTop10F’进行表达,重组蛋白经纯化后,分别用ELISA和动物保护性实验鉴定其生物学活性。结果VH和VL编码区基因全长分别为396bp和378bp,ScFv基因能在大肠埃希菌中高效表达,表达产物的分子量为34000,用NiSO4亲和层析柱成功纯化。功能性实验表明纯化的重组蛋白可以与Stx2毒素有效结合,能保护动物抵御毒素分子的攻击。结论成功地克隆S2C4单抗可变区基因,并构建、表达其单链抗体ScFv,为下一步进行该抗体人源化奠定实验基础。     

Expression of antibodies using single‐open reading frame vector design and polyprotein processing from mammalian cells

We describe a novel polyprotein precursor‐based approach to express antibodies from mammalian cells. Rather than expressing heavy and light chain proteins from separate expression units, the antibody heavy and light chains are contained in one single‐open reading frame (sORF) separated by an intein gene fused in frame. Inside mammalian cells this ORF is transcribed into a single mRNA, and translated into one polypeptide. The antibody heavy and light chains are separated posttranslationally, assembled into the functional antibody molecule, and secreted into culture medium. It is demonstrated that Pol I intein from P. horikoshii mediates protein splicing and cleavage reactions in mammalian cells, in the context of antibody heavy and light chain amino acid sequences. To allow the separation of antibody heavy chain, light chain, and the intein, we investigated a number of intein mutations designed to inhibit intein‐mediated splicing but preserve cleavage reactions. We have also designed constructs in which the signal peptide downstream from intein has altered hydrophobicity. The use of some of these mutant constructs resulted in more efficient antibody secretion, highlighting areas that can be further explored in improving such an expression system. An antibody secreted using one of the sORF constructs was characterized. This antibody has correct N‐terminal sequences for both of its heavy and light chains, correct heavy and light chain MW as well as intact MW as measured by mass spectrometry. Its affinity to antigen, as measured by surface plasmon resonance (SPR), is indistinguishable from that of the same antibody produced using conventional method. © 2009 American Institute of Chemical Engineers Biotechnol. Prog., 2009     

VIRS: A visual tool for identifying restriction sites in multiple DNA sequences

VIRS (A visual tool for identifying restriction sites in multiple DNA sequences) is an interactive web‐based program designed for restriction endonuclease cut sites prediction and visualization. It can afford to analyze multiple DNA sequences simultaneously and produce visual restriction maps with several useful options intended for users' customization. These options also perform in‐depth analysis of the restriction maps, such as providing virtual electrophoretic result for digested fragments. Different from other analytical tools, VIRS not only displays visual outputs but also provides the detailed properties of restriction endonucleases that are commercially available. All the information of these enzymes is stored in our internal database, which is updated monthly from the manufacturers' web pages. It is freely available online at http://bis.zju.edu.cn/virs/index.html . © 2009 American Institute of Chemical Engineers Biotechnol. Prog., 2009     

A computational approach for studying antibody–antigen interactions without prior structural information: the anti‐testosterone binding antibody as a case study

Given the increasing exploitation of antibodies in different contexts such as molecular diagnostics and therapeutics, it would be beneficial to unravel the atomistic level properties of antibody‐antigen complexes with the help of computational modeling. Thus, here we have studied the feasibility of computational tools to gather atomic scale information regarding the antibody‐antigen complexes solely starting from an amino acid sequence. First, we constructed a homology model for the anti‐testosterone binding antibody based on the knowledge based classification of complementary determining regions (CDRs) and implicit solvent molecular dynamics simulations. To further examine whether the generated homology model is suitable for studying antibody‐antigen interactions, docking calculations were carried out followed by binding free‐energy simulations. Our results indicate that with the antibody modeling approach presented here it is possible to construct accurate homology models for antibodies which correctly describes the antibody‐antigen interactions, and produces absolute binding free‐energies that are comparable with experimental values. In addition, our simulations suggest that the conformations of complementary determining regions (CDRs) may considerably change from the X‐ray configuration upon solvation. In conclusion, here we have introduced an antibody modeling workflow that can be used in studying the interactions between antibody and antigen solely based on an amino acid sequence, which in turn provides novel opportunities to tune the properties of antibodies in different applications. Proteins 2017; 85:322–331. © 2016 Wiley Periodicals, Inc.     

Resolving self-association of a therapeutic antibody by formulation optimization and molecular approaches

A common challenge encountered during development of high concentration monoclonal antibody formulations is preventing self-association. Depending on the antibody and its formulation, self-association can be seen as aggregation, precipitation, opalescence or phase separation. Here we report on an unusual manifestation of self-association, formation of a semi-solid gel or “gelation." Therapeutic monoclonal antibody C4 was isolated from human B cells based on its strong potency in neutralizing bacterial toxin in animal models. The purified antibody possessed the unusual property of forming a firm, opaque white gel when it was formulated at concentrations >30 mg/mL and the temperature was <6°C. Gel formation was reversible with temperature. Gelation was affected by salt concentration or pH, suggesting an electrostatic interaction between IgG monomers. A comparison of the C4 amino acid sequences to consensus germline sequences revealed differences in framework regions. A C4 variant in which the framework sequence was restored to the consensus germline sequence did not gel at 100 mg/mL at temperatures as low as 1°C. Additional genetic analysis was used to predict the key residue(s) involved in the gelation. Strikingly, a single substitution in the native antibody, replacing heavy chain glutamate 23 with lysine (E23K), was sufficient to prevent gelation. These results indicate that the framework region is involved in intermolecular interactions. The temperature dependence of gelation may be related to conformational changes near glutamate 23 or the regions it interacts with. Molecular engineering of the framework can be an effective approach to resolve the solubility issues of therapeutic antibodies.     

甲型H1N1流感病毒血凝素单克隆抗体轻链和重链可变区基因的巢式PCR扩增及序列分析

目的：采用巢式PCR对甲型H1N1流感病毒血凝素单克隆抗体的轻链和重链基因进行扩增,对获得的基因进行序列分析,并找出克隆鼠Igκ轻链和重链可变区基因的通用方法。方法：设计22对扩增鼠Igκ轻链可变区和重链可变区基因的引物,对6株鼠抗人甲型H1N1流感病毒血凝素单克隆抗体的轻链和重链可变区基因进行克隆并测序,与NCBI公布的鼠免疫球蛋白序列比对分析。结果：巢式PCR方法可以有效避免单克隆抗体克隆过程的假基因,并且得到的单克隆抗体的氨基酸序列均符合鼠免疫球蛋白可变区特征。结论：建立了克隆鼠免疫球蛋白轻链和重链可变区基因的通用方法,为后期克隆鼠源性单克隆抗体的可变区基因提供了基础,并为研究甲型H1N1流感病毒血凝素与抗体的结合位点提供了实验数据。     

蛋白质内含子在特定氨基酸序列中的剪接

基因工程技术已经被广泛应用于抗体的生产。但是由于抗体的分子量较大,导致合成抗体较为困难。蛋白质内含子是前体蛋白质中的一段氨基酸序列,能够将自身剪切出来,并将两端的外显子连接形成成熟的蛋白质。将抗体的Fab(antigen binding fragment)和Fc(crystalline fragment)分别与蛋白质内含子(intein) 的N端(IN)和C端(IC)融合表达,利用蛋白质内含子的剪接功能,可形成完整的抗体分子。KSCDKTH是存在于抗体铰链区(hinge region)的一段氨基酸序列,如果在KSCDKTH序列中筛选到高效剪接的蛋白质内含子,即可通过蛋白质剪接,将抗体分子的Fab和Fc剪接形成完整抗体。本文筛选发现,Ssp DnaX的3种断裂蛋白质内含子(S0, S1, S11)具有在KSCDKTH序列中高效剪接的能力,这一研究结果为抗体的剪接合成提供了可行性。     

Intracellular antibodies and cancer: New technologies offer therapeutic opportunities

Since the realisation that the antigen‐binding regions of antibodies, the variable (V) regions, can be uncoupled from the rest of the molecule to create fragments that recognise and abrogate particular protein functions in cells, the use of antibody fragments inside cells has become an important tool in bioscience. Diverse libraries of antibody fragments plus in vivo screening can be used to isolate single chain variable fragments comprising VH and VL segments or single V‐region domains. Some of these are interfering antibody fragments that compete with protein‐protein interactions, providing lead molecules for drug interactions that until now have been considered difficult or undruggable. It may be possible to deliver or express antibody fragments in target cells as macrodrugs per se. In future incarnations of intracellular antibodies, however, the structural information of the interaction interface of target and antibody fragment should facilitate development of binding site mimics as small drug‐like molecules. This is a new dawn for intracellular antibody fragments both as macrodrugs and as precursors of drugs to treat human diseases and should finally lead to the removal of the epithet of the ‘undruggable’ protein‐protein interactions.     

A humanized antibody inhibitor for cathepsin L

Cathepsin L (CTSL) is a cysteine protease involved in a variety of physiological and pathological processes. Potent inhibitors against CTSL have long been sought for drug development. Due to insufficient specificity and suboptimal pharmacological properties for current CTSL inhibitors, novel agents are still required for selectively blocking CTSL activity. Here we generated a humanized antibody inhibitor of CTSL by genetically fusing the inhibitory propeptide of procathepsin L to the N‐terminus of the light chain of a humanized antibody. The resulting antibody fusion could be stably expressed and displays highly potent inhibition activity and specificity toward CTSL. This work demonstrates a new approach for the rapid generation of antibody inhibitors of CTSL. It can possibly be extended to create inhibitory antibodies targeting other cathepsin proteases, providing novel research and therapeutic tools.     

Structural prediction of antibody‐APRIL complexes by computational docking constrained by antigen saturation mutagenesis library data

IgA nephropathy (IgAN) is the most prevalent cause of primary glomerular disease worldwide, and the cytokine A PRoliferation‐Inducing Ligand (APRIL) is emerging as a key player in IgAN pathogenesis and disease progression. For a panel of anti‐human APRIL antibodies with known antagonistic activity, we sought to define their structural mode of engagement to understand molecular mechanisms of action and aid rational antibody engineering. Reliable computational prediction of antibody‐antigen complexes remains challenging, and experimental methods such as X‐ray co‐crystallography and cryoEM have considerable technical, resource, and throughput barriers. To overcome these limitations, we implemented an integrated and accessible experimental‐computational workflow to more accurately predict structures of antibody‐APRIL complexes. Specifically, a yeast surface display library encoding site‐saturation mutagenized surface positions of APRIL was screened against a panel of anti‐APRIL antibodies to rapidly obtain a comprehensive biochemical profile of mutational impact on binding for each antibody. The experimentally derived mutational profile data were used as quantitative constraints in a computational docking workflow optimized for antibodies, resulting in robust structural models of antibody‐antigen complexes. The model results were confirmed by solving the cocrystal structure of one antibody‐APRIL complex, which revealed strong agreement with our model. The models were used to rationally select and engineer one antibody for cross‐species APRIL binding, which quite often aids further testing in relevant animal models. Collectively, we demonstrate a rapid, integrated computational‐experimental approach to robustly predict antibody‐antigen structures information, which can aid rational antibody engineering and provide insights into molecular mechanisms of action.     

Resolving self-association of a therapeutic antibody by formulation optimization and molecular approaches

A common challenge encountered during development of high concentration monoclonal antibody formulations is preventing self-association. Depending on the antibody and its formulation, self-association can be seen as aggregation, precipitation, opalescence or phase separation. Here we report on an unusual manifestation of self-association, formation of a semi-solid gel or “gelation." Therapeutic monoclonal antibody C4 was isolated from human B cells based on its strong potency in neutralizing bacterial toxin in animal models. The purified antibody possessed the unusual property of forming a firm, opaque white gel when it was formulated at concentrations >30 mg/mL and the temperature was <6°C. Gel formation was reversible with temperature. Gelation was affected by salt concentration or pH, suggesting an electrostatic interaction between IgG monomers. A comparison of the C4 amino acid sequences to consensus germline sequences revealed differences in framework regions. A C4 variant in which the framework sequence was restored to the consensus germline sequence did not gel at 100 mg/mL at temperatures as low as 1°C. Additional genetic analysis was used to predict the key residue(s) involved in the gelation. Strikingly, a single substitution in the native antibody, replacing heavy chain glutamate 23 with lysine (E23K), was sufficient to prevent gelation. These results indicate that the framework region is involved in intermolecular interactions. The temperature dependence of gelation may be related to conformational changes near glutamate 23 or the regions it interacts with. Molecular engineering of the framework can be an effective approach to resolve the solubility issues of therapeutic antibodies.     

Multi‐constraint computational design suggests that native sequences of germline antibody H3 loops are nearly optimal for conformational flexibility

The limited size of the germline antibody repertoire has to recognize a far larger number of potential antigens. The ability of a single antibody to bind multiple ligands due to conformational flexibility in the antigen‐binding site can significantly enlarge the repertoire. Among the six complementarity determining regions (CDRs) that generally comprise the binding site, the CDR H3 loop is particularly variable. Computational protein design studies showed that predicted low energy sequences compatible with a given backbone structure often have considerable similarity to the corresponding native sequences of naturally occurring proteins, indicating that native protein sequences are close to optimal for their structures. Here, we take a step forward to determine whether conformational flexibility, believed to play a key functional role in germline antibodies, is also central in shaping their native sequence. In particular, we use a multi‐constraint computational design strategy, along with the Rosetta scoring function, to propose that the native sequences of CDR H3 loops from germline antibodies are nearly optimal for conformational flexibility. Moreover, we find that antibody maturation may lead to sequences with a higher degree of optimization for a single conformation, while disfavoring sequences that are intrinsically flexible. In addition, this computational strategy allows us to predict mutations in the CDR H3 loop to stabilize the antigen‐bound conformation, a computational mimic of affinity maturation, that may increase antigen binding affinity by preorganizing the antigen binding loop. In vivo affinity maturation data are consistent with our predictions. The method described here can be useful to design antibodies with higher selectivity and affinity by reducing conformational diversity. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Toward high‐resolution homology modeling of antibody Fv regions and application to antibody–antigen docking

High‐resolution homology models are useful in structure‐based protein engineering applications, especially when a crystallographic structure is unavailable. Here, we report the development and implementation of RosettaAntibody, a protocol for homology modeling of antibody variable regions. The protocol combines comparative modeling of canonical complementarity determining region (CDR) loop conformations and de novo loop modeling of CDR H3 conformation with simultaneous optimization of V_L‐V_H rigid‐body orientation and CDR backbone and side‐chain conformations. The protocol was tested on a benchmark of 54 antibody crystal structures. The median root mean square deviation (rmsd) of the antigen binding pocket comprised of all the CDR residues was 1.5 Å with 80% of the targets having an rmsd lower than 2.0 Å. The median backbone heavy atom global rmsd of the CDR H3 loop prediction was 1.6, 1.9, 2.4, 3.1, and 6.0 Å for very short (4–6 residues), short (7–9), medium (10–11), long (12–14) and very long (17–22) loops, respectively. When the set of ten top‐scoring antibody homology models are used in local ensemble docking to antigen, a moderate‐to‐high accuracy docking prediction was achieved in seven of fifteen targets. This success in computational docking with high‐resolution homology models is encouraging, but challenges still remain in modeling antibody structures for sequences with long H3 loops. This first large‐scale antibody–antigen docking study using homology models reveals the level of “functional accuracy” of these structural models toward protein engineering applications. Proteins 2009; 74:497–514. © 2008 Wiley‐Liss, Inc.     

ScaffoldSeq: Software for characterization of directed evolution populations

ScaffoldSeq is software designed for the numerous applications—including directed evolution analysis—in which a user generates a population of DNA sequences encoding for partially diverse proteins with related functions and would like to characterize the single site and pairwise amino acid frequencies across the population. A common scenario for enzyme maturation, antibody screening, and alternative scaffold engineering involves naïve and evolved populations that contain diversified regions, varying in both sequence and length, within a conserved framework. Analyzing the diversified regions of such populations is facilitated by high‐throughput sequencing platforms; however, length variability within these regions (e.g., antibody CDRs) encumbers the alignment process. To overcome this challenge, the ScaffoldSeq algorithm takes advantage of conserved framework sequences to quickly identify diverse regions. Beyond this, unintended biases in sequence frequency are generated throughout the experimental workflow required to evolve and isolate clones of interest prior to DNA sequencing. ScaffoldSeq software uniquely handles this issue by providing tools to quantify and remove background sequences, cluster similar protein families, and dampen the impact of dominant clones. The software produces graphical and tabular summaries for each region of interest, allowing users to evaluate diversity in a site‐specific manner as well as identify epistatic pairwise interactions. The code and detailed information are freely available at http://research.cems.umn.edu/hackel . Proteins 2016; 84:869–874. © 2016 Wiley Periodicals, Inc.     

锌螯合物抗体可变区的克隆鉴定、真核表达及重组抗体的三维模建

筛选得到能分泌抗锌(Ⅱ)-二乙烯三胺五乙酸( Zn(H)-DTPA)的单克隆抗体杂交瘤细胞株2E8E3.经测定,此株细胞分泌IgM,轻链K型抗体.为了得到该抗体可变区序列,分别在小鼠抗体的信号肽区域和抗体的第一恒定区半胱氨酸残基(Cys)附近选择合适引物,通过RT-PCR方法,得到重(VH)、轻(V1)链可变区序列.将序列进行BLAST分析、SignalP分析、IMGT-V/QUEST分析,结果显示,克隆到的VH-2E8E3( 465 bp)和VL-2E8E3( 417 bp),序列与小鼠μ重链,κ轻链抗体可变区序列具有同源性(相似度＞97％),前端序列符合信号肽特征,抗体骨架区( Framework region,FR)及第一恒定区均含有保守半胱氨酸(Cys)残基.在此基础上,构建了重组抗体的真核表达载体,转染293T细胞并用间接ELISA法初步验证其具有表达活性.对重组抗体进行了三维模拟,并推测了主要氨基酸残基.试验结果表明,克隆到了正确的2E8E3抗体可变区序列并初步鉴定为有活性.     

An enriched structural kinase database to enable kinome‐wide structure‐based analyses and drug discovery

The development of a kinase structural database, the kinase knowledge base (KKB), is described. It covers all human kinase domain structures that have been deposited in the Protein Data Bank. All structures are renumbered using a common scheme, which enables efficient cross‐comparisons and multiple queries of interest to the kinase field. The common numbering scheme is also used to automatically annotate conserved residues and motifs, and conformationally classify the structures based on the DFG‐loop and Helix C. Analyses of residue conservation in the ATP binding site using the full human‐kinome–sequence alignment lead to the identification of a conserved hydrogen bond between the hinge region backbone and a glycine in the specificity surface. Furthermore, 90% of kinases are found to have at least one stabilizing interaction for the hinge region, which has not been described before.     

Kabat database and its applications: 30 years after the first variability plot

The Kabat Database was initially started in 1970 to determine the combining site of antibodies based on the available amino acid sequences at that time. Bence Jones proteins, mostly from human, were aligned, using the now-known Kabat numbering system, and a quantitative measure, variability, was calculated for every position. Three peaks, at positions 24-34, 50-56 and 89-97, were identified and proposed to form the complementarity determining regions (CDR) of light chains. Subsequently, antibody heavy chain amino acid sequences were also aligned using a different numbering system, since the locations of their CDRs (31-35B, 50-65 and 95-102) are different from those of the light chains. CDRL1 starts right after the first invariant Cys 23 of light chains, while CDRH1 is eight amino acid residues away from the first invariant Cys 22 of heavy chains. During the past 30 years, the Kabat database has grown to include nucleotide sequences, sequences of T cell receptors for antigens (TCR), major histocompatibility complex (MHC) class I and II molecules and other proteins of immunological interest. It has been used extensively by immunologists to derive useful structural and functional information from the primary sequences of these proteins. An overall view of the Kabat Database and its various applications are summarized here. The Kabat Database is freely available at http://immuno.bme.nwu.edu