QMEANDisCo-distance constraints applied on model quality estimation

Phage-Choline Kinase Inhibitor Combination to Control Pseudomonas aeruginosa: A Promising Combo

BACKGROUND

Pseudomonas aeruginosa is one of the most prevalent opportunistic pathogens in humans that has thrived and proved to be difficult to control in this "post-antibiotic era." Antibiotic alternatives are necessary for fighting against this resilient bacterium. Even though phages might not be "the wonder drug" that solves everything, they still provide a viable option to combat P. aeruginosa and curb the threat it imposes.

MAIN FINDINGS

The combination of antibiotics with phages, however, poses a propitious treatment option for P. aeruginosa. Choline kinase (ChoK) is the enzyme that synthesizes phosphorylcholine subsequently incorporated into lipopolysaccharide located at the outer membrane of gram-negative bacteria. Recently, inhibition of ChoKs has been proposed as a promising antibacterial strategy. Successful docking of Hemicholinium-3, a choline kinase inhibitor, to the model structure of P. aeruginosa ChoK also supports the use of this inhibitor or its derivatives to inhibit the growth of this microorganism.

CONCLUSION

Therefore, the combination of the novel antimicrobial "choline kinase inhibitors (ChoKIs)" with a phage cocktail or synthetic phages as a potential treatment for P. aeruginosa infection has been proposed.

Genotype-phenotype analysis, and assessment of the importance of the zinc-binding site in PHEX in Japanese patients with X-linked hypophosphatemic rickets using 3D structure modeling

X-linked hypophosphatemic rickets (XLH) is an inheritable type of rickets caused by inactivating variants in the phosphate regulating endopeptidase homolog X-linked (PHEX) gene, which results in the overproduction of fibroblast growth factor 23 (FGF23). The mechanism by which PHEX impairment leads to FGF23 overproduction is unknown. Because little is known regarding the genotype-phenotype correlation in Japanese XLH, we summarized the available clinical and genetic data and analyzed the genotype-phenotype relationships using 3-dimensional (3D) structure modeling to clarify the XLH pathophysiology. We retrospectively reviewed the clinical features and performed genetic analysis of 39 Japanese patients with XLH from 28 unrelated pedigrees carrying any known or novel PHEX variant. To predict changes in the 3D structure of mutant PHEX, we constructed a putative 3D model of each mutant and evaluated the effect of structural alteration by genotype-phenotype correlation analysis. Genetic analysis revealed 23 PHEX variants, including eight novel variants. They were associated with high i-FGF23 levels, hypophosphatemia, phosphaturia, high alkaline phosphatase levels, and short stature. No gene dosage effect or genotype-phenotype correlation was observed when truncating and non-truncating variants were compared. However, the conservation of the zinc-binding site and cavity in PHEX had an impact on the elevation of i-FGF23 levels. Via genotype-phenotype relationship analysis using 3D modeling, we showed that the zinc-binding site and cavity in PHEX can play a critical role in its function. These findings provide new genetic clues for investigating the function of PHEX and the pathogenesis of XLH.

Old pesticide, new use: Smart and safe enantiomer of isocarbophos in locust control

Locust plagues are still worldwide problems. Selecting active enantiomers from current chiral insecticides is necessary for controlling locusts and mitigating the pesticide pollution in agricultural lands. Herein, two enantiomers of isocarbophos (ICP) were separated and the enantioselectivity in insecticidal activity against the pest Locusta migratoria manilensis (L. migratoria) and mechanisms were investigated. The significant difference of LD₅₀ between (+)-ICP (0.609 mg/kg bw) and (-)-ICP (79.412 mg/kg bw) demonstrated that (+)-ICP was a more effective enantiomer. The enantioselectivity in insecticidal activity of ICP enantiomers could be attributed to the selective affinity to acetylcholinesterase (AChE). Results of in vivo and in vitro assays suggested that AChE was more sensitive to (+)-ICP. In addition, molecular docking showed that the -CDOKER energies of (+)-ICP and (-)-ICP were 25.6652 and 24.4169, respectively, which suggested a stronger affinity between (+)-ICP and AChE. Significant selectivity also occurred in detoxifying enzymes activities (carboxylesterases (CarEs) and glutathione S-transferases (GSTs)) and related gene expressions. Suppression of detoxifying enzymes activities with (+)-ICP treatment suggested that (-)-ICP may induce the detoxifying enzyme-mediated ICP resistance. A more comprehensive understanding of the enantioselectivity of ICP is necessary for improving regulation and risk assessment of ICP.

Physicochemical and structural characterization, epitope mapping and vaccine potential investigation of a new protein containing Tetratrico Peptide Repeats of Acinetobacter baumannii: An in-silico and in-vivo approach

Acinetobacter baumannii is an opportunistic multidrug-resistant pathogen that causes a significant mortality rate. The proteins containing Tetratrico Peptide Repeats (TPRs) are involved in the pathogenicity and virulence of bacteria and have different roles such as transfer of bacterial virulence factors to host cells, binding to the host cells and inhibition of phagolysosomal maturation. So, in this study, physicochemical properties of a new protein containing TPRs in A. baumannii which was named PcTPRs1 by this study were characterized and its 3D structure was predicted by in-silico tools. The protein B and T cell epitopes were mapped and its vaccine potential was in-silico and in-vivo investigated. Domain analysis indicated that the protein contains the Flp pilus assembly protein TadD domain which has three TPRs. The helix is dominant in the protein structure, and this protein is an outer membrane antigen which, is extremely conserved among A. baumannii strains; thus, has good properties to be applied as a recombinant vaccine. The best-predicted and refined model was applied in ligand-binding sites and conformational epitopes prediction. Based on epitope mapping results, several epitopes were characterized which could stimulate both immune systems. BLAST results showed the introduced epitopes are completely conserved among A. baumannii strains. The in-vivo analysis indicates that a 101 amino acid fragment of the protein which contains the best selected epitope, can produce a good protectivity against A. baumannii as well as the whole TPR protein and thus could be investigated as an effective subunit and potential vaccines.

Discovery of a novel and a rare Kristen rat sarcoma viral oncogene homolog (KRAS) gene mutation in colorectal cancer patients

Colorectal cancer (CRC) is one of the most important causes of morbidity and mortality in the developed world and is gradually more frequent in the developing world including Saudi Arabia. According to the Saudi Cancer Registry report 2015, CRC is the most common cancer in men (14.9%) and the second most prevalent cancer. Oncogenic mutations in the KRAS gene play a central role in tumorigenesis and are mutated in 30-40% of all CRC patients. To explore the prevalence of KRAS gene mutations in the Saudi population, we collected 80 CRC tumor tissues and sequenced the KRAS gene using automated sequencing technologies. The chromatograms presented mutations in 26 patients (32.5%) in four different codons, that is, 12, 13, 17, and 31. Most of the mutations were identified in codon 12 in 16 patients (61.5% of all mutations). We identified a novel mutation c.51 G>A in codon 17, where serine was substituted by arginine (S17R) in four patients. We also identified a very rare mutation, c.91 G>A, in which glutamic acid was replaced by lysine (E31K) in three patients. In conclusion, our findings further the knowledge about KRAS mutations in different ethnic groups is indispensable to fully understand their role in the development and progression of CRC.

Characterization of a MOB1 Homolog in the Apicomplexan Parasite Toxoplasma gondii

Simple Summary

Monopolar spindle One Binder1 (MOB1) proteins regulate key cellular functions, namely cell multiplication and cell division. The unicellular parasite Toxoplasma gondii transitions between several morphological stages, with the need to control the number of parasites in its cellular environment. We hypothesized that MOB1 proteins could participate in the regulation of the T. gondii life cycle, having identified one MOB1 protein (TgMOB1) coded in its genome. However, this study shows that TgMOB1 presents divergent features. While in organisms studied to date the lack of MOB1 has led to cell division defects, this did not occur in T. gondii in vitro cultures where mob1 was not an essential gene. Additionally, the identification of TgMOB1 proximity interacting partners detected novel MOB1 interactors. Still, TgMOB1 localizes to the region between the new-forming nuclei during cell division, and T. gondii parasites multiply slower with TgMOB1 overexpression and faster when there is a lack of TgMOB1, indicating an intricate role for TgMOB1 in T. gondii. This study uncovers new features of the T. gondii biology, a zoonotic parasite and model organism for the phylum Apicomplexa, and highlights the complex roles MOB1 proteins may assume, with possible implications for disease processes.

Abstract

Monopolar spindle One Binder1 (MOB1) proteins are conserved components of the tumor-suppressing Hippo pathway, regulating cellular processes such as cytokinesis. Apicomplexan parasites present a life cycle that relies on the parasites’ ability to differentiate between stages and regulate their proliferation; thus, Hippo signaling pathways could play an important role in the regulation of the apicomplexan life cycle. Here, we report the identification of one MOB1 protein in the apicomplexan Toxoplasma gondii. To characterize the function of MOB1, we generated gain-of-function transgenic lines with a ligand-controlled destabilization domain, and loss-of-function clonal lines obtained through CRISPR/Cas9 technology. Contrary to what has been characterized in other eukaryotes, MOB1 is not essential for cytokinesis in T. gondii. However, this picture is complex since we found MOB1 localized between the newly individualized daughter nuclei at the end of mitosis. Moreover, we detected a significant delay in the replication of overexpressing tachyzoites, contrasting with increased replication rates in knockout tachyzoites. Finally, using the proximity-biotinylation method, BioID, we identified novel members of the MOB1 interactome, a probable consequence of the observed lack of conservation of some key amino acid residues. Altogether, the results point to a complex evolutionary history of MOB1 roles in apicomplexans, sharing properties with other eukaryotes but also with divergent features, possibly associated with their complex life cycle.

Homology Modeling and Molecular Dynamics Simulations of Trypanosoma cruzi Phosphodiesterase b1

Cyclic nucleotide phosphodiesterases have been implicated in the proliferation, differentiation and osmotic regulation of trypanosomatids; in some trypanosomatid species, they have been validated as molecular targets for the development of new therapeutic agents. Because the experimental structure of Trypanosoma cruzi PDEb1 (TcrPDEb1) has not been solved so far, an homology model of the target was created using the structure of Trypanosoma brucei PDEb1 (TbrPDEb1) as a template. The model was refined by extensive enhanced sampling molecular dynamics simulations, and representative snapshots were extracted from the trajectory by combined clustering analysis. This structural ensemble was used to develop a structure-based docking model of the target. The docking accuracy of the model was validated by redocking and cross-docking experiments using all available crystal structures of TbrPDEb1, whereas the scoring accuracy was validated through a retrospective screen, using a carefully curated dataset of compounds assayed against TbrPDEb1 and/or TcrPDEb1. Considering the results from in silico validations, the model may be applied in prospective virtual screening campaigns to identify novel hits, as well as to guide the rational design of potent and selective inhibitors targeting this enzyme.

Characterisation of a putative M23-domain containing protein in Mycobacterium tuberculosis

Mycobacterium tuberculosis, the causative agent of tuberculosis remains a global health concern, further compounded by the high rates of HIV-TB co-infection and emergence of multi- and extensive drug resistant TB, all of which have hampered efforts to eradicate this disease. As a result, novel anti-tubercular interventions are urgently required, with the peptidoglycan component of the M. tuberculosis cell wall emerging as an attractive drug target. Peptidoglycan M23 endopeptidases can function as active cell wall hydrolases or degenerate activators of hydrolases in a variety of bacteria, contributing to important processes such as bacterial growth, division and virulence. Herein, we investigate the function of the Rv0950-encoded putative M23 endopeptidase in M. tuberculosis. In silico analysis revealed that this protein is conserved in mycobacteria, with a zinc-binding catalytic site predictive of hydrolytic activity. Transcript analysis indicated that expression of Rv0950c was elevated during lag and log phases of growth and reduced in stationary phase. Deletion of Rv0950c yielded no defects in growth, colony morphology, antibiotic susceptibility or intracellular survival but caused a reduction in cell length. Staining with a monopeptide-derived fluorescent D-amino acid, which spatially reports on sites of active PG biosynthesis or repair, revealed an overall reduction in uptake of the probe in ΔRv0950c. When stained with a dipeptide probe in the presence of cell wall damaging agents, the ΔRv0950c mutant displayed reduced sidewall labelling. As bacterial peptidoglycan metabolism is important for survival and pathogenesis, the role of Rv0950c and other putative M23 endopeptidases in M. tuberculosis should be explored further.

A Comparative Study between Spanish and British SARS-CoV-2 Variants

The study of the interaction between the SARS-CoV-2 spike protein and the angiotensin-converting enzyme 2 (ACE2) receptor is key to understanding binding affinity and stability. In the present report, we sought to investigate the differences between two already sequenced genome variants (Spanish and British) of SARS-CoV-2. Methods: In silico model evaluating the homology, identity and similarity in the genome sequence and the structure and alignment of the predictive spike by computational docking methods. Results: The identity results between the Spanish and British variants of the Spike protein were 28.67%. This close correspondence in the results between the Spanish and British SARS-CoV-2 variants shows that they are very similar (99.99%). The alignment obtained results in four deletions. There were 23 nucleotide substitutions also predicted which could affect the functionality of the proteins produced from this sequence. The interaction between the binding receptor domain from the spike protein and the ACE2 receptor produces some of the mutations found and, therefore, the energy of this ligand varies. However, the estimated antigenicity of the British variant is higher than its Spanish counterpart. Conclusions: Our results indicate that minimal mutations could interfere in the infectivity of the virus due to changes in the fitness between host cell recognition and interaction proteins. In particular, the N501Y substitution, situated in the RBD of the spike of the British variant, might be the reason for its extraordinary infective potential.

Global diversity of the gene encoding the Pfs25 protein-a Plasmodium falciparum transmission-blocking vaccine candidate

Background

Vaccines against the sexual stages of the malarial parasite Plasmodium falciparum are indispensable for controlling malaria and abrogating the spread of drug-resistant parasites. Pfs25, a surface antigen of the sexual stage of P. falciparum, is a leading candidate for transmission-blocking vaccine development. While clinical trials have reported that Pfs25-based vaccines are safe and effective in inducing transmission-blocking antibodies, the extent of the genetic diversity of Pfs25 in malaria endemic populations has rarely been studied. Thus, this study aimed to investigate the global diversity of Pfs25 in P. falciparum populations.

Methods

A database of 307 Pfs25 sequences of P. falciparum was established. Population genetic analyses were performed to evaluate haplotype and nucleotide diversity, analyze haplotypic distribution patterns of Pfs25 in different geographical populations, and construct a haplotype network. Neutrality tests were conducted to determine evidence of natural selection. Homology models of the Pfs25 haplotypes were constructed, subjected to molecular dynamics (MD), and analyzed in terms of flexibility and percentages of secondary structures.

Results

The Pfs25 gene of P. falciparum was found to have 11 unique haplotypes. Of these, haplotype 1 (H1) and H2, the major haplotypes, represented 70% and 22% of the population, respectively, and were dominant in Asia, whereas only H1 was dominant in Africa, Central America, and South America. Other haplotypes were rare and region-specific, resulting in unique distribution patterns in different geographical populations. The diversity in Pfs25 originated from ten single-nucleotide polymorphism (SNP) loci located in the epidermal growth factor (EGF)-like domains and anchor domain. Of these, an SNP at position 392 (GGA/GCA), resulting in amino acid substitution 131 (Gly/Ala), defined the two major haplotypes. The MD results showed that the structures of H1 and H2 variants were relatively similar. Limited polymorphism in Pfs25 could likely be due to negative selection.

Conclusions

The study successfully established a Pfs25 sequence database that can become an essential tool for monitoring vaccine efficacy, designing assays for detecting malaria carriers, and conducting epidemiological studies of P. falciparum. The discovery of the two major haplotypes, H1 and H2, and their conserved structures suggests that the current Pfs25-based vaccines could be used globally for malaria control.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-05078-6.

ABC transporter slr0982 affects response of Synechocystis sp. PCC 6803 to oxidative stress caused by methyl viologen

The exposure of methyl viologen (a bipyridine salt) can lead to the production of reactive oxygen species, causing oxidative stress to organisms. ABC transporters have been reported to be involved in multi-drug resistance and have a role in MV detoxification. Here, we performed a protein structure simulation of the Slr0982 protein encoding ABC transporters, and confirmed that the region from Phe57 to Gln257 was the ABC transporter-type domain of the Slr0982 protein. The results of protein sequence alignment showed that Slr0982 protein was similar to Slr2108 protein (polysialic acid transport ATP-binding protein) and Slr0354 protein (ABC transporter). We reported that the mutation of slr0982 reduced the tolerance of Synechocystis sp. PCC 6803 to oxidative stress induced by methyl viologen. The deletion of slr0982 reduced the ability of cells to resist oxidative stress. Our data confirmed that the deletion of slr0982 could affect the concentration of exopolysaccharide and the expression of some genes related to carbohydrate metabolism, thus decreasing polysaccharide transport. Importantly, the exogenous addition of exopolysaccharide extracted from wild type can effectively reduce the oxidative damage to Δslr0982 by methyl viologen. This study expands the role of ABC transporters in MV-induced oxidative stress and provides an insight into the further analysis of the response of cyanobacteria to oxidative stress.

Current Approaches in Supersecondary Structures Investigation

Proteins expressed during the cell cycle determine cell function, topology, and responses to environmental influences. The development and improvement of experimental methods in the field of structural biology provide valuable information about the structure and functions of individual proteins. This work is devoted to the study of supersecondary structures of proteins and determination of their structural motifs, description of experimental methods for their detection, databases, and repositories for storage, as well as methods of molecular dynamics research. The interest in the study of supersecondary structures in proteins is due to their autonomous stability outside the protein globule, which makes it possible to study folding processes, conformational changes in protein isoforms, and aberrant proteins with high productivity.

In silico evaluation of molecular virus-virus interactions taking place between Cotton leaf curl Kokhran virus- Burewala strain and Tomato leaf curl New Delhi virus

Background

Cotton leaf curl disease (CLCuD) is a disease of cotton caused by begomoviruses, leading to a drastic loss in the annual yield of the crop. Pakistan has suffered two epidemics of this disease leading to the loss of billions in annual exports. The speculation that a third epidemic of CLCuD may result as consequence of the frequent occurrence of Tomato leaf curl New Delhi virus (ToLCNDV) and Cotton leaf curl Kokhran Virus-Burewala Strain (CLCuKoV-Bu) in CLCuD infected samples, demand that the interactions taking between the two viruses be properly evaluated. This study is designed to assess virus-virus interactions at the molecular level and determine the type of co-infection taking place.

Methods

Based on the amino acid sequences of the gene products of both CLCuKoV-Bu and ToLCNDV, protein structures were generated using different software, i.e., MODELLER, I-TASSER, QUARKS, LOMETS and RAPTORX. A consensus model for each protein was selected after model quality assessment using ERRAT, QMEANDisCo, PROCHECK Z-Score and Ramachandran plot analysis. The active and passive residues in the protein structures were identified using the CPORT server. Protein–Protein Docking was done using the HADDOCK webserver, and 169 Protein–Protein Interaction (PPIs) were performed between the proteins of the two viruses. The docked complexes were submitted to the PRODIGY server to identify the interacting residues between the complexes. The strongest interactions were determined based on the HADDOCK Score, Desolvation energy, Van der Waals Energy, Restraint Violation Energy, Electrostatic Energy, Buried Surface Area and Restraint Violation Energy, Binding Affinity and Dissociation constant (K_d). A total of 50 ns Molecular Dynamic simulations were performed on complexes that exhibited the strongest affinity in order to validate the stability of the complexes, and to remove any steric hindrances that may exist within the structures.

Results

Our results indicate significant interactions taking place between the proteins of the two viruses. Out of all the interactions, the strongest were observed between the Replication Initiation protein (Rep) of CLCuKoV-Bu with the Movement protein (MP), Nuclear Shuttle Protein (NSP) of ToLCNDV (DNA-B), while the weakest were seen between the Replication Enhancer protein (REn) of CLCuKoV-Bu with the REn protein of ToLCNDV. The residues identified to be taking a part in interaction belonged to domains having a pivotal role in the viral life cycle and pathogenicity. It maybe deduced that the two viruses exhibit antagonistic behavior towards each other, and the type of infection may be categorised as a type of Super Infection Exclusion (SIE) or homologous interference. However, further experimentation, in the form of transient expression analysis, is needed to confirm the nature of these interactions and increase our understanding of the direct interactions taking place between two viruses.

Biochemical evaluation and molecular docking assessment of glucosamines from Neocarya macrophylla fruits against Naja nigricollis venom

Two new glucosamines, Microphyllose A and B were isolated from the chloroform fraction of Neocarya macrophylla fruit using flash column chromatography. The structures of these compounds were elucidated based on chemical tests and the analysis of their spectral data (IR, 1D- & 2D-NMR). The compounds have demonstrated significant (p < 0.05) antivenom activity against Naja nigricollis venom with 60 and 80% protection, respectively. When subjected to molecular docking, the compounds have demonstrated different binding affinities against three toxins (phospholipase A2, neurotoxin and cardiotoxin) from Naja nigricollis venom and they were further screened for ADMET analysis based on Lipinski's and Veber's rule and the compounds have failed absorptivity for oral medications. To the best of our knowledge, this is the first report of isolation and molecular docking analysis of these compounds from medicinal plants.

Discovery of new 2-phenyl-1H-benzo[d]imidazole core-based potent α-glucosidase inhibitors: Synthesis, kinetic study, molecular docking, and in vivo anti-hyperglycemic evaluation

In the present study, a series of 2-phenyl-1H-benzo[d]imidazole-based α-glucosidase inhibitors were synthesized and evaluated for their in vitro and in vivo anti-diabetic potential. Screening of an in-house library revealed a moderated α-glucosidase inhibitor, 6a with 3-(1H-benzo[d]imidazol-2-yl)aniline core, and then the structural optimization was performed to obtain more efficient derivatives. Most of these derivatives showed increased activity than 6a, and the most promising inhibitors were found to be compounds 15o and 22d with IC₅₀ values of 2.09 ± 0.04 and 0.71 ± 0.02 µM, respectively. Fluorescence quenching experiment confirmed the direct binding of compounds 15o and 22d with α-glucosidase. Kinetic study revealed that both compounds were non-competitive inhibitors, that was consistent with the result of molecular docking studies where they located at the allosteric site of the enzyme. Cell viability evaluation demonstrated the non-cytotoxicity of 15o and 22d against LO2 cells. Furthermore, the in vivo pharmacodynamic study revealed that compound 15o showed significant hypoglycemic activity and improved oral sucrose tolerance, comparable to the positive control acarbose.

Modulating Glycoside Hydrolase Activity between Hydrolysis and Transfer Reactions Using an Evolutionary Approach

The proteins within the CAZy glycoside hydrolase family GH13 catalyze the hydrolysis of polysaccharides such as glycogen and starch. Many of these enzymes also perform transglycosylation in various degrees, ranging from secondary to predominant reactions. Identifying structural determinants associated with GH13 family reaction specificity is key to modifying and designing enzymes with increased specificity towards individual reactions for further applications in industrial, chemical, or biomedical fields. This work proposes a computational approach for decoding the determinant structural composition defining the reaction specificity. This method is based on the conservation of coevolving residues in spatial contacts associated with reaction specificity. To evaluate the algorithm, mutants of α-amylase (TmAmyA) and glucanotransferase (TmGTase) from Thermotoga maritima were constructed to modify the reaction specificity. The K98P/D99A/H222Q variant from TmAmyA doubled the transglycosydation/hydrolysis (T/H) ratio while the M279N variant from TmGTase increased the hydrolysis/transglycosidation ratio five-fold. Molecular dynamic simulations of the variants indicated changes in flexibility that can account for the modified T/H ratio. An essential contribution of the presented computational approach is its capacity to identify residues outside of the active center that affect the reaction specificity.

Human Parainfluenza Virus 3 Phosphoprotein Is a Tetramer and Shares Structural and Interaction Features with Ebola Phosphoprotein VP35

The human parainfluenza virus 3 (HPIV3) poses a risk for pneumonia development in young children and immunocompromised patients. To investigate mechanisms of HPIV3 pathogenesis, we characterized the association state and host protein interactions of HPIV3 phosphoprotein (HPIV3 P), an indispensable viral polymerase cofactor. Sequence analysis and homology modeling predict that HPIV3 P possesses a long, disordered N-terminal tail (PTAIL) a coiled-coil multimerization domain (PMD), similar to the well-characterized paramyxovirus phosphoproteins from measles and Sendai viruses. Using a recombinantly expressed and purified construct of PMD and PTAIL, we show that HPIV3 P in solution is primarily an alpha-helical tetramer that is stable up to 60 °C. Pulldown and isothermal titration calorimetry experiments revealed that HPIV3 P binds the host hub protein LC8, and turbidity experiments demonstrated a new role for LC8 in increasing the solubility of HPIV3 P in the presence of crowding agents such as RNA. For comparison, we show that the multimerization domain of the Zaire Ebola virus phosphoprotein VP35 is also a tetramer and binds LC8 but with significantly higher affinity. Comparative analysis of the domain architecture of various virus phosphoproteins in the order Mononegavirales show multiple predicted and verified LC8 binding motifs, suggesting its prevalence and importance in regulating viral phosphoprotein structures. Our work provides evidence for LC8 binding to phosphoproteins with multiple association states, either tetrameric, as in the HPIV3 and Ebola phosphoproteins shown here, or dimeric as in rabies virus phosphoprotein. Taken together the data suggest that the association states of a virus-specific phosphoprotein and the complex formed by binding of the phosphoprotein to host LC8 are important regulators of viral function.

Identifying carbohydrate-active enzymes of Cutaneotrichosporon oleaginosus using systems biology

Background

The oleaginous yeast Cutaneotrichosporon oleaginosus represents one of the most promising microbial platforms for resource-efficient and scalable lipid production, with the capacity to accept a wide range of carbohydrates encapsulated in complex biomass waste or lignocellulosic hydrolysates. Currently, data related to molecular aspects of the metabolic utilisation of oligomeric carbohydrates are sparse. In addition, comprehensive proteomic information for C. oleaginosus focusing on carbohydrate metabolism is not available.

Results

In this study, we conducted a systematic analysis of carbohydrate intake and utilisation by C. oleaginosus and investigated the influence of different di- and trisaccharide as carbon sources. Changes in the cellular growth and morphology could be observed, depending on the selected carbon source. The greatest changes in morphology were observed in media containing trehalose. A comprehensive proteomic analysis of secreted, cell wall-associated, and cytoplasmatic proteins was performed, which highlighted differences in the composition and quantity of secreted proteins, when grown on different disaccharides. Based on the proteomic data, we performed a relative quantitative analysis of the identified proteins (using glucose as the reference carbon source) and observed carbohydrate-specific protein distributions. When using cellobiose or lactose as the carbon source, we detected three- and five-fold higher diversity in terms of the respective hydrolases released. Furthermore, the analysis of the secreted enzymes enabled identification of the motif with the consensus sequence LALL[LA]L[LA][LA]AAAAAAA as a potential signal peptide.

Conclusions

Relative quantification of spectral intensities from crude proteomic datasets enabled the identification of new enzymes and provided new insights into protein secretion, as well as the molecular mechanisms of carbo-hydrolases involved in the cleavage of the selected carbon oligomers. These insights can help unlock new substrate sources for C. oleaginosus, such as low-cost by-products containing difficult to utilize carbohydrates. In addition, information regarding the carbo-hydrolytic potential of C. oleaginosus facilitates a more precise engineering approach when using targeted genetic approaches. This information could be used to find new and more cost-effective carbon sources for microbial lipid production by the oleaginous yeast C. oleaginosus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-021-01692-2.

Natural Phytochemicals, Luteolin and Isoginkgetin, Inhibit 3C Protease and Infection of FMDV, In Silico and In Vitro

Foot-and-mouth-disease virus (FMDV) is a picornavirus that causes a highly contagious disease of cloven-hoofed animals resulting in economic losses worldwide. The 3C protease (3Cpro) is the main protease essential in the picornavirus life cycle, which is an attractive antiviral target. Here, we used computer-aided virtual screening to filter potential anti-FMDV agents from the natural phytochemical compound libraries. The top 23 filtered compounds were examined for anti-FMDV activities by a cell-based assay, two of which possessed antiviral effects. In the viral and post-viral entry experiments, luteolin and isoginkgetin could significantly block FMDV growth with low 50% effective concentrations (EC50). Moreover, these flavonoids could reduce the viral load as determined by RT-qPCR. However, their prophylactic activities were less effective. Both the cell-based and the fluorescence resonance energy transfer (FRET)-based protease assays confirmed that isoginkgetin was a potent FMDV 3Cpro inhibitor with a 50% inhibition concentration (IC50) of 39.03 ± 0.05 and 65.3 ± 1.7 μM, respectively, whereas luteolin was less effective. Analyses of the protein-ligand interactions revealed that both compounds fit in the substrate-binding pocket and reacted to the key enzymatic residues of the 3Cpro. Our findings suggested that luteolin and isoginkgetin are promising antiviral agents for FMDV and other picornaviruses.

Genetic and Phenotypic Characterization of the Novel Metallo-β-Lactamase NDM-29 From Escherichia coli

Objectives: The New Delhi metallo-β-lactamase (NDM) can hydrolyze almost all clinically available β-lactam antibiotics and has widely spread all over the world. NDM-29, a novel carbapenemase, was discovered in an Escherichia coli (19NC225) isolated from a patient with biliary tract infection in 2019 in China.

Methods: Conjugation, transformation, cloning test, fitness cost, PacBio Sequel, and Illumina sequencing were performed to analyze the genetic and phenotypic characterization of bla_NDM–29.

Results: The susceptibility testing results showed 19NC225 was resistant to cephalosporins, carbapenems, combinations of β-lactam and β-lactamase inhibitors, and levofloxacin. Conjugation and transformation were performed to verify the transferability of NDM-29-encoding plasmid, and cloning test was conducted to prove the function of bla_NDM–29 to increase carbapenem resistance. Furthermore, fitness cost test confirmed that the presence of NDM-29 exerts no survival pressure on bacteria. PacBio Sequel and Illumina sequencing were performed to analyze the genetic characterization of 19NC225, which contains two plasmids (pNC225-TEM1B and pNC225-NDM-29). pNC225-NDM-29, exhibiting 99.96% identity and 100% coverage with pNDM-BTR (an IncN1 plasmid from an E. coli in urine specimen from Beijing in 2013), showed responsibility for the multidrug-resistant (MDR) phenotype. Compared with bla_NDM–1, bla_NDM–29, located on pNC225-NDM-29, carries a G388A (D130N) mutation. The region harboring bla_NDM–29 is located in an ISKpn19-based transposon, and two Tn6292 remnants are symmetrically located upstream and downstream of the transposon. The sequence results also indicated several important virulence genes.

Conclusion: The findings of the novel carbapenemase NDM-29 could pose a threat to the control of antimicrobial resistance and arouse attention about the mutation of bacteria.

Kinetic characterization and thermostability of C. elegans cytoplasmic and mitochondrial malate dehydrogenases

Malate dehydrogenase (MDH) catalyzes the conversion of NAD⁺ and malate to NADH and oxaloacetate in the citric acid cycle. Eukaryotes have one MDH isozyme that is imported into the mitochondria and one in the cytoplasm. We overexpressed and purified Caenorhabditis elegans cytoplasmic MDH-1 and mitochondrial MDH-2 in E. coli. Our goal was to compare the kinetic and structural properties of these enzymes because C. elegans can survive adverse environmental conditions, such as lack of food and elevated temperatures. In steady-state enzyme kinetics assays, we measured K_M values for oxaloacetate of 54 and 52 μM and K_M values for NADH of 61 and 107 μM for MDH-1 and MDH-2, respectively. We partially purified endogenous MDH-1 and MDH-2 from a mixed population of worms and separated them using anion exchange chromatography. Both endogenous enzymes had a K_M for oxaloacetate similar to that of the corresponding recombinant enzyme. Recombinant MDH-1 and MDH-2 had maximum activity at 40 °C and 35 °C, respectively. In a thermotolerance assay, MDH-1 was much more thermostable than MDH-2. Protein homology modeling predicted that MDH-1 had more intersubunit salt-bridges than mammalian MDH1 enzymes, and these ionic interactions may contribute to its thermostability. In contrast, the MDH-2 homology model predicted fewer intersubunit ionic interactions compared to mammalian MDH2 enzymes. These results suggest that the increased stability of MDH-1 may facilitate its ability to remain active in adverse environmental conditions. In contrast, MDH-2 may use other strategies, such as protein binding partners, to function under similar conditions.

Mutational spectrum of BRAF gene in colorectal cancer patients in Saudi Arabia

Colorectal cancer (CRC) is one of the topmost causes of death in males in Saudi Arabia. In females, it was also within the top five cancer types. CRC is heterogeneous in terms of pathogenicity and molecular genetic pathways. It is very important to determine the genetic causes of CRC in the Saudi population. BRAF is one of the major genes involved in cancers, it participates in transmitting chemical signals from outside the cells into the nucleus of the cells and it is also shown to participate in cell growth. In this study, we mapped the spectrum of BRAF mutations in 100 Saudi patients with CRC. We collected tissue samples from colorectal cancer patients, sequenced the BRAF gene to identify gene alterations, and analyzed the data using different bioinformatics tools. We designed a three-dimensional (3D) homology model of the BRAF protein using the Swiss Model automated homology modeling platform to study the structural impact of these mutations using the Missense3D algorithm. We found six mutations in 14 patients with CRC. Four of these mutations are being reported for the first time. The novel frameshift mutations observed in CRC patients, such as c.1758delA (E586E), c.1826insT (Q609L), c.1860insA and c.1860insA/C (M620I), led to truncated proteins of 589, 610, and 629 amino acids, respectively, and potentially affected the structure and the normal functions of BRAF. These findings provide insights into the molecular etiology of CRC in general and to the Saudi population. BRAF genetic testing may also guide treatment modalities, and the treatment may be optimized based on personalized gene variations.

Molecular and genetic characterization of a large Brazilian cohort presenting hearing loss

Hearing loss is one of the most common sensory defects, affecting 5.5% of the worldwide population and significantly impacting health and social life. It is mainly attributed to genetic causes, but their relative contribution reflects the geographical region's socio-economic development. Extreme genetic heterogeneity with hundreds of deafness genes involved poses challenges for molecular diagnosis. Here we report the investigation of 542 hearing-impaired subjects from all Brazilian regions to search for genetic causes. Biallelic GJB2/GJB6 causative variants were identified in 12.9% (the lowest frequency was found in the Northern region, 7.7%), 0.4% carried GJB2 dominant variants, and 0.6% had the m.1555A > G variant (one aminoglycoside-related). In addition, other genetic screenings, employed in selected probands according to clinical presentation and presumptive inheritance patterns, identified causative variants in 2.4%. Ear malformations and auditory neuropathy were diagnosed in 10.8% and 3.5% of probands, respectively. In 3.8% of prelingual/perilingual cases, Waardenburg syndrome was clinically diagnosed, and in 71.4%, these diagnoses were confirmed with pathogenic variants revealed; seven out of them were novel, including one CNV. All these genetic screening strategies revealed causative variants in 16.2% of the cases. Based on causative variants in the molecular diagnosis and genealogy analyses, a probable genetic etiology was found in ~ 50% of the cases. The present study highlights the relevance of GJB2/GJB6 as a cause of hearing loss in all Brazilian regions and the importance of screening unselected samples for estimating frequencies. Moreover, when a comprehensive screening is not available, molecular diagnosis can be enhanced by selecting probands for specific screenings.

Characterization of Inducible HSP70 Genes in an Antarctic Yeast, Glaciozyma antarctica PI12, in Response to Thermal Stress

The induction of highly conserved heat shock protein 70 (HSP70) is often related to a cellular response due to harmful stress or adverse life conditions. In this study, we determined the expression of Hsp70 genes in the Antarctic yeast, Glaciozyma antarctica, under different several thermal treatments for several exposure periods. The main aims of the present study were (1) to determine if stress-induced Hsp70 could be used to monitor the exposure of the yeast species G. antarctica to various types of thermal stress; (2) to analyze the structures of the G. antarctica HSP70 proteins using comparative modeling; and (3) to evaluate the relationship between the function and structure of HSP70 in G. antarctica. In this study, we managed to amplify and clone 2 Hsp70 genes from G. antarctica named GaHsp70-1 and GaHsp70-2. The cells of G. antarctica expressed significantly inducible Hsp70 genes after the heat and cold shock treatments. Interestingly, GaHsp70-1 showed 2–6-fold higher expression than GaHsp70-2 after the heat and cold exposure. ATP hydrolysis analysis on both G. antarctica HSP70s proved that these psychrophilic chaperones can perform activities in a wide range of temperatures, such as at 37, 25, 15, and 4 °C. The 3D structures of both HSP70s revealed several interesting findings, such as the substitution of a β-sheet to loop in the N-terminal ATPase binding domain and some modest residue substitutions, which gave the proteins the flexibility to function at low temperatures and retain their functional activity at ambient temperatures. In conclusion, both analyzed HSP70s played important roles in the physiological adaptation of G. antarctica.

Potent Inhibition of HIF1α and p300 Interaction by a Constrained Peptide Derived from CITED2

Disrupting the interaction between HIF1α and p300 is a promising strategy to modulate the hypoxia response of tumor cells. Herein, we designed a constrained peptide inhibitor derived from the CITED2/p300 complex to disturb the HIF1α/p300 interaction. Through truncation/mutation screening and a terminal aspartic acid-stabilized strategy, a constrained peptide was constructed with outstanding biochemical/biophysical properties, especially in binding affinity, cell penetration, and serum stability. To date, our study was the first one to showcase that stabilized peptides derived from CITED2 using helix-stabilizing methods acted as a promising candidate for modulating hypoxia-inducible signaling.

Identification and Characterization of a Rhodopsin Kinase Gene in the Suckers of Octopus vulgaris: Looking around Using Arms?

Simple Summary

Octopus arms are a fascinating and evolutionarily unique sensory organ, with hundreds of motile suckers, each with thousands of sensory cells, lining eight highly flexible arms. Scientifically, there are many open questions regarding the sensory capabilities of the arms and specifically the highly innervated suckers. In our present work, we used a multidisciplinary approach to fully characterize the light-sensing molecule, Ov-GRK1, in the suckers, skin and retina of Octopus vulgaris. We sequenced the O. vulgaris GRK1 gene, defining a phylogenetic tree and performing a 3D structure model prediction. We found differences in the relative expression of mRNA in different sucker types at several locations along the arm, which might indicate a functional difference. Using labeling methods, we localized the expression to the highly sensitive sucker rim. Our findings indicate that octopus suckers, in specific areas of the arm, might have the ability for light sensing. We therefore suggest that suckers are tactile, chemical and light sensors.

Abstract

In their foraging behavior octopuses rely on arm search movements outside the visual field of the eyes. In these movements the environment is explored primarily by the suckers that line the entire length of the octopus arm. In this study, for the first time, we report the complete characterization of a light-sensing molecule, Ov-GRK1, in the suckers, skin and retina of Octopus vulgaris. We sequenced the O. vulgaris GRK1 gene, defining a phylogenetic tree and performing a 3D structure model prediction. Furthermore, we found differences in relative mRNA expression in different sucker types at several arm levels, and localized it through in situ hybridization. Our findings suggest that the suckers in octopus arms are much more multimodal than was previously shown, adding the potential for light sensing to the already known mechanical and chemical sensing abilities.

Identification and Characterization of an Antennae-Specific Glutathione S-Transferase From the Indian Meal Moth

Insect glutathione-S-transferases (GSTs) play essential roles in metabolizing endogenous and exogenous compounds. GSTs that are uniquely expressed in antennae are assumed to function as scavengers of pheromones and host volatiles in the odorant detection system. Based on this assumption, antennae-specific GSTs have been identified and functionally characterized in increasing number of insect species. In the present study, 17 putative GSTs were identified from the antennal transcriptomic dataset of the Indian meal moth, Plodia interpunctella, a severe stored-grain pest worldwide. Among the GSTs, only PiGSTd1 is antennae-specific according to both Fragments Per Kilobase Million (FPKM) and quantitative real-time PCR (qRT-PCR) analysis. Sequence analysis revealed that PiGSTd1 has a similar identity as many delta GSTs from other moths. Enzyme kinetic assays using 1-chloro-2,4-dinitrobenzene (CDNB) as substrates showed that the recombinant PiGSTd1 gave a K_m of 0.2292 ± 0.01805 mM and a V_max of 14.02 ± 0.2545 μmol·mg⁻¹·min⁻¹ under the optimal catalytic conditions (35°C and pH = 7.5). Further analysis revealed that the recombinant PiGSTd1 could efficiently degrade the sex pheromone component Z9-12:Ac (75.63 ± 5.52%), as well as aldehyde volatiles, including hexanal (89.10 ± 2.21%), heptanal (63.19 ± 5.36%), (E)-2-octenal (73.58 ± 3.92%), (E)-2-nonenal (75.81 ± 1.90%), and (E)-2-decenal (61.13 ± 5.24%). Taken together, our findings suggest that PiGSTd1 may play essential roles in degrading and inactivating a variety of odorants, especially sex pheromones and host volatiles of P. interpunctella.

Modeling human glucose-6-phosphate dehydrogenase mutations using C. elegans GSPD-1

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is an X-linked, recessive condition that causes intermittent jaundice or hemolytic anemia because of low NADPH levels in red blood cells. We performed steady-state enzyme kinetics with the recombinant C. elegans ortholog of human G6PD, GSPD-1, and two mutants containing amino acid changes found in human patients. The K M values for glucose-6-phosphate were 100 ± 27 µM, 80 ± 22 µM, and 1000 ± 300 µM for the wild-type, D60N, and R252L GSPD-1 enzymes, respectively. The specific activities of the D60N and R252L mutants were 59% and 11%, respectively, of the wild-type value. Protein homology modeling suggested that the R252L mutation was more severe because the mutation caused a shift in the position of some active site residues. The D60N mutation may have affected the conformation of an outer loop of the enzyme. These data demonstrate that GSPD-1 is a promising model for human G6PD deficiencies, with the advantage that potential treatments could be studied in vivo in C. elegans.

Structural basis for membrane recruitment of ATG16L1 by WIPI2 in autophagy

Autophagy is a cellular process that degrades cytoplasmic cargo by engulfing it in a double-membrane vesicle, known as the autophagosome, and delivering it to the lysosome. The ATG12-5-16L1 complex is responsible for conjugating members of the ubiquitin-like ATG8 protein family to phosphatidylethanolamine in the growing autophagosomal membrane, known as the phagophore. ATG12-5-16L1 is recruited to the phagophore by a subset of the phosphatidylinositol 3-phosphate-binding seven-bladedß -propeller WIPI proteins. We determined the crystal structure of WIPI2d in complex with the WIPI2 interacting region (W2IR) of ATG16L1 comprising residues 207-230 at 1.85 Å resolution. The structure shows that the ATG16L1 W2IR adopts an alpha helical conformation and binds in an electropositive and hydrophobic groove between WIPI2 ß-propeller blades 2 and 3. Mutation of residues at the interface reduces or blocks the recruitment of ATG12-5-16 L1 and the conjugation of the ATG8 protein LC3B to synthetic membranes. Interface mutants show a decrease in starvation-induced autophagy. Comparisons across the four human WIPIs suggest that WIPI1 and 2 belong to a W2IR-binding subclass responsible for localizing ATG12-5-16 L1 and driving ATG8 lipidation, whilst WIPI3 and 4 belong to a second W34IR-binding subclass responsible for localizing ATG2, and so directing lipid supply to the nascent phagophore. The structure provides a framework for understanding the regulatory node connecting two central events in autophagy initiation, the action of the autophagic PI 3-kinase complex on the one hand and ATG8 lipidation on the other.

Revealing the molecular interplay of curcumin as Culex pipiens Acetylcholine esterase 1 (AChE1) inhibitor

Emergence of vector borne diseases has continued to take toll on millions of lives since its inception. The use of insecticides began as vector control strategy in the early 1900's but the menace of insects is still prevalent. Additionally, the inadequate use of organophosphates and carbamates which target acetylcholine esterase (AChE), are known to develop resistance amongst vectors of transmission and are toxic to humans. In this study, extensive computational screening was performed using homology modelling, molecular docking, molecular dynamics (MD) simulation and free energy change calculation, which highlighted curcumin as a lead molecule out of ~ 1700 phytochemicals against Culex pipiens AChE. In vivo larvicidal activity was carried out along with in vivo and in vitro AChE inhibition assay to determine the biochemical efficacy of curcumin. Our study reveals that curcumin induces mortality in Cx. pipiens at an early stage of its life cycle by AChE inhibition. This also underlines the use of curcumin as a coming-age natural product insecticide.

Thymidylate kinase (TMK) of the photosynthetic, nitrogen-fixing cyanobacterium Nostoc sp. strain PCC7120: Biophysical, biochemical and physiological characterisation

Thymidylate kinase (TMK/TMPK) is an important enzyme in DNA biosynthesis and catalyses the conversion of dTMP to dTDP. Due to its therapeutic potential, the focus has been on characterizing the TMK proteins of pathogens and human origin, with very little information available on the TMK proteins of photosynthetic organisms and agriculturally important nitrogen-fixing organisms. In this work we report the characterisation of TMK in an evolutionarily ancient organism, cyanobacteria. The TMK protein of the photosynthetic, nitrogen-fixing cyanobacterium Nostoc PCC7120 (AnTMK) was found to have low conformational stability, which related to its low Tm of ~46 °C confirmed by Differential Scanning Fluorimetry (DSF) and Differential Scanning Calorimetry (DSC) techniques. The AnTMK protein exhibited substrate specificity for dTMP and ATP with K_m of 20.74 ± 1.47 μM and 20.17 ± 2.96 μM respectively. The enzyme kinetics data and the positive co-operativity observed between dTMP and ATP binding correlated well with the data obtained from Isothermal Titration Calorimetry (ITC). Homology model of the enzyme suggested that the binding mode of substrate nucleotides to the enzyme is conserved. When overexpressed constitutively in Nostoc PCC7120 (Antmk⁺), it supported faster growth measured in terms of chlorophyll a content under normal growth conditions, but exhibited lower photosynthetic efficiency. Compared to the vector control recombinant Nostoc AnpAM, the Antmk ⁺ cells exhibited higher photoinhibition at higher light irradiance with more open reaction centres and lower dissipation of heat, indicative of damage to photosynthetic machinery. This indicated that the TMK is likely to have a significant role in photosynthetic organisms.

Effects of temperature and pH on the structure of a metalloprotease from Lactobacillus fermentum R6 isolated from Harbin dry sausages and molecular docking between protease and meat protein

BACKGROUND

Microbial protease can interact with meat protein in fermented meat products at a certain pH and temperature. To investigate the effects of various pH values and temperatures on the structural characteristics of Lactobacillus fermentum R6 protease, which was isolated from Harbin dry sausages, spectroscopy techniques and molecular dynamics were utilized to evaluate structural changes.

RESULTS

The protease exhibited a stable spatial structure at pH 7 and 40 °C, and the extension of the protease structure was also promoted. Although the structure of the protease could be changed or destroyed by pH 8 and 70 °C, it was mainly determined by the changes of secondary and tertiary structures such as α-helix, β-sheet, β-turn and random coil. In addition, carbonyl vibration, -NH vibration, C-H stretching vibration and disulphide bonds were present in L. fermentum R6 protease under various pH and temperature conditions. Molecular docking showed that the protease can interact with myosin light chain, myosin heavy chain, actin and myoglobin.

CONCLUSION

The protease can maintain stable structure and interact with meat protein, which reflected certain application prospects in the fermentation of Harbin dry sausages. © 2021 Society of Chemical Industry.

A new P450 involved in the furanocoumarin pathway underlies a recent case of convergent evolution

Furanocoumarins are phytoalexins often cited as an example to illustrate the arms race between plants and herbivorous insects. They are distributed in a limited number of phylogenetically distant plant lineages, but synthesized through a similar pathway, which raised the question of a unique or multiple emergence in higher plants. The furanocoumarin pathway was investigated in the fig tree (Ficus carica, Moraceae). Transcriptomic and metabolomic approaches led to the identification of CYP76F112, a cytochrome P450 catalyzing an original reaction. CYP76F112 emergence was inquired using phylogenetics combined with in silico modeling and site-directed mutagenesis. CYP76F112 was found to convert demethylsuberosin into marmesin with a very high affinity. This atypical cyclization reaction represents a key step within the polyphenol biosynthesis pathway. CYP76F112 evolutionary patterns suggests that the marmesin synthase activity appeared recently in the Moraceae family, through a lineage-specific expansion and diversification. The characterization of CYP76F112 as the first known marmesin synthase opens new prospects for the use of the furanocoumarin pathway. It also supports the multiple acquisition of furanocoumarin in angiosperms by convergent evolution, and opens new perspectives regarding the ability of cytochromes P450 to evolve new functions related to plant adaptation to their environment.

Determination of the interaction between the receptor binding domain of 2019-nCoV spike protein, TMPRSS2, cathepsin B and cathepsin L, and glycosidic and aglycon forms of some flavonols

The novel coronavirus (COVID-19, SARS-CoV-2) is a rapidly spreading disease with a high mortality. In this research, the interactions between specific flavonols and the 2019-nCoV receptor binding domain (RBD), transmembrane protease, serine 2 (TMPRSS2), and cathepsins (CatB and CatL) were analyzed. According to the relative binding capacity index (RBCI) calculated based on the free energy of binding and calculated inhibition constants, it was determined that robinin (ROB) and gossypetin (GOS) were the most effective flavonols on all targets. While the binding free energy of ROB with the spike glycoprotein RBD, TMPRSS2, CatB, and CatL were -5.02, -7.57, -10.10, and -6.11 kcal/mol, the values for GOS were -4.67, -5.24, -8.31, and -6.76, respectively. Furthermore, both compounds maintained their stability for at least 170 ns on respective targets in molecular dynamics simulations. The molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) calculations also corroborated these data. Considering Lipinski's rule of five, ROB and GOS exhibited 3 (MW>500, N or O>10, NH or OH>5), and 1 (NH or OH>5) violations, respectively. Neither ROB nor GOS showed AMES toxicity or hepatotoxicity. The LD50 of these compounds in rats were 2.482 and 2.527 mol/kg, respectively. Therefore, we conclude that these compounds could be considered as alternative therapeutic agents in the treatment of COVID-19. However, the possible inhibitory effects of these compounds on cytochromes (CYPs) should be verified by in vitro or in vivo tests and their adverse effects on cellular energy metabolism should be minimized by performing molecular modifications if necessary.

Potent Human Single-Domain Antibodies Specific for a Novel Prefusion Epitope of Respiratory Syncytial Virus F Glycoprotein

Respiratory syncytial virus (RSV) poses great health threats to humans. However, there are no licensed vaccines or therapeutic drugs to date. Only one humanized monoclonal antibody, palivizumab, is available on the market, but it is used prophylactically and is limited to infants with high risk. With advances in antibody engineering, it has been found that a single-domain antibody (sdAb) can be therapeutically administered by inhalation, which would be more efficient for respiratory diseases. Here, we identified two human sdAbs, m17 and m35, by phage display technology. They specifically bind to RSV fusion glycoprotein (F protein) in the prefusion state with subnanomolar affinity and potently neutralize both RSV subtypes A and B with 50% inhibitory concentration (IC₅₀) values ranging from pM to nM. Interestingly, these sdAbs recognize a novel epitope, termed VI, that is unique to the prefusion state. This epitope is located at the C terminus of the F1 subunit, close to the viral membrane, and might be sterically restricted. We further find that m17 and m35 neutralize RSV by preventing the prefusion F conformational arrangement, thus inhibiting membrane fusion. These two sdAbs have the potential to be further developed as therapeutic candidates and may also provide novel insight for developing other antiviral reagents against RSV. IMPORTANCE Because respiratory syncytial virus (RSV) can cause serious respiratory disease in immunodeficient groups, including infants and seniors, the development of vaccines and therapeutic drugs, such as neutralizing antibodies, is urgently needed. Compared to the conventional full-length antibody, a single-domain antibody (sdAb) has been demonstrated to be efficient for respiratory diseases when administered by inhalation, thereby potentially introducing a kind of novel therapeutic agent in the market. Here, we discovered two potent neutralizing human sdAbs against RSV that recognized a novel prefusion epitope, termed VI, and prevented conformational arrangement during the fusion process. Our work provides not only therapeutic candidates but also novel targets for new drug and vaccine development.

Drugs repurposing against SARS-CoV2 and the new variant B.1.1.7 (alpha strain) targeting the spike protein: molecular docking and simulation studies

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV2) is responsible for the global COVID-19 pandemic and millions of deaths worldwide. In December 2020, a new alpha strain of SARS-CoV2 was identified in the United Kingdom. It was referred to as VUI 202012/01 (Alpha strain modelled under investigation, 2020, month 12, number 01). The interaction between spike protein and ACE2 receptor is a prerequisite for entering virion into the host cell. The present study is focussed on the spike protein of the SARS-COV 2, involving the comparison of binding affinity of new alpha strain modelled spike with previous strain spike (PDB ID:7DDN) using in silico molecular docking, dynamics and simulation studies. The molecular docking studies of the alpha strain modelled spike protein confirmed its higher affinity for the ACE2 receptor than the spike protein of the dominant strain. Similar computational approaches have also been used to investigate the potency of FDA approved drugs from the ZINC Database against the spike protein of new alpha strain modelled and old ones. The drug molecules which showed strong affinity for both the spike proteins are then subjected to ADME analysis. The overall binding energy of Conivaptan (-107.503 kJ/mol) and Trosec (-94.029 kJ/mol) is indicative of their strong binding affinities, well supported by interactions with critical residues.

•

We investigated the potential FDA drugs for repurposing against the spike protein of alpha strain modelled of SARS-CoV-2.

•

Spike protein of alpha strain modelled of SARS-CoV-2 has more affinity for the ACE2 receptor of host cell than previous strains and, therefore, is more contagious.

•

Conivaptan, Ecamsule and Trosec are common drugs that bind strongly with spike protein of both the strains .

•

Molecular Docking and simulation studies show that Conivaptan and Trosec emerged as potential inhibitors of the alpha strain modelled spike protein.

Fucosidosis in Tunisian patients: mutational analysis and homology-based modeling of FUCA1 enzyme

BACKGROUND

Fucosidosis is an autosomal recessive lysosomal storage disease caused by defective alpha-L-fucosidase (FUCA1) activity, leading to the accumulation of fucose-containing glycolipids and glycoproteins in various tissues. Clinical features include angiokeratoma, progressive psychomotor retardation, neurologic signs, coarse facial features, and dysostosis multiplex.

METHODS

All exons and flanking intron regions of FUCA1 were screened by direct sequencing to identify mutations and polymorphisms in three unrelated families with fucosidosis. Bioinformatics tools were then used to predict the impacts of novel alterations on the structure and function of proteins. Furthermore, the identified mutations were localized onto a 3D structure model using the DeepView Swiss-PdbViewer 4.1 software, which established a function-structure relationship of the FUCA1 proteins.

RESULTS

Four novel mutations were identified in this study. Two patients (P1 and P2) in Families 1 and 2 who had the severe phenotype were homoallelic for the two identified frameshift mutations p.K57Sfs*75 and p.F77Sfs*55, respectively. The affected patient (P3) from Family 3, who had the milder phenotype, was heterozygous for the novel missense mutation p.G332E and the novel splice site mutation c.662+5g>c. We verified that this sequence variation did not correspond to a polymorphism by testing 50 unrelated individuals. Additionally, 16 FUCA1 polymorphisms were identified. The structure prediction analysis showed that the missense mutation p.G332E would probably lead to a significant conformational change, thereby preventing the expression of the FUCA1 protein indeed; the 3D structural model of the FUCA1 protein reveals that the glycine at position 332 is located near a catalytic nucleophilic residue. This makes it likely that the enzymatic function of the protein with p.G332E is severely impaired.

CONCLUSION

These are the first FUCA1 mutations identified in Tunisia that cause the fucosidosis disease. Bioinformatics analysis allowed us to establish an approximate structure-function relationship for the FUCA1 protein, thereby providing better genotype/phenotype correlation knowledge.

Novel tetrahydrobenzo[b]thiophen-2-yl)urea derivatives as novel α-glucosidase inhibitors: Synthesis, kinetics study, molecular docking, and in vivo anti-hyperglycemic evaluation

α-Glucosidase inhibitors, which can inhibit the digestion of carbohydrates into glucose, are one of important groups of anti-type 2 diabetic drugs. In the present study, we report our effort on the discovery and optimization of α-glucosidase inhibitors with tetrahydrobenzo[b]thiophen-2-yl)urea core. Screening of an in-house library revealed a moderated α-glucosidase inhibitors, 5a, and then the following structural optimization was performed to obtain more efficient derivatives. Most of these derivatives showed increased inhibitory activity against α-glucosidase than the parental compound 5a (IC₅₀ of 26.71 ± 1.80 μM) and the positive control acarbose (IC₅₀ of 258.53 ± 1.27 μM). Among them, compounds 8r (IC₅₀ = 0.59 ± 0.02 μM) and 8s (IC₅₀ = 0.65 ± 0.03 μM) were the most potent inhibitors, and showed selectivity over α-amylase. The direct binding of both compounds with α-glucosidase was confirmed by fluorescence quenching experiments. Kinetics study revealed that these compounds were non-competitive inhibitors, which was consistent with the molecular docking results that compounds 8r and 8s showed high preference to bind to the allosteric site instead of the active site of α-glucosidase. In addition, compounds 8r and 8s were not toxic (IC₅₀ > 100 μM) towards LO2 and HepG2 cells. Finally, the in vivo anti-hyperglycaemic activity assay results indicated that compounds 8r could significantly decrease the level of plasma glucose and improve glucose tolerance in SD rats treated with sucrose. The present study provided the tetrahydrobenzo[b]thiophen-2-yl)urea chemotype for developing novel α-glucosidase inhibitors against type 2 diabetes.

Discovery of a new light-driven Li+/Na+-pumping rhodopsin with DTG motif

Microbial pumping rhodopsin is a seven-transmembrane retinal binding protein, which is light-driven ion pump with a functional key motif. Ion-pumping with the key motif and charged amino acids in the rhodopsin is biochemically important. The rhodopsins with DTG motif have been discovered in various eubacteria, and they function as H⁺ pump. Especially, the DTG motif rhodopsins transported H⁺ despite the replacement of a proton donor by Gly. We investigated Methylobacterium populi rhodopsin (MpR) in one of the DTG motif rhodopsin clades. To determine which ions the MpR transport, we tested with various monovalent ion solutions and determined that MpR transports Li⁺/Na⁺. By replacing the three negatively charged residues residues which are located in helix B, Glu32, Glu33, and Asp35, we concluded that the residues play a critical role in the transport of Li⁺/Na⁺. The MpR E33Q transported H⁺ in place of Li⁺/Na⁺, suggesting that Glu33 is a Li⁺/Na⁺ binding site on the cytoplasmic side. Gly93 in MpR was replaced by Asp to convert from the Li⁺/Na⁺ pump to the H⁺ pump, resulting in MpR G93D transporting H⁺. Dissociation constant (K_d) values of Na⁺ for MpR WT and E33Q were determined to be 4.0 and 72.5 mM, respectively. These results indicated the mechanism by which MpR E33Q transports H⁺. Up to now, various ion-pumping rhodopsins have been discovered, and Li⁺/Na⁺-pumping rhodopsins were only found in the NDQ motif in NaR. Here, we report a new light-driven Na⁺ pump MpR and have determined the important residues required for Li⁺/Na⁺-pumping different from previously known NaR.

Docking-generated multiple ligand poses for bootstrapping bioactivity classifying Machine Learning: Repurposing covalent inhibitors for COVID-19-related TMPRSS2 as case study

In the present work we introduce the use of multiple docked poses for bootstrapping machine learning-based QSAR modelling. Ligand-receptor contact fingerprints are implemented as descriptor variables. We implemented this method for the discovery of potential inhibitors of the serine protease enzyme TMPRSS2 involved the infectivity of coronaviruses. Several machine learners were scanned, however, Xgboost, support vector machines (SVM) and random forests (RF) were the best with testing set accuracies reaching 90%. Three potential hits were identified upon using the method to scan known untested FDA approved drugs against TMPRSS2. Subsequent molecular dynamics simulation and covalent docking supported the results of the new computational approach.

Characterization of Dynamic Regulatory Gene and Protein Networks in Wheat Roots Upon Perceiving Water Deficit Through Comparative Transcriptomics Survey

A well-developed root system benefits host plants by optimizing water absorption and nutrient uptake and thereby increases plant productivity. In this study we have characterized the root transcriptome using RNA-seq and subsequential functional analysis in a set of drought tolerant and susceptible genotypes. The goal of the study was to elucidate and characterize water deficit-responsive genes in wheat landraces that had been through long-term field and biochemical screening for drought tolerance. The results confirm genotype differences in water-deficit tolerance in line with earlier results from field trials. The transcriptomics survey highlighted a total of 14,187 differentially expressed genes (DEGs) that responded to water deficit. The characterization of these genes shows that all chromosomes contribute to water-deficit tolerance, but to different degrees, and the B genome showed higher involvement than the A and D genomes. The DEGs were mainly mapped to flavonoid, phenylpropanoid, and diterpenoid biosynthesis pathways, as well as glutathione metabolism and hormone signaling. Furthermore, extracellular region, apoplast, cell periphery, and external encapsulating structure were the main water deficit-responsive cellular components in roots. A total of 1,377 DEGs were also predicted to function as transcription factors (TFs) from different families regulating downstream cascades. TFs from the AP2/ERF-ERF, MYB-related, B3, WRKY, Tify, and NAC families were the main genotype-specific regulatory factors. To further characterize the dynamic biosynthetic pathways, protein-protein interaction (PPI) networks were constructed using significant KEGG proteins and putative TFs. In PPIs, enzymes from the CYP450, TaABA8OH2, PAL, and GST families play important roles in water-deficit tolerance in connection with MYB13-1, MADS-box, and NAC transcription factors.

In Silico Study of the Structure and Ligand Preference of Pyruvate Kinases from Cyanobacterium Synechocystis sp. PCC 6803

Finding reliable cheap sources for producing chemicals and materials is always challenging. During recent decades, photosynthetic organisms such as cyanobacteria, which used CO₂ as a carbon source for making products, have attracted a great deal of attention. Among cyanobacteria, Synechocystis sp. PCC 6803 has been considered as a model strain and has some desirable features that make it suitable for use as an industrial strain. Pyruvate kinase (PK) catalyzes the transformation of phosphoenolpyruvate (PEP) to pyruvate in the last step of glycolysis that is an essential enzyme to produce adenosine triphosphate (ATP) in all organisms. Therefore, it plays a critical role in regulating cell metabolism. However, active and allosteric sites of PK and allosteric mechanisms governing PK activity are poorly understood in many bacteria. This study was aimed to provide more insight into PKs of Synechocystis sp. PCC 6803, using in silico methods. The results indicated that predicted structures of PKs from Synechocystis sp. PCC 6803 are reliable and can be considered for further studies. Molecular docking studies suggested that for predicted structures of sll0587 and sll1275, respectively, there are three and two possible active or allosteric sites. Furthermore, molecular interaction analysis of modeled structures proposes that sll0587 is strongly inhibited by ATP and when ATP concentration is low, this isoenzyme is active.

Bioinformatic prediction and identification of immunogenic epitopes of the antigenic 14-3-3 protein of Echinococcus multilocularis

INTRODUCTION

Alveolar echinococcosis is a high-risk parasitic disease caused by the larval stage of Echinococcus multilocularis. The study aimed to predict and identify the dominant Th1/Th2 and B cell epitopes of the antigen protein 14-3-3 beta:alpha from Echinococcus multilocularis.

METHODS

A comparison of the four amino acid sequences of 14-3-3 beta:alpha was respectively derived from Echinococcus multilocularis, Rattus norvegicus, Canis lupus familiaris, and Homo sapiens was carried out by CLUSTALW to provide a basis for excluding similar epitopes. The amino acid sequence information was analyzed by SOPMA and the homology model was established by Swiss-Model. IEDB and SYFPEITHI were used to predict T cell epitopes. According to the Bcepred and ABCpred, the B cell epitopes were comprehensively predicted and analyzed. The dominant epitopes were validated by Lymphocyte Proliferation, ELISA, ELISpot, and Flow cytometry.

RESULTS

Eight potential epitopes of 14-3-3 from Echinococcus multilocularis were screened according to the results of prediction and analysis: 14-3-3_1-15, 14-3-3_6-21, 14-3-3_71-86, 14-3-3_144-157, 14-3-3_145-166, 14-3-3_146-160, 14-3-3_153-161, and 14-3-3_164-177. The 3D structure model of the protein was constructed and the location distribution of potential epitope was ascertained. Respectively, the epitopes of the dominant antigen of B cells were validated as 14-3-3_145-166 and 14-3-3_164-177; the Th1 dominant antigen epitopes were 14-3-3_6-21, 14-3-3_145-166; and the Th2 dominant epitopes was 14-3-3_145-166.

CONCLUSION

In this study, two dominant antigen epitopes of B cells, two Th1 dominant antigen epitopes, and one Th2 dominant antigen epitope were validated. Our work provides a basis for the subsequent development of efficient and safe vaccines targeting epitopes of Echinococcus multilocularis.

Effects of changes in glycan composition on glycoprotein dynamics: example of N-glycans on insulin receptor

Glycosylation is among the most common post-translational modifications in proteins, although it is observed in only about 10% of all the protein structures in protein data bank (PDB). Modifications of sugar composition in glycoproteins profoundly impact the overall physiology of the organism. One such example is the development of insulin resistance, which has been attributed to the removal of sialic acid residues from N-glycans of insulin receptor (IR) from various experimental studies. How such modifications affect the glycan-glycoprotein dynamics, and ultimately their function is not clearly understood to date. In this study, we performed molecular dynamics simulations of glycans in different environments. We studied the effects of removal of sialic acid on the glycan, as well as on the dynamics of leucine-rich repeat L1 domain of the IR ectodomain. We observed perturbations in L1 domain dynamics as a result of the removal of sialic acid. The perturbations include an increase in the flexibility of insulin-binding residues, which may affect insulin binding with IR. These changes are accompanied by perturbations in glycan-protein interactions and perturbation of long-range allosteric dynamics. Our observations will further aid in understanding the role of sugars in maintaining homeostasis and how changes in glycan composition may lead to perturbations in homeostasis, ultimately leading to conditions such as insulin resistance.

An astonishing wealth of new proteasome homologs

Motivation

The proteasome is the main proteolytic machine for targeted protein degradation in archaea and eukaryotes. While some bacteria also possess the proteasome, most of them contain a simpler and more specialized homolog, the heat shock locus V protease. In recent years, three further homologs of the proteasome core subunits have been characterized in prokaryotes: Anbu, BPH and connectase. With the inclusion of these members, the family of proteasome-like proteins now exhibits a range of architectural and functional forms, from the canonical proteasome, a barrel-shaped protease without pronounced intrinsic substrate specificity, to the monomeric connectase, a highly specific protein ligase.

Results

We employed systematic sequence searches to show that we have only seen the tip of the iceberg so far and that beyond the hitherto known proteasome homologs lies a wealth of distantly related, uncharacterized homologs. We describe a total of 22 novel proteasome homologs in bacteria and archaea. Using sequence and structure analysis, we analyze their evolutionary history and assess structural differences that may modulate their function. With this initial description, we aim to stimulate the experimental investigation of these novel proteasome-like family members.

Availability and implementation

The protein sequences in this study are searchable in the MPI Bioinformatics Toolkit (https://toolkit.tuebingen.mpg.de) with ProtBLAST/PSI-BLAST and with HHpred (database ‘proteasome_homologs’). The following data are available at https://data.mendeley.com/datasets/t48yhff7hs/3: (i) sequence alignments for each proteasome-like homolog, (ii) the coordinates for their structural models and (iii) a cluster-map file, which can be navigated interactively in CLANS and gives direct access to all the sequences in this study.

Supplementary information

Supplementary data are available at Bioinformatics online.

Binding affinity and mechanisms of SARS-CoV-2 variants

•

Theoretical investigation of the binding affinity of the SARS-CoV-2 variants.

•

The 501Y.V1 variant yielded the highest enhancements in binding affinity.

•

The binding affinity mainly enhanced by forming hydrogen bonds and salt bridges.

During the rapid worldwide spread of SARS-CoV-2, the viral genome has been undergoing numerous mutations, especially in the spike (S) glycoprotein gene that encode a type-I fusion protein, which plays an important role in the infectivity and transmissibility of the virus into the host cell. In this work, we studied the effect of S glycoprotein residue mutations on the binding affinity and mechanisms of SARS-CoV-2 using molecular dynamics simulations and sequence analysis. We quantitatively determined the degrees of binding affinity caused by different S glycoprotein mutations, and the result indicated that the 501Y.V1 variant yielded the highest enhancements in binding affinity (increased by 36.8%), followed by the N439K variant (increased by 29.5%) and the 501Y.V2 variant (increased by 19.6%). We further studied the structures, chemical bonds, binding free energies (enthalpy and entropy), and residue contribution decompositions of these variants to provide physical explanations for the changes in SARS-CoV-2 binding affinity caused by these residue mutations. This research identified the binding affinity differences of the SARS-CoV-2 variants and provides a basis for further surveillance, diagnosis, and evaluation of mutated viruses.