Construction and analysis of a detailed three-dimensional model of the ligand binding domain of the human B cell receptor CD40

The interaction between the human B cell receptor CD40 and its ligand on T cells is critical for B cell proliferation and the regulation of humoral immune responses. CD40 is a member of the tumor necrosis factor receptor (TNFR) family. We report here the construction and analysis of a detailed three-dimensional model of the TNFR-homologous extracellular region of CD40. This study provides an example for structure-based model building in the presence of low sequence similarity. The assessment of model quality and sequence-structure compatibility is emphasized, and limitations of the model are discussed. The current CD40 model predicts structural details beyond the backbone level. Features of the CD40 ligand binding site are discussed in conjunction with the results of a previous mutagenesis study. Proteins 27:59–70 © 1997 Wiley-Liss, Inc.     

The role of polar interactions in the molecular recognition of CD40L with its receptor CD40.

CD40 Ligand (CD40L) is transiently expressed on the surface of T-cells and binds to CD40, which is expressed on the surface of B-cells. This binding event leads to the differentiation, proliferation, and isotype switching of the B-cells. The physiological importance of CD40L has been demonstrated by the fact that expression of defective CD40L protein causes an immunodeficiency state characterized by high IgM and low IgG serum levels, indicating faulty T-cell dependent B-cell activation. To understand the structural basis for CD40L/CD40 association, we have used a combination of molecular modeling, mutagenesis, and X-ray crystallography. The structure of the extracellular region of CD40L was determined by protein crystallography, while the CD40 receptor was built using homology modeling based upon a novel alignment of the TNF receptor superfamily, and using the X-ray structure of the TNF receptor as a template. The model shows that the interface of the complex is composed of charged residues, with CD40L presenting basic side chains (K143, R203, R207), and CD40 presenting acidic side chains (D84, E114, E117). These residues were studied experimentally through site-directed mutagenesis, and also theoretically using electrostatic calculations with the program Delphi. The mutagenesis data explored the role of the charged residues in both CD40L and CD40 by switching to Ala (K143A, R203A, R207A of CD40L, and E74A, D84A, E114A, E117A of CD40), charge reversal (K143E, R203E, R207E of CD40L, and D84R, E114R, E117R of CD40), mutation to a polar residue (K143N, R207N, R207Q of CD40L, and D84N, E117N of CD40), and for the basic side chains in CD40L, isosteric substitution to a hydrophobic side chain (R203M, R207M). All the charge-reversal mutants and the majority of the Met and Ala substitutions led to loss of binding, suggesting that charged interactions stabilize the complex. This was supported by the Delphi calculations which confirmed that the CD40/CD40L residue pairs E74-R203, D84-R207, and E117-R207 had a net stabilizing effect on the complex. However, the substitution of hydrophilic side chains at several of the positions was tolerated, which suggests that although charged interactions stabilize the complex, charge per se is not crucial at all positions. Finally, we compared the electrostatic surface of TNF/TNFR with CD40L/CD40 and have identified a set of polar interactions surrounded by a wall of hydrophobic residues that appear to be similar but inverted between the two complexes.     

Immunoglobulin fold characteristics of B7-1 (CD80) and B7-2 (CD86).

B7-1 and B7-2 are expressed on antigen-presenting cells and bind to the CD28 and CTLA-4 receptors on T cells. These interactions trigger a costimulatory pathway that is essential for T-cell activation. B7-1 and B7-2 are members of the immunoglobulin superfamily (IgSF) and, despite sharing common function, have only limited sequence similarity. The B7-1 extracellular region was previously subdivided into 2 IgSF domains, an N-terminal V(ariable)-like domain, followed by a C(onstant)-like domain. We recently reported that the V-like domains of B7-1 and B7-2 share some significant sequence similarities with 3 major histocompatibility complex (MHC)-encoded members of the IgSF. We have now applied inverse folding methodology to assess the compatibility of the B7-1 and B7-2 extracellular region sequences with currently available 3-dimensional structures. In these calculations, the sequences of the N-terminal (V-like) domains in B7-1 and B7-2 were not compatible with known structures, including the IgSF V-set. In contrast, the sequences of the C-like domains were compatible with IgSF C-set structures and were best recognized by the beta 2-microglobulin (beta 2m) domain of MHC Class I. A sequence comparison of the C-like domains in the B7 molecules showed that 11 of 17 rigorously conserved residues in B7-1 and B7-2 are not IgSF C-1 set consensus residues. When mapped onto the corresponding positions of the beta 2m structure, the conserved residues in B7 cluster on the surface, where they may interact with the B7 V-like domain or other molecules.     

CD6 recognizes the neural adhesion molecule BEN.

CD6 and its ligand activated leukocyte cell adhesion molecule (ALCAM, CD166) have been detected on various immune cells and in the brain. CD6-ligand interactions have been implicated in the regulation of T cell function. ALCAM shares the same extracellular domain organization and significant sequence homology with the chicken neural adhesion molecule BEN. Although ALCAM's CD6 binding site is only partially conserved in BEN, CD6 specifically binds BEN, albeit with approximately 10-fold lower avidity than ALCAM. Differences in binding avidity are not detected when ALCAM and BEN fusion proteins containing the full-length extracellular regions are tested. Homotypic interactions between full-length forms are likely to account for these observations. The identified cross-species interaction between CD6 and BEN suggests that CD6-ligand interactions are highly conserved.     

The counterreceptor binding site of human CD2 exhibits an extended surface patch with multiple conformations fluctuating with millisecond to microsecond motions.

We have used 15N NMR relaxation experiments to probe, for the glycosylated human CD2 adhesion domain, the overall molecular motion, as well as very fast nanosecond-picosecond (ns-ps) and slow millisecond-microsecond (ms-microsecond) internal motions. Using a novel analysis method that considers all residues, we obtained a correlation time for the overall motion of 9.5 +/- 0.3 ns. Surprisingly, we found a large contiguous patch of residues in the counterreceptor (CD58) binding site of human CD2 exhibiting slow conformational exchange motions (ms-microsecond). On the other hand, almost none of the residues of the CD58 binding side display fast (ns-ps) internal motions of amplitudes larger than what is seen for well-ordered regions of the structure. Residues close to the N-glycosylation site, and the first N-acetylglucosamine of the high mannose glycan are as rigid as the protein core. Residues conserved in the immunoglobulin superfamily V-set domain are generally very rigid.     

Three-dimensional solution structure and conformational plasticity of the N-terminal scavenger receptor cysteine-rich domain of human CD5

The lymphocyte receptor CD5 influences cell activation by modifying the strength of the intracellular response initiated by antigen engagement. Regulation through CD5 involves the interaction of one or more of its three scavenger receptor cysteine-rich domains present in the extracellular region. Here, we present the 3D solution structure of a non-glycosylated double mutant of the N-terminal domain of human CD5 expressed in Escherichia coli (eCD5d1m), which has enhanced solubility compared to the non-glycosylated wild-type (eCD5d1). In common with a glycosylated form expressed in Pichia pastoris, the [¹⁵N,¹H]-correlation spectra of both eCD5d1 and eCD5d1m exhibit non-uniform temperature-dependent signal intensities, indicating extensive conformational fluctuations on the micro-millisecond timescale. Although approximately one half of the signals expected for the domain are absent at 298 K, essentially complete resonance assignments and a solution structure could be obtained at 318 K. Because of the sparse nature of the experimental restraint data and the potentially important contribution of conformational exchange to the nuclear Overhauser effect peak intensity, we applied inferential structure determination to calculate the eCD5d1m structure. The inferential structure determination ensemble has similar features to that obtained by traditional simulated annealing methods, but displays superior definition and structural quality. The eCD5d1m structure is similar to other members of the scavenger receptor cysteine-rich superfamily, but the position of the lone α helix differs due to interactions with the unique N-terminal region of the domain. The availability of an experimentally tractable form of CD5d1, together with its 3D structure, provides new tools for further investigation of its function within intact CD5.     

Molecular and Functional Characterization of the Four-Transmembrane Molecule L6 in Epidermal Keratinocytes

In normal human epidermal keratinocytes (NHEK) proteolytic detachment from the substrate induces a complex activation cascade including expression of new proteins, morphological alterations, and the onset of migration for epidermal regeneration. By subtractive cloning we have shown that L6, a four-transmembrane protein, is newly expressed after proteolytic keratinocyte detachment. In this study, we have generated a novel anti-L6 antibody (clone HD-pKe#104-1.1) and investigated L6 expression regulation in vitro and in vivo as well as L6 function in keratinocyte migration. Dispase-mediated detachment induced L6 expression in NHEK at the mRNA and protein level. Immunohistology of skin biopsies displayed a strong expression of L6 in follicular epidermis and epidermolytic lesions of autoimmune bullous dermatoses (bullous pemphigoid, pemphigus vulgaris), but not in normal interfollicular epidermis. In contrast to normal keratinocytes, HaCaT cells showed constitutive L6 expression, indicating a constitutively active phenotype. After artificial wounding of confluent HaCaT cultures, anti-L6 antibody strongly impaired cell migration velocity and migratory reepithelization of the defect, indicating L6 involvement in keratinocyte migration. These findings suggest that L6 is an important activation-dependent regulator of keratinocyte function and epidermal tissue regeneration.     

A proposed structural model of domain 1 of fasciclin III neural cell adhesion protein based on an inverse folding algorithm.

Fasciclin III is an integral membrane protein expressed on a subset of axons in the developing Drosophila nervous system. It consists of an intracellular domain, a transmembrane region, and an extracellular region composed of three domains, each predicted to form an immunoglobulin-like fold. The most N-terminal of these domains is expected to be important in mediating cell-cell recognition events during nervous system development. To learn more about the structure/function relationships in this cellular recognition molecule, a model structure of this domain was built. A sequence-to-structure alignment algorithm was used to align the protein sequence of the fasciclin III first domain to the immunoglobulin McPC603 structure. Based on this alignment, a model of the domain was built using standard homology modeling techniques. Side-chain conformations were automatically modeled using a rotamer search algorithm and the model was minimized to relax atomic overlaps. The resulting model is compact and has chemical characteristics consistent with related globular protein structures. This model is a de novo test of the sequence-to-structure alignment algorithm and is currently being used as the basis for mutagenesis experiments to discern the parts of the fasciclin III protein that are necessary for homophilic molecular recognition in the developing Drosophila nervous system.     

Variants of the CD40 ligand gene are not associated with increased susceptibility to tuberculosis in West Africa

Evidence for linkage between tuberculosis and human chromosomal region Xq26 has previously been described. The costimulatory molecule CD40 ligand, encoded by TNFSF5 and located at Xq26.3, is a promising positional candidate. Interactions between CD40 ligand and CD40 are involved in the development of humoral- and cell-mediated immunity, as well as the activation of macrophages, which are the primary host and effector cells for Mycobacterium tuberculosis. We hypothesised that common variation within TNFSF5 might affect susceptibility to tuberculosis disease and, thus, might be responsible for the observed linkage to Xq26. Sequencing 32 chromosomes from a Gambian population identified nine common polymorphisms within the coding, 3 and 5 regulatory sequences of the gene. Six single nucleotide polymorphisms (SNPs) and a 3 microsatellite were genotyped in 121 tuberculosis patients and their available parents. No association with tuberculosis was detected for these variants using a transmission disequilibrium test, although one SNP at –726 showed some evidence of association in males. This finding, however, did not replicate in a separate case control study of over 1,200 West African individuals. We conclude that common genetic variation in TNFSF5 is not likely to affect tuberculosis susceptibility in West Africa and the linkage observed in this region is not due to variation in TNFSF5.Sadly, Professor Steve Bennett passed away in March 2003     

Characterization of a comparative model of the extracellular domain of the epidermal growth factor receptor

The Epidermal Growth Factor (EGF) receptor is a tyrosine kinase that mediates the biological effects of ligands such as EGF and transforming growth factor alpha. An understanding of the molecular basis of its action has been hindered by a lack of structural and mutational data on the receptor. We have constructed comparative models of the four extracellular domains of the EGF receptor that are based on the structure of the first three domains of the insulin-like growth factor-1 (IGF-1) receptor. The first and third domains of the EGF receptor, L1 and L2, are right-handed beta helices. The second and fourth domains of the EGF receptor, S1 and S2, consist of the modules held together by disulfide bonds, which, except for the first module of the S1 domain, form rod-like structures. The arrangement of the L1 and S1 domains of the model are similar to that of the first two domains of the IGF-1 receptor, whereas that of the L2 and S2 domains appear to be significantly different. Using the EGF receptor model and limited information from the literature, we have proposed a number of regions that may be involved in the functioning of the receptor. In particular, the faces containing the large beta sheets in the L1 and L2 domains have been suggested to be involved with ligand binding of EGF to its receptor.     

Binding site on human immunoglobulin G for the affinity ligand HWRGWV

Affinity ligand HWRGWV has demonstrated the ability to isolate human immunoglobulin G (hIgG) from mammalian cell culture media. The ligand specifically binds hIgG through its Fc portion. This work shows that deglycosylation of hIgG has no influence on its binding to the HWRGWV ligand and the ligand does not compete with Protein A or Protein G in binding hIgG. It is suggested by the mass spectrometry (MS) data and docking simulation that HWRGWV binds to the pFc portion of hIgG and interacts with the amino acids in the loop Ser383–Asn389 (SNGQPEN) located in the C_H3 domain. Subsequent modeling has suggested a possible three‐dimensional minimized solution structure for the interaction of hIgG and the HWRGWV ligand. The results support the fact that a peptide as small as a hexamer can have specific interactions with large proteins such as hIgG. Copyright © 2009 John Wiley & Sons, Ltd.     

Molecular dynamics simulation of a human thiopurine S-methyltransferase complexed with 6-mercaptopurine model

Human thiopurine S-methyltransferase (TPMT) is an essential protein in 6-mercaptopurine (6MP) drug metabolism. To understand the pharmacogenetics of TPMT and 6MP, X-ray co-crystal structures of TPMT complexes with S-adenosyl-L-methionine (AdoMet) and 6MP are required. However, the co-crystal structure of this complex has not been reported because 6MP is poorly water soluble. We used molecular dynamics (MD) simulation to predict the structure of the complex of human TPMT-AdoHcy(CH₂)6MP, where the sulfur atoms of AdoHcy and 6MP were linked by a CH₂ group. After 1300 picoseconds of MD simulation, the trajectory showed that 6MP was stabilized in the TPMT active site by formation of non-bonded interactions between 6MP and Phe40, Pro196 and Arg226 side chains of TPMT. The intersulfur distance between AdoHcy and 6MP as well as the binding modes and the interactions of our TPMT-AdoHcy model are consistent with those observed in the X-ray crystal structure of murine TPMT-AdoHcy-6MP complex. The predicted binding modes of AdoHcy and 6MP in our model are consistent with those observed in murine TPMT X-ray crystal structures, which provides structural insights into the interactions of TPMT, AdoHcy, and 6MP at the atomic level and may be used as a starting point for further study of thiopurine drug pharmacogenetics.     

Antagonist properties of monoclonal antibodies targeting human CD28

Antagonist antibodies targeting CD28 have been proposed as an alternative to the use of CD80/86 antagonists to modulate T cell responses in autoimmunity and transplantation. Advantages would be the blockade of CD28-mediated co-stimulatory signals without impeding the co-inhibitory signals dependent on CD80 interactions with CTLA-4 and PD-L1 that are important for the control of immune responses and for the function of regulatory T cells. Anti-CD28 antibodies are candidate antagonists only if they prevent access to the CD80/86 ligands without simultaneously stimulating CD28 itself, a process that is believed to depend on receptor multimerization. In this study, we evaluated the impact of different formats of a potentially antagonist anti-human CD28 antibody on T cell activation. In particular, we examined the role of valency and of the presence of an Fc domain, two components that might affect receptor multimerization either directly or in the presence of accessory cells expressing Fc receptors. Among monovalent (Fab’, scFv), divalent (Fab’2), monovalent-Fc (Fv-Fc) and divalent-Fc (IgG) formats, only the monovalent formats showed consistent absence of induced CD28 multimerization and absence of associated activation of phosphoinositol-3-kinase, and clear antagonist properties in T cell stimulation assays. In contrast, divalent antibodies showed agonist properties that resulted in cell proliferation and cytokine release in an Fc-independent manner. Conjugation of monovalent antibodies with polyethylene glycol, α-1-antitrypsin or an Fc domain significantly extended their in vivo half-life without modifying their antagonist properties. In conclusion, these data indicate that monovalency is mandatory for maintaining the antagonistic activity of anti-CD28 monoclonal antibodies.     

Extending the B7 (CD80) gene family.

B7-1 and B7-2 are members of the immunoglobulin superfamily (IgSF) and important regulators of T cell-mediated immune responses. Despite sharing only limited sequence identity, B7-1 and B7-2 bind common receptors, CD28 and CTLA-4, on T cells and have similar functional properties. We have found that the extracellular V (ariable)-like domains of B7-1 and B7-2 share significant sequence similarities with 3 major histocompatibility complex (MHC)-encoded members of the IgSF: butyrophilin, myelin/oligodendrocyte glycoprotein, and the chicken MHC molecule, B-G. This raises the question whether there is an evolutionary link between the MHC, which encodes molecules regulating the antigen specificity of T lymphocyte responses, and B7 molecules, which co-stimulate these responses in antigen-nonspecific fashion.     

Homology modeling and molecular dynamics simulations of the glycine receptor ligand binding domain

We present a homology based model of the ligand binding domain (LBD) of the homopentameric alpha1 glycine receptor (GlyR). The model is based on multiple sequence alignment with other members of the nicotinicoid ligand gated ion channel superfamily and two homologous acetylcholine binding proteins (AChBP) from the freshwater (Lymnaea stagnalis) and saltwater (Aplysia californica) snails with known high resolution structure. Using two template proteins with known structure to model three dimensional structure of a target protein is especially advantageous for sequences with low homology as in the case presented in this paper. The final model was cross-validated by critical evaluation of experimental and published mutagenesis, functional and other biochemical studies. In addition, a complex structure with strychnine antagonist in the putative binding site is proposed based on docking simulation using Autodock program. Molecular dynamics (MD) simulations with simulated annealing protocol are reported on the proposed LBD of GlyR, which is stable in 5 ns simulation in water, as well as for a deformed LBD structure modeled on the corresponding domain determined in low-resolution cryomicroscopy structure of the alpha subunit of the full-length acetylcholine receptor (AChR). Our simulations demonstrate that the beta-sandwich central core of the protein monomer is fairly rigid in the simulations and resistant to deformations in water.     

Three-dimensional Model Structure for the Extracellular Domains of Fibroblast Growth Factor Receptor - 1 (FGFR-1)

We have developed a model for the two immunoglobulin-like extracellular domains DII and DIII of the FGF receptor 1 (FGFR-1), giving a special attention to the determination of the appropriate Ig set. The DII domain was aligned with the C-terminal domain of myosin light chain kinase (telokin) of the I set, and the DIII domain with the variable domain of the Bence-Jones immunoglobulin of the V set. Two assemblies, corresponding to different propositions for the domains relative orientation, have been refined and compared.Electronic Supplementary Material available.     

Molecular modeling of protein function regions

Experimental protein structures often provide extensive insight into the mode and specificity of small molecule binding, and this information is useful for understanding protein function and for the design of drugs. We have performed an analysis of the reliability with which ligand-binding information can be deduced from computer model structures, as opposed to experimentally derived ones. Models produced as part of the CASP experiments are used. The accuracy of contacts between protein model atoms and experimentally determined ligand atom positions is the main criterion. Only comparative models are included (i.e., models based on a sequence relationship between the protein of interest and a known structure). We find that, as expected, contact errors increase with decreasing sequence identity used as a basis for modeling. Analysis of the causes of errors shows that sequence alignment errors between model and experimental template have the most deleterious effect. In general, good, but not perfect, insight into ligand binding can be obtained from models based on a sequence relationship, providing there are no alignment errors in the model. The results support a structural genomics strategy based on experimental sampling of structure space so that all protein domains can be modeled on the basis of 30% or higher sequence identity.     

Molecular insights into the binding selectivity of a synthetic ligand DAAG to Fc fragment of IgG

Affinity chromatography with synthetic ligands has been focused as the potential alternative to protein A‐based chromatography for antibody capture because of its comparable selectivity and efficiency. Better understanding on the molecular interactions between synthetic ligand and antibody is crucial for improving and designing novel ligands. In this work, the molecular interaction mechanism between Fc fragment of IgG and a synthetic ligand (DAAG) was studied with molecular docking and dynamics simulation. The docking results on the consensus binding site (CBS) indicated that DAAG could bind to the CBS with the favorable orientation like a tripod for the top‐ranked binding complexes. The ligand‐Fc fragment complexes were then tested by molecular dynamics simulation at neutral condition (pH 7.0) for 10 ns. The results indicated that the binding of DAAG on the CBS of Fc fragment was achieved by the multimodal interactions, combining the hydrophobic interaction, electrostatic interaction, hydrogen bond, and so on. It was also found that multiple secondary interactions endowed DAAG with an excellent selectivity to Fc fragment. In addition, molecular dynamics simulation conducted at acidic condition (pH 3.0) showed that the departure of DAAG ligand from the surface of Fc fragment was the result of reduced interaction energies. The binding modes between DAAG and CBS not only shed light on the molecular mechanisms of DAAG for antibody purification but also provide useful information for the improvement of ligand design. Copyright © 2014 John Wiley & Sons, Ltd.     

Structure-function analysis of human IL-6: identification of two distinct regions that are important for receptor binding.

Interleukin-6 (IL-6) is a multifunctional cytokine that plays an important role in host defense. It has been predicted that IL-6 may fold as a 4 alpha-helix bundle structure with up-up-down-down topology. Despite a high degree of sequence similarity (42%) the human and mouse IL-6 polypeptides display distinct species-specific activities. Although human IL-6 (hIL-6) is active in both human and mouse cell assays, mouse IL-6 (mIL-6) is not active on human cells. Previously, we demonstrated that the 5 C-terminal residues of mIL-6 are important for activity, conformation, and stability (Ward LD et al., 1993, Protein Sci 2:1472-1481). To further probe the structure-function relationship of this cytokine, we have constructed several human/mouse IL-6 hybrid molecules. Restriction endonuclease sites were introduced and used to ligate the human and mouse sequences at junction points situated at Leu-62 (Lys-65 in mIL-6) in the putative connecting loop AB between helices A and B, at Arg-113 (Val-117 in mIL-6) at the N-terminal end of helix C, at Lys-150 (Asp-152 in mIL-6) in the connecting loop CD between helices C and D, and at Leu-178 (Thr-180 in mIL-6) in helix D. Hybrid molecules consisting of various combinations of these fragments were constructed, expressed, and purified to homogeneity. The conformational integrity of the IL-6 hybrids was assessed by far-UV CD. Analysis of their biological activity in a human bioassay (using the HepG2 cell line), a mouse bioassay (using the 7TD1 cell line), and receptor binding properties indicates that at least 2 regions of hIL-6, residues 178-184 in helix D and residues 63-113 in the region incorporating part of the putative connecting loop AB through to the beginning of helix C, are critical for efficient binding to the human IL-6 receptor. For human IL-6, it would appear that interactions between residues Ala-180, Leu-181, and Met-184 and residues in the N-terminal region may be critical for maintaining the structure of the molecule; replacement of these residues with the corresponding 3 residues in mouse IL-6 correlated with a significant loss of alpha-helical content and a 200-fold reduction in activity in the mouse bioassay. A homology model of mIL-6 based on the X-ray structure of human granulocyte colony-stimulating factor is presented.