Protein Structure Modeling with MODELLER

The flavonoid resokaempferol improves insulin secretion from healthy and dysfunctional pancreatic β-cells

BACKGROUND AND PURPOSE

The pharmacology of flavonoids on β-cell function is largely undefined especially in the context of defective secretion of insulin. We sought to identify flavonoids that increased the insulin-secreting function of β-cells and to explore the underlying mechanisms.

EXPERIMENTAL APPROACH

INS-1 β-cells in culture and islets of Langerhans isolated from control and diabetic male rats were used for insulin secretion experiments. Pharmacological and electrophysiological approaches were used for mechanistic studies.

KEY RESULTS

Among a set of flavonoids, exposure of INS-1 β-cells to resokaempferol (ResoK) enhanced glucose-stimulated insulin secretion and therefore we further characterised its activity and its pharmacological mechanism. ResoK glucose-dependently enhanced insulin secretion in INS-1 β-cells and pancreatic islets isolated from rats. Mechanistically, whole cell patch clamp recordings in INS-1 cells showed that ResoK rapidly and dose-dependently enhanced the L-type Ca2+ current whereas it was inactive towards T-type Ca2+ current. Accordingly, pharmacological inhibition of L-type Ca2+ current but not T-type Ca2+ current blocked the effects of ResoK on glucose-stimulated insulin secretion. ResoK was still active on dysfunctional β-cells as it ameliorated glucose-stimulated insulin secretion in glucotoxicity-induced dysfunctional INS-1 cells and in pancreatic islets isolated from diabetic rats.

CONCLUSION AND IMPLICATIONS

ResoK is a glucose-dependent activator of insulin secretion. Our results indicated that the effects of ResoK on insulin secretion involved its capacity to stimulate L-type Ca2+ currents in cultured β-cells. As ResoK was also effective on dysfunctional β-cells, our work provides a new approach to stimulating insulin secretion, using compounds based on the structure of ResoK.

Advancing rational pesticide development against Drosophila suzukii: bioinformatics tools and applications-a systematic review

CONTEXT

Drosophila suzukii (Matsumura, 1931) is a widespread agricultural pest responsible for significant damage to various soft-skinned fruit hosts. The revolutionary potential of bioinformatics in agriculture emerges from its ability to provide extensive information on pests, fungi, chemical resistance, implications of non-target species, and other critical aspects. This wealth of information allows researchers to engage in projects and applied research in diverse agricultural domains that face these challenges. In this context, bioinformatics tools play a fundamental role. The negative impact of pests on crops, resulting in substantial economic losses, has highlighted the importance of in silico methods.

METHODS

To achieve this, we conducted a systematic search in scientific databases using as keywords "Drosophila suzukii," "biopesticides," "simulations computational," and "in-silico." After applying the filters of relevance and publication date, we organized the articles and prioritized those that directly addressed that matched the keywords and the use of bioinformatics tools. Additionally, we included studies focusing on in silico assays of biopesticides, such as molecular docking. Our review aimed to present a collection of recent literature on biopesticides against Drosophila suzukii, emphasizing bioinformatics methods. Through this work, we strive to contribute to the literature of new perspectives on the development and efficiency of biopesticides, along with to advance research that may improve pest control strategies.

RESULTS

In the results of the systematic review, we found 2734 articles related to the selected keywords. Six of these articles directly address Drosophila suzukii and the use of bioinformatics tools in the search for alternatives in pest control. In the selected studies, we observed that two articles tend to focus on phylogenetic approaches, searching for gene sequences, amino acids, and constructing phylogenetic trees. The other three articles used molecular modeling and docking of receptors such as GABA and TRP with plant-derived and synthetic compounds to study intermolecular interactions. However, we identified gaps in these studies that could lead to further research in the biorational development of biopesticides using bioinformatics tools.

AlphaFold-SFA: Accelerated sampling of cryptic pocket opening, protein-ligand binding and allostery by AlphaFold, slow feature analysis and metadynamics

Sampling rare events in proteins is crucial for comprehending complex phenomena like cryptic pocket opening, where transient structural changes expose new binding sites. Understanding these rare events also sheds light on protein-ligand binding and allosteric communications, where distant site interactions influence protein function. Traditional unbiased molecular dynamics simulations often fail to sample such rare events, as the free energy barrier between metastable states is large relative to the thermal energy. This renders these events inaccessible on the timescales typically simulated by unbiased molecular dynamics, limiting our understanding of these critical processes. In this paper, we proposed a novel unsupervised learning approach termed as slow feature analysis (SFA) which aims to extract slowly varying features from high-dimensional temporal data. SFA trained on small unbiased molecular dynamics simulations launched from AlphaFold generated conformational ensembles manages to capture rare events governing cryptic pocket opening, protein-ligand binding, and allosteric communications in a kinase. Metadynamics simulations using SFA as collective variables manage to sample 'deep' cryptic pocket opening within a few hundreds of nanoseconds which was beyond the reach of microsecond long unbiased molecular dynamics simulations. SFA augmented metadynamics also managed to capture conformational plasticity of protein upon ligand binding/unbinding and provided novel insights into allosteric communication in receptor-interacting protein kinase 2 (RIPK2) which dictates protein-protein interaction. Taken together, our results show how SFA acts as a dimensionality reduction tool which bridges the gap between AlphaFold, molecular dynamics simulation and metadynamics in context of capturing rare events in biomolecules, extending the scope of structure-based drug discovery in the era of AlphaFold.

Computational analysis of human gut microbial prolyl oligopeptidases (POPs) reveal candidate genes as therapeutics for celiac disease

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

Phosphorylation at the disordered N-end makes HuR accumulate and dimerize in the cytoplasm

Human antigen R (HuR) is an RNA binding protein mainly involved in maintaining the stability and controlling the translation of mRNAs, critical for immune response, cell survival, proliferation and apoptosis. Although HuR is a nuclear protein, its mRNA translational-related function occurs at the cytoplasm, where the oligomeric form of HuR is more abundant. However, the regulation of nucleo-cytoplasmic transport of HuR and its connection with protein oligomerization remain unclear. In this work, we describe the phosphorylation of Tyr5 as a new hallmark for HuR activation. Our biophysical, structural and computational assays using phosphorylated and phosphomimetic HuR proteins demonstrate that phosphorylation of Tyr5 at the disordered N-end stretch induces global changes on HuR dynamics and conformation, modifying the solvent accessible surface of the HuR nucleo-cytoplasmic shuttling (HNS) sequence and releasing regions implicated in HuR dimerization. These findings explain the preferential cytoplasmic accumulation of phosphorylated HuR in HeLa cells, aiding to comprehend the mechanisms underlying HuR nucleus-cytoplasm shuttling and its later dimerization, both of which are relevant in HuR-related pathogenesis.

Deciphering the Role of Trehalose in Chroococcidiopsis sp. 029's High-Desiccation Resistance: Sequence Determination, Structural Modelling and Simulative Analysis of the 30S Ribosomal Subunit

Desert strains of the genus Chroococcidiopsis are among the most desiccation-resistant cyanobacteria capable of anhydrobiosis. The accumulation of two sugars, sucrose and trehalose, facilitates the entrance of anhydrobiotes into a reversible state of dormancy by stabilizing cellular components upon water removal. This study aimed to evaluate, at the atomistic level, the role of trehalose in desiccation resistance by using as a model system the 30S ribosomal subunit of the desert cyanobacterium Chroococcidiopsis sp. 029. Molecular dynamic simulations provided atomistic evidence regarding its protective role on the 30S molecular structure. Trehalose forms an enveloping shell around the ribosomal subunit and stabilizes the structures through a network of direct interactions. The simulation confirmed that trehalose actively interacts with the 30S ribosomal subunit and that, by replacing water molecules, it ensures ribosomal structural integrity during desiccation, thus enabling protein synthesis to be carried out upon rehydration.

Molecular basis and functional consequences of the interaction between the base excision repair DNA glycosylase NEIL1 and RPA

NEIL1 is a DNA glycosylase that recognizes and initiates base excision repair of oxidized bases. The ubiquitous ssDNA binding scaffolding protein, replication protein A (RPA), modulates NEIL1 activity in a manner that depends on DNA structure. Interaction between NEIL1 and RPA has been reported, but the molecular basis of this interaction has yet to be investigated. Using a combination of NMR spectroscopy and isothermal titration calorimetry (ITC), we show that NEIL1 interacts with RPA through two contact points. An interaction with the RPA32C protein recruitment domain was mapped to a motif in the common interaction domain (CID) of NEIL1 and a dissociation constant (Kd) of 200 nM was measured. A substantially weaker secondary interaction with the tandem RPA70AB ssDNA binding domains was also mapped to the CID. Together these two contact points reveal NEIL1 has a high overall affinity (Kd ∼ 20 nM) for RPA. A homology model of the complex of RPA32C with the NEIL1 RPA binding motif in the CID was generated and used to design a set of mutations in NEIL1 to disrupt the interaction, which was confirmed by ITC. The mutant NEIL1 remains catalytically active against a thymine glycol lesion in duplex DNA in vitro. Testing the functional effect of disrupting the NEIL1-RPA interaction in vivo using a Fluorescence Multiplex-Host Cell Reactivation (FM-HCR) reporter assay revealed an unexpected role for NEIL1 in nucleotide excision repair. These findings are discussed in the context of the role of NEIL1 in replication-associated repair.

Melodia: A Python Library for Protein Structure Analysis

SUMMARY

Analysing protein structure similarities is an important step in protein engineering and drug discovery. Methodologies that are more advanced than simple RMSD are available but often require extensive mathematical or computational knowledge for implementation. Grouping and optimising such tools in an efficient open-source library increases accessibility and encourages the adoption of more advanced metrics. Melodia is a Python library with a complete set of components devised for describing, comparing and analysing the shape of protein structures using differential geometry of three-dimensional curves and knot theory. It can generate robust geometric descriptors for thousands of shapes in just a few minutes. Those descriptors are more sensitive to structural feature variation than RMSD deviation. Melodia also incorporates sequence structural annotation and three-dimensional visualisations.

AVAILABILITY AND IMPLEMENTATION

Melodia is an open-source Python library freely available on https://github.com/rwmontalvao/Melodia_py, along with interactive Jupyter Notebook tutorials.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Critical amino acid residues regulating TRPA1 Zn2+ response: A comparative study across species

Cellular zinc ions (Zn2+) are crucial for signal transduction in various cell types. The transient receptor potential (TRP) ankyrin 1 (TRPA1) channel, known for its sensitivity to intracellular Zn2+ ([Zn2+]i), has been a subject of limited understanding regarding its molecular mechanism. Here, we used metal ion-affinity prediction, three-dimensional structural modeling, and mutagenesis, utilizing data from the Protein Data Bank and AlphaFold database, to elucidate the [Zn2+]i binding domain (IZD) structure composed by specific AAs residues in human (hTRPA1) and chicken TRPA1 (gTRPA1). External Zn2+ induced activation in hTRPA1, while not in gTRPA1. Moreover, external Zn2+ elevated [Zn2+]i specifically in hTRPA1. Notably, both hTRPA1 and gTRPA1 exhibited inherent sensitivity to [Zn2+]i, as evidenced by their activation upon internal Zn2+ application. The critical AAs within IZDs, specifically histidine at 983/984, lysine at 711/717, tyrosine at 714/720, and glutamate at 987/988 in IZD1, and H983/H984, tryptophan at 710/716, E854/E855, and glutamine at 979/980 in IZD2, were identified in hTRPA1/gTRPA1. Furthermore, mutations, such as the substitution of arginine at 919 (R919) to H919, abrogated the response to external Zn2+ in hTRPA1. Among single-nucleotide polymorphisms (SNPs) at Y714 and a triple SNP at R919 in hTRPA1, we revealed that the Zn2+ responses were attenuated in mutants carrying the Y714 and R919 substitution to asparagine and proline, respectively. Overall, this study unveils the intrinsic sensitivity of hTRPA1 and gTRPA1 to [Zn2+]i mediated through IZDs. Furthermore, our findings suggest that specific SNP mutations can alter the responsiveness of hTRPA1 to extracellular and intracellular Zn2+.

Insights into the Dissociation Process and Binding Pattern of the BRCT7/8-PHF8 Complex

DNA topoisomerase 2-binding protein 1 (Topbp1) plays a crucial role in activating the ataxia-telangiectasia mutated and rad3-related (ATR) complex to initiate DNA damage repair responses. For this process to occur, it is necessary for PHF8 to dissociate from Topbp1. Topbp1 binds to the acidic patch sequence (APS) of PHF8 through its C-terminal BRCT7/8 domain, and disrupting this interaction could be a promising strategy for cancer treatment. To investigate the dissociation process and binding pattern of BRCT7/8-PHF8, we employed enhanced sampling techniques, such as steered molecular dynamics (SMD) simulations and accelerated molecular dynamics (aMD) simulations, along with self-organizing maps (SOM) and time-resolved force distribution analysis (TRFDA) methodologies. Our results demonstrate that the dissociation of PHF8 from BRCT7/8 starts from the N-terminus, leading to the unfolding of the N-terminal helix. Additionally, we identified critical residues that play a pivotal role in this dissociation process. These findings provide valuable insights into the disassociation of PHF8 from BRCT7/8, which could potentially guide the development of novel drugs targeting Topbp1 for cancer therapy.

Antibacterial activities of Miang extracts against selected pathogens and the potential of the tannin-free extracts in the growth inhibition of Streptococcus mutans

Bacterial pathogens have remained a major public health concern for several decades. This study investigated the antibacterial activities of Miang extracts (at non-neutral and neutral pH) against Bacillus cereus TISTR 747, Escherichia coli ATCC 22595, Salmonella enterica serovar Typhimurium TISTR 292 and Streptococcus mutans DMST 18777. The potential of Polyvinylpolypyrrolidone (PVPP)-precipitated tannin-free Miang extracts in growth-inhibition of the cariogenic Streptococcus mutans DMST 18777 and its biofilms was also evaluated. The tannin-rich fermented extracts had the best bacterial growth inhibition against S. mutans DMST 18777 with an MIC of 0.29 and 0.72 mg/mL for nonfilamentous fungi (NFP) Miang and filamentous-fungi-processed (FFP) Miang respectively. This observed anti-streptococcal activity still remained after PVPP-mediated precipitation of bioactive tannins especially, in NFP and FFP Miang. Characterization of the PVPP-treated extracts using High performance liquid chromatography quadrupole-time of flight-mass spectrometry (HPLC-QToF-MS) analysis, also offered an insight into probable compound classes responsible for the activities. In addition, Crystal violet-staining also showed better IC50 values for NFP Miang (4.30 ± 0.66 mg/mL) and FFP Miang (12.73 ± 0.11 mg/mL) against S. mutans DMST 18777 biofilms in vitro. Homology modeling and molecular docking analysis using HPLC-MS identified ligands in tannin-free Miang supernatants, was performed against modelled S. mutans DMST 18777 sortase A enzyme. The in silico analysis suggested that the inhibition by NFP and FFP Miang might be attributed to the presence of ellagic acid, flavonoid aglycones, and glycosides. Thus, these Miang extracts could be optimized and explored as natural active pharmaceutical ingredients (NAPIs) for applications in oral hygienic products.

BioProtIS: Streamlining protein-ligand interaction pipeline for analysis in genomic and transcriptomic exploration

The identification of protein-ligand interactions plays a pivotal role in elucidating biological processes and discovering potential bioproducts. Harnessing the capabilities of computational methods in drug discovery, we introduce an innovative Inverted Virtual Screening (IVS) pipeline. This pipeline Integrated molecular dynamics and docking analyses to ensure that protein structures are not only energetically favorable but also representative of stable conformations. The primary objective of this pipeline is to automate and streamline the analysis of protein-ligand interactions at both genomic and transcriptomic scales. In the contemporary post-genomic era, high-throughput computational screening for bioproducts, biological systems, and therapeutic drugs has become a cornerstone practice. This approach offers the promise of cost-effectiveness, time efficiency, and optimization of laboratory work. Nevertheless, a notable deficiency persists in the availability of efficient pipelines capable of automating the virtual screening process, seamlessly integrating input and output, and leveraging the full potential of open-source tools. To bridge this critical gap, we have developed a versatile pipeline known as BioProtIS. This tool seamlessly integrates a suite of state-of-the-art tools, including Modeller, AlphaFold, Gromacs, FPOCKET, and AutoDock Vina, thus facilitating the streamlined docking of ligands with an expansive repertoire of proteins sourced from genomes and transcriptomes, and substrates. To assess the pipeline's performance, we employed the transcriptomes of Cereus jamacaru (a cactus species) and Aspisoma lineatum (firefly), along with the genome of Homo sapiens. This integration not only improves the accuracy of ligand-protein interactions by minimizing replicability deviations but also optimizes the discovery process by enabling the simultaneous evaluation of multiple substrates. Furthermore, our pipeline accommodates distinct testing scenarios, such as blind docking or site-specific targeting, which are invaluable in applications ranging from drug repositioning to the exploration of new allosteric binding sites and toxicity assessments. BioProtIS has been designed with modularity at its core. This inherent flexibility empowers users to make custom modifications directly within the source code, tailoring the pipeline to their specific research needs. Moreover, it lays the foundation for seamless integration of diverse docking algorithms in future iterations, promising ongoing advancements in the field of computational biology. This pipeline is available for free distribution and can be download at: https://github.com/BBMDO/BioProtIS.

Exploring the anti-gout potential of sunflower receptacles alkaloids: A computational and pharmacological analysis

Gout, a painful condition marked by elevated uric acid levels often linked to the diet's high purine and alcohol content, finds a potential treatment target in xanthine oxidase (XO), a crucial enzyme for uric acid production. This study explores the therapeutic properties of alkaloids extracted from sunflower (Helianthus annuus L.) receptacles against gout. By leveraging computational chemistry and introducing a novel R-based clustering algorithm, "TriDimensional Hierarchical Fingerprint Clustering with Tanimoto Representative Selection (3DHFC-TRS)," we assessed 231 alkaloid molecules from sunflower receptacles. Our clustering analysis pinpointed six alkaloids with significant gout-targeting potential, particularly emphasizing the fifth cluster's XO inhibition capabilities. Through molecular docking and the BatchDTA prediction model, we identified three top compounds-2-naphthylalanine, medroxalol, and fenspiride-with the highest XO affinity. Further molecular dynamics simulations assessed their enzyme active site interactions and binding free energies, employing MM-PBSA calculations. This investigation not only highlights the discovery of promising compounds within sunflower receptacle alkaloids via LC-MS but also introduces medroxalol as a novel gout treatment candidate, showcasing the synergy of computational techniques and LC-MS in drug discovery.

Structural insights of Labeo catla (catla) myxovirus resistance protein,GTP binding recognition and constitutive expression induced with Poly I:C

The Myxovirus resistance (Mx) proteins are critical effectors belonging to the super-family of guanidine triphosphatase, often stimulated by type I interferon (IFN) and mediates antiviral responses to restrict the replication of numerous viral genes in fishes. In teleosts, Mx proteins display diverse and complicated antiviral activity in different species. The present investigation seeks to characterize the Mx gene from Labeo catla upon induction by double-stranded (ds) RNA, polyinosinic-polycytidylic acid, (poly I: C). Molecular modeling and all-atoms molecular dynamics (MD) simulations were employed to understand the architecture of the GTPase domain and its plausible mode of GTP recognition in Mx protein. The full-length L. catla Mx (LcMx) gene sequence (1821 bp nucleotides) encodes an open reading frame of 606 amino acids. Domain search indicated conserved tripartite domain architecture of LcMx and forms a major cluster with the Mx from other teleosts. The positively charged Arginine and polar Glutamine residues from helix 3 and 4 of stalk region LcMx aid in homo-oligomerization. MD simulation portrayed the role of conserved critical residues aid in GTP recognition by the GTPase domain which perfectly corroborates with experimental findings and prior MD studies. After injection of poly I:C, the temporal mRNA profile showed that LcMx expression was significantly elevated in the spleen, brain, kidney, liver, muscle, heart, intestine, and gill tissues. Collectively, these results suggest that the elevated expression of the major innate immune defense gene Mx was able to inhibit the poly I: C mediated virulence in fish.Communicated by Ramaswamy H. Sarma.

State-of-the-art in engineering small molecule biosensors and their applications in metabolic engineering

Genetically encoded biosensors are crucial for enhancing our understanding of how molecules regulate biological systems. Small molecule biosensors, in particular, help us understand the interaction between chemicals and biological processes. They also accelerate metabolic engineering by increasing screening throughput and eliminating the need for sample preparation through traditional chemical analysis. Additionally, they offer significantly higher spatial and temporal resolution in cellular analyte measurements. In this review, we discuss recent progress in in vivo biosensors and control systems-biosensor-based controllers-for metabolic engineering. We also specifically explore protein-based biosensors that utilize less commonly exploited signaling mechanisms, such as protein stability and induced degradation, compared to more prevalent transcription factor and allosteric regulation mechanism. We propose that these lesser-used mechanisms will be significant for engineering eukaryotic systems and slower-growing prokaryotic systems where protein turnover may facilitate more rapid and reliable measurement and regulation of the current cellular state. Lastly, we emphasize the utilization of cutting-edge and state-of-the-art techniques in the development of protein-based biosensors, achieved through rational design, directed evolution, and collaborative approaches.

GPR101: Modeling a constitutively active receptor linked to X-linked acrogigantism

GPR101 is a G protein-coupled receptor (GPCR) implicated in a rare form of genetic gigantism known as X-linked acrogigantism, or X-LAG. In particular, X-LAG patients harbor microduplications in the long arm of the X-chromosome that invariably include the GPR101 gene. Duplications of the GPR101 gene lead to the formation of a new chromatin domain that causes over-expression of the receptor in the pituitary tumors of the patients. Notably, GPR101 is a constitutively active receptor, which stimulates cells to produce the second messenger cyclic AMP (cAMP) in the absence of ligands. Moreover, GPR101 was recently reported to constitutively activate not only the cAMP pathway via Gs, but also other G protein subunits (Gq/11 and G12/13). Hence, chemicals that block the constitutive activity of GPR101, known as inverse agonists, have the potential to be useful for the development of pharmacological tools for the treatment of X-LAG. In this study, we provide structural insights into the putative structure of GPR101 based on in-house built homology models, as well as third party models based on the machine learning methods AlphaFold and AlphaFold-Multistate. Moreover, we report a molecular dynamics study, meant to further probe the constitutive activity of GPR101. Finally, we provide a structural comparison with the closest GPCRs, which suggests that GPR101 does not share their natural ligands. While this manuscript was under review, cryo-electron microscopy structures of GPR101 were reported. These structures are expected to enable computer-aided ligand discovery efforts targeting GPR101.

Delivery to, and Reactivation of, the p53 Pathway in Cancer Cells Using a Grafted Cyclotide Conjugated with a Cell-Penetrating Peptide

Peptides are promising drug modalities that can modulate protein-protein interactions, but their application is hampered by their limited ability to reach intracellular targets. Here, we improved the cytosolic delivery of a peptide blocking p53:MDM2/X interactions using a cyclotide as a stabilizing scaffold. We applied several design strategies to improve intracellular delivery and found that the conjugation of the lead cyclotide to the cyclic cell-penetrating peptide cR10 was the most effective. Conjugation allowed cell internalization at micromolar concentration and led to elevated intracellular p53 levels in A549, MCF7, and MCF10A cells, as well as inducing apoptosis in A549 cells without causing membrane disruption. The lead peptide had >35-fold improvement in inhibitory activity and increased cellular uptake compared to a previously reported cyclotide p53 activator. In summary, we demonstrated the delivery of a large polar cyclic peptide in the cytosol and confirmed its ability to modulate intracellular protein-protein interactions involved in cancer.

Impact of PmrB mutations on clinical Klebsiella pneumoniae with variable colistin-susceptibilities: Structural insights and potent therapeutic solutions

Carbapenem-resistant Klebsiella pneumoniae (CRKP) infections continue to impose high morbidity threats to hospitalized patients worldwide, limiting therapeutic options to last-resort antibiotics like colistin. However, the dynamic genomic landscape of colistin-resistant K. pneumoniae (COLR-Kp) invoked ardent exploration of underlying molecular signatures for therapeutic propositions/designs. We unveiled the structural impact of the widespread and emerging PmrB mutations involved in colistin resistance (COLR) in K. pneumoniae. In the present study, clinical isolates of K. pneumoniae expressed variable susceptibilities to colistin (>0.5 μg/mL for resistant and ≤0.25 μg/mL for susceptible) despite mutations such as T157P, G207D and T246A. The protein sequences extracted from in-house sequenced genomes were used to model mutant PmrB proteins and analyze the underlying structural alterations. The mutations were contrasted based on molecular dynamics simulation trajectories, free-energy landscapes and structural flexibility profiles. The altered backbone flexibilities can be an essential factor for mutant selection by COLR K. pneumoniae and can provide clues to deal with emerging mutants. Furthermore, PmrB having high druggability confidence (>0.99), was explored as a potential target for 1396 virtually screened FDA-approved drug candidates. Among the top-10 compounds (scores >70), amphotericin B was found to be potential candidate with high affinity (Binding energy <-8 kcal/mol) and stable interactions (RMSF <0.7 Å) against PmrB druggable pockets, despite the mutations, which encourages future adjunct therapeutic research against COLR-Kp.

Effect of lactic acid bacteria fermentation on cow milk allergenicity and antigenicity: A review

Cow milk is a major allergenic food. The potential prevention and treatment effects of lactic acid bacteria (LAB)-fermented dairy products on allergic symptoms have garnered considerable attention. Cow milk allergy (CMA) is mainly attributed to extracellular and/or cell envelope proteolytic enzymes with hydrolysis specificity. Numerous studies have demonstrated that LAB prevents the risk of allergies by modulating the development and regulation of the host immune system. Specifically, LAB and its effectors can enhance intestinal barrier function and affect immune cells by interfering with humoral and cellular immunity. Fermentation hydrolysis of allergenic epitopes is considered the main mechanism of reducing CMA. This article reviews the linear epitopes of allergens in cow milk and the effect of LAB on these allergens and provides insight into the means of predicting allergenic epitopes by conventional laboratory analysis methods combined with molecular simulation. Although LAB can reduce CMA in several ways, the mechanism of action remains partially clarified. Therefore, this review additionally attempts to summarize the main mechanism of LAB fermentation to provide guidance for establishing an effective preventive and treatment method for CMA and serve as a reference for the screening, research, and application of LAB-based intervention.

Structural and physical features that distinguish tumor-controlling from inactive cancer neoepitopes

Neoepitopes arising from amino acid substitutions due to single nucleotide polymorphisms are targets of T cell immune responses to cancer and are of significant interest in the development of cancer vaccines. However, understanding the characteristics of rare protective neoepitopes that truly control tumor growth has been a challenge, due to their scarcity as well as the challenge of verifying true, neoepitope-dependent tumor control in humans. Taking advantage of recent work in mouse models that circumvented these challenges, here, we compared the structural and physical properties of neoepitopes that range from fully protective to immunologically inactive. As neoepitopes are derived from self-peptides that can induce immune tolerance, we studied not only how the various neoepitopes differ from each other but also from their wild-type counterparts. We identified multiple features associated with protection, including features that describe how neoepitopes differ from self as well as features associated with recognition by diverse T cell receptor repertoires. We demonstrate both the promise and limitations of neoepitope structural analysis and predictive modeling and illustrate important aspects that can be incorporated into neoepitope prediction pipelines.

Dual function of Zika virus NS2B-NS3 protease

Zika virus (ZIKV) serine protease, indispensable for viral polyprotein processing and replication, is composed of the membrane-anchored NS2B polypeptide and the N-terminal domain of the NS3 polypeptide (NS3pro). The C-terminal domain of the NS3 polypeptide (NS3hel) is necessary for helicase activity and contains an ATP-binding site. We discovered that ZIKV NS2B-NS3pro binds single-stranded RNA with a Kd of ~0.3 μM, suggesting a novel function. We tested various structural modifications of NS2B-NS3pro and observed that constructs stabilized in the recently discovered "super-open" conformation do not bind RNA. Likewise, stabilizing NS2B-NS3pro in the "closed" (proteolytically active) conformation using substrate inhibitors abolished RNA binding. We posit that RNA binding occurs when ZIKV NS2B-NS3pro adopts the "open" conformation, which we modeled using highly homologous dengue NS2B-NS3pro crystallized in the open conformation. We identified two positively charged fork-like structures present only in the open conformation of NS3pro. These forks are conserved across Flaviviridae family and could be aligned with the positively charged grove on NS3hel, providing a contiguous binding surface for the negative RNA strand exiting helicase. We propose a "reverse inchworm" model for a tightly intertwined NS2B-NS3 helicase-protease machinery, which suggests that NS2B-NS3pro cycles between open and super-open conformations to bind and release RNA enabling long-range NS3hel processivity. The transition to the closed conformation, likely induced by the substrate, enables the classical protease activity of NS2B-NS3pro.

A multi-layered computational structural genomics approach enhances domain-specific interpretation of Kleefstra syndrome variants in EHMT1

This study investigates the functional significance of assorted variants of uncertain significance (VUS) in euchromatic histone lysine methyltransferase 1 (EHMT1), which is critical for early development and normal physiology. EHMT1 mutations cause Kleefstra syndrome and are linked to various human cancers. However, accurate functional interpretations of these variants are yet to be made, limiting diagnoses and future research. To overcome this, we integrate conventional tools for variant calling with computational biophysics and biochemistry to conduct multi-layered mechanistic analyses of the SET catalytic domain of EHMT1, which is critical for this protein function. We use molecular mechanics and molecular dynamics (MD)-based metrics to analyze the SET domain structure and functional motions resulting from 97 Kleefstra syndrome missense variants within the domain. Our approach allows us to classify the variants in a mechanistic manner into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). Our findings reveal that the damaging variants are mostly mapped around the active site, substrate binding site, and pre-SET regions. Overall, we report an improvement for this method over conventional tools for variant interpretation and simultaneously provide a molecular mechanism for variant dysfunction.

A Multi-Layered Computational Structural Genomics Approach Enhances Domain-Specific Interpretation of Kleefstra Syndrome Variants in EHMT1

This study investigates the functional significance of assorted variants of uncertain significance (VUS) in euchromatic histone lysine methyltransferase 1 (EHMT1), which is critical for early development and normal physiology. EHMT1 mutations cause Kleefstra syndrome and are linked to various human cancers. However, accurate functional interpretation of these variants are yet to be made, limiting diagnoses and future research. To overcome this, we integrate conventional tools for variant calling with computational biophysics and biochemistry to conduct multi-layered mechanistic analyses of the SET catalytic domain of EHMT1, which is critical for this protein function. We use molecular mechanics and molecular dynamics (MD)-based metrics to analyze the SET domain structure and functional motions resulting from 97 Kleefstra syndrome missense variants within this domain. Our approach allows us to classify the variants in a mechanistic manner into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). Our findings reveal that the damaging variants are mostly mapped around the active site, substrate binding site, and pre-SET regions. Overall, we report an improvement for this method over conventional tools for variant interpretation and simultaneously provide a molecular mechanism of variant dysfunction.

Integrative Genome-Based Survey of the SARS-CoV-2 Omicron XBB.1.16 Variant

The XBB.1.16 SARS-CoV-2 variant, also known as Arcturus, is a recent descendant lineage of the recombinant XBB (nicknamed Gryphon). Compared to its direct progenitor, XBB.1, XBB.1.16 carries additional spike mutations in key antigenic sites, potentially conferring an ability to evade the immune response compared to other circulating lineages. In this context, we conducted a comprehensive genome-based survey to gain a detailed understanding of the evolution and potential dangers of the XBB.1.16 variant, which became dominant in late June. Genetic data indicates that the XBB.1.16 variant exhibits an evolutionary background with limited diversification, unlike dangerous lineages known for rapid changes. The evolutionary rate of XBB.1.16, which amounts to 3.95 × 10-4 subs/site/year, is slightly slower than that of its direct progenitors, XBB and XBB.1.5, which have been circulating for several months. A Bayesian Skyline Plot reconstruction suggests that the peak of genetic variability was reached in early May 2023, and currently, it is in a plateau phase with a viral population size similar to the levels observed in early March. Structural analyses indicate that, overall, the XBB.1.16 variant does not possess structural characteristics markedly different from those of the parent lineages, and the theoretical affinity for ACE2 does not seem to change among the compared variants. In conclusion, the genetic and structural analyses of SARS-CoV-2 XBB.1.16 do not provide evidence of its exceptional danger or high expansion capability. Detected differences with previous lineages are probably due to genetic drift, which allows the virus constant adaptability to the host, but they are not necessarily connected to a greater danger. Nevertheless, continuous genome-based monitoring is essential for a better understanding of its descendants and other lineages.

The structural characterization and bioactivity assessment of nonspecific lipid transfer protein 1 (nsLTP1) from caraway (Carum carvi) seeds

BACKGROUND

Carum carvi (caraway) of the Apiaceae family has been used in many cultures as a cooking spice and part of the folk medicine. Previous reports primarily focus on the medicinal properties of caraway seed essential oil and the whole seeds extract. However, no effort has been made to study caraway proteins and their potential pharmacological properties, including nonspecific lipid transfer protein (nsLTP), necessitating further research. The current study aimed to characterize nonspecific lipid transfer protein 1 (nsLTP1) from caraway seed, determine its three-dimensional structure, and analyze protein-ligand complex interactions through docking studies. We also evaluated nsLTP1 in vitro cytotoxic effect and antioxidant capacity. Additionally, nsLTP1 thermal- and pH- stability were investigated.

METHODS

Caraway nsLTP1 was purified using two-dimensional chromatography. The complete amino acid sequence of nsLTP1 was achieved by intact protein sequence for the first 20 residues and the overlapping digested peptides. The three-dimensional structure was predicted using MODELLER. Autodock Vina software was employed for docking fatty acids against caraway nsLTP1. Assessment of nsLTP1 cytotoxic activity was achieved by MTS assay, and the Trolox equivalent antioxidant capacity (TAC) was determined. Thermal and pH stability of the nsLTP1 was examined by circular dichroism (CD) spectroscopy.

RESULTS

Caraway nsLTP1 is composed of 91 residues and weighs 9652 Da. The three-dimensional structure of caraway nsLTP1 sequence was constructed based on searching known structures in the PDB. We chose nsLTP of Solanum melongena (PDB ID: 5TVI) as the modeling template with the highest identity among all other homologous proteins. Docking linolenic acid with caraway protein showed a maximum binding score of -3.6 kcal/mol. A preliminary screening of caraway nsLTP1 suppressed the proliferation of human breast cancer cell lines MDA-MB-231 and MCF-7 in a dose‑dependent manner with an IC₅₀ value of 52.93 and 44.76 μM, respectively. Also, nsLTP1 (41.4 μM) showed TAC up to 750.4 μM Trolox equivalent. Assessment of nsLTP1 demonstrated high thermal/pH stability.

CONCLUSION

To the best of our knowledge, this is the first study carried out on nsLTP1 from caraway seeds. We hereby report the sequence of nsLTP1 from caraway seeds and its possible interaction with respective fatty acids using in silico approach. Our data indicated that the protein had anticancer and antioxidant activities and was thermally stable.

SARS-CoV-2 Recombinants: Genomic Comparison between XBF and Its Parental Lineages

Recombination events are very common and represent one of the primary drivers of RNA virus evolution. The XBF SARS-CoV-2 lineage is one of the most recently generated recombinants during the COVID-19 pandemic. It is a recombinant of BA.5.2.3 and BA.2.75.3, both descendants of lineages that caused many concerns (BA.5 and BA.2.75, respectively). Here, we performed a genomic survey focused on comparing the recombinant XBF with its parental lineages to provide a comprehensive assessment of the evolutionary potential, epidemiological trajectory, and potential risks. Genetic analyses indicated that although XBF initially showed the typical expansion depicted by a steep curve, causing several concerns, currently there is no indication of significant expansion potential or a contagion rate surpassing that of other currently active or previously prevalent lineages. BSP indicated that the peak has been reached around 19 October 2022 and then the genetic variability suffered slight oscillations until early 5 March 2023 when the population size reduced for the last time starting its last plateau that is still lasting. Structural analyses confirmed its reduced potential, also indicating that properties of NTDs and RBDs of XBF and its parental lineages present no significant difference. Of course, cautionary measures must still be taken and genome-based monitoring remains the best tool for detecting any important changes in viral genome composition.

Khosta: A Genetic and Structural Point of View of the Forgotten Virus

Bats are well-known to be natural reservoirs of various zoonotic coronaviruses, which have caused outbreaks of severe acute respiratory syndrome (SARS) and the COVID-19 pandemic in 2002 and 2019, respectively. In late 2020, two new Sarbecoviruses were found in Russia, isolated in Rhinolophus bats, i.e., Khosta-1 in R. ferrumequinum and Khosta-2 in R. hipposideros. The potential danger associated with these new species of Sarbecovirus is that Khosta-2 has been found to interact with the same entry receptor as SARS-CoV-2. Our multidisciplinary approach in this study demonstrates that Khosta-1 and -2 currently appear to be not dangerous with low risk of spillover, as confirmed by prevalence data and by phylogenomic reconstruction. In addition, the interaction between Khosta-1 and -2 with ACE2 appears weak, and furin cleavage sites are absent. While the possibility of a spillover event cannot be entirely excluded, it is currently highly unlikely. This research further emphasizes the importance of assessing the zoonotic potential of widely distributed batborne CoV in order to monitor changes in genomic composition of viruses and prevent spillover events (if any).

Molecular dynamic simulation reveals the molecular interactions of epidermal growth factor receptor with musk xylene are involved in the carcinogenicity

Musk xylene (MX), a kind of personal care product, has become a new type of environmental contaminant in recent years. Long-term exposure to MX is associated with a variety of cancers, but the mechanism is still unclear. Meanwhile, our previous research showed that MX exposure could lead to malignant transformation of human liver cells L02 and up-regulation of multi genes which are involved in the MAPK signaling pathway, such as the epidermal growth factor receptor (EGFR). These findings indicated that the MAPK signaling pathway might be involved in the malignant transformation caused by MX, but the mechanism is also unclear. In this study, the underlying interaction mechanisms between EGFR and MX were investigated using molecular dynamics (MD) simulation. Results revealed that MX bound to the ECD of EGFR in four binding sites, which was mainly driven by van der Waals and nonpolar interactions, and the affinity of MX toward ECD was sIII > sI > sII > sIV. Further analysis through MD simulation found that s III, the site with the strongest binding, was coincidentally located at the binding area of EGF, which is the natural ligand of EGFR. Therefore, we speculated that MX may activate the MAPK signaling pathway by binding to EGFR in a similar way to EGF, and finally lead to tumorigenesis. In addition, the MM/PBSA method could also be utilized to calculate the hot residues in each binding site. The prediction of hot residues would provide some theoretical guidance for further study of the carcinogenesis mechanisms of MX both in MD simulation and experimental research.

SARS-CoV CH.1.1 Variant: Genomic and Structural Insight

In early February 2023, the Omicron subvariant XBB.1.5, also known as "Kraken", accounted for more than 44% of new COVID-19 cases worldwide, whereas a relatively new Omicron subvariant named CH.1.1, deemed "Orthrus", accounted for less than 6% of new COVID-19 cases during the subsequent weeks. This emerging variant carries a mutation, L452R, previously observed in the highly pathogenic Delta and the highly transmissible BA.4 and BA.5 variants, necessitating a shift to active surveillance to assure adequate preparedness for likely future epidemic peaks. We provide a preliminary understanding of the global distribution of this emerging SARS-CoV-2 variant by combining genomic data with structural molecular modeling. In addition, we shield light on the number of specific point mutations in this lineage that may have functional significance, thereby increasing the risk of disease severity, vaccine resistance, and increased transmission. This variant shared about 73% of the mutations with Omicron-like strains. Our homology modeling analysis revealed that CH.1.1 may have a weakened interaction with ACE2 and that its electrostatic potential surface appears to be more positive than that of the reference ancestral virus. Finally, our phylogenetic analysis revealed that this likely-emerging variant was already cryptically circulating in European countries prior to its first detection, highlighting the importance of having access to whole genome sequences for detecting and controlling emerging viral strains.

Translational protein RpsE as an alternative target for novel nucleoside analogues to treat MDR Enterobacter cloacae ATCC 13047: network analysis and molecular dynamics study

The pathogenic Enterobacter cloacae subsp. cloacae str. ATCC 13047 has contemporarily emerged as a multi-drug resistant strain. To formulate an effective treatment option, alternative therapeutic methods need to be explored. The present study focused on Gene Interaction Network study of 46 antimicrobial resistance genes to reveal the densely interconnecting and functional hub genes in E. cloacae ATCC 13047. The AMR genes were subjected to clustering, topological and functional enrichment analysis, revealing rpsE (RpsE), acrA (AcrA) and arnT (ArnT) as novel therapeutic drug targets for hindering drug resistance in the pathogenic strain. Network topology further indicated translational protein RpsE to be exploited as a promising drug-target candidate for which the structure was predicted, optimized and validated through molecular dynamics simulations (MDS). Absorption, distribution, metabolism and excretion screening recognized ZINC5441082 (N-Isopentyladenosine) (Lead_1) and ZINC1319816 (cyclopentyl-aminopurinyl-hydroxymethyl-oxolanediol) (Lead_2) as orally bioavailable compounds against RpsE. Molecular docking and MDS confirmed the binding efficacy and protein-ligand complex stability. Furthermore, binding free energy (G_bind) calculations, principal component and free energy landscape analyses affirmed the predicted nucleoside analogues against RpsE protein to be comprehensively examined as a potential treatment strategy against E. cloacae ATCC 13047.

Production, homology modeling and mutagenesis studies on GlcH glucose transporter from Prochlorococcus sp. strain SS120

The marine cyanobacterium Prochlorococcus is one of the main primary producers on Earth, which can take up glucose by using the high affinity, multiphasic transporter GlcH. We report here the overexpression of glcH from Prochlorococcus marinus strain SS120 in Escherichia coli. Modeling studies of GlcH using the homologous MelB melibiose transporter from Salmonella enterica serovar Typhimurium showed high conservation at the overall fold. We observed that an important structural interaction, mediated by a strong hydrogen bond between D8 and R141, is conserved in Prochlorococcus, although the corresponding amino acids in MelB from Salmonella are different. Biased docking studies suggested that when glucose reaches the pocket of the transporter and interacts with D8 and R141, the hydrogen bond network in which these residues are involved could be disrupted, favoring a conformational change with the subsequent translocation of the glucose molecule towards the cytoplasmic region of the pmGlcH structure. Based on these theoretical predictions and on the conservation of N117 and W348 in other MelB structures, D8, N117, R141 and W348 were mutated to glycine residues. Their key role in glucose transport was evaluated by glucose uptake assays. N117G and W348G mutations led to 17 % decrease in glucose uptake, while D8G and R141G decreased the glucose transport by 66 % and 92 % respectively. Overall, our studies provide insights into the Prochlorococcus 3D-structure of GlcH, paving the way for further analysis to understand the features which are involved in the high affinity and multiphasic kinetics of this transporter.

Molecular In-Depth on the Epidemiological Expansion of SARS-CoV-2 XBB.1.5

Since the beginning of the pandemic, the generation of new variants periodically recurs. The XBB.1.5 SARS-CoV-2 variant is one of the most recent. This research was aimed at verifying the potential hazard of this new subvariant. To achieve this objective, we performed a genome-based integrative approach, integrating results from genetic variability/phylodynamics with structural and immunoinformatic analyses to obtain as comprehensive a viewpoint as possible. The Bayesian Skyline Plot (BSP) shows that the viral population size reached the plateau phase on 24 November 2022, and the number of lineages peaked at the same time. The evolutionary rate is relatively low, amounting to 6.9 × 10−4 subs/sites/years. The NTD domain is identical for XBB.1 and XBB.1.5 whereas their RBDs only differ for the mutations at position 486, where the Phe (in the original Wuhan) is replaced by a Ser in XBB and XBB.1, and by a Pro in XBB.1.5. The variant XBB.1.5 seems to spread more slowly than sub-variants that have caused concerns in 2022. The multidisciplinary molecular in-depth analyses on XBB.1.5 performed here does not provide evidence for a particularly high risk of viral expansion. Results indicate that XBB.1.5 does not possess features to become a new, global, public health threat. As of now, in its current molecular make-up, XBB.1.5 does not represent the most dangerous variant.

A unique network of attack, defence and competence on the outer membrane of the periodontitis pathogen Tannerella forsythia

Periodontopathogenic Tannerella forsythia uniquely secretes six peptidases of disparate catalytic classes and families that operate as virulence factors during infection of the gums, the KLIKK-peptidases. Their coding genes are immediately downstream of novel ORFs encoding the 98-132 residue potempins (Pot) A, B1, B2, C, D and E. These are outer-membrane-anchored lipoproteins that specifically and potently inhibit the respective downstream peptidase through stable complexes that protect the outer membrane of T. forsythia, as shown in vivo. Remarkably, PotA also contributes to bacterial fitness in vivo and specifically inhibits matrix metallopeptidase (MMP) 12, a major defence component of oral macrophages, thus featuring a novel and highly-specific physiological MMP inhibitor. Information from 11 structures and high-confidence homology models showed that the potempins are distinct β-barrels with either a five-stranded OB-fold (PotA, PotC and PotD) or an eight-stranded up-and-down fold (PotE, PotB1 and PotB2), which are novel for peptidase inhibitors. Particular loops insert like wedges into the active-site cleft of the genetically-linked peptidases to specifically block them either via a new "bilobal" or the classic "standard" mechanism of inhibition. These results discover a unique, tightly-regulated proteolytic armamentarium for virulence and competence, the KLIKK-peptidase/potempin system.

Benchmarking TriadAb using targets from the second antibody modeling assessment

Computational modeling and design of antibodies has become an integral part of today's research and development in antibody therapeutics. Here we describe the Triad Antibody Homology Modeling (TriadAb) package, a functionality of the Triad protein design platform that predicts the structure of any heavy and light chain sequences of an antibody Fv domain using template-based modeling. To gauge the performance of TriadAb, we benchmarked against the results of the Second Antibody Modeling Assessment (AMA-II). On average, TriadAb produced main-chain carbonyl root-mean-square deviations between models and experimentally determined structures at 1.10 Å, 1.45 Å, 1.41 Å, 3.04 Å, 1.47 Å, 1.27 Å, 1.63 Å in the framework and the six complementarity-determining regions (H1, H2, H3, L1, L2, L3), respectively. The inaugural results are comparable to those reported in AMA-II, corroborating with our internal bench-based experiences that models generated using TriadAb are sufficiently accurate and useful for antibody engineering using the sequence design capabilities provided by Triad.

Y12F mutation in Pseudomonas plecoglossicida S7 lipase enhances its thermal and pH stability for industrial applications: a combination of in silico and in vitro study

Appropriate amino acid substitutions are critical for protein engineering to redesign catalytic properties of industrially important enzymes like lipases. The present study aimed for improving the environmental stability of lipase from Pseudomonas plecoglossicida S7 through site-directed mutagenesis driven by computational studies. lipA gene was amplified and sequenced. Both wild type (WT) and mutant type (MT) lipase genes were expressed into the pET SUMO system. The expressed proteins were purified and characterized for pH and thermostability. The lipase gene belonged to subfamily I.1 lipase. Molecular dynamics revealed that Y12F-palmitic acid complex had a greater binding affinity (-6.3 Kcal/mol) than WT (-6.0 Kcal/mol) complex. Interestingly, MDS showed that the binding affinity of WT-complex (-130.314 ± 15.11 KJ/mol) was more than mutant complex (-108.405 ± 69.376 KJ/mol) with a marked increase in the electrostatic energy of mutant (-26.969 ± 12.646 KJ/mol) as compared to WT (-15.082 ± 13.802 KJ/mol). Y12F mutant yielded 1.27 folds increase in lipase activity at 55 °C as compared to the purified WT protein. Also, Y12F mutant showed increased activity (~ 1.2 folds each) at both pH 6 and 10. P. plecoglossicida S7. Y12F mutation altered the kinetic parameters of MT (K_m- 1.38 mM, V_max- 22.32 µM/min) as compared to WT (K_m- 1.52 mM, V_max- 29.76 µM/min) thus increasing the binding affinity of mutant lipase. Y12F mutant lipase with better pH and thermal stability can be used in biocatalysis.

Molecular dynamic simulation and functional analysis of pathogenic PTEN mutations in glioblastoma

PTEN, a dual-phosphatase and scaffold protein, is one of the most commonly mutated tumour suppressor gene across various cancer types in human. The aim of this study therefore was to investigate the stability, structural and functional effects, and pathogenicity of 12 missense PTEN mutations (R15S, E18G, G36R, N49I, Y68H, I101T, C105F, D109N, V133I, C136Y, R173C and N276S) found by next generation sequencing of the PTEN gene in tissue samples obtained from glioblastoma patients. Computational tools and molecular dynamic simulation programs were used to identify the deleterious effects of these mutations. Furthermore, PTEN mRNA and protein expression levels were evaluated by qRT-PCR, Western Blot, and immunohistochemistry staining methods. Various computational tools predicted strong deleterious effects for the G36R, C105F, C136Y and N276S mutations. Molecular dynamic simulation revealed a significant decrease in protein stability for the Y68H and N276S mutations when compared with the wild type protein; whereas, C105F, D109N, V133I and R173C showed partial stability reduction. Significant residual fluctuations were observed in the R15S, N49I and C136Y mutations and radius of gyration graphs revealed the most compact structure for D109N and least for C136Y. In summary, our study is the first one to show the presence of PTEN E18G, N49I, D109N and N276S mutations in glioblastoma patients; where, D109N is neutral and N276S is a damaging and disease-associated mutation.Communicated by Ramaswamy H. Sarma.

Theoretical 3D Modeling of NLRP3 Inflammasome Complex

The NOD-like receptor pyrin domain containing 3 (NLRP3) is a multidomain protein that plays a key role in innate immune response. Structures of NLRP3 in different conformational states and bound to cognate partners are available. In this chapter we present an approach to model the oligomeric structure of NLRP3 by homology modeling using multiple templates, symmetry, and refinement. The overall process presented here represents advanced exercise in structural modeling that provides unique insights into the biological role and activation of NLRP3 oligomer. Finally, the same approach can be easily adapted to the rest of the members of the NLRP family.

Protein Homology Modeling for Effective Drug Design

The effective drug design, especially for combating the multi-drug-resistant bacterial pathogens, requires more and more sophisticated procedures to obtain novel lead-like compounds. New classes of enzymes should be explored, especially those that help bacteria overcome existing treatments. The homology modeling is useful in obtaining the models of new enzymes; however, the active sites of them are sometimes present in closed conformations in the crystal structures, not suitable for drug design purposes. In such difficult cases, the combination of homology modeling, molecular dynamics simulations, and fragment screening can give satisfactory results.

The PANDORA Software for Anchor-Restrained Peptide:MHC Modeling

Major histocompatibility complexes (MHC) play a key role in the immune surveillance system in all jawed vertebrates. MHC class I molecules randomly sample cytosolic peptides from inside the cell, while MHC class II sample exogenous peptides. Both types of peptide:MHC complex are then presented on the cell surface for recognition by αβ T cells (CD8+ and CD4+, respectively). The three-dimensional structure of such complexes can give crucial insights in the presentation and recognition mechanisms. For this reason, softwares like PANDORA have been developed to rapidly and accurately generate peptide:MHC (pMHC) 3D structures. In this chapter, we describe the protocol of PANDORA. PANDORA exploits the structural knowledge on anchor pockets that MHC molecules use to dock peptides. PANDORA provides anchor positions as restraints to guide the modeling process. This allows PANDORA to generate twenty 3D models in just about 5 min. PANDORA is highly customizable, easy to install, supports parallel processing, and is suitable to provide large datasets for deep learning algorithms.

Specificities of Protein Homology Modeling for Allosteric Drug Design

The allosteric binding sites are usually located in the flexible areas of proteins, which are hardly visible in the crystal structures. However, there are notable exceptions like allosteric sites in receptors in class B and C of GPCRs, which are located within a well-defined bundle of transmembrane helices. Class B and C evolved from class A and even after swapping of orthosteric and allosteric sites the central binding site persisted and it can be used for easy design of allosteric drugs. However, studying the ligand binding to the allosteric sites in the most populated class A of GPCRs is still a challenge, since they are located mostly in unresolved parts of the receptor's structure, and especially N-terminus. This chapter provides an example of cannabinoid CB1 receptor N-terminal homology modeling, ligand-guided modeling of the allosteric site in GABA receptor, as well as C-linker modeling in the potassium ion channels where the allosteric phospholipid ligand PIP2 is bound.

Assessing the effects of β-triketone herbicides on HPPD from environmental bacteria using a combination of in silico and microbiological approaches

4-hydroxyphenylpyruvate dioxygenase (HPPD) is the molecular target of β-triketone herbicides in plants. This enzyme, involved in the tyrosine pathway, is also present in a wide range of living organisms, including microorganisms. Previous studies, focusing on a few strains and using high herbicide concentrations, showed that β-triketones are able to inhibit microbial HPPD. Here, we measured the effect of agronomical doses of β-triketone herbicides on soil bacterial strains. The HPPD activity of six bacterial strains was tested with 1× or 10× the recommended field dose of the herbicide sulcotrione. The selected strains were tested with 0.01× to 15× the recommended field dose of sulcotrione, mesotrione, and tembotrione. Molecular docking was also used to measure and model the binding mode of the three herbicides with the different bacterial HPPD. Our results show that responses to herbicides are strain-dependent with Pseudomonas fluorescens F113 HPPD activity not inhibited by any of the herbicide tested, when all three β-triketone herbicides inhibited HPPD in Bacillus cereus ATCC14579 and Shewanella oneidensis MR-1. These responses are also molecule-dependent with tembotrione harboring the strongest inhibitory effect. Molecular docking also reveals different binding potentials. This is the first time that the inhibitory effect of β-triketone herbicides is tested on environmental strains at agronomical doses, showing a potential effect of these molecules on the HPPD enzymatic activity of non-target microorganisms. The whole-cell assay developed in this study, coupled with molecular docking analysis, appears as an interesting way to have a first idea of the effect of herbicides on microbial communities, prior to setting up microcosm or even field experiments. This methodology could then largely be applied to other family of pesticides also targeting an enzyme present in microorganisms.

Genetic and Structural Data on the SARS-CoV-2 Omicron BQ.1 Variant Reveal Its Low Potential for Epidemiological Expansion

The BQ.1 SARS-CoV-2 variant, also known as Cerberus, is one of the most recent Omicron descendant lineages. Compared to its direct progenitor BA.5, BQ.1 has some additional spike mutations in some key antigenic sites, which confer further immune escape ability over other circulating lineages. In such a context, here, we perform a genome-based survey aimed at obtaining a complete-as-possible nuance of this rapidly evolving Omicron subvariant. Genetic data suggest that BQ.1 represents an evolutionary blind background, lacking the rapid diversification that is typical of a dangerous lineage. Indeed, the evolutionary rate of BQ.1 is very similar to that of BA.5 (7.6 × 10-4 and 7 × 10-4 subs/site/year, respectively), which has been circulating for several months. The Bayesian Skyline Plot reconstruction indicates a low level of genetic variability, suggesting that the peak was reached around 3 September 2022. Concerning the affinity for ACE2, structure analyses (also performed by comparing the properties of BQ.1 and BA.5 RBD) indicate that the impact of the BQ.1 mutations may be modest. Likewise, immunoinformatic analyses showed moderate differences between the BQ.1 and BA5 potential B-cell epitopes. In conclusion, genetic and structural analyses on SARS-CoV-2 BQ.1 suggest no evidence of a particularly dangerous or high expansion capability. Genome-based monitoring must continue uninterrupted for a better understanding of its descendants and all other lineages.

Potent Activation of Human but Not Mouse TRPA1 by JT010

Transient receptor potential (TRP) ankyrin repeat 1 (TRPA1), which is involved in inflammatory pain sensation, is activated by endogenous factors, such as intracellular Zn²⁺ and hydrogen peroxide, and by irritant chemical compounds. The synthetic compound JT010 potently and selectively activates human TRPA1 (hTRPA1) among the TRPs. Therefore, JT010 is a useful tool for analyzing TRPA1 functions in biological systems. Here, we show that JT010 is a potent activator of hTRPA1, but not mouse TRPA1 (mTRPA1) in human embryonic kidney (HEK) cells expressing hTRPA1 and mTRPA1. Application of 0.3-100 nM of JT010 to HEK cells with hTRPA1 induced large Ca²⁺ responses. However, in HEK cells with mTRPA1, the response was small. In contrast, both TRPA1s were effectively activated by allyl isothiocyanate (AITC) at 10-100 μM. Similar selective activation of hTRPA1 by JT010 was observed in electrophysiological experiments. Additionally, JT010 activated TRPA1 in human fibroblast-like synoviocytes with inflammation, but not TRPA1 in mouse dorsal root ganglion cells. As cysteine at 621 (C621) of hTRPA1, a critical cysteine for interaction with JT010, is conserved in mTRPA1, we applied JT010 to HEK cells with mutations in mTRPA1, where the different residue of mTRPA1 with tyrosine at 60 (Y60), with histidine at 1023 (H1023), and with asparagine at 1027 (N1027) were substituted with cysteine in hTRPA1. However, these mutants showed low sensitivity to JT010. In contrast, the mutation of hTRPA1 at position 669 from phenylalanine to methionine (F669M), comprising methionine at 670 in mTRPA1 (M670), significantly reduced the response to JT010. Moreover, the double mutant at S669 and M670 of mTRPA1 to S669E and M670F, respectively, induced slight but substantial sensitivity to 30 and 100 nM JT010. Taken together, our findings demonstrate that JT010 potently and selectively activates hTRPA1 but not mTRPA1.

Phenotype-based drug screening: An in vivo strategy to classify and identify the chemical compounds modulating zebrafish M-cell regeneration

Several disease-modulatory FDA-approved drugs are being used in patients with neurodegenerative diseases. However, information on their toxicity-related profiles is very limited. Therefore, measurement of drug toxicity is essential to increase the knowledge of their side effects. This study aimed to identify compounds that can modulate M-cell regeneration by causing neuro-protection and -toxicity. Here, we developed a simple and efficient in vivo assay using Tg (hsp: Gal4FF62A; UAS: nfsB-mCherry) transgenic zebrafish larvae. Interestingly, via the phenotype-based drug screening approach, we rapidly investigated 1,260 compounds from the United States drug collection and validated these in large numbers, including 14 compounds, that were obstructing this regeneration process. Next, 4 FDA-approved drugs out of 14 compounds were selected as the lead hits for in silico analysis to clarify their binding patterns with PTEN and SOCS3 signaling due to their significant potential in the inhibition of axon regeneration. Molecular docking studies indicated good binding affinity of all 4 drugs with the respective signaling molecules. This may point to PTEN and SOCS3 as the signaling molecules responsible for reducing axon regeneration. Moreover, the acute effect of compounds in reducing M-cell regeneration delineated their toxic effect. In conclusion, our in vivo along with in silico screening strategy will promote the rapid translation of new therapeutics to improve knowledge of the toxicity profile of approved/non-approved drugs efficiently.

In Vitro Rescue of the Bile Acid Transport Function of ABCB11 Variants by CFTR Potentiators

ABCB11 is responsible for biliary bile acid secretion at the canalicular membrane of hepatocytes. Variations in the ABCB11 gene cause a spectrum of rare liver diseases. The most severe form is progressive familial intrahepatic cholestasis type 2 (PFIC2). Current medical treatments have limited efficacy. Here, we report the in vitro study of Abcb11 missense variants identified in PFIC2 patients and their functional rescue using cystic fibrosis transmembrane conductance regulator potentiators. Three ABCB11 disease-causing variations identified in PFIC2 patients (i.e., A257V, T463I and G562D) were reproduced in a plasmid encoding an Abcb11-green fluorescent protein. After transfection, the expression and localization of the variants were studied in HepG2 cells. Taurocholate transport activity and the effect of potentiators were studied in Madin–Darby canine kidney (MDCK) clones coexpressing Abcb11 and the sodium taurocholate cotransporting polypeptide (Ntcp/Slc10A1). As predicted using three-dimensional structure analysis, the three variants were expressed at the canalicular membrane but showed a defective function. Ivacaftor, GLP1837, SBC040 and SBC219 potentiators increased the bile acid transport of A257V and T463I and to a lesser extent, of G562D Abcb11 missense variants. In addition, a synergic effect was observed when ivacaftor was combined with SBC040 or SBC219. Such potentiators could represent new pharmacological approaches for improving the condition of patients with ABCB11 deficiency due to missense variations affecting the function of the transporter.

Unraveling key interactions and the mechanism of demethylation during hAGT-mediated DNA repair via simulations

Alkylating agents pose the biggest threat to the genomic integrity of cells by damaging DNA bases through regular alkylation. Such damages are repaired by several automated types of machinery inside the cell. O6-alkylguanine-DNA alkyltransferase (AGT) is an enzyme that performs the direct repair of an alkylated guanine base by transferring the alkyl group to a cysteine residue. In the present study, using extensive MD simulations and hybrid QM/MM calculations, we have investigated the key interactions between the DNA lesion and the hAGT enzyme and elucidated the mechanisms of the demethylation of the guanine base. Our simulation shows that the DNA lesion is electrostatically stabilized by the enzyme and the Arg135 of hAGT enzyme provides the main driving force to flip the damaged base into the enzyme. The QM/MM calculations show demethylation of the damaged base as a three-step process in a thermodynamically feasible and irreversible manner. Our calculations show that the final product forms via Tyr114 in a facile way in contrast to the previously proposed Lys-mediated route.

Identification of key amino acid residues toward improving the catalytic activity and substrate specificity of plant-derived cytochrome P450 monooxygenases CYP716A subfamily enzyme for triterpenoid production in Saccharomyces cerevisiae

Triterpenoids constitute a group of specialized plant metabolites with wide structural diversity and high therapeutic value for human health. Cytochrome P450 monooxygenases (CYP) are a family of enzymes important for generating the structural diversity of triterpenoids by catalyzing the site-specific oxidization of the triterpene backbone. The CYP716 enzyme family has been isolated from various plant families as triterpenoid oxidases; however, their experimental crystal structures are not yet available and the detailed catalytic mechanism remains elusive. Here, we address this challenge by integrating bioinformatics approaches with data from other CYP families. Medicago truncatula CYP716A12, the first functionally characterized CYP716A subfamily enzyme, was chosen as the model for this study. We performed homology modeling, structural alignment, in silico site-directed mutagenesis, and molecular docking analysis to search and screen key amino acid residues relevant to the catalytic activity and substrate specificity of the CYP716A subfamily enzyme in triterpenoid biosynthesis. An in vivo functional analysis using engineered yeast that endogenously produced plant-derived triterpenes was performed to elucidate the results. When the amino acids in the signature region and substrate recognition sites (SRSs) were substituted, the product profile of CYP716A12 was modified. We identified amino acid residues that control the substrate contraction of the enzyme (D292) and engineered the enzyme to improve its catalytic activity and substrate specificity (D122, I212, and Q358) for triterpenoid biosynthesis. In addition, we demonstrated the versatility of this strategy by changing the properties of key residues in SRSs to improve the catalytic activity of Arabidopsis thaliana CYP716A1 (S356) and CYP716A2 (M206, F210) at C-28 on the triterpene backbone. This research has the potential to help in the production of desired triterpenoids in engineered yeast by increasing the catalytic activity and substrate specificity of plant CYP716A subfamily enzymes.

Structure based design of effective HtpG-derived vaccine antigens against M. tuberculosis

Vaccine development against Tuberculosis is a strong need, given the low efficacy of the sole vaccine hitherto used, the Bacillus Calmette–Guérin (BCG) vaccine. The chaperone-like protein HtpG_Mtb of M. tuberculosis is a large dimeric and multi-domain protein with promising antigenic properties. We here used biophysical and biochemical studies to improve our understanding of the structural basis of HtpG_Mtb functional role and immunogenicity, a precious information to engineer improved antigens. We showed that HtpG_Mtb is a dimeric nucleotide-binding protein and identified the dimerisation interface on the C-terminal domain of the protein. We also showed that the most immunoreactive regions of the molecule are located on the C-terminal and middle domains of the protein, whereas no role is played by the catalytic N-terminal domain in the elicitation of the immune response. Based on these observations, we experimentally validated our predictions in mice, using a plethora of immunological assays. As an outcome, we designed vaccine antigens with enhanced biophysical properties and ease of production, albeit conserved or enhanced antigenic properties. Our results prove the efficacy of structural vaccinology approaches in improving our understanding of the structural basis of immunogenicity, a precious information to engineer more stable, homogeneous, efficiently produced, and effective vaccine antigens.

Transcriptomic responses of the zearalenone (ZEN)-detoxifying yeast Apiotrichum mycotoxinivorans to ZEN exposure

Zearalenone (ZEN) is a potent oestrogenic mycotoxin that is mainly produced by Fusarium species and is a serious environmental pollutant in animal feeds. Apiotrichum mycotoxinivorans has been widely used as a feed additive to detoxify ZEN. However, the effects of ZEN on A. mycotoxinivorans and its detoxification mechanisms remain unclear. In this study, transcriptomic and bioinformatic analyses were used to investigate the molecular responses of A. mycotoxinivorans to ZEN exposure and the genetic basis of ZEN detoxification. We detected 1424 significantly differentially expressed genes (DEGs), of which 446 were upregulated and 978 were downregulated. Functional and enrichment analyses showed that ZEN-induced genes were significantly associated with xenobiotic metabolism, oxidative stress response, and active transport systems. However, ZEN-inhibited genes were mainly related to cell division, cell cycle, and fungal development. Subsequently, bioinformatic analysis identified candidate ZEN-detoxification enzymes. The Baeyer-Villiger monooxygenases and carboxylesterases, which are responsible for the formation and subsequent hydrolysis of a new ZEN lactone, respectively, were significantly upregulated. In addition, the expression levels of genes related to conjugation and transport involved in the xenobiotic detoxification pathway were significantly upregulated. Moreover, the expression levels of genes encoding enzymatic antioxidants and those related to growth and apoptosis were significantly upregulated and downregulated, respectively, which made it possible for A. mycotoxinivorans to survive in a highly toxic environment and efficiently detoxify ZEN. This is the first systematic report of ZEN tolerance and detoxification in A. mycotoxinivorans. We identified the metabolic enzymes that were potentially involved in detoxifying ZEN in the GMU1709 strain and found that ZEN-induced transcriptional regulation of genes is key to withstanding highly toxic environments. Hence, our results provide valuable information for developing enzymatic detoxification systems or engineering this detoxification pathway in other species.

Characterization of SARS-CoV-2 Escape Mutants to a Pair of Neutralizing Antibodies Targeting the RBD and the NTD

Mutations in the spike protein of SARS-CoV-2 can lead to evasion from neutralizing antibodies and affect the efficacy of passive and active immunization strategies. Immunization of mice harboring an entire set of human immunoglobulin variable region gene segments allowed to identify nine neutralizing monoclonal antibodies, which either belong to a cluster of clonally related RBD or NTD binding antibodies. To better understand the genetic barrier to emergence of SARS-CoV-2 variants resistant to these antibodies, escape mutants were selected in cell culture to one antibody from each cluster and a combination of the two antibodies. Three independently derived escape mutants to the RBD antibody harbored mutations in the RBD at the position T478 or S477. These mutations impaired the binding of the RBD antibodies to the spike protein and conferred resistance in a pseudotype neutralization assay. Although the binding of the NTD cluster antibodies were not affected by the RBD mutations, the RBD mutations also reduced the neutralization efficacy of the NTD cluster antibodies. The mutations found in the escape variants to the NTD antibody conferred resistance to the NTD, but not to the RBD cluster antibodies. A variant resistant to both antibodies was more difficult to select and only emerged after longer passages and higher inoculation volumes. VOC carrying the same mutations as the ones identified in the escape variants were also resistant to neutralization. This study further underlines the rapid emergence of escape mutants to neutralizing monoclonal antibodies in cell culture and indicates the need for thorough investigation of escape mutations to select the most potent combination of monoclonal antibodies for clinical use.