Mutations designed to destabilize the receptor-bound conformation increase MICA-NKG2D association rate and affinity

Preclinical characterization of Pan-NKG2D ligand-binding NKG2D receptor decoys

NKG2D and its ligands are critical regulators of protective immune responses controlling infections and cancer, defining a crucial immune signaling axis. Current therapeutic efforts targeting this axis almost exclusively aim at enhancing NKG2D-mediated effector functions. However, this axis can drive disease processes when dysregulated, in particular, driving stem-like cancer cell reprogramming and tumorigenesis through receptor/ligand self-stimulation on tumor cells. Despite complexities with its structure and biology, we developed multiple novel engineered proteins that functionally serve as axis-blocking NKG2D "decoys" and report biochemical, structural, in vitro, and in vivo evaluation of their functionality.

Targeting the NKG2D/NKG2D-L axis in acute myeloid leukemia

Natural killer group 2, member D (NKG2D) receptor is a crucial activating receptor in the immune recognition and eradication of abnormal cells by natural killer (NK) cells, and T lymphocytes. NKG2D can transmit activation signals and activate the immune system by recognizing the NKG2D ligands (NKG2D-L) on acute myeloid leukemia (AML) cells. Downregulation of NKG2D-L in AML can circumvent resistance to chemotherapy and immune recognition. Considering this effect, the exploration of targeting the NKG2D/NKG2D-L axis is considered to have tremendous potential for the discovery of novel biomacromolecule antibodies and pharmacological modulators in AML. This review was to outline the impact of NKG2D/NKG2D-L axis on intrinsic immunosurveillance and the development of AML. Furthermore, the NKG2D/NKG2D-L axis related modulators and progress in preclinical and clinical trials was also to be reviewed.

Alternative Splice Transcripts for MHC Class I-like MICA Encode Novel NKG2D Ligands with Agonist or Antagonist Functions

MHC class I chain-related proteins A and B (MICA and MICB) and UL16-binding proteins are ligands of the activating NKG2D receptor involved in cancer and immune surveillance of infection. Structurally, MICA/B proteins contain an α3 domain, whereas UL16-binding proteins do not. We identified novel alternative splice transcripts for MICA encoding five novel MICA isoforms: MICA-A, -B1, -B2, -C, and -D. Alternative splicing associates with MICA*015 and *017 and results from a point deletion (G) in the 5' splice donor site of MICA intron 4 leading to exon 3 and exon 4 skipping and/or deletions. These changes delete the α3 domain in all isoforms, and the α2 domain in the majority of isoforms (A, B1, C, and D). Endothelial and hematopoietic cells contained endogenous alternative splice transcripts and isoforms. MICA-B1, -B2, and -D bound NKG2D by surface plasmon resonance and were expressed at the cell surface. Functionally, MICA-B2 contains two extracellular domains (α1 and α2) and is a novel potent agonist ligand for NKG2D. We found that MICA-D is a new truncated form of MICA with weak affinity for NKG2D despite lacking α2 and α3 domains. MICA-D may functionally impair NKG2D activation by competing with full-length MICA or MICA-B2 for NKG2D engagement. Our study established NKG2D binding for recombinant MICA-B1 but found no function for this isoform. New truncated MICA isoforms exhibit a range of functions that may drive unexpected immune mechanisms and provide new tools for immunotherapy.

Genetics, genomics, and evolutionary biology of NKG2D ligands

Human and mouse NKG2D ligands (NKG2DLs) are absent or only poorly expressed by most normal cells but are upregulated by cell stress, hence, alerting the immune system in case of malignancy or infection. Although these ligands are numerous and highly variable (at genetic, genomic, structural, and biochemical levels), they all belong to the major histocompatibility complex class I gene superfamily and bind to a single, invariant, receptor: NKG2D. NKG2D (CD314) is an activating receptor expressed on NK cells and subsets of T cells that have a key role in the recognition and lysis of infected and tumor cells. Here, we review the molecular diversity of NKG2DLs, discuss the increasing appreciation of their roles in a variety of medical conditions, and propose several explanations for the evolutionary force(s) that seem to drive the multiplicity and diversity of NKG2DLs while maintaining their interaction with a single invariant receptor.

Biochemical Basis for Dominant Inheritance, Variable Penetrance, and Maternal Effects in RBP4 Congenital Eye Disease

Gestational vitamin A (retinol) deficiency poses a risk for ocular birth defects and blindness. We identified missense mutations in RBP4, encoding serum retinol binding protein, in three families with eye malformations of differing severity, including bilateral anophthalmia. The mutant phenotypes exhibit dominant inheritance, but incomplete penetrance. Maternal transmission significantly increases the probability of phenotypic expression. RBP normally delivers retinol from hepatic stores to peripheral tissues, including the placenta and fetal eye. The disease mutations greatly reduce retinol binding to RBP, yet paradoxically increase the affinity of RBP for its cell surface receptor, STRA6. By occupying STRA6 nonproductively, the dominant-negative proteins disrupt vitamin A delivery from wild-type proteins within the fetus, but also, in the case of maternal transmission, at the placenta. These findings establish a previously uncharacterized mode of maternal inheritance, distinct from imprinting and oocyte-derived mRNA, and define a group of hereditary disorders plausibly modulated by dietary vitamin A.

Increasing affinity of interferon-γ receptor 1 to interferon-γ by computer-aided design

We describe a computer-based protocol to design protein mutations increasing binding affinity between ligand and its receptor. The method was applied to mutate interferon-γ receptor 1 (IFN-γ-Rx) to increase its affinity to natural ligand IFN-γ, protein important for innate immunity. We analyzed all four available crystal structures of the IFN-γ-Rx/IFN-γ complex to identify 40 receptor residues forming the interface with IFN-γ. For these 40 residues, we performed computational mutation analysis by substituting each of the interface receptor residues by the remaining standard amino acids. The corresponding changes of the free energy were calculated by a protocol consisting of FoldX and molecular dynamics calculations. Based on the computed changes of the free energy and on sequence conservation criteria obtained by the analysis of 32 receptor sequences from 19 different species, we selected 14 receptor variants predicted to increase the receptor affinity to IFN-γ. These variants were expressed as recombinant proteins in Escherichia coli, and their affinities to IFN-γ were determined experimentally by surface plasmon resonance (SPR). The SPR measurements showed that the simple computational protocol succeeded in finding two receptor variants with affinity to IFN-γ increased about fivefold compared to the wild-type receptor.

Structure of NKp65 bound to its keratinocyte ligand reveals basis for genetically linked recognition in natural killer gene complex

The natural killer (NK) gene complex (NKC) encodes numerous C-type lectin-like receptors that govern the activity of NK cells. Although some of these receptors (Ly49s, NKG2D, CD94/NKG2A) recognize MHC or MHC-like molecules, others (Nkrp1, NKRP1A, NKp80, NKp65) instead bind C-type lectin-like ligands to which they are genetically linked in the NKC. To understand the basis for this recognition, we determined the structure of human NKp65, an activating receptor implicated in the immunosurveillance of skin, bound to its NKC-encoded ligand keratinocyte-associated C-type lectin (KACL). Whereas KACL forms a homodimer resembling other C-type lectin-like dimers, NKp65 is monomeric. The binding mode in the NKp65-KACL complex, in which a monomeric receptor engages a dimeric ligand, is completely distinct from those used by Ly49s, NKG2D, or CD94/NKG2A. The structure explains the exceptionally high affinity of the NKp65-KACL interaction compared with other cell-cell interaction pairs (KD = 6.7 × 10(-10) M), which may compensate for the monomeric nature of NKp65 to achieve cell activation. This previously unreported structure of an NKC-encoded receptor-ligand complex, coupled with mutational analysis of the interface, establishes a docking template that is directly applicable to other genetically linked pairs in the NKC, including Nkrp1-Clr, NKRP1A-LLT1, and NKp80-AICL.

Slow, reversible, coupled folding and binding of the spectrin tetramerization domain

Many intrinsically disordered proteins (IDPs) are significantly unstructured under physiological conditions. A number of these IDPs have been shown to undergo coupled folding and binding reactions whereby they can gain structure upon association with an appropriate partner protein. In general, these systems display weaker binding affinities than do systems with association between completely structured domains, with micromolar K(d) values appearing typical. One such system is the association between α- and β-spectrin, where two partially structured, incomplete domains associate to form a fully structured, three-helix bundle, the spectrin tetramerization domain. Here, we use this model system to demonstrate a method for fitting association and dissociation kinetic traces where, using typical biophysical concentrations, the association reactions are expected to be highly reversible. We elucidate the unusually slow, two-state kinetics of spectrin assembly in solution. The advantages of studying kinetics in this regime include the potential for gaining equilibrium constants as well as rate constants, and for performing experiments with low protein concentrations. We suggest that this approach would be particularly appropriate for high-throughput mutational analysis of two-state reversible binding processes.

Combining different design strategies for rational affinity maturation of the MICA-NKG2D interface

We redesigned residues on the surface of MICA, a protein that binds the homodimeric immunoreceptor NKG2D, to increase binding affinity with a series of rational, incremental changes. A fixed-backbone RosettaDesign protocol scored a set of initial mutations, which we tested by surface plasmon resonance for thermodynamics and kinetics of NKG2D binding, both singly and in combination. We combined the best four mutations at the surface with three affinity-enhancing mutations below the binding interface found with a previous design strategy. After curating design scores with three cross-validated tests, we found a linear relationship between free energy of binding and design score, and to a lesser extent, enthalpy and design score. Multiple mutants bound with substantial subadditivity, but in at least one case full additivity was observed when combining distant mutations. Altogether, combining the best mutations from the two strategies into a septuple mutant enhanced affinity by 50-fold, to 50 nM, demonstrating a simple, effective protocol for affinity enhancement.

Preparative Protein Production from Inclusion Bodies and Crystallization: A Seven-Week Biochemistry Sequence

We describe how to produce and purify proteins from E. coli inclusion bodies by adapting versatile, preparative-scale techniques to the undergraduate laboratory schedule. This seven-week sequence of experiments fits into an annual cycle of research activity in biochemistry courses. Recombinant proteins are expressed as inclusion bodies, which are collected, washed, then solubilized in urea. Stepwise dialysis to dilute urea over the course of a week produces refolded protein. Column chromatography is used to purify protein into fractions, which are then analyzed with gel electrophoresis and concentration assays. Students culminate the project by designing crystallization trials in sitting-drop trays. Student evaluation of the experience has been positive, listing 5-12 new techniques learned, which are transferrable to graduate research in academia and industry.

A free-energy landscape for coupled folding and binding of an intrinsically disordered protein in explicit solvent from detailed all-atom computations

The N-terminal repressor domain of neural restrictive silencer factor (NRSF) is an intrinsically disordered protein (IDP) that binds to the paired amphipathic helix (PAH) domain of mSin3. An NMR experiment revealed that the minimal binding unit of NRSF is a 15-residue segment that adopts a helical structure upon binding to a cleft of mSin3. We computed a free-energy landscape of this system by an enhanced conformational sampling method, all-atom multicanonical molecular dynamics. The simulation started from a configuration where the NRSF segment was fully disordered and distant from mSin3 in explicit solvent. In the absence of mSin3, the disordered NRSF segment thermally fluctuated between hairpins, helices, and bent structures. In the presence of mSin3, the segment bound to mSin3 by adopting the structures involved in the isolated state, and non-native and native complexes were formed. The free-energy landscape comprised three superclusters, and free-energy barriers separated the superclusters. The native complex was located at the center of the lowest free-energy cluster. When NRSF landed in the largest supercluster, the generated non-native complex moved on the landscape to fold into the native complex, by increasing the interfacial hydrophobic contacts and the helix content. When NRSF landed in other superclusters, the non-native complex overcame the free-energy barriers between the various segment orientations in the binding cleft of mSin3. Both population-shift and induced-fit (or induced-folding) mechanisms work cooperatively in the coupled folding and binding. The diverse structural adaptability of NRSF may be related to the hub properties of the IDP.

Systematic mutation and thermodynamic analysis of central tyrosine pairs in polyspecific NKG2D receptor interactions

The homodimeric, activating natural killer cell receptor NKG2D interacts with multiple monomeric ligands polyspecifically, yet without central conformational flexibility. Crystal structures of multiple NKG2D-ligand interactions have identified the NKG2D tyrosine pair Tyr 152 and Tyr 199 as forming multiple specific but diverse interactions with MICA and related proteins. Here we systematically altered each tyrosine to tryptophan, phenylalanine, isoleucine, leucine, valine, serine, and alanine to measure the effect of mutation on affinity and thermodynamics for binding a range of similar ligands: MICA, the higher-affinity ligand MICB, and MICdesign, a high-affinity version of MICA that shares all NKG2D contact residues with MICA. Affinity and residue size were related: tryptophan could often substitute for tyrosine without loss of affinity; loss of the tyrosine hydroxyl through mutation to phenylalanine was tolerated more at position 152 than 199; and the smallest residues coincide with lowest affinities in general. NKG2D mutant van't Hoff binding thermodynamics generally show that substitution of other residues for tyrosine causes a moderate positive or flat van't Hoff slope consistent with moderate loss of binding enthalpy. One set of NKG2D mutations caused MICA to adopt a positive van't Hoff slope corresponding to absorption of heat, and another set caused MICB to adopt a negative slope of greater heat release than wild-type. MICdesign shared one example of the first set with MICA and one of the second set with MICB. When the NKG2D mutation affinities were arranged according to change in nonpolar surface area and compared to results from specific antibody-antigen and protein-peptide interactions, it was found that hydrophobic surface loss in NKG2D reduced binding affinity less than reported in the other contexts. The hydrophobic effect at the center of the NKG2D binding appears more similar to that at the periphery of an antibody-antigen binding site than at its center. Therefore the polyspecific NKG2D binding site is more tolerant of structural alteration in general than either an antibody-antigen or protein-peptide binding site, and this tolerance may adapt NKG2D to a broad range of protein surfaces with micromolar affinity.

RosettaBackrub--a web server for flexible backbone protein structure modeling and design

The RosettaBackrub server (http://kortemmelab.ucsf.edu/backrub) implements the Backrub method, derived from observations of alternative conformations in high-resolution protein crystal structures, for flexible backbone protein modeling. Backrub modeling is applied to three related applications using the Rosetta program for structure prediction and design: (I) modeling of structures of point mutations, (II) generating protein conformational ensembles and designing sequences consistent with these conformations and (III) predicting tolerated sequences at protein–protein interfaces. The three protocols have been validated on experimental data. Starting from a user-provided single input protein structure in PDB format, the server generates near-native conformational ensembles. The predicted conformations and sequences can be used for different applications, such as to guide mutagenesis experiments, for ensemble-docking approaches or to generate sequence libraries for protein design.

MICA polymorphism identified by whole genome array associated with NKG2D-mediated cytotoxicity in T-cell large granular lymphocyte leukemia

BACKGROUND

Large granular lymphocyte leukemia is a semi-autonomous clonal proliferation of cytotoxic T cells accompanied by immune cytopenias and various autoimmune conditions. Due to the rarity of this disease and its association with autoimmune diseases, a theoretical germline or somatic mutation might have significant penetrance, thus enabling detection, even from samples of suboptimal size, through genome-wide association studies.

DESIGN AND METHODS

To investigate a non-mendelian genetic predisposition to large granular lymphocyte leukemia, we used a step-wise method for gene discovery. First, a modified 'random forests' technique was used for candidate gene identification: this was followed by traditional allele-specific polymerase chain reaction, sequencing modalities, and mechanistic assays.

RESULTS

Our analysis found an association with MICA, a non-peptide-presenting, tightly regulated, stress-induced MHC-like molecule and cognate receptor for NKG2D, found abundantly on large granular lymphocyte leukemia cells. Sequencing of germline DNA revealed a higher frequency of MICA*00801/A5.1 in patients with large granular lymphocyte leukemia than in matched controls (64% versus 41%, P<0.001, homozygous 40% versus 15%, P<0.001). Flow cytometry was employed to determine the expression of MICA within hematologic compartments, showing that the signal intensity of MICA was increased in granulocytes from neutropenic patients with large granular lymphocyte leukemia in comparison with that in controls (P=0.033). Furthermore, neutrophil counts were inversely correlated with MICA expression (R(2)=0.50, P=0.035). Finally, large granular lymphocyte leukemia cells were able to selectively kill MICA(+) Ba/F3 lymphocytes transfected with human MICA*019 in a dose-dependent manner compared to naïve cells (P<0.001), an effect mitigated by administration of an anti-NKG2D antibody (P=0.033).

CONCLUSIONS

Our results illustrate that MICA-NKG2D played a role in disease pathogenesis in the majority of patients in our cohort of cases of large granular lymphocyte leukemia and further investigation into this signaling axis may provide potent therapeutic targets.

Kinetic advantage of intrinsically disordered proteins in coupled folding-binding process: a critical assessment of the "fly-casting" mechanism

Intrinsically disordered proteins (IDPs) are recognized to play important roles in many biological functions such as transcription and translation regulation, cellular signal transduction, protein phosphorylation, and molecular assemblies. The coupling of folding with binding through a "fly-casting" mechanism has been proposed to account for the fast binding kinetics of IDPs. In this article, experimental data from the literature were collated to verify the kinetic advantages of IDPs, while molecular simulations were performed to clarify the origin of the kinetic advantages. The phosphorylated KID-kinase-inducible domain interacting domain (KIX) complex was used as an example in the simulations. By modifying a coarse-grained model with a native-centric Gō-like potential, we were able to continuously tune the degree of disorder of the phosphorylated KID domain and thus investigate the intrinsic role of chain flexibility in binding kinetics. The simulations show that the "fly-casting" effect is not only due to the greater capture radii of IDPs. The coupling of folding with binding of IDPs leads to a significant reduction in binding free-energy barrier. Such a reduction accelerates the binding process. Although the greater capture radius has been regarded as the main factor in promoting the binding rate of IDPs, we found that this parameter will also lead to the slower translational diffusion of IDPs when compared with ordered proteins. As a result, the capture rate of IDPs was found to be slower than that of ordered proteins. The main origin of the faster binding for IDPs are the fewer encounter times required before the formation of the final binding complex. The roles of the interchain native contacts fraction (Q(b)) and the mass-center distance (DeltaR) as reaction coordinates are also discussed.

Structure of natural killer cell receptor KLRG1 bound to E-cadherin reveals basis for MHC-independent missing self recognition

The cytolytic activity of natural killer (NK) cells is regulated by inhibitory receptors that detect the absence of self molecules on target cells. Structural studies of missing self recognition have focused on NK receptors that bind MHC. However, NK cells also possess inhibitory receptors specific for non-MHC ligands, notably cadherins, which are downregulated in metastatic tumors. We determined the structure of killer cell lectin-like receptor G1 (KLRG1) in complex with E-cadherin. KLRG1 mediates missing self recognition by binding to a highly conserved site on classical cadherins, enabling it to monitor expression of several cadherins (E-, N-, and R-) on target cells. This site overlaps the site responsible for cell-cell adhesion but is distinct from the integrin alpha(E)beta(7) binding site. We propose that E-cadherin may coengage KLRG1 and alpha(E)beta(7) and that KLRG1 overcomes its exceptionally weak affinity for cadherins through multipoint attachment to target cells, resulting in inhibitory signaling.

OptGraft: A computational procedure for transferring a binding site onto an existing protein scaffold

One of the many challenging tasks of protein design is the introduction of a completely new function into an existing protein scaffold. In this study, we introduce a new computational procedure OptGraft for placing a novel binding pocket onto a protein structure so as its geometry is minimally perturbed. This is accomplished by introducing a two-level procedure where we first identify where are the most appropriate locations to graft the new binding pocket into the protein fold by minimizing the departure from a set of geometric restraints using mixed-integer linear optimization. On identifying the suitable locations that can accommodate the new binding pocket, CHARMM energy calculations are employed to identify what mutations in the neighboring residues, if any, are needed to ensure that the minimum energy conformation of the binding pocket conserves the desired geometry. This computational framework is benchmarked against the results available in the literature for engineering a copper binding site into thioredoxin protein. Subsequently, OptGraft is used to guide the transfer of a calcium-binding pocket from thermitase protein (PDB: 1thm) into the first domain of CD2 protein (PDB:1hng). Experimental characterization of three de novo redesigned proteins with grafted calcium-binding centers demonstrated that they all exhibit high affinities for terbium (Kd) approximately 22, 38, and 55 microM) and can selectively bind calcium over magnesium.

Size-exclusion chromatography can identify faster-associating protein complexes and evaluate design strategies

Background

We previously developed a set of rationally designed mutant MICA protein ligands for the NKG2D immunoreceptor in which MICA was mutated at residues that do not contact NKG2D. Some of these MICA mutants, predicted by RosettaDesign to be destabilized, bound NKG2D with affinities enhanced by more than an order of magnitude when evaluated by surface plasmon resonance (SPR).

Findings

Small-zone size-exclusion chromatography (SEC) detected persistent high-affinity MICA mutant-NKG2D complexes in solution as early-eluting peaks. The SEC binding assay used standard protein purification instrumentation to evaluate complex stability, qualitatively paralleled the SPR results, and successfully discriminated among complexes that differed only in on-rates. We used the SEC binding assay, along with SPR, to assess the results of a follow-up design strategy targeting the non-interfacial redesigned region. Both SEC and SPR agreed that these mutations did not enhance affinity as much as previous mutants. When the SEC binding assay was run in 1 M urea, only the highest affinity complex was detected.

Conclusion

This SEC binding assay provides a correlation with SPR results for protein complex affinities, detecting changes in complex on-rates, and tunable to lower sensitivity with 1 M urea. The SEC binding assay is complementary to other protein design evaluation methods, can be adapted to the undergraduate research laboratory, and may provide additional structural information about changes in hydrodynamic radii from elution times. Our assay allowed us to conclude that further alteration of MICA at non-contacting residues is unlikely to further enhance NKG2D affinity.

Association of MICA-TM and MICB C1_2_A microsatellite polymorphisms with tumor progression in patients with colorectal cancer

Purpose

The major histocompatibility complex class I related A (MICA) and MICB molecules are ligands of NKG2D receptors on natural killer cells, gamma/delta T cells, and CD8aß T cells that mediate host antitumor immune response. The role of MICA-TM and MICB C1_2_A alleles in patients with colorectal cancer has not yet been investigated.

Methods

We have analyzed the MICA-TM and MICB C1_2_A polymorphisms in colorectal cancer patients (n = 79) by polymerase chain reaction amplification, subsequent electrophoresis, and sequencing in comparison to a previously analyzed cohort of healthy controls (n = 306). Allele frequencies obtained for MICA-TM and MICB C1_2_A were compared to histopathological data regarding tumor invasion, disease progression, microsatellite instability, and the presence of KRAS mutations (codon 12) and analyzed for possible impact on tumor-related survival (n = 61).

Results

Allele frequencies of MICA-TM and MICB C1_2_A polymorphisms were not different in patients with colorectal cancer in comparison to normal controls. In colorectal cancer patients, MICA-TM A4 allele was directly and MICA-TM A5 allele was inversely associated with lymph node involvement and advanced UICC stages. Tumor-related survival in colorectal cancer patients was significantly reduced in the presence of the MICA-TM A4 allele (p = 0.015). In patients with microsatellite stable tumors, survival was reduced in association with the MICA-TM A4 allele (p = 0.006) and MICA-TM A9 allele (p = 0.034), but increased in patients showing the MICA-TM A5 allele (p = 0.042).

Conclusions

Specific MICA-TM alleles seem to influence tumor progression and midterm survival of patients with colorectal cancer, indicating an important role of host innate immune predisposition involving NKG2D mediated antitumor response.

A survey of the year 2007 literature on applications of isothermal titration calorimetry

Elucidation of the energetic principles of binding affinity and specificity is a central task in many branches of current sciences: biology, medicine, pharmacology, chemistry, material sciences, etc. In biomedical research, integral approaches combining structural information with in-solution biophysical data have proved to be a powerful way toward understanding the physical basis of vital cellular phenomena. Isothermal titration calorimetry (ITC) is a valuable experimental tool facilitating quantification of the thermodynamic parameters that characterize recognition processes involving biomacromolecules. The method provides access to all relevant thermodynamic information by performing a few experiments. In particular, ITC experiments allow to by-pass tedious and (rarely precise) procedures aimed at determining the changes in enthalpy and entropy upon binding by van't Hoff analysis. Notwithstanding limitations, ITC has now the reputation of being the "gold standard" and ITC data are widely used to validate theoretical predictions of thermodynamic parameters, as well as to benchmark the results of novel binding assays. In this paper, we discuss several publications from 2007 reporting ITC results. The focus is on applications in biologically oriented fields. We do not intend a comprehensive coverage of all newly accumulated information. Rather, we emphasize work which has captured our attention with originality and far-reaching analysis, or else has provided ideas for expanding the potential of the method.