High affinity extremes in combinatorial libraries and repertoires

By generating a large diversity of molecules, the immune system selects antibodies that bind antigens. Sharing the same approach, combinatorial biotechnologies use a large library of compounds to screen for molecules of high affinity to a given target. Understanding the properties of the best binders in the pool aids the design of the library. In particular, how does the maximum affinity increase with the size of the library or repertoire? We consider two alternative models to examine the properties of extreme affinities. In the first model, affinities are distributed lognormally, while in the second, affinities are determined by the number of matches to a target sequence. The second model more explicitly models nucleic acids (DNA or RNA) and proteins such as antibodies. Using extreme value theory we show that the logarithm of the mean of the highest affinity in a combinatorial library grows linearly with the square root of the log of the library size. When there is an upper bound to affinity, this “absolute maximum” is also approached approximately linearly with root log library size, reaching the upper limit abruptly. The design of libraries may benefit from considering how this plateau is reached as the library size is increased.     

Proteome sampling by the HLA class I antigen processing pathway

The peptide repertoire that is presented by the set of HLA class I molecules of an individual is formed by the different players of the antigen processing pathway and the stringent binding environment of the HLA class I molecules. Peptide elution studies have shown that only a subset of the human proteome is sampled by the antigen processing machinery and represented on the cell surface. In our study, we quantified the role of each factor relevant in shaping the HLA class I peptide repertoire by combining peptide elution data, in silico predictions of antigen processing and presentation, and data on gene expression and protein abundance. Our results indicate that gene expression level, protein abundance, and rate of potential binding peptides per protein have a clear impact on sampling probability. Furthermore, once a protein is available for the antigen processing machinery in sufficient amounts, C-terminal processing efficiency and binding affinity to the HLA class I molecule determine the identity of the presented peptides. Having studied the impact of each of these factors separately, we subsequently combined all factors in a logistic regression model in order to quantify their relative impact. This model demonstrated the superiority of protein abundance over gene expression level in predicting sampling probability. Being able to discriminate between sampled and non-sampled proteins to a significant degree, our approach can potentially be used to predict the sampling probability of self proteins and of pathogen-derived proteins, which is of importance for the identification of autoimmune antigens and vaccination targets.     

Over-expression of RAD51 or RAD54 but not RAD51/4 enhances extra-chromosomal homologous recombination in the human sarcoma (HT-1080) cell line

RAD51 and RAD54, members of the RAD52 epistasis group, play key roles in homologous recombination (HR). The efficiency of homologous recombination (HR) can be increased by over-expression of either of them. A vector that allows co-expression of RAD51 and RAD54 was constructed to investigate interactions between the two proteins during extra-chromosomal HR. The efficiency of extra-chromosomal HR evaluated by GFP extra-chromosomal HR was enhanced (110-245%) in different transfected Human sarcoma (HT-1080) cell colonies. We observed that RAD51 clearly promotes extra-chromosomal HR; however, the actions of RAD54 in extra-chromosomal HR were weak. Our data suggest that RAD51 may function as a universal factor during HR, whereas RAD54 mainly functions in other types of HR (gene targeting or intra-chromosomal HR), which involves interaction with chromosomal DNA.     

Effect of decreasing the affinity of the class II-associated invariant chain peptide on the MHC class II peptide repertoire in the presence or absence of H-2M

The class II-associated invariant chain peptide (CLIP) region of the invariant chain (Ii) directly influences MHC class II presentation by occupying the MHC class II peptide-binding groove, thereby preventing premature loading of peptides. Different MHC class II alleles exhibit distinct affinities for CLIP, and a low affinity interaction has been associated with decreased dependence upon H-2M and increased susceptibility to rheumatoid arthritis, suggesting that decreased CLIP affinity alters the MHC class II-bound peptide repertoire, thereby promoting autoimmunity. To examine the role of CLIP affinity in determining the MHC class II peptide repertoire, we generated transgenic mice expressing either wild-type human Ii or human Ii containing a CLIP region of low affinity for MHC class II. Our data indicate that although degradation intermediates of Ii containing a CLIP region with decreased affinity for MHC class II do not remain associated with I-A(b), this does not substantially alter the peptide repertoire bound by MHC class II or increase autoimmune susceptibility in the mice. This implies that the affinity of the CLIP:MHC class II interaction is not a strong contributory factor in determining the probability of developing autoimmunity. In contrast, in the absence of H-2M, MHC class II peptide repertoire diversity is enhanced by decreasing the affinity of CLIP for MHC class II, although MHC class II cell surface expression is reduced. Thus, we show clearly, in vivo, the critical chaperone function of H-2M, which preserves MHC class II molecules for high affinity peptide binding upon dissociation of Ii degradation intermediates.     

Isolation of multiple cell-binding ligands from different phage displayed-peptide libraries

A technical challenge in the development of biosensor devices for cancer detection and diagnosis is the identification of ligands that recognize cancer cells with high affinity and specificity. Furthermore, it is unlikely that one cell-binding ligand will provide sufficient biological information, thus, multiple ligands for a given cancer type will be needed for confident clinical diagnosis. Biopanning of phage displayed peptide libraries is a route to isolation of specific cell-binding reagents. A potential approach towards isolation of multiple ligands for a single cell type is to pan against the same cell type using different peptide libraries. Here we report the synthesis of a new 20-mer peptide-phage library and its use to select a peptide that binds to the large cell lung carcinoma cell line, H1299. The isolated phage clone binds H1299 cells 80 times better than a control phage and can distinguish between H1299 and normal control cells. The phage clone also binds to the lung pleura epidermoid cell line, Calu-1 but not to all lung carcinoma cell lines. The peptide is functional outside the context of the phage and tetramerization of the peptide on a trilysine core improves the affinity of the peptide. The tetrameric peptide can be used to deliver a fluorescent quantum dot to H1299 cells. Unexpectedly, the peptide shares sequence similarity to a previously isolated H1299-binding peptide isolated from a different 20-mer peptide library. Data suggests that the two peptides target the same cellular receptor. Our results imply that cell-based biopanning can isolate cell-binding ligands that may be of utility for cancer diagnosis, and isolation of cell-targeting peptides from different peptide libraries can expand the repertoire of cell-binding reagents.     

秩——丰盛度分布是热带蛾类集合判别的有用工具吗?

制图并分析秩——丰盛度分布(RADs)是描述物种集合特征的有效方法。研究一直认为,来自不同生境的RADs秩——丰盛度分布倾向于符合不同的模型,如对数序列分布或对数正态分布。我们用从婆罗洲北部收集的63套尺蠖蛾样本数据,给这个想法一个严格的检验。通过对数序列分布和对数正态分布间的二项分布比较及外加另外三个生态位-空间分隔分布的模型间多项分布比较,发现生境干扰或海拔均无法显著地预测最适合的RAD模型。另一方面,样本多样性和在较小程度上的样本大小(个体数)与RAD模型显著相关。我们的结论是,尽管RAD模型拟合有许多其它的有趣应用,但在根据生境类型来判别物种集合方面是无效的。     

Identification of novel isoforms of human RAD52

In yeast, RAD52 has been shown to be essential for homologous recombination of DNA and to be involved in the repair of double-stranded DNA breaks. Recently, the human homologue of yeast RAD52, a 418-amino-acid protein, has been identified. In this study, we report three different isoforms of human RAD52 isolated from brain and testis cDNA libraries. cDNAs of these isoforms contain distinct insertions and encode truncated proteins due to translational frame-shifts. The three isoforms consist of 226-, 139-, and 118-amino-acid residues, and are designated as RAD52beta, gamma, and delta, respectively. The original RAD52 is termed as RAD52alpha in this paper. Messages of these isoforms have been detected in various human tissues. We found that the RAD52 isoforms were unable to interact with RAD52alpha because of partial defect of the self-interaction domain. Furthermore, like RAD52alpha, the isoforms have been shown to bind to both single-stranded and double-stranded DNA. These results suggest that RAD52beta, gamma, and delta might affect RAD52alpha function through their DNA-binding property and their inability to bind to RAD52alpha. Thus, these isoforms might act as dominant negative mutants or negative regulators of RAD52alpha.     

RAD51AP1 mediates RAD51 activity through nucleosome interaction

RAD51-associated protein 1 (RAD51AP1) is a key protein in the homologous recombination (HR) DNA repair pathway. Loss of RAD51AP1 leads to defective HR, genome instability, and telomere erosion. RAD51AP1 physically interacts with the RAD51 recombinase and promotes RAD51-mediated capture of donor DNA, synaptic complex assembly, and displacement-loop formation when tested with nucleosome-free DNA substrates. In cells, however, DNA is packaged into chromatin, posing an additional barrier to the complexities of the HR reaction. In this study, we show that RAD51AP1 binds to nucleosome core particles (NCPs), the minimum basic unit of chromatin in which approximately two superhelical turns of 147 bp double-stranded DNA are wrapped around one histone octamer with no free DNA ends remaining. We identified a C-terminal region in RAD51AP1, including its previously mapped DNA-binding domain, as critical for mediating the association between RAD51AP1 and both the NCP and the histone octamer. Using in vitro surrogate assays of HR activity, we show that RAD51AP1 is capable of promoting duplex DNA capture and initiating joint-molecule formation with the NCP and chromatinized template DNA, respectively. Together, our results suggest that RAD51AP1 directly assists in the RAD51-mediated search for donor DNA in chromatin. We present a model, in which RAD51AP1 anchors the DNA template through affinity for its nucleosomes to the RAD51-ssDNA nucleoprotein filament.     

Telomere biology: integrating chromosomal end protection with DNA damage response

Telomeres play the key protective role at chromosomes. Many studies indicate that loss of telomere function causes activation of DNA damage response. Here, we review evidence supporting interdependence between telomere maintenance and DNA damage response and present a model in which these two pathways are combined into a single mechanism for protecting chromosomal integrity. Proteins directly involved in telomere maintenance and DNA damage response include Ku, DNA-PKcs, RAD51D, PARP-2, WRN and RAD50/MRE11/NBS1 complex. Since most of these proteins participate in the repair of DNA double-strand breaks (DSBs), this was perceived by many authors as a paradox, given that telomeres function to conceal natural DNA ends from mechanisms that detect and repair DSBs. However, we argue here that the key function of one particular DSB protein, Ku, is to prevent or control access of telomerase, the enzyme that synthesises telomeric sequences, to both internal DSBs and natural chromosomal ends. This view is supported by observations that Ku has a high affinity for DNA ends; it acts as a negative regulator of telomerase and that telomerase itself can target internal DSBs. Ku then directs other DSB repair/telomere maintenance proteins to either repair DSBs at internal chromosomal sites or prevent uncontrolled elongation of telomeres by telomerase. This model eliminates the above paradox and provides a testable scenario in which the role of DSB repair proteins is to protect chromosomal integrity by balancing repair activities and telomere maintenance. In our model, a close association between telomeres and different DNA damage response factors is not an unexpected event, but rather a logical result of chromosomal integrity maintenance activities. Review related to the 15th International Chromosome Conference (ICC XV), held in September 2004, Brunel University, London, UK     

Phage display: Concept,innovations, applications and future

Phage display is the technology that allows expression of exogenous (poly)peptides on the surface of phage particles. The concept is simple in principle: a library of phage particles expressing a wide diversity of peptides is used to select those that bind the desired target. The filamentous phage M13 is the most commonly used vector to create random peptide display libraries. Several methods including recombinant techniques have been developed to increase the diversity of the library. On the other extreme, libraries with various biases can be created for specific purposes. For instance, when the sequence of the peptide that binds the target is known, its affinity and selectivity can be increased by screening libraries created with limited mutagenesis of the peptide. Phage libraries are screened for binding to synthetic or native targets. The initial screening of library by basic biopanning has been extended to column chromatography including negative screening and competition between selected phage clones to identify high affinity ligands with greater target specificity. The rapid isolation of specific ligands by phage display is advantageous in many applications including selection of inhibitors for the active and allosteric sites of the enzymes, receptor agonists and antagonists, and G-protein binding modulatory peptides. Phage display has been used in epitope mapping and analysis of protein-protein interactions. The specific ligands isolated from phage libraries can be used in therapeutic target validation, drug design and vaccine development. Phage display can also be used in conjunction with other methods. The past innovations and those to come promise a bright future for this field.     

Synthesis and Screening of an Indexed Motif-Library Containing Non-proteinogenic Amino Acids

In an effort to increase the probability of finding novel peptides in resin-bound combinatorial libraries displaying affinity to various macromolecular targets, we increased the diversity of a solid-phase library considerably by synthesizing multiple structures on each bead – a motif-library – including 45 building blocks. The building blocks consist of L -aa, D -aa and eight hydrophobic non-proteinogenic α-amino acids. A library with the format O-Z_0–1-O-Z_0–1-O-XX-resin was synthesized giving the four motifs OOOXX, OZOOXX, OOZOXX, OZOZOXX corresponding to 364.500 different motifs (45³×4 theoretical combinations). The positions O are defined amino acids while Z represents three mixtures Π, Ω, ϖ, where Π is a mixture of polar and charged residues, Ω is a mixture of aliphatic residues and ϖ is a mixture of aromatic residues. X represents a mixture of all 45 residues. The library was screened with the macromolecular target streptavidin which served as a model receptor. Binding peptides were sequenced by microsequencing. We included small amounts of norvaline and norleucine in the library, which served as index residues to be able to distinguish between LD -amino acids and other residues with the same retention time in the HPLC system. Beads that interact with the receptor were found, and the binding motifs that appeared had no homology to known binding motifs found in either L -aa or D -aa libraries, instead motifs with the non-proteinogenic residues L -phenylglycine, O-benzyl-L -hydroxyproline and O-benzyl-L -tyrosine dominated. The novel peptides inhibit binding of biotin to streptavidin but do not bind to avidin, and the affinity is higher than the peptides found in linear all L -aa peptide libraries. © 1997 European Peptide Society and John Wiley & Sons, Ltd     

Parameters and determinants of responses to selection in antibody libraries

The sequences of antibodies from a given repertoire are highly diverse at few sites located on the surface of a genome-encoded larger scaffold. The scaffold is often considered to play a lesser role than highly diverse, non-genome-encoded sites in controlling binding affinity and specificity. To gauge the impact of the scaffold, we carried out quantitative phage display experiments where we compare the response to selection for binding to four different targets of three different antibody libraries based on distinct scaffolds but harboring the same diversity at randomized sites. We first show that the response to selection of an antibody library may be captured by two measurable parameters. Second, we provide evidence that one of these parameters is determined by the degree of affinity maturation of the scaffold, affinity maturation being the process by which antibodies accumulate somatic mutations to evolve towards higher affinities during the natural immune response. In all cases, we find that libraries of antibodies built around maturated scaffolds have a lower response to selection to other arbitrary targets than libraries built around germline-based scaffolds. We thus propose that germline-encoded scaffolds have a higher selective potential than maturated ones as a consequence of a selection for this potential over the long-term evolution of germline antibody genes. Our results are a first step towards quantifying the evolutionary potential of biomolecules.     

Thermodynamic additivity of sequence variations: an algorithm for creating high affinity peptides without large libraries or structural information

Background

There is a significant need for affinity reagents with high target affinity/specificity that can be developed rapidly and inexpensively. Existing affinity reagent development approaches, including protein mutagenesis, directed evolution, and fragment-based design utilize large libraries and/or require structural information thereby adding time and expense. Until now, no systematic approach to affinity reagent development existed that could produce nanomolar affinity from small chemically synthesized peptide libraries without the aid of structural information.

Methodology/Principal Findings

Based on the principle of additivity, we have developed an algorithm for generating high affinity peptide ligands. In this algorithm, point-variations in a lead sequence are screened and combined in a systematic manner to achieve additive binding energies. To demonstrate this approach, low-affinity lead peptides for multiple protein targets were identified from sparse random sequence space and optimized to high affinity in just two chemical steps. In one example, a TNF-α binding peptide with K_d = 90 nM and high target specificity was generated. The changes in binding energy associated with each variation were generally additive upon combining variations, validating the basis of the algorithm. Interestingly, cooperativity between point-variations was not observed, and in a few specific cases, combinations were less than energetically additive.

Conclusions/Significance

By using this additivity algorithm, peptide ligands with high affinity for protein targets were generated. With this algorithm, one of the highest affinity TNF-α binding peptides reported to date was produced. Most importantly, high affinity was achieved from small, chemically-synthesized libraries without the need for structural information at any time during the process. This is significantly different than protein mutagenesis, directed evolution, or fragment-based design approaches, which rely on large libraries and/or structural guidance. With this algorithm, high affinity/specificity peptide ligands can be developed rapidly, inexpensively, and in an entirely chemical manner.     

Statistical analysis of antigen receptor spectratype data

MOTIVATION: The effectiveness of vertebrate adaptive immunity depends crucially on the establishment and maintenance of extreme diversity in the antigen receptor repertoire. Spectratype analysis is a method used in clinical and basic immunological settings in which antigen receptor length diversity is assessed as a surrogate for functional diversity. The purpose of this paper is to describe the systematic derivation and application of statistical methods for the analysis of spectratype data. RESULTS: The basic probability model used for spectratype analysis is the multinomial model with n, the total number of counts, indeterminate. We derive the appropriate statistics and statistical procedures for testing hypotheses regarding differences in antigen receptor distributions and variable repertoire diversity in different treatment groups.We then apply these methods to spectratype data obtained from several healthy donors to examine the differences between normal CD4+ and CD8+ T cell repertoires, and to data from a thymus transplant patient to examine the development of repertoire diversity following the transplant.     

RAD18 mediates resistance to ionizing radiation in human glioma cells

Radioresistance remains a major challenge in the treatment of glioblastoma multiforme (GBM). RAD18 a central regulator of translesion DNA synthesis (TLS), has been shown to play an important role in regulating genomic stability and DNA damage response. In the present study, we investigate the relationship between RAD18 and resistance to ionizing radiation (IR) and examined the expression levels of RAD18 in primary and recurrent GBM specimens. Our results showed that RAD18 is an important mediator of the IR-induced resistance in GBM. The expression level of RAD18 in glioma cells correlates with their resistance to IR. Ectopic expression of RAD18 in RAD18-low A172 glioma cells confers significant resistance to IR treatment. Conversely, depletion of endogenous RAD18 in RAD18-high glioma cells sensitized these cells to IR treatment. Moreover, RAD18 overexpression confers resistance to IR-mediated apoptosis in RAD18-low A172 glioma cells, whereas cells deficient in RAD18 exhibit increased apoptosis induced by IR. Furthermore, knockdown of RAD18 in RAD18-high glioma cells disrupts HR-mediated repair, resulting in increased accumulation of DSB. In addition, clinical data indicated that RAD18 was significantly higher in recurrent GBM samples that were exposed to IR compared with the corresponding primary GBM samples. Collectively, our findings reveal that RAD18 may serve as a key mediator of the IR response and may function as a potential target for circumventing IR resistance in human GBM.     

The RAD52 recombinational repair pathway is essential in pol30 (PCNA) mutants that accumulate small single-stranded DNA fragments during DNA synthesis.

To identify in vivo pathways that compensate for impaired proliferating cell nuclear antigen (PCNA or Pol30p in yeast) activity, we performed a synthetic lethal screen with the yeast pol30-104 mutation. We identified nine mutations that display synthetic lethality with pol30-104; three mutations affected the structural gene for the large subunit of replication factor C (rfc1), which loads PCNA onto DNA, and six mutations affected three members of the RAD52 epistasis group for DNA recombinational repair (rad50, rad52 and rad57). We also found that pol30-104 displayed synthetic lethality with mutations in other members of the RAD52 epistasis group (rad51 and rad54), but not with mutations in members of the RAD3 nor the RAD6 epistasis group. Analysis of nine different pol30 mutations shows that the requirement for the RAD52 pathway is correlated with a DNA replication defect but not with the relative DNA repair defect caused by pol30 mutations. In addition, mutants that require RAD52 for viability (pol30-100, pol30-104, rfc1-1 and rth1delta) accumulate small single-stranded DNA fragments during DNA replication in vivo. Taken together, these data suggest that the RAD52 pathway is required when there are defects in the maturation of Okazaki fragments.     

Assessment of flexible backbone protein design methods for sequence library prediction in the therapeutic antibody Herceptin-HER2 interface

Computational protein design methods can complement experimental screening and selection techniques by predicting libraries of low-energy sequences compatible with a desired structure and function. Incorporating backbone flexibility in computational design allows conformational adjustments that should broaden the range of predicted low-energy sequences. Here, we evaluate computational predictions of sequence libraries from different protocols for modeling backbone flexibility using the complex between the therapeutic antibody Herceptin and its target human epidermal growth factor receptor 2 (HER2) as a model system. Within the program RosettaDesign, three methods are compared: The first two use ensembles of structures generated by Monte Carlo protocols for near-native conformational sampling: kinematic closure (KIC) and backrub, and the third method uses snapshots from molecular dynamics (MD) simulations. KIC or backrub methods were better able to identify the amino acid residues experimentally observed by phage display in the Herceptin-HER2 interface than MD snapshots, which generated much larger conformational and sequence diversity. KIC and backrub, as well as fixed backbone simulations, captured the key mutation Asp98Trp in Herceptin, which leads to a further threefold affinity improvement of the already subnanomolar parental Herceptin-HER2 interface. Modeling subtle backbone conformational changes may assist in the design of sequence libraries for improving the affinity of antibody-antigen interfaces and could be suitable for other protein complexes for which structural information is available.     

Interactions among Trypanosoma brucei RAD51 paralogues in DNA repair and antigenic variation

Homologous recombination in Trypanosoma brucei is used for moving variant surface glycoprotein (VSG) genes into expression sites during immune evasion by antigenic variation. A major route for such VSG switching is gene conversion reactions in which RAD51, a universally conserved recombinase, catalyses homology-directed strand exchange. In any eukaryote, RAD51-directed strand exchange in vivo is mediated by further factors, including RAD51-related proteins termed Rad51 paralogues. These appear to be ubiquitously conserved, although their detailed roles in recombination remain unclear. In T. brucei, four putative RAD51 paralogue genes have been identified by sequence homology. Here we show that all four RAD51 paralogues act in DNA repair, recombination and RAD51 subnuclear dynamics, though not equivalently, while mutation of only one RAD51 paralogue gene significantly impedes VSG switching. We also show that the T. brucei RAD51 paralogues interact, and that the complexes they form may explain the distinct phenotypes of the mutants as well as observed expression interdependency. Finally, we document the Rad51 paralogues that are encoded by a wide range of protists, demonstrating that the Rad51 paralogue repertoire in T. brucei is unusually large among microbial eukaryotes and that one member of the protein family corresponds with a key, conserved eukaryotic Rad51 paralogue.     

Solid-state, dye-labeled DNA detects volatile compounds in the vapor phase

This paper demonstrates a previously unreported property of deoxyribonucleic acid—the ability of dye-labeled, solid-state DNA dried onto a surface to detect odors delivered in the vapor phase by changes in fluorescence. This property is useful for engineering systems to detect volatiles and provides a way for artificial sensors to emulate the way cross-reactive olfactory receptors respond to and encode single odorous compounds and mixtures. Recent studies show that the vertebrate olfactory receptor repertoire arises from an unusually large gene family and that the receptor types that have been tested so far show variable breadths of response. In designing biomimetic artificial noses, the challenge has been to generate a similarly large sensor repertoire that can be manufactured with exact chemical precision and reproducibility and that has the requisite combinatorial complexity to detect odors in the real world. Here we describe an approach for generating and screening large, diverse libraries of defined sensors using single-stranded, fluorescent dye–labeled DNA that has been dried onto a substrate and pulsed with brief exposures to different odors. These new solid-state DNA-based sensors are sensitive and show differential, sequence-dependent responses. Furthermore, we show that large DNA-based sensor libraries can be rapidly screened for odor response diversity using standard high-throughput microarray methods. These observations describe new properties of DNA and provide a generalized approach for producing explicitly tailored sensor arrays that can be rationally chosen for the detection of target volatiles with different chemical structures that include biologically derived odors, toxic chemicals, and explosives.