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.

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.

Nucleotide binding and oligomerization domain 1 (NOD1), a cytoplasmic pattern recognition receptor (PRR) and is a key component for modulating innate immunity and signaling. It is highly specific to γ-D-Glu-mDAP (iE-DAP), a cell wall component of Gram-negative and few Gram-positive bacteria. In the absence of the experimental structure of NOD1 leucine rich repeat (NOD1-LRR) domain, the NOD signaling cascade mediated through NOD1 and iE-DAP interaction is poorly understood. Herein, we modeled 3D structure of zebrafish NOD1-LRR (zNOD1-LRR) through a protein-threading approach and structural integrity of the model was assessed using molecular dynamics simulations. Molecular interaction analysis of iE-DAP and zNOD1-LRR, their complex stability and binding free energy studies were conducted to anticipate the ligand binding residues in zNOD1. Our study revealed that His775, Lys777, Asp803, Gly805, Trp807, Asn831, Ser833, Ile859 and Trp861 situated in the β-sheet region of zNOD1-LRR could be involved in iE-DAP recognition, which correlates the earlier findings in human. Comparison of binding free energies of native and mutant zNOD1-iE-DAP complexes delineated His775, Lys777, Asp803, Ser833 and Ile859 as the pivotal residues for energetic stability of NOD1 and iE-DAP interaction. This study provides the first comprehensive description of biophysical and biochemical parameters responsible for NOD1 and iE-DAP interaction in zebrafish, which is expected to shed more light on NOD1 signaling and therapeutic applications in other organisms.