PRDM9 is the sole hybrid sterility gene identified so far in vertebrates. PRDM9 gene encodes a protein with an immensely variable zinc-finger (ZF) domain that determines the site of meiotic recombination hotspots genome-wide. In this study, the terminal ZF domain of PRDM9 on bovine chromosome 1 and its paralog on chromosome 22 were characterized in 225 samples from five ruminant species (cattle, yak, mithun, sheep and goat). We found extraordinary variation in the number of PRDM9 zinc fingers (6 to 12). We sequenced PRDM9 ZF encoding region from 15 individuals (carrying the same ZF number in both copies) and found 43 different ZF domain sequences. Ruminant zinc fingers of PRDM9 were found to be diversifying under positive selection and concerted evolution, specifically at positions involved in defining their DNA-binding specificity, consistent with the reports from other vertebrates such as mice, humans, equids and chimpanzees. ZF-encoding regions of the PRDM7, a paralog of PRDM9 on bovine chromosome 22 and on unknown chromosomes in other studied species were found to contain 84 base repeat units as in PRDM9, but there were multiple disruptive mutations after the first repeat unit. The diversity of the ZFs suggests that PRDM9 may activate recombination hotspots that are largely unique to each ruminant species.

Middle East respiratory syndrome coronavirus (MERS-CoV) has evolved to navigate through the sophisticated network of a host's immune system. The immune evasion mechanism including type 1 interferon and protein kinase R-mediated antiviral stress responses has been recently attributed to the involvement of MERS-CoV protein 4a (p4a) that masks the viral dsRNA. However, the structural mechanism of how p4a recognizes and establishes contacts with dsRNA is not well explained. In this study, we report a dynamic mechanism deployed by p4a to engage the viral dsRNA and make it unavailable to the host immune system. Multiple variants of p4a-dsRNA were created and investigated through extensive molecular dynamics procedures to highlight crucial interfacial residues that may be used as potential pharmacophores for future drug development. The structural analysis revealed that p4a exhibits a typical αβββα fold structure, as found in other dsRNA-binding proteins. The α1 helix and the β1-β2 loop play a crucial role in recognizing and establishing contacts with the minor grooves of dsRNA. Further, mutational and binding free energy analyses suggested that in addition to K63 and K67, two other residues, K27 and W45, might also be crucial for p4a-dsRNA stability.

Direct inhibition of tumor necrosis factor-alpha (TNF-α) action is considered a promising way to prevent or treat TNF-α-associated diseases. The trimeric form of TNF-α binds to its receptor (TNFR) and activates the downstream signaling pathway. The interaction of TNF-α with molecular-grade dimethyl sulfoxide (DMSO) in an equal volumetric ratio renders TNF-α inert, in this state, TNF-α fails to activate TNFR. Here, we aimed to examine the inhibition of TNF-α function by various concentrations of DMSO. Its higher concentration led to stronger attenuation of TNF-α-induced cytokine secretion by fibroblasts, and of their death. We found that this inhibition was mediated by a perturbation in the formation of the functional TNF-α trimer. Molecular dynamics simulations revealed a transient interaction between DMSO molecules and the central hydrophobic cavity of the TNF-α homodimer, indicating that a brief interaction of DMSO with the TNF-α homodimer may disrupt the formation of the functional homotrimer. We also found that the sensitizing effect of actinomycin D on TNF-α-induced cell death depends upon the timing of these treatments and on the cell type. This study will help to select an appropriate concentration of DMSO as a working solvent for the screening of water-insoluble TNF-α inhibitors.

The Anti-Müllerian Hormone (AMH) is a member of the transforming growth factor beta (TGF-β) superfamily, playing a significant role in cell proliferation, differentiation and apoptosis. In females, AMH is secreted throughout their reproductive life span from ovaries, whereas in males it is secreted by gonadal cells at a very early stage of testicular development. AMH is a promising marker of ovarian reserve in women and can be used to measure the female reproductive lifespan. In the present study, we cloned and sequenced the GC rich AMH gene from Indian riverine buffalo (Bubalus bubalis) and goat (Capra hircus). Obtained sequences were compared to the AMH sequences of other mammals, and corresponding amino acid sequences revealed that the caprine and bovine AMH sequences are more closely related to each other than to those of other mammals. Furthermore, we analyzed the chromosomal localization of AMH genes in mammalian species to understand potential syntenic relationship. The AMH gene is localized between the sequences for the SF3A and JSRP1 genes and maintains this precise location in relation to other nearby genes. The dN/dS ratio of AMH gene did not indicate any pressure for either positive or negative selection; thus, the physiological function of the AMH gene in the reproduction of these two ruminant species remains very vital. Similar to other mammals, the AMH gene may be an important indicator for regulating female reproductive biology function in bovine, cetacean, caprine, and camelidae.

β-defensins are adsorbable on the sperm surface in the male reproductive tract (MRT) and enhance sperm functional characteristics. The beta-defensin 129 (DEFB129) antimicrobial peptide is involved in sperm maturation, motility, and fertilization. However, its role in bovine fertility has not been well investigated. This study examines the relationship between the bovine BBD129 gene and Bos indicus x Bos taurus bull fertility. The complete coding sequence of BBD129 mRNA was identified by RNA Ligase Mediated-Rapid Amplification of cDNA End (RLM-RACE) and Sanger sequencing methodologies. It consisted of 582 nucleotides (nts) including 5' untranslated region (UTR) (46nts) and 3'UTR (23nts). It conserves all beta-defensin-like features. The expression level of BBD129 was checked by RT-qPCR and maximal expression was detected in the corpus-epididymis region compared to other parts of MRT. Polymorphism in BBD129 was also confirmed by Sanger sequencing of 254 clones from 5 high fertile (HF) and 6 low fertile (LF) bulls at two positions, 169 T > G and 329A > G, which change the S57A and N110S in the protein sequence respectively. These two mutations give rise to four types of BBD129 haplotypes. The non-mutated TA-BBD129 (169 T/329A) haplotype was substantially more prevalent among high-fertile bulls (P < 0.005), while the double-site mutated GG-BBD129 (169 T > G/329A > G) haplotype was significantly more prevalent among low-fertile bulls (P < 0.005). The in silico analysis confirmed that the polymorphism in BBD129 results in changes in mRNA secondary structure, protein conformations, protein stability, extracellular-surface availability, post-translational modifications (O-glycosylation and phosphorylation), and affects antibacterial and immunomodulatory capabilities. In conclusion, the mRNA expression of BBD129 in the MRT indicates its region-specific dynamics in sperm maturation. BBD129 polymorphisms were identified as the deciding elements accountable for the changed proteins with impaired functionality, contributing to cross-bred bulls' poor fertility.

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 (NOD)-like receptors (NLRs) are innate immune receptors that recognize bacterial cell wall components and initiate host immune response. Structure and function of NLRs have been well studied in human and mice, but little information exists on genetic composition and role of these receptors in innate immune system of water buffalo--a species known for its exceptional disease resistance. Here, a comparative study on the functional domains of NOD1 and NOD2 was performed across different species. The NOD mediated in-vitro cellular responses were studied in buffalo peripheral blood mononuclear cells, resident macrophages, mammary epithelial, and fibroblast cells. Buffalo NOD1 (buNOD1) and buNOD2 showed conserved domain architectures as found in other mammals. The domains of buNOD1 and buNOD2 showed analogy in secondary and tertiary conformations. Constitutive expressions of NODs were ubiquitous in different tissues. Following treatment with NOD agonists, peripheral lymphocytes showed an IFN-γ response along-with production of pro-inflammatory cytokines. Alveolar macrophages and mammary epithelial cells showed NOD mediated in-vitro immune response through NF-κB dependent pathway. Fibroblasts showed pro-inflammatory cytokine response following agonist treatment. Our study demonstrates that both immune and non-immune cells could generate NOD-mediated responses to pathogens though the type and magnitude of response depend on the cell types. The structural basis of ligand recognition by buffalo NODs and knowledge of immune response by different cell types could be useful for development of non-infective innate immune modulators and next generation anti-inflammatory compounds.

ALRs (AIM2-like Receptors) are germline encoded PRRs that belong to PYHIN gene family of cytokines, which are having signature N-terminal PYD (Pyrin, PAAD or DAPIN) domain and C-terminal HIN-200 (hematopoietic, interferon-inducible nuclear protein with 200 amino acid repeat) domain joined by a linker region. The positively charged HIN-200 domain senses and binds with negatively charged phosphate groups of single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA) purely through electrostatic attractions. On the other hand, PYD domain interacts homotypically with a PYD domain of other mediators to pass the signals to effector molecules downwards the pathways for inflammatory responses. There is remarkable inter-specific diversity in the numbers of functional PYHIN genes e.g. one in cow, five in human, thirteen in mice etc., while there is a unique loss of PYHIN genes in the bat genomes which was revealed by Ahn et al. (2016) by studying genomes of ten different bat species belonging to sub-orders yinpterochiroptera and yangochiroptera. The conflicts between host and pathogen interfaces are compared with "Red queen's arms race" which is also described as binding seeking dynamics and binding avoidance dynamics. As a result of this never-ending rivalry, eukaryotes developed PRRs as antiviral mechanism while viruses developed counter mechanisms to evade host immune defense. The PYHIN receptors are directly engaged with pathogenic molecules, so these should have evolved under the influence of selection pressures. In the current study, we investigated the nature of selection pressure on different domain types of IFI16-like (IFI16-L) PYHIN genes in ruminants.

Goat milk, a choice of substitution to mother's milk for its composition, fulfils nutritional requirement of infants, pregnant mothers and older people. The present study was carried out to unravel the milk proteome profiles from geographically and genetically diverse goat breeds by gel based 2DE and nLC-MS/MS. A total of 1307 functional proteins comprising casein and other low abundance proteins were identified. Gene annotations revealed that the majority of the proteins were involved in binding function, catalytic activity and structural molecules and localised in nucleus and membrane. The distinguished proteins were involved in 144 KEGG pathways in information processing, metabolism, cellular process, organismal systems and diseases. The large number of proteins and peptides including bioactive peptides were reported from goat milk from diverse agro-climatic regions of India indicating their significant potential for human health applications. SIGNIFICANCE: Goat milk in India is used in various Ayurvedic formulations to treat a number of ailments and allergies as well as for nutraceutical formulations. The study identifies milk protein variants both at protein and DNA level and subsequent identification of proteins by 2DE and nLC-MS/MS resulting in a proteome comprising of 1307 proteins. The specific proteins and peptides having significant role in immune regulation, disease pathways, cellular growth and metabolism have been identified. The results contribute to goat milk protein and peptide database which is very limited. We identified proteins for specific functional categories and associated them with different pathways for studying functional diversity of goat milk proteins. The proteins and peptides identified can be used for multiple human health application.

Toll-like receptors (TLRs) play a fundamental role in the inflammatory response against invading pathogens. However, the dysregulation of TLR-signaling pathways is implicated in several autoimmune/inflammatory diseases. Here, we show that a novel small molecule TLR-inhibitor (TAC5) and its derivatives TAC5-a, TAC5-c, TAC5-d, and TAC5-e predominantly antagonized poly(I:C) (TLR3)-, imiquimod (TLR7)-, TL8-506 (TLR8)-, and CpG-oligodeoxynucleotide (TLR9)-induced signaling pathways. TAC5 and TAC5-a significantly hindered the activation of nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB), reduced the phosphorylation of mitogen-activated protein kinases, and inhibited the secretion of tumor necrosis factor-α (TNF-α) and interleukin-6. Besides, TAC5-a prevented the progression of psoriasis and systemic lupus erythematosus (SLE) in mice. Interestingly, TAC5 and TAC5-a did not affect Pam3CSK4 (TLR1/2)-, FSL-1 (TLR2/6)-, or lipopolysaccharide (TLR4)-induced TNF-α secretion, indicating their specificity towards endosomal TLRs (TLR3/7/8/9). Collectively, our data suggest that the TAC5 series of compounds are potential candidates for treating autoimmune diseases such as psoriasis or SLE.

A single nucleotide substitution in intron 3 of IGF2 in pigs abrogates a binding site for a repressor and leads to a 3-fold up-regulation of IGF2 in skeletal muscle. The mutation has major effects on muscle growth, size of the heart, and fat deposition. Here, we have identified the repressor and find that the protein, named ZBED6, is previously unknown, specific for placental mammals, and derived from an exapted DNA transposon. Silencing of Zbed6 in mouse C2C12 myoblasts affected Igf2 expression, cell proliferation, wound healing, and myotube formation. Chromatin immunoprecipitation (ChIP) sequencing using C2C12 cells identified about 2,500 ZBED6 binding sites in the genome, and the deduced consensus motif gave a perfect match with the established binding site in Igf2. Genes associated with ZBED6 binding sites showed a highly significant enrichment for certain Gene Ontology classifications, including development and transcriptional regulation. The phenotypic effects in mutant pigs and ZBED6-silenced C2C12 myoblasts, the extreme sequence conservation, its nucleolar localization, the broad tissue distribution, and the many target genes with essential biological functions suggest that ZBED6 is an important transcription factor in placental mammals, affecting development, cell proliferation, and growth.

The present study was conducted to study the diversity of MHC-DRB3 alleles in Indian cattle and buffalo breeds. Previously reported BoLA-DRB exon 2 alleles of Indian Zebu cattle, Bos taurus cattle, buffalo, sheep, and goats were analyzed for the identities and divergence among various allele sequences. Comparison of predicted amino acid residues of DRB3 exon 2 alleles with similar alleles from other ruminants revealed considerable congruence in amino acid substitution pattern. These alleles showed a high degree of nucleotide and amino acid polymorphism at positions forming peptide-binding regions. A higher rate of nonsynonymous substitution was detected at the peptide-binding regions, indicating that BoLA-DRB3 allelic sequence evolution was driven by positive selection.

Cathelicidins are an ancient class of antimicrobial peptides (AMPs) with broad spectrum bactericidal activities. In this study, we investigated the diversity and biological activity of cathelicidins of buffalo, a species known for its disease resistance. A series of new homologs of cathelicidin4 (CATHL4), which were structurally diverse in their antimicrobial domain, was identified in buffalo. AMPs of newly identified buffalo CATHL4s (buCATHL4s) displayed potent antimicrobial activity against selected Gram positive (G+) and Gram negative (G-) bacteria. These peptides were prompt to disrupt the membrane integrity of bacteria and induced specific changes such as blebing, budding, and pore like structure formation on bacterial membrane. The peptides assumed different secondary structure conformations in aqueous and membrane-mimicking environments. Simulation studies suggested that the amphipathic design of buCATHL4 was crucial for water permeation following membrane disruption. A great diversity, broad-spectrum antimicrobial action, and ability to induce an inflammatory response indicated the pleiotropic role of cathelicidins in innate immunity of buffalo. This study suggests short buffalo cathelicidin peptides with potent bactericidal properties and low cytotoxicity have potential translational applications for the development of novel antibiotics and antimicrobial peptidomimetics.

The NADPH-dependent HC-toxin reductases (HCTR1 and 2) encoded by enzymatic class of disease resistance homologous genes (Hm1 and Hm2) protect maize by detoxifying a cyclic tetrapeptide, HC-toxin, secreted by the fungus Cochliobolus carbonum race 1(CCR1). Unlike the other classes' resistance (R) genes, HCTR-mediated disease resistance is an inimitable mechanism where the avirulence (Avr) component from CCR1 is not involved in toxin degradation. In this study, we attempted to decipher cofactor (NADPH) recognition and mode of HC-toxin binding to HCTRs through molecular docking, molecular dynamics (MD) simulations and binding free energy calculation methods. The rationality and the stability of docked complexes were validated by 30-ns MD simulation. The binding free energy decomposition of enzyme-cofactor complex was calculated to find the driving force behind cofactor recognition. The overall binding free energies of HCTR1-NADPH and HCTR2-NADPH were found to be -616.989 and -16.9749 kJ mol-1 respectively. The binding free energy decomposition revealed that the binding of NADPH to the HCTR1 is mainly governed by van der Waals and nonpolar interactions, whereas electrostatic terms play dominant role in stabilizing the binding mode between HCTR2 and NADPH. Further, docking analysis of HC-toxin with HCTR-NADPH complexes showed a distinct mode of binding and the complexes were stabilized by a strong network of hydrogen bond and hydrophobic interactions. This study is the first in silico attempt to unravel the biophysical and biochemical basis of cofactor recognition in enzymatic class of R genes in cereal crop maize.

The cDNA sequence of leuteinizing hormone (LH) beta subunit was determined to understand the molecular basis for silent oestrus behavior and poor response to superovulation in buffalo. The LH-beta cDNA contains an open reading frame 426 bp long. The deduced sequence of the LH-beta is 141 amino acids in length. The amino acid sequences of the Indian river buffalo LH-beta subunit showed overall similarity to those of other mammals. The nucleotide sequence variability of LH-beta was studied in more than approximately 350 Indian buffaloes covering five different breeds. The results of the sequence analysis showed that the buffalo LH-beta gene is not highly conserved and non-synonymous mutations are not rare, at least in the samples collected randomly from five different breeds and buffalo populations. A total of seven different variants were obtained. In spite of its crucial role in reproduction, variation of the LH-beta gene was found present in this species. The polymorphisms found were unique in the Indian river buffalo population.

The innate immune system is activated upon virus invasion of a host cell by recognizing viral component, such as dsRNA through specific receptors, resulting in the production of type- I IFNs, which confer an antiviral state within the invaded as well as surrounding cells. In the present study, fibroblast, monocyte and macrophage cells derived from water Buffalo (Bubalus bubalis) were exposed to a synthetic dsRNA analogue, poly I:C to mimic viral invasion in each cell type. Recognition of poly I:C through cytosolic helicase receptors RIG-I and MDA5 molecule lead to the activation of the RLR pathway, subsequently activating the MAVS-IRF3/7 cascade and the production of antiviral effector molecule like IFNβ and ISGs. Within the different cell types, we identified variability in RLR receptor and IFNβ expression after poly I:C administration. Fibroblasts responded quickly and strongly with IFNβ production, followed by macrophages and monocytes. Despite absolute expression variability among different cell types the expression trend of RLRs pathway genes were similar. Length of poly I:C molecule also influence IFNβ expression in response of RLR pathway. Short (LMW) poly I:C induce stronger IFN-β expression in myeloid (macrophage and monocyte) cells. In contrast long (HMW) poly I:C preferably elicit higher IFNβ expression in non-myeloid (fibroblast) cell. Therefore, MDA5 and RIG-1 plays an indispensable role in eliciting antiviral response in non- immune (fibroblast) host cell. Thus, stimulation of RLR pathway with suitable and potentially cell-type specific agonist molecules successfully elicit antiviral state in the host animal, with fibroblasts conferring a stronger antiviral state compared with the monocytes and macrophages.

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.

Heat shock protein 70 (HSP70) is an important chaperone, involved in protein folding, refolding, translocation and complex remodeling reactions under normal as well as stress conditions. However, expression of HSPA1A gene in heat and cold stress conditions associates with other chaperons and perform its function. Experimental structure for Camel HSP70 protein (cHSP70) has not been reported so far. Hence, we constructed 3D models of cHSP70 through multi- template comparative modeling with HSP110 protein of S. cerevisiae (open state) and with HSP70 protein of E. coli 70kDa DnaK (close state) and relaxed them for 100 nanoseconds (ns) using all-atom Molecular Dynamics (MD) Simulation. Two stable conformations of cHSP70 with Substrate Binding Domain (SBD) in open and close states were obtained. The collective mode analysis of different transitions of open state to close state and vice versa was examined via Principal Component Analysis (PCA) and Minimum Distance Matrix (MDM). The results provide mechanistic representation of the communication between Nucleotide Binding Domain (NBD) and SBD to identify the role of sub domains in conformational change mechanism, which leads the chaperone cycle of cHSP70. Further, residues present in the chaperon functioning site were also identified through protein-peptide docking. This study provides an overall insight into the inter domain communication mechanism and identification of the chaperon binding cavity, which explains the underlying mechanism involved during heat and cold stress conditions in camel.

RIG1 and MDA5 have emerged as important intracellular innate pattern recognition receptors that recognize viral RNA and mediate cellular signals controlling Type I interferon (IFN-I) response. Buffalo RIG1 and MDA5 genes were investigated to understand the mechanism of receptor induced antiviral response. Sequence analysis revealed that RIG1 and MDA5 maintain a domain arrangement that is common in mammals. Critical binding site residues of the receptors are evolutionary conserved among mammals. Molecular dynamics simulations suggested that RIG1 and MDA5 follow a similar, if not identical, dsRNA binding pattern that has been previously reported in human. Moreover, binding free energy calculation revealed that MDA5 had a greater affinity towards dsRNA compared to RIG1. Constitutive expressions of RLR genes were ubiquitous in different tissues without being specific to immune organs. Poly I:C stimulation induced elevated expressions of IFN-β and IFN-stimulated genes (ISGs) through interferon regulatory factors (IRFs) mediated pathway in buffalo foetal fibroblast cells. The present study provides crucial insights into the structure and function of RIG1 and MDA5 receptors in buffalo.

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.