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.     

Structures and recognition modes of toll‐like receptors

Toll‐like receptors (TLRs) recognize common structural patterns in diverse microbial molecules and play central roles in the innate immune response. The structures of extracellular domains and their ligand complexes of several TLRs have been determined by X‐ray crystallography. Here, we discuss recent advances on structures and activation mechanisms of TLRs. Despite the differences in interaction areas of ligand with TLRs, the extracellular domains of TLRs all adopt horseshoe‐shaped structures and the overall M‐shape of the TLR–ligand complexes is strikingly similar. The structural rearrangement information of TLRs sheds new light on their ligand‐recognition and ‐activation mechanisms. Proteins 2016; 85:3–9. © 2016 Wiley Periodicals, Inc.     

The evolution of bat nucleic acid‐sensing Toll‐like receptors

We characterized the nucleic acid‐sensing Toll‐like receptors (TLR) of a New World bat species, the common vampire bat (Desmodus rotundus), and through a comparative molecular evolutionary approach searched for general adaptation patterns among the nucleic acid‐sensing TLRs of eight different bats species belonging to three families (Pteropodidae, Vespertilionidae and Phyllostomidae). We found that the bat TLRs are evolving slowly and mostly under purifying selection and that the divergence pattern of such receptors is overall congruent with the species tree, consistent with the evolution of many other mammalian nuclear genes. However, the chiropteran TLRs exhibited unique mutations fixed in ligand‐binding sites, some of which involved nonconservative amino acid changes and/or targets of positive selection. Such changes could potentially modify protein function and ligand‐binding properties, as some changes were predicted to alter nucleic acid binding motifs in TLR 9. Moreover, evidence for episodic diversifying selection acting specifically upon the bat lineage and sublineages was detected. Thus, the long‐term adaptation of chiropterans to a wide variety of environments and ecological niches with different pathogen profiles is likely to have shaped the evolution of the bat TLRs in an order‐specific manner. The observed evolutionary patterns provide evidence for potential functional differences between bat and other mammalian TLRs in terms of resistance to specific pathogens or recognition of nucleic acids in general.     

Crystal structure of the C‐terminal domain of mouse TLR9

Toll‐like receptors (TLRs) are important pattern recognition receptors that function in innate immunity. Elucidating the structure and signaling mechanisms of TLR9, a sensor of foreign and endogenous DNA, is essential for understanding its key role in immunity against microbial pathogens as well as in autoimmunity. Abundant evidence suggests that the TLR9‐CTD (C‐terminal domain) by itself is capable of DNA binding and signaling. The crystal structure of unliganded mouse TLR9‐CTD is presented. TLR9‐CTD exhibits one unique feature, a cluster of stacked aromatic and arginine side chains on its concave face. Overall, its structure is most related to the TLR8‐CTD, suggesting a similar mode of ligand binding and signaling. Proteins 2014; 82:2874–2878. © 2014 Wiley Periodicals, Inc.     

A leucine-rich repeat assembly approach for homology modeling of the human TLR5-10 and mouse TLR11-13 ectodomains

So far, 13 groups of mammalian Toll-like receptors (TLRs) have been identified. Most TLRs have been shown to recognize pathogen-associated molecular patterns from a wide range of invading agents and initiate both innate and adaptive immune responses. The TLR ectodomains are composed of varying numbers and types of leucine-rich repeats (LRRs). As the crystal structures are currently missing for most TLR ligand-binding ectodomains, homology modeling enables first predictions of their three-dimensional structures on the basis of the determined crystal structures of TLR ectodomains. However, the quality of the predicted models that are generated from full-length templates can be limited due to low sequence identity between the target and templates. To obtain better templates for modeling, we have developed an LRR template assembly approach. Individual LRR templates that are locally optimal for the target sequence are assembled into multiple templates. This method was validated through the comparison of a predicted model with the crystal structure of mouse TLR3. With this method, we also constructed ectodomain models of human TLR5, TLR6, TLR7, TLR8, TLR9, and TLR10 and mouse TLR11, TLR12, and TLR13 that can be used as first passes for a computational simulation of ligand docking or to design mutation experiments. This template assembly approach can be extended to other repetitive proteins.     

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.     

Interleukin‐8 and CXCL10 expression in oral keratinocytes and fibroblasts via Toll‐like receptors

Oral keratinocytes and fibroblasts may be the first line of host defense against oral microorganisms. Here, the contention that oral keratinocytes and fibroblasts recognize microbial components via Toll‐like receptors (TLRs) and participate in development of oral inflammation was examined. It was found that immortalized oral keratinocytes (RT7), fibroblasts (GT1) and primary cells express mRNA of TLRs 1–10. Interleukin‐8 (IL‐8) production by RT7 cells was induced by treatment with TLRs 1–9 with the exception of TLR7 agonist, whereas GT1 cells were induced to produce IL‐8 by all TLR agonists tested except for TLR7 and TLR9. GT1 cells showed increased CXCL10 production following treatment with agonists for TLR1/2, TLR3, TLR4, and TLR5, whereas only those for TLR3 and TLR5 increased CXCL10 production in RT7 cells. Moreover, TLR agonists differentially regulated tumor necrosis factor‐alpha‐induced IL‐8 and CXCL10 production by the tested cell types. These findings suggest that recognition of pathogenic microorganisms in oral keratinocytes and fibroblasts by TLRs may have important roles in orchestrating host immune responses via production of various chemokines.     

Signatures of diversifying selection and convergence acting on passerine Toll‐like receptor 4 in an evolutionary context

Positive selection acting on Toll‐like receptors (TLRs) has been recently investigated to reveal evolutionary mechanisms of host–pathogen molecular co‐adaptation. Much of this research, however, has focused mainly on the identification of sites predicted to be under positive selection, bringing little insight into the functional differences and similarities among species and a limited understanding of convergent evolution in the innate immune molecules. In this study, we provide evidence of phenotypic variability in the avian TLR4 ligand‐binding region (LBR), the direct interface between host and pathogen molecular structures. We show that 55 passerine species vary substantially in the distribution of electrostatic potential on the surface of the receptor, and based on these distinct patterns, we identified four species clusters. Seven of the 34 evolutionarily nonconservative and positively selected residues correspond topologically to sites previously identified as being important for lipopolysaccharide, lipid IVa or MD‐2 binding. Five of these positions codetermine the identity of the charge clusters. Groups of species that host‐related communities of pathogens were predicted to cluster based on their TLR4 LBR charge. Despite some evidence for convergence among taxa, there were no clear associations between the TLR4 LBR charge distribution and any of the general ecological characteristics compared (migration, latitudinal distribution and diet). Closely related species, however, mostly belonged to the same surface charge cluster indicating that phylogenetic constraints are key determinants shaping TLR4 adaptive evolution. Our results suggest that host innate immune evolution is consistent with Fahrenholz's rule on the cospeciation of hosts and their parasites.     

Activation of Toll‐like receptor 9 inhibits lipopolysaccharide‐induced receptor activator of nuclear factor kappa‐ B ligand expression in rat B lymphocytes

B lymphocytes express multiple TLRs that regulate their cytokine production. We investigated the effect of TLR4 and TLR9 activation on receptor activator of NF‐κB ligand (RANKL) expression by rat spleen B cells. Splenocytes or purified spleen B cells from Rowett rats were cultured with TLR4 ligand Escherichia coli LPS and/or TLR9 ligand CpG‐oligodeoxynucleotide (CpG‐ODN) for 2 days. RANKL mRNA expression and the percentage of RANKL‐positive B cells were increased in rat splenocytes challenged by E. coli LPS alone. The increases were less pronounced when cells were treated with both CpG‐ODN and E. coli LPS. Microarray analysis showed that expressions of multiple cyclin‐dependent kinase (CDK) pathway‐related genes were up‐regulated only in cells treated with both E. coli LPS and CpG‐ODN. This study suggests that CpG‐ODN inhibits LPS‐induced RANKL expression in rat B cells via regulation of the CDK pathway.     

Comparative analysis of species-specific ligand recognition in Toll-like receptor 8 signaling: a hypothesis

Toll-like receptors (TLRs) play a central role in the innate immune response by recognizing conserved structural patterns in a variety of microbes. TLRs are classified into six families, of which TLR7 family members include TLR7, 8, and 9, which are localized to endolysosomal compartments recognizing viral infection in the form of foreign nucleic acids. In our current study, we focused on TLR8, which has been shown to recognize different types of ligands such as viral or bacterial ssRNA as well as small synthetic molecules. The primary sequences of rodent and non-rodent TLR8s are similar, but the antiviral compound (R848) that activates the TLR8 pathway is species-specific. Moreover, the factors underlying the receptor's species-specificity remain unknown. To this end, comparative homology modeling, molecular dynamics simulations refinement, automated docking and computational mutagenesis studies were employed to probe the intermolecular interactions between this anti-viral compound and TLR8. Furthermore, comparative analyses of modeled TLR8 (rodent and non-rodent) structures have shown that the variation mainly occurs at LRR14-15 (undefined region); hence, we hypothesized that this variation may be the primary reason for the exhibited species-specificity. Our hypothesis was further bolstered by our docking studies, which clearly showed that this undefined region was in close proximity to the ligand-binding site and thus may play a key role in ligand recognition. In addition, the interface between the ligand and TLR8s varied depending upon the amino acid charges, free energy of binding, and interaction surface. Therefore, our current work provides a hypothesis for previous in vivo studies in the context of TLR signaling.     

Up‐regulated TLR4 in cardiomyocytes exacerbates heart failure after long‐term myocardial infarction

It remains unclear whether and how cardiomyocytes contribute to the inflammation in chronic heart failure (CHF). We recently reviewed the capacity of cardiomyocytes to initiate inflammation, by means of expressing certain immune receptors such as toll‐like receptors (TLRs) that respond to pathogen‐ and damage‐associated molecular patterns (PAMP and DAMP). Previous studies observed TLR4‐mediated inflammation within days of myocardial infarction (MI). This study examined TLR4 expression and function in cardiomyocytes of failing hearts after 4 weeks of MI in rats. The increases of TLR4 mRNA and proteins, as well as inflammatory cytokine production, were observed in both the infarct and remote myocardium. Enhanced immunostaining for TLR4 was observed in cardiomyocytes but not infiltrating leucocytes. The injection of lentivirus shRNA against TLR4 into the infarcted heart decreased inflammatory cytokine production and improved heart function in vivo. Accordingly, in cardiomyocytes isolated from CHF hearts, increases of TLR4 mRNA and proteins were detected. More robust binding of TLR4 with lipopolysaccharide (LPS), a PAMP ligand for TLR4, and heat shock protein 60 (HSP60), a DAMP ligand for TLR4, was observed in CHF cardiomyocytes under a confocal microscope. The maximum binding capacity (B_max) of TLR4 was increased for LPS and HSP60, whereas the binding affinity (Kd) was not significantly changed. Furthermore, both LPS and HSP60 induced more robust production of inflammatory cytokines in CHF cardiomyocytes, which was reduced by TLR4‐blocking antibodies. We conclude that the expression, ligand‐binding capacity and pro‐inflammatory function of cardiomyocyte TLR4 are up‐regulated after long‐term MI, which promote inflammation and exacerbate heart failure.     

The functional effects of physical interactions among Toll-like receptors 7, 8, and 9

Toll-like receptors (TLRs) TLR1, TLR2, TLR4, and TLR6 are evolutionarily conserved, highly homologous, and localized to plasma membranes of host cells and recognize pathogen-associated molecular patterns (PAMPs) derived from bacterial membranes. These receptors cooperate in a pairwise combination to elicit or inhibit the inflammatory signals in response to certain PAMPs. The other TLRs that are evolutionarily closely related and highly homologous are TLR7, TLR8, and TLR9. They are all confined to the membranes of endosomes and recognize similar molecular structures, the oligonucleotide-based PAMPs. However, the cooperative interactions among these receptors that may modulate the inflammatory signaling in response to their cognate agonists are not reported. We report here for the first time the functional effects of one TLR on the other among TLR7, TLR8, and TLR9. The results indicate that TLR8 inhibits TLR7 and TLR9, and TLR9 inhibits TLR7 but not vice versa in HEK293 cells transfected with TLRs in a pairwise combination. This is concluded by selectively activating one TLR over the other by using small molecule TLR agonists. We also show that these inhibitory interactions are the result of direct or indirect physical interactions between the TLRs. The murine TLR8 that does not respond to any known human TLR8 agonists also inhibits both murine and human TLR7. The implications of the inhibitory interactions among these TLRs in host-pathogen recognition and subsequent inflammatory responses are not obvious. However, given the complexity in expression pattern in a particular cell type and the variation in distribution and response to different pathogens and stress signals in different cell types, the inhibitory physical interactions among these TLRs may play a role in balancing the inflammatory outcome from a given cell type to a specific challenge.     

TLR stimulation modifies BLyS receptor expression in follicular and marginal zone B cells

Through their differential interactions with B lymphocyte stimulator (BLyS) and a proliferation-inducing ligand (APRIL), the three BLyS family receptors play central roles in B cell survival and differentiation. Recent evidence indicates BLyS receptor levels shift following BCR ligation, suggesting that activation cues can alter overall BLyS receptor profiles and thus ligand sensitivity. In this study, we show that TLR stimuli also alter BLyS receptor expression, but in contrast to BCR ligation, TLR9 and TLR4 signals, preferentially increase transmembrane activator calcium modulator and cyclophilin ligand interactor (TACI) expression. Although both of these TLRs act through MyD88-dependent mechanisms to increase TACI expression, they differ in terms of their downstream mediators and the B cell subset affected. Surprisingly, only TLR4 relies on c-Rel and p50 to augment TACI expression, whereas TLR9 does not. Furthermore, although all follicular and marginal zone B cells up-regulate TACI in response to TLR9 stimulation, only marginal zone B cells and a subset of follicular B cells respond to TLR4. Finally, we find that both BLyS and APRIL enhance viability among quiescent and BCR-stimulated B cells. However, although BLyS enhances viability among TLR stimulated B cells, APRIL does not, suggesting that TACI but not BLyS receptor 3 may share survival promoting pathways with TLRs.     

Positive selection and convergent evolution shape molecular phenotypic traits of innate immunity receptors in tits (Paridae)

Despite widespread variability and redundancy abounding animal immunity, little is currently known about the rate of evolutionary convergence (functionally analogous traits not inherited from a common ancestor) in host molecular adaptations to parasite selective pressures. Toll‐like receptors (TLRs) provide the molecular interface allowing hosts to recognize pathogenic structures and trigger early danger signals initiating an immune response. Using a novel combination of bioinformatic approaches, here we explore genetic variation in ligand‐binding regions of bacteria‐sensing TLR4 and TLR5 in 29 species belonging to the tit family of passerine birds (Aves: Paridae). Three out of the four consensual positively selected sites in TLR4 and six out of 14 positively selected positions in TLR5 were located on the receptor surface near the functionally important sites, and based on the phylogenetic pattern evolved in a convergent (parallel) manner. This type of evolution was also seen at one N‐glycosylation site and two positively selected phosphorylation sites, providing the first evidence of convergence in post‐translational modifications in evolutionary immunology. Finally, the overall mismatch between phylogeny and the clustering of surface charge distribution demonstrates that convergence is common in overall TLR4 and TLR5 molecular phenotypes involved in ligand binding. Our analysis did not reveal any broad ecological traits explaining the convergence observed in electrostatic potentials, suggesting that information on microbial symbionts may be needed to explain TLR evolution. Adopting state‐of‐the‐art predictive structural bionformatics, we have outlined a new broadly applicable methodological approach to estimate the functional significance of positively selected variation and test for the adaptive molecular convergence in protein‐coding polymorphisms.     

UNC93B1 physically associates with human TLR8 and regulates TLR8-mediated signaling

Toll-like receptors (TLRs) 3, 7, 8, and 9 are localized to intracellular compartments where they encounter foreign or self nucleic acids and activate innate and adaptive immune responses. The endoplasmic reticulum (ER)-resident membrane protein, UNC93B1, is essential for intracellular trafficking and endolysosomal targeting of TLR7 and TLR9. TLR8 is phylogenetically and structurally related to TLR7 and TLR9, but little is known about its localization or function. In this study, we demonstrate that TLR8 localized to the early endosome and the ER but not to the late endosome or lysosome in human monocytes and HeLa transfectants. UNC93B1 physically associated with human TLR8, similar to TLRs 3, 7, and 9, and played a critical role in TLR8-mediated signaling. Localization analyses of TLR8 tail-truncated mutants revealed that the transmembrane domain and the Toll/interleukin-1 receptor domain were required for proper targeting of TLR8 to the early endosome. Hence, although UNC93B1 participates in intracellular trafficking and signaling for all nucleotide-sensing TLRs, the mode of regulation of TLR localization differs for each TLR.     

TLR3 engagement induces IRF‐3‐dependent apoptosis in androgen‐sensitive prostate cancer cells and inhibits tumour growth in vivo

Toll‐like receptors (TLRs) are a family of highly conserved transmembrane proteins expressed in epithelial and immune cells that recognize pathogen associated molecular patterns. Besides their role in immune response against infections, numerous studies have shown an important role of different TLRs in cancer, indicating these receptors as potential targets for cancer therapy. We previously demonstrated that the activation of TLR3 by the synthetic double‐stranded RNA analogue poly I:C induces apoptosis of androgen‐sensitive prostate cancer (PCa) LNCaP cells and, much less efficiently, of the more aggressive PC3 cell line. Therefore, in this study we selected LNCaP cells to investigate the mechanism of TLR3‐mediated apoptosis and the in vivo efficacy of poly I:C‐based therapy. We show that interferon regulatory factor‐3 (IRF‐3) signalling plays an essential role in TLR3‐mediated apoptosis in LNCaP cells through the activation of the intrinsic and extrinsic apoptotic pathways. Interestingly, hardly any apoptosis was induced by poly I:C in normal prostate epithelial cells RWPE‐1. We also demonstrate for the first time the direct anticancer effect of poly I:C as a single therapeutic agent in a well‐established human androgen‐sensitive PCa xenograft model, by showing that tumour growth is highly impaired in poly I:C‐treated immunodeficient mice. Immunohistochemical analysis of PCa xenografts highlights the antitumour role of poly I:C in vivo both on cancer cells and, indirectly, on endothelial cells. Notably, we show the presence of TLR3 and IRF‐3 in both human normal and PCa clinical samples, potentially envisaging poly I:C‐based therapy for PCa.     

Homology modeling and structural comparison of leucine rich repeats of toll like receptors 1-10 of ruminants

Toll-like receptors (TLRs) are transmembrane receptors composed of extra cellular leucine rich repeats (LRRs) that identify specific pathogen associated molecular patterns triggering a innate immune cascade. The LRR regions of TLR 1–10 proteins of goat (Capra hircus), sheep (Ovis aries), buffalo (Bubalus bubalis) and bovine (Bos taurus) were modeled using MODELLER 9v7 tool and validated. The similarities and variations of these 10 TLRs extracellular regions of each species were compared using online servers like FATCAT, SSM and SSAP. It was evident that the LRRs of TLRs like 1, 2, 3 and 6 showed structural convergence with <1 % RMSD deviation while TLRs like 5, 7, 8 and 9 had high divergence. Docking analysis showed that TLR 2, 3 and 7 of all the selected four ruminant species were able to bind with their corresponding ligands like Peptidoglycan (PGN), Poly I:C, Resiquimod (R-848) and Imiquimod. However, there were variations in the active site regions, interacting residues and the number of bonded interactions. Variations seen among TLR structures and their ligand binding characteristics is likely to be responsible for species and breed specific genetic resistance observed among species or breeds.     

Toll‐like receptor signalling cross‐activates the autophagic pathway to restrict Salmonella Typhimurium growth in macrophages

It has been long recognised that activation of toll‐like receptors (TLRs) induces autophagy to restrict intracellular bacterial growth. However, the mechanisms of TLR‐induced autophagy are incompletely understood. Salmonella Typhimurium is an intracellular pathogen that causes food poisoning and gastroenteritis in humans. Whether TLR activation contributes to S. Typhimurium‐induced autophagy has not been investigated. Here, we report that S. Typhimurium and TLRs shared a common pathway to induce autophagy in macrophages. We first showed that S. Typhimurium‐induced autophagy in a RAW264.7 murine macrophage cell line was mediated by the AMP‐activated protein kinase (AMPK) through activation of the TGF‐β‐activated kinase (TAK1), a kinase activated by multiple TLRs. AMPK activation led to increased phosphorylation of Unc‐51‐like autophagy activating kinase (ULK1) at S317 and S555. ULK1 phosphorylation at these two sites in S. Typhimurium‐infected macrophages overrode the inhibitory effect of mTOR on ULK1 activity due to mTOR‐mediated ULK1 phosphorylation at S757. Lipopolysaccharide (LPS), flagellin, and CpG oligodeoxynucleotide, which activate TLR4, TLR5, and TLR9, respectively, increased TAK1 and AMPK phosphorylation and induced autophagy in RAW264.7 cells and in bone marrow‐derived macrophages. However, LPS was unable to induce TAK1 and AMPK phosphorylation and autophagy in TLR4‐deficient macrophages. TAK1 and AMPK‐specific inhibitors blocked S. Typhimurium‐induced autophagy and xenophagy and increased the bacterial growth in RAW264.7 cells. These observations collectively suggest that activation of the TAK1–AMPK axis through TLRs is essential for S. Typhimurium‐induced autophagy and that TLR signalling cross‐activates the autophagic pathway to clear intracellular bacteria.