Nucleotide binding and oligomerization domain (NOD2) is a key component of innate immunity that is highly specific for muramyl dipeptide (MDP)-a peptidoglycan component of bacterial cell wall. MDP recognition by NOD2-leucine rich repeat (LRR) domain activates NF-κB signaling through a protein-protein interaction between caspase activating and recruitment domains (CARDs) of NOD2 and downstream receptor interacting and activating protein kinase 2 (RIP2). Due to the lack of crystal/NMR structures, MDP recognition and CARD-CARD interaction are poorly understood. Herein, we have predicted the probable MDP and CARD-CARD binding surfaces in zebrafish NOD2 (zNOD2) using various in silico methodologies. The results show that the conserved residues Phe819, Phe871, Trp875, Trp929, Trp899, and Arg845 located at the concave face of zNOD2-LRR confer MDP recognition by hydrophobic and hydrogen bond (H-bond) interactions. Molecular dynamics simulations reveal a stable association between the electropositive surface on zNOD2-CARDa and the electronegative surface on zRIP2-CARD reinforced mostly by H-bonds and electrostatic interactions. Importantly, a 3.5 Å salt bridge is observed between Arg60 of zNOD2-CARDa and Asp494 of zRIP2-CARD. Arg11 and Lys53 of zNOD2-CARDa and Ser498 and Glu508 of zRIP2-CARD are critical residues for CARD-CARD interaction and NOD2 signaling. The 2.7 Å H-bond between Lys104 of the linker and Glu508 of zRIP2-CARD suggests a possible role of the linker for shaping CARD-CARD interaction. These findings are consistent with existing mutagenesis data. We provide first insight into MDP recognition and CARD-CARD interaction in the zebrafish that will be useful to understand the molecular basis of NOD signaling in a broader perspective.

The nucleotide-binding and oligomerization domain (NOD)-containing protein 1 (NOD1) plays the pivotal role in host-pathogen interface of innate immunity and triggers immune signalling pathways for the maturation and release of pro-inflammatory cytokines. Upon the recognition of iE-DAP, NOD1 self-oligomerizes in an ATP-dependent fashion and interacts with adaptor molecule receptor-interacting protein 2 (RIP2) for the propagation of innate immune signalling and initiation of pro-inflammatory immune responses. This interaction (mediated by NOD1 and RIP2) helps in transmitting the downstream signals for the activation of NF-κB signalling pathway, and has been arbitrated by respective caspase-recruitment domains (CARDs). The so-called CARD-CARD interaction still remained contradictory due to inconsistent results. Henceforth, to understand the mode and the nature of the interaction, structural bioinformatics approaches were employed. MD simulation of modelled 1:1 heterodimeric complexes revealed that the type-Ia interface of NOD1CARD and the type-Ib interface of RIP2CARD might be the suitable interfaces for the said interaction. Moreover, we perceived three dynamically stable heterotrimeric complexes with an NOD1:RIP2 ratio of 1:2 (two numbers) and 2:1. Out of which, in the first trimeric complex, a type-I NOD1-RIP2 heterodimer was found interacting with an RIP2CARD using their type-IIa and IIIa interfaces. However, in the second and third heterotrimer, we observed type-I homodimers of NOD1 and RIP2 CARDs were interacting individually with RIP2CARD and NOD1CARD (in type-II and type-III interface), respectively. Overall, this study provides structural and dynamic insights into the NOD1-RIP2 oligomer formation, which will be crucial in understanding the molecular basis of NOD1-mediated CARD-CARD interaction in higher and lower eukaryotes.

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.

Viral infections are one of the major challenges in aquaculture production, and considered as the potential threat for fish farming. Toll-like receptor (TLR) 3 and TLR22 are highly specialized innate immune receptors that recognize double-stranded (ds)-RNA of viruses resulting in the induction of innate immunity. The existence of TLR3 and TLR22 only in aquatic animals indicates their distinctive characteristics in viral infection; however, the studies in exploring their structural features and dsRNA binding mechanism are still elusive. Here, we studied the structural and functional differentiations of TLR3 and TLR22 in zebrafish by employing comparative modeling and molecular dynamics simulation. Comparative structural analysis revealed a distinct spatial arrangement of TLR22 ectodomain with a flattened horseshoe-shape conformation as compared to other TLRs. Essential dynamics studies showed that unlike TLR3, TLR22 possessed a prominent motion, elasticity and twisting at both terminus separated by a distance equivalent to the length of a short-sized dsRNA. Interaction analysis of polyinosinic:polycytidylic acid (poly I:C) and dsRNA depicted leucine-rich-repeats (LRR)2-3 and LRR18-19 (in TLR3) and LRRNT-LRR3 and LRR22-24 (in TLR22) as the potential binding sites. The short-sized dsRNA binds tightly across its full-length with TLR22-monomer, and suggested that TLR22 dimer may sense long-sized dsRNA. Binding energy (BE) calculation using MM/PBSA method from the TLR3- and TLR22-ligand complexes revealed an adequate binding affinity between TLR22-monomer and dsRNA as like as TLR3-dimer-dsRNA complex. Mutagenesis and BE computation of key residues suggested their involvement in dsRNA recognition. These findings can be helpful for therapeutic applications against viral diseases in fish.

Nucleotide-binding and oligomerization domain (NOD)-like receptors (NLRs), the first line of defense, are the cytosolic pattern recognition receptors (PRRs) that regulate the inflammatory activity in response to invading pathogens. NLRs are the members of AAA+ ATPase superfamily that comprises of N-terminal EBD(s), a centrally positioned NOD/NACHT and varying range of LRRs towards the C-terminal end. Due to the lack of structural data, the functional aspects of NLRP-signaling mechanism, which includes pathogen recognition, nucleotide-binding, and sensor-adaptor-effector interactions, are not fully understood. In this study, we implemented structural bioinformatics approaches including protein modeling, docking, and molecular dynamics simulations to explore the structural-dynamic features of ADP-/ATP-Mg2+ binding in NLRPNACHT models. Our results indicate a similar mode of ATP-Mg2+ binding in all NLRPNACHT models and the interacting residues are found consistent with reported mutagenesis data. Accompanied by the key amino acids (proposed to be crucial for ATP-Mg2+ coordination), we further have noticed that some additional conserved residues (including ‘Trp’ of the PhhCW motif, and ‘Phe’ and ‘Tyr’ of the GFxxxxRxxYF motif) are potentially interacting with ATP during dynamics; which require further experimentation for legitimacy. Overall, this study will help in understanding the ADP-/ATP-Mg2+ binding mechanisms in NLRPs in a broader perspective and the proposed ATP-binding pocket will aid in designing novel inhibitors for the regulation of inflammasome activity.