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.     

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.     

Energy-Based Analysis and Prediction of the Orientation between Light- and Heavy-Chain Antibody Variable Domains

Diversity in antibody structure is crucial to the ability of the adaptive immune system to recognize the tremendously diverse set of potential antigens. The diversity in structure is most apparent in the six hypervariable loops of the complementarity-determining regions. However, given that these loops occur at the interface of the heavy- and light-chain variable domains and form the antigen-binding site, the relative orientation of the heavy- and light-chain variable domains can create another source of structural diversity leading to changes in antigen binding. Here, we first reexamine the diversity of V_L:V_H orientations in existing antibody crystal structures using 153 nonredundant sequences, demonstrating that the variation in V_L:V_H orientation is greater than that expected from effects of crystal packing, antigen binding, or the presence of antibody constant regions and increases, on average, as sequence similarity decreases for residues in the interface between the domains. We developed a tool for predicting the relative orientations of the heavy- and light-chain variable domains using side-chain rotamer sampling in the interface and molecular-mechanics-based energy calculations. When using variable domain backbones from the crystal structures, the predicted orientation is very close (< 1 Å RMSD) to the crystallographically observed orientation in most cases, confirming that the V_L:V_H orientation is determined by the antibody sequence and suggesting an approach to predicting the relative orientation of the variable domains when building homology models of antibodies. When applied to antibody homology models generated from templates with 55-75% sequence identity, we predict the V_L:V_H orientation of 20 antibodies with an average/median RMSD of 2.1/1.6 Å to the crystal structures.     

Structure-Function Relationship of Tumor Necrosis Factor (TNF) and Its Receptor Interaction Based on 3D Structural Analysis of a Fully Active TNFR1-Selective TNF Mutant

Tumor necrosis factor (TNF) is an important cytokine that suppresses carcinogenesis and excludes infectious pathogens to maintain homeostasis. TNF activates its two receptors [TNF receptor (TNFR) 1 and TNFR2], but the contribution of each receptor to various host defense functions and immunologic surveillance is not yet clear. Here, we used phage display techniques to generate receptor-selective TNF mutants that activate only one TNFR. These TNF mutants will be useful in the functional analysis of TNFR.Six amino acids in the receptor binding interface (near TNF residues 30, 80, and 140) were randomly mutated by polymerase chain reaction. Two phage libraries comprising over 5 million TNF mutants were constructed. By selecting the mutants without affinity for TNFR1 or TNFR2, we successfully isolated 4 TNFR2-selective candidates and 16 TNFR1-selective candidates, respectively. The TNFR1-selective candidates were highly mutated near residue 30, whereas TNFR2-selective candidates were highly mutated near residue 140, although both had conserved sequences near residues 140 and 30, respectively. This finding suggested that the phage display technique was suitable for identifying important regions for the TNF interaction with TNFR1 and TNFR2. Purified clone R1-6, a TNFR1-selective candidate, remained fully bioactive and had full affinity for TNFR1 without activating TNFR2, indicating the usefulness of the R1-6 TNF mutant in analyzing TNFR1 receptor function.To further elucidate the receptor selectivity of R1-6, we examined the structure of R1-6 by X-ray crystallography. The results suggested that R31A and R32G mutations strongly influenced electrostatic interaction with TNFR2, and that L29K mutation contributed to the binding of R1-6 to TNFR1. This phage display technique can be used to efficiently construct functional mutants for analysis of the TNF structure-function relationship, which might facilitate in silico drug design based on receptor selectivity.     

Plasticity in the Contribution of T Cell Receptor Variable Region Residues to Binding of Peptide–HLA-A2 Complexes

One hypothesis accounting for major histocompatibility complex (MHC) restriction by T cell receptors (TCRs) holds that there are several evolutionary conserved residues in TCR variable regions that contact MHC. While this “germline codon” hypothesis is supported by various lines of evidence, it has been difficult to test. The difficulty stems in part from the fact that TCRs exhibit low affinities for pep/MHC, thus limiting the range of binding energies that can be assigned to these key interactions using mutational analyses. To measure the magnitude of binding energies involved, here we used high-affinity TCRs engineered by mutagenesis of CDR3. The TCRs included a high-affinity, MART-1/HLA-A2-specific single-chain TCR and two other high-affinity TCRs that all contain the same Vα region and recognize the same MHC allele (HLA-A2), with different peptides and Vβ regions. Mutational analysis of residues in CDR1 and CDR2 of the three Vα2 regions showed the importance of the key germline codon residue Y51. However, two other proposed key residues showed significant differences among the TCRs in their relative contributions to binding. With the use of single-position, yeast-display libraries in two of the key residues, MART-1/HLA-A2 selections also revealed strong preferences for wild-type germline codon residues, but several alternative residues could also accommodate binding and, hence, MHC restriction. Thus, although a single residue (Y51) could account for a proportion of the energy associated with positive selection (i.e., MHC restriction), there is significant plasticity in requirements for particular side chains in CDR1 and CDR2 and in their relative binding contributions among different TCRs.     

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.”     

Crystal structure of vascular endothelial growth factor-B in complex with a neutralising antibody Fab fragment

Vascular endothelial growth factor (VEGF) B effects blood vessel formation by binding to VEGF receptor 1. To study the specifics of the biological profile of VEGF-B in both physiological and pathological angiogenesis, a neutralising anti-VEGF-B antibody (2H10) that functions by inhibiting the binding of VEGF-B to VEGF receptor 1 was developed. Here, we present the structural features of the ‘highly ordered’ interaction of the Fab fragment of this antibody (Fab-2H10) with VEGF-B. Two molecules of Fab-2H10 bind to symmetrical binding sites located at each pole of the VEGF-B homodimer, giving a unique U-shaped topology to the complex that has not been previously observed in the VEGF family. VEGF-B residues essential for binding to the antibody are contributed by both monomers of the cytokine. Our detailed analysis reveals that the neutralising effect of the antibody occurs by virtue of the steric hindrance of the receptor-binding interface. These findings suggest that functional complementarity between VEGF-B and 2H10 can be harnessed both in analysing the therapeutic potential of VEGF-B and as an antagonist of receptor activation.     

Engineering a high-affinity anti-IL-15 antibody: crystal structure reveals an α-helix in VH CDR3 as key component of paratope

Interleukin (IL) 15 is an inflammatory cytokine that plays an essential role in the activation, proliferation, and maintenance of specific natural killer cell and T-cell populations, and has been implicated as a mediator of inflammatory diseases. An anti-IL-15 antibody that blocked IL-15-dependent cellular responses was isolated by phage display and optimised via mutagenesis of the third complementarity-determining regions (CDRs) of variable heavy (VH) and variable light chains. Entire repertoires of improved variants were recombined with each other to explore the maximum potential sequence space. DISC0280, the most potent antibody isolated using this comprehensive strategy, exhibits a 228-fold increase in affinity and a striking 40,000-fold increase in cellular potency compared to its parent. Such a wholesale recombination strategy therefore represents a useful method for exploiting synergistic potency gains as part of future antibody engineering efforts. The crystal structure of DISC0280 Fab (fragment antigen binding), in complex with human IL-15, was determined in order to map the structural epitope and paratope. The most remarkable feature revealed lies within the paratope and is a novel six-amino-acid α-helix that sits within the VH CDR3 loop at the center of the antigen binding site. This is the first report to describe an α-helix as a principal component of a naturally derived VH CDR3 following affinity maturation.     

A dominant conformational role for amino acid diversity in minimalist protein-protein interfaces

Recent studies have shown that highly simplified interaction surfaces consisting of combinations of just two amino acids, Tyr and Ser, exhibit high affinity and specificity. The high functional levels of such minimalist interfaces might thus indicate small contributions of greater amino acid diversity seen in natural interfaces. Toward addressing this issue, we have produced a pair of binding proteins built on the fibronectin type III scaffold, termed “monobodies.” One monobody contains the Tyr/Ser binary-code interface (termed YS) and the other contains an expanded amino acid diversity interface (YSX), but both bind to an identical target, maltose-binding protein. The YSX monobody bound with higher affinity, a slower off rate and a more favorable enthalpic contribution than the YS monobody. High-resolution X-ray crystal structures revealed that both proteins bound to an essentially identical epitope, providing a unique opportunity to directly investigate the role of amino acid diversity in a protein interaction interface. Surprisingly, Tyr still dominates the YSX paratope and the additional amino acid types are primarily used to conformationally optimize contacts made by tyrosines. Scanning mutagenesis showed that while all contacting Tyr side chains are essential in the YS monobody, the YSX interface was more tolerant to mutations. These results suggest that the conformational, not chemical, diversity of additional types of amino acids provided higher functionality and evolutionary robustness, supporting the dominant role of Tyr and the importance of conformational diversity in forming protein interaction interfaces.     

Reconstitution and Engineering of Apoptotic Protein Interactions on the Bacterial Cell Surface

The interactions between pro- and anti-apoptotic members of the Bcl-2 class of proteins control whether a cell lives or dies, and the study of these protein-protein interactions has been an area of intense research. In this report, we describe a new tool for the study and engineering of apoptotic protein interactions that is based on the flow cytometric detection of these interactions on the surface of Escherichia coli. After validation of the assay with the well-studied interaction between the Bak(72-87) peptide and the anti-apoptotic protein Bcl-x_L, the effect of both increasing and decreasing Bak peptide length on Bcl-x_L binding was investigated. Previous work demonstrated that the Bak(72-87) peptide also binds to the anti-apoptotic protein Bcl-2, albeit with lower binding affinity compared to Bcl-x_L. Here, we demonstrate that a slightly longer Bak peptide corresponding to amino acids 72-89 of Bak binds Bcl-x_L and Bcl-2 equally well. Approximate binding affinity calculations on these peptide-protein complexes confirm the experimental observations. The flow cytometric assay was also used to screen a saturation mutagenesis library of Bak(72-87) variants for improved affinity to Bcl-x_L. The best variants obtained from this library exhibit an apparent K_d to Bcl-x_L 4-fold lower than that of wild-type Bak(72-87).     

Crystal structure of the LG3 domain of endorepellin, an angiogenesis inhibitor

Endorepellin, the C-terminal region of perlecan, inhibits angiogenesis by disrupting actin cytoskeleton and focal adhesions. The C-terminal laminin-like globular domain (LG3) of endorepellin directs most of this antiangiogenic activity. To investigate the angiostatic mechanism and to identify structural determinants, we have solved crystal structures of the LG3 domain in both apo- and calcium-bound forms at resolutions of 1.5 Å and 2.8 Å, respectively. The conserved core has the jellyroll fold characteristic of LG domains. The calcium-induced structural changes seem very restricted, and the calcium binding site appears to be preformed, suggesting that the bound calcium ion, rather than structural rearrangements, contributes to antiangiogenesis. We have identified H4268 on the EF loop as a key residue for the biochemical function of LG3, since its mutation abolishes antiangiogenic activity, and mutant LG3 can no longer form a direct interaction with integrin. Taken together, we propose that these two distinct structural elements contribute to the angiostatic effect of endorepellin.     

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.     

The role of RNA sequence and structure in RNA--protein interactions

We investigate the sequence and structural properties of RNA-protein interaction sites in 211 RNA-protein chain pairs, the largest set of RNA-protein complexes analyzed to date. Statistical analysis confirms and extends earlier analyses made on smaller data sets. There are 24.6% of hydrogen bonds between RNA and protein that are nucleobase specific, indicating the importance of both nucleobase-specific and -nonspecific interactions. While there is no significant difference between RNA base frequencies in protein-binding and non-binding regions, distinct preferences for RNA bases, RNA structural states, protein residues, and protein secondary structure emerge when nucleobase-specific and -nonspecific interactions are considered separately. Guanine nucleobase and unpaired RNA structural states are significantly preferred in nucleobase-specific interactions; however, nonspecific interactions disfavor guanine, while still favoring unpaired RNA structural states. The opposite preferences of nucleobase-specific and -nonspecific interactions for guanine may explain discrepancies between earlier studies with regard to base preferences in RNA-protein interaction regions. Preferences for amino acid residues differ significantly between nucleobase-specific and -nonspecific interactions, with nonspecific interactions showing the expected bias towards positively charged residues. Irregular protein structures are strongly favored in interactions with the protein backbone, whereas there is little preference for specific protein secondary structure in either nucleobase-specific interaction or -nonspecific interaction. Overall, this study shows strong preferences for both RNA bases and RNA structural states in protein-RNA interactions, indicating their mutual importance in protein recognition.     

A disulfide-free single-domain V(L) intrabody with blocking activity towards huntingtin reveals a novel mode of epitope recognition

We present the crystal structure and biophysical characterization of a human V_L [variable domain immunoglobulin (Ig) light chain] single-domain intrabody that binds to the huntingtin (Htt) protein and has been engineered for antigen recognition in the absence of its intradomain disulfide bond, otherwise conserved in the Ig fold. Analytical ultracentrifugation demonstrated that the αHtt-V_L 12.3 domain is a stable monomer under physiological conditions even at concentrations > 20 μM. Using peptide SPOT arrays, we identified the minimal binding epitope to be EKLMKAFESLKSFQ, comprising the N-terminal residues 5-18 of Htt and including the first residue of the poly-Gln stretch. X-ray structural analysis of αHtt-V_L both as apo protein and in the presence of the epitope peptide revealed several interesting insights: first, the role of mutations acquired during the combinatorial selection process of the αHtt-V_L 12.3 domain—initially starting from a single-chain Fv fragment—that are responsible for its stability as an individually soluble Ig domain, also lacking the disulfide bridge, and second, a previously unknown mode of antigen recognition, revealing a novel paratope. The Htt epitope peptide adopts a purely α-helical structure in the complex with αHtt-V_L and is bound at the base of the complementarity-determining regions (CDRs) at the concave β-sheet that normally gives rise to the interface between the V_L domain and its paired V_H (variable domain Ig heavy chain) domain, while only few interactions with CDR-L1 and CDR-L3 are formed. Notably, this noncanonical mode of antigen binding may occur more widely in the area of in vitro selected antibody fragments, including other Ig-like scaffolds, possibly even if a V_H domain is present.     

Structural and thermodynamic signatures of DNA recognition by Mycobacterium tuberculosis DnaA

An essential protein, DnaA, binds to 9-bp DNA sites within the origin of replication oriC. These binding events are prerequisite to forming an enigmatic nucleoprotein scaffold that initiates replication. The number, sequences, positions, and orientations of these short DNA sites, or DnaA boxes, within the oriCs of different bacteria vary considerably. To investigate features of DnaA boxes that are important for binding Mycobacterium tuberculosis DnaA (MtDnaA), we have determined the crystal structures of the DNA binding domain (DBD) of MtDnaA bound to a cognate MtDnaA-box (at 2.0 Å resolution) and to a consensus Escherichia coli DnaA-box (at 2.3 Å). These structures, complemented by calorimetric equilibrium binding studies of MtDnaA DBD in a series of DnaA-box variants, reveal the main determinants of DNA recognition and establish the [T/C][T/A][G/A]TCCACA sequence as a high-affinity MtDnaA-box. Bioinformatic and calorimetric analyses indicate that DnaA-box sequences in mycobacterial oriCs generally differ from the optimal binding sequence. This sequence variation occurs commonly at the first 2 bp, making an in vivo mycobacterial DnaA-box effectively a 7-mer and not a 9-mer. We demonstrate that the decrease in the affinity of these MtDnaA-box variants for MtDnaA DBD relative to that of the highest-affinity box TTGTCCACA is less than 10-fold. The understanding of DnaA-box recognition by MtDnaA and E. coli DnaA enables one to map DnaA-box sequences in the genomes of M. tuberculosis and other eubacteria.     

RCSB Protein Data bank: Tools for visualizing and understanding biological macromolecules in 3D

Now in its 52nd year of continuous operations, the Protein Data Bank (PDB) is the premiere open‐access global archive housing three‐dimensional (3D) biomolecular structure data. It is jointly managed by the Worldwide Protein Data Bank (wwPDB) partnership. The Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB) is funded by the National Science Foundation, National Institutes of Health, and US Department of Energy and serves as the US data center for the wwPDB. RCSB PDB is also responsible for the security of PDB data in its role as wwPDB‐designated Archive Keeper. Every year, RCSB PDB serves tens of thousands of depositors of 3D macromolecular structure data (coming from macromolecular crystallography, nuclear magnetic resonance spectroscopy, electron microscopy, and micro‐electron diffraction). The RCSB PDB research‐focused web portal (RCSB.org) makes PDB data available at no charge and without usage restrictions to many millions of PDB data consumers around the world. The RCSB PDB training, outreach, and education web portal (PDB101.RCSB.org) serves nearly 700 K educators, students, and members of the public worldwide. This invited Tools Issue contribution describes how RCSB PDB (i) is organized; (ii) works with wwPDB partners to process new depositions; (iii) serves as the wwPDB‐designated Archive Keeper; (iv) enables exploration and 3D visualization of PDB data via RCSB.org; and (v) supports training, outreach, and education via PDB101.RCSB.org. New tools and features at RCSB.org are presented using examples drawn from high‐resolution structural studies of proteins relevant to treatment of human cancers by targeting immune checkpoints.