MIB2: metal ion-binding site prediction and modeling server

Strontium Attenuates LPS-Induced Inflammation via TLR4/MyD88/NF-κB Pathway in Bovine Ruminal Epithelial Cells

Subacute ruminal acidosis (SARA) is a common nutritional metabolic disease in ruminants that causes significant economic losses to dairy farming. Strontium (Sr) is known to be involved in bone metabolism and exhibits potent anti-inflammatory effects. To evaluate the effect of Sr on inflammation in bovine ruminal epithelial cells, a model of LPS-induced inflammation was established in this study, and the cell viability of bovine ruminal epithelial cells was measured using CCK-8. The production of pro-inflammatory cytokines was measured by ELISA and real-time PCR, respectively. The related proteins of the TLR4/MyD88/NF-κB pathway were assayed through Western blotting, and the fluorescence of p-p65 and p-IκB were assayed by immunofluorescence. Molecular docking of Sr and TLR4/MyD88/NF-κB pathway-related proteins was performed using MIB2 ( http://bioinfo.cmu.edu.tw/MIB2/ ). Results showed that after treatment for 24 h, the cell viability was decreased at the high concentration of Sr (≥ 10 mmol/L). Sr significantly decreased the production of TNF-α, IL-1β, and IL-6, downregulated the related proteins expression of the TLR4/MyD88/NF-κB pathway, and reduced the fluorescence levels of p-p65 and p-IκB. The NF-κB pathway inhibitor PDTC and molecular docking further revealed that Sr reduced LPS-induced pro-inflammatory cytokines production via the TLR4/MyD88/NF-κB pathway. These results suggest that Sr reduces LPS-induced pro-inflammatory cytokines production via the TLR4/MyD88/NF-κB pathway, thereby exerting an anti-inflammatory effect in bovine ruminal epithelial cells, providing a basis for Sr in the treatment of bovine rumen acidosis disease.

Bioinspired mp20 mimicking uricase in ZIF-8: Metal ion dependent for controllable activity

Mini protein mimicking uricase (mp20) has shown significant potential as a replacement for natural enzymes in the development of uric acid biosensors. However, the design of mp20 has resulted to an inactive form of peptide, causing of loss their catalytic activity. Herein, this paper delineates the impact of various metal cofactors on the catalytic activity of mp20. The metal ion-binding site prediction and docking (MIB) web server was employed to identify the metal ion binding sites and their affinities towards mp20 residues. Among the tested metal ions, Cu2+ displayed the highest docking score, indicating its preference for interaction with Thr16 and Asp17 residues of mp20. To assess the catalytic activity of mp20 in the presence of metal ions, uric acid assays was monitored using a colorimetric method. The presence of Cu2+ in the assays promotes the activation of mp20, resulting in a color change based on quinoid production. Furthermore, the encapsulation of the mp20 within zeolitic imidazolate framework-8 (ZIF-8) notably improved the stability of the biomolecule. In comparison to the naked mp20, the encapsulated ZIFs biocomposite (mp20@ZIF-8) demonstrates superior stability, selectivity and sensitivity. ZIF's porous shells provides excellent protection, broad detection (3-100 μM) with a low limit (4.4 μM), and optimal function across pH (3.4-11.4) and temperature (20-100°C) ranges. Cost-effective and stable mp20@ZIF-8 surpasses native uricase, marking a significant biosensor technology breakthrough. This integration of metal cofactor optimization and robust encapsulation sets new standards for biosensing applications.

Conserved cysteine residues in Kaposi's sarcoma herpesvirus ORF34 are necessary for viral production and viral pre-initiation complex formation

Kaposi's sarcoma herpesvirus (KSHV) ORF34 plays a significant role as a component of the viral pre-initiation complex (vPIC), which is indispensable for late gene expression across beta- and gammaherpesviruses. Although the key role of ORF34 within the vPIC and its function as a hub protein have been recognized, further clarification regarding its specific contribution to vPIC functionality and interactions with other components is required. This study employed a deep learning algorithm-assisted structural model of ORF34, revealing highly conserved amino acid residues across human beta- and gammaherpesviruses localized in structured domains. Thus, we engineered ORF34 alanine-scanning mutants by substituting conserved residues with alanine. These mutants were evaluated for their ability to interact with other vPIC factors and restore viral production in cells harboring the ORF34-deficient KSHV-BAC. Our experimental results highlight the crucial role of the four cysteine residues conserved in ORF34: a tetrahedral arrangement consisting of a pair of C-Xn-C consensus motifs. This suggests the potential incorporation of metal cations in interacting with ORF24 and ORF66 vPIC components, facilitating late gene transcription, and promoting overall virus production by capturing metal cations. In summary, our findings underline the essential role of conserved cysteines in KSHV ORF34 for effective vPIC assembly and viral replication, thereby enhancing our understanding of the complex interplay between the vPIC components.

IMPORTANCE

The initiation of late gene transcription is universally conserved across the beta- and gammaherpesvirus families. This process employs a viral pre-initiation complex (vPIC), which is analogous to a cellular PIC. Although KSHV ORF34 is a critical factor for viral replication and is a component of the vPIC, the specifics of vPIC formation and the essential domains crucial for its function remain unclear. Structural predictions suggest that the four conserved cysteines (C170, C175, C256, and C259) form a tetrahedron that coordinates the metal cation. We investigated the role of these conserved amino acids in interactions with other vPIC components, late gene expression, and virus production to demonstrate for the first time that these cysteines are pivotal for such functions. This discovery not only deepens our comprehensive understanding of ORF34 and vPIC dynamics but also lays the groundwork for more detailed studies on herpesvirus replication mechanisms in future research.

Deciphering retention effect of extracellular polymeric substances to typical heavy metals and their interaction with key inner enzymes of Candidatus Kuenenia

This study employed spectroscopy, metagenomics, and molecular simulation to investigate the inhibitory effects of Cd(II) and Cu(II) on the anammox system, examining both intracellular and extracellular effects. At concentrations of 5 mg/L, Cd(II) and Cu(II) significantly reduced nitrogen removal efficiency by 41.46 % and 62.03 %, respectively. Additionally, elevated metal concentrations were correlated with decreased extracellular polymeric substances (EPS), thereby reducing their capacity to absorb heavy metals, particularly Cu(II), which decreased from 76.47 % to 14.67 %. Spectral analysis revealed alterations in the secondary structures of EPS induced by Cd(II) and Cu(II), decreasing the ratio of extracellular protein α-helix to (β-sheet + random coil), which resulted in looser extracellular protein configurations. The results of the metagenomics study showed that the abundance of Candidatus Kuenenia and its genes encoding nitrogen removal-related enzymes was reduced. The abundance of hzs-γ was reduced by 35.09 % at a concentration of 5 mg/L Cu(II). Conversely, genes associated with metal efflux enzymes, like czcR, increased by 54.86 % at 2 mg/L Cd(II). Molecular docking revealed robust bindings of Cd(II) to HZS-α (-342.299 ± 218.165 kJ/mol) and Cu(II) to HZS-γ (-880.934 ± 55.526 kJ/mol). This study elucidated the inhibitory mechanisms of Cd(II) and Cu(II) on the anammox system, providing insights into the resistance of anammox bacteria to heavy metals.

Weak magnetic field promotes denitrification by stimulating ferromagnetic ion-containing metalloprotein expression

Weak magnetic field (WMF) has been recognized to promote biological denitrification processes; however, the underlying mechanisms remain largely unexplored, hindering the optimization of its effectiveness. Here, we systematically investigated the effects of WMF on denitrification performance, enzyme activity, microbial community, and metaproteome in packed bed bioreactors treating high nitrate wastewater under different WMF intensities and C:N ratios. Results showed that WMFs significantly promoted denitrification by consistently stimulating the activities of denitrifying reductases and NAD+/NADH biosynthesis across decreasing C:N ratios. Reductases and electron transfer enzymes involved in denitrification were overproduced due to the significantly enriched overexpression of ferromagnetic ion-containing (FIC) metalloproteins. We also observed WMFs' intensity-dependent selective pressure on microbial community structures despite the effects being limited compared to those caused by changing C:N ratios. By coupling genome-centric metaproteomics and structure prediction, we found the dominant denitrifier, Halomonas, was outcompeted by Pseudomonas and Azoarcus under WMFs, likely due to its structural deficiencies in iron uptake, suggesting that advantageous ferromagnetic ion acquisition capacity was necessary to satisfy the substrate demand for FIC metalloprotein overproduction. This study advances our understanding of the biomagnetic effects in the context of complex communities and highlights WMF's potential for manipulating FIC protein-associated metabolism and fine-tuning community structure.

An integrated multitool analysis contributes elements to interpreting unclassified factor IX missense variants associated with hemophilia B

BACKGROUND

Dissection of genotype-phenotype relationships in hemophilia B (HB) is particularly relevant for challenging (mild HB) or for HB-associated but unclassified factor (F)IX missense variants.

OBJECTIVE

To contribute elements to interpret unclassified HB-associated FIX missense variants by a multiple-level approach upon identification of a reported, but uncharacterized, FIX missense variant associated with mild HB.

METHODS

Molecular modeling of wild-type and V92A FIX variants, expression studies in HEK293 cells with evaluation of protein (ELISA, western blotting) and activity (activated partial thromboplastin time-based/chromogenic assays) levels after recombinant expression, and multiple prediction tools.

RESULTS

The F9(NM_000133.4):c.275T>C (p.V92A) variant was found in a mild HB patient (antigen, 45.4 U/dL; coagulant activity, 23.6 IU/dL; specific activity, 0.52). Newly generated molecular models showed alterations in Gla/EGF1-EGF2 domain conformation impacting Ca++ affinity and protein-protein interactions with activated factor XI (FXIa). Multitool analysis indicated a moderate impact on protein structure/function of the valine-to-alanine substitution, in accordance with patient and modeling data. Expression studies on the V92A variant showed a specific activity (0.49 ± 0.07; wild-type, 1.0 ± 0.1) recapitulating that of the natural variant, and pointed toward a moderate activation impairment as the main determinant underlying the p.V92A defect. The validated multitool approach, integrated with evidence-based data, was challenged on a panel (n = 9) of unclassified FIX missense variants, which resulted in inferred protein (secretion/function) outputs and HB severity.

CONCLUSION

The rational integration of multitool and multiparameter analyses contributed elements to interpret genotype/phenotype relationships of unclassified FIX missense variants, with implications for diagnosis, management, and treatment of HB patients, and potentially translatable into other human disorders.

Designable immobilization of D-allulose 3-epimerase on bimetallic organic frameworks based on metal ion compatibility for enhanced D-allulose production

D-allulose, a low-calorie rare sugar catalyzed by D-allulose 3-epimerase (DAE), is highly sought after for its potential health benefits. However, poor reusability and stability of DAE limited its popularization in industrial applications. Although metal-organic frameworks (MOFs) offer a promising enzyme platform for enzyme immobilization, developing customized strategies for MOF immobilization of enzymes remains challenging. In this study, we introduce a designable strategy involving the construction of bimetal-organic frameworks (ZnCo-MOF) based on metal ions compatibility. The DAE@MOFs materials were prepared and characterized, and the immobilization of DAE and the enzymatic characteristics of the MOF-immobilized DAE were subsequently evaluated. Remarkably, DAE@ZnCo-MOF exhibited superior recyclability which could maintain 95 % relative activity after 8 consecutive cycles. The storage stability is significantly improved compared to the free form, with a relative activity of 116 % remaining after 30 days. Molecular docking was also employed to investigate the interaction between DAE and the components of MOFs synthesis. The results demonstrate that the DAE@ZnCo-MOF exhibited enhanced catalytic efficiency and increased stability. This study introduces a viable and adaptable MOF-based immobilization strategy for enzymes, which holds the potential to expand the implementation of enzyme biocatalysts in a multitude of disciplines.

Conserved cysteine residues in Kaposi's sarcoma herpesvirus ORF34 are necessary for viral production and viral pre-initiation complex formation

Kaposi's sarcoma herpesvirus (KSHV) ORF34 plays a significant role as a component of the viral pre-initiation complex (vPIC), which is indispensable for late gene expression across beta and gamma herpesviruses. Although the key role of ORF34 within the vPIC and its function as a hub protein have been recognized, further clarification regarding its specific contribution to vPIC functionality and interactions with other components is required. This study employed a deep-learning algorithm-assisted structural model of ORF34, revealing highly conserved amino acid residues across human beta- and gamma-herpesviruses localized in structured domains. Thus, we engineered ORF34 alanine-scanning mutants by substituting conserved residues with alanine. These mutants were evaluated for their ability to interact with other vPIC factors and restore viral production in cells harboring the ORF34-deficient KSHV-BAC. Our experimental results highlight the crucial role of the 4 cysteine residues conserved in ORF34: a tetrahedral arrangement consisting of a pair of C-Xn-C consensus motifs. This suggests the potential incorporation of metal cations in interacting with ORF24 and ORF66 vPIC components, facilitating late gene transcription, and promoting overall virus production by capturing metal cations. In summary, our findings underline the essential role of conserved cysteines in KSHV ORF34 for effective vPIC assembly and viral replication, thereby enhancing our understanding of the complex interplay between the vPIC components.

Explorations on the antiviral potential of zinc and magnesium salts against chikungunya virus: implications for therapeutics

Background

Chikungunya virus (CHIKV), which causes chikungunya fever, is an arbovirus of public health concern with no approved antiviral therapies. A significant proportion of patients develop chronic arthritis after an infection. Zinc and magnesium salts help the immune system respond effectively against viral infections. This study explored the antiviral potential of zinc sulphate, zinc acetate, and magnesium sulphate against CHIKV infection.

Methods

The highest non-toxic concentration of the salts (100 µM) was used to assess the prophylactic, virucidal, and therapeutic anti-CHIKV activities. Dose-dependent antiviral effects were investigated to find out the 50% inhibitory concentration of the salts. Entry bypass assay was conducted to find out whether the salts affect virus entry or post entry stages. Virus output in all these experiments was estimated using a focus-forming unit assay, real-time RT-PCR, and immunofluorescence assay.

Results

Different time- and temperature-dependent assays revealed the therapeutic antiviral activity of zinc and magnesium salts against CHIKV. A minimum exposure of 4 hours and treatment initiation within 1 to 2 hours of infection are required for inhibition of CHIKV. Entry assays revealed that zinc salt affected virus-entry. Entry bypass assays suggested that both salts affected post-entry stages of CHIKV. In infected C57BL6 mice orally fed with zinc and magnesium salts, a reduction in viral RNA copy number was observed.

Conclusion

The study results suggest zinc salts exert anti-CHIKV activity at entry and post entry stages of the virus life cycle, while magnesium salt affect CHIKV at post entry stages. Overall, the study highlights the significant antiviral potential of zinc sulphate, zinc acetate, and magnesium sulphate against CHIKV, which can be exploited in designing potential therapeutic strategies for early treatment of chikungunya patients, thereby reducing the virus-associated persistent arthritis.

Multi-omics reveals new links between Fructosamine-3-Kinase (FN3K) and core metabolic pathways

Fructosamine-3-kinases (FN3Ks) are a conserved family of repair enzymes that phosphorylate reactive sugars attached to lysine residues in peptides and proteins. Although FN3Ks are present across the Tree of Life and share detectable sequence similarity to eukaryotic protein kinases, the biological processes regulated by these kinases are largely unknown. To address this knowledge gap, we leveraged the FN3K CRISPR Knock-Out (KO) HepG2 cell line alongside an integrative multi-omics study combining transcriptomics, metabolomics, and interactomics to place these enzymes in a pathway context. The integrative analyses revealed the enrichment of pathways related to oxidative stress response, lipid biosynthesis (cholesterol and fatty acids), and carbon and co-factor metabolism. Moreover, enrichment of nicotinamide adenine dinucleotide (NAD) binding proteins and localization of human FN3K (HsFN3K) to mitochondria suggests potential links between FN3K and NAD-mediated energy metabolism and redox balance. We report specific binding of HsFN3K to NAD compounds in a metal and concentration-dependent manner and provide insight into their binding mode using modeling and experimental site-directed mutagenesis. Our studies provide a framework for targeting these understudied kinases in diabetic complications and metabolic disorders where redox balance and NAD-dependent metabolic processes are altered.

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+.

Structural basis of substrate recognition and allosteric activation of the proapoptotic mitochondrial HtrA2 protease

The mitochondrial serine protease HtrA2 is a human homolog of the Escherichia coli Deg-proteins exhibiting chaperone and proteolytic roles. HtrA2 is involved in both apoptotic regulation via its ability to degrade inhibitor-of-apoptosis proteins (IAPs), as well as in cellular maintenance as part of the cellular protein quality control machinery, by preventing the possible toxic accumulation of aggregated proteins. In this study, we use advanced solution NMR spectroscopy methods combined with biophysical characterization and biochemical assays to elucidate the crucial role of the substrate recognizing PDZ domain. This domain regulates the protease activity of HtrA2 by triggering an intricate allosteric network involving the regulatory loops of the protease domain. We further show that divalent metal ions can both positively and negatively modulate the activity of HtrA2, leading to a refined model of HtrA2 regulation within the apoptotic pathway.

Application of SUMO fusion technology for the enhancement of stability and activity of lysophospholipase from Pyrococcus abyssi

Heterologous production of proteins in Escherichia coli has raised several challenges including soluble production of target proteins, high levels of expression and purification. Fusion tags can serve as the important tools to overcome these challenges. SUMO (small ubiquitin-related modifier) is one of these tags whose fusion to native protein sequence can enhance its solubility and stability. In current research, a simple, efficient and cost-effective method is being discussed for the construction of pET28a-SUMO vector. In order to improve the stability and activity of lysophospholipase from Pyrococcus abyssi (Pa-LPL), a 6xHis-SUMO tag was fused to N-terminal of Pa-LPL by using pET28a-SUMO vector. Recombinant SUMO-fused enzyme (6 H-S-PaLPL) works optimally at 35 °C and pH 6.5 with remarkable thermostability at 35-95 °C. Thermo-inactivation kinetics of 6 H-S-PaLPL were also studied at 35-95 °C with first order rate constant (kIN) of 5.58 × 10- 2 h-1 and half-life of 12 ± 0 h at 95 °C. Km and Vmax for the hydrolysis of 4-nitrophenyl butyrate were calculated to be 2 ± 0.015 mM and 3882 ± 22.368 U/mg, respectively. 2.4-fold increase in Vmax of Pa-LPL was observed after fusion of 6xHis-SUMO tag to its N-terminal. It is the first report on the utilization of SUMO fusion tag to enhance the overall stability and activity of Pa-LPL. Fusion of 6xHis-SUMO tag not only aided in the purification process but also played a crucial role in increasing the thermostability and activity of the enzyme. SUMO-fused enzyme, thus generated, can serve as an important candidate for degumming of vegetable oils at industrial scale.

Purification and characterization of glutamate dehydrogenase from rainbow trout (Oncorhynchus mykiss) liver and molecular docking studies

Glutamate dehydrogenase (GDH) participates in the energy metabolism of proteins and the synthesis of metabolites important for the organism. In this study, GDH enzyme was purified from the liver of rainbow trout (Oncorhynchus mykiss) by 2',5'-ADP Sepharose 4B affinity chromatography in one step. As a result of this purification process, GDH enzyme was purified 171-fold with 5.83 U/mg protein-specific activity. The characterization experiments presented that the storage stability of the purified GDH enzyme was determined as -80°C; optimum temperature 40°C; it was determined that the optimum ionic strength was 100 mM phosphate buffer and the optimum pH was 8.00. Sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and PAGE studies showed that the natural molar mass of the purified GDH enzyme was 346.74 kDa, and the molar mass of its subunits was 53.71 kDa. Km and Vmax values for substrates and coenzymes of GDH enzyme purified from rainbow trout liver were calculated, and the lowest Km value was found in NAD+ (1.86 mM) and the highest Vmax value in NH4 + (1.79 U/mL). The effects of some metal ions, vitamins, and solvents on the activity of the purified GDH enzyme were investigated and also IC50 values and inhibition types. The metal ion with the lowest IC50 value is Ag+ (8.65 ± 1.68 μM), and the vitamin is B6 (0.77 ± 0.04 mM). The binding affinities of inhibitors were investigated with molecular docking, based on the conformational state of GDH.

Methotrexate for Drug Repurposing as an Anti-Aggregatory Agent to Mercuric Treated α-Chymotrypsinogen-A

Protein aggregation is related to numerous pathological conditions like Alzheimer's and Parkinson's disease. In our study, we have shown that an already existing FDA-approved drug; methotrexate (MTX) can be reprofiled on preformed α-chymotrypsinogen A (α-Cgn A) aggregates. The zymogen showed formation of aggregates upon interaction with mercuric ions, with increasing concentration of Hg2Cl2 (0-150 µM). The hike in ThT and ANS fluorescence concomitant with blue shift, bathochromic shift and the hyperchromic effect in the CR absorbance, RLS and turbidity measurements, substantiate the zymogen β-rich aggregate formation. The secondary structural alterations of α- Cgn A as analyzed by CD measurements, FTIR and Raman spectra showed the transformation of native β-barrel conformation to β-inter-molecular rich aggregates. The native α- Cgn A have about 30% α-helical content which was found to be about 3% in presence of mercuric ions suggesting the formation of aggregates. The amorphous aggregates were visualized by SEM. On incubation of Hg2Cl2 treated α- Cgn A with increasing concentration of the MTX resulted in reversing aggregates to the native-like structure. These results were supported by remarkable decrease in ThT and ANS fluorescence intensities and CR absorbance and also consistent with CD, FTIR, and Raman spectroscopy data. MTX was found to increase the α-helical content of the zymogen from 3 to 15% proposing that drug is efficient in disrupting the β-inter-molecular rich aggregates and reverting it to native like structure. The SEM images are in accordance with CD data showing the disintegration of aggregates. The most effective concentration of the drug was found to be 120 µM. Molecular docking analysis showed that MTX molecule was surrounded by the hydrophobic residues including Phe39, His40, Arg145, Tyr146, Thr151, Gly193, Ser195, and Gly216 and conventional hydrogen bonds, including Gln73 (bond length: 2.67Å), Gly142 (2.59Å), Thr144 (2.81Å), Asn150 (2.73Å), Asp153 (2.71Å), and Cys191 (2.53Å). This investigation will help to find the use of already existing drugs to cure protein misfolding-related abnormalities.

Electrophysiological properties and structural prediction of the SARS-CoV-2 viroprotein E

COVID-19, the infectious disease caused by the most recently discovered coronavirus SARS- CoV-2, has caused millions of sick people and thousands of deaths all over the world. The viral positive-sense single-stranded RNA encodes 31 proteins among which the spike (S) is undoubtedly the best known. Recently, protein E has been reputed as a potential pharmacological target as well. It is essential for the assembly and release of the virions in the cell. Literature describes protein E as a voltage-dependent channel with preference towards monovalent cations whose intracellular expression, though, alters Ca2+ homeostasis and promotes the activation of the proinflammatory cascades. Due to the extremely high sequence identity of SARS-CoV-2 protein E (E-2) with the previously characterized E-1 (i.e., protein E from SARS-CoV) many data obtained for E-1 were simply adapted to the other. Recent solid state NMR structure revealed that the transmembrane domain (TMD) of E-2 self-assembles into a homo-pentamer, albeit the oligomeric status has not been validated with the full-length protein. Prompted by the lack of a common agreement on the proper structural and functional features of E-2, we investigated the specific mechanism/s of pore-gating and the detailed molecular structure of the most cryptic protein of SARS-CoV-2 by means of MD simulations of the E-2 structure and by expressing, refolding and analyzing the electrophysiological activity of the transmembrane moiety of the protein E-2, in its full length. Our results show a clear agreement between experimental and predictive studies and foresee a mechanism of activity based on Ca2+ affinity.

A comprehensive review of protein-centric predictors for biomolecular interactions: from proteins to nucleic acids and beyond

Proteins interact with diverse ligands to perform a large number of biological functions, such as gene expression and signal transduction. Accurate identification of these protein-ligand interactions is crucial to the understanding of molecular mechanisms and the development of new drugs. However, traditional biological experiments are time-consuming and expensive. With the development of high-throughput technologies, an increasing amount of protein data is available. In the past decades, many computational methods have been developed to predict protein-ligand interactions. Here, we review a comprehensive set of over 160 protein-ligand interaction predictors, which cover protein-protein, protein-nucleic acid, protein-peptide and protein-other ligands (nucleotide, heme, ion) interactions. We have carried out a comprehensive analysis of the above four types of predictors from several significant perspectives, including their inputs, feature profiles, models, availability, etc. The current methods primarily rely on protein sequences, especially utilizing evolutionary information. The significant improvement in predictions is attributed to deep learning methods. Additionally, sequence-based pretrained models and structure-based approaches are emerging as new trends.

Kinetic properties of glucose 6-phosphate dehydrogenase and inhibition effects of several metal ions on enzymatic activity in vitro and cells

Due to the non-degradable and persistent nature of metal ions in the environment, they are released into water bodies, where they accumulate in fish. In order to assess pollution in fish, the enzyme, glucose 6-phosphate dehydrogenase (G6PD), has been employed as a biomarker due to sensitivity to various ions. This study investigates the kinetic properties of the G6PD enzyme in yellow catfish (Pelteobagrus fulvidraco), and analyzes the effects of these metal ions on the G6PD enzyme activity in the ovarian cell line (CCO) of channel catfish (Ictalurus punctatus). IC50 values and inhibition types of G6PD were determined in the metal ions Cu2+, Al3+, Zn2+, and Cd2+. While, the inhibition types of Cu2+ and Al3+ were the competitive inhibition, Zn2+ and Cd2+ were the linear mixed noncompetitive and linear mixed competitive, respectively. In vitro experiments revealed an inverse correlation between G6PD activity and metal ion concentration, mRNA levels and enzyme activity of G6PD increased at the lower metal ion concentration and decreased at the higher concentration. Our findings suggest that metal ions pose a significant threat to G6PD activity even at low concentrations, potentially playing a crucial role in the toxicity mechanism of metal ion pollution. This information contributes to the development of a biomonitoring tool for assessing metal ion contamination in aquatic species.

An ancient bacterial zinc acquisition system identified from a cyanobacterial exoproteome

Bacteria have developed fine-tuned responses to cope with potential zinc limitation. The Zur protein is a key player in coordinating this response in most species. Comparative proteomics conducted on the cyanobacterium Anabaena highlighted the more abundant proteins in a zur mutant compared to the wild type. Experimental evidence showed that the exoprotein ZepA mediates zinc uptake. Genomic context of the zepA gene and protein structure prediction provided additional insights on the regulation and putative function of ZepA homologs. Phylogenetic analysis suggests that ZepA represents a primordial system for zinc acquisition that has been conserved for billions of years in a handful of species from distant bacterial lineages. Furthermore, these results show that Zur may have been one of the first regulators of the FUR family to evolve, consistent with the scarcity of zinc in the ecosystems of the Archean eon.

PinMyMetal: A hybrid learning system to accurately model metal binding sites in macromolecules

Metal ions are vital components in many proteins for the inference and engineering of protein function, with coordination complexity linked to structural (4-residue predominate), catalytic (3-residue predominate), or regulatory (2-residue predominate) roles. Computational tools for modeling metal ions in protein structures, especially for transient, reversible, and concentration-dependent regulatory sites, remain immature. We present PinMyMetal (PMM), a sophisticated hybrid machine learning system for predicting zinc ion localization and environment in macromolecular structures. Compared to other predictors, PMM excels in predicting regulatory sites (median deviation of 0.34 Å), demonstrating superior accuracy in locating catalytic sites (median deviation of 0.27 Å) and structural sites (median deviation of 0.14 Å). PMM assigns a certainty score to each predicted site based on local structural and physicochemical features independent of homolog presence. Interactive validation through our server, CheckMyMetal, expands PMM's scope, enabling it to pinpoint and validates diverse functional zinc sites from different structure sources (predicted structures, cryo-EM and crystallography). This facilitates residue-wise assessment and robust metal binding site design. The lightweight PMM system demands minimal computing resources and is available at https://PMM.biocloud.top. While currently trained on zinc, the PMM workflow can easily adapt to other metals through expanded training data.

Unraveling potential mechanism of different metal ions effect on anammox through big data analysis, molecular docking and molecular dynamics simulation

Heavy metals (HMs) have been widely reported to pose an adverse effect on anaerobic ammonia oxidation (anammox) bacteria, yet the underlying mechanisms remain unclear. This study provides new insights into the potential mechanisms of interaction between HMs and functional enzymes through big date analysis, molecular docking and molecular dynamics simulation. The statistical analysis indicated that 10 mg/L Cu(II) and Cd(II) reduced nitrogen removal rate (NRR) by 85% and 43%, while 5 mg/L Fe(II) enhanced NRR by 29%. Additionally, the results of molecular simulations provided a microscopic interpretation for these macroscopic data. Molecular docking revealed that Hg(II) formed a distinctive binding site on ferritin, while other HMs resided at iron oxidation sites. Furthermore, HMs exhibited distinct binding sites on hydrazine dehydrogenase. Concurrently, the molecular dynamics simulation results further substantiated their capacity to form complexes. Cu(II) displayed the strongest binding affinity with ferritin for -1576 ± 79 kJ/mol in binding free energy calculation. Moreover, Cd(II) bound to ferritin and HDH for -1052.67 ± 58.49 kJ/mol, -290.02 ± 49.68 kJ/mol, respectively. This research addressed a crucial knowledge gap, shedding light on potential applications for remediating heavy metal-laden industrial wastewater.

Structural Plasticity Allows Calumenin-1 to Moonlight as a Ca2+-Independent Chaperone: Pb2+ Enables Probing Alternate Inhibitory Conformation

Human calumenin-1 (HsCalu-1) is an endoplasmic reticulum (ER) and Golgi-resident Ca2+-binding protein of the hepta-EF-hand superfamily that plays a vital role in maintaining the cytoplasmic Ca2+ concentration below toxic levels by interacting with Sarco/endoplasmic reticulum Ca2+-ATPase (SERCA) and ryanodine receptors (RyR), indicating its role in Ca2+ homeostasis in the ER. HsCalu-1 seems to be able to exhibit structural plasticity to achieve its plethora of functions. In this study, we demonstrate that HsCalu-1 acts as a chaperone in both its intrinsically disordered state (apo form) and the structured state (Ca2+-bound form). HsCalu-1 chaperone activity is independent of Ca2+ and Pb2+ binding attenuating its chaperone-like activity. Incidentally, Pb2+ binds to HsCalu-1 with lower affinity (KD = 38.46 μM) (compared to Ca2+-binding), leading to the formation of a less-stable conformation as observed by a sharp drop in its melting temperature Tm from 67 °C in the Ca2+-bound form to 43 °C in the presence of Pb2+. The binding site for Pb2+ was mapped as being in the EF-Hand-234 domain of HsCalu-1, a region that overlaps with the Ca2+-dependent initiator of its functional fold. A change in the secondary and tertiary structure, leading to a less-stable but compact conformation upon Pb2+ binding, is the mechanism by which the chaperone-like activity of HsCalu-1 is diminished. Our results not only demonstrate the chaperone activity by a protein in its disordered state but also explain, using Pb2+ as a probe, that the multiple functions of calumenin are due to its ability to adopt a quasi-stable conformation.

M-Ionic: prediction of metal-ion-binding sites from sequence using residue embeddings

MOTIVATION

Understanding metal-protein interaction can provide structural and functional insights into cellular processes. As the number of protein sequences increases, developing fast yet precise computational approaches to predict and annotate metal-binding sites becomes imperative. Quick and resource-efficient pre-trained protein language model (pLM) embeddings have successfully predicted binding sites from protein sequences despite not using structural or evolutionary features (multiple sequence alignments). Using residue-level embeddings from the pLMs, we have developed a sequence-based method (M-Ionic) to identify metal-binding proteins and predict residues involved in metal binding.

RESULTS

On independent validation of recent proteins, M-Ionic reports an area under the curve (AUROC) of 0.83 (recall = 84.6%) in distinguishing metal binding from non-binding proteins compared to AUROC of 0.74 (recall = 61.8%) of the next best method. In addition to comparable performance to the state-of-the-art method for identifying metal-binding residues (Ca2+, Mg2+, Mn2+, Zn2+), M-Ionic provides binding probabilities for six additional ions (i.e. Cu2+, Po43-, So42-, Fe2+, Fe3+, Co2+). We show that the pLM embedding of a single residue contains sufficient information about its neighbours to predict its binding properties.

AVAILABILITY AND IMPLEMENTATION

M-Ionic can be used on your protein of interest using a Google Colab Notebook (https://bit.ly/40FrRbK). The GitHub repository (https://github.com/TeamSundar/m-ionic) contains all code and data.

Exploring the potential of Halalkalibacterium halodurans laccase for endosulfan and chlorophacinone degradation: insights from molecular docking and molecular dynamics simulations

Pesticides are widely used in agriculture but at the same time, a majority of them are known to cause serious harm to health and the environment. In the recent past, laccases have been reported as key enzymes having the ability to degrade pollutants by converting them into less toxic forms. In this investigation, laccase from polyextremophilic bacterium Halalkalibacterium halodurans C-125 was analyzed for its structural, physicochemical, and functional characterization using in silico approaches. The 3D model of the said enzyme is unknown; therefore, the model was generated by template-independent modeling using ROBETTA, I-TASSER, and Alphafold server. The best-generated model from Alphafold with a confidence of 0.95 was validated from ERRAT and Verify 3D scores of 89.95 and 91.80%, respectively. The Ramachandran plot generated using the PROCHECK server further predicted the accuracy of the model with 93.7% and 5.9% of residues present in most favored and additional allowed regions of the plot respectively. The active sites, ion binding sites, and subcellular localization of laccase were also predicted. The generated model was docked with 121 pollutants (pesticides, insecticides, herbicides, fungicides, and rodenticides) for its degradation potential towards these pollutants. Two ligands chlorophacinone (based on the highest binding energy) and endosulfan (based on agricultural uses) were selected for molecular dynamic simulation studies. Endosulfan as a pesticide is banned but in some countries governments allow its use for special purposes which need serious consideration on developing bioremediation approaches for endosulfan degradation. MD simulation studies revealed that both chlorophacinone and endosulfan form hydrogen bonds and hydrophobic bonds with the active site of laccase and chlorophacinone-laccase complex were more stable in comparison to endosulfan. The present investigation provides insight into the structural features of laccase and its potential for the degradation of pesticides which can be further validated by experimental data.Communicated by Ramaswamy H. Sarma.

Hunting down zinc(II)-binding sites in proteins with distance matrices

MOTIVATION

In recent years, high-throughput sequencing technologies have made available the genome sequences of a huge variety of organisms. However, the functional annotation of the encoded proteins often still relies on low-throughput and costly experimental studies. Bioinformatics approaches offer a promising alternative to accelerate this process. In this work, we focus on the binding of zinc(II) ions, which is needed for 5%-10% of any organism's proteins to achieve their physiologically relevant form.

RESULTS

To implement a predictor of zinc(II)-binding sites in the 3D structures of proteins, we used a neural network, followed by a filter of the network output against the local structure of all known sites. The latter was implemented as a function comparing the distance matrices of the Cα and Cβ atoms of the sites. We called the resulting tool Master of Metals (MOM). The structural models for the entire proteome of an organism generated by AlphaFold can be used as input to our tool in order to achieve annotation at the whole organism level within a few hours. To demonstrate this, we applied MOM to the yeast proteome, obtaining a precision of about 76%, based on data for homologous proteins.

AVAILABILITY AND IMPLEMENTATION

Master of Metals has been implemented in Python and is available at https://github.com/cerm-cirmmp/Master-of-metals.

Evolutionarily conserved heat shock protein, HtpX, as an adjunct target against antibiotic-resistant Neisseria gonorrhoeae

The emergence of multiple drug resistance and extreme drug resistance pathogens necessitates the continuous evaluation of the pathogenic genome to identify conserved molecular targets and their respective inhibitors. In this study, we mapped the global mutational landscape of Neisseria gonorrhoeae (an intracellular pathogen notoriously known to cause the sexually transmitted disease gonorrhoea). We identified highly variable amino acid positions in the antibiotic target genes like the penA, ponA, 23s rRNA, rpoB, gyrA, parC, mtrR and porB. Some variations are directly reported to confer resistance to the currently used front-line drugs like ceftriaxone, cefixime, azithromycin and ciprofloxacin. Further, by whole genome comparison and Shannon entropy analysis, we identified a completely conserved protein HtpX in the drug-resistant as well as susceptible isolates of N. gonorrhoeae (NgHtpX). Comparison with the only available information of Escherichia coli HtpX suggested it to be a transmembrane metalloprotease having a role in stress response. The critical zinc-binding residue of NgHtpX was mapped to E141. By applying composite high throughput screening followed by MD simulations, we identified pemirolast and thalidomide as high-energy binding ligands of NgHtpX. Following cloning and expression of the purified metal-binding domain of NgHtpX (NgHtpXd), its Zn2+ -binding (Kd  = 0.4 µM) and drug-binding (pemirolast, Kd  = 3.47 µM; and thalidomide, Kd  = 1.04 µM) potentials were determined using in-vitro fluorescence quenching experiment. When tested on N. gonorrhoeae cultures, both the ligands imposed a dose-dependent reduction in viability. Overall, our results provide high entropy positions in the targets of presently used antibiotics, which can be further explored to understand the AMR mechanism. Additionally, HtpX and its specific inhibitors identified can be utilised effectively in managing gonococcal infections.

Phosphatase and Pseudo-Phosphatase Functions of Phosphatase of Regenerating Liver 3 (PRL-3) Are Insensitive to Divalent Metals In Vitro

Phosphatase of regenerating liver 3 (PRL-3) is associated with cancer metastasis and has been shown to interact with the cyclin and CBS domain divalent metal cation transport mediator (CNNM) family of proteins to regulate the intracellular concentration of magnesium and other divalent metals. Despite PRL-3's importance in cancer, factors that regulate PRL-3's phosphatase activity and its interactions with CNNM proteins remain unknown. Here, we show that divalent metal ions, including magnesium, calcium, and manganese, have no impact on PRL-3's structure, stability, phosphatase activity, or CNNM binding capacity, indicating that PRL-3 does not act as a metal sensor, despite its interaction with CNNM metal transporters. In vitro approaches found that PRL-3 is a broad but not indiscriminate phosphatase, with activity toward di- and tri-nucleotides, phosphoinositols, and NADPH but not other common metabolites. Although calcium, magnesium, manganese, and zinc-binding sites were predicted near the PRL-3 active site, these divalent metals did not specifically alter PRL-3's phosphatase activity toward a generic substrate, its transition from an inactive phospho-cysteine intermediate state, or its direct binding with the CBS domain of CNNM. PRL-3's insensitivity to metal cations negates the possibility of its role as an intracellular metal content sensor for regulating CNNM activity. Further investigation is warranted to define the regulatory mechanisms governing PRL-3's phosphatase activity and CNNM interactions, as these findings could hold potential therapeutic implications in cancer treatment.

In Silico Simulations Reveal Molecular Mechanism of Uranyl Ion Toxicity towards DNA-Binding Domain of PARP-1 Protein

The molecular toxicity of the uranyl ion (UO22+) in living cells is primarily determined by its high affinity to both native and potential metal-binding sites that commonly occur in the structure of biomolecules. Recent advances in computational and experimental research have shed light on the structural properties and functional impacts of uranyl binding to proteins, organic ligands, nucleic acids, and their complexes. In the present work, we report the results of the computational investigation of the uranyl-mediated loss of DNA-binding activity of PARP-1, a eukaryotic enzyme that participates in DNA repair, cell differentiation, and the induction of inflammation. The latest experimental studies have shown that the uranyl ion directly interacts with its DNA-binding subdomains, zinc fingers Zn1 and Zn2, and alters their tertiary structure. Here, we propose an atomistic mechanism underlying this process and compute the free energy change along the suggested pathway. Our Quantum Mechanics/Molecular Mechanics (QM/MM) simulations of the Zn2-UO22+ complex indicate that the uranyl ion replaces zinc in its native binding site. However, the resulting state is destroyed due to the spontaneous internal hydrolysis of the U-Cys162 coordination bond. Despite the enthalpy of hydrolysis being +2.8 kcal/mol, the overall reaction free energy change is -0.6 kcal/mol, which is attributed to the loss of domain's native tertiary structure originally maintained by a zinc ion. The subsequent reorganization of the binding site includes the association of the uranyl ion with the Glu190/Asp191 acidic cluster and significant perturbations in the domain's tertiary structure driven by a further decrease in the free energy by 6.8 kcal/mol. The disruption of the DNA-binding interface revealed in our study is consistent with previous experimental findings and explains the loss of PARP-like zinc fingers' affinity for nucleic acids.

Molecular docking in organic, inorganic, and hybrid systems: a tutorial review

Molecular docking simulation is a very popular and well-established computational approach and has been extensively used to understand molecular interactions between a natural organic molecule (ideally taken as a receptor) such as an enzyme, protein, DNA, RNA and a natural or synthetic organic/inorganic molecule (considered as a ligand). But the implementation of docking ideas to synthetic organic, inorganic, or hybrid systems is very limited with respect to their use as a receptor despite their huge popularity in different experimental systems. In this context, molecular docking can be an efficient computational tool for understanding the role of intermolecular interactions in hybrid systems that can help in designing materials on mesoscale for different applications. The current review focuses on the implementation of the docking method in organic, inorganic, and hybrid systems along with examples from different case studies. We describe different resources, including databases and tools required in the docking study and applications. The concept of docking techniques, types of docking models, and the role of different intermolecular interactions involved in the docking process to understand the binding mechanisms are explained. Finally, the challenges and limitations of dockings are also discussed in this review.

Graphical abstract

Zn2+ and Cu2+ Interaction with the Recognition Interface of ACE2 for SARS-CoV-2 Spike Protein

The spike protein (S) of SARS-CoV-2 is able to bind to the human angiotensin-converting enzyme 2 (ACE2) receptor with a much higher affinity compared to other coronaviruses. The binding interface between the ACE2 receptor and the spike protein plays a critical role in the entry mechanism of the SARS-CoV-2 virus. There are specific amino acids involved in the interaction between the S protein and the ACE2 receptor. This specificity is critical for the virus to establish a systemic infection and cause COVID-19 disease. In the ACE2 receptor, the largest number of amino acids playing a crucial role in the mechanism of interaction and recognition with the S protein is located in the C-terminal part, which represents the main binding region between ACE2 and S. This fragment is abundant in coordination residues such as aspartates, glutamates, and histidine that could be targeted by metal ions. Zn²⁺ ions bind to the ACE2 receptor in its catalytic site and modulate its activity, but it could also contribute to the structural stability of the entire protein. The ability of the human ACE2 receptor to coordinate metal ions, such as Zn²⁺, in the same region where it binds to the S protein could have a crucial impact on the mechanism of recognition and interaction of ACE2-S, with consequences on their binding affinity that deserve to be investigated. To test this possibility, this study aims to characterize the coordination ability of Zn²⁺, and also Cu²⁺ for comparison, with selected peptide models of the ACE2 binding interface using spectroscopic and potentiometric techniques.

Metal3D: a general deep learning framework for accurate metal ion location prediction in proteins

Metal ions are essential cofactors for many proteins and play a crucial role in many applications such as enzyme design or design of protein-protein interactions because they are biologically abundant, tether to the protein using strong interactions, and have favorable catalytic properties. Computational design of metalloproteins is however hampered by the complex electronic structure of many biologically relevant metals such as zinc . In this work, we develop two tools - Metal3D (based on 3D convolutional neural networks) and Metal1D (solely based on geometric criteria) to improve the location prediction of zinc ions in protein structures. Comparison with other currently available tools shows that Metal3D is the most accurate zinc ion location predictor to date with predictions within 0.70 ± 0.64 Å of experimental locations. Metal3D outputs a confidence metric for each predicted site and works on proteins with few homologes in the protein data bank. Metal3D predicts a global zinc density that can be used for annotation of computationally predicted structures and a per residue zinc density that can be used in protein design workflows. Currently trained on zinc, the framework of Metal3D is readily extensible to other metals by modifying the training data.

Characterization of a novel GH26 β-mannanase from Paenibacillus polymyxa and its application in the production of mannooligosaccharides

A novel glycoside hydrolase family 26 β-mannanase gene ppman26a was cloned from Paenibacillus polymyxa KF-1. The full-length enzyme PpMan26A and its truncated products CBM35pp (aa 35-328) and PpMan26A-Δ205 (aa 206-656) were overexpressed in Escherichia coli. PpMan26A hydrolyzed locust bean gum, guar gum, konjac gum and ivory nut mannan, with the highest specific activity toward konjac gum. The K_m and k_cat values for konjac gum were 2.13 mg/mL and 416.66 s^-1, respectively. The oligosaccharides fraction obtained from the hydrolysis of konjac gum by PpMan26A was analyzed by matrix-assisted laser desorption ionization-time-of-flight mass spectrometer (MALDI-TOF-MS). The degradation products were mainly mannooligosaccharides with a degree of polymerization of 3-8. CBM35pp exerted strong binding activity toward mannans but without β-mannanase activity. PpMan26A-Δ205, with the deletion of the N-terminal CBM domain, showed lower substrate binding capacity, resulting in reduced enzymatic activity and thermostability. This study complements our understanding of GH26 β-mannanases and expands the potential industrial application of PpMan26A.

Metalloproteome of human-infective RNA viruses: a study towards understanding the role of metal ions in virology

Metalloproteins and metal-based inhibitors have been shown to effectively combat infectious diseases, particularly those caused by RNA viruses. In this study, a diverse set of bioinformatics methods was employed to identify metal-binding proteins of human RNA viruses. Seventy-three viral proteins with a high probability of being metal-binding proteins were identified. These proteins included 40 zinc-, 47 magnesium- and 14 manganese-binding proteins belonging to 29 viral species and eight significant viral families, including Coronaviridae, Flaviviridae and Retroviridae. Further functional characterization has revealed that these proteins play a critical role in several viral processes, including viral replication, fusion and host viral entry. They fall under the essential categories of viral proteins, including polymerase and protease enzymes. Magnesium ion is abundantly predicted to interact with these viral enzymes, followed by zinc. In addition, this study also examined the evolutionary aspects of predicted viral metalloproteins, offering essential insights into the metal utilization patterns among different viral species. The analysis indicates that the metal utilization patterns are conserved within the functional classes of the proteins. In conclusion, the findings of this study provide significant knowledge on viral metalloproteins that can serve as a valuable foundation for future research in this area.