Homology modeling of human Toll‐like receptors TLR7, 8, and 9 ligand‐binding domains

Toll‐like receptors (TLRs) play a key role in the innate immune system. The TLR7, 8, and 9 compose a family of intracellularly localized TLRs that signal in response to pathogen‐derived nucleic acids. So far, there are no crystallographic structures for TLR7, 8, and 9. For this reason, their ligand‐binding mechanisms are poorly understood. To enable first predictions of the receptor–ligand interaction sites, we developed three‐dimensional structures for the leucine‐rich repeat ectodomains of human TLR7, 8, and 9 based on homology modeling. To achieve a high sequence similarity between targets and templates, structural segments from all known TLR ectodomain structures (human TLR1/2/3/4 and mouse TLR3/4) were used as candidate templates for the modeling. The resulting models support previously reported essential ligand‐binding residues. They also provide a basis to identify three potential receptor dimerization mechanisms. Additionally, potential ligand‐binding residues are identified using combined procedures. We suggest further investigations of these residues through mutation experiments. Our modeling approach can be extended to other members of the TLR family or other repetitive proteins.     

DNA binding to proteolytically activated TLR9 is sequence-independent and enhanced by DNA curvature

Toll-like receptor 9 (TLR9) recognizes microbial DNA in endolysosomal compartments. The ectodomain of TLR9 must be proteolytically cleaved by endosomal proteases to produce the active receptor capable of inducing an innate immune signal. We show that the cleaved TLR9 ectodomain is a monomer in solution and that DNA ligands with phosphodiester backbones induce TLR9 dimerization in a sequence-independent manner. Ligands with phosphorothioate (PS) backbones induce the formation of large TLR9-DNA aggregates, possibly due to the propensity of PS ligands to self-associate. DNA curvature-inducing proteins including high-mobility group box 1 and histones H2A and H2B significantly enhance TLR9 binding, suggesting that TLR9 preferentially recognizes curved DNA backbones. Our work sheds light on the molecular mechanism of TLR9 activation by endogenous protein-nucleic acid complexes, which are associated with autoimmune diseases including systemic lupus erythematosus.     

Recognition highlights: Toll-like receptors

As part of our first line of defense against invading organisms, Toll-like receptors recognize bacterial and viral molecules, and launch an inflammatory response.     

C‐terminal deletion of leucine‐rich repeats from YopM reveals a heterogeneous distribution of stability in a cooperatively folded protein

Terminal deletions of units from α‐helical repeat proteins have provided insight into the physical origins of their cooperativity. To test if the same principles governing cooperativity apply to β‐sheet‐containing repeat proteins, we have created a series of C‐terminal deletion constructs from a large leucine‐rich repeat (LRR) protein, YopM. We have examined the structure and stability of the resulting deletion constructs by a combination of solution spectroscopy, equilibrium denaturation studies, and limited proteolysis. Surprisingly, a high degree of nonuniformity was found in the stability distribution of YopM. Unlike previously studied repeat proteins, we identified several key LRR that on deletion disrupt nearby structure, at distances as far away as up to three repeats, in YopM. This partial unfolding model is supported by limited proteolysis studies and by point substitution in repeats predicted to be disordered as a result of deletion of adjacent repeats. We show that key internal‐ and terminal‐caps must be present to maintain the structural integrity in adjacent regions (roughly four LRRs long) of decreased stability. The finding that full‐length YopM maintains a high level of cooperativity in equilibrium unfolding underscores the importance of interfacial interactions in stabilizing locally unstable regions of structure.     

Leucine‐rich repeat receptor‐like kinase II phylogenetics reveals five main clades throughout the plant kingdom

Receptor‐like kinases (RLKs) represent the largest group of cell surface receptors in plants. The monophyletic leucine‐rich repeat (LRR)‐RLK subfamily II is considered to contain the somatic embryogenesis receptor kinases (SERKs) and NSP‐interacting kinases known to be involved in developmental processes and cellular immunity in plants. There are only a few published studies on the phylogenetics of LRR‐RLKII; unfortunately these suffer from poor taxon/gene sampling. Hence, it is not clear how many and what main clades this family contains, let alone what structure–function relationships exist. We used 1342 protein sequences annotated as ‘SERK’ and ‘SERK‐like’ plus related sequences in order to estimate phylogeny within the LRR‐RLKII clade, using the nematode protein kinase Pelle as an outgroup. We reconstruct five main clades (LRR‐RLKII 1–5), in each of which the main pattern of land plant relationships re‐occurs, confirming previous hypotheses that duplication events happened in this gene subfamily prior to divergence among land plant lineages. We show that domain structures and intron–exon boundaries within the five clades are well conserved in evolution. Furthermore, phylogenetic patterns based on the separate LRR and kinase parts of LRR‐RLKs are incongruent: whereas the LRR part supports a LRR‐RLKII 2/3 sister group relationship, the kinase part supports clades 1/2. We infer that the kinase part includes few ‘radical’ amino acid changes compared with the LRR part. Finally, our results confirm that amino acids involved in each LRR‐RLKII–receptor complex interaction are located at N‐capping residues, and that the short amino acid motifs of this interaction domain are highly conserved throughout evolution within the five LRR‐RLKII clades.     

A mutational analysis of the cytosolic domain of the tomato Cf‐9 disease‐resistance protein shows that membrane‐proximal residues are important for Avr9‐dependent necrosis

The tomato Cf‐9 gene encodes a membrane‐anchored glycoprotein that imparts race‐specific resistance against the tomato leaf mould fungus Cladosporium fulvum in response to the avirulence protein Avr9. Although the N‐terminal half of the extracellular leucine‐rich repeat (eLRR) domain of the Cf‐9 protein determines its specificity for Avr9, the C‐terminal half, including its small cytosolic domain, is postulated to be involved in signalling. The cytosolic domain of Cf‐9 carries several residues that are potential sites for ubiquitinylation or phosphorylation, or signals for endocytic uptake. A targeted mutagenesis approach was employed to investigate the roles of these residues and cellular processes in Avr9‐dependent necrosis triggered by Cf‐9. Our results indicate that the membrane‐proximal region of the cytosolic domain of Cf‐9 plays an important role in Cf‐9‐mediated necrosis, and two amino acids within this region, a threonine (T835) and a proline (P838), are particularly important for Cf‐9 function. An alanine mutation of T835 had no effect on Cf‐9 function, but an aspartic acid mutation, which mimics phosphorylation, reduced Cf‐9 function. We therefore postulate that phosphorylation/de‐phosphorylation of T835 could act as a molecular switch to determine whether Cf‐9 is in a primed or inactive state. Yeast two‐hybrid analysis was used to show that the cytosolic domain of Cf‐9 interacts with the cytosolic domain of tomato VAP27. This interaction could be disrupted by an alanine mutation of P838, whereas interaction with CITRX remained unaffected. We therefore postulate that a proline‐induced kink in the membrane‐proximal region of the cytosolic domain of Cf‐9 may be important for interaction with VAP27, which may, in turn, be important for Cf‐9 function.     

Neutralization of pro‐inflammatory monocytes by targeting TLR2 dimerization ameliorates colitis

Monocytes have emerged as critical driving force of acute inflammation. Here, we show that inhibition of Toll‐like receptor 2(TLR2) dimerization by a TLR2 transmembrane peptide (TLR2‐p) ameliorated DSS‐induced colitis by interfering specifically with the activation of Ly6C⁺ monocytes without affecting their recruitment to the colon. We report that TLR2‐p directly interacts with TLR2 within the membrane, leading to inhibition of TLR2–TLR6/1 assembly induced by natural ligands. This was associated with decreased levels of extracellular signal‐regulated kinases (ERK) signaling and reduced secretion of pro‐inflammatory cytokines, such as interleukin (IL)‐6, IL‐23, IL‐12, and IL‐1β. Altogether, our study provides insights into the essential role of TLR2 dimerization in the activation of pathogenic pro‐inflammatory Ly6C^hi monocytes and suggests that inhibition of this aggregation by TLR2‐p might have therapeutic potential in the treatment of acute gut inflammation.     

Comprehensive modeling and functional analysis of Toll‐like receptor ligand‐recognition domains

Toll‐like receptors (TLRs) are innate immune pattern‐recognition receptors endowed with the capacity to detect microbial pathogens based on pathogen‐associated molecular patterns. The understanding of the molecular principles of ligand recognition by TLRs has been greatly accelerated by recent structural information, in particular the crystal structures of leucine‐rich repeat‐containing ectodomains of TLR2, 3, and 4 in complex with their cognate ligands. Unfortunately, for other family members such as TLR7, 8, and 9, no experimental structural information is currently available. Methods such as X‐ray crystallography or nuclear magnetic resonance are not applicable to all proteins. Homology modeling in combination with molecular dynamics may provide a straightforward yet powerful alternative to obtain structural information in the absence of experimental (structural) data, provided that the generated three‐dimensional models adequately approximate what is found in nature. Here, we report the development of modeling procedures tailored to the structural analysis of the extracellular domains of TLRs. We comprehensively compared secondary structure, torsion angles, accessibility for glycosylation, surface charge, and solvent accessibility between published crystal structures and independently built TLR2, 3, and 4 homology models. Finding that models and crystal structures were in good agreement, we extended our modeling approach to the remaining members of the TLR family from human and mouse, including TLR7, 8, and 9.     

Direct measurement of the interaction force between immunostimulatory CpG-DNA and TLR9 fusion protein

The specific interaction between human Toll-like receptor 9 (TLR9)-ectodomain (ECD)-fusion protein and immunostimulatory CpG-DNA was measured using force spectroscopy. Flexible tethers were used between receptors and surface as well as between DNA and atomic force microscope tip to make efficient recognition studies possible. The molecular recognition forces detected are in the range of 50 to 150 ± 20 pN at the used force-loading rates, and the molecular interaction probability was much reduced when the receptors were blocked with free CpG-DNA. A linear increase of the unbinding force with the logarithm of the loading rate was found over the range 0.1 to 30 nN/s. This indicates a single potential barrier characterizing the energy landscape and no intermediate state for the unbinding pathway of CpG-DNA from the TLR9-ECD. Two important kinetic parameters for CpG-DNA interaction with TLR9-ECD were determined from the force-loading rate dependency: an off-rate of k(off) = 0.14 ± 0.10 s(-1) and a binding interaction length of x(β) = 0.30 ± 0.03 nm, which are consistent with literature values. Various models for the molecular interaction of this innate immune receptor binding to CpG-DNA are discussed.     

Metal allergens nickel and cobalt facilitate TLR4 homodimerization independently of MD2

Development of contact allergy requires cooperation of adaptive and innate immunity. Ni²⁺ stimulates innate immunity via TLR4/MD2, the bacterial LPS receptor. This likely involves receptor dimerization, but direct proof is pending and it is unclear if related haptens share this mechanism. We reveal Co²⁺ as second metal stimulating TLR4 and confirm necessity of H456/H458 therein. Experiments with a new TLR4 dimerization mutant established dimerization as a mechanism of metal‐ and LPS‐induced TLR4 activation. Yet, in interaction studies only LPS‐ but not metal‐induced dimerization required MD2. Consistently, soluble TLR4 expressed without MD2 inhibited metal‐ but not LPS‐induced responses, opening new therapeutic perspectives.     

Crystal structure of the N‐terminal methyltransferase‐like domain of anamorsin

Anamorsin is a recently identified molecule that inhibits apoptosis during hematopoiesis. It contains an N‐terminal methyltransferase‐like domain and a C‐terminal Fe‐S cluster motif. Not much is known about the function of the protein. To better understand the function of anamorsin, we have solved the crystal structure of the N‐terminal domain at 1.8 Å resolution. Although the overall structure resembles a typical S‐adenosylmethionine (SAM) dependent methyltransferase fold, it lacks one α‐helix and one β‐strand. As a result, the N‐terminal domain as well as the full‐length anamorsin did not show S‐adenosyl‐l ‐methionine (AdoMet) dependent methyltransferase activity. Structural comparisons with known AdoMet dependent methyltransferases reveals subtle differences in the SAM binding pocket that preclude the N‐terminal domain from binding to AdoMet. The N‐terminal methyltransferase‐like domain of anamorsin probably functions as a structural scaffold to inhibit methyl transfers by out‐competing other AdoMet dependant methyltransferases or acts as bait for protein–protein interactions.Proteins 2014; 82:1066–1071. © 2013 Wiley Periodicals, Inc.     

Study of Q224K,V152G double mutation in bean PGIP2, an LRR protein for plant defense—An in silico approach

Polygalacturonase inhibiting proteins (PGIPs) are leucine‐rich repeat (LRR) proteins from plants that are organized into multigene families. They act as specific inhibitors against Polygalacturonases (PGs) from phytopathogens and share high sequence identity within species. We performed in silico mutation (Q224K and V152G) in PGIP2 from Phaseolus vulgaris to corresponding residues of another member, PGIP1. This mutation is known to cause 100% loss of inhibition against the PG of fungus Fusarium phyllophilum (Fp). A comparative analysis between PGIP2 and the double mutant, using 50 ns molecular dynamics simulations explored structural difference affecting PG binding properties. Simulations revealed that the mutation at 224, strains this residue which acts as a lock for the PGIP‐PG complex through main chain H‐bond. Changes in secondary structural elements and strain in the bend region along the convex face of the solenoidal protein affected the flexibility of the mutant protein. At the concave interacting face of the mutant, subtle changes in the sidechain behavior of the PG‐binding residues occurred in a concerted manner revealing flipping of aromatic rings to be crucial to avoid steric clash with FpPG in PGIP2. Docking PGIP2 and the mutant protein individually to FpPG illustrated the inability of the latter to inhibit FpPG leaving its active site free. Our study demonstrates that the effect of mutation affects the flexibility of the protein along the convex face, while binding specificity is altered through the concave face imparting minimal change in the typical structure supported by the LRRs. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

The evolutionary characteristics and structural biology of Gallus toll‐like receptor 21

Toll‐like receptors (TLRs) are an important part of the innate immune system, acting as a first line of defense against many invading pathogens. The ligand known to bind Gallus toll‐like receptor 21 (gTLR21) is the unmethylated cytosine phosphate guanine dideoxy nucleotide motif; however, the evolutionary characteristics and structural biology of gTLR21 are poorly elaborated. Our results suggest that gTLR21 is phylogenetically and evolutionarily related to the TLR11 family and is perhaps a close ortholog of the Mus TLR13. Structural biology of homology modeling of the gTLR21 ectodomain structure suggests that it has no Z‐loop like that seen in Mus TLR9. The cytosolic toll‐IL‐1 receptor region of gTLR21 contains a central 4‐stranded parallel β‐sheet (βA‐βD) surrounded by 5 α‐helices (αA‐αE) on both sides, a highly conserved structure also seen in other TLRs. Molecular docking analysis reveals that the gTLR21 ectodomain has the potential to distinguish between different ligands. Homodimer analysis results also suggest that Phe842 and Pro844 of the BB loop and Cys876 of the αC helix in gTLR21 are conserved in other cytosolic toll‐IL‐1 receptor domains of other TLRs and may contribute to the docking of homodimers. Our study on the evolutionary characteristics and structural biology of gTLR21 reveals that the molecule may have a broader role to play in innate immune system; however, further experimental validation is required to confirm our findings.     

A threshold level of TLR9 mRNA predicts cellular responsiveness to CpG-ODN in haematological and non-haematological tumour cell lines

The human toll like receptor 9 (TLR9) detects differences between microbial and host DNA, based on unmethylated deoxycytidyl deoxyguanosine dinucleotide (CpG) motifs, leading to activation of both innate and adaptive immune mechanisms. The synthetic TLR9 agonist, CpG-ODN, can substitute for microbial DNA in these responses, and is in clinical trials as an immunomodulatory agent in diseases as diverse as infections, cancer and allergic disorders. Human TLR9 is expressed on cells of haematopoietic origin (principally plasmacytoid dendritic cells and B cells), but has also been described as being expressed on a number of other cell types. In order to clarify the expression and function of TLR9 in a range of cells of both haematopoietic and non-haematopoietic origin, we investigated the level of expression of TLR9 mRNA, and the ability of the cells to respond to CpG-ODN by upregulation of cell surface markers, cytokine production, cellular proliferation and activation of NFκB. Our data show that the cellular response to CpG-ODN depended on a threshold level of expression of TLR9. TLR9 was widely expressed amongst B cell tumours (with the exception of myeloma cell lines), but we did not find either threshold levels of expression of TLR9 or responses to CpG-ODN in several myeloma or myeloid tumour cell lines or any non-haematological tumour cell lines tested in our study. TLR9-positive cells varied significantly in their responses to CpG-ODN, and the level of TLR9 expression beyond the threshold did not correlate with the magnitude of the response to CpG-ODN. Finally, CpG-ODN induced NFκB activation and increased cellular proliferation in Hek293 cells that had been stably transfected with hTLR9, but did not affect the expression of surface markers or synthesis of IL-6, IL-10 or TNF-α. Thus both haematological and non-haematological cells expressing appropriate levels of TLR9 respond to CpG-ODN, but the nature of the TLR9-mediated response is dependent on cell type.