Search

Nine viruses of the bean common mosaic virus subgroup of potyviruses are major international crop pathogens, but their phylogenetically closest relatives from non-crop plants have mostly been found only in SE Asia and Oceania, which is thus likely to be their "centre of emergence". We have compared over 700 of the complete genomic ORFs of the crop pandemic and the non-crop viruses in various ways. Only one-third of crop virus genomes are non-recombinant, but more than half the non-crop virus genomes are. Four of the viruses were from crops domesticated in the Old World (Africa to SE Asia), and the other five were from New World crops. There was a temporal signal in only three of the crop virus datasets, but it confirmed that the most recent common ancestors of all the crop viruses were before inter-continental marine trade started after 1492 CE, whereas all the crown clusters of the phylogenies are from after that date. The non-crop virus datasets are genetically more diverse than those of the crop viruses, and Tajima's D analyses showed that their populations were contracting, and only one of the crop viruses had a significantly expanding population. dN/dS analyses showed that most of the genes and codons in all the viruses were under significant negative selection, and the few that were under significant positive selection were mostly in the PIPO-encoding region of the P3 protein, or the PIPO protein itself. Interestingly, more positively selected codons were found in non-crop than in crop viruses, and, as the hosts of the former were taxonomically more diverse than the latter, this may indicate that the positively selected codons are involved in host range determination; AlphaFold3 modelling was used to investigate this possibility.

Eotaxin-3 is a key chemokine with a relevant role in eosinophilic esophagitis, a rare chronic immune/antigen-mediated inflammatory disorder. Eotaxin-3 is a potent activator of eosinophil emergence and migration, which may lead to allergic airway inflammation. We investigated, using bioinformatics tools, the protein structure and the possible effects of the known variations reported in public databases. Following a procedure already established, we created a 3D model of the whole protein and modeled the structure of 105 protein variants due to known point mutations. The effects of the amino acid substitution at the level of impact on protein structure, stability, and possibly function were detected by the bioinformatics procedure and described in detail. A web application was implemented to browse the results of the analysis and visualize the 3D models, with the opportunity of downloading the models and analyzing them using their own software. Among 105 amino acid substitutions investigated, the study evidenced in 44 cases at least one change in any of the investigated structural parameters. Other six variations are also relevant, although a structural effect was not detected by our analysis, because they affected amino acids highly conserved, which suggests a possible function role. All these variations should be the object of particular attention, as they may induce a loss of functionality in the protein.

The origin of flight and laryngeal echolocation in bats is likely to have been accompanied by evolutionary changes in other aspects of their sensory biology. Of all sensory modalities in bats, olfaction is perhaps the least well understood. Olfactory receptors (ORs) function in recognizing odour molecules, with crucial roles in evaluating food, as well as in processing social information. Here we compare OR repertoire sizes across taxa and apply a new pipeline that integrates comparative genome data with protein structure modelling and then we employ molecular docking techniques with small molecules to analyse OR functionality based on binding energies. Our results suggest a sharp contraction in odorant recognition of the functional OR repertoire during the origin of bats, consistent with a reduced dependence on olfaction. We also compared bat lineages with contrasting different ecological characteristics and found evidence of differences in OR gene expansion and contraction, and in the composition of ORs with different tuning breadths. The strongest binding energies of ORs in non-echolocating fruit-eating bats were seen to correspond to ester odorants, although we did not detect a quantitative advantage of functional OR repertoires in these bats compared with echolocating insectivorous species. Overall, our findings based on molecular modelling and computational docking suggest that bats have undergone olfactory evolution linked to dietary adaptation. Our results from extant and ancestral bats help to lay the groundwork for targeted experimental functional tests in the future., (© 2024 John Wiley & Sons Ltd.)

Unlike animals, variability in transcription factors (TFs) and their binding regions (TFBRs) across the plants species is a major problem that most of the existing TFBR finding software fail to tackle, rendering them hardly of any use. This limitation has resulted into underdevelopment of plant regulatory research and rampant use of Arabidopsis-like model species, generating misleading results. Here, we report a revolutionary transformers-based deep-learning approach, PTFSpot, which learns from TF structures and their binding regions' co-variability to bring a universal TF-DNA interaction model to detect TFBR with complete freedom from TF and species-specific models' limitations. During a series of extensive benchmarking studies over multiple experimentally validated data, it not only outperformed the existing software by >30% lead but also delivered consistently >90% accuracy even for those species and TF families that were never encountered during the model-building process. PTFSpot makes it possible now to accurately annotate TFBRs across any plant genome even in the total lack of any TF information, completely free from the bottlenecks of species and TF-specific models., (© The Author(s) 2024. Published by Oxford University Press.)

Microscopic, photosynthetic prokaryotes and eukaryotes, collectively referred to as microalgae, are widely studied to improve our understanding of key metabolic pathways (e.g., photosynthesis) and for the development of biotechnological applications. Omics technologies, which are now common tools in biological research, have been shown to be critical in microalgal research. In the past decade, significant technological advancements have allowed omics technologies to become more affordable and efficient, with huge datasets being generated. In particular, where studies focused on a single or few proteins decades ago, it is now possible to study the whole proteome of a microalgae. The development of mass spectrometry-based methods has provided this leap forward with the high-throughput identification and quantification of proteins. This review specifically provides an overview of the use of proteomics in fundamental (e.g., photosynthesis) and applied (e.g., lipid production for biofuel) microalgal research, and presents future research directions in this field.

Lipoxygenases make several biological functions in cells, based on the products of the catalyzed reactions. In diatoms, microalgae ubiquitous in aquatic ecosystems, lipoxygenases have been noted for the oxygenation of fatty acids with the production of oxylipins, which are involved in many physiological and pathological processes in marine organisms. The interest in diatoms' lipoxygenases and oxylipins has increased due to their possible biotechnological applications, ranging from ecology to medicine. We investigated using bioinformatics and molecular docking tools the lipoxygenases of diatoms and the possible interaction with substrates. A large-scale analysis of sequence resources allowed us to retrieve 45 sequences of lipoxygenases from diatoms. We compared and analyzed the sequences by multiple alignments and phylogenetic trees, suggesting the possible clustering in phylogenetic groups. Then, we modelled the 3D structure of representative enzymes from the different groups and investigated in detail the structural and functional properties by docking simulations with possible substrates. The results allowed us to propose a classification of the lipoxygenases from diatoms based on their sequence features, which may be reflected in specific structural differences and possible substrate specificity.

The majority of picornaviral 3C proteinases (3C pro ) cleavage sites possess glutamine at the P1 position. Plant nepovirus 3C-like proteinases (3CL pro ) show however much broader specificity, cleaving not only after glutamine, but also after several basic and hydrophobic residues. To investigate this difference, we employed AlphaFold to generate structural models of twelve selected 3CL pro , representing six substrate specificities. Generally, we observed favorable correlations between the architecture and charge of nepovirus proteinase S1 subsites and their ability to accept or restrict larger residues. The models identified a conserved aspartate residue close to the P1 residue in the S1 subsites of all nepovirus proteinases examined, consistent with the observed strong bias against negatively-charged residues at the P1 position of nepovirus cleavage sites. Finally, a cramped S4 subsite along with the presence of two unique histidine and serine residues explains the strict requirement of the grapevine fanleaf virus proteinase for serine at the P4 position., 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 The Authors. Published by Elsevier Inc. All rights reserved.)

Introduction: The sirtuin (Silent mating type information regulation 2 homolog)1(SIRT1) protein plays a vital role in many disorders such as diabetes, cancer, obesity, inflammation, and neurodegenerative and cardiovascular diseases. The objective of this in silico analysis of SIRT1's functional single nucleotide polymorphisms (SNPs) was to gain valuable insight into the harmful effects of non-synonymous SNPs (nsSNPs) on the protein. The objective of the study was to use bioinformatics methods to investigate the genetic variations and modifications that may have an impact on the SIRT1 gene's expression and function., Methods: nsSNPs of SIRT1 protein were collected from the dbSNP site, from its three (3) different protein accession IDs. These were then fed to various bioinformatic tools such as SIFT, Provean, and I- Mutant to find the most deleterious ones. Functional and structural effects were examined using the HOPE server and I-Tasser. Gene interactions were predicted by STRING software. The SIFT, Provean, and I-Mutant tools detected the most deleterious three nsSNPs (rs769519031, rs778184510, and rs199983221)., Results: Out of 252 nsSNPs, SIFT analysis showed that 94 were deleterious, Provean listed 67 dangerous, and I-Mutant found 58 nsSNPs resulting in lowered stability of proteins. HOPE modelling of rs199983221 and rs769519031 suggested reduced hydrophobicity due to Ile 4Thr and Ile223Ser resulting in decreased hydrophobic interactions. In contrast, on modelling rs778184510, the mutant protein had a higher hydrophobicity than the wild type., Conclusions: Our study reports that three nsSNPs (D357A, I223S, I4T) are the most damaging mutations of the SIRT1 gene. Mutations may result in altered protein structure and functions. Such altered protein may be the basis for various disorders. Our findings may be a crucial guide in establishing the pathogenesis of various disorders., Competing Interests: No competing interests were disclosed., (Copyright: © 2024 Thayyil Menambath D et al.)

The outer membrane lipase (oml) gene, encoding a novel autotransporter-dependent lipase from Pseudomonas guariconensis, was cloned and sequenced. The oml gene has an open reading frame of 1866 bp. It encodes the 621 amino acid autotransporter-dependent GDSL lipase (OML), which has the highest sequence similarity (64.08 %) with the EstA of Pseudomonas aeruginosa (PDB:3kvn.1. A). OML was expressed and purified, which showed a purified band of approximately 70 kDa. The purified enzyme showed maximum activity at pH 9 and 40 °C. Substrate specificity studies and kinetic study by Lineweaver-Burk plot of purified OML showed K m of 1.27 mM and Vmax of 333.33 U/mL with p-nitrophenyl palmitate. The purified enzyme showed good stability in the presence of hexane, methanol, and ethanol, while the presence of the metal ion Mg 2+ showed maximum lipase activity. Bioinformatics analysis supported the in vitro findings by predicting enzyme substrate specificity towards long-chain fatty acids and fatty acids with shorter chain lengths. The stability of the interaction of the protein-ligand complex (OML-ricinoleic acid) was confirmed using MDS and castor oil bioconversion using purified OML was confirmed using High-Performance Liquid Chromatography (HPLC)., Competing Interests: Declaration of competing interest The authors declare no conflict of interest in publishing this article in this journal., (Copyright © 2023 Elsevier B.V. All rights reserved.)

The results of tertiary structure assessment at CASP15 are reported. For the first time, recognizing the outstanding performance of AlphaFold 2 (AF2) at CASP14, all single-chain predictions were assessed together, irrespective of whether a template was available. At CASP15, there was no single stand-out group, with most of the best-scoring groups-led by PEZYFoldings, UM-TBM, and Yang Server-employing AF2 in one way or another. Many top groups paid special attention to generating deep Multiple Sequence Alignments (MSAs) and testing variant MSAs, thereby allowing them to successfully address some of the hardest targets. Such difficult targets, as well as lacking templates, were typically proteins with few homologues. Local divergence between prediction and target correlated with localization at crystal lattice or chain interfaces, and with regions exhibiting high B-factor factors in crystal structure targets, and should not necessarily be considered as representing error in the prediction. However, analysis of exposed and buried side chain accuracy showed room for improvement even in the latter. Nevertheless, a majority of groups produced high-quality predictions for most targets, which are valuable for experimental structure determination, functional analysis, and many other tasks across biology. These include those applying methods similar to those used to generate major resources such as the AlphaFold Protein Structure Database and the ESM Metagenomic atlas: the confidence estimates of the former were also notably accurate., (© 2023 The Authors. Proteins: Structure, Function, and Bioinformatics published by Wiley Periodicals LLC.)

The affinity of specific phenolic compounds (PCs) and capsaicinoids (CAPs) present in three Capsicum annuum varieties (Friariello, Cayenne and Dzuljunska Sipka) to the transient receptor potential vanilloid member 1 (TRPV1) was investigated by integrating an analytic approach for the simultaneous extraction and analysis through high-performance liquid chromatography coupled with ion trap mass spectrometry (HPLC/ITMS) and UV detection (HPLC-UV) of PCs and CAPs and structural bioinformatics based on the protein modelling and molecular simulations of protein-ligand docking. Overall, a total of 35 compounds were identified in the different samples and CAPs were quantified. The highest content of total polyphenols was recorded in the pungent Dzuljunska Sipka variety (8.91 ± 0.05 gGAE/Kg DW) while the lowest was found in the non-pungent variety Friariello (3.58 ± 0.02 gGAE/Kg DW). Protein modelling generated for the first time a complete model of the homotetrameric human TRPV1, and it was used for docking simulations with the compounds detected via the analytic approach, as well as with other compounds, as an inhibitor reference. The simulations indicate that different capsaicinoids can interact with the receptor, providing details on the molecular interaction, with similar predicted binding energy values. These results offer new insights into the interaction of capsaicinoids with TRPV1 and their possible actions.

Protein bioinformatics has been applied to a myriad of opportunities in biocatalysis from enzyme engineering to enzyme discovery, but its application in enzyme immobilization is still very limited. Enzyme immobilization brings clear advantages in terms of sustainability and cost-efficiency but is still limited in its implementation. This, because it is a technique that remains tied to a quasi-blind protocol of trial and error, and therefore, is regarded as a time-intensive and costly approach. Here, we present the use of a set of bioinformatic tools to rationalize the results of protein immobilization that have been previously described. The study of proteins with these new tools allows the discovery of key driving forces in the process of immobilization that explain the obtained results, moving us a step closer to the final goal: predictive enzyme immobilization protocols., (© 2023 The Authors. ChemBioChem published by Wiley-VCH GmbH.)

We have used an integrated computational approach to explore the role of vitamin C and vitamin D in preventing aggregation of Fibrinogen A alpha-chain (FGA cter ) protein responsible for renal amyloidosis. We modelled structures of E524K / E526K mutants of FGA cter protein and examined the potential interactions of these mutants with vitamin C and vitamin D3. Interaction of these vitamins at the amyloidogenic site may prevent the intermolecular interaction required for amyloid formation. The binding free energy values of vitamin C and vitamin D3 for E524K FGA cter and E526K FGA cter are - 67.12 ± 30.46 kJ/mole and - 79.45 ± 26.12 kJ/mol, respectively. Experimental studies using Congo red absorption, aggregation index studies and AFM imaging show encouraging results. The AFM images of E526K FGA cter contained more extensive and higher protofibril aggregates, whereas, in the presence of vitamin D3, small monomeric and oligomeric aggregates were observed. Overall, the works provide interesting results about vitamin C and D role in preventing renal amyloidosis., Competing Interests: Declaration of Interest Statement The authors declared no conflict of interest., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Objective: Idiopathic hypogonadotropic hypogonadism (IHH) is classified into two groups-Kalman syndrome and normosmic IHH (nIHH). Half of all cases can be explained by mutations in >50 genes. Targeted gene panel testing with nexrt generation sequencing (NGS) is required for patients without typical phenotypic findings. The aim was to determine the genetic etiologies of patients with IHH using NGS, including 54 IHH-associated genes, and to present protein homology modeling and protein stability analyzes of the detected variations., Methods: Clinical and demographic data of 16 patients (eight female), aged between 11.6-17.8 years, from different families were assessed. All patients were followed up for a diagnosis of nIHH, had normal cranial imaging, were without anterior pituitary hormone deficiency other than gonadotropins, had no sex chromosome anomaly, had no additional disease, and underwent genetic analysis with NGS between the years 2008-2021. Rare variants were classified according to the variant interpretation framework of the American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology. Changes in protein structure caused by variations were modeled using RoseTTAFold and changes in protein stability resulting from variation were analyzed., Results: Half of the 16 had no detectable variation. Three (18.75%) had a homozygous (pathogenic) variant in the GNRHR gene, one (6.25%) had a compound heterozygous [likely pathogenic-variants of uncertain significance (VUS)] variant in PROK2 and four (25%) each had a heterozygous (VUS) variant in HESX1, FGF8, FLRT3 and DMXL2 . Protein models showed that variants interpreted as VUS according to ACMG could account for the clinical IHH., Conclusion: The frequency of variation detection was similar to the literature. Modelling showed that the variant in five different genes, interpreted as VUS according to ACMG, could explain the clinical IHH., (©Copyright 2023 by Turkish Society for Pediatric Endocrinology and Diabetes | The Journal of Clinical Research in Pediatric Endocrinology published by Galenos Publishing House.)

Meiotic recombination 11 ( MRE11 ) is a component of the tri-molecular MRE11 -RAD50-NBS1 (MRN) complex, which functions as an exonuclease and endonuclease which is involved in identifying, signalling, protecting and repairing double-strand breaks in DNA (DSBs). Ataxia-telangiectasia-like disorder (ATLD) 1 and Nijmegen breakage syndrome (NBS)-like disorder are MRE11 associated diseases. In the present study, we used an integrated computational approach to identify the most deleterious SNPs and their structural and functional impact on human MRE11 . Five of the 68 observed non-synonymous SNP (nsSNPs; I162T, S273C, W210C, D311Y and R364L) should be worked on due to their strong possible pathogenicity and the risk of changing protein properties. All the nsSNPs were highly conserved and decrease the protein stability located in the MRE11 nuclease and MRE11 DNA binding presumed domain. R364L and I162T were predicted to be involved in post-translational modification (PTM) sites. Furthermore, we also analysed the regulatory effect of noncoding SNPs on MRE11 gene regulation in which 6 SNPs were found to affect gene regulation. All six noncoding SNPs predicted chromatin interactive site whereas only one SNP was noted its association with miRNA binding site which disrupts 5 miRNA conserved site. These findings help future studies to get more insights into the role of these variants in the alteration of the MRE11 function.Communicated by Ramaswamy H. Sarma.

The refinement of predicted 3D models aims to bring them closer to the native structure by fixing errors including unusual bonds and torsion angles and irregular hydrogen bonding patterns. Refinement approaches based on molecular dynamics (MD) simulations using different types of restraints have performed well since CASP10. ReFOLD, developed by the McGuffin group, was one of the many MD-based refinement approaches, which were tested in CASP 12. When the performance of the ReFOLD method in CASP12 was evaluated, it was observed that ReFOLD suffered from the absence of a reliable guidance mechanism to reach consistent improvement for the quality of predicted 3D models, particularly in the case of template-based modelling (TBM) targets. Therefore, here we propose to utilize the local quality assessment score produced by ModFOLD6 to guide the MD-based refinement approach to further increase the accuracy of the predicted 3D models. The relative performance of the new local quality assessment guided MD-based refinement protocol and the original MD-based protocol ReFOLD are compared utilizing many different official scoring methods. By using the per-residue accuracy (or local quality) score to guide the refinement process, we are able to prevent the refined models from undesired structural deviations, thereby leading to more consistent improvements. This chapter will include a detailed analysis of the performance of the local quality assessment guided MD-based protocol versus that deployed in the original ReFOLD method., (© 2023. Springer Science+Business Media, LLC, part of Springer Nature.)

Barth syndrome is a rare disease affecting mitochondria structure and function in males. In our previous study, we have shown a new mutation (c.83T>A, p.Val28Glu) in the TAZ gene in two affected patients with congenital cardiomyopathy. Furthermore, women in this family had no mutations in their blood cells, whereas they only had mutations in the oral epithelial cells. The objective of the project was to evaluate the effect of intertissue mosaicisms on the Barth syndrome phenotypes, searching for another disease-related loci on chromosome X and finally to assess the consequences of the mutation. We conducted the advanced genetic study including cytogenetic research (constitutional karyotyping in blood and fibroblasts), NGS sequencing (with custom chromosome X sequencing together with the evaluation of loss of heterozygosity (LOH) and aberrations (CNV) in the whole genome) in four different tissues and sequencing of tafazzin and deoxyribonuclease 1 like 1 transcripts. The presence of deletions within the 5'untranslated region of the TAZ gene and/or the noncoding regions of the DNASE1L1 gene were detected in several tissues. Whereas, there is no intertissue mosaicism regarding point mutation in TAZ gene in all investigated tissues in female carriers. Only the male patient presented biochemical markers and neurological symptoms of Barth syndrome. All the female carriers are healthy and have normal tafazzin and deoxyribonuclease 1 like 1 transcripts in 2 analyzed tissues. The conclusion of this study is that we cannot rule out or confirm mosaicism in the noncoding regions of TAZ or DNASE1L1 genes, but this is not clinically relevant in female carriers because they are healthy. Finally, it has been proven that mutation (c.83T>A, p.Val28Glu) is responsible for disease in males in this family., Competing Interests: Conflicts of interest The authors declare no conflict of interest., (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

Aquaporins are an ancient family of membrane channel proteins present in all eukaryotes and most prokaryotes, and apart from water, allow the transport of neutral solutes and organic compounds through the pore. These proteins are essential role differentially expressed during ripening in Fragaria vesca fruits. Fv PIP2-1a is intensively expressed in fruits, inclusive several other proteins member are differentially expressed in fruit but also in other plant tissues. Phylogenetic analysis shows that Fv PIP2-1a grouped with other Fragaria proteins and far apart from other F. vesca PIP proteins. A structural model for Fv PIP2-1a protein was built by comparative modelling methodology, which was validated and refined by molecular dynamics simulation. Fv PIP2-1a structure consists of 6 transmembrane regions and two NPA domains. The mobilization of water was analyzed by molecular docking simulations in wildtype and two mutants. Interestingly, the mutant FvPIP2-1a_H214G allowed the prediction of an increment in the flux of water molecules. On contrary, structural analysis predicted that H214E mutation blocked passage of water associated to constriction of the pore., 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. Raul Herrera reports financial support was provided by ANID- Fondecyt and ANID-Anillo., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Since its discovery, bovine theileriosis has caused major socioeconomic losses in sub-Saharan Africa. Acaricide resistance of the intermediate host, paucity of therapeutics, and lack of sufficiently cross-protective vaccines increase the risk of parasite spread due to global warming. Here, we highlight three important areas that require investigation to develop next-generation vaccines., Competing Interests: Declaration of interests The authors have no competing interests to declare., (Copyright © 2022 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Visual opsins of vertebrates and invertebrates diversified independently and converged to detect ultraviolet to long wavelengths (LW) of green or red light. In both groups, colour vision largely derives from opsin number, expression patterns and changes in amino acids interacting with the chromophore. Functional insights regarding invertebrate opsin evolution have lagged behind those for vertebrates because of the disparity in genomic resources and the lack of robust in vitro systems to characterize spectral sensitivities. Here, we review bioinformatic approaches to identify and model functional variation in opsins as well as recently developed assays to measure spectral phenotypes. In particular, we discuss how transgenic lines, cAMP-spectroscopy and sensitive heterologous expression platforms are starting to decouple genotype-phenotype relationships of LW opsins to complement the classical physiological-behavioural-phylogenetic toolbox of invertebrate visual sensory studies. We illustrate the use of one heterologous method by characterizing novel LW Gq opsins from 10 species, including diurnal and nocturnal Lepidoptera, a terrestrial dragonfly and an aquatic crustacean, expressing them in HEK293T cells, and showing that their maximum absorbance spectra ( λ max ) range from 518 to 611 nm. We discuss the advantages of molecular approaches for arthropods with complications such as restricted availability, lateral filters, specialized photochemistry and/or electrophysiological constraints. This article is part of the theme issue 'Understanding colour vision: molecular, physiological, neuronal and behavioural studies in arthropods'.

Isoflavonoid synthase (IFS) is a critical enzyme for the biosynthesis of over 2400 isoflavonoids. Isoflavonoids are an important class of plant secondary metabolites that have a range of pharmaceutical and nutraceutical properties. With growing interest in isoflavonoids from both research and industrial perspectives, efforts are being forwarded to enhance isoflavonoid production in-planta and ex-planta; therefore, in-silico analysis and characterisation of available IFS protein sequences are needed. The present study is the first-ever attempt toward phylogenetic analysis and protein modelling of available IFS protein sequences. Phylogenetic analysis has shown that IFS amino acid sequences have 86.4% pairwise identity and 26.5% identical sites, and the sequences were grouped into six different clades. The presence of a β-hairpin and extra loop at catalytic sites of Trifolium pratense , Beta vulgaris and Medicago truncatula , respectively, compared with Glycyrrhiza echinata are critical structural differences that may affect catalytic function. Protein docking highlighted the preference of selected IFS for liquiritigenin compared with naringenin and has listed T. pratense as the most efficient candidate for heterologous biosynthesis of isoflavonoids. The in-silico characterisation of IFS represented in this study is vital in realising the new bioengineering endeavours and will help in the characterisation and selection of IFS candidate enzymes for heterologous biosynthesis of isoflavonoids.

Main Conclusion: Genetic variation in seed protein composition, seed protein gene expression and predictions of seed protein physiochemical properties were documented in C. sativa and other Camelina species. Seed protein diversity was examined in six Camelina species (C. hispida, C. laxa, C. microcarpa, C. neglecta, C. rumelica and C. sativa). Differences were observed in seed protein electrophoretic profiles, total seed protein content and amino acid composition between the species. Genes encoding major seed proteins (cruciferins, napins, oleosins and vicilins) were catalogued for C. sativa and RNA-Seq analysis established the expression patterns of these and other genes in developing seed from anthesis through to maturation. Examination of 187 C. sativa accessions revealed limited variation in seed protein electrophoretic profiles, though sufficient to group the majority into classes based on high MW protein profiles corresponding to the cruciferin region. C. sativa possessed four distinct types of cruciferins, named CsCRA, CsCRB, CsCRC and CsCRD, which corresponded to orthologues in Arabidopsis thaliana with members of each type encoded by homeologous genes on the three C. sativa sub-genomes. Total protein content and amino acid composition varied only slightly; however, RNA-Seq analysis revealed that CsCRA and CsCRB genes contributed > 95% of the cruciferin transcripts in most lines, whereas CsCRC genes were the most highly expressed cruciferin genes in others, including the type cultivar DH55. This was confirmed by proteomics analyses. Cruciferin is the most abundant seed protein and contributes the most to functionality. Modelling of the C. sativa cruciferins indicated that each type possesses different physiochemical attributes that were predicted to impart unique functional properties. As such, opportunities exist to create C. sativa cultivars with seed protein profiles tailored to specific technical applications., (© 2022. Crown.)

Multi-drug resistant (MDR) Shigella flexneri 2a, one of the leading bacterial agents of diarrhoeal mortality, has posed challenges in treatment strategies. The present study was conducted to identify potential therapeutic biomarkers using gene interaction network (GIN) in order to understand the cellular and molecular level interactions of both antimicrobial resistance (AMR) and virulence genes through topological and clustering metrics. Statistically significant differential gene expression (DGE), structural chemistry and dynamics were incorporated to elucidate biomarker for sustainable therapeutic regimen against MDR S. flexneri . Functional enrichments and topological metrics revealed evgS, ybjZ, tolC, gyrA, parC and their direct interactors to be associated with diverse AMR mechanisms. Histidine kinase EvgS was considered as the hub protein due to its highest prevalence in the molecular interactome profiles of both the AMR (71.6%) and virulence (45.8%) clusters interconnecting several genes concerning two-component system (TCS). DGE profiles of ΔPhoPQ (deleted regulatory PhoP and sensor PhoQ) led to the upregulation of TCS comprising EvgSA thereby validating EvgS as a promising therapeutic biomarker. Druggability and structural stability of EvgS was assessed through thermal shifts, backbone stability and coarse dynamics refinement. Structure-function relationship was established revealing the C-terminal extracellular domain as the drug-binding site which was further validated through molecular dynamics simulation. Structure elucidation of identified biomarker followed by secondary and tertiary structural validation would prove pivotal for future therapeutic interventions against subverting both AMR and virulence posed by this strain., Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-022-03325-w., Competing Interests: Conflict of interestThe authors have no competing financial or non-financial interests to declare that are relevant to the content of this article., (© King Abdulaziz City for Science and Technology 2022, Springer Nature or its licensor holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.)

Aluminum (Al) toxicity poses a significant challenge for the yield improvement of chickpea, which is an economically important legume crop with high nutritional value in human diets. The genetic basis of Al-tolerance in chickpea remains unclear. Here, we assessed the Al-tolerance of 8 wild Cicer and one cultivated chickpea (PBA Pistol) accessions by measuring the root elongation in solution culture under control (0 μM Al 3+ ) and Al treatments (15, 30 μM Al 3+ ). Compared to PBA Pistol, the wild Cicer accessions displayed both tolerant and sensitive phenotypes, supporting wild Cicer as a potential genetic pool for Al-tolerance improvement. To identify potential genes related to Al-tolerance in chickpea, genome-wide screening of multidrug and toxic compound extrusion (MATE) encoding genes was performed. Fifty-six MATE genes were identified in total, which can be divided into 4 major phylogenetic groups. Four chickpea MATE genes ( CaMATE1-4 ) were clustered with the previously characterized citrate transporters MtMATE66 and MtMATE69 in Medicago truncatula . Transcriptome data showed that CaMATE1-4 have diverse expression profiles, with CaMATE2 being root-specific. qRT-PCR analyses confirmed that CaMATE2 and CaMATE4 were highly expressed in root tips and were up-regulated upon Al treatment in all chickpea lines. Further measurement of carboxylic acids showed that malonic acid, instead of malate or citrate, is the major extruded acid by Cicer spp. root. Protein structural modeling analyses revealed that CaMATE2 has a divergent substrate-binding cavity from Arabidopsis AtFRD3, which may explain the different acid-secretion profile for chickpea. Pangenome survey showed that CaMATE1-4 have much higher genetic diversity in wild Cicer than that in cultivated chickpea. This first identification of CaMATE2 and CaMATE4 responsive to Al 3+ treatment in Cicer paves the way for future functional characterization of MATE genes in Cicer spp., and to facilitate future design of gene-specific markers for Al-tolerant line selection in chickpea breeding programs., 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 Jia, Pradeep, Vance, Zhang, Weir, Wei, Deng, Zhang, Xu, Zhao, Berger, Bell and Li.)

Protease inhibitors are widely studied since the unrestricted activity of proteases can cause extensive organ lesions. In particular, elastase activity is involved in the pathophysiology of acute lung injury, for example during SARS-CoV-2 infection, while serine proteases and thrombin-like proteases are involved in the development and/or pathology of the nervous system. Natural protease inhibitors have the advantage to be reversible and with few side effects and thus are increasingly considered as new drugs. Kunitz-type protease inhibitors (KTPIs), reported in the venom of various organisms, such as wasps, spiders, scorpions, and snakes, have been studied for their potent anticoagulant activity and widespread protease inhibitor activity. Putative KTPI anticoagulants have been identified in transcriptomic resources obtained for two blister beetle species, Lydus trimaculatus and Mylabris variabilis . The KTPIs of L. trimaculatus and M. variabilis were characterized by combined transcriptomic and bioinformatics methodologies. The full-length mRNA sequences were divided on the base of the sequence of the active sites of the putative proteins. In silico protein structure analyses of each group of translational products show the biochemical features of the active sites and the potential protease targets. Validation of these genes is the first step for considering these molecules as new drugs for use in medicine.

Background: Aldehyde dehydrogenases (ALDHs) are a family of NAD(P) + dependent enzymes that detoxify aldehydes by promoting their oxidation to respective carboxylic acids. The role of ALDH enzymes in various plant species has been extensively studied, revealing their critical role in salinity, drought, heat, and heavy metal stress tolerance. Despite their physiological significance, ALDH genes in Sorghum bicolor have yet to be studied thoroughly., Results: In this study, a total of 19 ALDH genes have been identified that have been grouped into ten families based on the criteria of the ALDH gene nomenclature committee. Segmental duplication assisted more in the enhancement of SbALDH gene family members than tandem duplication. All the identified SbALDH members made a cluster with monocot rice and maize in the phylogenetic tree rather than dicot species, suggesting the pre-eudicot-monocot separation of the ALDH superfamily members. The gene structure and protein domain were found to be mostly conserved in separate phylogenetic classes, indicating that each family played an important role in evolution. Expression analysis revealed that several SbALDHs were expressed in various tissues, developmental stages, and in response to abiotic stresses, indicating that they can play roles in plant growth, development, or stress adaptation. Interestingly, the majority of the SbALDH genes were found to be highly responsive to drought stress, and the SbALDH18B1 transcript showed maximum enhancement in all the stress conditions. The presence of cis-acting elements (mainly ABRE and MBS) in the promoter region of these genes might have a significant role in drought tolerance., Conclusions: Our findings add to the current understanding, evolutionary history, and contribution of SbALDHs in stress tolerance, and smooth the path of further functional validation of these genes., (© 2022. The Author(s).)

Male infertility is a global problem affecting a considerable part of the male population. Current guidelines and practices aimed at diagnosing the cause of this problem still have low diagnostic yield. As novel candidate genes for infertility emerge, their functional role needs to be investigated in patient populations. The present study aimed to investigate testis‑specific serine kinase 1B ( TSSK1B ), which was discovered in a previously diagnosed patient. Sanger sequencing of the coding regions and exon borders of TSSK1B was performed in a cohort of 100 male Bulgarian patients with unresolved infertility causes. Missense mutations were discovered in 10% of patients and were associated with clinical data on sperm dysmorphology. Two previously unreported mutations were discovered, p.3D>N and p.52F>L. All mutations were scored via in silico predictors and protein modelling using AlphaFold2. The present findings indicated an association between TSSK1B mutations and asthenoteratozoospermia, with further missense mutations in patients with azoospermia and teratozoospermia. Mutations in TSSK1B may be a cause of undiagnosed cases of male infertility and should be considered when molecular diagnostics are warranted.

Antifreeze proteins (AFPs) from the model crop, Brachypodium distachyon , allow freeze survival and attenuate pathogen-mediated ice nucleation. Intriguingly, Brachypodium AFP genes encode two proteins, an autonomous AFP and a leucine-rich repeat (LRR). We present structural models which indicate that ice-binding motifs on the ~13 kDa AFPs can "spoil" nucleating arrays on the ~120 kDa bacterial ice nucleating proteins used to form ice at high sub-zero temperatures. These models are consistent with the experimentally demonstrated decreases in ice nucleating activity by lysates from wildtype compared to transgenic Brachypodium lines. Additionally, the expression of Brachypodium LRRs in transgenic Arabidopsis inhibited an immune response to pathogen flagella peptides (flg22). Structural models suggested that this was due to the affinity of the LRR domains to flg22. Overall, it is remarkable that the Brachypodium genes play multiple distinctive roles in connecting freeze survival and anti-pathogenic systems via their encoded proteins' ability to adsorb to ice as well as to attenuate bacterial ice nucleation and the host immune response.

Hepatitis C virus (HCV), which infected 71 million worldwide and about 5%-6% are from Pakistan, is an ssRNA virus, responsible for end-stage liver disease. To date, no effective therapy is available to cure this disease. Hence, it is important to study the most prevalent genotypes infecting human population and design novel vaccine or small molecule inhibitors to control the infections associated with HCV. Therefore, in this study clinical samples ( n  = 35; HCV-3a) from HCV patients were subjected to Sanger sequencing method. The sequencing of the core gene, which is generally considered as conserved, involved in the detection, quantitation and genotyping of HCV was performed. Multiple mutations, that is, R46C, R70Q, L91C, G60E, N/S105A, P108A, N110I, S116V, G90S, A77G and G145R that could be linked with response to antiviral therapies were detected. Phylogenetic analysis suggests emerging viral isolates are circulating in Pakistan. Using ab initio modelling technique, we predicted the 3D structure of core protein and subjected to molecular dynamics simulation to extract the most stable conformation of the structure for further analysis. Immunoinformatic approaches were used to propose a multi-epitopes vaccine against HCV by using core protein. The vaccine constructs consist of nine CTL and three HTL epitopes joined by different linkers were docked against the two reported Toll-like receptors (TLR-3 and TLR-8). Docking of vaccine construct with TLR-3 and TLR-8 shows proper binding and in silico expression of the vaccine resulted in a CAI value of 0.93. These analyses suggest that specific immune responses may be produced by the proposed vaccine.Communicated by Ramaswamy H. Sarma.

Mutations in the RYR1 gene, encoding ryanodine receptor 1 (RyR1), are a well-known cause of Central Core Disease (CCD) and Multi-minicore Disease (MmD). We screened a cohort of 153 patients carrying an histopathological diagnosis of core myopathy (cores and minicores) for RYR1 mutation. At least one RYR1 mutation was identified in 69 of them and these patients were further studied. Clinical and histopathological features were collected. Clinical phenotype was highly heterogeneous ranging from asymptomatic or paucisymptomatic hyperCKemia to severe muscle weakness and skeletal deformity with loss of ambulation. Sixty-eight RYR1 mutations, generally missense, were identified, of which 16 were novel. The combined analysis of the clinical presentation, disease progression and the structural bioinformatic analyses of RYR1 allowed to associate some phenotypes to mutations in specific domains. In addition, this study highlighted the structural bioinformatics potential in the prediction of the pathogenicity of RYR1 mutations. Further improvement in the comprehension of genotype-phenotype relationship of core myopathies can be expected in the next future: the actual lack of the human RyR1 crystal structure paired with the presence of large intrinsically disordered regions in RyR1, and the frequent presence of more than one RYR1 mutation in core myopathy patients, require designing novel investigation strategies to completely address RyR1 mutation effect., (© 2022. The Author(s).)

G555F mutant of Fibrinogen A alpha-chain (FGA) is reported to be associated with kidney amyloidosis. In the current study, we have modelled the G555F mutant and examined the mutation's effect on the structural and functional level. We have also docked Vitamin C and D3 on the mutant's amyloidogenic region to identify if these vitamins can bind amyloidogenic regions. Further, we analyzed if they could prevent or modulate amyloid formation by stopping critical interactions in amyloidogenic regions in FGA. We used the wild type FGA model protein as a control. Our docking and molecular dynamics simulation results indicate stronger Vitamin D3 binding than Vitamin C to the amyloidogenic region of the mutant protein. The RMSD, radius of gyration, and RMSF values were higher for the G555F mutant than the FGA wild type protein. Overall, the results support these vitamins' potential as a therapeutic and anti-amyloidogenic agent for FGA renal amyloidosis., (© 2021. The Author(s), under exclusive licence to Springer Nature Switzerland AG.)

Human menin is a nuclear protein that participates in many cellular processes, as transcriptional regulation, DNA damage repair, cell signaling, cell division, proliferation, and migration, by interacting with many other proteins. Mutations of the gene encoding menin cause multiple endocrine neoplasia type 1 (MEN1), a rare autosomal dominant disorder associated with tumors of the endocrine glands. In order to characterize the structural and functional effects at protein level of the hundreds of missense variations, we investigated by computational methods the wild-type menin and more than 200 variants, predicting the amino acid variations that change secondary structure, solvent accessibility, salt-bridge and H-bond interactions, protein thermostability, and altering the capability to bind known protein interactors. The structural analyses are freely accessible online by means of a web interface that integrates also a 3D visualization of the structure of the wild-type and variant proteins. The results of the study offer insight into the effects of the amino acid variations in view of a more complete understanding of their pathological role.

Background: Recently, we have shown that seven genes, namely GBP5 , IRS2 , KRT4 , LINCOO707 , MRPL55 , RRS1 and SLC4A11, have prognostic power for the overall survival in ovarian cancer (OC)., Methods: We present an analysis on the association of these genes with any phenotypes and mutations indicative of involvement in female cancers and predict the structural and functional consequences of those SNPS using in silico tools., Results: These seven genes present with 976 SNPs/mutations that are associated with human cancers, out of which 284 related to female cancers. We have then analysed the mutation impact on amino acid polarity, charge and water affinity, leading to the identification of 30 mutations in gynaecological cancers where amino acid (aa) changes lead to opposite polarity, charges and water affinity. Out of these 30 mutations identified, only a missense mutation (i.e., R831C/R804C in uterine corpus endometrial carcinomas, UCEC) was suggestive of structural damage on the SLC4A11 protein., Conclusions: We demonstrate that the R831C/R804C mutation is deleterious and the predicted ΔΔG values suggest that the mutation reduces the stability of the protein. Future in vitro studies should provide further insight into the role of this transporter protein in UCEC.

Oomycete and fungal interactions with plants can be neutral, symbiotic or pathogenic with different impact on plant health and fitness. Both fungi and oomycetes can generate so-called effector proteins in order to successfully colonize the host plant. These proteins modify stress pathways, developmental processes and the innate immune system to the microbes' benefit, with a very different outcome for the plant. Investigating the biological and functional roles of effectors during plant-microbe interactions are accessible through bioinformatics and experimental approaches. The next generation protein modeling software RoseTTafold and AlphaFold2 have made significant progress in defining the 3D-structure of proteins by utilizing novel machine-learning algorithms using amino acid sequences as their only input. As these two methods rely on super computers, Google Colabfold alternatives have received significant attention, making the approaches more accessible to users. Here, we focus on current structural biology, sequence motif and domain knowledge of effector proteins from filamentous microbes and discuss the broader use of novel modelling strategies, namely AlphaFold2 and RoseTTafold, in the field of effector biology. Finally, we compare the original programs and their Colab versions to assess current strengths, ease of access, limitations and future applications.

The rapid developments in gene sequencing technologies achieved in the recent decades, along with the expansion of knowledge on the three-dimensional structures of proteins, have enabled the construction of proteome-scale databases of protein models such as the Genome3D and ModBase. Nevertheless, although gene products are usually expressed as individual polypeptide chains, most biological processes are associated with either transient or stable oligomerisation. In the PDB databank, for example, ~40% of the deposited structures contain at least one homo-oligomeric interface. Unfortunately, databases of protein models are generally devoid of multimeric structures. To tackle this particular issue, we have developed ProtCHOIR, a tool that is able to generate homo-oligomeric structures in an automated fashion, providing detailed information for the input protein and output complex. ProtCHOIR requires input of either a sequence or a protomeric structure that is queried against a pre-constructed local database of homo-oligomeric structures, then extensively analyzed using well-established tools such as PSI-Blast, MAFFT, PISA and Molprobity. Finally, MODELLER is employed to achieve the construction of the homo-oligomers. The output complex is thoroughly analyzed taking into account its stereochemical quality, interfacial stabilities, hydrophobicity and conservation profile. All these data are then summarized in a user-friendly HTML report that can be saved or printed as a PDF file. The software is easily parallelizable and also outputs a comma-separated file with summary statistics that can straightforwardly be concatenated as a spreadsheet-like document for large-scale data analyses. As a proof-of-concept, we built oligomeric models for the Mabellini Mycobacterium abscessus structural proteome database. ProtCHOIR can be run as a web-service and the code can be obtained free-of-charge at http://lmdm.biof.ufrj.br/protchoir., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

The recently emerged SARS-CoV2 caused a major pandemic of coronavirus disease (COVID-19). Non structural protein 1 (nsp1) is found in all beta coronavirus that cause severe respiratory disease. This protein is considered as a virulence factor and has an important role in pathogenesis. This study aims to elucidate the structural conformations of nsp1 to aid in the prediction of epitope sites and identification of important residues for targeted therapy against COVID-19. In this study, molecular modelling coupled with molecular dynamics simulations were performed to analyse the conformational landscape of nsp1 homologs of SARS-CoV1, SARS-CoV2 and MERS-CoV. Principal component analysis escorted by free energy landscape revealed that SARS-CoV2 nsp1 protein shows greater flexibility compared to SARS-CoV1 and MERS-CoV nsp1. Sequence comparison reveals that 28 mutations are present in SARS-CoV2 nsp1 protein compared to SARS-CoV1 nsp1. Several B-cell and T-cell epitopes were identified by an immunoinformatics approach. SARS-CoV2 nsp1 protein binds with the interface region of the palm and finger domain of POLA1 via hydrogen bonding and salt bridge interactions. Taken together, these in silico findings may help in the development of therapeutics specific against COVID-19., 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., (© 2021 Elsevier B.V. All rights reserved.)

Background: Boltzmann machines are energy-based models that have been shown to provide an accurate statistical description of domains of evolutionary-related protein and RNA families. They are parametrized in terms of local biases accounting for residue conservation, and pairwise terms to model epistatic coevolution between residues. From the model parameters, it is possible to extract an accurate prediction of the three-dimensional contact map of the target domain. More recently, the accuracy of these models has been also assessed in terms of their ability in predicting mutational effects and generating in silico functional sequences., Results: Our adaptive implementation of Boltzmann machine learning, adabmDCA, can be generally applied to both protein and RNA families and accomplishes several learning set-ups, depending on the complexity of the input data and on the user requirements. The code is fully available at https://github.com/anna-pa-m/adabmDCA . As an example, we have performed the learning of three Boltzmann machines modeling the Kunitz and Beta-lactamase2 protein domains and TPP-riboswitch RNA domain., Conclusions: The models learned by adabmDCA are comparable to those obtained by state-of-the-art techniques for this task, in terms of the quality of the inferred contact map as well as of the synthetically generated sequences. In addition, the code implements both equilibrium and out-of-equilibrium learning, which allows for an accurate and lossless training when the equilibrium one is prohibitive in terms of computational time, and allows for pruning irrelevant parameters using an information-based criterion., (© 2021. The Author(s).)

Lactic acid bacteria (LAB) are important microorganisms in food fermentation. In the food industry, bacteriophages (phages or bacterial viruses) may cause the disruption of LAB-dependent processes with product inconsistencies and economic losses. LAB phages use diverse adhesion devices to infect their host, yet the overall picture of host-binding mechanisms remains incomplete. Here, we aimed to determine the structure and topology of the adhesion devices of two lytic siphophages, OE33PA and Vinitor162, infecting the wine bacteria Oenococcus oeni . These phages possess adhesion devices with a distinct composition and morphology and likely use different infection mechanisms. We primarily used AlphaFold2, an algorithm that can predict protein structure with unprecedented accuracy, to obtain a 3D model of the adhesion devices' components. Using our prior knowledge of the architecture of the LAB phage host-binding machineries, we also reconstituted the topology of OE33PA and Vinitor162 adhesion devices. While OE33PA exhibits original structures in the assembly of its bulky adhesion device, Vinitor162 harbors several carbohydrate-binding modules throughout its long and extended adhesion device. Overall, these results highlight the ability of AlphaFold2 to predict protein structures and illustrate its great potential in the study of phage structures and host-binding mechanisms.

The metabolism of polyphenolic polymers is essential to the development and response to environmental changes of organisms from all kingdoms of life, but shows particular diversity in plants. In contrast to other biopolymers, whose polymerisation is catalysed by homologous gene families, polyphenolic metabolism depends on phenoloxidases, a group of heterogeneous oxidases that share little beyond the eponymous common substrate. In this review, we provide an overview of the differences and similarities between phenoloxidases in their protein structure, reaction mechanism, substrate specificity, and functional roles. Using the example of laccases (LACs), we also performed a meta-analysis of enzyme kinetics, a comprehensive phylogenetic analysis and machine-learning based protein structure modelling to link functions, evolution, and structures in this group of phenoloxidases. With these approaches, we generated a framework to explain the reported functional differences between paralogs, while also hinting at the likely diversity of yet undescribed LAC functions. Altogether, this review provides a basis to better understand the functional overlaps and specificities between and within the three major families of phenoloxidases, their evolutionary trajectories, and their importance for plant primary and secondary metabolism., 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 Blaschek and Pesquet.)

Purple acid phosphatases belong to metallo-phosphatase family. Intracellular phosphatases are crucial for phosphorus (P) distribution in the cell and are highly induced in phosphorus-deprived conditions in the soil. Disparate PAP isoforms exist within discrete subcellular compartments in Setaria italica and their expression in P deprived conditions fosters phosphorus amelioration. We isolated the SiPAP18 gene and developed the homology SiPAP18 protein model based on the crystal structure of the Kidney bean PvPAP (PDB ID: 2QFP) as template (sequence similarity 42.7%) using Modeller 9.12 with adequate validation. Structure model analysis shows the significance of five conserved signatures with seven metal-paired amino acid residues during P-deprivation induced phosphorus amelioration., (© 2021 Biomedical Informatics.)

Parasitic nematodes constitute one of the major threats to human health, causing diseases of major socioeconomic importance worldwide. Recent estimates indicate that more than 1 billion people are infected with parasitic nematodes around the world. Current measures to combat parasitic nematode infections include anthelmintic drugs. However, heavy exposure to anthelmintics has selected populations of livestock parasitic nematodes that are no longer susceptible to the drugs, rendering several anthelmintics useless for parasitic nematode control in many areas of the world. The rapidity with which anthelmintic resistance developed in response to these drugs suggests that increasing the selective pressure on human parasitic nematodes will also rapidly generate resistant worm populations. Therefore, development of new anthelmintics is of major importance before resistance becomes widespread in human parasitic nematode populations. G-Protein Coupled Receptors (GPCRs) represent an important target for many pharmacological interventions due to their ubiquitous expression in various cell types. GPCRs contribute to numerous physiological processes, and their ligand binding sites located on cell surfaces make them accessible targets and attractive substrates in terms of druggability. In fact, ∼35 % of Food and Drug Administration (FDA) and European Medicines Agency (EMA) approved drugs target GPCRs and their associated proteins, with over 300 additional drugs targeting GPCRs at the clinical trial stage. Nematode Chemosensory GPCRs (NemChRs) are unique to nematodes, and therefore represent ideal substrates for target-based drug discovery. Here we set out to identify NemChRs that are transcriptionally active inside the host, and to use these NemChRs in a reverse pharmacological screen to impede parasitic development. Our data identified several NemChRs, and we focused on one that was expressed in neuronal cells and exhibited the highest fold change in transcription after host activation. Next, we performed homology modelling and molecular dynamics simulations of this NemChR in order to conduct a virtual screening campaign to identify candidate drug targets which were ranked and selected for experimental testing in bioassays. Taken together, our results identify and characterize a candidate NemChR drug target, and provide a chemogenomic pipeline for identifying nematicide substrates., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

Background: Animal mitochondria play a central role in energy production in the cells through the oxidative phosphorylation (OXPHOS) pathway. Recent studies of selection on different mitochondrial OXPHOS genes have revealed the adaptive implications of amino acid changes in these subunits. In hares, climatic variation and/or introgression were suggested to be at the origin of such adaptation. Here we looked for evidence of positive selection in three mitochondrial OXPHOS genes, using tests of selection, protein structure modelling and effects of amino acid substitutions on the protein function and stability. We also used statistical models to test for climate and introgression effects on sites under positive selection., Results: Our results revealed seven sites under positive selection in ND4 and three sites in Cytb. However, no sites under positive selection were observed in the COX1 gene. All three subunits presented a high number of codons under negative selection. Sites under positive selection were mapped on the tridimensional structure of the predicted models for the respective mitochondrial subunit. Of the ten amino acid replacements inferred to have evolved under positive selection for both subunits, six were located in the transmembrane domain. On the other hand, three codons were identified as sites lining proton translocation channels. Furthermore, four codons were identified as destabilizing with a significant variation of Δ vibrational entropy energy between wild and mutant type. Moreover, our PROVEAN analysis suggested that among all positively selected sites two fixed amino acid replacements altered the protein functioning. Our statistical models indicated significant effects of climate on the presence of ND4 and Cytb protein variants, but no effect by trans-specific mitochondrial DNA introgression, which is not uncommon in a number of hare species., Conclusions: Positive selection was observed in several codons in two OXPHOS genes. We found that substitutions in the positively selected codons have structural and functional impacts on the encoded proteins. Our results are concordantly suggesting that adaptations have strongly affected the evolution of mtDNA of hare species with potential effects on the protein function. Environmental/climatic changes appear to be a major trigger of this adaptation, whereas trans-specific introgressive hybridization seems to play no major role for the occurrence of protein variants.

Objective: In current study we performed sequencing of palm domain of HCV-NS5B gene, its ancestral analysis along with amino acids substitution analysis. These analysis were done to find the molecular basis of the viral resistance against Sofosbuvir drug., Methods: Blood samples from individuals with chronic Hepatitis C infection that were resistant to Sofosbuvir were collected. The samples were processed for their molecular characterization that included RNA extraction, Complementary DNA (cDNA) synthesize, nested PCR, gel elution, Sequencing, ancestral and 3D structure analysis., Results: Evolutionary analysis revealed that current study sequences (QAU-01, QAU-02) clustered with a previously studied sequence, KY971494.1. Moreover, we reports multiple novel amino acid substitutions in the palm domain of NS5B gene such as Ile116Val, Asn117Gly, Glu246Ala, Val252Ala, Glu258Gln, Cys262Leu, Ser269Arg, Ala272Thr, Ile293Leu, Lys304Arg, Asn307Gly, Ala338Val and Arg345Gly in our query sequence (QAU-01). At 246 and 269 position in (QAU-02), no substitution was observed., Conclusions: We have noticed that the current sequences are relatively emerging and could have been originated from aforementioned sequence recently. Based on the current results, we suggests that these substitutions could be associated with structural or functional impairment of protein and could also be may be considered as resistance associated substitutions (RAS) to Sofosbuvir drug., (Copyright © 2021. Published by Elsevier Ltd.)

The Paget disease (PDB; OMIM is 167250) is a chronic bone disease caused by pathogenic mutations in Sequestome1/p62 (SQSTM1) gene. This study has aimed to interpret the relationship of PDB linked SQSTM1 mutations with protein structure and its molecular dynamic features. The disease causative missense mutations were initially collected, and then analyzed for their, exonic and domain distribution, impact on secondary and tertiary structures, and their ability on protein-ligand interactions, using a combination of systems biology approaches. Our results show that most PDB linked SQSTM1 missense mutations affect amino acid residues clustered within or near the UBA domain (aa 389-434), which participates in the ubiquitination of substrates. We also report that the majority mutations occurred in α-helices over β-strands but their effects on the secondary structure were mostly neutral. Global tertiary structure deviations were minimal; however, at amino acid residue level minor structural changes were evident. The molecular dynamics simulation analysis showed that both PB1 and UBA domains were under constant structural fluctuations resulting in closed form conformation of SQSMT1 protein structure, when it is bound to PRKCI ligand. We also found salt bridge conformation changes in the UBA domain of SQSTM1 mutants when they bound to the PRKCI interactor protein. This finding suggests the possibility that mutations in SQSTM1 could impair its ability to ubiquitinate the substrates, eventually affecting autophagy and apoptosis, especially in mature osteoclasts. This study presents the additional insight into structure and function relationship between SQSTM1 mutations and PDB pathogenesis. Communicated by Ramaswamy H. Sarma.

Mitochondria was long thought to be an "end function" organelle that regulated the metabolic flux and apoptosis in the cell. However, with the discovery of the mitochondrial peptide (MDP) humanin (HN/MTRNR2), the cytoprotective and pro-survival applications of MDPs have taken the forefront of therapeutic and diagnostic research. However, the regulation of humanin-like MDPs in natural model systems that can tolerate lethal environmental and cytotoxic insults remains to be investigated. Red-eared sliders are champion anaerobes that can withstand three continuous months of anoxia followed by rapid bouts of oxygen reperfusion without incurring cellular damage. Freshwater turtles employ extensive physiological and biochemical strategies to combat anoxia, with metabolic rate depression and a global enhancement of antioxidant and cytoprotective pathways being the two most important contributors. The main aim of this study was to uncover and characterize the humanin-homologue in freshwater turtles as well as investigate the differential regulation of humanin in response to short and long-term oxygen deprivation. In this study we have used de novo and homology-based protein modelling to elucidate the putative structure of humanin in red-eared sliders as well as an ELISA and western immunoblotting to confirm the protein abundance in the turtle brain and six peripheral tissues during control, 5 h, and 20 h anoxia (n = 4/group). We found that a humanin-homologue (TSE-humanin) is present in red-eared sliders and it may play a cytoprotective role against oxidative damage.

Congenital hypofibrinogenemia is a rare bleeding disorder characterized by a proportional decrease of functional and antigenic fibrinogen levels. Hypofibrinogenemia can be considered the phenotypic expression of heterozygous loss of function mutations occurring within one of the three fibrinogen genes ( FGA , FGB , and FGG ). Clinical manifestations are highly variable; most patients are usually asymptomatic, but may appear with mild to severe bleeding or thrombotic complications. We have sequenced all exons of the FGA , FGB , and FGG genes using the DNA isolated from the peripheral blood in two unrelated probands with mild hypofibrinogenemia. Coagulation screening, global hemostasis, and functional analysis tests were performed. Molecular modeling was used to predict the defect of synthesis and structural changes of the identified mutation. DNA sequencing revealed a novel heterozygous variant c.1421G>A in exon 8 of the FGB gene encoding a Bβ chain (p.Trp474Ter) in both patients. Clinical data from patients showed bleeding episodes. Protein modelling confirmed changes in the secondary structure of the molecule, with the loss of three β sheet arrangements. As expected by the low fibrinogen levels, turbidity analyses showed a reduced fibrin polymerisation and imaging difference in thickness fibrin fibers. We have to emphasize that our patients have a quantitative fibrinogen disorder; therefore, the reduced function is due to the reduced concentration of fibrinogen, since the Bβ chains carrying the mutation predicted to be retained inside the cell. The study of fibrinogen molecules using protein modelling may help us to understand causality and effect of novel genetic mutations.

The insect chemosensory repertoires of Odorant Receptors (ORs) and Gustatory Receptors (GRs) together represent one of the largest families of ligand-gated ion channels. Previous analyses have identified homologous 'Gustatory Receptor-Like' (GRL) proteins across Animalia, but the evolutionary origin of this novel class of ion channels is unknown. We describe a survey of unicellular eukaryotic genomes for GRLs, identifying several candidates in fungi, protists and algae that contain many structural features characteristic of animal GRLs. The existence of these proteins in unicellular eukaryotes, together with ab initio protein structure predictions, provide evidence for homology between GRLs and a family of uncharacterized plant proteins containing the DUF3537 domain. Together, our analyses suggest an origin of this protein superfamily in the last common eukaryotic ancestor., Competing Interests: RB, CD, DM No competing interests declared, (© 2020, Benton et al.)

Background: Preeclampsia (PE) is a frequently occurring multisystemic disease affecting ~5% of pregnancies. PE patients may develop HELLP syndrome (haemolysis, elevated liver enzymes, and low platelet), a mother and foetus life-threatening condition. Research into HELLP's genetic origin has been relatively unsuccessful, mainly because normal placental function and blood pressure regulation involve the fine-regulation of hundreds of genes., Objective: To identify new genes and mutations constituting potential biomarkers for HELLP syndrome., Study Design: The present case-control study involved whole-exome sequencing of 79 unrelated HELLP women. Candidate variants were screened in a control population constituted by 176 individuals. Stringent bioinformatics filters were used for selecting potentially etiological sequence variants in a subset of 487 genes. We used robust in silico mutation modelling for predicting the potential effect on protein structure., Results: We identified numerous sequence variants in genes related to angiogenesis/coagulation/blood pressure regulation, cell differentiation/communication/adhesion, cell cycle and transcriptional gene regulation, extracellular matrix biology, lipid metabolism and immunological response. Five sequence variants generated premature stop codons in genes playing an essential role in placental physiology (STOX1, PDGFD, IGF2, MMP1 and DNAH11). Six variants (ERAP1- p.Ile915Thr, ERAP2- p.Leu837Ser, COMT-p.His192Gln, CSAD-p.Pro418Ser, CDH1- p.Ala298Thr and CCR2-p.Met249Lys) led to destabilisation of protein structure as they had significant energy and residue interaction-related changes. We identified at least two mutations in 57% of patients, arguing in favour of a polygenic origin for the HELLP syndrome., Conclusion: Our results provide novel evidence regarding PE/HELLP's genetic origin, leading to new biomarkers, having potential clinical usefulness, being proposed., (Copyright © 2020 International Society for the Study of Hypertension in Pregnancy. Published by Elsevier B.V. All rights reserved.)