Modeling protein quaternary structure of homo- and hetero-oligomers beyond binary interactions by homology

Molecular dynamics of SARS-CoV-2 omicron variants from Philippine isolates against hesperidin as spike protein inhibitor

SARS-CoV-2 remains a global threat with new sublineages posing challenges, particularly in the Philippines. Hesperidin (HD) is being studied as a potential prophylactic for COVID-19. However, the virus's rapid evolution could alter how HD binds to it, affecting its effectiveness. Here, we study the mutation-induced variabilities of HD dynamics and their effects on molecular energetics in SARS-CoV-2 spike receptor complex systems. We considered eight different point mutations present in the Omicron variant. Root-mean-square deviation and binding energy analysis showed that S477N and Omicron did not eject HD throughout the simulation. Hydrogen bond distribution analysis highlighted the involvement of hydrogen bonding in mutant-HD stabilization, especially for S477N and Omicron. Root-mean-square fluctuation analysis revealed evidence of Y505H destabilization on complex systems, while distal-end loop mutations increased loop flexibility for all models bearing the three mutations. Per-residue energy decomposition demonstrated that Q493R substitution increased HD interaction. Free energy landscape and essential dynamics through principal component analysis provided insights into the conformational subspace distribution of mutant model molecular dynamics trajectories. In conclusion, significant mutations contributed to the HD interaction in different ways. S477N has shown significant binding contributions through favorable ligand interaction, while other mutations contribute via conformational modifications, increased affinity due to sidechain mutations, and increased loop flexibility.

Strategies to enhance the hydrolytic activity of Escherichia coli BL21 penicillin G acylase based on heterologous expression and targeted mutagenesis

Penicillin G acylase (PGA) serves as a critical biocatalyst for the hydrolysis of penicillin G, yielding 6-aminopenicillanic acid, a vital precursor for β-lactam semi-synthetic antibiotics. The catalytic efficiency of PGA, however, remains suboptimal in native Escherichia coli strains. To improve this, E. coli BL21 was engineered as a microbial cell factory via heterologous expression and site-directed mutagenesis to enhance PGA activity. The heterologous pga gene from Providencia rettgeri was integrated into E. coli BL21 (DE3) for the biosynthesis of PGA, achieving a PGA activity of 253 ± 2 U/mL after 16 hours of fermentation. The N167 site underwent mutation, producing the sites N167A and N167I. Plasmids carrying these mutations were introduced into E. coli BL21(DE3), and the enzymatic activities were recorded as 293 ± 3 U/mL for the N167A mutant and 238 ± 2 U/mL for the N167I mutant. This study not only introduces a novel approach to enhancing PGA activity but also illustrates the potential for catalytic optimization through targeted modifications of the enzyme's active site.

Elucidating specific interactions for designing novel pyrrolamide derivatives as potential GyrB inhibitors based on ab initio fragment molecular orbital calculations

Tuberculosis (TB), the second leading infectious killer, causes serious public health problems worldwide. To develop novel anti-TB agents, many biochemical studies have targeted the subunit B of DNA gyrase (GyrB), which captures a second DNA segment and responses for ATP hydrolysis. Here, we investigated specific interactions between GyrB residues and existing pyrrolamide derivatives at an electronic level using ab initio fragment molecular orbital (FMO) calculations and designed potent inhibitors against GyrB. The evaluated binding affinities between GyrB and pyrrolamides were confirmed to be consistent with the IC50 values obtained from previous experiments. Thus, we employed the most potent pyrrolamide (compound 1) as a lead compound and proposed novel pyrrolamide derivatives. The specific interactions between GyrB and these derivatives were investigated using molecular mechanic optimizations and FMO calculations. The results revealed that our proposed derivatives had strong hydrogen bonds with Asp79 and Arg141 and exhibited electrostatic interactions with Glu56 and Ile84 of GyrB. In addition, the binding affinity between GyrB and compound 1 was enhanced significantly by the replacement at the R3 site of compound 1. The present results may provide structural concepts for the rational design of potent GyrB inhibitors as anti-TB agents.Communicated by Ramaswamy H. Sarma.

Different receptor models show differences in ligand binding strength and location: a computational drug screening for the tick-borne encephalitis virus

The tick-borne encephalitis virus (TBE) is a neurotrophic disease that has spread more rapidly throughout Europe and Asia in the past few years. At the same time, no cure or specific therapy is known to battle the illness apart from vaccination. To find a pharmacologically relevant drug, a computer-aided drug screening was initiated. Such a procedure probes a possible binding of a drug to the RNA Polymerase of TBE. The crystal structure of the receptor, however, includes missing and partially modeled regions, which rendered the structure incomplete and of questionable use for a thorough drug screening procedure. The quality of the receptor model was addressed by studying three putative structures created. We show that the choice of receptor models greatly influences the binding affinity of potential drug molecules and that the binding location could also be significantly impacted. We demonstrate that some drug candidates are unsuitable for one model but show decent results for another. Without any prejudice on the three employed receptor models, the study reveals the imperative need to investigate the receptor structure before drug binding is probed whether experimentally or computationally.

Conformation and binding of 12 Microcystin (MC) congeners to PPP1 using molecular dynamics simulations: A potential approach in support of an improved MC risk assessment

Microcystins (MCs) occur frequently during cyanobacterial blooms worldwide, representing a group of currently about 300 known MC congeners, which are structurally highly similar. Human exposure to MCs via contaminated water, food or dietary supplements can lead to severe intoxications with ensuing high morbidity and in some cases mortality. Currently, one MC congener (MC-LR) is almost exclusively considered for risk assessment (RA) by the WHO. Many MC congeners co-occur during bloom events, of which MC-LR is not the most toxic. Indeed, MC congeners differ dramatically in their inherent toxicity, consequently raising question about the reliability of the WHO RA and the derived guidance values. Molecular dynamics (MD) simulation can aid in understanding differences in toxicity, as experimental validation for all known MC congeners is not feasible. Therefore, we present MD simulations of a total of twelve MC congeners, of which eight MC congeners were simulated for the first time. We show that depending on their structure and toxicity class, MCs adapt to different backbone conformations. These backbone conformations are specific to certain MC congeners and can change or shift to other conformations upon binding to PPP1, affecting the stability of the binding. Analysis of the interactions with PPP1 demonstrated that there are frequently occurring patterns for individual MC congeners, and that published PPP interactions could be reproduced. In addition, common but also unique patterns were found for individual MC congeners, suggesting differences in binding behaviour. The MD simulations presented here therefore enhance our understanding of MC congener-specific differences and demonstrated that congener-specific investigations are prerequisite for allowing characterisation of yet untested or even unknown MC congeners, thereby allowing for a novel potential approach in support of an improved RA of microcystins in humans.

GlnA3Mt is able to glutamylate spermine but it is not essential for the detoxification of spermine in Mycobacterium tuberculosis

Mycobacterium tuberculosis is well adapted to survive and persist in the infected host, escaping the host's immune response. Since polyamines such as spermine, which are synthesized by infected macrophages, are able to inhibit the growth of M. tuberculosis, the pathogen needs strategies to cope with these toxic metabolites. The actinomycete Streptomyces coelicolor, a close relative of M. tuberculosis, makes use of a gamma-glutamylation pathway to functionally neutralize spermine. We therefore considered whether a similar pathway would be functional in M. tuberculosis. In the current study, we demonstrated that M. tuberculosis growth was inhibited by the polyamine spermine. Using in vitro enzymatic assays we determined that GlnA3Mt (Rv1878) possesses genuine gamma-glutamylspermine synthetase catalytic activity. We further showed that purified His-Strep-GlnA3Mt, as well as native GlnA3Mt, prefer spermine as a substrate over putrescine, cadaverine, spermidine, or other monoamines and amino acids, suggesting that GlnA3Mt may play a specific role in the detoxification of the polyamine spermine. However, the deletion of the glnA3 gene in M. tuberculosis did not result in growth inhibition or enhanced sensitivity of M. tuberculosis in the presence of high spermine concentrations. Gene expression analysis of spermine-treated M. tuberculosis revealed no difference in the level of glnA3Mt expression relative to untreated cells, whereas a gene encoding a previously characterized efflux pump (Mmr; rv3065) was significantly upregulated. This suggests that bacterial survival under elevated spermine concentrations can not only be achieved by detoxification of spermine itself but also by mechanisms resulting in decreased spermine levels in the bacteria.

IMPORTANCE

Upon Mycobacterium tuberculosis infection macrophages synthesize the polyamine spermine, which at elevated concentrations is toxic for M. tuberculosis. Based on our investigations of spermine resistance in the closely related actinomycete Streptomyces coelicolor, we hypothesized that the glutamylspermine synthetase GlnA3 may be responsible for the resistance of M. tuberculosis against toxic spermine. Here we show that GlnA3Mt can indeed covalently modify spermine via glutamylation. However, GlnA3Mt is probably not the only resistance mechanism since a glnA3 null mutant of M. tuberculosis can survive under spermine stress. Gene expression studies suggest that an efflux pump may participate in resistance. Thus a combination of GlnA3Mt and specific efflux pumps acting as putative spermine transporters may constitute an active spermine-detoxification system in M. tuberculosis.

Functional Characterization of an Aldol Condensation Synthase PheG for the Formation of Hispidin from Phellinus Igniarius

Hispidin (1) is a polyphenolic compound with a wide range of pharmacological activities that is distributed in both plants and fungi. In addition to natural extraction, hispidin can be obtained by chemical or enzymatic synthesis. In this study, the identification and characterization of an undescribed enzyme, PheG, from Phellinus igniarius (P. igniarius), which catalyzes the construction of a key C─C bond in the enzymatic synthesis of hispidin are reported. It is demonstrated in vitro that PheG generates hispidin by catalyzing C─C bond formation in the aldol condensation reaction. Based on these results, a plausible pathway for hispidin biosynthesis is proposed by utilizing the primary triacetic acid lactone (TAL, 2) and 3,4-dihydroxybenzaldehyde (3). The mechanisms for the aldol condensation reaction of PheG are investigated using molecular dynamics (MD) simulations, molecular mechanics/generalized Born surface area (MM/GBSA) binding free energy calculations, density functional theory, and site-specific mutations. The locations of the key amino acid residues that catalyze the conversion of substrates 2 and 3 to hispidin at the active site of PheG-1 are identified. This study provides a new method for preparing hispidin with high efficiency and low cost.

The mutational spectrum of NRAS gene discovers a novel frameshift mutation (E49R) in Saudi colorectal cancer patients

Colorectal cancer (CRC) is a major health problem the world face currently and one of the leading causes of death worldwide. CRC is genetically heterogeneous and multiple genetic aberrations may appear on course of the disease throughout patient's lifetime. Genetic biomarkers such as BRAF, KRAS, and NRAS may provide early precision treatment options that are crucial for patient survival and well-being. The aim of this study was to identify pathogenic mutations in the NRAS gene causing colorectal cancer in the Saudi population. We enrolled 80 CRC tumor tissue samples and performed molecular analyses to establish the mutation spectrum status in the western region of Saudi Arabia. We identified 5 different mutations in 10 patients, 4 of whom were reported previously (G10R, E37K, Q61K, and Q61*) in the literature while we discovered one novel lethal insertion mutation (E49R). A novel identified insertion mutation was present in the third codon of the NRAS gene [c.145 insA (p.Glu49ArgTer85)], causing a frameshift in the amino acid sequence of the protein, and leading to an aberrant and truncated protein of 85 amino acids. Subsequent bioinformatics analysis showed that the mutation was highly deleterious and affected protein function to a greater extent. This identification may further improve the prognosis of CRC and benefit subsequent treatment choices.

The OnSPN2 from the nipa palm hispid beetle Octodonta nipae is a multipurpose defense tool against proteases from different peptidase families

Serpins (serine protease inhibitors) constitute a superfamily of proteins with functional diversity and unusual conformational flexibility. In insects, serpins act as multiple inhibitors, by forming inactive acyl-enzyme complexes, in regulating Spätzles activation, phenoloxidases (POs) activity, and other cytokines. In this study, we present the cloning and characterization of Octodonta nipae serpin2 (OnSPN2), a 415 residues protein homologous to Tenebrio molitor 42Dd-like. Notably, OnSPN2 features an arginine residue (R364) at the P1 position, and additional arginine residues (R362, R367) at the P3 and P3' positions, respectively which is crucial for protease inhibition. Immunohistochemistry (IHC) and Western blot analyses revealed that OnSPN2 is primarily synthesized in plasmatocytes and then released into the plasma to exert its function. RNA interference results indicated that OnSPN2 knockdown may depress serine protease in melanization and remarkably increase the transcript level of Attacin in hemolymph, but its messenger RNA levels were not changed upon immune induction. Reciprocal co-immunoprecipitation assay results confirmed that OnSPN2 binds to OnPPAF1 and OnSP8, indicating its role as a negative regulator in the PO and AMP pathway. Intriguingly, several cathepsin-L isoforms were identified in the OnSPN2 immunoprecipitated samples. The cathepsin-L inhibition assays and protein-protein docking results, identified cathepsin-L as a potential target of OnSPN2. These results indicate that OnSPN2 is produced as an intracellular resident and additionally is associated with the PO and AMP pathway. OnSPN2 represents a multiple defense tool that may provide multiple antiproteolytic functions.

Accurate Prediction of Protein Complex Stoichiometry by Integrating AlphaFold3 and Template Information

Protein structure prediction methods require stoichiometry information (i.e., subunit counts) to predict the quaternary structure of protein complexes. However, this information is often unavailable, making stoichiometry prediction crucial for complexes with unknown stoichiometry. Despite its importance, few computational methods address this challenge. In this study, we present an approach that integrates AlphaFold3 structure predictions with homologous template data to predict stoichiometry. The method generates candidate stoichiometries, builds structural models for them using AlphaFold3, ranks them based on AlphaFold3 scores, and further refine predictions with template-based information when available. In the 16th community-wide Critical Assessment of Techniques for Protein Structure Prediction (CASP16), our method achieved 71.4% top-1 accuracy and 92.9% top-3 accuracy, outperforming other predictors in terms of the overall performance. This demonstrates the complementary strengths of AlphaFold3- and template-based predictions and highlights its applicability for uncharacterized protein complexes lacking stoichiometry data.

Mechanism of low molecular weight impurity formation in an IgG1 monoclonal antibody formulation

Formulation robustness study was performed for a biosimilar monoclonal antibody (IgG1) manufactured at Dr. Reddy's Laboratory, where the pH and concentration level of excipients in the drug product formulation were systematically varied from the target formulation. It was observed that the IgG1 formulation having relatively low pH and high citrate (buffer salt) concentration were predisposed to the formation of low molecular weight impurities. Mass spectrometry analysis of the mAb1 fragments detected the pyroglutamate species from LC-LC dimer and fragmentation in the -DKTH- amino acid sequence of the heavy chain. Blind docking indicated binding of citrate with Lysine 222 residue in the proximity of Cys224 could have potentially fragmented IgG1.

Exploring 1-alkene biosynthesis in bacterial antagonists and Jeotgalicoccus sp. ATCC 8456

Terminal olefins are important platform chemicals, drop-in compatible hydrocarbons and also play an important role as biocontrol agents of plant pathogens. Currently, 1-alkenes are derived from petroleum, although microbial biosynthetic routes are known. Jeotgalicoccus sp. ATCC 8456 produces 1-alkenes via the fatty acid decarboxylase OleTJE. UndA and UndB are recently identified non-heme iron oxidases converting medium-chain fatty acids into terminal alkenes. Our knowledge about the diversity and natural function of OleTJE, UndA, and UndB homologs is scarce. We applied a combined screening strategy-solid-phase microextraction coupled with gas chromatography-mass spectrometry (SPME GC-MS) and polymerase chain reaction (PCR)-based amplification-to survey an environmental strain collection for microbial 1-alkene producers and their corresponding enzymes. Our results reinforce the high level of conservation of UndA and UndB genes across the genus Pseudomonas. In vivo production of defined 1-alkenes (C9-C13; C15; C19) was directed by targeted feeding of fatty acids. Lauric acid feeding enabled 1-undecene production to a concentration of 3.05 mg l-1 in Jeotgalicoccus sp. ATCC 8456 and enhanced its production by 105% in Pseudomonas putida 1T1 (1.10 mg l-1). Besides, whole genome sequencing of Jeotgalicoccus sp. ATCC 8456 enabled reconstruction of the 1-alkene biosynthetic pathway. These results advance our understanding of microbial 1-alkene synthesis and the underlying genetic basis.

Differential Inhibition by Cenobamate of Canonical Human Nav1.5 Ion Channels and Several Point Mutants

Cenobamate is a new and highly effective antiseizure compound used for the treatment of adults with focal onset seizures and particularly for epilepsy resistant to other antiepileptic drugs. It acts on multiple targets, as it is a positive allosteric activator of γ-aminobutyric acid type A (GABAA) receptors and an inhibitor of neuronal sodium channels, particularly of the late or persistent Na+ current. We recently evidenced the inhibitory effects of cenobamate on the peak and late current component of the human cardiac isoform hNav1.5. The determined apparent IC50 values of 87.6 µM (peak) and 46.5 µM (late current) are within a clinically relevant range of concentrations (the maximal plasma therapeutic effective concentration for a daily dose of 400 mg in humans is 170 µM). In this study, we built a 3D model of the canonical hNav1.5 channel (UniProt Q14524-1) in open conformation using AlphaFold2, embedded it in a DPPC lipid bilayer, corrected the residue protonation state (pH 7.2) with H++, and added 2 Na+ ions in the selectivity filter. By molecular docking, we found the cenobamate binding site in the central cavity. We identified 10-point mutant variants in the binding site region and explored them via docking and MD. Mutants N1462K/Y (rs1064795922, rs199473614) and M1765R (rs752476527) (by docking) and N932S (rs2061582195) (by MD) featured higher predicted affinity than wild-type.

Insights into electron transfer and bifurcation of the Synechocystis sp. PCC6803 hydrogenase reductase module

The NAD+-reducing soluble [NiFe] hydrogenase (SH) is the key enzyme for production and consumption of molecular hydrogen (H2) in Synechocystis sp. PCC6803. In this study, we focused on the reductase module of the SynSH and investigated the structural and functional aspects of its subunits, particularly the so far elusive role of HoxE. We demonstrated the importance of HoxE for enzyme functionality, suggesting a regulatory role in maintaining enzyme activity and electron supply. Spectroscopic analysis confirmed that HoxE and HoxF each contain one [2Fe2S] cluster with an almost identical electronic structure. Structure predictions, alongside experimental evidence for ferredoxin interactions, revealed a remarkable similarity between SynSH and bifurcating hydrogenases, suggesting a related functional mechanism. Our study unveiled the subunit arrangement and cofactor composition essential for biological electron transfer. These findings enhance our understanding of NAD+-reducing [NiFe] hydrogenases in terms of their physiological function and structural requirements for biotechnologically relevant modifications.

Elucidating the potential of bioactive of Trichoderma sp.. in combating pathogenesis by Fusarium sp.. by targeting pectin lyases: a bioinformatics approach

Pectin lyase is an industrially important enzyme, predominately used in fruit juice clarification and retting of fibers. It also promotes pathogenesis via the degradation of the pectin. The phytopathogen, Fusarium infects various crops and causes several diseases. Trichoderma sp. is a promising biocontrol agent that is vital in maintaining plant health and disease prevention. In the current study, a computational approach utilizing structure prediction, molecular docking, molecular dynamics, and MM-PBSA analysis was used to analyze the potential role of bioactive compounds secreted by Trichoderma sp. in inhibiting the pectin lyase enzyme from Fusarium proliferatum, F. fujikuroi, F. graminearum, F. oxysporum and F. verticillioides. Molecular docking with secondary metabolites revealed that Viridiofungin A secreted by Trichoderma harzianum and Virone secreted by T. virens are bioactive compounds with immense potential to inhibit PNLs of Fusarium species. Further, the rigidity of the structure and stability of the docked complex were confirmed via Molecular dynamic simulations assessed through multiple parameters from the simulation trajectory data. Dual culture assay of T. harzianum and T. virens with F. proliferatum, F. fujikuroi, F. graminearum, F. oxysporum, and F. verticillioides showed variable mycelial inhibition. The research provides insight into the potential of the bioactive compounds secreted by Trichoderma species as an effective agent for the inhibition of pectin lyases produced by phytopathogens, especially Fusarium species. The proposed research can be used to develop bioformulations that function as biopesticides, offering a sustainable replacement for chemical products.

Exploring the Structure-Function Relationships in a 5-Aminolevulinic Acid Synthase and the Use of Protein Engineering to Expand its Substrate Range

5-Aminolevulinate synthase (ALAS) is a PLP-dependent enzyme that catalyzes the production of 5-aminolevulinate from succinyl-CoA and glycine. Its ability to catalyze the essentially irreversible C-C bond formation has significant potential in chemoenzymatic synthesis of α-amino ketones. Native ALAS, unfortunately, is extremely substrate-selective, and this seriously limits its synthetic utility. Here, we have used three different protein engineering strategies to overcome this problem for the acyl-CoA substrate. By combining previously reported mutation results and structural analysis, a series of site-saturation mutagenesis/screening efforts were focused on R21, T82, N84, and T362 of Rhodopseudomonas palustris ALAS. These yielded single, double, and triple mutants with significantly improved substrate ranges. The steady-state kinetic parameters of several key variants were determined. These data were analyzed in the framework of the ALAS catalytic mechanism to identify the steps that may have been impacted. The most active variant was used in a larger-scale reaction to demonstrate its synthetic potential. Taken together, our results show how ALAS might become a useful biocatalyst for α-amino ketone synthesis and have also allowed us to comment on the relative merits of each the three protein engineering strategies utilized.

Recent advances and challenges in protein complex model accuracy estimation

Estimation of model accuracy plays a crucial role in protein structure prediction, aiming to evaluate the quality of predicted protein structure models accurately and objectively. This process is not only key to screening candidate models that are close to the real structure, but also provides guidance for further optimization of protein structures. With the significant advancements made by AlphaFold2 in monomer structure, the problem of single-domain protein structure prediction has been widely solved. Correspondingly, the importance of assessing the quality of single-domain protein models decreased, and the research focus has shifted to estimation of model accuracy of protein complexes. In this review, our goal is to provide a comprehensive overview of the reference and statistical metrics, as well as representative methods, and the current challenges within four distinct facets (Topology Global Score, Interface Total Score, Interface Residue-Wise Score, and Tertiary Residue-Wise Score) in the field of complex EMA.

New bacteriophage-derived lysins, LysJ and LysF, with the potential to control Bacillus anthracis

Bacillus anthracis is an etiological agent of anthrax, a severe zoonotic disease that can be transmitted to people and cause high mortalities. Bacteriophages and their lytic enzymes, endolysins, have potential therapeutic value in treating infections caused by this bacterium as alternatives or complements to antibiotic therapy. They can also be used to identify and detect B. anthracis. Endolysins of two B. anthracis Wbetavirus phages, J5a and F16Ba which were described by us recently, differ significantly from the best-known B. anthracis phage endolysin PlyG from Wbetavirus genus bacteriophage Gamma and a few other Wbetavirus genus phages. They are larger than PlyG (351 vs. 233 amino acid residues), contain a signal peptide at their N-termini, and, by prediction, have a different fold of cell binding domain suggesting different structural basis of cell epitope recognition. We purified in a soluble form the modified versions of these endolysins, designated by us LysJ and LysF, respectively, and depleted of signal peptides. Both modified endolysins could lyse the B. anthracis cell wall in zymogram assays. Their activity against the living cells of B. anthracis and other species of Bacillus genus was tested by spotting on the layers of bacteria in soft agar and by assessing the reduction of optical density of bacterial suspensions. Both methods proved the effectiveness of LysJ and LysF in killing the anthrax bacilli, although the results obtained by each method differed. Additionally, the lytic efficiency of both proteins was different, which apparently correlates with differences in their amino acid sequence. KEY POINTS: • LysJ and LysF are B. anthracis-targeting lysins differing from lysins studied so far • LysJ and LysF could be overproduced in E. coli in soluble and active forms • LysJ and LysF are active in killing cells of B. anthracis virulent strains.

Trans-Placental Transfer Mechanisms of Aromatic Amine Antioxidants (AAs) and p-Phenylenediamine Quinones (PPD-Qs): Evidence from Human Gestation Exposure and the Rat Uterine Perfusion Model

Aromatic amine antioxidants (AAs), as rubber additives, and their ozone photochemical oxidation products of p-phenylenediamine quinone (PPD-Qs) have attracted great attention recently due to their wide environmental occurrences and toxicity. However, there is currently no research on the exposure risks during pregnancy and their trans-placental transfer mechanisms. Herein, 20 AAs and six PPD-Qs were analyzed in 60 maternal urine and fifty-six amniotic fluid samples (n = 53 pairs). ΣAAs (median: 8.57 and 15.4 ng/mL) and ΣPPD-Qs (0.236 and 2.29 ng/mL) were both observed, where the median concentration of PPD-Qs was significantly (p < 0.05) higher than that of the parent PPDs (0.130 and 0.092 ng/mL) in the maternal urine and amniotic fluid samples, respectively. The result of the self-established rat uterine perfusion model and molecular docking analysis suggested that passive diffusion and active transport patterns were involved in the trans-placental transfer. This study will raise concerns regarding intrauterine exposure and the trans-placental transfer mechanisms to AAs/PPD-Qs during pregnancy.

(+)-3,6-Epoxymaaliane: A Novel Derivative of (+)-Bicyclogermacrene Oxidation Catalyzed by CYP450 BM3-139-3 and Its Variants

(+)-Bicyclogermacrene is a sesquiterpene compound found in various plant essential oils and serves as a crucial precursor for multiple biologically active compounds. Many derivatives of (+)-bicyclogermacrene have been shown to exhibit valuable bioactivities. Cytochrome P450 BM3 from Bacillus megaterium can catalyze a variety of substrates and different types of oxidation reactions, making it become a powerful tool for oxidizing terpenes. In this study, we employed P450 BM3-139-3 variant for in vitro enzymatic oxidation of (+)-bicyclogermacrene, identifying a novel oxidized derivative of (+)-bicyclogermacrene, named (+)-3,6-epoxymaaliane, and an unknown sesquiterpenoid in a ratio of 70 : 30 (by GC peak area). (+)-3,6-Epoxymaaliane showed demonstrated antibacterial activities toward Escherichia coli and Staphylococcus aureus. To obtain a better variant of the monooxygenase with a high selectivity to form (+)-3,6-epoxymaaliane, we combined alanine scanning with the "Focused Rational Iterative Site-Specific Mutagenesis" (FRISM) strategy to modify the closest residues within 5 Å radius surrounding the substrate to create a small-but-smart library of mutants. Consequently, it gave an optimal variant with 1.6-fold improvement, in a turnover number (TON) of up to 964 toward (+)-3,6-epoxymaaliane production with a higher product selectivity.

Drug Discovery in the Age of Artificial Intelligence: Transformative Target-Based Approaches

The complexities inherent in drug development are multi-faceted and often hamper accuracy, speed and efficiency, thereby limiting success. This review explores how recent developments in machine learning (ML) are significantly impacting target-based drug discovery, particularly in small-molecule approaches. The Simplified Molecular Input Line Entry System (SMILES), which translates a chemical compound's three-dimensional structure into a string of symbols, is now widely used in drug design, mining, and repurposing. Utilizing ML and natural language processing techniques, SMILES has revolutionized lead identification, high-throughput screening and virtual screening. ML models enhance the accuracy of predicting binding affinity and selectivity, reducing the need for extensive experimental screening. Additionally, deep learning, with its strengths in analyzing spatial and sequential data through convolutional neural networks (CNNs) and recurrent neural networks (RNNs), shows promise for virtual screening, target identification, and de novo drug design. Fragment-based approaches also benefit from ML algorithms and techniques like generative adversarial networks (GANs), which predict fragment properties and binding affinities, aiding in hit selection and design optimization. Structure-based drug design, which relies on high-resolution protein structures, leverages ML models for accurate predictions of binding interactions. While challenges such as interpretability and data quality remain, ML's transformative impact accelerates target-based drug discovery, increasing efficiency and innovation. Its potential to deliver new and improved treatments for various diseases is significant.

Detectability of cytokine and chemokine using ELISA, following sample-inactivation using Triton X-100 or heat

Clinical samples are routinely inactivated before molecular assays to prevent pathogen transmission. Antibody-based assays are sensitive to changes in analyte conformation, but the impact of inactivation on the analyte detectability has been overlooked. This study assessed the effects of commonly used inactivation-methods, Triton X-100 (0.5%) and heat (60 °C, 1 h), on cytokine/chemokine detection in plasma, lung aspirates, and nasopharyngeal samples. Heat significantly reduced analyte detectability in plasma (IL-12p40, IL-15, IL-16, VEGF, IL-7, TNF-β) by 33-99% (p ≤ 0.02), while Triton X-100 minimally affected analytes in plasma and nasopharyngeal samples (11-37%, p ≤ 0.04) and had no significant impact on lung aspirates. Structural analysis revealed that cytokines affected by heat had more hydrophobic residues and higher instability-indices. As the protein-detectability was affected differently in different sample types, the sample environment could also influence protein stability. This underscores the importance of selecting the most suitable inactivation methods for clinical samples to ensure accurate cytokine/chemokine analysis in both clinical and research settings.

A new silicon phthalocyanine dye induces pyroptosis in prostate cancer cells during photoimmunotherapy

Photoimmunotherapy (PIT) combines the specificity of antibodies with the cytotoxicity of light activatable photosensitizers (PS) and is a promising new cancer therapy. We designed and synthesized, in a highly convergent manner, the silicon phthalocyanine dye WB692-CB2, which is novel for being the first light-activatable PS that can be directly conjugated via a maleimide linker to cysteines. In the present study we conjugated WB692-CB2 to a humanized antibody with engineered cysteines in the heavy chains that specifically targets the prostate-specific membrane antigen (PSMA). The resulting antibody dye conjugate revealed high affinity and specificity towards PSMA-expressing prostate cancer cells and induced cell death after irradiation with red light. Treated cells exhibited morphological characteristics associated with pyroptosis. Mechanistic studies revealed the generation of reactive oxygen species, triggering a cascade of intracellular events involving lipid peroxidation, caspase-1 activation, gasdermin D cleavage and membrane rupture followed by release of pro-inflammatory cellular contents. In first in vivo experiments, PIT with our antibody dye conjugate led to a significant reduction of tumor growth and enhanced overall survival in mice bearing subcutaneous prostate tumor xenografts. Our study highlights the future potential of the new phthalocyanine dye WB692-CB2 as PS for the fluorescence-based detection and PIT of cancer, including local prostate tumor lesions, and systemic activation of anti-tumor immune responses by the induction of pyroptosis.

Two Novel Variants in the CHRNA2 and SCN2A Genes in Italian Patients with Febrile Seizures

BACKGROUND

Febrile seizures (FSs) are the most common form of epilepsy in children aged between six months and five years. The exact cause is unknown, but several studies have demonstrated the importance of genetic predisposition, with increasing involvement of receptors and ion channels. The present study aims to identify novel pathogenic variants in Italian patients with FSs.

METHODS

We performed targeted panel sequencing in a cohort of 21 patients with FSs. In silico analysis was performed to predict the pathogenic role of the resulting variants.

RESULTS

We found two novel variants segregating in two families with FSs: c.1021C>G (p.Leu341Val) in the CHRNA2 gene and c.140A>G (p.Glu47Gly) in SCN2A.

CONCLUSIONS

The c.1021C>G (p.Leu341Val) variant leads to a codon change of highly conserved leucine to valine at position 341 and is located in segments M3 of the subunit, which is important for channel gating. The c.140A>G (p.Glu47Gly) variant causes a substitution of glutamic acid with glycine at position 47 of the protein, which is highly conserved across the species. Moreover, it is located in the N-terminal domain, a region commonly affected in ASD, which impacts the inactivation kinetics and voltage dependence of steady-state activation. Further analyses are needed to better explain the role of CHRNA2 and SCN2A in the development of febrile seizures.

Novel bi-allelic DNAH3 variants cause oligoasthenoteratozoospermia

Background

Oligoasthenoteratozoospermia (OAT) is a widespread cause of male infertility. One of the usual clinical manifestations of OAT is multiple morphological abnormalities of the sperm flagella (MMAF), which are frequently associated with mutations and defects in the dynein family. However, the relationship between the newly identified Dynein Axonemal Heavy Chain 3 (DNAH3) mutation and oligonasthenospermia in humans has not yet been established.

Methods

Whole exome sequencing, pathogenicity analysis, and species conservation analysis of mutation sites were conducted on two patients from different unrelated families with DNAH3 mutations. We identified representative mutation sites and predicted the protein structure following these mutations. The sperm characteristics of the two patients with DNAH3 mutations were verified using Papanicolaou staining, scanning electron microscopy, and transmission electron microscopy. Additionally, mRNA and protein levels were assessed through RT-qPCR and Western blotting.

Results

The biallelic mutations in the first progenitor included a heterozygous deletion and insertion, c.6535_6536 delinsAC (to infect mutation (p.Asp2179Thr), and stop codon premutation, c.3249G > A (p.Trp1083Ter). In Family II, the patient (P2) harbored a DNAH3 heterozygous missense mutation, c. 10439G> A(p.Arg3480Gln), along with a stop codon premutation, (c.10260G > A; p.Trp3420Ter). Patients with premature termination of transcription or translation due to DNAH3 mutations exhibit OAT phenotypes, including fibrous sheath dysplasia and multiple tail malformations. We identified the representative sites after mutation, predicted the protein structure, and assessed changes in the protein levels of DNAH3 and related genes following mutations. Notably,a significant reduction in DNAH3 protein expression was validated in these patients. We may explore in the future how DNAH3 affects sperm motility and quality through regulatory mechanisms involving protein structural changes.

Conclusion

Novel biallelic mutations in DNAH3, especially those resulting in a premature stop codon, may alter protein expression, structure, and active site, leading to spermatogenic failure and potentially inducing OAT. The discovery of new mutations in DNAH3 may be the key to the diagnosis and treatment of OAT.

Repetitive transcranial magnetic stimulation alleviates motor impairment in Parkinson's disease: association with peripheral inflammatory regulatory T-cells and SYT6

BACKGROUND

Repetitive transcranial magnetic stimulation (rTMS) has been used to treat various neurological disorders. However, the molecular mechanism underlying the therapeutic effect of rTMS on Parkinson's disease (PD) has not been fully elucidated. Neuroinflammation like regulatory T-cells (Tregs) appears to be a key modulator of disease progression in PD. If rTMS affects the peripheral Tregs in PD remains unknown.

METHODS

Here, we conducted a prospective clinical study (Chinese ClinicalTrials. gov: ChiCTR 2100051140) involving 54 PD patients who received 10-day rTMS (10 Hz) stimulation on the primary motor cortex (M1) region or sham treatment. Clinical and function assessment as well as flow cytology study were undertaken in 54 PD patients who were consecutively recruited from the department of neurology at Zhujiang hospital between September 2021 and January 2022. Subsequently, we implemented flow cytometry analysis to examine the Tregs population in spleen of MPTP-induced PD mice that received rTMS or sham treatment, along with quantitative proteomic approach reveal novel molecular targets for Parkinson's disease, and finally, the RNA interference method verifies the role of these new molecular targets in the treatment of PD.

RESULTS

We demonstrated that a 10-day rTMS treatment on the M1 motor cortex significantly improved motor dysfunction in PD patients. The beneficial effects persisted for up to 40 days, and were associated with an increase in peripheral Tregs. There was a positive correlation between Tregs and motor improvements in PD cases. Similarly, a 10-day rTMS treatment on the brains of MPTP-induced PD mice significantly ameliorated motor symptoms. rTMS reversed the downregulation of circulating Tregs and tyrosine hydroxylase neurons in these mice. It also increased anti-inflammatory mediators, deactivated microglia, and decreased inflammatory cytokines. These effects were blocked by administration of a Treg inhibitor anti-CD25 antibody in MPTP-induced PD mice. Quantitative proteomic analysis identified TLR4, TH, Slc6a3 and especially Syt6 as the hub node proteins related to Tregs and rTMS therapy. Lastly, we validated the role of Treg and rTMS-related protein syt6 in MPTP mice using the virus interference method.

CONCLUSIONS

Our clinical and experimental studies suggest that rTMS improves motor function by modulating the function of Tregs and suppressing toxic neuroinflammation. Hub node proteins (especially Syt6) may be potential therapeutic targets.

TRIAL REGISTRATION

Chinese ClinicalTrials, ChiCTR2100051140. Registered 15 December 2021, https://www.chictr.org.cn/bin/project/edit?pid=133691.

Genome-wide analysis of the SWEET gene family and its response to powdery mildew and leaf spot infection in the common oat (Avena sativa L.)

The nutritional quality and yield of oats (Avena sativa) are often compromised by plant diseases such as red leaf, powdery mildew, and leaf spot. Sugars Will Eventually be Exported Transporters (SWEETs) are newly identified sugar transporters involved in regulating plant growth and stress responses. However, the roles of SWEET genes in biotic stress responses remain uncharacterized in oats. In this study, 13 AsSWEET genes were identified across nine chromosomes of the oat genome, all of which were predicted to contain seven transmembrane regions. Phylogenetic analysis revealed four clades of AsSWEET proteins, with high homology to SWEET proteins in the Poaceae family. Collinearity analysis demonstrated strong relationships between oat and Zea mays SWEETs. Using subcellular localization prediction tools, AsSWEET proteins were predicted to localize to the plasma membrane. Promoter analysis revealed cis-acting elements associated with light response, growth, and stress regulation. Six AsSWEET proteins were predicted to interact in a network centered on AsSWEET1a and AsSWEET11. Gene expression analysis of two oat varieties, 'ForagePlus' and 'Molasses', indicated significant expression differences in several AsSWEET genes following infection with powdery mildew or leaf spot, including AsSWEET1a, AsSWEET1b, AsSWEET2b, AsSWEET3a, AsSWEET11, and AsSWEET16. These SWEET genes are potential candidates for disease resistance in oats. This study provides a foundation for understanding the regulatory mechanisms of AsSWEET genes, particularly in response to powdery mildew and leaf spot, and offers insights for enhancing oat molecular breeding.

In silico analysis and structural vaccinology prediction of Toxoplasma gondii ROP41 gene via immunoinformatics methods as a vaccine candidate

INTRODUCTION

Toxoplasma gondii (T. gondii) infects all warm-blooded animals, including humans. Currently, no effective treatments exist to prevent the generation of chronic tissue cysts in infected hosts. Therefore, developing a vaccine to protect to deal with toxoplasmosis is a promising strategy, as a single immunization could provide lifelong protective immunity. Rhoptry proteins (ROPs) play a vital role for the parasite's survival within host cells and perform critical functions during different phases of parasite invasion. Little is known about ROP41 gene. Nevertheless, Understanding the characteristics of ROP41 will enhance diagnostic and vaccine research.

MATERIALS AND METHODS

The current article provides a comprehensive analysis of the essential components of the ROP41 protein, including its transmembrane domain, physico-chemical properties, subcellular location, tertiary and secondary structures, and potential T- and B-cell epitopes. These features were determined by many bioinformatics approaches to identify possible epitopes for developing a highly effective vaccine.

RESULTS

ROP41 protein showed 36 possible post-translational modification regions. The ROP41 protein secondary structure contains 17.35 % extended strand, 33.47 % alpha-helix, and 49.18 % random coil. Also, ROP41 showed many possible B- and T-cell epitopes. According to the Ramachandran plot, 90.78 % of amino acid residues had been placed in favored, 3.28 % in outlier, and 5.94 % in allowed areas. Also, the allergenicity and antigenicity evaluation indicated that ROP41 is non-allergenic and immunogenic.

CONCLUSION

The current study offered critical basic and conceptual information on ROP41 to increase a successful vaccine in opposition to continual and acute toxoplasmosis for in addition in vivo assessments. Further research is necessary for the development of vaccines utilizing ROP41 alone or combined with various antigens.

DockQ v2: improved automatic quality measure for protein multimers, nucleic acids, and small molecules

MOTIVATION

It is important to assess the quality of modeled biomolecules to benchmark and assess the performance of different prediction methods. DockQ has emerged as the standard tool for assessing the quality of protein interfaces in model structures against given references. However, as predictions of large multimers with multiple chains become more common, DockQ needs to be updated with more functionality for robustness and speed. Moreover, as the field progresses and more methods are released to predict interactions between proteins and other types of molecules, such as nucleic acids and small molecules, it becomes necessary to have a tool that can assess all types of interactions.

RESULTS

Here, we present a complete reimplementation of DockQ in pure Python. The updated version of DockQ is more portable, faster and introduces novel functionalities, such as automatic DockQ calculations for multiple interfaces and automatic chain mapping with multi-threading. These enhancements are designed to facilitate comparative analyses of protein complexes, particularly large multi-chain complexes. Furthermore, DockQ is now also able to score interfaces between proteins, nucleic acids, and small molecules.

AVAILABILITY AND IMPLEMENTATION

DockQ v2 is available online at: https://wallnerlab.org/DockQ.

Barley as a production platform for oral vaccines in sustainable fish aquaculture

Vaccination is the most effective measure to prevent disease outbreaks in fish aquaculture, with oral vaccine administration emerging as the most practical approach. However, oral vaccines face a notable limitation due to insufficient stimulation of the complex gut-associated lymphoid tissue caused by factors such as vaccine degradation, poor absorption, and recognition by the immune cells. An innovative solution to these limitations lies in the plant-based production of recombinant vaccines. Plant cells enable the production and targeted storage of recombinant vaccines in specific cell organelles which ensure superior protection from degradation and contain natural compounds acting as adjuvants. Our study explores the potential of barley (Hordeum vulgare), a globally significant cereal crop, for producing orally administered subunit vaccines against viral infections affecting economically important fish species in the Salmonidae and Cyprinidae families. Through Agrobacterium-mediated transformation of immature barley embryos, we have generated homozygous T2 generation of transgenic barley expressing recombinant antigens of spring viremia of carp virus and infectious salmon anaemia virus. The expression of these plant-based recombinant vaccines was confirmed by immunodetection, which was supported by fluorescence observation, specifically in the seed endosperm. The antigenicity of transgenic plant material containing recombinant antigens was evaluated using an intubation model of common carp (Cyprinus carpio), revealing a substantial upregulation of the immunoglobulin transcripts in both systemic and mucosal tissues over a period of 28 days following a single dose of transgenic antigens. Collectively, these results underscore the potential of barley-based recombinant vaccines for disease prevention in fish aquaculture.

Identifying Rab2 Protein as a Key Interactor of Centrin1 Essential for Leishmania donovani Growth

Previously, we have demonstrated that deletion of a growth-regulating gene (LdCen1) in the Leishmania donovani parasite (LdCen1-/-) attenuated the parasite's intracellular amastigote growth but not the growth of extracellular promastigotes. LdCen1-/- parasites were found to be safe and efficacious against homologous and heterologous Leishmania species as a vaccine candidate in animal models. The reason for the differential growth of LdCen1-/- between the two stages of the parasite needed investigation. Here, we report that LdCen1 interacts with a novel Ras-associated binding protein in L. donovani (LdRab2) to compensate for the growth of LdCen1-/- promastigotes. LdRab2 was isolated by protein pull-down from the parasite lysate, followed by nano-LC-MS/MS identification. The RAB domain sequence and the functional binding partners of the LdRab2 protein were predicted via Search Tool for the Retrieval of Interacting Proteins (STRING) analysis. The closeness of the LdRab2 protein to other reported centrin-binding proteins with different functions in other organisms was analyzed via phylogenetic analysis. Furthermore, in vitro and in silico analyses revealed that LdRab2 also interacts with other L. donovani centrins 3-5. Since centrin is a calcium-binding protein, we further investigated calcium-based interactions and found that the binding of LdRab2 to LdCen1 and LdCen4 is calcium-independent, whereas the interactions with LdCen3 and LdCen5 are calcium-dependent. The colocalization of LdCen1 and LdRab2 at the cellular basal-body region by immunofluorescence supports their possible functional association. The elevated expression of the LdRab2 protein in the mutant promastigotes suggested a probable role in compensating for the promastigote growth of this mutant strain, probably in association with other parasite centrins.

The fork protection complex promotes parental histone recycling and epigenetic memory

The inheritance of parental histones across the replication fork is thought to mediate epigenetic memory. Here, we reveal that fission yeast Mrc1 (CLASPIN in humans) binds H3-H4 tetramers and operates as a central coordinator of symmetric parental histone inheritance. Mrc1 mutants in a key connector domain disrupted segregation of parental histones to the lagging strand comparable to Mcm2 histone-binding mutants. Both mutants showed clonal and asymmetric loss of H3K9me-mediated gene silencing. AlphaFold predicted co-chaperoning of H3-H4 tetramers by Mrc1 and Mcm2, with the Mrc1 connector domain bridging histone and Mcm2 binding. Biochemical and functional analysis validated this model and revealed a duality in Mrc1 function: disabling histone binding in the connector domain disrupted lagging-strand recycling while another histone-binding mutation impaired leading strand recycling. We propose that Mrc1 toggles histones between the lagging and leading strand recycling pathways, in part by intra-replisome co-chaperoning, to ensure epigenetic transmission to both daughter cells.

Melastatin subfamily Transient Receptor Potential channels support spermatogenesis in planarian flatworms

The Transient Receptor Potential superfamily of proteins (TRPs) form cation channels that are abundant in animal sensory systems. Amongst TRPs, the Melastatin-related subfamily (TRPMs) is composed of members that respond to temperature, pH, sex hormones, and various other stimuli. Some TRPMs exhibit enriched expression in gonads of vertebrate and invertebrate species, but their contributions to germline development remain to be determined. We identified twenty-one potential TRPMs in the planarian flatworm Schmidtea mediterranea and analyzed their anatomical distribution of expression by whole-mount in situ hybridization. Enriched expression of two TRPMs (Smed-TRPM-c and Smed-TRPM-l) was detected in testis, whereas eight TRPM genes had detectable expression in patterns representative of neuronal and/or sensory cell types. Functional analysis of TRPM homologs by RNA-interference (RNAi) revealed that disruption of Smed-TRPM-c expression results in reduced sperm development, indicating a role for this receptor in supporting spermatogenesis. Smed-TRPM-l RNAi did not result in a detectable phenotype, but it increased sperm development deficiencies when combined with Smed-TRPM-c RNAi. Fluorescence in situ hybridization revealed expression of Smed-TRPM-c in early spermatogenic cells within testes, suggesting cell-autonomous regulatory functions in germ cells for this gene. In addition, Smed-TRPM-c RNAi resulted in reduced numbers of presumptive germline stem cell clusters in asexual planarians, suggesting that Smed-TRPM-c supports establishment, maintenance, and/or expansion of spermatogonial germline stem cells. While further research is needed to identify the factors that trigger Smed-TRPM-c activity, these findings reveal one of few known examples for TRPM function in direct regulation of sperm development.

Comparative Genomic and Functional Analysis of c-di-GMP Metabolism and Regulatory Proteins in Bacillus velezensis LQ-3

Bacillus velezensis is a promising candidate for biocontrol applications. A common second messenger molecule, bis-(3,5)-cyclic-dimeric-guanosine monophosphate (c-di-GMP), has the ability to regulate a range of physiological functions that impact the effectiveness of biocontrol. However, the status of the c-di-GMP signaling pathway in biocontrol strain LQ-3 remains unknown. Strain LQ-3, which was isolated from wheat rhizosphere soil, has shown effective control of wheat sharp eyespot and has been identified as B. velezensis through whole-genome sequencing analyses. In this study, we investigated the intracellular c-di-GMP signaling pathway of LQ-3 and further performed a comparative genomic analysis of LQ-3 and 29 other B. velezensis strains. The results revealed the presence of four proteins containing the GGDEF domain, which is the conserved domain for c-di-GMP synthesis enzymes. Additionally, two proteins were identified with the EAL domain, which represents the conserved domain for c-di-GMP degradation enzymes. Furthermore, one protein was found to possess a PilZ domain, indicative of the conserved domain for c-di-GMP receptors in LQ-3. These proteins are called DgcK, DgcP, YybT, YdaK, PdeH, YkuI, and DgrA, respectively; they are distributed in a similar manner across the strains and belong to the signal transduction family. We selected five genes from the aforementioned seven genes for further study, excluding YybT and DgrA. They all play a role in regulating the motility, biofilm formation, and colonization of LQ-3. This study reveals the c-di-GMP signaling pathway associated with biocontrol features in B. velezensis LQ-3, providing guidance for the prevention and control of wheat sharp eyespot by LQ-3.

In silico analysis of balsaminol as anti-viral agents targeting SARS-CoV-2 main protease, spike receptor binding domain and papain-like protease receptors

Plant-derived phytochemicals from medicinal plants are becoming increasingly attractive natural sources of antimicrobial and antiviral agents due to their therapeutic value, mechanism of action, level of toxicity and bioavailability. The continued emergence of more immune-evasive strains and the rate of resistance to current antiviral drugs have created a need to identify new antiviral agents against SARS-CoV-2. This study investigated the antiviral potential of balsaminol, a bioactive compound from Momordica balsamina, and its inhibitory activities against SARS-CoV-2 receptor proteins. In this study, three Food and Drug Administration (FDA) COVID-19 approved drugs namely; nirmatrelvir, ritonavir and remdesivir were used as positive control. Molecular docking was performed to determine the predominant binding mode (most negative Gibbs free energy of binding/ΔG) and inhibitory activity of balsaminol against SARS-CoV-2 receptor proteins. The pharmacokinetics, toxicity, physicochemical and drug-like properties of balsaminol were evaluated to determine its potential as an active oral drug candidate as well as its non-toxicity in humans. The results show that balsaminol E has the highest binding affinity to the SARS CoV-2 papain-like protease (7CMD) with a free binding energy of - 8.7 kcal/mol, followed by balsaminol A interacting with the spike receptor binding domain (6VW1) with - 8.5 kcal/mol and balsaminol C had a binding energy of - 8.1 kcal/mol with the main protease (6LU7) comparable to the standard drugs namely ritonavir, nirmatrelvir and remdesivir. However, the ADMET and drug-like profile of balsaminol F favours it as a better potential drug candidate and inhibitor of the docked SARS-CoV-2 receptor proteins. Further preclinical studies are therefore recommended.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-024-00241-0.

Development of New Multi-Glycosylation Routes to Facilitate the Biosynthesis of Sweetener Mogrosides from Bitter Immature Siraitia Grosvenorii Using Engineered Escherichia coli

Mogrosides, which have various pharmacological activities, are mainly extracted from Siraitia grosvenorii (Luo Han Guo) and are widely used as natural zero-calorie sweeteners. Unfortunately, the difficult cultivation and long maturation time of Luo Han Guo have contributed to a shortage of mogrosides. To overcome this obstacle, we developed a highly efficient biosynthetic method using engineered Escherichia coli to synthesize sweet mogrosides from bitter mogrosides. Three UDP-glycosyltransferase (UGT) genes with primary/branched glycosylation catalytic activity at the C3/C24 sites of mogrosides were screened and tested. Mutant M3, which could catalyze the glycosylation of nine types of mogrosides, was obtained through enhanced catalytic activity. This improvement in β-(1,6)-glycosidic bond formation was achieved through single nucleotide polymorphisms and direct evolution, guided by 3D structural analysis. A new multienzyme system combining three UGTs and UDP-glucose (UDPG) regeneration was developed to avoid the use of expensive UDPG. Finally, the content of sweet mogrosides in the immature Luo Han Guo extract increased significantly from 57% to 95%. This study not only established a new multienzyme system for the highly efficient production of sweet mogrosides from immature Luo Han Guo but also provided a guideline for the high-value utilization of rich bitter mogrosides from agricultural waste and residues.

The 26 S proteasome in Entamoeba histolytica: divergence of the substrate binding pockets from host proteasomes

OBJECTIVE

Proteasomes are conserved proteases crucial for proteostasis in eukaryotes and are promising drug targets for protozoan parasites. Yet, the proteasomes of Entamoeba histolytica remain understudied. The study's objective was to analyse the differences in the substrate binding pockets of amoeba proteasomes from those of host, and computational modelling of β5 catalytic subunit, with the goal of finding selective inhibitors.

RESULTS

Comparative sequence analysis revealed differences in substrate binding sites of E. histolytica proteasomes, especially in the S1 and S3 pockets of the catalytic beta subunits, implying differences in substrate preference and susceptibility to inhibitors from host proteasomes. This was strongly supported by significantly lower sensitivity to MG132 mediated inhibition of amoebic proteasome β5 subunit's chymotryptic activity compared to human proteasomes, also reflected in lower sensitivity of E. histolytica to MG132 for inhibition of proliferation. Computational models of β4 and β5 subunits, and a docked β4-β5 model revealed a binding pocket between β4-β5, similar to that of Leishmania tarentolae. Selective inhibitors for visceral leishmaniasis, LXE408 and compound 8, docked well to this pocket. This functional and sequence-based analysis predicts differences between amoebic and host proteasomes that can be utilized to develop rationally designed, selective inhibitors against E. histolytica.

Genotype-phenotype association and functional analysis of hnRNPA1 mutations in amyotrophic lateral sclerosis

BACKGROUND

Pathogenic variants in hnRNPA1 have been reported in amyotrophic lateral sclerosis (ALS) patients. However, studies on hnRNPA1 mutant spectrum and pathogenicity of variants were rare.

METHODS

We performed whole exome sequencing of ALS-associated genes and subsequent verification of rare variants in hnRNPA1 in our ALS patients. The hnRNPA1 mutations reported in literature were reviewed and combined with our results to determine the genotype-phenotype relationship. Functional analysis of the novel variant p.G195A was performed in vitro by transfection of mutant hnRNPA1 into 293T cell.

RESULTS

Among 207 ALS patients recruited, 3 rare hnRNPA1 variants were identified (mutant frequency 1.45%), including two recurrent mutations (p.P340S and p.G283R), and a novel rare variant p.G195A. In combination with previous reports, there are 27 ALS patients with 15 hnRNPA1 mutations identified. Disease onset age was 47.90 ± 1.52 years with predominant limb onset. The p.P340S mutation caused flail arm syndrome (FAS) in two independent families with extended life expectancy. The newly identified p.G195A mutation, lying at the start of the PrLD ("prion-like" domain)/LCD (low-complexity domain), causes local structural changes in 3D protein prediction. Upon sodium arsenite exposure, mutant hnRNPA1 retained in the nucleus but deficit of cytoplasmic G3BP1-positive stress granule clearance was observed. This is different from the p.P340S mutation which caused both cytoplasmic translocation and stress granule formation. No cytoplasmic TDP-43 translocation was observed.

CONCLUSION

Mutations in hnRNPA1 are overall minor in ALS patients. The p.P340S mutation is associated with manifestation of FAS. Mutations in LCD of hnRNPA1 cause stress granule misprocessing.

AI-Driven Deep Learning Techniques in Protein Structure Prediction

Protein structure prediction is important for understanding their function and behavior. This review study presents a comprehensive review of the computational models used in predicting protein structure. It covers the progression from established protein modeling to state-of-the-art artificial intelligence (AI) frameworks. The paper will start with a brief introduction to protein structures, protein modeling, and AI. The section on established protein modeling will discuss homology modeling, ab initio modeling, and threading. The next section is deep learning-based models. It introduces some state-of-the-art AI models, such as AlphaFold (AlphaFold, AlphaFold2, AlphaFold3), RoseTTAFold, ProteinBERT, etc. This section also discusses how AI techniques have been integrated into established frameworks like Swiss-Model, Rosetta, and I-TASSER. The model performance is compared using the rankings of CASP14 (Critical Assessment of Structure Prediction) and CASP15. CASP16 is ongoing, and its results are not included in this review. Continuous Automated Model EvaluatiOn (CAMEO) complements the biennial CASP experiment. Template modeling score (TM-score), global distance test total score (GDT_TS), and Local Distance Difference Test (lDDT) score are discussed too. This paper then acknowledges the ongoing difficulties in predicting protein structure and emphasizes the necessity of additional searches like dynamic protein behavior, conformational changes, and protein-protein interactions. In the application section, this paper introduces some applications in various fields like drug design, industry, education, and novel protein development. In summary, this paper provides a comprehensive overview of the latest advancements in established protein modeling and deep learning-based models for protein structure predictions. It emphasizes the significant advancements achieved by AI and identifies potential areas for further investigation.

Genome-Wide Identification and Expression Analysis of SlNRAMP Genes in Tomato under Nutrient Deficiency and Cadmium Stress during Arbuscular Mycorrhizal Symbiosis

Arbuscular mycorrhizal (AM) fungi are well known for enhancing phosphorus uptake in plants; however, their regulating roles in cation transporting gene family, such as natural resistance-associated macrophage protein (NRAMP), are still limited. Here, we performed bioinformatics analysis and quantitative expression assays of tomato SlNRAMP 1 to 5 genes under nutrient deficiency and cadmium (Cd) stress in response to AM symbiosis. These five SlNRAMP members are mainly located in the plasma or vacuolar membrane and can be divided into two subfamilies. Cis-element analysis revealed several motifs involved in phytohormonal and abiotic regulation in their promoters. SlNRAMP2 was downregulated by iron deficiency, while SlNRAMP1, SlNRAMP3, SlNRAMP4, and SlNRAMP5 responded positively to copper-, zinc-, and manganese-deficient conditions. AM colonization reduced Cd accumulation and expression of SlNRAMP3 but enhanced SlNRAMP1, SlNRAMP2, and SlNRMAP4 in plants under Cd stress. These findings provide valuable genetic information for improving tomato resilience to nutrient deficiency and heavy metal stress by developing AM symbiosis.

Structure based functional identification of an uncharacterized protein from Coxiella burnetii involved in adipogenesis

Coxiella burnetii, the causative agent of Q fever, is an intracellular pathogen posing a significant global public health threat. There is a pressing need for dependable and effective treatments, alongside an urgency for further research into the molecular characterization of its genome. Within the genomic landscape of Coxiella burnetii, numerous hypothetical proteins remain unidentified, underscoring the necessity for in-depth study. In this study, we conducted comprehensive in silico analyses to identify and prioritize potential hypothetical protein of Coxiella burnetii, aiming to elucidate the structure and function of uncharacterized protein. Furthermore, we delved into the physicochemical properties, localization, and molecular dynamics and simulations, and assessed the primary, secondary, and tertiary structures employing a variety of bioinformatics tools. The in-silico analysis revealed that the uncharacterized protein contains a conserved Mth938-like domain, suggesting a role in preadipocyte differentiation and adipogenesis. Subcellular localization predictions indicated its presence in the cytoplasm, implicating a significant role in cellular processes. Virtual screening identified ligands with high binding affinities, suggesting the protein's potential as a drug target against Q fever. Molecular dynamics simulations confirmed the stability of these complexes, indicating their therapeutic relevance. The findings provide a structural and functional overview of an uncharacterized protein from C. burnetii, implicating it in adipogenesis. This study underscores the power of in-silico approaches in uncovering the biological roles of uncharacterized proteins and facilitating the discovery of new therapeutic strategies. The findings provide valuable preliminary data for further investigation into the protein's role in adipogenesis.

The structure assessment web server: for proteins, complexes and more

The 'structure assessment' web server is a one-stop shop for interactive evaluation and benchmarking of structural models of macromolecular complexes including proteins and nucleic acids. A user-friendly web dashboard links sequence with structure information and results from a variety of state-of-the-art tools, which facilitates the visual exploration and evaluation of structure models. The dashboard integrates stereochemistry information, secondary structure information, global and local model quality assessment of the tertiary structure of comparative protein models, as well as prediction of membrane location. In addition, a benchmarking mode is available where a model can be compared to a reference structure, providing easy access to scores that have been used in recent CASP experiments and CAMEO. The structure assessment web server is available at https://swissmodel.expasy.org/assess.

Temperature effect in the inhibition of PLA2 activity of Bothrops brazili venom by Rosmarinic and Chlorogenic acids, experimental and computational approaches

Myotoxicity caused by snakebite envenoming emerges as one of the main problems of ophidic accidents as it is not well neutralized by the current serum therapy. A promising alternative is to search for efficient small molecule inhibitors that can act against multiple venom components. Phospholipase A2 (PLA2) is frequently found in snake venom and is usually associated with myotoxicity. Thus it represents an excellent target for the search of new treatments. This work reports the effect of temperature in the inhibition of catalytic properties of PLA2 from Bothrops brazili venom by Rosmarinic (RSM) and Chlorogenic (CHL) acids through experimental and computational approaches. Three temperatures were evaluated (25, 37 and 50 °C). In the experimental section, enzymatic assays showed that RSM is a better inhibitor in all three temperatures. At 50 °C, the inhibition efficiency decayed significantly for both acids. Docking studies revealed that both ligands bind to the hydrophobic channel of the protein dimer where the phospholipid binds in the catalytic process, interacting with several functional residues. In this context, RSM presents better interaction energies due to stronger interactions with chain B of the dimer. Molecular dynamics simulations showed that RSM can establish selective interactions with ARG112B of PLA2, which is located next to residues of the putative Membrane Disruption Site in PLA2-like structures. The affinity of RSM and CHL acids towards PLA2 is mainly driven by electrostatic interactions, especially salt bridge interactions established with residues ARG33B (for CHL) and ARG112B (RSM) and hydrogen bonds with residue ASP89A. The inability of CHL to establish a stable interaction with ARG112B was identified as the reason for its lower inhibition efficiency compared to RSM at the three temperatures. Furthermore, extensive structural analysis was performed to explain the lower inhibition efficiency at 50 °C for both ligands. The analysis performed in this work provides important information for the future design of new inhibitors.Communicated by Ramaswamy H. Sarma.

Genome-Wide Analysis of the SRPP/REF Gene Family in Taraxacum kok-saghyz Provides Insights into Its Expression Patterns in Response to Ethylene and Methyl Jasmonate Treatments

Taraxacum kok-saghyz (TKS) is a model plant and a potential rubber-producing crop for the study of natural rubber (NR) biosynthesis. The precise analysis of the NR biosynthesis mechanism is an important theoretical basis for improving rubber yield. The small rubber particle protein (SRPP) and rubber elongation factor (REF) are located in the membrane of rubber particles and play crucial roles in rubber biosynthesis. However, the specific functions of the SRPP/REF gene family in the rubber biosynthesis mechanism have not been fully resolved. In this study, we performed a genome-wide identification of the 10 TkSRPP and 2 TkREF genes' family members of Russian dandelion and a comprehensive investigation on the evolution of the ethylene/methyl jasmonate-induced expression of the SRPP/REF gene family in TKS. Based on phylogenetic analysis, 12 TkSRPP/REFs proteins were divided into five subclades. Our study revealed one functional domain and 10 motifs in these proteins. The SRPP/REF protein sequences all contain typical REF structural domains and belong to the same superfamily. Members of this family are most closely related to the orthologous species T. mongolicum and share the same distribution pattern of SRPP/REF genes in T. mongolicum and L. sativa, both of which belong to the family Asteraceae. Collinearity analysis showed that segmental duplication events played a key role in the expansion of the TkSRPP/REFs gene family. The expression levels of most TkSRPP/REF members were significantly increased in different tissues of T. kok-saghyz after induction with ethylene and methyl jasmonate. These results will provide a theoretical basis for the selection of candidate genes for the molecular breeding of T. kok-saghyz and the precise resolution of the mechanism of natural rubber production.

EGG: Accuracy Estimation of Individual Multimeric Protein Models Using Deep Energy-Based Models and Graph Neural Networks

Reliable and accurate methods of estimating the accuracy of predicted protein models are vital to understanding their respective utility. Discerning how the quaternary structure conforms can significantly improve our collective understanding of cell biology, systems biology, disease formation, and disease treatment. Accurately determining the quality of multimeric protein models is still computationally challenging, as the space of possible conformations is significantly larger when proteins form in complex with one another. Here, we present EGG (energy and graph-based architectures) to assess the accuracy of predicted multimeric protein models. We implemented message-passing and transformer layers to infer the overall fold and interface accuracy scores of predicted multimeric protein models. When evaluated with CASP15 targets, our methods achieved promising results against single model predictors: fourth and third place for determining the highest-quality model when estimating overall fold accuracy and overall interface accuracy, respectively, and first place for determining the top three highest quality models when estimating both overall fold accuracy and overall interface accuracy.

Characterization of Sf9 cell clones with differential susceptibilities to Sf-rhabdovirus X+3.7 and Sf-rhabdovirus X- replication

Our laboratory previously discovered a novel rhabdovirus in the Spodoptera frugiperda Sf9 insect cell line that was designated as Sf-rhabdovirus. Using limiting dilution, this cell line was found to be a mixed population of cells infected by Sf-rhabdovirus variants containing either the full length X accessory gene with a 3.7 kb internal duplication (designated as Sf-rhabdovirus X+3.7) or lacking the duplication and part of the X gene (designated as Sf-rhabdovirus X-), and cells that were negative for Sf-rhabdovirus. In this paper, we found that the Sf-rhabdovirus negative cell clones had sub-populations with different susceptibilities to the replication of Sf-rhabdovirus X+3.7 and X- variants: cell clone Sf9-13F12 was more sensitive to replication by both virus variants compared to Sf9-3003; moreover, Sf9-3003 showed more resistance to X+3.7 replication than to X- replication. RNA-Seq analysis indicated significant differentially expressed genes in the Sf9-13F12 and Sf9-3003 cell clones further supporting that distinct sub-populations of virus-negative cells co-exist in the parent Sf9 cell line.

Development of transgenic Paulownia trees expressing antimicrobial thionin genes for enhanced resistance to fungal infections using chitosan nanoparticles

There is an increasing focus on genetically altering Paulownia trees to enhance their resistance against fungal infections, given their rapid growth and quality wood production. The aim of this research was to establish a technique for incorporating two antimicrobial thionin genes, namely thionin-60 (thio-60) and thionin-63 (thio-63), into Paulownia tomentosa and Paulownia hybrid 9501 through the utilization of chitosan nanoparticles. The outcomes revealed the successful gene transfer into Paulownia trees utilizing chitosan nanoparticles. The effectiveness of thionin proteins against plant pathogens Fusarium and Aspergillus was examined, with a specific focus on Fusarium equiseti due to limited available data. In non-transgenic Paulownia species, the leaf weight inhibition percentage varied from 25 to 36 %, whereas in transgenic species, it ranged from 22 to 7 %. In general, Paulownia species expressing thio-60 displayed increased resistance to F. equiseti, while those expressing thio-63 exhibited heightened resistance to A. niger infection. The thionin proteins displayed a strong affinity for the phospholipid bilayer of the fungal cell membrane, demonstrating their capability to disrupt its structure. The transgenic plants created through this technique showed increased resistance to fungal infections. Thionin-60 demonstrated superior antifungal properties in comparison to thio-63, being more effective at disturbing the fungal cell membrane. These findings indicate that thio-60 holds potential as a novel antifungal agent and presents a promising approach for enhancing the antimicrobial traits of genetically modified Paulownia trees.

A Survey of Deep Learning Methods for Estimating the Accuracy of Protein Quaternary Structure Models

The quality prediction of quaternary structure models of a protein complex, in the absence of its true structure, is known as the Estimation of Model Accuracy (EMA). EMA is useful for ranking predicted protein complex structures and using them appropriately in biomedical research, such as protein-protein interaction studies, protein design, and drug discovery. With the advent of more accurate protein complex (multimer) prediction tools, such as AlphaFold2-Multimer and ESMFold, the estimation of the accuracy of protein complex structures has attracted increasing attention. Many deep learning methods have been developed to tackle this problem; however, there is a noticeable absence of a comprehensive overview of these methods to facilitate future development. Addressing this gap, we present a review of deep learning EMA methods for protein complex structures developed in the past several years, analyzing their methodologies, data and feature construction. We also provide a prospective summary of some potential new developments for further improving the accuracy of the EMA methods.

Copper Resistance Mechanism and Copper Response Genes in Corynebacterium crenatum

Heavy metal resistance mechanisms and heavy metal response genes are crucial for microbial utilization in heavy metal remediation. Here, Corynebacterium crenatum was proven to possess good tolerance in resistance to copper. Then, the transcriptomic responses to copper stress were investigated, and the vital pathways and genes involved in copper resistance of C. crenatum were determined. Based on transcriptome analysis results, a total of nine significantly upregulated DEGs related to metal ion transport were selected for further study. Among them, GY20_RS0100790 and GY20_RS0110535 belong to transcription factors, and GY20_RS0110270, GY20_RS0100790, and GY20_RS0110545 belong to copper-binding peptides. The two transcription factors were studied for the function of regulatory gene expression. The three copper-binding peptides were displayed on the C. crenatum surface for a copper adsorption test. Furthermore, the nine related metal ion transport genes were deleted to investigate the effect on growth in copper stress. This investigation provided the basis for utilizing C. crenatum in copper bioremediation.

Novel benzimidazole angiotensin receptor blockers with anti-SARS-CoV-2 activity equipotent to that of nirmatrelvir: computational and enzymatic studies

BACKGROUND

Hypertension worsens outcomes in SARS-CoV-2 patients. Sartans, a type of antihypertensive angiotensin receptor blocker-(ARB), reduce COVID-19 morbidity and mortality by targeting angiotensin-converting enzyme-2 (ACE2). This study aimed to evaluate the antiviral and antihypertensive effects of nirmatrelvir, commercial sartans (candesartan, losartan, and losartan carboxylic (Exp3174)), and newly synthesized sartans (benzimidazole-N-biphenyl carboxyl (ACC519C) and benzimidazole-N-biphenyl tetrazole (ACC519T)), compared to nirmatrelvir, the antiviral component of Paxlovid.

RESEARCH DESIGN AND METHODS

Surface plasmon resonance (SPR) and enzymatic studies assessed drug effects on ACE2. Antiviral abilities were tested with SARS-CoV-2-infected Vero E6 cells, and antihypertensive effects were evaluated using angiotensin II-contracted rabbit iliac arteries.

RESULTS

Benzimidazole-based candesartan and ACC519C showed antiviral activity comparable to nirmatrelvir (95% inhibition). Imidazole-based losartan, Exp3174, and ACC519T were less potent (75%-80% and 50%, respectively), with Exp3174 being the least effective. SPR analysis indicated high sartans-ACE2 binding affinity. Candesartan and nirmatrelvir combined had greater inhibitory and cytopathic effects (3.96%) than individually (6.10% and 5.08%). ACE2 enzymatic assays showed varying effects of novel sartans on ACE2. ACC519T significantly reduced angiotensin II-mediated contraction, unlike nirmatrelvir and ACC519T(2).

CONCLUSION

This study reports the discovery of a new class of benzimidazole-based sartans that significantly inhibit SARS-CoV-2, likely due to their interaction with ACE2.