A review on structure-function mechanism and signaling pathway of serine/threonine protein PIM kinases as a therapeutic target.

The proviral integration for the Moloney murine leukemia virus (PIM) kinases, belonging to serine/threonine kinase family, have been found to be overexpressed in various types of cancers, such as prostate, breast, colon, endometrial, gastric, and pancreatic cancer. The three isoforms PIM kinases i.e., PIM1, PIM2, and PIM3 share a high degree of sequence and structural similarity and phosphorylate substrates controlling tumorigenic phenotypes like proliferation and cell survival. Targeting short-lived PIM kinases presents an intriguing strategy as in vivo knock-down studies result in non-lethal phenotypes, indicating that clinical inhibition of PIM might have fewer adverse effects. The ATP binding site (hinge region) possesses distinctive attributes, which led to the development of novel small molecule scaffolds that target either one or all three PIM isoforms. Machine learning and structure-based approaches have been at the forefront of developing novel and effective chemical therapeutics against PIM in preclinical and clinical settings, and none have yet received approval for cancer treatment. The stability of PIM isoforms is maintained by PIM kinase activity, which leads to resistance against PIM inhibitors and chemotherapy; thus, to overcome such effects, PIM proteolysis targeting chimeras (PROTACs) are now being developed that specifically degrade PIM proteins. In this review, we recapitulate an overview of the oncogenic functions of PIM kinases, their structure, function, and crucial signaling network in different types of cancer, and the potential of pharmacological small-molecule inhibitors. Further, our comprehensive review also provides valuable insights for developing novel antitumor drugs that specifically target PIM kinases in the future. In conclusion, we provide insights into the benefits of degrading PIM kinases as opposed to blocking their catalytic activity to address the oncogenic potential of PIM kinases.

A novel heterozygote allele in caprine melanocortin 1 receptor (MC1R) gene: an association with heat stress traits.

Climate change negatively influences the productive and reproductive abilities of goats. There is a need to understand the relationship between heat stress and genes that may aid in the development of climate-resilient goats. Melanism variation in goats plays a role in thermoregulation, in which the melanogenic genes have a pleiotropic effect on the regulation of physiological responses and behavior that are altered due to heat stress in the animals. Thus, the present study was conducted to establish a possible association between the coat color gene (MC1R) and heat stress characteristics. The physiological responses and cortisol levels were recorded in forty different coat-colored goats. The genotyping of the animals revealed four SNPs at the 183rd (C/T), 332nd (C/G), 748th (G/T), and 801st (C/G) positions, among which the black and brown goat populations had novel SNPs at the 332nd position. Eight haplotypes were constructed, and an association study revealed that haplotypes (CCGG, TCGG, and CCTC) that were linked to white animals had lower cortisol values, rectal temperature, skin temperature, and respiration rate. The multivariate and cluster analyses revealed that the white goats were distinct from the rest of the goats. In addition, the docking results revealed the residues that were forming the interaction complex, which could play a role in melanogenesis in the animals and, in turn, the heat stress ability of the goats. Altogether, the results of the present study could pave the way for more research into coat color genes and their relationship with heat stress traits.

GIpred: a computational tool for prediction of GIGANTEA proteins using machine learning algorithm.

In plants, GIGANTEA (GI) protein plays different biological functions including carbon and sucrose metabolism, cell wall deposition, transpiration and hypocotyl elongation. This suggests that GI is an important class of proteins. So far, the resource-intensive experimental methods have been mostly utilized for identification of GI proteins. Thus, we made an attempt in this study to develop a computational model for fast and accurate prediction of GI proteins. Ten different supervised learning algorithms i.e., SVM, RF, JRIP, J48, LMT, IBK, NB, PART, BAGG and LGB were employed for prediction, where the amino acid composition (AAC), FASGAI features and physico-chemical (PHYC) properties were used as numerical inputs for the learning algorithms. Higher accuracies i.e., 96.75% of AUC-ROC and 86.7% of AUC-PR were observed for SVM coupled with AAC + PHYC feature combination, while evaluated with five-fold cross validation. With leave-one-out cross validation, 97.29% of AUC-ROC and 87.89% of AUC-PR were respectively achieved. While the performance of the model was evaluated with an independent dataset of 18 GI sequences, 17 were observed as correctly predicted. We have also performed proteome-wide identification of GI proteins in wheat, followed by functional annotation using Gene Ontology terms. A prediction server "GIpred" is freely accessible at http://cabgrid.res.in:8080/gipred/ for proteome-wide recognition of GI proteins. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s12298-022-01130-6.

Identification and futuristic role of novel polymorphism of caprine PrP gene.

The Caprine Prion Protein (PrP) gene polymorphism in three different native Indian goat populations of Southern Odisha, namely Ganjam (a registered breed of India), Ghumusari and Raighar was studied. The 876 bp amplified segment of PrP gene contains full length coding sequence of 771 bp. In Ganjam and Ghumusari goats, any difference of nucleotide sequence was not identified. However, the comparison of nucleotide sequences of Raighar goats and goats of other locality revealed a change in nucleotide at five different positions (G190A, G724A, A727T, C775G and C800T) which includes two non-synonymous nucleotide changes. The non-synonymous nucleotide change resulted a change in amino acid at two different positions (Ser234Cys and Lys246Phe) in mature polypeptide which were not reported earlier and therefore, considered as novel. On the basis of these variants of PrP gene phylogenetic tree was constructed which showed that Ganjam and Raighar goats appeared in different clade. Since any occurrence of Scrapie infection in goats of Odisha was not reported, it can be proposed that these changes in amino acid may be responsible as resistance allele.

PredCRG: A computational method for recognition of plant circadian genes by employing support vector machine with Laplace kernel.

BACKGROUND: Circadian rhythms regulate several physiological and developmental processes of plants. Hence, the identification of genes with the underlying circadian rhythmic features is pivotal. Though computational methods have been developed for the identification of circadian genes, all these methods are based on gene expression datasets. In other words, we failed to search any sequence-based model, and that motivated us to deploy the present computational method to identify the proteins encoded by the circadian genes. RESULTS: Support vector machine (SVM) with seven kernels, i.e., linear, polynomial, radial, sigmoid, hyperbolic, Bessel and Laplace was utilized for prediction by employing compositional, transitional and physico-chemical features. Higher accuracy of 62.48% was achieved with the Laplace kernel, following the fivefold cross- validation approach. The developed model further secured 62.96% accuracy with an independent dataset. The SVM also outperformed other state-of-art machine learning algorithms, i.e., Random Forest, Bagging, AdaBoost, XGBoost and LASSO. We also performed proteome-wide identification of circadian proteins in two cereal crops namely, Oryza sativa and Sorghum bicolor, followed by the functional annotation of the predicted circadian proteins with Gene Ontology (GO) terms. CONCLUSIONS: To the best of our knowledge, this is the first computational method to identify the circadian genes with the sequence data. Based on the proposed method, we have developed an R-package PredCRG ( https://cran.r-project.org/web/packages/PredCRG/index.html ) for the scientific community for proteome-wide identification of circadian genes. The present study supplements the existing computational methods as well as wet-lab experiments for the recognition of circadian genes.

An in vivo and in silico analysis of novel variation in TMBIM6 gene affecting cardiopulmonary traits of Indian goats.

Transmembrane Bax Inhibitor Motif-containing 6 (TMBIM6) gene acts as calcium leak channel and negatively regulates autophagy and autophagosome formation. The TMBIM6 gene was amplified and searched for variation in three different goat populations (i.e. Black Bengal, Ganjam and Raighar) of Odisha state of the India. The result indicated two substitutions i.e. 55th position (C55T) and 95th position (C95A) in the amplified region of the gene resulting in change of amino acids (Leu > Phe and Thr > Asn). The identified SNPs were combined to form haplotypes and animals were grouped accordingly. Structural analysis showed minor changes (5%) in between mutant and wild TMBIM6 protein structures. However, any functional variation could not be identified with respect to the calcium ligand and open pore state. But an alteration of calcium binding site was found. The binding interaction of calcium with the TMBIM6 protein was hydrophobic in nature in closed state whereas hydrophilic in open pore stage. The stress releasing function was the result of calcium leakage controlled by amino acids coded by exon 4 and exon 5 regions of TMBIM6 gene. The effect of breed and haplotype on cardiopulmonary traits was studied. The data on cardiopulmonary traits of body i.e. rectal temperature, skin temperature, heart rate and respiration rate were recorded when ambient temperature usually remained the highest. The statistical analysis showed, significant difference in rectal temperature, skin temperature and respiration rate among these goat populations. The haplotypes (CC and TA) were found to have a significant (P < 0.05) effect on rectal temperature, skin temperature and respiration rate. However, any such significant effect could not be identified in recorded heart rate. The objective of the present study to identify the genetic variations in TMBIM6 gene having significant effect on cardiopulmonary traits which can be further uses as the molecular markers to improve heat tolerance mechanism in goats.

Arabinosyltransferase C enzyme of Mycobacterium tuberculosis, a potential drug target: An insight from molecular docking study.

Multi-drug resistant in Mycobacterium tuberculosis (M.tb) is considered as major bottleneck in the treatment and cure of tuberculosis (TB). Several anti-tubercular drugs fail in its efficacy due to drug-resistant M.tb developed mechanism for resistance. So, research across globe has been carried out to develop effective anti-TB drugs to improve the treatment of these strains. Traditional drug development methods have been proved unsuccessful as it fails to develop a broad-spectrum drug due to lack of structure based approach. Several studies have been conducted in this regard and identified several drug target sites that influence drug-resistant M.tb strains. In this study, the attempt was to study the interaction between the protein Arabinosyltransferase C with the two existing drugs (Ethambutol and Isoniazid) and five modified molecules derived from Ethambutol by calculating their binding affinity and mode of binding through molecular docking study using AutoDock 4. From the comparison study of the existing drug (EMB and INH) and the five proposed modified molecules (Emb1, Emb2, Emb3, Emb4 and Emb5), it is analysed that Emb1 and Emb3 with binding affinities -5.77 kcal/mol and -5.13 kcal/mol respectively can be considered as potential inhibitors of Arabinosyltransferase C in Mycobacterium tuberculosis which is responsible for cell wall synthesis. The facts provided may be further verified experimentally for future drug discovery process to make a stand against tuberculosis and contribute an advance research for worthy antimycobacterial strategies.

Functional diversity and metabolic profile of microbial community of mine soils with different levels of chromium contamination.

Microbial communities provide useful information about any chemical and physical changes in the environment and play an essential role in maintaining soil fertility. Biolog® eco-plates method was used to study the functional diversity of microbial communities, and their correlation with soil organic carbon (OC), microbial biomass and activities, under three different soil conditions of Sukinda chromite mining area of Odisha, India during August 2016. The OC, available nitrogen, available phosphorus and available potash were significantly (p < 0.05) lower in in situ and overburden soils as compared to forest soil. The average development rate of average well color development values decreased with incubation time in all soil conditions. The utilization of six categories of carbon sources by soil microbes decreased with the increase in chromium load and biplot analysis suggested that carbohydrate, polymer and amino acid utilizing microbes were dominant in mining soils. The ecotoxicological status of chromite mine soil would be useful for formulating strategies of possible bioremediation program.

Production, characterization and application of thermostable, alkaline α-amylase (AA11) from Bacillus cereus strain SP-CH11 isolated from Chilika Lake.

An alkaliphile bacterial strain designated as CH11 was isolated from the sediments of Chilika Lake, Odisha. The isolate showed stupendous growth and production of α-amylase at pH 10.0. Through 16S rRNA gene based molecular technique this isolate was identified as Bacillus cereus strain SP-CH11 having GenBank Accession No. KT992791. Homogenous ~55 kDa extracellular α-amylase was extracted with 241.304, 26.26 and 3.2-fold acceleration in specific activity, purification fold and yield respectively. The alkaline α-amylase AA11 was further characterized. At pH 9.0 the purified enzyme AA11 was highly stable while retaining 88-100% functional viability at temperature range from 35 to 65 °C, confirming its thermostability nature. It showed stability with powdered and liquid detergents at 7 mg/mL and 100-fold dilutions respectively. AA11 efficiently removed the starch stain from cotton fabrics. The findings of this study indicate that the isolate CH11 is a source of novel alkaline α-amylase that has promising application in food and detergent industries.

Elucidating the molecular interaction of Zebrafish (Danio rerio) peptidoglycan recognition protein 2 with diaminopimelic acid and lysine type peptidoglycans using in silico approaches.

Peptidoglycan recognition proteins (PGRPs) belong to the family of pattern recognition receptor, represent the major constituent of innate immunity. Although PGRPs are structurally conserved through evolution, their involvement in innate immunity is different in vertebrates and invertebrates. They are highly specific towards recognition of ligands and can hydrolyze bacterial peptidoglycans (PGNs). Zebrafish PGRPs (zPGRPs) have both peptidoglycans lytic amidase activity and broad-spectrum bactericidal activity, but far less is known about how these receptors recognize these microbial ligands. Such studies are hindered due to lack of structural and functional configuration of zPGRPs. Therefore, in this study, we predicted the three-dimensional structure of zPGRP2 through theoretical modeling, investigated the conformational and dynamic properties through molecular dynamics simulations. Molecular docking study revealed the microbial ligands, that is, muramyl pentapeptide-DAP , muramyl pentapeptide-LYS, muramyl tripeptide-DAP, muramyl tripeptide-Lys, muramyl tetrapeptide-DAP, muramyl tetrapeptide-LYS and tracheal cytotoxin interacts with the conserved amino acids of the ligand recognition site comprised of ß1, α2, α4, ß4 and loops connecting ß1 - α2, α2 - ß2, ß3 - ß4 and α4 - α5. Conserved His31, His32, Ala34, Ile35, Pro36, Lys38, Asp60, Trp61, Trp63, Ala89, His90, Asp106, His143 and Arg144 are predicted to essential for binding and provides stability to these zPGRP-PGN complexes. Our study provides basic molecular information for further research on the immune mechanisms of PGRP's in Zebrafish. The plasticity of the zPGRP's binding site revealed by these microbial ligands suggests an intrinsic capacity of the innate immune system to rapidly evolve specificities to meet new microbial challenges in the future.Communicated by Ramaswamy H. Sarma.

Bioactivity guided isolation of antidiabetic and antioxidant compound from Xylocarpus granatum J. Koenig bark.

The present study was designed to identify antidiabetic and antioxidant constituents from ethanol bark extract of a medicinally important mangrove plant Xylocarpus granatum J. Koenig, using in vitro bioactivity-guided fractionation. The repeated fractionation of X. granatum ethanol bark extract (XGEB) by silica gel column chromatography yielded a compound with strong antidiabetic and antioxidant potential. The bioactive compounds likely to be present in the XGEB fraction were identified by FT-IR, 1H & 13C NMR and MS analysis and determined as a limonoid derivative Xyloccensin-I, by comparing spectral data with the literature reports. The isolated compound demonstrated excellent in vitro antidiabetic potential IC50 values of 0.25 and 0.16 mg/ml, respectively for α-amylase and α-glucosidase inhibition study. The antioxidant potential assayed by DPPH, ABTS, superoxide and hydrogen peroxide scavenging studies exhibited that the isolated compound could scavenge these free radicals with IC50 values of 0.041, 0.039, 0.096 and 0.235 mg/ml, respectively. Further, in silico study was performed to find the antidiabetic activity of Xyloccensin-I by docking it against α-glucosidase enzyme. The study demonstrated that Xyloccensin-I have satisfactory interactions and binding energies when docked into target which further confirms the possible mode of antidiabetic action of the isolated compound. The bioactivity assays further asserts the antidiabetic and antioxidant efficacy of the isolated compound which strongly suggests that Xyloccensin-I holds promise in the pharmaceutical industry. The results from this study provide new mechanistic evidence justifying, at least in part, the traditional use of X. granatum extract for antidiabetic and antioxidants activities.

Insights into the mode of flavin mononucleotide binding and catalytic mechanism of bacterial chromate reductases: A molecular dynamics simulation study.

Enzymes from natural sources protect the environment via complex biological mechanisms, which aid in reductive immobilization of toxic metals including chromium. Nevertheless, progress was being made in elucidating high-resolution crystal structures of reductases and their binding with flavin mononucleotide (FMN) to understand the underlying mechanism of chromate reduction. Therefore, herein, we employed molecular dynamics (MD) simulations, principal component analysis (PCA), and binding free energy calculations to understand the dynamics behavior of these enzymes with FMN. Six representative chromate reductases in monomeric and dimeric forms were selected to study the mode, dynamics, and energetic component that drive the FMN binding process. As evidenced by MD simulation, FMN prefers to bind the cervix formed between the catalytic domain surrounded by strong conserved hydrogen bonding, electrostatic, and hydrophobic contacts. The slight movement and reorientation of FMN resulted in breakage of some crucial H-bonds and other nonbonded contacts, which were well compensated with newly formed H-bonds, electrostatic, and hydrophobic interactions. The critical residues aiding in tight anchoring of FMN within dimer were found to be strongly conserved in the bacterial system. The molecular mechanics combined with the Poisson-Boltzmann surface area binding free energy of the monomer portrayed that the van der Waals and electrostatic energy contribute significantly to the total free energy, where, the polar solvation energy opposes the binding of FMN. The proposed proximity relationships between enzyme and FMN binding site presented in this study will open up better avenues to engineer enzymes with optimized chromate reductase activity for sustainable bioremediation of heavy metals.

Interactions between Ibuprofen and Silicified-MCC: Characterization, Drug Release and Modeling Approaches.

Analysis of the binding interactions of ibuprofen and silicified-microcrystalline cellulose (SMCC) has been undertaken. Co-processing of ibuprofen with SMCC was carried out by solid state ball milling, and aqueous state equilibration followed by freeze drying to investigate the effect of silicified-microcrystalline cellulose on ligand. Molecular docking study revealed that ibuprofen formed complex through hydrogen bond with microcrystalline cellulose (MCC) and silicon dioxide (SiO2); the binding energy between MCC and SiO2, and ibuprofen and SMCC were found as -1.11 and -1.73 kcal/mol respectively. The hydrogen bond lengths were varying from 2.028 to 2.056 Å. Interaction of Si atom of SMCC molecule with Pi-Orbital of ibuprofen has shown the bond length of 4.263 Å. Significant improvement in dissolution of ibuprofen has been observed as a result of interaction. Binary and ternary interactions revealed more stabilizing interactions with ibuprofen and SMCC compared to SMCC formation.

Drug-in-mucoadhesive type film for ocular anti-inflammatory potential of amlodipine: Effect of sulphobutyl-ether-beta-cyclodextrin on permeation and molecular docking characterization.

Mucoadhesive type ocular film has been prepared for studying the anti-inflammatory potential of amlodipine (AML) on carrageenan-induced rabbit model and the effect of sulphobutyl-ether-beta-cyclodextrin on corneal permeation was tested. Hydroxypropyl methylcellulose (HPMC) ocular film was prepared after complexation of amlodipine with ß-cyclodextrin, (BCD), hydroxypropyl ß-cyclodextrin (HPCD), and sulfobutylether ß-cyclodextrin (SBCD). The film without cyclodextrin showed a maximum swelling, and erosion to the highest extent. Both drug release and permeation were highly diffusion controlled and highest improvement was observed with SBCD due to increased dissolution, compared to other formulations with or without cyclodextrin. Highest binding energy and highest extent of amorphization were noticed in the SBCD film formulation. Improved amlodipine release in-vitro and ocular permeation were found by the HPMC film formulation after complexation of the drug with cyclodextrins wherein SBCD exhibited both to the highest extent. Binary and ternary systems molecular docking studies confirmed the lowest energy of binding between amlodipine and BCD compared to HBCD and SBCD. Signs of acute inflammation were mitigated within 2 h of film application in the cul-de-sac. Presence of sulphobutyl-ether ß-cyclodextrin in the amlodipine-HPMC film can improve ocular permeation significantly and could be utilized as mucoadhesive type formulation for anti-inflammatory activity.

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.

In Silico Study of miRNA Based Gene Regulation, Involved in Solid Cancer, by the Assistance of Argonaute Protein.

Solid tumor is generally observed in tissues of epithelial or endothelial cells of lung, breast, prostate, pancreases, colorectal, stomach, and bladder, where several genes transcription is regulated by the microRNAs (miRNAs). Argonaute (AGO) protein is a family of protein which assists in miRNAs to bind with mRNAs of the target genes. Hence, study of the binding mechanism between AGO protein and miRNAs, and also with miRNAs-mRNAs duplex is crucial for understanding the RNA silencing mechanism. In the current work, 64 genes and 23 miRNAs have been selected from literatures, whose deregulation is well established in seven types of solid cancer like lung, breast, prostate, pancreases, colorectal, stomach, and bladder cancer. In silico study reveals, miRNAs namely, miR-106a, miR-21, and miR-29b-2 have a strong binding affinity towards PTEN, TGFBR2, and VEGFA genes, respectively, suggested as important factors in RNA silencing mechanism. Furthermore, interaction between AGO protein (PDB ID-3F73, chain A) with selected miRNAs and with miRNAs-mRNAs duplex were studied computationally to understand their binding at molecular level. The residual interaction and hydrogen bonding are inspected in Discovery Studio 3.5 suites. The current investigation throws light on understanding miRNAs based gene silencing mechanism in solid cancer.

Drug Target Identification and Elucidation of Natural Inhibitors for Bordetella petrii: An In Silico Study.

Environmental microbes like Bordetella petrii has been established as a causative agent for various infectious diseases in human. Again, development of drug resistance in B. petrii challenged to combat against the infection. Identification of potential drug target and proposing a novel lead compound against the pathogen has a great aid and value. In this study, bioinformatics tools and technology have been applied to suggest a potential drug target by screening the proteome information of B. petrii DSM 12804 (accession No. PRJNA28135) from genome database of National Centre for Biotechnology information. In this regards, the inhibitory effect of nine natural compounds like ajoene (Allium sativum), allicin (A. sativum), cinnamaldehyde (Cinnamomum cassia), curcumin (Curcuma longa), gallotannin (active component of green tea and red wine), isoorientin (Anthopterus wardii), isovitexin (A. wardii), neral (Melissa officinalis), and vitexin (A. wardii) have been acknowledged with anti-bacterial properties and hence tested against identified drug target of B. petrii by implicating computational approach. The in silico studies revealed the hypothesis that lpxD could be a potential drug target and with recommendation of a strong inhibitory effect of selected natural compounds against infection caused due to B. petrii, would be further validated through in vitro experiments.

Structural and functional insights into CARDs of zebrafish (Danio rerio) NOD1 and NOD2, and their interaction with adaptor protein RIP2.

Nucleotide-binding and oligomerization domain-containing protein 1 (NOD1) and NOD2 are cytosolic pattern-recognition receptors (PRRs) composed of an N-terminal caspase activation and recruitment domain (CARD), a central NACHT domain and C-terminal leucine-rich repeats (LRRs). They play a vital role in innate immune signaling by activating the NF-κB pathway via recognition of peptidoglycans by LRRs, and ATP-dependent self-oligomerization of NACHT followed by downstream signaling. After oligomerization, CARD/s play a crucial role in activating downstream signaling via the adaptor molecule, RIP2. Due to the inadequacy of experimental 3D structures of CARD/s of NOD2 and RIP2, and results from differential experimental setups, the RIP2-mediated CARD-CARD interaction has remained as a contradictory statement. We employed a combinatorial approach involving protein modeling, docking, molecular dynamics simulation, and binding free energy calculation to illuminate the molecular mechanism that shows the possible involvement of either the acidic or basic patch of zebrafish NOD1/2-CARD/a and RIP2-CARD in CARD-CARD interaction. Herein, we have hypothesized 'type-I' mode of CARD-CARD interaction in NOD1 and NOD2, where NOD1/2-CARD/a involve their acidic surfaces to interact with RIP2. Asp37 and Glu51 (of NOD1) and Arg477, Arg521 and Arg529 (of RIP2) were identified to be crucial for NOD1-RIP2 interaction. However, in NOD2-RIP2, Asp32 (of NOD2) and Arg477 and Arg521 (of RIP2) were anticipated to be significant for downstream signaling. Furthermore, we found that strong electrostatic contacts and salt bridges are crucial for protein-protein interactions. Altogether, our study has provided novel insights into the RIP2-mediated CARD-CARD interaction in zebrafish NOD1 and NOD2, which will be helpful to understand the molecular basis of the NOD1/2 signaling mechanism.

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.