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 and functional investigation of zebrafish (Danio rerio) NOD1 leucine rich repeat domain and its interaction with iE-DAP.

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.

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.

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.