Mahesh Chandra Patra

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

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.

Diversity, Antimicrobial Action and Structure-Activity Relationship of Buffalo Cathelicidins.

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.

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

Comparative Genomic Analysis of Buffalo (Bubalus bubalis) NOD1 and NOD2 Receptors and Their Functional Role in In-Vitro Cellular Immune Response

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.

Insight into buffalo (Bubalus bubalis) RIG1 and MDA5 receptors: a comparative study on dsRNA recognition and in-vitro antiviral response.

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.

Comparative modeling of human kappa opioid receptor and docking analysis with the peptide YFa

The kappa opioid receptor belongs to the super family of G protein - coupled receptors that are of utmost significance in the development of potent analgesic drugs for the treatment of severe pain. An accurate evaluation of the ligand binding pathways into this receptor at molecular level may play a key role in the design of new molecules with more desirable properties and reduced side effects. In this study, homology model of the human kappa opioid receptor was developed by MODELLER using the X-ray crystal structure of bovine rhodopsin as template. Initial structure of the receptor was refined computationally with energy minimization and molecular dynamics simulation at 300 K in a pre-equilibrated phospholipid bilayer by GROMACS. The Met-enkaphalin-Arg-Phe based opioid peptide YFa (YGGFMKKKFMRF) designed and characterized by our laboratory was docked into the optimized model and the critical amino acids responsible for binding were identified. A number of low energy binding poses of YFa with the receptor were assessed after the molecular docking in which the peptide was observed to interact with the receptors extracellular amino acids through hydrogen bonds. The human kappa opioid receptor model optimized in a phospholipid bilayer should provide a good starting point for further characterization of the binding modes of other opioid ligands. Furthermore, the biologically favorable molecular interactions between YFa and human kappa opioid receptor observed by our study might be able to justify the specificity of this peptide.

Structural Analysis of Respirasomes in Electron Transfer Pathway of Acidithiobacillus ferrooxidans: A Computer-Aided Molecular Designing Study

Acidithiobacillus ferrooxidans obtains its metabolic energy by reducing extracellular ferrous iron with either downhill or uphill electron transfer pathway. The downhill electron transfer pathway has been substantially explored in recent years to underpin the mechanism of iron respiration but, there exists a wide gap in our present understanding on how these proteins are organized as a supercomplex and what sort of atomic level interactions governs their stability in the iron respiratory chain. In the present study, we aimed at unraveling the structural basis of supermolecular association of respirasomes using protein threading, protein-protein docking, and molecular dynamics (MD) simulation protocols. Our results revealed that Phe312 of outer membrane cytochrome c plays a crucial role in diffusing electrons from heme C group to Asp73 of rusticyanin. In line with the previous experimental results, His143 of rusticyanin was found to have a stable interaction with Glu121 of periplasmic cytochrome c4. Cytochrome c4 interacts with subunit B of cytochrome c oxidase through Lys146 and Thr148 of the conserved hydrophobic/aromatic motif 145-WKWTFSY-151 to attain stability during simulation. Phe468 of cytochrome c oxidase was found indispensable for stabilizing heme aa3 during MD simulation. Taken together, we conclude that the molecular interactions of charged and hydrophobic amino acids present on the surface of each respirasome form a hypothetical electron wire in the iron respiratory supercomplex of A. ferrooxidans.

Comparison of Copy Number of HSF Genes in Two Buffalo Genomes.

ABSTRACT The copy number variation (CNV) is the number of copies of a particular gene in the genotype of an individual. Recent evidences show that the CNVs can vary in frequency and occurrence between breeds. These variations reportedly allowed different breeds to adapt to different environments. As copy number variations follow Mendelian pattern of inheritance, identification and distribution of these variants between populations can be used to infer the evolutionary history of the species. In this study, we have examined the absolute copy number of four Heat shock factor genes viz. HSF-1, 2, 4, and 5 in two different breeds of buffalo species using real-time PCR. Here, we report that the absolute copy number of HSF2 varies between the two breeds. In contrast no significant difference was observed in the copy number for HSF-1, 4, and 5 between the two breeds. Our results provide evidence for the presence of breed specific differences in HSF2 genomic copy number. This seems to be the first step in delineating the genetic factors underlying environmental adaptation between the two breeds. Nevertheless, a more detailed study is needed to characterize the functional consequence of this variation.