Identification of MDP (muramyl dipeptide)-binding key domains in NOD2 (nucleotide-binding and oligomerization domain-2) receptor of Labeo rohita

Functional characterization of RIP2 in large yellow croaker (Larimichthys crocea), a protein involved in the host antiviral responses via NF-κΒ, IRF3/7 related signaling

As an adaptor protein functions essentially in the activation of NF-κΒ and MAPK signaling pathways mediated by NOD1 and NOD2, RIP2 plays important roles in the host innate immune responses. In the present study, the RIP2 ortholog termed Lc-RIP2 was identified and characterized in large yellow croaker (Larimichthys crocea). It was revealed that Lc-RIP2 is consisted of an open reading frame (ORF) of 1695 bp, encoding a protein of 564 aa, with an N-terminal kinase domain and a C-terminal caspase activation and recruitment domain (CARD). Subcellular localization assays demonstrated that Lc-RIP2 was a cytosolic protein, which was broadly distributed in the examined tissues/organs, and could be induced in response to poly I:C, LPS, PGN, and Pseudomonas plecoglossicida stimulations in vivo according to qRT-PCR analysis. Notably, Lc-RIP2 overexpression in vitro was sufficient to abolish SVCV proliferation in EPC cells, and could significantly induce the activation of NF-κB, IRF3, IRF7, and IFN1 promoters. In addition, luciferase assays found that Lc-RIP2 could cooperate with Lc-MAVS, Lc-TRAF3, Lc-TRAF6, Lc-IRF3, and Lc-IRF7 in NF-κB activation, associate with Lc-TRIF, Lc-MAVS, Lc-TRAF3, Lc-IRF3, and Lc-IRF7 in IRF3 activation, enhance Lc-TRIF, Lc-MAVS, Lc-TRAF3, and Lc-TRAF6 mediated IRF7 activation, and Lc-IRF3 mediated IFN1 activation, whereas suppress NF-κB activation when co-expressed with Lc-TRIF. Co-immunoprecipitation (Co-IP) assays also demonstrated that Lc-RIP2 interacts separately with Lc-TRIF, Lc-MAVS, Lc-TRAF3, Lc-TRAF6, Lc-IRF3, and Lc-IRF7. It is thus collectively indicated that Lc-RIP2 function dominantly in the regulation of the host innate immune signaling.

Structure of fish Toll-like receptors (TLR) and NOD-like receptors (NLR)

Innate immunity driven by pattern recognition receptor (PRR) protects the host from invading pathogens. Aquatic animals like fish where the adaptive immunity is poorly developed majorly rely on their innate immunity modulated by PRRs like toll-like receptors (TLR) and NOD-like receptors (NLR). However, current development to improve the fish immunity via TLR/NLR signaling is affected by a poor understanding of its mechanistic and structural features. This review discusses the structure of fish TLRs/NLRs and its interaction with pathogen associated molecular patterns (PAMPs) and downstream signaling molecules. Over the past one decade, significant progress has been done in studying the structure of TLRs/NLRs in higher eukaryotes; however, structural studies on fish innate immune receptors are undermined. Several novel TLR genes are identified in fish that are absent in higher eukaryotes, but the function is still poorly understood. Unlike the fundamental progress achieved in developing antagonist/agonist to modulate human innate immunity, analogous studies in fish are nearly lacking due to structural inadequacy. This underlies the importance of exploring the structural and mechanistic details of fish TLRs/NLRs at an atomic and molecular level. This review outlined the mechanistic and structural basis of fish TLR and NLR activation.

Involvement of ayu NOD2 in NF-κB and MAPK signaling pathways: Insights into functional conservation of NOD2 in antibacterial innate immunity

Nucleotide oligomerization domain 2 (NOD2) is a major cytoplasmic sensor for pathogens and is critical for the clearance of cytosolic bacteria in mammals. However, studies regarding NOD2, especially the initiated signaling pathways, are scarce in teleost species. In this study, we identified a NOD2 molecule (PaNOD2) from ayu (Plecoglossus altivelis). Bioinformatics analysis showed the structure of NOD2 to be highly conserved during vertebrate evolution. Dual-luciferase reporter assays examined the activation of NF-κB signaling and Western blotting analysis detected the phosphorylation of three MAP kinases (p-38, Erk1/2, and JNK1/2). Functional study revealed that, like its mammalian counterparts, PaNOD2 was the receptor of the bacterial cell wall component muramyl dipeptide (MDP), and the leucine-rich repeat motif was responsible for the recognition and binding of PaNOD2 with the ligand. Overexpression of PaNOD2 activated the NF-κB signaling pathway, leading to the upregulation of inflammatory cytokines, including TNF-α and IL-1β in HEK293T cells and ayu head kidney-derived monocytes/macrophages (MO/MΦ). Particularly, we found that PaNOD2 activated the MAPK signaling pathways, as indicated by the increased phosphorylation of p-38, Erk1/2, and JNK1/2, which have not been characterized in any teleost species previously. Our findings proved that the NOD2 molecule and initiated pathways are conserved between mammals and ayu. Therefore, ayu could be used as an animal model to investigate NOD2-based diseases and therapeutic applications.

Molecular characterization, expression and functional analysis of NOD1, NOD2 and NLRC3 in Nile tilapia (Oreochromis niloticus)

The nucleotide-binding oligomerization domain proteins NOD1, NOD2 and NLRC3 are cytoplasmic pattern recognition receptors (PRRs) of the Nod-like receptor (NLR) family. In the present study, the Nile tilapia (Oreochromis niloticus) NOD1 (ntNOD1), NOD2 (ntNOD2) and NLRC3 (ntNLRC3) genes were cloned and characterized. The full-length ntNOD1, ntNOD2 and ntNLRC3 genes were 3924, 3886 and 4574 bp, encoding 941, 986 and 1130 amino acids, respectively. The three Nod-like receptors have a NACHT domain and a C-terminal leucine-rich repeat (LRR) domain. In addition, ntNOD1 and ntNOD2 have a N-terminal CARD domain (ntNOD2 has two). Phylogenetic analysis showed that the three NLRs are highly conserved. Tissue expression analysis of the three receptors revealed that the highest mRNA and protein levels of ntNOD1, ntNOD2 and ntNLRC3 were in the spleen. The expression patterns of NLRs during embryonic development showed that the expression levels of ntNOD2 and ntNLRC3 significantly increased from 2 to 8 days post-fertilization (dpf). The expression levels of ntNOD1 significantly increased from 2 to 6 dpf, decreased at 7 dpf and then increased at 8 dpf. Upon stimulation with an intraperitoneal injection of Streptococcus agalactiae, expression levels of the ntNOD1, ntNOD2 and ntNLRC3 mRNA and protein were clearly altered in the blood, spleen, kidney, intestine and gill. Furthermore, after cotransfection with an NF-κB reporter plasmid, NF-κB activation in ntNOD1-overexpressing 293T cells significantly increased compared with that in control cells, before or after i-EDPA-stimulation. By contrast, compared with control, ntNOD2 and ntNLRC3 had no effect on NF-κB activation in 293T cells, when their potential ligands were not stimulated. However, after MDP-stimulation, ntNOD2 and ntNLRC3 overexpression increased NF-κB activation in 293T cells. NOD1 and NLRC3 were uniformly distributed throughout the cytoplasm in 293T cells, whereas NOD2 was distributed throughout the cytoplasm and nucleus. Our results indicate that the three Nod-like receptors are functionally conserved and may play pivotal roles in defense against pathogens such as Streptococcus agalactiae.

Understanding the molecular differential recognition of muramyl peptide ligands by LRR domains of human NOD receptors

Human nucleotide-binding oligomerization domain proteins, hNOD1 and hNOD2, are host intracellular receptors with C-terminal leucine-rich repeat (LRR) domains, which recognize specific bacterial peptidoglycan (PG) fragments as their ligands. The specificity of this recognition is dependent on the third amino acid of the stem peptide of the PG ligand, which is usually meso-diaminopimelic acid (mesoDAP) or l-lysine (l-Lys). Since the LRR domains of hNOD receptors had been experimentally shown to confer the PG ligand-sensing specificity, we developed three-dimensional structures of hNOD1-LRR and the hNOD2-LRR to understand the mechanism of differential recognition of muramyl peptide ligands by hNOD receptors. The hNOD1-LRR and hNOD2-LRR receptor models exhibited right-handed curved solenoid shape. The hot-spot residues experimentally proved to be critical for ligand recognition were located in the concavity of the NOD-LRR and formed the recognition site. Our molecular docking analyses and molecular electrostatic potential mapping studies explain the activation of hNOD-LRRs, in response to effective molecular interactions of PG ligands at the recognition site; and conversely, the inability of certain PG ligands to activate hNOD-LRRs, by deviations from the recognition site. Based on molecular docking studies using PG ligands, we propose few residues — G825, D826 and N850 in hNOD1-LRR and L904, G905, W931, L932 and S933 in hNOD2-LRR, evolutionarily conserved across different host species, which may play a major role in ligand recognition. Thus, our integrated experimental and computational approach elucidates the molecular basis underlying the differential recognition of PG ligands by hNOD receptors.

Mutual Regulation of NOD2 and RIG-I in Zebrafish Provides Insights into the Coordination between Innate Antibacterial and Antiviral Signaling Pathways

Nucleotide-binding oligomerization domain-containing protein 2 (NOD2) and retinoic acid-inducible gene I (RIG-I) are two important cytosolic pattern recognition receptors (PRRs) in the recognition of pathogen-associated molecular patterns (PAMPs), initiating innate antibacterial and antiviral signaling pathways. However, the relationship between these PRRs, especially in teleost fish models, is rarely reported. In this article, we describe the mutual regulation of zebrafish NOD2 (DrNOD2) and RIG-I (DrRIG-I) in innate immune responses. Luciferase assays were conducted to determine the activation of NF-κB and interferon signaling. Morpholino-mediated knockdown and mRNA-mediated rescue were performed to further confirm the regulatory roles between DrNOD2 and DrRIG-I. Results showed that DrNOD2 and DrRIG-I shared conserved structural hallmarks with their mammalian counterparts, and activated DrRIG-I signaling can induce DrNOD2 production. Surprisingly, DrNOD2-initiated signaling can also induce DrRIG-I expression, indicating that a mutual regulatory mechanism may exist between them. Studies conducted using HEK293T cells and zebrafish embryos showed that DrRIG-I could negatively regulate DrNOD2-activated NF-κB signaling, and DrNOD2 could inhibit DrRIG-I-induced IFN signaling. Moreover, knocking down DrRIG-I expression by morpholino could enhance DrNOD2-initiated NF-κB activation, and vice versa, which could be rescued by their corresponding mRNAs. Results revealed a mutual feedback regulatory mechanism underlying NOD2 and RIG-I signaling pathways in teleosts. This mechanism reflects the coordination between cytosolic antibacterial and antiviral PRRs in the complex network of innate immunity.

Casting the critical regions in nucleotide binding oligomerization domain 2 protein: a signature mediated structural dynamics approach

Nucleotide binding oligomerization domain 2 (NOD2), a protein involved in the first line defence mechanism has a pivotal role in innate immunity. Impaired function of this protein is implicated in disorders such as Blau syndrome and Crohn's disease. Since an altered function is linked to protein's structure, we framed a systematic strategy to interpret the structure-function relationship of the protein. Initiated with mutation-based pattern prediction and identified a distant ortholog (DO) of NOD2 from which the intra-residue interaction network was elucidated. The network was used to identify hotspots that serve as critical points to maintain the stable architecture of the protein. Structural comparison of NOD2 domains with a DO revealed the minimal number of intra-protein interactions required by the protein to maintain the structural fold. In addition, the conventional molecular dynamics simulation emphasized the conformational transitions at hot spot residues between native NOD2 domains and its respective mutants (G116R, R42W and R54A) structures. The analysis of intra-protein interactions globally and the displacement of residues locally around the mutational site revealed loss of several critical bonds and residues vital for the protein's function. Conclusively we report, about 10 residues in leucine-rich repeat, 13 residues in NOD and 6 residues in CARD domain are required by the NOD2 to maintain its function. This protocol will help the researchers to achieve for more prospective studies to attest druggable site utility in discovering novel drug candidates.

A genome-wide survey of expansive NLR-C subfamily in miiuy croaker and characterization of the NLR-B30.2 genes

NOD-like receptors (NLRs) are essential intracellular pattern-recognition receptors that respond to pathogens and regulate innate immunity. NLRs include three distinct subfamilies: NLR-A, NLR-B and NLR-C, thereinto, NLR-C as a large subfamily is unique to bony fish and little research about it has been done. In the current study, we identified the members of NLR-B and NLR-C subfamilies containing 2 and 48 genes respectively in miiuy croaker. Compared with other teleosts except for zebrafish, NLR-C subfamily genes occurred expansion in miiuy croaker. The gene expansions of NLR-C subfamily may illustrate adaptive genome evolution in response to specific aquatic environments. Structural analysis showed that the N-terminus of NLR-C subfamily receptors has different characteristics of the domains including RING domain, FISNA domain or PYRIN domain. Interestingly, the C-terminus of 18 NLR-C subfamily members contains an extra B30.2 domain (named NLR-B30.2 genes) which plays an important role in antiviral immune recognition. Simultaneously, molecular evolutionary analysis indicated that the positively sites in miiuy croaker are mainly located in NACHT domain which was the vital region for signal transduction in immune response. Significantly, pathogens challenge in spleen and macrophages demonstrated that NLR-B30.2 genes exhibited more sensitive response to virus than bacteria, suggesting these genes play enhanced roles in innate antiviral immunity, which may represent a new family used for antiviral infection.

Structural and functional insights into CARDs of zebrafish (Danio rerio) NOD1 and NOD2, and their interaction with adaptor protein RIP2

Nucleotide-binding and oligomerization domain-containing protein 1 (NOD1) and NOD2 are cytosolic pattern-recognition receptors (PRRs) composed of an N-terminal caspase activation and recruitment domain (CARD), a central NACHT domain and C-terminal leucine-rich repeats (LRRs). They play a vital role in innate immune signaling by activating the NF-κB pathway via recognition of peptidoglycans by LRRs, and ATP-dependent self-oligomerization of NACHT followed by downstream signaling. After oligomerization, CARD/s play a crucial role in activating downstream signaling via the adaptor molecule, RIP2. Due to the inadequacy of experimental 3D structures of CARD/s of NOD2 and RIP2, and results from differential experimental setups, the RIP2-mediated CARD-CARD interaction has remained as a contradictory statement. We employed a combinatorial approach involving protein modeling, docking, molecular dynamics simulation, and binding free energy calculation to illuminate the molecular mechanism that shows the possible involvement of either the acidic or basic patch of zebrafish NOD1/2-CARD/a and RIP2-CARD in CARD-CARD interaction. Herein, we have hypothesized 'type-I' mode of CARD-CARD interaction in NOD1 and NOD2, where NOD1/2-CARD/a involve their acidic surfaces to interact with RIP2. Asp37 and Glu51 (of NOD1) and Arg477, Arg521 and Arg529 (of RIP2) were identified to be crucial for NOD1-RIP2 interaction. However, in NOD2-RIP2, Asp32 (of NOD2) and Arg477 and Arg521 (of RIP2) were anticipated to be significant for downstream signaling. Furthermore, we found that strong electrostatic contacts and salt bridges are crucial for protein-protein interactions. Altogether, our study has provided novel insights into the RIP2-mediated CARD-CARD interaction in zebrafish NOD1 and NOD2, which will be helpful to understand the molecular basis of the NOD1/2 signaling mechanism.

Structural models of zebrafish (Danio rerio) NOD1 and NOD2 NACHT domains suggest differential ATP binding orientations: insights from computational modeling, docking and molecular dynamics simulations

Nucleotide-binding oligomerization domain-containing protein 1 (NOD1) and NOD2 are cytosolic pattern recognition receptors playing pivotal roles in innate immune signaling. NOD1 and NOD2 recognize bacterial peptidoglycan derivatives iE-DAP and MDP, respectively and undergoes conformational alternation and ATP-dependent self-oligomerization of NACHT domain followed by downstream signaling. Lack of structural adequacy of NACHT domain confines our understanding about the NOD-mediated signaling mechanism. Here, we predicted the structure of NACHT domain of both NOD1 and NOD2 from model organism zebrafish (Danio rerio) using computational methods. Our study highlighted the differential ATP binding modes in NOD1 and NOD2. In NOD1, γ-phosphate of ATP faced toward the central nucleotide binding cavity like NLRC4, whereas in NOD2 the cavity was occupied by adenine moiety. The conserved ‘Lysine’ at Walker A formed hydrogen bonds (H-bonds) and Aspartic acid (Walker B) formed electrostatic interaction with ATP. At Sensor 1, Arg328 of NOD1 exhibited an H-bond with ATP, whereas corresponding Arg404 of NOD2 did not. ‘Proline’ of GxP motif (Pro386 of NOD1 and Pro464 of NOD2) interacted with adenine moiety and His511 at Sensor 2 of NOD1 interacted with γ-phosphate group of ATP. In contrast, His579 of NOD2 interacted with the adenine moiety having a relatively inverted orientation. Our findings are well supplemented with the molecular interaction of ATP with NLRC4, and consistent with mutagenesis data reported for human, which indicates evolutionary shared NOD signaling mechanism. Together, this study provides novel insights into ATP binding mechanism, and highlights the differential ATP binding modes in zebrafish NOD1 and NOD2.

A conformational analysis of mouse Nalp3 domain structures by molecular dynamics simulations, and binding site analysis

Scrutinizing various nucleotide-binding oligomerization domain (NOD)-like receptor (NLR) genes in higher eukaryotes is very important for understanding the intriguing mechanism of the host defense against pathogens. The nucleotide-binding domain (NACHT), leucine-rich repeat (LRR), and pyrin domains (PYD)-containing protein 3 (Nalp3), is an intracellular innate immune receptor and is associated with several immune system related disorders. Despite Nalp3's protective role during a pathogenic invasion, the molecular features and structural organization of this crucial protein is poorly understood. Using comparative modeling and molecular dynamics simulations, we have studied the structural architecture of Nalp3 domains, and characterized the dynamic and energetic parameters of adenosine triphosphate (ATP) binding in NACHT, and pathogen-derived ligands muramyl dipeptide (MDP) and imidazoquinoline with LRR domains. The results suggested that walker A, B and extended walker B motifs were the key ATP binding regions in NACHT that mediate self-oligomerization. The analysis of the binding sites of MDP and imidazoquinoline revealed LRR 7-9 to be the most energetically favored site for imidazoquinoline interaction. However, the binding free energy calculations using the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) method indicated that MDP is incompatible for activating the Nalp3 molecule in its monomeric form, and suggest its complex interaction with NOD2 or other NLRs accounts for MDP recognition. The high binding affinity of ATP with NACHT was correlated to the experimental data for human NLRs. Our binding site prediction for imidazoquinoline in LRR warrants further investigation via in vivo models. This is the first study that provides ligand recognition in mouse Nalp3 and its spatial structural arrangements.

LRRsearch: An asynchronous server-based application for the prediction of leucine-rich repeat motifs and an integrative database of NOD-like receptors

The leucine-rich repeat (LRR) motifs of the nucleotide-binding oligomerization domain like receptors (NLRs) play key roles in recognizing and binding various pathogen associated molecular patterns (PAMPs) resulting in the activation of downstream signaling and innate immunity. Therefore, identification of LRR motifs is very important to study ligand-receptor interaction. To date, available resources pose restrictions including both false negative and false positive prediction of LRR motifs from the primary protein sequence as their algorithms are relied either only on sequence based comparison or alignment techniques or are over biased for a particular LRR containing protein family. Therefore, to minimize the error (≤5%) and to identify a maximum number of LRR motifs in the wide range of proteins, we have developed "LRRsearch" web-server using position specific scoring matrix (PSSM) of 11 residue LRR-HCS (highly conserved segment) which are frequently observed motifs in the most divergent classes of LRR containing proteins. A data library of 421 proteins, distributed among five known NLR families has also been integrated with the "LRRsearch" for the rich user experience. The access to the "LRRsearch" program is freely available at http://www.lrrsearch.com/.

Exploration of the binding modes of buffalo PGRP1 receptor complexed with meso-diaminopimelic acid and lysine-type peptidoglycans by molecular dynamics simulation and free energy calculation

The peptidoglycan recognition proteins (PGRPs) are the key components of innate-immunity, and are highly specific for the recognition of bacterial peptidoglycans (PGN). Among different mammalian PGRPs, the PGRP1 binds to murein PGN of Gram-positive bacteria (lysine-type) and also have bactericidal activity towards Gram-negative bacteria (diaminopimelic acid or Dap-type). Buffaloes are the major sources of milk and meat in Asian sub-continents and are highly exposed to bacterial infections. The PGRP activates the innate-immune signaling, but their studies has been confined to limited species due to lack of structural and functional information. So, to understand the structural constituents, 3D model of buffalo PGRP1 (bfPGRP1) was constructed and conformational and dynamics properties of bfPGRP1 was studied. The bfPGRP1 model highly resembled human and camel PGRP structure, and shared a highly flexible N-terminus and centrally placed L-shaped cleft. Docking simulation of muramyl-tripeptide, tetrapeptide, pentapeptide-Dap-(MTP-Dap, MTrP-Dap and MPP-Dap) and lysine-type (MTP-Lys, MTrP-Lys and MPP-Lys) in AutoDock 4.2 and ArgusLab 4.0.1 anticipated β1, α2, α4, β4, and loops connecting β1-α2, α2-β2, β3-β4 and α4-α5 as the key interacting domains. The bfPGRP1-ligand complex molecular dynamics simulation followed by free binding energy (BE) computation conceded BE values of -18.30, -35.53, -41.80, -25.03, -24.62 and -22.30 kJ mol(-1) for MTP-Dap, MTrP-Dap, MPP-Dap, MTP-Lys, MTrP-Lys and MPP-Lys, respectively. The groove-surface and key binding residues involved in PGN-Dap and Lys-type interaction intended by the molecular docking, and were also accompanied by significant BE values directed their importance in pharmacogenomics, and warrants further in vivo studies for drug targeting and immune signaling pathways exploration.