A new clustering of antibody CDR loop conformations

Previous analyses of the complementarity-determining regions (CDRs) of antibodies have focused on a small number of “canonical” conformations for each loop. This is primarily the result of the work of Chothia and coworkers, most recently in 1997. Because of the widespread utility of antibodies, we have revisited the clustering of conformations of the six CDR loops with the much larger amount of structural information currently available. In this work, we were careful to use a high-quality data set by eliminating low-resolution structures and CDRs with high B-factors or high conformational energies. We used a distance function based on directional statistics and an effective clustering algorithm with affinity propagation. With this data set of over 300 nonredundant antibody structures, we were able to cover 28 CDR-length combinations (e.g., L1 length 11, or “L1-11” in our CDR-length nomenclature) for L1, L2, L3, H1, and H2. The Chothia analysis covered only 20 CDR-lengths. Only four of these had more than one conformational cluster, of which two could easily be distinguished by gene source (mouse/human; κ/λ) and one could easily be distinguished purely by the presence and the positions of Pro residues (L3-9). Thus, using the Chothia analysis does not require the complicated set of “structure-determining residues” that is often assumed. Of our 28 CDR-lengths, 15 have multiple conformational clusters, including 10 for which the Chothia analysis had only one canonical class. We have a total of 72 clusters for non-H3 CDRs; approximately 85% of the non-H3 sequences can be assigned to a conformational cluster based on gene source and/or sequence. We found that earlier predictions of “bulged” versus “nonbulged” conformations based on the presence or the absence of anchor residues Arg/Lys94 and Asp101 of H3 have not held up, since all four combinations lead to a majority of conformations that are bulged. Thus, the earlier analyses have been significantly enhanced by the increased data. We believe that the new classification will lead to improved methods for antibody structure prediction and design.     

Backbone Flexibility of CDR3 and Immune Recognition of Antigens

Conformational entropy is an important component of protein–protein interactions; however, there is no reliable method for computing this parameter. We have developed a statistical measure of residual backbone entropy in folded proteins by using the ?–ψ distributions of the 20 amino acids in common secondary structures. The backbone entropy patterns of amino acids within helix, sheet or coil form clusters that recapitulate the branching and hydrogen bonding properties of the side chains in the secondary structure type. The same types of residues in coil and sheet have identical backbone entropies, while helix residues have much smaller conformational entropies. We estimated the backbone entropy change for immunoglobulin complementarity-determining regions (CDRs) from the crystal structures of 34 low-affinity T-cell receptors and 40 high-affinity Fabs as a result of the formation of protein complexes. Surprisingly, we discovered that the computed backbone entropy loss of only the CDR3, but not all CDRs, correlated significantly with the kinetic and affinity constants of the 74 selected complexes. Consequently, we propose a simple algorithm to introduce proline mutations that restrict the conformational flexibility of CDRs and enhance the kinetics and affinity of immunoglobulin interactions. Combining the proline mutations with rationally designed mutants from a previous study led to 2400-fold increase in the affinity of the A6 T-cell receptor for Tax-HLAA2. However, this mutational scheme failed to induce significant binding changes in the already-high-affinity C225–Fab/huEGFR interface. Our results will serve as a roadmap to formulate more effective target functions to design immune complexes with improved biological functions.     

Exploring the capacity of minimalist protein interfaces: interface energetics and affinity maturation to picomolar KD of a single-domain antibody with a flat paratope

A major architectural class in engineered binding proteins ("antibody mimics") involves the presentation of recognition loops off a single-domain scaffold. This class of binding proteins, both natural and synthetic, has a strong tendency to bind a preformed cleft using a convex binding interface (paratope). To explore their capacity to produce high-affinity interfaces with diverse shape and topography, we examined the interface energetics and explored the affinity limit achievable with a flat paratope. We chose a minimalist paratope limited to two loops found in a natural camelid heavy-chain antibody (VHH) that binds to ribonuclease A. Ala scanning of the VHH revealed only three "hot spot" side chains and additional four residues important for supporting backbone-mediated interactions. The small number of critical residues suggested that this is not an optimized paratope. Using selection from synthetic combinatorial libraries, we enhanced its affinity by >100-fold, resulting in variants with Kd as low as 180 pM with no detectable loss of binding specificity. High-resolution crystal structures revealed that the mutations induced only subtle structural changes but extended the network of interactions. This resulted in an expanded hot spot region including four additional residues located at the periphery of the paratope with a concomitant loss of the so-called "O-ring" arrangement of energetically inert residues. These results suggest that this class of simple, single-domain scaffolds is capable of generating high-performance binding interfaces with diverse shape. More generally, they suggest that highly functional interfaces can be designed without closely mimicking natural interfaces.     

Distinct roles of overlapping and non-overlapping regions of hub protein interfaces in recognition of multiple partners

Cellular functions of an organism are maintained by protein-protein interactions. Those proteins that bind multiple partners asynchronously (date hub proteins) are important to make the interaction network coordinated. It is known that many date hub proteins bind different partners at overlapping (OV) interfaces. To understand how OV interfaces of date hub proteins can recognize multiple partners, we analyzed the difference between OV and non-overlapping (Non-OV) regions of interfaces involved in the binding of different partners. By using the structures of 16 date hub proteins with various interaction partners (ranging from 5 to 33), we compared buried surface area, compositions of amino acid residues and secondary structures, and side-chain orientations. It was found that buried interface residues are important for recognizing multiple partners, while exposed interface residues are important for determining specificity to a particular ligand. In addition, our analyses reveal that residue compositions in OV and Non-OV regions are different and that residues in OV region show diverse side-chain torsion angles to accommodate binding to multiple targets.     

Structural Analysis of Peroxide-Soaked MnSOD Crystals Reveals Side-On Binding of Peroxide to Active-Site Manganese

The superoxide dismutase (SOD) enzymes are important antioxidant agents that protect cells from reactive oxygen species. The SOD family is responsible for catalyzing the disproportionation of superoxide radical to oxygen and hydrogen peroxide. Manganese- and iron-containing SOD exhibit product inhibition whereas Cu/ZnSOD does not. Here, we report the crystal structure of Escherichia coli MnSOD with hydrogen peroxide cryotrapped in the active site. Crystallographic refinement to 1.55 Å and close inspection revealed electron density for hydrogen peroxide in three of the four active sites in the asymmetric unit. The hydrogen peroxide molecules are in the position opposite His26 that is normally assumed by water in the trigonal bipyramidal resting state of the enzyme. Hydrogen peroxide is present in active sites B, C, and D and is side-on coordinated to the active-site manganese. In chains B and D, the peroxide is oriented in the plane formed by manganese and ligands Asp167 and His26. In chain C, the peroxide is bound, making a 70° angle to the plane. Comparison of the peroxide-bound active site with the hydroxide-bound octahedral form shows a shifting of residue Tyr34 towards the active site when peroxide is bound. Comparison with peroxide-soaked Cu/ZnSOD indicates end-on binding of peroxide when the SOD does not exhibit inhibition by peroxide and side-on binding of peroxide in the product-inhibited state of MnSOD.     

An Alternative Conformation of the T-Cell Receptor α Constant Region

αβ T-cell receptors (TcRs) play a central role in cellular immune response. They are members of the Ig superfamily, with extracellular regions of the α and β chains each comprising a V-type domain and a C-type domain. We have determined the ectodomain structure of an αβ TcR, which recognizes the autoantigen myelin basic protein. The 2.0-Å-resolution structure reveals canonical main-chain conformations for the V^α, V^β, and C^β domains, but the C^α domain exhibits a main-chain conformation remarkably different from those previously reported for TcR crystal structures. The global IgC-like fold is maintained, but a piston-like rearrangement between BC and DE β-turns results in β-strand slippage. This substantial conformational change may represent a signaling intermediate. Our structure is the first example for the Ig fold of the increasingly recognized concept of “metamorphic proteins.”     

Systematic Mutational Analysis of Peptide Inhibition of the p53-MDM2/MDMX Interactions

Inhibition of the interaction between the tumor suppressor protein p53 and its negative regulators MDM2 and MDMX is of great interest in cancer biology and drug design. We previously reported a potent duodecimal peptide inhibitor, termed PMI (TSFAEYWNLLSP), of the p53-MDM2 and -MDMX interactions. PMI competes with p53 for MDM2 and MDMX binding at an affinity roughly 2 orders of magnitude higher than that of ^17-28p53 (ETFSDLWKLLPE) of the same length; both peptides adopt nearly identical α-helical conformations in the complexes, where the three highlighted hydrophobic residues Phe, Trp, and Leu dominate PMI or ^17-28p53 binding to MDM2 and MDMX. To elucidate the molecular determinants for PMI activity and specificity, we performed a systematic Ala scanning mutational analysis of PMI and ^17-28p53. The binding affinities for MDM2 and MDMX of a total of 35 peptides including 10 truncation analogs were quantified, affording a complete dissection of energetic contributions of individual residues of PMI and ^17-28p53 to MDM2 and MDMX association. Importantly, the N8A mutation turned PMI into the most potent dual-specific antagonist of MDM2 and MDMX reported to date, registering respective K_d values of 490 pM and 2.4 nM. The co-crystal structure of N8A-PMI-^25-109MDM2 was determined at 1.95 Å, affirming that high-affinity peptide binding to MDM2/MDMX necessitates, in addition to optimized intermolecular interactions, enhanced helix stability or propensity contributed by non-contact residues. The powerful empirical binding data and crystal structures present a unique opportunity for computational studies of peptide inhibition of the p53-MDM2/MDMX interactions.     

Structural and Thermodynamic Characterization of Pre- and Postpolymerization States in the F4 Fimbrial Subunit FaeG

Enterotoxigenic Escherichia coli expressing F4 fimbriae are the major cause of porcine colibacillosis and are responsible for significant death and morbidity in neonatal and postweaned piglets. Via the chaperone-usher pathway, F4 fimbriae are assembled into thin, flexible polymers mainly composed of the single-domain adhesin FaeG. The F4 fimbrial system has been labeled eccentric because the F4 pilins show some features distinct from the features of pilins of other chaperone-usher-assembled structures. In particular, FaeG is much larger than other pilins (27  versus ∼ 17 kDa), grafting an additional carbohydrate binding domain on the common immunoglobulin-like core. Structural data of FaeG during different stages of the F4 fimbrial biogenesis process, combined with differential scanning calorimetry measurements, confirm the general principles of the donor strand complementation/exchange mechanisms taking place during pilus biogenesis via the chaperone-usher pathway.     

Crystal structure of scytalidoglutamic peptidase with its first potent inhibitor provides insights into substrate specificity and catalysis

Scytalidoglutamic peptidase (SGP) from Scytalidium lignicolum is the founding member of the newly discovered\ family of peptidases, G1, so far found exclusively in fungi. The crystal structure of SGP revealed a previously undescribed fold for peptidases and a unique catalytic dyad of residues Gln53 and Glu136. Surprisingly, the beta-sandwich structure of SGP is strikingly similar to members of the carbohydrate-binding concanavalin A-like lectins/glucanases superfamily. By analogy with the active sites of aspartic peptidases, a mechanism employing nucleophillic attack by a water molecule activated by the general base functionality of Glu136 has been proposed. Here, we report the first crystal structures of SGP in complex with two transition state peptide analogs designed to mimic the tetrahedral intermediate of the proteolytic reaction. Of these two analogs, the one containing a central S-hydroxyl group is a potent sub-nanomolar inhibitor of SGP. The inhibitor binds non-covalently to the concave surface of the upper beta-sheet and enables delineation of the S4 to S3' substrate specificity pockets of the enzyme. Structural differences in these pockets account for the unique substrate preferences of SGP among peptidases having an acidic pH optimum. Inhibitor binding is accompanied by a structuring of the region comprising residues Tyr71-Gly80 from being mostly disordered in the apoenzyme and leading to positioning of crucial active site residues for establishing enzyme-inhibitor contacts. In addition, conformational rearrangements are seen in a disulfide bridged surface loop (Cys141-Cys148), which moves inwards, partially closing the open substrate binding cleft of the native enzyme. The non-hydrolysable scissile bond analog of the inhibitor is located in the active site forming close contacts with Gln53 and Glu136. The nucleophilic water molecule is displaced and a unique mode of binding is observed with the S-OH of the inhibitor occupying the oxyanion binding site of the proposed tetrahedral intermediate. Details of the enzyme-inhibitor interactions and mechanistic interpretations are discussed.     

Structural details of HIV-1 recognition by the broadly neutralizing monoclonal antibody 2F5: epitope conformation, antigen-recognition loop mobility, and anion-binding site

2F5 is a monoclonal antibody with potent and broadly neutralizing activity against HIV-1. It targets the membrane-proximal external region (MPER) of the gp41 subunit of the envelope glycoprotein and interferes with the process of fusion between viral and host cell membranes. This study presents eight 2F5 F_ab′ crystal structures in complex with various gp41 peptide epitopes. These structures reveal several key features of this antibody-antigen interaction. (1) Whenever free of contacts caused by crystal artifacts, the extended complementarity-determining region H3 loop is mobile; this is true for ligand-free and epitope-bound forms. (2) The interaction between the antibody and the gp41 ELDKWA epitope core is absolutely critical, and there are also close and specific contacts with residues located N-terminal to the epitope core. (3) Residues located at the C-terminus of the gp41 ELDKWA core do not interact as tightly with the antibody. However, in the presence of a larger peptide containing the gp41 fusion peptide segment, these residues adopt a conformation consistent with the start of an α-helix. (4) At high sulfate concentrations, the electron density maps of 2F5 F_ab′-peptide complexes contain a peak that may mark a binding site for phosphate groups of negatively charged lipid headgroups. The refined atomic-level details of 2F5 paratope-epitope interactions revealed here should contribute to a better understanding of the mechanism of 2F5-based virus neutralization, in general, and prove important for the design of potential vaccine candidates intended to elicit 2F5-like antibody production.     

DRoP: A Water Analysis Program Identifies Ras-GTP-Specific Pathway of Communication between Membrane-Interacting Regions and the Active Site

Ras GTPase mediates several cellular signal transduction pathways and is found mutated in a large number of cancers. It is active in the GTP-bound state, where it interacts with effector proteins, and at rest in the GDP-bound state. The catalytic domain is tethered to the membrane, with which it interacts in a nucleotide-dependent manner. Here we present the program Detection of Related Solvent Positions (DRoP) for crystallographic water analysis on protein surfaces and use it to study Ras. DRoP reads and superimposes multiple Protein Data Bank coordinates, transfers symmetry-related water molecules to the position closest to the protein surface, and ranks the waters according to how well conserved and tightly clustered they are in the set of structures. Coloring according to this rank allows visualization of the results. The effector-binding region of Ras is hydrated with highly conserved water molecules at the interface between the P-loop, switch I, and switch II, as well as at the Raf-RBD binding pocket. Furthermore, we discovered a new conserved water-mediated H-bonding network present in Ras-GTP, but not in Ras-GDP, that links the nucleotide sensor residues R161 and R164 on helix 5 to the active site. The double mutant RasN85A/N86A, where the final link between helix 5 and the nucleotide is not possible, is a severely impaired enzyme, while the single mutant RasN86A, with partial connection to the active site, has a wild-type hydrolysis rate. DRoP was instrumental in determining the water-mediated connectivity networks that link two lobes of the catalytic domain in Ras.     

Structure-Based Prediction of the Peptide Sequence Space Recognized by Natural and Synthetic PDZ Domains

Protein-protein recognition, frequently mediated by members of large families of interaction domains, is one of the cornerstones of biological function. Here, we present a computational, structure-based method to predict the sequence space of peptides recognized by PDZ domains, one of the largest families of recognition proteins. As a test set, we use a considerable amount of recent phage display data that describe the peptide recognition preferences for 169 naturally occurring and engineered PDZ domains. For both wild-type PDZ domains and single point mutants, we find that 70-80% of the most frequently observed amino acids by phage display are predicted within the top five ranked amino acids. Phage display frequently identified recognition preferences for amino acids different from those present in the original crystal structure. Notably, in about half of these cases, our algorithm correctly captures these preferences, indicating that it can predict mutations that increase binding affinity relative to the starting structure. We also find that we can computationally recapitulate specificity changes upon mutation, a key test for successful forward design of protein-protein interface specificity. Across all evaluated data sets, we find that incorporation backbone sampling improves accuracy substantially, irrespective of using a crystal or NMR structure as the starting conformation. Finally, we report successful prediction of several amino acid specificity changes from blind tests in the DREAM4 peptide recognition domain specificity prediction challenge. Because the foundational methods developed here are structure based, these results suggest that the approach can be more generally applied to specificity prediction and redesign of other protein-protein interfaces that have structural information but lack phage display data.     

Structural Basis of Conserved Cysteine in the Fibroblast Growth Factor Family: Evidence for a Vestigial Half-Cystine

The 22 members of the mouse/human fibroblast growth factor (FGF) family of proteins contain a conserved cysteine residue at position 83 (numbering scheme of the 140-residue form of FGF-1). Sequence and structure information suggests that this position is a free cysteine in 16 members and participates as a half-cystine in at least 3 (and perhaps as many as 6) other members. While a structural role as a half-cystine provides a stability basis for possible selective pressure, it is less clear why this residue is conserved as a free cysteine (although free buried thiols can limit protein functional half-life). To probe the structural role of the free cysteine at position 83 in FGF-1, we constructed Ala, Ser, Thr, Val, and Ile mutations and determined their effects on structure and stability. These results show that position 83 in FGF-1 is thermodynamically optimized to accept a free cysteine. A second cysteine mutation was introduced into wild-type FGF-1 at adjacent position Ala66, which is known to participate as a half-cystine with position 83 in FGF-8, FGF-19, and FGF-23. Results show that, unlike position 83, a free cysteine at position 66 destabilizes FGF-1; however, upon oxidation, a near-optimal disulfide bond is formed between Cys66 and Cys83, resulting in ∼ 14 kJ/mol of increased thermostability. Thus, while the conserved free cysteine at position 83 in the majority of the FGF proteins may have a principal role in limiting functional half-life, evidence suggests that it is a vestigial half-cystine.     

Rapid Structural Characterization of Human Antibody-Antigen Complexes through Experimentally Validated Computational Docking

If we understand the structural rules governing antibody (Ab)-antigen (Ag) interactions in a given virus, then we have the molecular basis to attempt to design and synthesize new epitopes to be used as vaccines or optimize the antibodies themselves for passive immunization. Comparing the binding of several different antibodies to related Ags should also further our understanding of general principles of recognition.To obtain and compare the three-dimensional structure of a large number of different complexes, however, we need a faster method than traditional experimental techniques. While biocomputational docking is fast, its results might not be accurate. Combining experimental validation with computational prediction may be a solution.As a proof of concept, here we isolated a monoclonal Ab from the blood of a human donor recovered from dengue virus infection, characterized its immunological properties, and identified its epitope on domain III of dengue virus E protein through simple and rapid NMR chemical shift mapping experiments. We then obtained the three-dimensional structure of the Ab/Ag complex by computational docking, using the NMR data to drive and validate the results. In an attempt to represent the multiple conformations available to flexible Ab loops, we docked several different starting models and present the result as an ensemble of models equally agreeing with the experimental data. The Ab was shown to bind a region accessible only in part on the viral surface, explaining why it cannot effectively neutralize the virus.     

Determination of the transition state ensemble for the folding of ubiquitin from a combination of Phi and Psi analyses

Protein engineering techniques have emerged as powerful tools for characterizing transition states (TSs) for protein folding. Recently, the Ψ analysis, in which double-histidine mutations create the possibility of reversible crosslinking in the native state, has been proposed as an additional approach to the well-established Φ analysis. We present here a combination of these two procedures for defining the structure of the TS of ubiquitin, a small α/β protein that has been used extensively as a model system for both experimental and computational studies of the protein-folding process. We performed a series of molecular dynamics simulations in which Φ and Ψ values were used as ensemble-averaged structural restraints to determine an ensemble of structures representing the TS of ubiquitin. Although the available Ψ values for ubiquitin did not, by themselves, generate well-defined TS ensembles, the inclusion of the restricted set of zero or unity values, but not fractional ones, provided useful complementary information to the Φ analysis. Our results show that the TS of ubiquitin is formed by a relatively narrow ensemble of structures exhibiting an overall native-like topology in which the N-terminal and C-terminal regions are in close proximity.