Jitendra Maharana

Structural insights of rohu TLR3, its binding site analysis with fish reovirus dsRNA, poly I:C and zebrafish TRIF

In response to double stranded RNA (dsRNA) viruses, toll-like receptor 3 (TLR3) in fish activates signaling like human, and induces innate immunity. This suggested the existence of dsRNA binding domains in fish TLR3 as reported in higher vertebrates. In in silico analysis, leucine rich repeat (LRR) regions (4-6, 13-14, 20-22), and LRR (8-15, 17-24) were identified as key domains in rohu TLR3 as poly I:C and dsRNA of fish reovirus (AGCRV,VHSV and IHNV) binding regions. 3D-models of rohu TLR3-TIR and zebrafish TRIF were generated by homology and ab initio modeling respectively, and their interacting domains were predicted. This is the first report of TLR3 modeling in fish.

Structural insights into the MDP binding and CARD–CARD interaction in zebrafish (Danio rerio) NOD2: a molecular dynamics approach

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. Copyright

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

In lower eukaryotes-like fish, innate immunity contributed by various pattern recognition receptor (PRR) plays an essential role in protection against diseases. Nucleotide-binding and oligomerization domain (NOD)-2 is a cytoplasmic PRR that recognizes MDP (muramyl dipeptide) of the Gram positive and Gram negative bacteria as ligand and activates signalling to induce innate immunity. Hypothesizing a similar NOD2 signalling pathway of higher eukaryotes, the peripheral blood leucocytes (PBLs) of rohu (Labeo rohita) was stimulated with MDP. The data of quantitative real-time PCR (qRT-PCR) revealed MDP-mediated inductive expression of NOD2 and its down-stream molecule RICK/RIP2 (receptor-interacting serine-threonine protein kinase-2). This observation suggested the existence of MDP-binding sites in rohu NOD2 (rNOD2). To investigate it, 3D model of ligand-binding leucine-rich repeat (LRR) region of rNOD2 (rNOD2-LRR) was constructed following ab initio and threading approaches in I-TASSER web server. Structural refinement of the model was performed by energy minimization, and MD (molecular dynamics) simulation was performed in GROMACS (Groningen Machine for Chemical Simulations). The refined model of rNOD2-LRR was validated through SAVES, ProSA, ProQ, WHAT IF and MolProbity servers, and molecular docking with MDP was carried out in GOLD 4.1. The result of docking identified LRR3-7 comprising Lys820, Phe821, Asn822, Arg847, Gly849, Trp877, Trp901 and Trp931 as MDP-binding critical amino acids in rNOD2. This is the first study in fish to provide an insight into the 3D structure of NOD2-LRR region and its important motifs that are expected to be engaged in MDP binding and innate immunity.

Understanding the distinguishable structural and functional features in zebrafish TLR3 and TLR22, and their binding modes with fish dsRNA viruses: an exploratory structural model analysis

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.

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.

Structural and dynamic investigation of bovine folate receptor alpha (FOLR1), and role of ultra-high temperature processing on conformational and thermodynamic characteristics of FOLR1–folate complex

The folate receptor alpha (FOLR1) present in milk has widely been studied to investigate the effects of pasteurization, ultra-high temperature (UHT) processing and fermentation on net folate concentration. However, the folate binding mechanism with FOLR1, and effect of temperature on FOLR1-folate complex is poorly explored till now in bovine milk which is a chief resource of folate. Despite of enormous importance of folic acid and the routine intake of bovine milk, folic acid deficiency diseases are common in human race. To understand the folate deficiency in milk after processing, in absence of experimental structure, 3D model of bovine FOLR1 (bvFOLR1) was built followed by 40ns molecular dynamics (MD) simulation. The folate and its derivatives binding sites in bvFOLR1 were anticipated by molecular docking using AutoDock 4.2. Essential MD studies suggested the presence of a longer signal peptide (22 residues) and a short propeptide (7 residues) at the C-terminus that may cleaved during post-translational modification. MD analysis of bvFOLR1-folate complex at 298, 323, 353, 373 and 408K followed by binding energy (BE) calculation showed maximum binding affinity at ∼353K. However, at 373K and UHT (408K), the folate BE is significantly decreased with substantial conformational alteration. Heating at UHT followed by cooling within 298-408K range demoed no structural reformation with temperature reduction, and the folate was displaced from the active site. This study presented the disintegration of folate from bvFOLR1 during high temperature processing and revealed a lower folate concentration in UHT milk and dairy products.

Molecular recognition of avirulence protein (avrxa5) by eukaryotic transcription factor xa5 of rice (Oryza sativa L.): Insights from molecular dynamics simulations

The avirulence gene avrxa5 of bacterial blight pathogen Xanthomonas oryzae pv. oryzae (Xoo) recognized by the resistant rice lines having corresponding resistance (xa5) gene in a gene-for-gene manner. We used a combinatorial approach involving protein-protein docking, molecular dynamics (MD) simulations and binding free energy calculations to gain novel insights into the gene-for-gene mechanism that governs the direct interaction of R-Avr protein. From the best three binding poses predicted by molecular docking, MD simulations were performed to explore the dynamic binding mechanism of xa5 and avrxa5. Molecular Mechanics/Poisson Boltzmann Surface Area (MM/PBSA) techniques were employed to calculate the binding free energy and to uncover the thriving force behind the molecular recognition of avrxa5 by eukaryotic transcription factor xa5. Binding free energy analysis revealed van der Waals term as the most constructive component that favors the xa5 and avrxa5 interaction. In addition, hydrogen bonds (H-bonds) and essential electrostatic interactions analysis highlighted amino acid residues Lys54/Asp870, Lys56/Ala868, Lys56/Ala866, Lys56/Glu871, Ile59/His862, Gly61/Phe858, His62/Arg841, His62/Leu856, Ser101/Ala872 and Ser105/Asp870 plays pivotal role for the energetically stability of the R-Avr complex. Insights gained from the present study are expected to unveil the molecular mechanisms that define the transcriptional activator mediated transcriptome modification in host plants.

Molecular dynamics simulation of neuropeptide B and neuropeptide W in the dipalmitoylphosphatidylcholine membrane bilayer

GPR7 and GPR8 are recently deorphanized G-protein-coupled receptors that are implicated in the regulation of neuroendocrine function, feeding behavior, and energy homeostasis. Neuropeptide B (NPB) and neuropeptide W (NPW) are two membrane-bound hypothalamic peptides, which specifically antagonize GPR7 and GPR8. Despite years of research, an accurate estimation of structure and molecular recognition of these neuropeptide systems still remains elusive. Herein, we investigated the structure, orientation, and interaction of NPB and NPW in a dipalmitoylphosphatidylcholine bilayer using long-range molecular dynamics (MD) simulation. During 30-ns simulation, membrane-embedded helical axes of NPB and NPW tilted 30 and 15°, respectively, from the membrane normal in order to overcome possible hydrophobic mismatch with the lipid bilayer. The calculation of various structural parameters indicated that NPW is more rigid and compact as compared to NPB. Qualitatively, the peptides exhibited flexible N-terminal (residues 1–12) and rigid C-terminal α-helical parts (residues 13–21), confirming previous NMR data. A strong electrostatic attraction between C-termini and headgroup atoms caused translocation of the peptides towards lower leaflet of the bilayer. The stabilizing hydrogen bonds (H-bonds) between phosphate groups and Trp1, Lys3, and Arg15 of the peptides played important roles for membrane anchoring. MD simulations of Alanine (Ala) mutants revealed that WYK->Ala variant of NPB/NPW lacked crucial H-bond interactions with phospholipid headgroups and also caused severe misfolding in NPB. Altogether, the knowledge of preferred structural fold and interaction of neuropeptides within the membrane bilayer will be useful to develop synthetic agonist or antagonist peptides for GPR7 and GPR8.

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.

NOD1CARD Might Be Using Multiple Interfaces for RIP2-Mediated CARD-CARD Interaction: Insights from Molecular Dynamics Simulation

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.