Search

Abstract: The common feature of receptors and antibodies against the ligand is that both display very specific, high affinity binding towards the ligand. Therefore, it can be hypothesized that the paratope of antibodies may exhibit homology with distinct domains of the receptor. By locating the hormone epitopes and determining the structure of the paratopes, it should be possible to identify the contact points between the ligand and the receptor. This hypothesis has been tested using hCG monoclonal antibodies (MAbs) recognizing different epitopes and having different effects on hormone binding and response. The β subunit and heterodimer specific antibodies inhibited both hormone binding and response, while the α subunit specific antibodies inhibited response without affecting binding. The single chain fragment variables (ScFvs) produced from these antibodies also retained the properties of the parent antibodies. The amino acid sequences of the ScFvs exhibited homology to different regions of the receptor; the β subunit specific antibody being homologous to the concave surface of the leucine rich repeats (LRR) of the receptor, particularly the concave surface of the LRRs, while the heterodimer specific antibody showed homology to the hinge region. The α subunit specific antibody showed homology to the transmembrane domain of the receptor. The exact locations of the epitopes of the monoclonal antibodies in the hormone molecule have also been identified. The data presented here also support the model of glycoprotein hormone–receptor interaction in which the hormone binds to the extracellular domain through the β subunit and then the α subunit is brought in contact with the transmembrane domain leading to signal transduction. [Copyright &y& Elsevier]

Background: Bcl-2 family proteins are key regulators of mitochondrial integrity and comprise both pro-and anti-apoptotic proteins. Bax a pro-apoptotic member localizes as monomers in the cytosol of healthy cells and accumulates as oligomers in mitochondria of apoptotic cells. The Bcl-2 homology-3 (BH3) domain regulates interactions within the family, but regions other than BH3 are also critical for Bax function. Thus, the N-terminus has been variously implicated in targeting to mitochondria, interactions with BH3-only proteins as well as conformational changes linked to Bax activation. The transmembrane (TM) domains (a5-a6 helices in the core and a9 helix in the C-terminus) in Bax are implicated in localization to mitochondria and triggering cytotoxicity. Here we have investigated N-terminus modulation of TM function in the context of regulation by the anti-apoptotic protein Bcl-xL. Results: Deletion of 29 amino acids in the Bax N-terminus (Bax 30-192) caused constitutive accumulation at mitochondria and triggered high levels of cytotoxicity, not inhibited by Bcl-xL. Removal of the TM domains (Bax 30-105) abrogated mitochondrial localization but resulted in Bcl-xL regulated activation of endogenous Bax and Bax-Bak dependent apoptosis. Inclusion of the a5-a6 helices/TMI domain (Bax 30-146) phenocopied Bax 30-192 as it restored mitochondrial localization, Bcl-xL independent cytotoxicity and was not dependent on endogenous Bax-Bak. Inhibition of function and localization by Bcl-xL was restored in Bax 1-146, which included the TM1 domain. Regardless of regulation by Bcl-xL, all N-terminal deleted constructs immunoprecipitated Bcl-xLand converged on caspase-9 dependent apoptosis consistent with mitochondrial involvement in the apoptotic cascade. Sub-optimal sequence alignments of Bax and Bcl-xL indicated a sequence similarity between the a5-a6 helices of Bax and Bcl-xL. Alanine substitutions of three residues (T14A-S15A-S16A) in the N-terminus (Bax-Ala3) attenuated regulation by the serine-threonine kinase Akt/PKB but not by Bcl-xL indicative of distinct regulatory mechanisms. Conclusion: Collectively, the analysis of Bax deletion constructs indicates that the N-terminus drives conformational changes facilitating inhibition of cytotoxicity by Bcl-xL. We speculate that the TM1 helices may serve as 'structural antagonists' for BH3-Bcl-xL interactions, with this function being regulated by the N-terminus in the intact protein. [ABSTRACT FROM AUTHOR]

Resource-poor areas with moisture deficit lands following aerobic and direct seeded rice (DSR) methods of cultivation face severe problems of iron deficiency. In this study, Bengal and Assam Aus rice panel was phenotyped at the seedling stage using an iron-deprived hydroponic medium for various shoot and root traits. A novel iron deficiency scoring scale was used to classify the tolerance reaction and could range anywhere between 0 and 9, indicating the most tolerant and susceptible, respectively. The GWAS results identified four putative candidate genes; OsFLA for number of leaves and shoot length, OsBIDK1 for root traits; average diameter, volume, biomass, projected area, and surface area, OsHPL3 for chlorophyll index of the third leaf and AKR2B (XBOS252) was for Fe score, (which was earlier reported in relation to Xa21). The nsSNP (nsSNPs) variations in these gene sequences were used to group the panel and identify superior haplotypes and donors. BR16 was identified as a superior donor, with higher chlorophyll index and shoot length than RA23, also higher values for root traits like root average diameter, root volume, root projected area and root surface area followed by Shete Bhado. The impact of identified nsSNPs on protein structure and stability was investigated. The conserved domains detected in the mutated proteins of the superior haplotypes are very informative, highlighting that natural selection favors abiotic stress tolerant variants in resource poor areas. Thus, justifying our choice of Aus landraces for association mapping of Fe deficiency tolerant genes in rice., (© 2024. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Toll-like receptors (TLRs) are major players in the innate immune system-recognizing pathogens and differentiating self/non-self components of immunity. These proteins are present either on the plasma membrane or endosome and recognize pathogens at their extracellular domains. They are characterized by a single transmembrane helix and an intracellular toll-interleukin-1 receptor (TIR) domain. Few TIRs directly invoke downstream signaling, while others require other TIR domains of adaptors like TIR domain-containing adaptor-inducing interferon-β (TRIF) and TRIF-related adaptor molecule (TRAM). On recognizing pathogenic lipopolysaccharides, TLR4 dimerises and interacts with the intracellular TRAM dimer through the TIR domain to recruit a downstream signaling adaptor (TRIF). We have performed an in-depth study of the structural effect of two mutations (P116H and C117H) at the dimeric interface of the adaptor TRAM, which are known to abrogate downstream signaling. We modeled the structure and performed molecular dynamics studies in order to decipher the structural basis of this effect. We observed that these mutations led to an increased radius of gyration of the complex and resulted in several changes to the interaction energy values when compared against the wild type (WT) and positive control mutants. We identified highly interacting residues as hubs in the WT dimer, and a few such hubs that were lost in the mutant dimers. Changes in the protein residue path, hampering the information flow between the crucial A86/E87/D88/D89 and T155/S156 sites, were observed for the mutants. Overall, we show that such residue changes can have subtle but long-distance effects, impacting the signaling path allosterically., (© 2024 The Author(s). Proteins: Structure, Function, and Bioinformatics published by Wiley Periodicals LLC.)

Single stranded RNA viruses are primary causative agents for pandemics, causing extensive morbidity and mortality worldwide. A pivotal question in pandemic preparedness and therapeutic intervention is what are the specific mutations which are more likely to emerge during such global health crises? This study aims to identify markers for mutations with the highest probability of emergence in these pandemics, focusing on the SARS-CoV-2 spike protein, an essential and therapeutically significant viral protein, starting from sequence information from the onset of the pandemic until July 2022. Quite consistently, we observed that emerged mutations tended to demonstrate a high genetic score, which reflects high similarity of the type of codon required for translation between an amino acid and to the mutated one. Further, this pattern is also observed in therapeutically significant proteins of other ssRNA pandemic viruses, including influenza (HA, NA), spike proteins of Ebola, envelope of Dengue and Chikungunya. We propose that the genetic score serves as an initial indicator, preceding the actual impact of the mutation on viral fitness. Finally, we developed a comprehensive computational pipeline to further explore and predict the subsequent effects of mutations on viral fitness. We believe that our pipeline can narrow down and predict future mutations in therapeutically important viral proteins during a pandemic., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Introduction: Black pepper (Piper nigrum L.) is a non-model spice crop of significant agricultural and biological importance. The 'quick wilt' disease caused by the oomycete Phytophthora capsici is a major threat, leading to substantial crop loss. The molecular mechanisms governing the plant immune responses to this pathogen remain unclear. This study employs RNA sequencing and transcriptome analysis to explore the defense mechanisms of P. nigrum against P. capsici., Results: Two-month-old P. nigrum plantlets were subjected to infection with P. capsici, and leaf samples were collected at 6- and 12-hours post-inoculation. RNA was extracted, sequenced, and the resulting data were processed and assembled. Differential gene expression analysis was conducted to identify genes responding to the infection. Additionally, the study investigated the involvement of Salicylic acid (SA), Jasmonic acid (JA), and Ethylene (ET) signalling pathways. Our transcriptome assembly comprised 64,667 transcripts with 96.7% completeness, providing valuable insights into the P. nigrum transcriptome. Annotation of these transcripts identified functional categories and domains, provided details on molecular processes. Gene expression analysis identified 4,714 transcripts at 6 h post-infection (hpi) and 9,416 at 12 hpi as differentially expressed, revealing dynamic regulation of immune-related genes. Furthermore, the study investigated key genes involved in biosynthesis pathways of Salicylic acid, Jasmonic acid, and Ethylene signalling. Notably, we found differential regulation of critical genes associated with these pathways while comparing data before and after infection, thereby shedding light on their roles in defense mechanism in P. nigrum defense., Conclusions: This comprehensive transcriptome analysis of P. nigrum response to P. capsici attack provides valuable insights into the plant defense mechanisms. The dynamic regulation of innate immunity and the involvement of key signalling pathways highlight the complexity of the plant-pathogen interaction. This study contributes to our understanding of plant immunity and offers potential strategies for enhancing P. nigrum resistance to this harmful pathogen., (© 2024. The Author(s).)

An essential adaptive strategy in insects is the evolution of olfactory receptors (ORs) to recognize important volatile environmental chemical cues. Our model species, Ceratosolen fusciceps, a specialist wasp pollinator of Ficus racemosa, likely possesses an OR repertoire that allows it to distinguish fig-specific volatiles in highly variable environments. Using a newly assembled genome-guided transcriptome, we annotated 63 ORs in the species and reconstructed the phylogeny of Ceratosolen ORs in conjunction with other hymenopteran species. Expression analysis showed that though ORs were mainly expressed in the female antennae, 20% were also expressed in nonantennal tissues such as the head, thorax, abdomen, legs, wings, and ovipositor. Specific upregulated expression was observed in OR30C in the head and OR60C in the wings. We identified OR expression from all major body parts of female C. fusciceps, suggesting novel roles of ORs throughout the body. Further examination of the OR expression of C. fusciceps in widely separated geographical locations, that is, South (urban) and Northeast (rural) India, revealed distinct OR expression levels in different locations. This discrepancy likely parallels the observed variation in fig volatiles between these regions and provides new insights into the evolution of insect ORs and their expression across geographical locations and tissues., (© 2024 Institute of Zoology, Chinese Academy of Sciences.)

Celiac disease (CD) is a common autoimmune disorder in which the patients are unable to digest gluten, which is present in foods made up of wheat, barley and rye. Whilst diagnosis happens late in 80% of the cases, avoidance of such foods appears to be the common solution. Alternative management strategies are required for the patients and their families since CD is also genetically carried over. Probiotic therapeutics and the consumption of appropriate enzymes, such as prolyloligopeptidases (POPs), from gut-friendly bacteria could reduce the disease burden and provide a better lifestyle for CD patients. We have examined around 5000 gut bacterial genomes and identified nearly 4000 non-redundant putative POPs. A select set of 10 gut bacterial POP sequences were subject to three-dimensional modelling, ligand docking and molecular dynamics simulations where stable interactions were observed between the POPs and gluten peptides. Our study provides sequence and structural analysis of potential POP enzymes in gut bacterial genomes, which form a strong basis to offer probiotic solutions to CD patients. In particular, these enzymes could be lead future therapeutics for this disease., (© 2024. The Author(s).)

Phosphorylation, a fundamental biochemical switch, intricately regulates protein function and signaling pathways. Our study employs extensive computational structural analyses on a curated data set of phosphorylated and unphosphorylated protein pairs to explore the multifaceted impact of phosphorylation on protein conformation. Using normal mode analysis (NMA), we investigated changes in protein flexibility post-phosphorylation, highlighting an enhanced level of structural dynamism. Our findings reveal that phosphorylation induces not only local changes at the phosphorylation site but also extensive alterations in distant regions, showcasing its far-reaching influence on protein structure-dynamics. Through in-depth case studies on polyubiquitin B and glycogen synthase kinase-3 beta, we elucidate how phosphorylation at distinct sites leads to variable structural and dynamic modifications, potentially dictating functional outcomes. While phosphorylation largely preserves the residue motion correlation, it significantly disrupts low-frequency global modes, presenting a dualistic impact on protein dynamics. We also explored alterations in the total accessible surface area (ASA), emphasizing region-specific changes around phosphorylation sites. This study sheds light on phosphorylation-induced conformational changes, dynamic modulation, and surface accessibility alterations, leveraging an integrated computational approach with RMSD, NMA, and ASA, thereby contributing to a comprehensive understanding of cellular regulation and suggesting promising avenues for therapeutic interventions., Competing Interests: The authors declare no competing financial interest., (© 2024 The Authors. Published by American Chemical Society.)

Interactions between proteins are vital in almost all biological processes. The characterization of protein-protein interactions helps us understand the mechanistic basis of biological processes, thereby enabling the manipulation of proteins for biotechnological and clinical purposes. The interface residues of a protein-protein complex are assumed to have the following two properties: (a) they always interact with a residue of a partner protein, which forms the basis for distance-based interface residue identification methods, and (b) they are solvent-exposed in the isolated form of the protein and become buried in the complex form, which forms the basis for Accessible Surface Area (ASA)-based methods. The study interrogates this popular assumption by recognizing interface residues in protein-protein complexes through these two methods. The results show that a few residues are identified uniquely by each method, and the extent of conservation, propensities, and their contribution to the stability of protein-protein interaction varies substantially between these residues. The case study analyses showed that interface residues, unique to distance, participate in crucial interactions that hold the proteins together, whereas the interface residues unique to the ASA method have a potential role in the recognition, dynamics, and specificity of the complex and can also be a hotspot. Overall, the study recommends applying both distance and ASA methods so that some interface residues missed by either method but crucial to the stability, recognition, dynamics, and function of protein-protein complexes are identified in a complementary manner., (© 2023 Wiley Periodicals LLC.)

Proteins such as enzymes perform their function by predominant non-covalent bond interactions between transiently interacting units. There is an impact on the overall structural topology of the protein, albeit transient nature of such interactions, that enable proteins to deactivate or activate. This aspect of the alteration of the structural topology is studied by employing protein structural networks, which are node-edge representative models of protein structure, reported as a robust tool for capturing interactions between residues. Several methods have been optimized to collect meaningful, functionally relevant information by studying alteration of structural networks. In this article, different methods of comparing protein structural networks are employed, along with spectral decomposition of graphs to study the subtle impact of protein-protein interactions. A detailed analysis of the structural network of interacting partners is performed across a dataset of around 900 pairs of bound complexes and corresponding unbound protein structures. The variation in network parameters at, around, and far away from the interface are analyzed. Finally, we present interesting case studies, where an allosteric mechanism of structural impact is understood from communication-path detection methods. The results of this analysis are beneficial in understanding protein stability, for future engineering, and docking studies., (© 2023 Wiley Periodicals LLC.)

Protein domains are structural, functional, and evolutionary units. These domains bring out the diversity of functionality by means of interactions with other co-existing domains and provide stability. Hence, it is important to study intra-protein inter-domain interactions from the perspective of types of interactions. Domains within a chain could interact over short timeframes or permanently, rather like protein-protein interactions (PPIs). However, no systematic study has been carried out between two classes, namely permanent and transient domain-domain interactions. In this work, we studied 263 two-domain proteins, belonging to either of these classes and their interfaces on the basis of several factors, such as interface area and details of interactions (number, strength, and types of interactions). We also characterized them based on residue conservation at the interface, correlation of residue motions across domains, its involvement in repeat formation, and their involvement in particular molecular processes. Finally, we could analyze the interactions arising from domains in two-domain monomeric proteins, and we observed significant differences between these two classes of domain interactions and a few similarities. This study will help to obtain a better understanding of structure-function and folding principles of multi-domain proteins., (© 2023 Wiley Periodicals LLC.)

Insect Odorant Binding Proteins (OBPs) constitute important components of their olfactory apparatus, as they are essential for odor recognition. OBPs undergo conformational changes upon pH change, altering their interactions with odorants. Moreover, they can form heterodimers with novel binding characteristics. Anopheles gambiae OBP1 and OBP4 were found capable of forming heterodimers possibly involved in the specific perception of the attractant indole. In order to understand how these OBPs interact in the presence of indole and to investigate the likelihood of a pH-dependent heterodimerization mechanism, the crystal structures of OBP4 at pH 4.6 and 8.5 were determined. Structural comparison to each other and with the OBP4-indole complex (3Q8I, pH 6.85) revealed a flexible N-terminus and conformational changes in the α4-loop-α5 region at acidic pH. Fluorescence competition assays showed a weak binding of indole to OBP4 that becomes further impaired at acidic pH. Additional Molecular Dynamic and Differential Scanning Calorimetry studies displayed that the influence of pH on OBP4 stability is significant compared to the modest effect of indole. Furthermore, OBP1-OBP4 heterodimeric models were generated at pH 4.5, 6.5, and 8.5, and compared concerning their interface energy and cross-correlated motions in the absence and presence of indole. The results indicate that the increase in pH may induce the stabilization of OBP4 by increasing its helicity, thereby enabling indole binding at neutral pH that further stabilizes the protein and possibly promotes the creation of a binding site for OBP1. A decrease in interface stability and loss of correlated motions upon transition to acidic pH may provoke the heterodimeric dissociation allowing indole release. Finally, we propose a potential OBP1-OBP4 heterodimer formation/disruption mechanism induced by pH change and indole binding., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier B.V. All rights reserved.)

Toll-like receptors (TLRs) are pattern recognition receptors present on the surface of cells playing a crucial role in innate immunity. One of the TLRs, TLR4, recognizes LPS (Lipopolysaccharide) as its ligand leading to the release of anti-inflammatory mediators as well as pro-inflammatory cytokines through signal transduction and domain recruitment. TLR4 homodimerizes at its intracellular TIR (Toll/interleukin-1 receptor) domain that helps in the recruitment of the TRAM/TICAM2 (TIR domain-containing adaptor molecule 2) molecule. TRAM also contains TIR domain which in turn, dimerizes and functions as an adapter protein to further recruit TRIF/TICAM1 (TIR domain-containing adaptor molecule 1) protein for mediating downstream signaling. Apart from LPS, TLR4 also recognizes endogenous ligands like fibrinogen, HMGB1, and hyaluronan in autoimmune conditions and sepsis. We employed computational approaches to target TRAM and recognize small molecule inhibitors from small molecules of natural origin, as contained in the Super Natural II database. Finally, cell reporter assays and NMR studies enabled the identification of promising lead compounds. Hence, this study aims to attenuate the signaling of the TLR4-TRAM-TRIF cascade in these auto-inflammatory conditions., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (© 2023 The Authors. Published by Elsevier B.V. on behalf of Research Network of Computational and Structural Biotechnology.)

Introduction: Cardiomyopathies are complex heart diseases with significant prevalence around the world. Among these, primary forms are the major contributors to heart failure and sudden cardiac death. As a high-energy demanding engine, the heart utilizes fatty acids, glucose, amino acid, lactate and ketone bodies for energy to meet its requirement. However, continuous myocardial stress and cardiomyopathies drive towards metabolic impairment that advances heart failure (HF) pathogenesis. So far, metabolic profile correlation across different cardiomyopathies remains poorly understood., Methods: In this study, we systematically explore metabolic differences amongst primary cardiomyopathies. By assessing the metabolic gene expression of all primary cardiomyopathies, we highlight the significantly shared and distinct metabolic pathways that may represent specialized adaptations to unique cellular demands. We utilized publicly available RNA-seq datasets to profile global changes in the above diseases (| log2FC | ≥ 0.28 and BH  adjusted p-val 0.1) and performed gene set analysis (GSA) using the PAGE statistics on KEGG pathways., Results: Our analysis demonstrates that genes in arachidonic acid metabolism (AA) are significantly perturbed across cardiomyopathies. In particular, the arachidonic acid metabolism gene PLA2G2A interacts with fibroblast marker genes and can potentially influence fibrosis during cardiomyopathy., Conclusion: The profound significance of AA metabolism within the cardiovascular system renders it a key player in modulating the phenotypes of cardiomyopathies., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (© 2023 Chauhan and Sowdhamini.)

Actin filaments help in maintaining the cell structure and coordinating cellular movements and cargo transport within the cell. Actin participates in the interaction with several proteins and also with itself to form the helical filamentous actin (F-actin). Actin-binding proteins (ABPs) and actin-associated proteins (AAPs) coordinate the actin filament assembly and processing, regulate the flux between globular G-actin and F-actin in the cell, and help maintain the cellular structure and integrity. We have used protein-protein interaction data available through multiple sources (STRING, BioGRID, mentha, and a few others), functional annotation, and classical actin-binding domains to identify actin-binding and actin-associated proteins in the human proteome. Here, we report 2482 AAPs and present an analysis of their structural and sequential domains, functions, evolutionary conservation, cellular localization, abundance, and tissue-specific expression patterns. This analysis provides a base for the characterization of proteins involved in actin dynamics and turnover in the cell., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Meenakshi S, Rao, Mayor and Sowdhamini.)

One of the key proteins that are present in the Z-disc of cardiac tissues, CSRP3, has been implicated in dilated and hypertrophic cardiomyopathy leading to heart failure. Although multiple cardiomyopathy-related mutations have been reported to reside on the two LIM domains and the disordered regions connecting the domains in this protein, the exact role of the disordered linker region is not clear. The linker harbors a few post-translational modification sites and is expected to be a regulatory site. We have carried out evolutionary studies on 5614 homologs spanning across taxa. We also performed molecular dynamics simulations of full-length CSRP3 to show that the length variations and conformational flexibility of the disordered linker could provide additional levels of functional modulation. Finally, we show that the CSRP3 homologs with widely different lengths of the linker regions could display diversity in their functional specifications. The present study provides a useful perspective to our understanding of the evolution of the disordered region between CSRP3 LIM domains., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Chauhan and Sowdhamini.)

Myotubularin or MTM1 is a lipid phosphatase that regulates vesicular trafficking in the cell. The MTM1 gene is mutated in a severe form of muscular disease, X-linked myotubular myopathy or XLMTM, affecting 1 in 50,000 newborn males worldwide. There have been several studies on the disease pathology of XLMTM, but the structural effects of missense mutations of MTM1 are underexplored due to the unavailability of a crystal structure. MTM1 consists of three domains-a lipid-binding N-terminal GRAM domain, the phosphatase domain and a coiled-coil domain which aids dimerisation of Myotubularin homologs. While most mutations reported to date map to the phosphatase domain of MTM1, the other two domains on the sequence are also frequently mutated in XLMTM. To understand the overall structural and functional effects of missense mutations on MTM1, we curated several missense mutations and performed in silico and in vitro studies. Apart from significantly impaired binding to substrate, abrogation of phosphatase activity was observed for a few mutants. Possible long-range effects of mutations from non-catalytic domains on phosphatase activity were observed as well. Coiled-coil domain mutants have been characterised here for the first time in XLMTM literature., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2023 The Authors.)

Odorant receptors (ORs) are important class of proteins involved in olfactory behaviour of insects. These are GPCR-like heptahelical transmembrane proteins with inverted topology compared to GPCR and require a co-receptor (ORco) for their function. OR function can be modulated through small molecules and negative modulation can be beneficial in case of disease vectors like Aedes aegypti . OR4 of A. aegypti is implicated in host recognition through human odour. Aedes aegypti is a vector for viruses that spread diseases like dengue, Zika and Chikungunya. In this study, we have attempted to model the full-length structure of OR4 and the ORco of A. aegypti due to lack of experimental structure. Further, we have screened a library of natural compounds (>0.3 million) along with known repellent molecules against ORco and OR4. Many natural compounds, including those from plants like Ocimum tenuiflorum (Holy Basil) and Piper nigrum (Black pepper) , were found to have better binding affinity towards ORco compared to known repellents like DEET providing an alternative to existing repellent molecules. For specific inhibitor of OR4, several natural compounds (including those from plant like Mulberry) were identified. Further, we have utilized multiple docking approaches and conservation analysis to understand the interaction between OR4 and ORco. It was observed that the residues from the seventh transmembrane helix of OR4 and pore forming helix of ORco could play an important role along with known intracellular loop 3 residues in mediating the heteromer formation of OR and ORco., Competing Interests: The authors declare no competing interests., (© 2023 The Authors.)

Non-structural protein (NS1) is a 350 amino acid long conserved protein in the dengue virus. Conservation of NS1 is expected due to its importance in dengue pathogenesis. The protein is known to exist in dimeric and hexameric states. The dimeric state is involved in its interaction with host proteins and viral replication, and the hexameric state is involved in viral invasion. In this work, we performed extensive structure and sequence analysis of NS1 protein, and uncovered the role of NS1 quaternary states in its evolution. A three-dimensional modeling of unresolved loop regions in NS1 structure is performed. "Conserved" and "Variable" regions within NS1 protein were identified from sequences obtained from patient samples and the role of compensatory mutations in selecting destabilizing mutations were identified. Molecular dynamics (MD) simulations were performed to extensively study the effect of a few mutations on NS1 structure stability and compensatory mutations. Virtual saturation mutagenesis, predicting the effect of every individual amino acid substitution on NS1 stability sequentially, revealed virtual-conserved and variable sites. The increase in number of observed and virtual-conserved regions across NS1 quaternary states suggest the role of higher order structure formation in its evolutionary conservation. Our sequence and structure analysis could enable in identifying possible protein-protein interfaces and druggable sites. Virtual screening of nearly 10,000 small molecules, including FDA-approved drugs, permitted us to recognize six drug-like molecules targeting the dimeric sites. These molecules could be promising due to their stable interactions with NS1 throughout the simulation., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Toll like receptors (TLRs) play a pivotal role in innate and adaptive immunity. There are 10 TLRs in the human genome, of which TLR10 is the least characterized. Genetic polymorphism of TLR10 has been shown to be associated with multiple diseases including tuberculosis and rheumatoid arthritis. TLR10 consists of an extracellular domain (ECD), a single-pass transmembrane (TM) helix and intracellular TIR (Toll/Interleukin-1 receptor) domain. ECD is employed for ligand recognition and the intracellular domain interacts with other TIR domain-containing adapter proteins for signal transduction. Experimental structure of ECD or TM domain is not available for TLR10. In this study, we have modelled multiple forms of TLR10-ECD dimers, such as closed and open forms, starting from available structures of homologues. Subsequently, multiple full-length TLR10 homodimer models were generated by utilizing homology modelling and protein-protein docking. The dynamics of these models in membrane-aqueous environment revealed the global motion of ECD and TIR domain towards membrane bilayer. The TIR domain residues exhibited high root mean square fluctuation compared to ECD. The 'closed form' model was observed to be energetically more favorable than 'open form' model. The evaluation of persistent interchain interactions, along with their conservation score, unveiled critical residues for each model. Further, the binding of dsRNA to TLR10 was modelled by defined and blind docking approaches. Differential binding of dsRNA to the protomers of TLR10 was observed upon simulation that could provide clues on ligand disassociation. Dynamic network analysis revealed that the 'open form' model can be the functional form while 'closed form' model can be the apo form of TLR10., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2023 The Authors.)

Background: Moringa concanensis Nimmo (MC), a plant that resembles Moringa oleifera Lam. (MO), has less scientific information but has traditionally been used as a medicinal plant. Moringa species have long been known for their medicinal qualities, which include antioxidant, anti-inflammatory, anticancer, and antihyperglycemic effects. We investigated the antidiabetic potential of MC and MO species in this study by using transcriptome profiling, metabolite analysis, and in vitro assay studies., Results: Our transcriptome analysis revealed the expression of enzymes involved in the biosynthesis of quercetin, chlorogenic acid, and benzylamine, all of which have previously been shown to have antidiabetic activity. We compared the expression patterns of five different tissues from MC and MO and it was found that the key enzymes involved in the biosynthesis of these compounds were highly expressed in leaf tissue. The expression estimated by MC transcriptome data in different tissues was verified using RT-qPCR analysis. The amount of these compounds was further quantified in the crude leaf extract of both species and found that MC had a higher abundance of quercetin and chlorogenic acid than MO. The crude leaf extract from both MC and MO were further tested in vitro, and the results demonstrated strong inhibitory activity for α-glucosidase and DPP-IV enzymes. Our findings suggest that compounds in leaf tissue, such as quercetin, benzylamine, and chlorogenic acid, could play a significant role in this antidiabetic activity. In addition, when comparing MO plants, we found that MC had a slightly higher effect in expression, abundance, and inhibitory activity., Conclusions: This study presents the first report of MC transcriptome data, as well as a comparison of its anti-diabetic activity to MO. Our analysis discussed the significance of leaf tissue in antidiabetic activity compared to other tissues of both species. Overall, this study not only provides transcriptome resources for Moringa species, but also sheds light on antidiabetic potential of both species., (© 2022. The Author(s).)

Homologous proteins can display high structural variation due to evolutionary divergence at low sequence identity. This classical inverse relationship between sequence identity and structural similarity, established many years ago, has remained true between homologous proteins of known structure over time. However, a large number of heteromeric proteins also exist in the structural data bank, where the interacting subunits belong to the same fold and maintain low sequence identity between themselves. It is not clear if there is any selection pressure to deviate from the inverse sequence-structure relationship for such interacting distant homologs, in comparison to pairs of homologs which are not known to interact. We examined 12,824 fold pairs of interacting homologs of known structure, which includes both heteromers and multi-domain proteins. These were compared with monomeric proteins, resulting in 26,082 fold pairs as a dataset of non-interacting homologous systems. Interacting homologs were found to retain higher structural similarity than non-interacting homologs at diminishing sequence identity in a statistically significant manner. Interacting homologs are more similar in their 3D structures than non-interacting homologs and have a preference towards symmetric association. There appears to be a structural constraint between remote homologs due to this commitment., (© 2022 The Authors. FEBS Open Bio published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.)

Cardiomyopathies are progressive disease conditions that give rise to an abnormal heart phenotype and are a leading cause of heart failures in the general population. These are complex diseases that show co-morbidity with other diseases. The molecular interaction network in the localised disease neighbourhood is an important step toward deciphering molecular mechanisms underlying these complex conditions. In this pursuit, we employed network medicine techniques to systematically investigate cardiomyopathy's genetic interplay with other diseases and uncover the molecular players underlying these associations. We predicted a set of candidate genes in cardiomyopathy by exploring the DIAMOnD algorithm on the human interactome. We next revealed how these candidate genes form association across different diseases and highlighted the predominant association with brain, cancer and metabolic diseases. Through integrative systems analysis of molecular pathways, heart-specific mouse knockout data and disease tissue-specific transcriptomic data, we screened and ascertained prominent candidates that show abnormal heart phenotype, including NOS3, MMP2 and SIRT1. Our computational analysis broadens the understanding of the genetic associations of cardiomyopathies with other diseases and holds great potential in cardiomyopathy research., (© 2022. The Author(s).)

RNA is the key player in many cellular processes such as signal transduction, replication, transport, cell division, transcription, and translation. These diverse functions are accomplished through interactions of RNA with proteins. However, protein-RNA interactions are still poorly derstood in contrast to protein-protein and protein-DNA interactions. This knowledge gap can be attributed to the limited availability of protein-RNA structures along with the experimental difficulties in studying these complexes. Recent progress in computational resources has expanded the number of tools available for studying protein-RNA interactions at various molecular levels. These include tools for predicting interacting residues from primary sequences, modelling of protein-RNA complexes, predicting hotspots in these complexes and insights into derstanding in the dynamics of their interactions. Each of these tools has its strengths and limitations, which makes it significant to select an optimal approach for the question of interest. Here we present a mini review of computational tools to study different aspects of protein-RNA interactions, with focus on overall application, development of the field and the future perspectives., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Bheemireddy, Sandhya, Srinivasan and Sowdhamini.)

Lantana camara L. is an invasive species of global concern. An ornamental plant originating from central America, it has now spread across natural and human-dominated habitats across tropical and subtropical regions globally. Understanding the population and evolutionary genetics of this species could help gain deeper insights into invasion biology, and provide tools for more effective management. Such investigation would require a relatively good quality genome assembly. While there have been reports of a transcriptome, it has been challenging to construct the genome assembly because of the large genome size. We present here the first draft genome assembly of Lantana camara L. which has an N50 value of 62 Kb, genome completeness of 99.3% and genome coverage of 74.3%. We hope that such an assembly will help researchers study colonization history, the genetic basis of adaptation and invasiveness, and help design strategies to contain the invasiveness of this plant, allowing biodiversity recovery in several parts of the globe., (© 2022 Biomedical Informatics.)

Aggregation of β-amyloid (Aβ) peptide is one of the hallmarks of Alzheimer's disease (AD) which results in chronic and progressive neurodegeneration of the brain. A recent study by our group have shown the ability of cyclic disulfide-rich peptides ("cyclotides") isolated from a medicinal plant, Clitoria ternatea , to inhibit the aggregation of Aβ peptides and reduce oxidative stress caused by reactive oxygen species using in vivo models of transgenic Caenorhabditis elegans . In the present study, through extensive computational docking and multi-ns molecular dynamics (MD) simulation, we evaluated if cyclotides can stably bind to Aβ molecules and/or destabilize the Aβ fibril by preventing conformational changes from α-helical to β-sheet rich structures. We demonstrate that cyclotides bind effectively and stably to different forms of Aβ structures via hydrogen bonding and hydrophobic interactions. One of the conserved hydrophobic interface residues, Tyr10 was mutated to Ala and the impact of this virtual mutation was estimated by additional MD simulations for the wild-type (WT) and mutant protein-peptide complexes. A detailed MD simulation analyses revealed that cyclotides form hydrogen bonds with the toxic amyloid assemblies thereby weakening the inter-strand hydrogen bonds between the Aβ peptide. The φ-ѱ distribution map of residues in the cyclotide binding pocket that ideally adopt β-sheet conformation show deviation towards right-handed ɑ-helical (ɑ R ) conformation. This effect was similar to that observed for the Tyr10Ala mutant and doubly so, for the cyclotide bound form. It is therefore possible to hypothesise that the opening up of amyloid β-sheet is due to an unfolding process occurring in the Aβ caused by cyclotide binding and inhibition. Our current findings provide novel structural insights on the mode of interaction between cyclotides and Aβ fibrils and describe their anti-amyloid aggregation potential. This sheds light on the future of cyclotide-based drug design against protein aggregation, a hallmark event in many neurodegenerative diseases., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Kalmankar, Gehi and Sowdhamini.)

Toll/Interleukin-1 receptor (TIR) domains are cytoplasmic domain that mediates receptor signalling. These domains are present in proteins like Toll-like receptors (TLR), its signaling adaptors and Interleukins, that form a major part of the immune system. These TIR domain containing signaling adaptors binds to the TLRs and interacts with their TIR domains for downstream signaling. We have examined the evolutionary divergence across the tree of life of two of these TIR domain containing adaptor molecules (TICAM) i.e., TIR domain-containing adapter-inducing interferon-β (TRIF/TICAM1) and TIR domain containing adaptor molecule2 (TRAM/TICAM2), by using computational approaches. We studied their orthologs, domain architecture, conserved motifs, and amino acid variations. Our study also adds a timeframe to infer the duplication of TICAM protein from Leptocardii and later divergence into TICAM1/TRIF and TICAM2/TRAM. More evidence of TRIF proteins was seen, but the absence of conserved co-existing domains such as TRIF-NTD, TIR, and RHIM domains in distant relatives hints on diversification and adaptation to different biological functions. TRAM was lost in Actinopteri and has conserved domain architecture of TIR across species except in Aves. An additional isoform of TRAM, TAG (TRAM adaptor with the GOLD domain), could be identified in species in the Mesozoic era. Finally, the Hypothesis based Likelihood ratio test was applied to look for selection pressure amongst orthologues of TRIF and TRAM to search for positively selected sites. These residues were mostly seen in the non-structural region of the proteins. Overall, this study unravels evolutionary information on the adaptors TRAM and TRIF and how well they had duplicated to perform diverse functions by changes in their domain architecture across lineages., (© 2022. The Author(s).)

Fenugreek (Trigonella foenum-graecum L.) is a self-pollinated leguminous crop belonging to the Fabaceae family. It is a multipurpose crop used as herb, spice, vegetable and forage. It is a traditional medicinal plant in India attributed with several nutritional and medicinal properties including antidiabetic and anticancer. We have performed a combined transcriptome assembly from RNA sequencing data derived from leaf, stem and root tissues. Around 209,831 transcripts were deciphered from the assembly of 92% completeness and an N50 of 1382 bases. Whilst secondary metabolites of medicinal value, such as trigonelline, diosgenin, 4-hydroxyisoleucine and quercetin, are distributed in several tissues, we report transcripts that bear sequence signatures of enzymes involved in the biosynthesis of such metabolites and are highly expressed in leaves, stem and roots. One of the antidiabetic alkaloid, trigonelline and its biosynthesising enzyme, is highly abundant in leaves. These findings are of value to nutritional and the pharmaceutical industry., (© 2022. The Author(s).)

Planarians have a remarkable ability to undergo whole-body regeneration. Successful regeneration outcome is determined by processes like polarity establishment at the wound site, which is followed by pole (organizer) specification. Interestingly, these determinants are almost exclusively expressed by muscles in these animals. However, the molecular toolkit that enables the functional versatility of planarian muscles remains poorly understood. Here we report that SMED_DDX24, a D-E-A-D Box RNA helicase, is necessary for planarian survival and regeneration. We found that DDX24 is enriched in muscles and its knockdown disrupts muscle fiber organization. This leads to defective pole specification, which in turn results in misregulation of many positional control genes specifically during regeneration. ddx24 RNAi also upregulates wound-induced Wnt signalling. Suppressing this ectopic Wnt activity rescues the knockdown phenotype by enabling better anterior pole regeneration. To summarize, our work highlights the role of an RNA helicase in muscle fiber organization, and modulating amputation-induced wnt levels, both of which seem critical for pole re-organization, thereby regulating whole-body regeneration., Competing Interests: Declaration of competing interest None declared., (Copyright © 2022. Published by Elsevier Inc.)

Background: Plant species from Rosaceae family are economically important. One of the major environmental factors impacting those species is cold stress. Although several Rosaceae plant genomes have recently been sequenced, there have been very few research conducted on cold upregulated genes and their promoter binding sites. In this study, we used computational approaches to identify and analyse potential cold stress response genes across ten Rosaceae family members., Results: Cold stress upregulated gene data from apple and strawberry were used to identify syntelogs in other Rosaceae species. Gene duplication analysis was carried out to better understand the distribution of these syntelog genes in different Rosaceae members. A total of 11,145 popular abiotic stress transcription factor-binding sites were identified in the upstream region of these potential cold-responsive genes, which were subsequently categorised into distinct transcription factor (TF) classes. MYB classes of transcription factor binding site (TFBS) were abundant, followed by bHLH, WRKY, and AP2/ERF. TFBS patterns in the promoter regions were compared among these species and gene families, found to be quite different even amongst functionally related syntelogs. A case study on important cold stress responsive transcription factor family, AP2/ERF showed less conservation in TFBS patterns in the promoter regions. This indicates that syntelogs from the same group may be comparable at the gene level but not at the level of cis-regulatory elements. Therefore, for such genes from the same family, different repertoire of TFs could be recruited for regulation and expression. Duplication events must have played a significant role in the similarity of TFBS patterns amongst few syntelogs of closely related species., Conclusions: Our study overall suggests that, despite being from the same gene family, different combinations of TFs may play a role in their regulation and expression. The findings of this study will provide information about potential genes involved in the cold stress response, which will aid future functional research of these gene families involved in many important biological processes., (© 2022. The Author(s).)

A dominant rust resistance gene, VG 9514-Rgene was isolated through map-based cloning. Sequence analysis revealed non-synonymous mutations in the TIR, NBS and LRR region of the R-protein. Candidate gene-based markers from these SNPs revealed complete co-segregation of the isolated VG 9514-Rgene with rust resistance in a RIL population and confirmed their map position in between FRS 72 and SSR_GO340445 markers in arahy03 chromosome. Blastp search of VG 9514-Rprotein detected Arahy.T6DCA5 with >80.0% identity that localized at 142,544,745.0.142,549,184 in arahy03 chromosome. K a /K s calculation revealed that VG 9514-Rgene had undergone positive selection compared to four homologous genes in the groundnut genome. Homology based structure modelling of this R-protein revealed a typical consensus three-dimensional folding of TIR-NBS-LRR protein. Non-synonymous mutations in susceptible version of R-protein were mapped and found E268Q mutation in hhGRExE motif, Y309F in RNBS-A motif and I579T in MHD motif of NB-ARC domain are probable candidates for loss of function., (Copyright © 2022 Elsevier B.V. All rights reserved.)

The olfactory system is capable of detecting and distinguishing thousands of environmental odorants that play a key role in reproduction, social behaviours including pheromones influenced classical events. Membrane secretary odorant binding proteins (OBPs) are soluble lipocalins, localized in the nasal membrane of mammals. They bind and carry odorants within the nasal epithelium to putative olfactory transmembrane receptors (ORs). OBP has not yet been exploited to develop a suitable technique to detect oestrus which is being reported as a difficult task in buffalo. In the present study, using molecular biology and protein engineering approaches, we have cloned six novel OBP isoforms from buffalo nasal epithelium odorant-binding proteins (bnOBPs). Furthermore, 3 D models were developed and molecular-docking, dynamics experiments were performed by in silico approaches. In particular, we found four residues (Phe104, Phe134, Phe69 and Asn118) in OBP1a, which contributed to favourable interactions towards two sex pheromones, specifically oleic acid and p-cresol. We expressed this protein in Escherichia coli from female buffalo urine and validated through fluorescence quenching studies to show similar strong binding affinities of OBP1a to oleic acid and p-cresol. By using structural data, the binding specificity was also verified by site-directed mutagenesis of the four residues followed by in vitro binding assays. Our results enable us to better understand the functions of different nasal epithelium OBP isoforms in buffaloes. They also lead to improved understanding of the interaction between olfactory proteins and odorants to develop highly selective biosensing devices for non-invasive detection of oestrus in buffaloes. Communicated by Ramaswamy H. Sarma.

The ongoing COVID-19 pandemic highlights the need to tackle viral variants, expand the number of antigens, and assess diverse delivery systems for vaccines against emerging viruses. In the present study, a DNA vaccine candidate was generated by combining in tandem envelope protein domain III (EDIII) of dengue virus serotypes 1-4 and a dengue virus (DENV)-2 non-structural protein 1 (NS1) protein-coding region. Each domain was designed as a serotype-specific consensus coding sequence derived from different genotypes based on the whole genome sequencing of clinical isolates in India and complemented with data from Africa. This sequence was further optimized for protein expression. In silico structural analysis of the EDIII consensus sequence revealed that epitopes are structurally conserved and immunogenic. The vaccination of mice with this construct induced pan-serotype neutralizing antibodies and antigen-specific T cell responses. Assaying intracellular interferon (IFN)-γ staining, immunoglobulin IgG2(a/c)/IgG1 ratios, and immune gene profiling suggests a strong Th1-dominant immune response. Finally, the passive transfer of immune sera protected AG129 mice challenged with a virulent, non-mouse-adapted DENV-2 strain. Our findings collectively suggest an alternative strategy for dengue vaccine design by offering a novel vaccine candidate with a possible broad-spectrum protection and a successful clinical translation either as a stand alone or in a mix and match strategy., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Sequence alignments are models that capture the structural, functional and evolutionary relationships between proteins. Structure-guided sequence alignments are helpful in the case of distantly related proteins with poor sequence identity, thus rendering routine sequence alignment methods ineffective. Protein Alignment organized as Structural Superfamilies or PASS2 database provides such sequence alignments of protein domains within a superfamily as per the Structural Classification of Proteins extended (SCOPe) database. The current update of PASS2 (i.e. PASS2.7) is following the latest release of SCOPe (2.07) and we provide data for 14 323 protein domains that are <40% identical and are organized into 2024 superfamilies. Several useful features derived from the alignments, such as conserved secondary structural motifs, HMMs and residues conserved across the superfamily, are also reported. Protein domains that are deviant from the rest of the members of a superfamily may compromise the quality of the alignment, and we found this to be the case in ∼7% of the total superfamilies we considered. To improve the alignment by objectively identifying such 'outliers', in this update, we have used a k-means-based unsupervised machine learning method for clustering superfamily members, where features provided were length of domains aligned, Cα-RMSD derived from the rigid-body superposition of all members and gaps contributed to the alignment by each domain. In a few cases, we have split the superfamily as per the clusters predicted and provided complete data for each cluster. A new feature included in this update is absolutely conserved interactions (ACIs) between residue backbones and side chains, which are obtained by aligning protein structure networks using structure-guided sequence alignments of superfamilies. ACIs provide valuable information about functionally important residues and the structure-function relationships of proteins. The ACIs and the corresponding conserved networks for backbone and sidechain have been marked on the superimposed structure separately., Database Url: The updated version of the PASS2 database is available at http://caps.ncbs.res.in/pass2/., (© The Author(s) 2022. Published by Oxford University Press.)

Cardiomyopathies are a severe and chronic cardiovascular burden worldwide, affecting a large cohort in the general population. Cysteine and glycine-rich protein 3 (CSRP3) is one of key proteins implicated in dominant dilated cardiomyopathy (DCM) and hypertrophic cardiomyopathy (HCM). In this study, we device a rapid in silico screening protocol that creates a mutational landscape map for all possible allowed and disallowed substitutions in the protein of interest. This map provides the structural and functional insights on the stability of LIM domains of CSRP3. Further, the sequence analysis delineates the eukaryotic CSRP3 protein orthologs which complements the mutational map, but provide limited information of amino acid exchanges. Next, we also evaluated the effect of HCM/DCM mutations on these domains. One of highly destabilising mutations-L44P (also disease causing) and a neutral mutation-L44M were further subjected to molecular dynamics (MD) simulations. The results establish that L44P substitution affects the LIM domain structure by altering secondary structure and due to loss of hydrophobic interaction with Phenylananine 35. The present study provides a useful perspective to our understanding of the role of mutations in the CSRP3 LIM domains and their evolution. This study provides a novel computational screening method for quick identification of key mutation sites for specific protein structures that can reduce the burden on experimental research., (© 2022. The Author(s).)

Disulphide bonds are stabilizing crosslinks in proteins and serve to enhance their thermal stability. In proteins that are small and rich in disulphide bonds, they could be the major determining factor for the choice of conformational state since their constraints on appropriate backbone conformation can be substantial. Such crosslinks and their positional conservation could itself enable protein family and functional association. Despite the importance of the field, there is no comprehensive database on disulphide crosslinks that is available to the public. Herein we provide information on disulphides in DSDBASE2.0, an updated and significantly expanded database that is freely available, fully annotated and manually curated database on native and modelled disulphides. The web interface also provides several useful computational tools that have been specifically developed for proteins containing disulphide crosslinks. The modelling of disulphide crosslinks is performed using stereochemical criteria, coded within our Modelling of Disulphides in Proteins (MODIP) algorithm. The inclusion of modelled disulphides potentially enhances the loop database substantially, thereby permitting the recognition of compatible polypeptide segments that could serve as templates for immediate modelling. The DSDBASE2.0 database has been updated to include 153,944 PDB entries, 216,096 native and 20,153,850 modelled disulphide bond segments from PDB January 2021 release. The current database also provides a resource to user-friendly search for multiple disulphide bond containing loops, along with annotation of their function using GO and subcellular localization of the query. Furthermore, it is possible to obtain the three-dimensional models of disulphide-rich small proteins using an independent algorithm, RANMOD, that generates and examines random, but allowed backbone conformations of the polypeptide. DSDBASE2.0 still remains the largest open-access repository that organizes all disulphide bonds of proteins on a single platform. The database can be accessed from http://caps.ncbs.res.in/dsdbase2., (© The Author(s) 2022. Published by Oxford University Press.)

Saffron (Crocus sativus L.) is the low yielding plant of medicinal and economic importance. Therefore, it is of interest to report the draft genome sequence of C. sativus. The draft genome of C. sativus has been assembled using Illumina sequencing and is 3.01 Gb long covering 84.24% of genome. C. sativus genome annotation identified 53,546 functional genes (including 5726 transcription factors), 862,275 repeats and 964,231 SSR markers. The genes involved in the apocarotenoids biosynthesis pathway (crocin, crocetin, picrocrocin, and safranal) were found in the draft genome analysis., Competing Interests: The authors declare no conflict of Interest., (© 2022 Biomedical Informatics.)

G-protein-coupled receptors (GPCRs) are membrane proteins which play an important role in many cellular processes and are excellent drug targets. Despite the existence of several US Food and Drug Administration (FDA)-approved GPCR-targeting drugs, there is a continuing challenge of side effects owing to the nonspecific nature of drug binding. We have investigated the diversity of the ligand binding site for this class of proteins against their cognate ligands using computational docking, even if their structures are known already in the ligand-complexed form. The cognate ligand of some of these receptors dock at allosteric binding site with better score than the binding at the conservative site. Interestingly, amino acid residues at such allosteric binding site are not conserved across GPCR subfamilies. Such a computational approach can assist in the prediction of specific allosteric binders for GPCRs., Competing Interests: Declaration of Conflicting Interests: The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article., (© The Author(s) 2021.)

The age of genomics has given us a wealth of information and the tools to study whole genomes. This, in turn, has facilitated genome-wide studies among organisms that were relatively less studied in the pre-genomic era or are non-model organisms. This paves the way to the discovery of interesting evolutionary patterns, which are brought to light by genome-wide surveys of protein superfamilies. Phosphorylation is a post-translational modification that is utilised across all clades of life, and acts as an important signalling switch, regulating several cellular processes. Tyrosine phosphatases, which are found predominantly in eukaryotes, act on phosphorylated tyrosine residues and sometimes on other substrates. Extending on our previous effort to look for tyrosine phosphatases in the human genome, we have looked for sequences of the cysteine-based tyrosine phosphatase superfamily in thirty mammalian genomes from all across Mammalia and validated the sequences with the presence of the signature catalytic motif. Domain architecture annotation, followed by in-depth analysis, revealed interesting taxon-specific patterns such as subtle differences between the protein families in marsupials and early mammals versus placental mammals. Finally, we discuss an interesting case of loss of the tyrosine phosphatase domain from a gene product in the course of eutherian evolution., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Background: Africa has one of the highest incidences of gonorrhea. Neisseria gonorrhoeae is gaining resistance to most of the available antibiotics, compromising treatment across the world. Whole-genome sequencing (WGS) is an efficient way of predicting AMR determinants and their spread in the population. Recent advances in next-generation sequencing technologies like Oxford Nanopore Technology (ONT) have helped in the generation of longer reads of DNA in a shorter duration with lower cost. Increasing accuracy of base-calling algorithms, high throughput, error-correction strategies, and ease of using the mobile sequencer MinION in remote areas lead to its adoption for routine microbial genome sequencing. To investigate whether MinION-only sequencing is sufficient for WGS and downstream analysis in resource-limited settings, we sequenced the genomes of 14 suspected N. gonorrhoeae isolates from Nairobi, Kenya., Methods: Using WGS, the isolates were confirmed to be cases of N. gonorrhoeae ( n = 9), and there were three co-occurrences of N. gonorrhoeae with Moraxella osloensis and N. meningitidis ( n = 2). N. meningitidis has been implicated in sexually transmitted infections in recent years. The near-complete N. gonorrhoeae genomes ( n = 10) were analyzed further for mutations/factors causing AMR using an in-house database of mutations curated from the literature., Results: We observe that ciprofloxacin resistance is associated with multiple mutations in both gyrA and parC. Mutations conferring tetracycline ( rpsJ ) and sulfonamide ( folP ) resistance and plasmids encoding beta-lactamase were seen in all the strains, and tet(M) -containing plasmids were identified in nine strains. Phylogenetic analysis clustered the 10 isolates into clades containing previously sequenced genomes from Kenya and countries across the world. Based on homology modeling of AMR targets, we see that the mutations in GyrA and ParC disrupt the hydrogen bonding with quinolone drugs and mutations in FolP may affect interaction with the antibiotic., Conclusion: Here, we demonstrate the utility of mobile DNA sequencing technology in producing a consensus genome for sequence typing and detection of genetic determinants of AMR. The workflow followed in the study, including AMR mutation dataset creation and the genome identification, assembly, and analysis, can be used for any clinical isolate. Further studies are required to determine the utility of real-time sequencing in outbreak investigations, diagnosis, and management of infections, especially in resource-limited settings., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Juma, Sankaradoss, Ndombi, Mwaura, Damodar, Nazir, Pandit, Khurana, Masika, Chirchir, Gachie, Krishna, Sowdhamini, Anzala and Meenakshi.)

Protein-ligand binding prediction has extensive biological significance. Binding affinity helps in understanding the degree of protein-ligand interactions and is a useful measure in drug design. Protein-ligand docking using virtual screening and molecular dynamic simulations are required to predict the binding affinity of a ligand to its cognate receptor. Performing such analyses to cover the entire chemical space of small molecules requires intense computational power. Recent developments using deep learning have enabled us to make sense of massive amounts of complex data sets where the ability of the model to "learn" intrinsic patterns in a complex plane of data is the strength of the approach. Here, we have incorporated convolutional neural networks to find spatial relationships among data to help us predict affinity of binding of proteins in whole superfamilies toward a diverse set of ligands without the need of a docked pose or complex as user input. The models were trained and validated using a stringent methodology for feature extraction. Our model performs better in comparison to some existing methods used widely and is suitable for predictions on high-resolution protein crystal (⩽2.5 Å) and nonpeptide ligand as individual inputs. Our approach to network construction and training on protein-ligand data set prepared in-house has yielded significant insights. We have also tested DEELIG on few COVID-19 main protease-inhibitor complexes relevant to the current public health scenario. DEELIG-based predictions can be incorporated in existing databases including RSCB PDB, PDBMoad, and PDBbind in filling missing binding affinity data for protein-ligand complexes., Competing Interests: Declaration of conflicting interests: The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article., (© The Author(s) 2021.)

Neurotoxic aggregation of β-amyloid (Aβ) peptides is a hallmark of Alzheimer's disease and increased reactive oxygen species (ROS) is an associated process. In the present study, we report the neuroprotective effects of disulfide-rich, circular peptides from Clitoria ternatea ( C. ternatea ) (butterfly pea) on Aβ-induced toxicity in transgenic Caenorhabditis elegans . Cyclotides (∼30 amino acids long) are a special class of cyclic cysteine knot peptides. We show that cyclotide-rich fractions from different plant tissues delay Aβ-induced paralysis in the transgenic CL4176 strain expressing the human muscle-specific Aβ 1-42 gene. They also improved Aβ-induced chemotaxis defects in CL2355 strain expressing Aβ 1-42 in the neuronal cells. ROS assay suggests that this protection is likely mediated by the inhibition of Aβ oligomerization. Furthermore, Aβ deposits were reduced in the CL2006 strain treated with the fractions. The study shows that cyclotides from C. ternatea could be a source of a novel pharmacophore scaffold against neurodegenerative diseases.

Understanding exposures to low doses of ionizing radiation are relevant since most environmental, diagnostic radiology and occupational exposures lie in this region. However, the molecular mechanisms that drive cellular responses at these doses, and the subsequent health outcomes, remain unclear. A local monazite-rich high level natural radiation area (HLNRA) in the state of Kerala on the south-west coast of Indian subcontinent show radiation doses extending from ≤ 1 to ≥ 45 mGy/y and thus, serve as a model resource to understand low dose mechanisms directly on healthy humans. We performed quantitative discovery proteomics based on multiplexed isobaric tags (iTRAQ) coupled with LC-MS/MS on human peripheral blood mononuclear cells from HLNRA individuals. Several proteins involved in diverse biological processes such as DNA repair, RNA processing, chromatin modifications and cytoskeletal organization showed distinct expression in HLNRA individuals, suggestive of both recovery and adaptation to low dose radiation. In protein-protein interaction (PPI) networks, YWHAZ (14-3-3ζ) emerged as the top-most hub protein that may direct phosphorylation driven pro-survival cellular processes against radiation stress. PPI networks also identified an integral role for the cytoskeletal protein ACTB, signaling protein PRKACA; and the molecular chaperone HSPA8. The data will allow better integration of radiation biology and epidemiology for risk assessment [Data are available via ProteomeXchange with identifier PXD022380].

This protocol describes a stepwise process to identify proteins of interest from a query proteome derived from NGS data. We implemented this protocol on Moringa oleifera transcriptome to identify proteins involved in secondary metabolite and vitamin biosynthesis and ion transport. This knowledge-driven protocol identifies proteins using an integrated approach involving sensitive sequence search and evolutionary relationships. We make use of functionally important residues (FIR) specific for the query protein family identified through its homologous sequences and literature. We screen protein hits based on the clustering with true homologues through phylogenetic tree reconstruction complemented with the FIR mapping. The protocol was validated for the protein hits through qRT-PCR and transcriptome quantification. Our protocol demonstrated a higher specificity as compared to other methods, particularly in distinguishing cross-family hits. This protocol was effective in transcriptome data analysis of M. oleifera as described in Pasha et al.•Knowledge-driven protocol to identify secondary metabolite synthesizing protein in a highly specific manner.•Use of functionally important residues for screening of true hits.•Beneficial for metabolite pathway reconstruction in any (species, metagenomics) NGS data., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2020 The Authors. Published by Elsevier B.V.)

N-acetylglucosamine 6-phosphate deacetylase (NagA) catalyzes the conversion of N-acetylglucosamine-6-phosphate to glucosamine-6-phosphate in amino sugar catabolism. This conversion is an essential step in the catabolism of sialic acid in several pathogenic bacteria, including Pasteurella multocida, and thus NagA is identified as a potential drug target. Here, we report the unique structural features of NagA from P. multocida (PmNagA) resolved to 1.95 Å. PmNagA displays an altered quaternary architecture with unique interface interactions compared to its close homolog, the Escherichia coli NagA (EcNagA). We confirmed that the altered quaternary structure is not a crystallographic artifact using single particle electron cryo-microscopy. Analysis of the determined crystal structure reveals a set of hot-spot residues involved in novel interactions at the dimer-dimer interface. PmNagA binds to one Zn 2+ ion in the active site and demonstrates kinetic parameters comparable to other bacterial homologs. Kinetic studies reveal that at high substrate concentrations (~10-fold the K M ), the tetrameric PmNagA displays hysteresis similar to its distant neighbor, the dimeric Staphylococcus aureus NagA (SaNagA). Our findings provide key information on structural and functional properties of NagA in P. multocida that could be utilized to design novel antibacterials., (© 2020 Wiley Periodicals LLC.)