Identification of common genetic factors and immune-related pathways associating more than two autoimmune disorders: implications on risk, diagnosis, and treatment.

Autoimmune disorders (ADs) are chronic conditions resulting from failure or breakdown of immunological tolerance, resulting in the host immune system attacking its cells or tissues. Recent studies report shared effects, mechanisms, and evolutionary origins among ADs; however, the possible factors connecting them are unknown. This study attempts to identify gene signatures commonly shared between different autoimmune disorders and elucidate their molecular pathways linking the pathogenesis of these ADs using an integrated gene expression approach. We employed differential gene expression analysis across 19 datasets of whole blood/peripheral blood cell samples with five different autoimmune disorders (rheumatoid arthritis, multiple sclerosis, systemic lupus erythematosus, Crohn's disease, and type 1 diabetes) to get nine key genes-EGR1, RUNX3, SMAD7, NAMPT, S100A9, S100A8, CYBB, GATA2, and MCEMP1 that were primarily involved in cell and leukocyte activation, leukocyte mediated immunity, IL-17, AGE-RAGE signaling in diabetic complications, prion disease, and NOD-like receptor signaling confirming its role in immune-related pathways. Combined with biological interpretations such as gene ontology (GO), pathway enrichment, and protein-protein interaction (PPI) network, our current study sheds light on the in-depth research on early detection, diagnosis, and prognosis of different ADs.

Identification of Potential Genes and Critical Pathways in Postoperative Recurrence of Crohn's Disease by Machine Learning And WGCNA Network Analysis.

Background: Crohn's disease (CD) is a chronic idiopathic inflammatory bowel disease affecting the entire gastrointestinal tract from the mouth to the anus. These patients often experience a period of symptomatic relapse and remission. A 20 - 30% symptomatic recurrence rate is reported in the first year after surgery, with a 10% increase each subsequent year. Thus, surgery is done only to relieve symptoms and not for the complete cure of the disease. The determinants and the genetic factors of this disease recurrence are also not well-defined. Therefore, enhanced diagnostic efficiency and prognostic outcome are critical for confronting CD recurrence. Methods: We analysed ileal mucosa samples collected from neo-terminal ileum six months after surgery (M6=121 samples) from Crohn's disease dataset (GSE186582). The primary aim of this study is to identify the potential genes and critical pathways in post-operative recurrence of Crohn's disease. We combined the differential gene expression analysis with Recursive feature elimination (RFE), a machine learning approach to get five critical genes for the postoperative recurrence of Crohn's disease. The features (genes) selected by different methods were validated using five binary classifiers for recurrence and remission samples: Logistic Regression (LR), Decision tree classifier (DT), Support Vector Machine (SVM), Random Forest classifier (RF), and K-nearest neighbor (KNN) with 10-fold cross-validation. We also performed weighted gene co-expression network analysis (WGCNA) to select specific modules and feature genes associated with Crohn's disease postoperative recurrence, smoking, and biological sex. Combined with other biological interpretations, including Gene Ontology (GO) analysis, pathway enrichment, and protein-protein interaction (PPI) network analysis, our current study sheds light on the in-depth research of CD diagnosis and prognosis in postoperative recurrence. Results: PLOD2, ZNF165, BOK, CX3CR1, and ARMCX4, are the important genes identified from the machine learning approach. These genes are reported to be involved in the viral protein interaction with cytokine and cytokine receptors, lysine degradation, and apoptosis. They are also linked with various cellular and molecular functions such as Peptidyl-lysine hydroxylation, Central nervous system maturation, G protein-coupled chemoattractant receptor activity, BCL-2 homology (BH) domain binding, Gliogenesis and negative regulation of mitochondrial depolarization. WGCNA identified a gene co-expression module that was primarily involved in mitochondrial translational elongation, mitochondrial translational termination, mitochondrial translation, mitochondrial respiratory chain complex, mRNA splicing via spliceosome pathways, etc.; Both the analysis result emphasizes that the mitochondrial depolarization pathway is linked with CD recurrence leading to oxidative stress in promoting inflammation in CD patients. Conclusion: These key genes serve as the novel diagnostic biomarker for the postoperative recurrence of Crohn's disease. Thus, among other treatment options present until now, these biomarkers would provide success in both diagnosis and prognosis, aiming for a long-lasting remission to prevent further complications in CD.

M Protein from Dengue virus oligomerizes to pentameric channel protein: in silico analysis study.

The Dengue virus M protein is a 75 amino acid polypeptide with two helical transmembranes (TM). The TM domain oligomerizes to form an ion channel, facilitating viral release from the host cells. The M protein has a critical role in the virus entry and life cycle, making it a potent drug target. The oligomerization of the monomeric protein was studied using ab initio modeling and molecular dynamics (MD) simulation in an implicit membrane environment. The representative structures obtained showed pentamer as the most stable oligomeric state, resembling an ion channel. Glutamic acid, threonine, serine, tryptophan, alanine, isoleucine form the pore-lining residues of the pentameric channel, conferring an overall negative charge to the channel with approximate length of 51.9 Å. Residue interaction analysis (RIN) for M protein shows that Ala94, Leu95, Ser112, Glu124, and Phe155 are the central hub residues representing the physicochemical interactions between domains. The virtual screening with 165 different ion channel inhibitors from the ion channel library shows monovalent ion channel blockers, namely lumacaftor, glipizide, gliquidone, glisoxepide, and azelnidipine to be the inhibitors with high docking scores. Understanding the three-dimensional structure of M protein will help design therapeutics and vaccines for Dengue infection.

Recursive Feature Elimination-based Biomarker Identification for Open Neural Tube Defects.

Background: Open spina bifida (myelomeningocele) is the result of the failure of spinal cord closing completely and is the second most common and severe birth defect. Open neural tube defects are multifactorial, and the exact molecular mechanism of the pathogenesis is not clear due to disease complexity for which prenatal treatment options remain limited worldwide. Artificial intelligence techniques like machine learning tools have been increasingly used in precision diagnosis. Objective: The primary objective of this study is to identify key genes for open neural tube defects using a machine learning approach that provides additional information about myelomeningocele in order to obtain a more accurate diagnosis. Materials and Methods: Our study reports differential gene expression analysis from multiple datasets (GSE4182 and GSE101141) of amniotic fluid samples with open neural tube defects. The sample outliers in the datasets were detected using principal component analysis (PCA). We report a combination of the differential gene expression analysis with recursive feature elimination (RFE), a machine learning approach to get 4 key genes for open neural tube defects. The features selected were validated using five binary classifiers for diseased and healthy samples: Logistic Regression (LR), Decision tree classifier (DT), Support Vector Machine (SVM), Random Forest classifier (RF), and K-nearest neighbour (KNN) with 5-fold cross-validation. Results: Growth Associated Protein 43 (GAP43), Glial fibrillary acidic protein (GFAP), Repetin (RPTN), and CD44 are the important genes identified in the study. These genes are known to be involved in axon growth, astrocyte differentiation in the central nervous system, post-traumatic brain repair, neuroinflammation, and inflammation-linked neuronal injuries. These key genes represent a promising tool for further studies in the diagnosis and early detection of open neural tube defects. Conclusion: These key biomarkers help in the diagnosis and early detection of open neural tube defects, thus evaluating the progress and seriousness in diseases condition. This study strengthens previous literature sources of confirming these biomarkers linked with open NTD's. Thus, among other prenatal treatment options present until now, these biomarkers help in the early detection of open neural tube defects, which provides success in both treatment and prevention of these defects in the advanced stage.

Amino Acid and Peptide-Based Liquid Crystals: An Overview.

The role of amino acids and peptides has found remarkable usage in both living systems and nonliving materials, which have enabled its utility by virtue of crafting molecular architectures through covalent bonds and non-covalent interactions. In material chemistry, the role of peptides in Liquid Crystals (LCs) is profound, especially in the rapid construction of supramolecular hierarchical networks. The importance of LCs for a variety of societal needs leads to the synthesis of innumerable LCs by conventional mesogenic strategy and nonconventional molecular design principles. For example, electronic appliances, including flat panel TV displays, electronic notebooks, digital cameras, domestic devices, use LCs as an integral component for such applications. In addition, LCs are useful in biological systems, including stem cell research, sensors for bacteria, virus, and proteins. These accomplishments are possible mostly due to the non-conventional molecular design principles for crafting LCs using smaller molecular motifs. The usage of amino acids and peptides in LCs facilitates many intrinsic characteristics, including side-chain diversity, chirality, directionality, reversibility, electro-optical, columnar axis, stimuli-responsive complex molecular architectures. The next essential criteria for any LCs design for useful applications are room temperature LC (RT-LC); therefore, the quest for such LCs system remains highly significant. Evidently, there are around half a million liquid crystalline molecules; only a handful of RTLCs has been found, as there is no simple, precise strategy or molecular design principles to obtain RT-LC systems. The smaller molecular motifs of amino acids and linear peptides as a structural part of mesogenic molecules led to many LC phases with properties, including lyotropic, thermotropic, and its applications in different realms. Therefore, this review serves as a compilation of Small Peptide-based LCs (SPLCs) exhibiting lyotropic and thermotropic phases with applications in the recent advancements.

Molecular Interaction between Distal C-Terminal Domain of the CB1 Cannabinoid Receptor and Cannabinoid Receptor Interacting Proteins (CRIP1a/CRIP1b).

We have investigated the structure of the distal C-terminal domain of the of the CB1 cannabinoid receptor (CB1R) to study its interactions with CRIP1a and CRIP1b using computational techniques. The amino acid sequence from the distal C-terminal domain of CB1R (G417-L472) was found to be unique, as it does not share sequence similarity with other protein structures, so the structure was predicted using ab initio modeling. The computed model of the distal C-terminal region of CB1R has a helical region between positions 441 and 455. The CRIP1a and CRIP1b were modeled using Rho-GDI 2 as a template. The three-dimensional model of the distal C-terminal domain of the CB1R was docked with both CRIP1a as well as CRIP1b to study the crucial interactions between CB1R and CRIP1a/b. The last nine residues of CB1R (S464TDTSAEAL4722) are known to be a CRIP1a/b binding site. The majority of the key interactions were identified in this region, but notable interactions were also observed beyond theses nine residues. The multiple interactions between Thr418 (CB1R) and Asn61 (CRIP1a) as well as Asp430 (CB1R) and Lys76 (CRIP1a) indicate their importance in the CB1R-CRIP1a interaction. In the case of CRIP1b, multiple hydrogen bond interactions between Asn437 (CB1R) and Glu77 (CRIP1b) were observed. These interactions can be critical for CB1R's interaction with CRIP1a/b, and targeting them for further experimental studies can advance information about CRIP1a/b functionality.

CB1 cannabinoid receptor-phosphorylated fourth intracellular loop structure-function relationships.

A peptide comprising the juxtamembrane C-terminal intracellular loop 4 (IL4) of the CB1 cannabinoid receptor possesses three Serine residues (Ser402, Ser411 and Ser415). Here we report the effect of Ser phosphorylation on the CB1 IL4 peptide conformation and cellular signaling functions using nuclear magnetic resonance spectroscopy, circular dichroism, G protein activation and cAMP production. Circular dichroism studies indicated that phosphorylation at various Ser residues induced helical structure in different environments. NMR data indicates that helical content varies in the order of IL4pSer411 > IL4pSer415 > IL4 > IL4pSer402. The efficacy of phosphorylated IL4 peptides in activating Go and Gi3 ([35S]GTPγS binding) and inhibiting cAMP accumulation in N18TG2 cells were correlated with helicity changes. Treatment of cells with bradykinin, which activates PKC, augmented CB1-mediated inhibition of cAMP accumulation, and this was reversed by a PKC inhibitor, suggesting that phosphorylation of serine might be a physiologically relevant modification in vivo. We conclude that phosphorylation-dependent alterations of helicity of CB1 IL4 peptides can increase efficacy of G protein signaling.

In silico interaction analysis of cannabinoid receptor interacting protein 1b (CRIP1b) - CB1 cannabinoid receptor.

Cannabinoid Receptor Interacting Protein isoform 1b (CRIP1b) is known to interact with the CB1 receptor. Alternative splicing of the CNRIP1 gene produces CRIP1a and CRIP1b with a difference in the third exon only. Exons 1 and 2 encode for a functional domain in both proteins. CRIP1a is involved in regulating CB1 receptor internalization, but the function of CRIP1b is not very well characterized. Since there are significant identities in functional domains of these proteins, CRIP1b is a potential target for drug discovery. We report here predicted structure of CRIP1b followed by its interaction analysis with CB1 receptor by in-silico methods A number of complementary computational techniques, including, homology modeling, ab-initio and protein threading, were applied to generate three-dimensional molecular models for CRIP1b. The computed model of CRIP1b was refined, followed by docking with C terminus of CB1 receptor to generate a model for the CRIP1b- CB1 receptor interaction. The structure of CRIP1b obtained by homology modelling using RHO_GDI-2 as template is a sandwich fold structure having beta sheets connected by loops, similar to predicted CRIP1a structure. The best scoring refined model of CRIP1b in complex with the CB1 receptor C terminus peptide showed favourable polar interactions. The overall binding pocket of CRIP1b was found to be overlapping to that of CRIP1a. The Arg82 and Cys126 of CRIP1b are involved in the majority of hydrogen bond interactions with the CB1 receptor and are possible key residues required for interactions between the CB1 receptor and CRIP1b.

Membrane-induced structure of novel human tachykinin hemokinin-1 (hHK1).

PPT-C encoded hemokinin-1(hHK-1) of Homo sapiens (TGKASQFFGLM) is a structurally distinct neuropeptide among the tachykinin family that participate in the NK-1 receptor downstream signaling processes. Subsequently, signal transduction leads to execution of various effector functions which includes aging, immunological, and central nervous system (CNS) regulatory actions. However the conformational pattern of ligand receptor binding is unclear. The three-dimensional structure of the hemokinin-1 in aqueous and micellar environment has been studied by one and two-dimensional proton nuclear magnetic resonance (2D 1H-NMR spectroscopy) and distance geometry calculations. Data shows that hemokinin-1 was unstructured in aqueous environment; anionic detergent SDS induces α-helix formation. Proton NMR assignments have been carried out with the aid of correlation spectroscopy (gradient-COSY and TOCSY) and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The inter proton distances and dihedral angle constraints obtained from the NMR data have been used in torsion angle dynamics algorithm for NMR applications (CYANA) to generate a family of structures, which have been refined using restrained energy minimization and dynamics. Helical conformation is observed from residue K3-M11. The conformational range of the peptide revealed by NMR studies has been analyzed in terms of characteristic secondary features. Observed conformational features have been compared to that of Substance P potent NK1 agonist. Thus the report provides a structural insight to study hHK-1-NK1 interaction that is essential for hHK1 based signaling events.

Effect of methylene group insertions on the structural rigidity of Aib containing helices.

Nonprotein amino acids are being extensively used in the design of synthetic peptides to create new structure mimics. In this study we report the effect of methylene group insertions in a heptapeptide Boc-Ala1-Leu2-Aib3-Xxx4-Ala5-Leu6-Aib7-OMe which nicely folds into a mixed 310 -/α-helical structure when Xxx= Ala. Analogs of this peptide have been made and studied by replacing central Xxx4 residue with Glycine (α-residue), ß-Alanine (ß-Αla), γ-aminobutyric acid (Gaba), and ε-aminocaproic acid (ε-Aca). NMR and circular dichroism were used to study the solution structure of these peptides. Crystals of the peptides containing alanine, ß-Αla, and Gaba reveal that increasing the number of central methylene (-CH2 -) groups introduces local perturbations even as the helical structure is retained.

Proteomic profiling of serum samples from chikungunya-infected patients provides insights into host response.

BACKGROUND: Chikungunya is a highly debilitating febrile illness caused by Chikungunya virus, a single-stranded RNA virus, which is transmitted by Aedes aegypti or Aedes albopictus mosquito species. The pathogenesis and host responses in individuals infected with the chikungunya virus are not well understood at the molecular level. We carried out proteomic profiling of serum samples from chikungunya patients in order to identify molecules associated with the host response to infection by this virus. RESULTS: Proteomic profiling of serum obtained from the infected individuals resulted in identification of 569 proteins. Of these, 63 proteins were found to be differentially expressed (≥ 2-fold) in patient as compared to control sera. These differentially expressed proteins were involved in various processes such as lipid metabolism, immune response, transport, signal transduction and apoptosis. CONCLUSIONS: This is the first report providing a global proteomic profile of serum samples from individuals infected with the chikungunya virus. Our data provide an insight into the proteins that are involved as host response factors during an infection. These proteins include clusterin, apolipoproteins and S100A family of proteins.

Molecular recognition of tachykinin receptor selective agonists: insights from structural studies.

This Review deals essentially with the elucidation of structural features of Tachykinin family of neuropeptides, which are known to interact through three distinct GPCR subtypes, namely NK1 (Neurokinin 1), NK2 (Neurokinin 2) and NK3 (Neurokinin 3) receptors. In mammals, Tachykinins have been shown to elicit a wide array of activities such as powerful vasodilatation, hypertensive action and stimulation of extravascular smooth muscle and are known to be involved in a variety of clinical conditions including chronic pain, Parkinson's disease, Alzheimer's disease, depression, rheumatoid arthritis, irritable bowel syndrome and asthma. This broad spectrum of action of Tachykinins is attributed to the lack of selectivity of tachykinins to their receptors. All tachykinins interact with all the three-receptor subtypes with SP preferring NK1, NKA preferring NK2 and NKB preferring NK3. This lack of specificity can be accounted for by the conformational flexibility of these short, linear peptides. Hence, identification of structural features of the agonists important for receptor binding and biological activity is of great significance in unraveling the molecular mechanisms involved in tachykinin receptor activation and also in rational design of novel therapeutic agents. Understanding structure of the ligand-receptor complex and analysis of topography of the binding pocket of the tachykinin receptor is also crucial in rational design of drugs.

Proteomic analysis of purified protein derivative of Mycobacterium tuberculosis.

BACKGROUND: Purified protein derivative (PPD) has been used for more than half a century as an antigen for the diagnosis of tuberculosis infection based on delayed type hypersensitivity. Although designated as "purified," in reality, the composition of PPD is highly complex and remains ill-defined. In this report, high resolution mass spectrometry was applied to understand the complexity of its constituent components. A comparative proteomic analysis of various PPD preparations and their functional characterization is likely to help in short-listing the relevant antigens required to prepare a less complex and more potent reagent for diagnostic purposes. RESULTS: Proteomic analysis of Connaught Tuberculin 68 (PPD-CT68), a tuberculin preparation generated from M. tuberculosis, was carried out in this study. PPD-CT68 is the protein component of a commercially available tuberculin preparation, Tubersol, which is used for tuberculin skin testing. Using a high resolution LTQ-Orbitrap Velos mass spectrometer, we identified 265 different proteins. The identified proteins were compared with those identified from PPD M. bovis, PPD M. avium and PPD-S2 from previous mass spectrometry-based studies. In all, 142 proteins were found to be shared between PPD-CT68 and PPD-S2 preparations. Out of the 354 proteins from M. tuberculosis-derived PPDs (i.e. proteins in either PPD-CT68 or PPD-S2), 37 proteins were found to be shared with M. avium PPD and 80 were shared with M. bovis PPD. Alignment of PPD-CT68 proteins with proteins encoded by 24 lung infecting bacteria revealed a number of similar proteins (206 bacterial proteins shared epitopes with 47 PPD-CT68 proteins), which could potentially be involved in causing cross-reactivity. The data have been deposited to the ProteomeXchange with identifier PXD000377. CONCLUSIONS: Proteomic and bioinformatics analysis of different PPD preparations revealed commonly and differentially represented proteins. This information could help in delineating the relevant antigens represented in various PPDs, which could further lead to development of a lesser complex and better defined skin test antigen with a higher specificity and sensitivity.

Theoretical and in vitro studies of a C-terminal peptide from PGKC of Leishmania mexicana mexicana.

Trypanosomatids cause deadly diseases in humans. Of the various biochemical pathways in trypanosomatids, glycolysis, has received special attention because of being sequestered in peroxisome like organelles critical for the survival of the parasites. This study focuses on phosphoglycerate kinase (PGK) from Leishmania spp. which, exists in two isoforms, the cytoplasmic PGKB and glycosomal PGKC differing in their biochemical properties. Computational analysis predicted the likelihood of a transmembrane helix only in the glycosomal isoform PGKC, of approximate length 20 residues in the 62-residue extension, ending at, arginine residues R471 and R472. From experimental studies using circular dichroism and NMR with deuterated sodium dodecyl sulfate, we find that the transmembrane helix spans residues 448±2 to 476 in Leishmania mexicana PGKC. The significance of this observation is discussed in the context of glycosomal transport and substrate tunneling.

Molecular modeling of neurokinin B and tachykinin NK3 receptor complex.

The tachykinin receptor NK3 is a member of the rhodopsin family of G-protein coupled receptors. The NK3 receptor has been regarded as an important drug target due to diverse physiological functions and its possible role in the pathophysiology of psychiatric disorders, including schizophrenia. The NK3 receptor is primarily activated by the tachykinin peptide hormone neurokinin B (NKB) which is the most potent natural agonist for the NK3 receptor. NKB has been reported to play a vital role in the normal human reproduction pathway and in potentially life threatening diseases such as pre-eclampsia and as a neuroprotective agent in the case of neurodegenerative diseases. Agonist binding to the receptor is a critical event in initiating signaling, and therefore a characterization of the structural features of the agonists can reveal the molecular basis of receptor activation and help in rational design of novel therapeutics. In this study a molecular model for the interaction of the primary ligand NKB with its G-protein coupled receptor NK3 has been developed. A three-dimensional model for the NK3 receptor has been generated by homology modeling using rhodopsin as a template. A knowledge based docking of the NMR derived bioactive conformation of NKB to the receptor has been performed utilizing limited ligand binding data obtained from photoaffinity labeling and site-directed mutagenesis studies. A molecular model for the NKB-NK3 receptor complex obtained sheds light on the topographical features of the binding pocket of the receptor and provides insight into the biochemical data currently available for the receptor.

Three-dimensional structure of Phyllomedusin, a NK1 receptor agonist bound to dodecylphosphocholine micelles.

Phyllomedusin, an amphibian tachykinin decapeptide, has been shown to be selective for Neurokinin 1 receptor. Because the micelle-associated structure may be relevant to the Phyllomedusin-receptor interaction, the three-dimensional structure of the Phyllomedusin in aqueous and micellar environments has been studied by two-dimensional proton nuclear magnetic resonance (2D (1)H NMR spectroscopy) and distance geometry calculations. Sequence specific resonance assignments of protons have been made from correlation spectroscopy (TOCSY, DQF-COSY) and NOESY spectroscopy. The interproton distance constraints and dihedral angle constraints have been utilized to generate a family of structures using DYANA. The CD and NMR results show that, while in water Phyllomedusin prefers to be in an extended chain conformation, whereas in the presence of dodecylphosphocholine micelles, a membrane model system, a partial helical conformation is induced. Analysis of NMR data indicates that the global fold of Phyllomedusin can be explained in terms of equilibrium between 3(10)-helix and alpha-helix from residue 4 to 10. An extended highly flexible N-terminus displays some degree of order and a possible turn structure. A comparison between the conformational features of Phyllomedusin and different Neurokinin 1 receptor agonist indicates several common features in the distribution of hydrophobic and hydrophilic residues. The conformational similarities suggest that the molecules interact with receptor in an analogous manner.