TLR4 in complex with MD2 senses the presence of lipid A (LA) and initiates a signaling cascade that curb the infection. This complex is evolutionarily conserved and can initiate the immune system in response to a variety of LAs. In this study, molecular dynamics simulation (25 ns) was performed to elucidate the differential behavior of TLR4/MD2 complex in response to Rhodobacter sphaeroides lipid A (RsLA). Penta-acyl chain-containing RsLA is at the verge of agonist (6 acyl-chains) and antagonist (4 acyl-chains) structure, and activates the TLR4 pathway in horses and hamsters, while inhibiting in humans and murine. In the time-evolved coordinates, the promising factors that dictated the differential response included the local and global mobility pattern of complexes, solvent-accessible surface area of ligand, and surface charge distributions of TLR4 and MD2. We showed that the GlcN1-GlcN2 backbone acquires agonist (3FXI)-like configurations in horses and hamsters, while acquiring antagonist (2E59)-like configurations in humans and murine systems. Moreover, analysis of F126 behavior in the MD2 F126 loop (amino acids 123-129) and loop EF (81-89) suggested that certain sequence variations also contribute to species-specific response. This study underlines the TLR4 signaling mechanism and provides new therapeutic opportunities.

Toll-like receptors (TLRs) activate a potent immunostimulatory response. There is clear evidence that overactivation of TLRs leads to infectious and inflammatory diseases. Recent biochemical studies have shown that the membrane-bound form of ST2 (ST2L), a member of the Toll-like/IL-1 receptor superfamily, negatively regulates MyD88-dependent TLR signaling pathways by sequestrating the adapters MyD88 and Mal (TIRAP). Specifically, ST2L attenuates the recruitment of Mal and MyD88 adapters to their receptors through its intracellular TIR domain. Thus, ST2L is a potent molecule that acts as a key regulator of endotoxin tolerance and modulates innate immunity. So far, the inhibitory mechanism of ST2L at the molecular level remains elusive. To develop a working hypothesis for the interactions between ST2L, TLRs (TLR1, 2, 4, and 6), and adapter molecules (MyD88 and Mal), we constructed three-dimensional models of the TIR domains of TLR4, 6, Mal, and ST2L based on homology modeling. Since the crystal structures of the TIR domains of TLR1, 2 as well as the NMR solution structure of MyD88 are known, we utilized these structures in our analysis. The TIR domains of TLR1, 2, 4, 6, MyD88, Mal and ST2L were subjected to molecular dynamics (MD) simulations in an explicit solvent environment. The refined structures obtained from the MD simulations were subsequently used in molecular docking studies to probe for potential sites of interactions. Through protein-protein docking analysis, models of the essential complexes involved in TLR2 and 4 signaling and ST2L inhibiting processes were developed. Our results suggest that ST2L may exert its inhibitory effect by blocking the molecular interface of Mal and MyD88 adapters mainly through its BB-loop region. Our predicted oligomeric signaling models may provide a basis for the understanding of the assembly process of TIR domain interactions, which has thus far proven to be difficult via in vivo studies.

Activating transcription factor 3 (ATF3) is induced by inflammatory responses, cell death, cytokines, and oxidative stress conditions. ATF3 is a negative regulator in the Toll-like receptor 4 signalling pathway. The principal molecule in this pathway is nuclear factor κB (NF-κB) that translocates into the nucleus to initiate the transcription of inflammatory mediators. However, scarce data are available regarding the interaction of ATF3 and p65, a part of the NF-κB dimer. Therefore, we studied the mechanism of regulation of p65 by ATF3 in RAW 264.7 cells. First, LPS-mediated NF-κB activation was confirmed, and then the direct interaction of ATF3 and p65 was observed through immunoprecipitation (IP). The presence of histone deacetylase 1 (HDAC1) was also detected in the complex. In ATF3 deficient cells, NF-κB activity was up-regulated and HDAC1 was not detected by IP. These observations suggest that p65 is attenuated by ATF3 such that ATF3 recruits HDAC1 to the ATF3/p65 complex and facilitates the deacetylation of p65. Likewise, inflammatory response genes were induced by translocated NF-κB in ATF3-deficient cells. Cumulatively, we uncovered a novel mechanism for the negative regulation of NF-κB by ATF3 via direct interaction with p65.

Streptococcus pneumoniae is a polysaccharide-encapsulated bacterium. The capsule thickens during blood invasion compared with the thinner capsules observed in asymptomatic nasopharyngeal colonization. However, the underlying mechanism regulating differential CPS expression remains unclear. CPS synthesis requires energy that is supplied by ATP. Previously, we demonstrated a correlation between ATP levels and adenylate kinase in S. pneumoniae (SpAdK). A dose-dependent induction of SpAdK in serum was also reported. To meet medical needs, this study aimed to elucidate the role of SpAdK in the regulation of CPS production. CPS levels in S. pneumoniae type 2 (D39) increased proportionally with SpAdK levels, but they were not related to pneumococcal autolysis. Moreover, increased SpAdK levels resulted in increased total tyrosine kinase Cps2D levels and phosphorylated Cps2D, which is a regulator of CPS synthesis in the D39 strain. Our results also indicated that the SpAdK and Cps2D proteins interact in the presence of Mg-ATP. In addition, in silico analysis uncovered the mechanism behind this protein-protein interaction, suggesting that SpAdK binds with the Cps2D dimer. This established the importance of the ATP-binding domain of Cps2D. Taken together, the biophysical interaction between SpAdK and Cps2D plays an important role in enhancing Cps2D phosphorylation, which results in increased CPS synthesis.

A variety of ligands differ in their capacity to bind the receptor, elicit gene expression, and modulate physiological responses. Such receptors include Toll-like receptors (TLRs), which recognize various patterns of pathogens and lead to primary innate immune activation against invaders, and G-protein coupled receptors (GPCRs), whose interaction with their cognate ligands activates heterotrimeric G proteins and regulates specific downstream effectors, including immuno-stimulating molecules. Once TLRs are activated, they lead to the expression of hundreds of genes together and bridge the arm of innate and adaptive immune responses. We characterized the gene expression profile of Toll-like receptor 4 (TLR4) in RAW 264.7 cells when it bound with its ligand, 2-keto-3-deoxyoctonate (KDO), the active part of lipopolysaccharide. In addition, to determine the network communications among the TLR, Janus kinase (JAK)/signal transducer and activator of transcription (STAT), and GPCR, we tested RAW 264.7 cells with KDO, interferon-β, or cAMP analog 8-Br. The ligands were also administered as a pair of double and triple combinations.

The interleukin-1 receptor-associated kinase (IRAK) family comprises critical signaling mediators of the TLR/IL-1R signaling pathways. IRAKs are Ser/Thr kinases. There are 4 members in the vertebrate genome (IRAK1, IRAK2, IRAKM, and IRAK4) and an IRAK homolog, Pelle, in insects. IRAK family members are highly conserved in vertebrates, but the evolutionary relationship between IRAKs in vertebrates and insects is not clear. To investigate the evolutionary history and functional divergence of IRAK members, we performed extensive bioinformatics analysis. The phylogenetic relationship between IRAK sequences suggests that gene duplication events occurred in the evolutionary lineage, leading to early vertebrates. A comparative phylogenetic analysis with insect homologs of IRAKs suggests that the Tube protein is a homolog of IRAK4, unlike the anticipated protein, Pelle. Furthermore, the analysis supports that an IRAK4-like kinase is an ancestral protein in the metazoan lineage of the IRAK family. Through functional analysis, several potentially diverged sites were identified in the common death domain and kinase domain. These sites have been constrained during evolution by strong purifying selection, suggesting their functional importance within IRAKs. In summary, our study highlighted the molecular evolution of the IRAK family, predicted the amino acids that contributed to functional divergence, and identified structural variations among the IRAK paralogs that may provide a starting point for further experimental investigations.

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.

Ginseng (Panax ginseng Meyer) has been shown to have anti-aging effects in animal and clinical studies. However, the molecular mechanisms by which ginseng exerts these effects remain unknown. Here, the anti-aging effect of Korean red ginseng (KRG) in rat testes was examined by system biology analysis. KRG water extract prepared in feed pellets was administered orally into 12 month old rats for 4 months, and gene expression in testes was determined by microarray analysis. Microarray analysis identified 33 genes that significantly changed. Compared to the 2 month old young rats, 13 genes (Rps9, Cyp11a1, RT1-A2, LOC365778, Sv2b, RGD1565959, RGD1304748, etc.) were up-regulated and 20 genes (RT1-Db1, Cldn5, Svs5, Degs1, Vdac3, Hbb, LOC684355, Svs5, Tmem97, Orai1, Insl3, LOC497959, etc.) were down-regulated by KRG in the older rats. Ingenuity Pathway Analysis of untreated aged rats versus aged rats treated with KRG showed that the affected most was Cyp11a1, responsible for C21-steroid hormone metabolism, and the top molecular and cellular functions are organ morphology and reproductive system development and function. When genes in young rat were compared with those in the aged rat, sperm capacitation related genes were down-regulated in the old rat. However, when genes in the old rat were compared with those in the old rat treated with KRG, KRG treatment up-regulated C21-steroid hormone metabolism. Taken together, Cyp11a1 expression is decreased in the aged rat, however, it is up-regulated by KRG suggesting that KRG seems enhance testes function via Cyp11a1.

Toll-like receptors (TLRs) are pattern recognition receptors that recognize pathogens based on distinct molecular signatures. The human (h)TLR1, 2, 6 and 10 belong to the hTLR1 subfamilies, which are localized in the extracellular regions and activated in response to diverse ligand molecules. Due to the unavailability of the hTLR10 crystal structure, the understanding of its homo and heterodimerization with hTLR2 and hTLR1 and the ligand responsible for its activation is limited. To improve our understanding of the TLR10 receptor-ligand interaction, we used homology modeling to construct a three dimensional (3D) structure of hTLR10 and refined the model through molecular dynamics (MD) simulations. We utilized the optimized structures for the molecular docking in order to identify the potential site of interactions between the homo and heterodimer (hTLR10/2 and hTLR10/1). The docked complexes were then used for interaction with ligands (Pam(3)CSK(4) and PamCysPamSK(4)) using MOE-Dock and ASEDock. Our docking studies have shown the binding orientations of hTLR10 heterodimer to be similar with other TLR2 family members. However, the binding orientation of hTLR10 homodimer is different from the heterodimer due to the presence of negative charged surfaces at the LRR11-14, thereby providing a specific cavity for ligand binding. Moreover, the multiple protein-ligand docking approach revealed that Pam(3)CSK(4) might be the ligand for the hTLR10/2 complex and PamCysPamSK(4,) a di-acylated peptide, might activate hTLR10/1 hetero and hTLR10 homodimer. Therefore, the current modeled complexes can be a useful tool for further experimental studies on TLR biology.

Activation of naïve B lymphocytes by extracellular ligands, e.g. antigen, lipopolysaccharide (LPS) and CD40 ligand, induces a combination of common and ligand-specific phenotypic changes through complex signal transduction pathways. For example, although all three of these ligands induce proliferation, only stimulation through the B cell antigen receptor (BCR) induces apoptosis in resting splenic B cells. In order to define the common and unique biological responses to ligand stimulation, we compared the gene expression changes induced in normal primary B cells by a panel of ligands using cDNA microarrays and a statistical approach, CLASSIFI (Cluster Assignment for Biological Inference), which identifies significant co-clustering of genes with similar Gene Ontology annotation.

Large-scale gene expression studies have mainly focused on highly expressed and 'discriminatory' genes to decipher key regulatory processes. Biological responses are consequence of the concerted action of gene regulatory network, thus, limiting our attention to genes having the most significant variations is insufficient for a thorough understanding of emergent whole genome response. Here we comprehensively analyzed the temporal oligonucleotide microarray data of lipopolysaccharide (LPS) stimulated macrophages in 4 genotypes; wildtype, Myeloid Differentiation factor 88 (MyD88) knockout (KO), TIR-domain-containing adapter-inducing interferon-beta (TRIF) KO and MyD88/TRIF double KO (DKO). Pearson correlations computed on the whole genome expression between different genotypes are extremely high (>0.98), indicating a strong co-regulation of the entire expression network. Further correlation analyses reveal genome-wide response is biphasic, i) acute-stochastic mode consisting of small number of sharply induced immune-related genes and ii) collective mode consisting of majority of weakly induced genes of diverse cellular processes which collectively adjust their expression level. Notably, temporal correlations of a small number of randomly selected genes from collective mode show scalability. Furthermore, in collective mode, the transition from large scatter in expression distributions for single ORFs to smooth linear lines emerges as an organizing principle when grouping of 50 ORFs and above. With this emergent behavior, the role of MyD88, TRIF and novel MyD88, TRIF-independent processes for gene induction can be linearly superposed to decipher quantitative whole genome differential control of transcriptional and mRNA decay machineries. Our work demonstrates genome-wide co-regulated responses subsequent to specific innate immune stimulus which have been largely neglected.

Negative regulation of Toll-like receptor-4 (TLR4) is anticipated to control the pathogen-induced exaggerated immune response. However, effective TLR4 antagonists with scarce off-target effects are yet to be developed. To fill this void, we sought to design small peptide-inhibitors of the TLR4/MD2-LPS interaction. Here we report novel TLR4-antagonistic peptides (TAP), identified through phage display, endowed with the LPS-induced proinflammation inhibition, and confirmed in mice. TAPs-attributed TLR4-antagonism were initially evaluated through NF-κB inhibition in HEK-blue hTLR4 and RAW264.7 cells, and further reinforced by the downregulation of MAPKs (mitogen-activated protein kinases), NF-κB, interleukin 6, and suppression of the oxidative-stress products and iNOS in macrophages and human peripheral blood mononuclear cells (hPBMCs). Among these, TAP2 specifically halted the TLR4, but not other TLRs signaling, which was further confirmed by the biophysical kinetic assay. Finally, TAP2 diminished LPS-elicited systemic cytokine response in vivo, suggesting that TAPs, specifically TAP2, have the potential to treat TLR4-mediated immune ailments.

Despite being crucial for combating microbes, paradoxical Toll-like receptors (TLRs) signaling have been associated with the aggravation of multiple immune disorders such as systemic lupus erythematosus, psoriasis, rheumatoid arthritis, and nonalcoholic steatohepatitis. The stoichiometry and precise arrangement of the interaction of adapters (via their Toll/interleukin-1 receptor [TIR] domains) are indispensable for the activation of TLRs and of downstream signaling cascades. Among adapters, plasma membrane-anchored MyD88 adaptor-like (MAL) has the potential for BB-loop-mediated self-oligomerization and interacts with other TIR domain-containing adaptors through αC and αD helices. Here, we used information on the MAL-αC interface to exploit its pharmacophores and to design a decoy peptide (MIP2) with broad-range TLR-inhibitory abilities. MIP2 abrogated MyD88- and TRIF-dependent lipopolysaccharide (LPS)-induced TLR4 signaling in murine and human cell lines and manifested a therapeutic potential in models of psoriasis, systemic lupus erythematosus, nonalcoholic steatohepatitis, and sepsis. Levels of hallmark serological and histological biomarkers were significantly restored and the disease symptoms were substantially ameliorated by MIP2 treatment of the animals. Collectively, our biophysical, in vitro, and in vivo findings suggest that MIP2 has broad specificity for TLRs and may be effective in modulating autoimmune complications caused by microbial or environmental factors.

Mutations affect gene functions related to cancer behavior, including cell growth, metastasis, and drug responses. Genome-wide profiling of cancer mutations and drug responses has identified actionable targets that can be utilized for the management of cancer patients. Here, the recapitulation of pharmacogenomic data revealed that the mutation of EPHB6 is associated with paclitaxel resistance in cancer cells. Experimental data confirmed that the EPHB6 mutation induces paclitaxel resistance in various cancer types, including lung, skin, and liver cancers. EPHB6 mutation-induced paclitaxel resistance was mediated by an interaction with EPHA2, which promotes c-Jun N-terminal kinase (JNK)-mediated cadherin 11 (CDH11) expression. We demonstrated that EPHB6-mutated cells acquire cell adhesion-mediated drug resistance (CAM-DR) in association with CDH11 expression and RhoA/focal adhesion kinase (FAK) activation. Targeted inhibition of EPHA2 or CDH11 reversed the acquired paclitaxel resistance, suggesting its potential clinical utility. The present results suggest that the EPHB6 mutation and its downstream EPHA2/JNK/CDH11/RhoA/FAK signaling axis are novel diagnostic and therapeutic targets for overcoming paclitaxel resistance in cancer patients.

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 novel coronavirus disease, caused by severe acute respiratory coronavirus 2 (SARS-CoV-2), rapidly spreading around the world, poses a major threat to the global public health. Herein, we demonstrated the binding mechanism of PF-07321332, α-ketoamide, lopinavir, and ritonavir to the coronavirus 3-chymotrypsin-like-protease (3CLpro) by means of docking and molecular dynamic (MD) simulations. The analysis of MD trajectories of 3CLpro with PF-07321332, α-ketoamide, lopinavir, and ritonavir revealed that 3CLpro-PF-07321332 and 3CLpro-α-ketoamide complexes remained stable compared with 3CLpro-ritonavir and 3CLpro-lopinavir. Investigating the dynamic behavior of ligand-protein interaction, ligands PF-07321332 and α-ketoamide showed stronger bonding via making interactions with catalytic dyad residues His41-Cys145 of 3CLpro. Lopinavir and ritonavir were unable to disrupt the catalytic dyad, as illustrated by increased bond length during the MD simulation. To decipher the ligand binding mode and affinity, ligand interactions with SARS-CoV-2 proteases and binding energy were calculated. The binding energy of the bespoke antiviral PF-07321332 clinical candidate was two times higher than that of α-ketoamide and three times than that of lopinavir and ritonavir. Our study elucidated in detail the binding mechanism of the potent PF-07321332 to 3CLpro along with the low potency of lopinavir and ritonavir due to weak binding affinity demonstrated by the binding energy data. This study will be helpful for the development and optimization of more specific compounds to combat coronavirus disease.

Systemic lupus erythematosus (SLE) is a common chronic autoimmune inflammatory disease According to recent studies, signaling through Toll-like receptor (TLR) protein, which promotes the production of inflammatory cytokines, leads to the development of SLE TLR-inhibitory peptide 1 (TIP1) has been newly identified for the treatment of autoimmune diseases.

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.

The potential toxicity of nanoparticles, particularly to neurons, is a major concern. In this study, we assessed the cytotoxicity of silica-coated magnetic nanoparticles containing rhodamine B isothiocyanate dye (MNPs@SiO2(RITC)) in HEK293 cells, SH-SY5Y cells, and rat primary cortical and dopaminergic neurons. In cells treated with 1.0 μg/μl MNPs@SiO2(RITC), the expression of several genes related to the proteasome pathway was altered, and proteasome activity was significantly reduced, compared with control and with 0.1 μg/μl MNPs@SiO2(RITC)-treated cells. Due to the reduction of proteasome activity, formation of cytoplasmic inclusions increased significantly in HEK293 cells over-expressing the α-synuclein interacting protein synphilin-1 as well as in primary cortical and dopaminergic neurons. Primary neurons, particularly dopaminergic neurons, were more vulnerable to MNPs@SiO2(RITC) than SH-SY5Y cells. Cellular polyamines, which are associated with protein aggregation, were significantly altered in SH-SY5Y cells treated with MNPs@SiO2(RITC). These findings highlight the mechanisms of neurotoxicity incurred by nanoparticles.