Structural Elucidation of Inter-CARD Interfaces involved in NOD2 Tandem CARD Association and RIP2 Recognition

Nucleotide-binding and oligomerization domain-containing protein 2 (NOD2) recognizes the muramyl dipeptide and activates the NF-κB signaling cascade following its interaction with receptor-interacting protein 2 (RIP2) via caspase recruitment domains (CARDs). The NOD2-RIP2 interaction is not understood well due to inadequate structural information. Using comparative modeling and multimicrosecond timescale molecular dynamics simulations, we have demonstrated the association of NOD2-CARDs (CARDa-CARDb) and their interaction with RIP2^CARD. Our results suggest that a negatively charged interface of NOD2^CARDa and positively charged type-Ia interface of NOD2^CARDb are crucial for CARDa-CARDb association and the type-Ia interface of NOD2^CARDa and type-Ib interface of RIP2^CARD predicted to be involved in 1:1 CARD-CARD interaction. Moreover, the direct interaction of NOD2^CARDb with RIP2^CARD signifies the importance of both CARDs of NOD2 in RIP2-mediated CARD-CARD interaction. Altogether, the structural results could help in understanding the underlying molecular details of the NOD2-RIP2 association in higher and lower eukaryotes.

Molecular and functional characterization of spotted snakehead NOD1 with an emphasis on structural insights into iE-DAP binding motifs employing advanced bioinformatic tools

Nucleotide-binding and oligomerization domain (NOD)-like receptors (NLRs) are cytosolic receptors implicated in recognition of intracellular pathogen associated molecular patterns (PAMPs) and danger associated molecular patterns (DAMPs). Depending upon their effector binding domain (EBD) at the C-terminal, the NLRs are categorized into NLRA, NLRB, NLRC, NLRP and NLRX. NOD1 is a pivotal player in immune responses against bacterial and viral invasions and interacts with pathogens via C-terminal leucine rich repeat (LRR) domain. This study aims at characterizing NOD1 in an economically important teleost of the Indian subcontinent, spotted snakehead Channa punctata. The understanding of pathogen-receptor interaction in teleosts is still obscure. In light of this, combinatorial approach involving protein modeling, docking, MD simulation and binding free energy calculation were employed to identify key motifs involved in binding iE-DAP. In silico analysis revealed that NOD1 consists of 943 amino acids comprising of one caspase recruitment domain (CARD) at N-terminal, one central NACHT domain and nine leucine rich repeat (LRR) regions at C-terminal. Structural dynamics study showed that the C-terminal β-sheet LRR4-7 region is involved in iE-DAP binding. NOD1 was ubiquitously and constitutively expressed in all tissues studied. Differential expression profile of NOD1 induced by Aeromonas hydrophila infection was also investigated. Lymphoid organs and phagocytes of infected spotted snakehead showed significant downregulation of NOD1 expression. The current study thus gives an insight into structural and functional dynamics of NOD1 which might have future prospect for structure-based drug designing in teleosts.Communicated by Ramaswamy H. Sarma.

Recognition of DAP and activation of NF-κB by cytosolic sensor NOD1 in Oreochromis niloticus

As a lower vertebrate, the immune defense mechanism of fish mainly depends on the innate immune system. Nucleotide-binding oligomerization domain-like receptors (NLRs) are an important class of pattern recognition receptors in the innate immune system. In this study, NOD1 gene was cloned and characterized in Nile tilapia (Oreochromis niloticus). The ORF of Nile tilapia NOD1 gene was 2826 bp long and encoded 941 amino acid residues with a structure of CARD-NACHT-LRRs that was similar to the other counterparts in mammals and fishes. Phylogenetic and synteny analysis showed that NOD1 was conserved among different fishes and existed at least in the early stage of fish evolution. Expression pattern revealed that NOD1 mRNA was constitutively expressed in the tested tissues, while had high expression level in main immune organs and mucosal immune tissues (liver, head kidney, spleen, blood, gill, and intestine). Following Streptococcus agalactiae challenge, Nile tilapia NOD1 mRNA expression levels were altered in immune organs (liver, head kidney, spleen, blood), and the expression pattern was similar in liver, spleen and blood. Furthermore, the ligand recognition and signaling pathway of Nile tilapia NOD1 were also analyzed, it showed that NOD1 could recognize Tri-DAP intracellularly and activated NF-κB signaling pathway. In summary, our results indicated that the Nile tilapia NOD1 may play an important role in innate immune system and provided a basis for the functional study of NOD1 in teleost.

Insights into structure and dynamics of extracellular domain of Toll-like receptor 5 in Cirrhinus mrigala (mrigala): A molecular dynamics simulation approach

The toll-like receptor 5 (TLR5) is the most conserved important pattern recognition receptors (PRRs) often stimulated by bacterial flagellins and plays a major role in the first-line defense against invading pathogenic bacteria and in immune homeostasis. Experimental crystallographic studies have shown that the extracellular domain (ECD) of TLR5 recognizes flagellin of bacteria and functions as a homodimer in model organism zebrafish. However, no structural information is available on TLR5 functionality in the major carp Cirrhinus mrigala (mrigala) and its interaction with bacterial flagellins. Therefore, the present study was undertaken to unravel the structural basis of TLR5-flagellin recognition in mrigala using structural homodimeric TLR5-flagellin complex of zebrafish as reference. Integrative structural modeling and molecular dynamics simulations were employed to explore the structural and mechanistic details of TLR5 recognition. Results from structural snapshots of MD simulation revealed that TLR5 consistently formed close interactions with the three helices of the D1 domain in flagellin on its lateral side mediated by several conserved amino acids. Results from the intermolecular contact analysis perfectly substantiate with the findings of per residue-free energy decomposition analysis. The differential recognition mediated by flagellin to TLR5 in mrigala involves charged residues at the interface of binding as compared to the zebrafish complex. Overall our results shows TLR5 of mrigala involved in innate immunity specifically recognized a conserved site on flagellin which advocates the scientific community to explore host-specific differences in receptor activation.

Sodium valproate attenuates the iE-DAP induced inflammatory response by inhibiting the NOD1-NF-κB pathway and histone modifications in bovine mammary epithelial cells

The anti-inflammatory effects of sodium valproate (VPA) in vivo and in vitro have been demonstrated in recent studies. The aim of this study was to evaluate whether VPA can suppress inflammation in bovine mammary epithelial cells (BMECs) stimulated by γ-D-glutamyl-meso-diaminopimelic acid (iE-DAP). First, the concentration and treatment points of iE-DAP and VPA were optimized. Then, BMECs were cultured in complete media and separated into four groups: untreated control cells (CON group), cells stimulated by 10 μg/mL iE-DAP for 6 h (DAP group), cells stimulated by 0.5 mmol/L VPA for 6 h (VPA group), and cells pretreated with VPA (0.5 mmol/L) for 6 h followed by 10 μg/mL of iE-DAP for 6 h (VD group). The results showed that the level of  interleukin-6 (IL-6) and tumor necrosis factor-α (TNF-α) in the culture medium increased in the iE-DAP-treated cells and that pretreatment with VPA reversed this increase. iE-DAP increased both mRNA and protein expression levels of nucleotide-binding oligomerization domain-containing protein 1 (NOD1) and receptor-interacting protein kinas (RIPK2) and activated inhibitor of NF-κB (IκB) and nuclear factor-kappa B p65 (NF-κB p65) through phosphorylation. Upon activation of the NF-κB pathway, the expression of the pro-inflammatory cytokines IL-6, interleukin-8 (IL-8) and interleukin-1β (IL-1β), the acute phase protein serum amyloid A 3 (SAA3) and the lingual antimicrobial peptide (LAP) but not  haptoglobi (HP) or bovine neutrophil beta defensing 5 (BNBD5) were increased in the DAP group. The VPA pretreatment induced the acetylation of signal transducers and activators of transcription(STAT1) and histone 3 (H3) by inhibiting histone deacetylase (HDAC) and then suppressed the NF-κB pathway. Moreover, VPA induced autophagy and reduced apoptosis in BMECs in the VD group. These results suggested that VPA treatment can attenuate the inflammatory response induced by iE-DAP.

Cloning, characterization and ontogenetic expression profile of RAG-2 and IgM in angelfish (Pterophyllum scalare)

Early larval developmental stages of fish are highly susceptible to opportunistic pathogens until the complete maturation of the lymphoid organs. Knowledge of the expression pattern of important markers of adaptive immune system during the ontogenetic development is essential before vaccinating the fish. In the present study, Pterophyllum scalare (angelfish) was taken to explore the relative expression profile of developmental markers of adaptive immunity, recombination activating gene-2 (RAG-2) and immunoglobulin M (IgM). The fishes were bred and early developmental stages (0-45 days post-hatched) were used to assess the expression profile. The genes, RAG-2 and IgM were cloned and sequenced with the base pair lengths of 1958 bp and 225 bp respectively. The mRNA expression of RAG-2 appeared at insignificant level at the first day of hatching, but the expression was significantly increased from 24 dph (days post-hatching) onwards and reached its peak at 27 dph. The results proved that the maturation of lymphoid organs was completed at 27 dph as the respective protein is involved in the V(D)J recombination, important for the maturation of lymphoid organs. A similar trend was also observed in the mRNA transcript levels of IgM gene and a significantly high expression was detected from 27 dph onwards. The present study suggested that the suitable time for vaccination in P. scalare could be taken at 27 dph, as the maturation and development of lymphoid organs is completed thus helps in stimulating the adaptive response of immunity against any pathogen.

Molecular characterization, constitutive expression and GTP binding mechanism of Cirrhinus mrigala (Hamilton, 1822) Myxovirus resistance (Mx) protein

Myxovirus resistance (Mx) proteins represents the subclass of the dynamin superfamily of large Guanosine triphosphates (GTPases), play esential role in intracellular vesicle trafficking, endocytosis, organelle homeostasis and mitochondria distribution. These proteins are key players of the vertebrate immune system, induced by type-I and type-III interferons (IFN) of infected host and inhibit viral replication by sequestering its nucleoprotein. In the present study, we report the sequencing and characterization of Cirrhinus mrigala Mx protein (CmMx) for the first time and observed its constitutive expression in different tissues for a period of fourteen days. The synthetic peptide, LSGVALPRGTGI, was dissolved in PBS and injected into a rabbit and the antibody raised against CmMx was used to study the level of its expression. The full length of the CmMx cDNA is 2244 bp with a molecular mass of 70.9 kDa and a predicted isoelectric point of 8.25. The 627 amino acids polypeptide formed of three main functional domains: N-terminal GTPase domain (GD), a middle domain (MD) and GTPase effector domain (GED) with carboxy terminal leucine zipper motif. The 3D models of CmMx protein was modeled based on available close structural homologs and further validated through molecular dynamics (MD) simulations. MD study revealed the importance of G-domain responsible for recognition of GTP, which perfectly corroborate with earlier studies. MM/PBSA binding free energy analysis displayed that van der Waals and electrostatic energy were the key driving force behind molecular recognition of GTP by CmMx protein. The results from this study will illuminate more lights into the ongoing research on myxovirus resistance protein and its role in inhibition of viral replication in other eukaryotic system as well.

Bioactivity guided isolation of antidiabetic and antioxidant compound from Xylocarpus granatum J. Koenig bark

The present study was designed to identify antidiabetic and antioxidant constituents from ethanol bark extract of a medicinally important mangrove plant Xylocarpus granatum J. Koenig, using in vitro bioactivity-guided fractionation. The repeated fractionation of X. granatum ethanol bark extract (XGEB) by silica gel column chromatography yielded a compound with strong antidiabetic and antioxidant potential. The bioactive compounds likely to be present in the XGEB fraction were identified by FT-IR, ¹H & ¹³C NMR and MS analysis and determined as a limonoid derivative Xyloccensin-I, by comparing spectral data with the literature reports. The isolated compound demonstrated excellent in vitro antidiabetic potential IC₅₀ values of 0.25 and 0.16 mg/ml, respectively for α-amylase and α-glucosidase inhibition study. The antioxidant potential assayed by DPPH, ABTS, superoxide and hydrogen peroxide scavenging studies exhibited that the isolated compound could scavenge these free radicals with IC₅₀ values of 0.041, 0.039, 0.096 and 0.235 mg/ml, respectively. Further, in silico study was performed to find the antidiabetic activity of Xyloccensin-I by docking it against α-glucosidase enzyme. The study demonstrated that Xyloccensin-I have satisfactory interactions and binding energies when docked into target which further confirms the possible mode of antidiabetic action of the isolated compound. The bioactivity assays further asserts the antidiabetic and antioxidant efficacy of the isolated compound which strongly suggests that Xyloccensin-I holds promise in the pharmaceutical industry. The results from this study provide new mechanistic evidence justifying, at least in part, the traditional use of X. granatum extract for antidiabetic and antioxidants activities.

Molecular cloning, GTP recognition mechanism and tissue-specific expression profiling of myxovirus resistance (Mx) protein in Labeo rohita (Hamilton) after Poly I:C induction

The myxovirus resistance (Mx) proteins belong to interferon-induced dynamin GTPase and play pivotal role in the inhibition of replication of numerous viruses. These antiviral proteins are released in usual or diseased condition to prevent the viral attack and to carry regular cellular activities like endocytosis and trafficking of nucleoproteins into the nucleus. The invasion of virus up-regulates the expression of Mx transcripts and double-stranded RNA mimic like polyinosinic polycytidyilic acid (Poly I:C). To understand the tissue-specific expression profiling and mechanism of GTP recognition of Mx protein from Labeo rohita (rohu), the full-length gene was cloned, sequenced and characterized through various Bioinformatics tools for the first time. The Mx cDNA was comprised of 2297 bp, and the open reading frame of 1938 bp encodes polypeptide of 631 amino acids. The coding sequence of Mx protein possess the signature motif of dynamin superfamily, LPRG(S/K)GIVTR, the tripartite guanosine-5/triphosphate (GTP)-binding motif (GXXXSGKS/T, DXXG and T/NKXD) and the leucine zipper motifs at the C-terminal end, well conserved in all interferon-induced Mx protein in vertebrates. Western blotting confirmed the molecular weight of Mx protein to be 72 kDa. After the intraperitoneal challenge of L. rohita with a Poly I:C, up-regulation of Mx protein was observed in brain, spleen, liver, kidney, intestine, heart, muscle, and gill. Ontogeny study displayed pronounced expression of Mx protein in all stages of the developmental of Rohu after Poly I:C induction. However a persistent expression of Mx transcript was also observed in Rohu egg as well as milt without induction with Poly I:C. Higher expression of Mx gene was observed on 96 h where it was 6.4 folds higher than the control. The computational modelling of Mx protein portrayed the tripartite N-terminal G-domain that binds to GTP, the bundle-signaling element (BSE) which interconnects the G-domain to the elongated stalk domain and C-terminal helical stalk domain. In agreement with the experimental studies, a series of conserved residues viz., Gln52, Ser53, Ser54, Leu68, Pro69, Gly71, Gly73, Thr76, Asp151, Gly154, Thr220, Lys221, Val251, Cys253, Arg254, and Gly255 were computed to be indispensable for tight anchoring of GTP within binding cavity of G-domain. The binding free energy calculation study depicted that the van der Waals and electrostatic terms contributs significantly to molecular recognition of GTP. Collectively, our study provides mechanistic insights into the tissue-specific expression profiling and GTP binding mechanism of Mx protein from Labeo rohita, which is expected to drive further research on several cellular events including viral resistance and endocytosis in the near future.

Mitogen-activated protein kinase kinase 6 is involved in the immune response to bacterial di-/tripeptide challenge in grass carp Ctenopharyngodon idella

Mitogen-activated protein kinase kinase 6 (MKK6) is an essential component of the p38MAPK signaling pathway, which is involved in the modulation of inflammation, cell apoptosis and survival responses in mammals. However, the function of MKK6s in teleosts is still unclear. In this study, a fish MKK6 homolog (CiMKK6) was first identified from the grass carp (Ctenopharyngodon idella), a freshwater fish. CiMKK6 cDNA encodes a putative protein of 357 amino acids that contains conserved structural characteristics of the MKK6 family, including the S_TKc domain, SVAKT motif and DVD site. The deduced CiMKK6 protein exhibits high sequence homology with other reported fish MKK6s and shares the closest relationship with MKK6 from Danio rerio. Quantitative real-time PCR (qRT-PCR) analysis revealed that CiMKK6 mRNA was widely expressed in all tested tissues and stages of embryonic development. Additionally, the transcript levels of CiMKK6 in the intestine were significantly upregulated in response to bacterial muramyl dipeptide (MDP) and L-Ala-γ-D-Glu-meso-diaminopimelic acid (Tri-DAP) stimulation. Moreover, subcellular localization analysis indicated that CiMKK6 was distributed in both the cytoplasm and the nucleus of HEK293T cells. Finally, overexpression of CiMKK6 significantly enhanced the transcriptional activity of the AP-1 reporter gene in HEK293T cells. Overall, these findings may help better clarify the immune function of teleost MKK6s and provide new insight into the immune defense mechanisms of grass carp.

Deciphering the ATP-binding mechanism(s) in NLRP-NACHT 3D models using structural bioinformatics approaches

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.

Structural bioinformatics insights into ATP binding mechanism in zebrafish (Danio rerio) cyclin-dependent kinase-like 5 (zCDKL5) protein

In mammalian systems, the conserved cyclin-dependent protein kinases (CDKs) control the process of cell division and curb the transcription mechanism in response to diverse signaling events that are essential for the catalytic activity. In zebrafish, zCDKL5 portrays differential expression profiling in several tissues and presumed to play a vital role in the neuronal development. In this present study, the sequence-structure relationship and mode of ATP binding in zCDKL5 was unveiled through theoretical modeling, molecular docking, and MD simulations. Like human CDKs, the modeled zCDKL5 was found to be bipartite in nature, where, ATP binds to the central cavity of the catalytic domain through a strong network of H-bonding, electrostatic, and hydrophobic interactions. MD simulation portrayed that conserved residues, viz, Ile10, Gly11, Glu12, Val18, Val64, Glu81, Cys143, and Asp144 were indispensable for tight anchoring of ATP and contribute to the stability of the zCDKL5-ATP complex. MM/PBSA binding free energy analysis displayed that van der Waal energy (ΔG _vwd ) and Electrostatic energy (ΔG _ele ) were the major contributors towards the overall binding free energy. Thus, the comparative structural bioinformatics approach has shed new insights into the dynamics and ATP binding mechanism of zCDKL5. The results from the study will help to undertake further research on the role of phosphorylated CDKL5 in the onset of neurodevelopmental disorders caused by mutations in higher eukaryotic systems.

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.

Elucidating the interfaces involved in CARD-CARD interactions mediated by NLRP1 and Caspase-1 using molecular dynamics simulation

Inflammasomes are the multi-protein caspase-activating complexes mainly assembled by the sensor proteins (NLRs/ALRs), adaptor molecule (ASC) and effector molecule pro-caspase-1 for the production and release of proinflammatory cytokines, IL-1β and IL-18. NLRP1 is the first NLR known to assemble the multi-protein complex. Unlike NLRP3, NLRP1 has an additional effector binding domain (CARD) at the carboxyl-terminal, which is reported to interact with pro-caspase-1 (precluding the recruitment of ASC) for the transmission of danger signals. So far no direct interaction has been observed between the NLRP1 and CASP1 at the structural level. In this study, an attempt has been made to elucidate the possible mode of interaction(s) between CASP1 and NLRP1 CARDs using structural bioinformatics approaches. The results revealed that the type-Ia patch of CASP1^CARD (R10, K11, and R55) is probably the favorable interface for 1:1 interaction. Moreover, the interactions mediated in the type-II and/(or) type-III interfaces of counter CARDs can also be not ruled out altogether. Overall, the findings of this study can be beneficial in understanding the underlying molecular mechanism(s) associated with NLRP1-mediated inflammasome.

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.

POP1 might be recruiting its type-Ia interface for NLRP3-mediated PYD-PYD interaction: Insights from MD simulation

Inflammasomes are multiprotein caspase-activating complexes that enhance the maturation and release of proinflammatory cytokines (IL-1β and IL-18) in response to the invading pathogen and/or host-derived cellular stress. These are assembled by the sensory proteins (viz NLRC4, NLRP1, NLRP3, and AIM-2), adaptor protein (ASC), and effector molecule procaspase-1. In NLRP3-mediated inflammasome activation, ASC acts as a mediator between NLRP3 and procaspase-1 for the transmission of signals. A series of homotypic protein-protein interactions (NLRP3^PYD :ASC^PYD and ASC^CARD :CASP1^CARD ) propagates the downstream signaling for the production of proinflammatory cytokines. Pyrin-only protein 1 (POP1) is known to act as the regulator of inflammasome. It modulates the ASC-mediated inflammasome assembly by interacting with pyrin domain (PYD) of ASC. However, despite similar electrostatic surface potential, the interaction of POP1 with NLRP3^PYD is obscured till date. Herein, to explore the possible PYD-PYD interactions between NLRP3^PYD and POP1, a combined approach of protein-protein docking and molecular dynamics simulation was adapted. The current study revealed that POP1's type-Ia interface and type-Ib interface of NLRP3^PYD might be crucial for 1:1 PYD-PYD interaction. In addition to type-I mode of interaction, we also observed type-II and type-III interaction modes in two different dynamically stable heterotrimeric complexes (POP1-NLRP3-NLRP3 and POP1-NLRP3-POP1). The inter-residual/atomic distance calculation exposed several critical residues that possibly govern the said interaction, which need further investigation. Overall, the findings of this study will shed new light on hitherto concealed molecular mechanisms underlying NLRP3-mediated inflammasome, which will have strong future therapeutic implications.

Bacterial peptidoglycan with amidated meso-diaminopimelic acid evades NOD1 recognition: an insight into NOD1 structure–recognition

Nucleotide-binding oligomerization domain-containing protein 1 (NOD1) is an intracellular pattern recognition receptor that recognizes bacterial peptidoglycan (PG) containing meso-diaminopimelic acid (mesoDAP) and activates the innate immune system. Interestingly, a few pathogenic and commensal bacteria modify their PG stem peptide by amidation of mesoDAP (mesoDAPNH2). In the present study, NOD1 stimulation assays were performed using bacterial PG containing mesoDAP (PGDAP) and mesoDAPNH2 (PGDAPNH2) to understand the differences in their biomolecular recognition mechanism. PGDAP was effectively recognized, whereas PGDAPNH2 showed reduced recognition by the NOD1 receptor. Restimulation of the NOD1 receptor, which was initially stimulated with PGDAP using PGDAPNH2, did not show any further NOD1 activation levels than with PGDAP alone. But the NOD1 receptor initially stimulated with PGDAPNH2 responded effectively to restimulation with PGDAP. The biomolecular structure–recognition relationship of the ligand-sensing leucine-rich repeat (LRR) domain of human NOD1 (NOD1–LRR) with PGDAP and PGDAPNH2 was studied by different computational techniques to further understand the molecular basis of our experimental observations. The d-Glu–mesoDAP motif of GMTPDAP, which is the minimum essential motif for NOD1 activation, was found involved in specific interactions at the recognition site, but the interactions of the corresponding d-Glu–mesoDAP motif of PGDAPNH2 occur away from the recognition site of the NOD1 receptor. Hot-spot residues identified for effective PG recognition by NOD1–LRR include W820, G821, D826 and N850, which are evolutionarily conserved across different host species. These integrated results thus successfully provided the atomic level and biochemical insights on how PGs containing mesoDAPNH2 evade NOD1–LRR receptor recognition.

InteractoMIX: a suite of computational tools to exploit interactomes in biological and clinical research

Virtually all the biological processes that occur inside or outside cells are mediated by protein–protein interactions (PPIs). Hence, the charting and description of the PPI network, initially in organisms, the interactome, but more recently in specific tissues, is essential to fully understand cellular processes both in health and disease. The study of PPIs is also at the heart of renewed efforts in the medical and biotechnological arena in the quest of new therapeutic targets and drugs. Here, we present a mini review of 11 computational tools and resources tools developed by us to address different aspects of PPIs: from interactome level to their atomic 3D structural details. We provided details on each specific resource, aims and purpose and compare with equivalent tools in the literature. All the tools are presented in a centralized, one-stop, web site: InteractoMIX (http://interactomix.com).

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.