"Redirecting an anti-IL-1ß antibody to bind a new, unrelated and computationally predicted epitope on hIL-17A".

Antibody engineering technology is at the forefront of therapeutic antibody development. The primary goal for engineering a therapeutic antibody is the generation of an antibody with a desired specificity, affinity, function, and developability profile. Mature antibodies are considered antigen specific, which may preclude their use as a starting point for antibody engineering. Here, we explore the plasticity of mature antibodies by engineering novel specificity and function to a pre-selected antibody template. Using a small, focused library, we engineered AAL160, an anti-IL-1ß antibody, to bind the unrelated antigen IL-17A, with the introduction of seven mutations. The final redesigned antibody, 11.003, retains favorable biophysical properties, binds IL-17A with sub-nanomolar affinity, inhibits IL-17A binding to its cognate receptor and is functional in a cell-based assay. The epitope of the engineered antibody can be computationally predicted based on the sequence of the template antibody, as is confirmed by the crystal structure of the 11.003/IL-17A complex. The structures of the 11.003/IL-17A and the AAL160/IL-1ß complexes highlight the contribution of germline residues to the paratopes of both the template and re-designed antibody. This case study suggests that the inherent plasticity of antibodies allows for re-engineering of mature antibodies to new targets, while maintaining desirable developability profiles.

Accurate profiling of full-length Fv in highly homologous antibody libraries using UMI tagged short reads.

Deep parallel sequencing (NGS) is a viable tool for monitoring scFv and Fab library dynamics in many antibody engineering high-throughput screening efforts. Although very useful, the commonly used Illumina NGS platform cannot handle the entire sequence of scFv or Fab in a single read, usually focusing on specific CDRs or resorting to sequencing VH and VL variable domains separately, thus limiting its utility in comprehensive monitoring of selection dynamics. Here we present a simple and robust method for deep sequencing repertoires of full length scFv, Fab and Fv antibody sequences. This process utilizes standard molecular procedures and unique molecular identifiers (UMI) to pair separately sequenced VH and VL. We show that UMI assisted VH-VL matching allows for a comprehensive and highly accurate mapping of full length Fv clonal dynamics in large highly homologous antibody libraries, as well as identification of rare variants. In addition to its utility in synthetic antibody discovery processes, our method can be instrumental in generating large datasets for machine learning (ML) applications, which in the field of antibody engineering has been hampered by conspicuous paucity of large scale full length Fv data.

Computational Design of Epitope-Specific Functional Antibodies.

The ultimate goal of protein design is to introduce new biological activity. We propose a computational approach for designing functional antibodies by focusing on functional epitopes, integrating large-scale statistical analysis with multiple structural models. Machine learning is used to analyze these models and predict specific residue-residue contacts. We use this approach to design a functional antibody to counter the proinflammatory effect of the cytokine interleukin-17A (IL-17A). X-ray crystallography confirms that the designed antibody binds the targeted epitope and the interaction is mediated by the designed contacts. Cell-based assays confirm that the antibody is functional. Importantly, this approach does not rely on a high-quality 3D model of the designed complex or even a solved structure of the target. As demonstrated here, this approach can be used to design biologically active antibodies, removing some of the main hurdles in antibody design and in drug discovery.

Computational design of antibodies.

Antibody design aims to create new antibodies with biological activity that can be used in therapy and research. Traditional methods for antibody discovery, such as animal immunization and large-scale library screening, generate antibodies that bind to the target of interest, but do not necessarily have the desired functional effect. Computational methods can be utilized as a means to guide the search for biologically relevant antibodies, focusing on specificity and affinity determinants to target a particular region of the antigen. Such an approach would allow for the design of epitope-specific antibodies that will have the desired effect on the function of the targeted protein.

Small molecule correctors of F508del-CFTR discovered by structure-based virtual screening.

Folding correctors of F508del-CFTR were discovered by in silico structure-based screening utilizing homology models of CFTR. The intracellular segment of CFTR was modeled and three cavities were identified at inter-domain interfaces: (1) Interface between the two Nucleotide Binding Domains (NBDs); (2) Interface between NBD1 and Intracellular Loop (ICL) 4, in the region of the F508 deletion; (3) multi-domain interface between NBD1:2:ICL1:2:4. We hypothesized that compounds binding at these interfaces may improve the stability of the protein, potentially affecting the folding yield or surface stability. In silico structure-based screening was performed at the putative binding-sites and a total of 496 candidate compounds from all three sites were tested in functional assays. A total of 15 compounds, representing diverse chemotypes, were identified as F508del folding correctors. This corresponds to a 3% hit rate, ~tenfold higher than hit rates obtained in corresponding high-throughput screening campaigns. The same binding sites also yielded potentiators and, most notably, compounds with a dual corrector-potentiator activity (dual-acting). Compounds harboring both activity types may prove to be better leads for the development of CF therapeutics than either pure correctors or pure potentiators. To the best of our knowledge this is the first report of structure-based discovery of CFTR modulators.

The cystic fibrosis-causing mutation deltaF508 affects multiple steps in cystic fibrosis transmembrane conductance regulator biogenesis.

The deletion of phenylalanine 508 in the first nucleotide binding domain of the cystic fibrosis transmembrane conductance regulator is directly associated with >90% of cystic fibrosis cases. This mutant protein fails to traffic out of the endoplasmic reticulum and is subsequently degraded by the proteasome. The effects of this mutation may be partially reversed by the application of exogenous osmolytes, expression at low temperature, and the introduction of second site suppressor mutations. However, the specific steps of folding and assembly of full-length cystic fibrosis transmembrane conductance regulator (CFTR) directly altered by the disease-causing mutation are unclear. To elucidate the effects of the ΔF508 mutation, on various steps in CFTR folding, a series of misfolding and suppressor mutations in the nucleotide binding and transmembrane domains were evaluated for effects on the folding and maturation of the protein. The results indicate that the isolated NBD1 responds to both the ΔF508 mutation and intradomain suppressors of this mutation. In addition, identification of a novel second site suppressor of the defect within the second transmembrane domain suggests that ΔF508 also effects interdomain interactions critical for later steps in the biosynthesis of CFTR.

Solution structure of Archaeglobus fulgidis peptidyl-tRNA hydrolase (Pth2) provides evidence for an extensive conserved family of Pth2 enzymes in archea, bacteria, and eukaryotes.

The solution structure of protein AF2095 from the thermophilic archaea Archaeglobus fulgidis, a 123-residue (13.6-kDa) protein, has been determined by NMR methods. The structure of AF2095 is comprised of four alpha-helices and a mixed beta-sheet consisting of four parallel and anti-parallel beta-strands, where the alpha-helices sandwich the beta-sheet. Sequence and structural comparison of AF2095 with proteins from Homo sapiens, Methanocaldococcus jannaschii, and Sulfolobus solfataricus reveals that AF2095 is a peptidyl-tRNA hydrolase (Pth2). This structural comparison also identifies putative catalytic residues and a tRNA interaction region for AF2095. The structure of AF2095 is also similar to the structure of protein TA0108 from archaea Thermoplasma acidophilum, which is deposited in the Protein Data Bank but not functionally annotated. The NMR structure of AF2095 has been further leveraged to obtain good-quality structural models for 55 other proteins. Although earlier studies have proposed that the Pth2 protein family is restricted to archeal and eukaryotic organisms, the similarity of the AF2095 structure to human Pth2, the conservation of key active-site residues, and the good quality of the resulting homology models demonstrate a large family of homologous Pth2 proteins that are conserved in eukaryotic, archaeal, and bacterial organisms, providing novel insights in the evolution of the Pth and Pth2 enzyme families.

The SufE sulfur-acceptor protein contains a conserved core structure that mediates interdomain interactions in a variety of redox protein complexes.

The isc and suf operons in Escherichia coli represent alternative genetic systems optimized to mediate the essential metabolic process of iron-sulfur cluster (Fe-S) assembly under basal or oxidative-stress conditions, respectively. Some of the proteins in these two operons share strong sequence homology, e.g. the cysteine desulfurases IscS and SufS, and presumably play the same role in the oxygen-sensitive assembly process. However, other proteins in these operons share no significant homology and occur in a mutually exclusive manner in Fe-S assembly operons in other organisms (e.g. IscU and SufE). These latter proteins presumably play distinct roles adapted to the different assembly mechanisms used by the two systems. IscU has three invariant cysteine residues that function as a template for Fe-S assembly while accepting a sulfur atom from IscS. SufE, in contrast, does not function as an Fe-S assembly template but has been suggested to function as a shuttle protein that uses a persulfide linkage to a single invariant cysteine residue to transfer a sulfur atom from SufS to an alternative Fe-S assembly template. Here, we present and analyze the 2.0A crystal structure of E.coli SufE. The structure shows that the persulfide-forming cysteine occurs at the tip of a loop with elevated B-factors, where its side-chain is buried from solvent exposure in a hydrophobic cavity located beneath a highly conserved surface. Despite the lack of sequence homology, the core of SufE shows strong structural similarity to IscU, and the sulfur-acceptor site in SufE coincides with the location of the cysteine residues mediating Fe-S cluster assembly in IscU. Thus, a conserved core structure is implicated in mediating the interactions of both SufE and IscU with the mutually homologous cysteine desulfurase enzymes present in their respective operons. A similar core structure is observed in a domain found in a variety of Fe-S cluster containing flavoenzymes including xanthine dehydrogenase, where it also mediates interdomain interactions. Therefore, the core fold of SufE/IscU has been adapted to mediate interdomain interactions in diverse redox protein systems in the course of evolution.

Solution NMR structure of the iron-sulfur cluster assembly protein U (IscU) with zinc bound at the active site.

IscU is a highly conserved protein that serves as the scaffold for IscS-mediated assembly of iron-sulfur ([Fe-S]) clusters. We report the NMR solution structure of monomeric Haemophilus influenzae IscU with zinc bound at the [Fe-S] cluster assembly site. The compact core of the globular structure has an alpha-beta sandwich architecture with a three-stranded antiparallel beta-sheet and four alpha-helices. A nascent helix is located N-terminal to the core structure. The zinc is ligated by three cysteine residues and one histidine residue that are located in and near conformationally dynamic loops at one end of the IscU structure. Removal of the zinc metal by chelation results in widespread loss of structure in the apo form. The zinc-bound IscU may be a good model for iron-loaded IscU and may demonstrate structural features found in the [Fe-S] cluster bound form. Structural and functional similarities, genomic context in operons containing other homologous genes, and distributions of conserved surface residues support the hypothesis that IscU protein domains are homologous (i.e. derived from a common ancestor) with the SufE/YgdK family of [Fe-S] cluster assembly proteins.

Solution NMR structure of the 30S ribosomal protein S28E from Pyrococcus horikoshii.

We report NMR assignments and solution structure of the 71-residue 30S ribosomal protein S28E from the archaean Pyrococcus horikoshii, target JR19 of the Northeast Structural Genomics Consortium. The structure, determined rapidly with the aid of automated backbone resonance assignment (AutoAssign) and automated structure determination (AutoStructure) software, is characterized by a four-stranded beta-sheet with a classic Greek-key topology and an oligonucleotide/oligosaccharide beta-barrel (OB) fold. The electrostatic surface of S28E exhibits positive and negative patches on opposite sides, the former constituting a putative binding site for RNA. The 13 C-terminal residues of the protein contain a consensus sequence motif constituting the signature of the S28E protein family. Surprisingly, this C-terminal segment is unstructured in solution.

Using multiple structure alignments, fast model building, and energetic analysis in fold recognition and homology modeling.

We participated in the fold recognition and homology sections of CASP5 using primarily in-house software. The central feature of our structure prediction strategy involved the ability to generate good sequence-to-structure alignments and to quickly transform them into models that could be evaluated both with energy-based methods and manually. The in-house tools we used include: a) HMAP (Hybrid Multidimensional Alignment Profile)-a profile-to-profile alignment method that is derived from sequence-enhanced multiple structure alignments in core regions, and sequence motifs in non-structurally conserved regions. b) NEST-a fast model building program that applies an "artificial evolution" algorithm to construct a model from a given template and alignment. c) GRASP2-a new structure and alignment visualization program incorporating multiple structure superposition and domain database scanning modules. These methods were combined with model evaluation based on all atom and simplified physical-chemical energy functions. All of these methods were under development during CASP5 and consequently a great deal of manual analysis was carried out at each stage of the prediction process. This interactive model building procedure has several advantages and suggests important ways in which our and other methods can be improved, examples of which are provided.

Structural genomics: computational methods for structure analysis.

The success of structural genomics initiatives requires the development and application of tools for structure analysis, prediction, and annotation. In this paper we review recent developments in these areas; specifically structure alignment, the detection of remote homologs and analogs, homology modeling and the use of structures to predict function. We also discuss various rationales for structural genomics initiatives. These include the structure-based clustering of sequence space and genome-wide function assignment. It is also argued that structural genomics can be integrated into more traditional biological research if specific biological questions are included in target selection strategies.

Solution structure of Vibrio cholerae protein VC0424: a variation of the ferredoxin-like fold.

The structure of Vibrio cholerae protein VC0424 was determined by NMR spectroscopy. VC0424 belongs to a conserved family of bacterial proteins of unknown function (COG 3076). The structure has an alpha-beta sandwich architecture consisting of two layers: a four-stranded antiparallel beta-sheet and three side-by-side alpha-helices. The secondary structure elements have the order alphabetaalphabetabetaalphabeta along the sequence. This fold is the same as the ferredoxin-like fold, except with an additional long N-terminal helix, making it a variation on this common motif. A cluster of conserved surface residues on the beta-sheet side of the protein forms a pocket that may be important for the biological function of this conserved family of proteins.