Polarization and charge-transfer interactions play an important role in ligand-receptor complexes containing metals, and only quantum mechanics methods can adequately describe their contribution to the binding energy. In this work, we selected a set of benzenesulfonamide ligands of human Carbonic Anhydrase II (hCA II)-an important druggable target containing a Zn2+ ion in the active site-as a case study to predict the binding free energy in metalloprotein-ligand complexes and designed specialized computational methods that combine the ab initio fragment molecular orbital (FMO) method and GRID approach. To reproduce the experimental binding free energy in these systems, we adopted a machine-learning approach, here named formula generator (FG), considering different FMO energy terms, the hydrophobic interaction energy (computed by GRID) and logP. The main advantage of the FG approach is that it can find nonlinear relations between the energy terms used to predict the binding free energy, explicitly showing their mathematical relation. This work showed the effectiveness of the FG approach, and therefore, it might represent an important tool for the development of new scoring functions. Indeed, our scoring function showed a high correlation with the experimental binding free energy (R2 = 0.76-0.95, RMSE = 0.34-0.18), revealing a nonlinear relation between energy terms and highlighting the relevant role played by hydrophobic contacts. These results, along with the FMO characterization of ligand-receptor interactions, represent important information to support the design of new and potent hCA II inhibitors.

Since fish metalloproteins are still not thoroughly characterized, the aim of this study was to investigate the acidic/basic nature of biomolecules involved in the sequestration of twelve selected metals in the soluble hepatic fraction of an important aquatic bioindicator organism, namely the fish species northern pike (Esox lucius). For this purpose, the hyphenated system HPLC-ICP-MS was applied, with chromatographic separation based on anion/cation-exchange principle at physiological pH (7.4). The results indicated predominant acidic nature of metal-binding peptides/proteins in the studied hepatic fraction. More than 90 % of Ag, Cd, Co, Cu, Fe, Mo, and Pb were eluted with negatively charged biomolecules, and >70 % of Bi, Mn, and Zn. Thallium was revealed to bind equally to negatively and positively charged biomolecules, and Cs predominantly to positively charged ones. The majority of acidic (negatively charged) metalloproteins/peptides were coeluted within the elution time range of applied standard proteins, having pIs clustered around 4-6. Furthermore, binding of several metals (Ag, Cd, Cu, Zn) to two MT-isoforms was assumed, with Cd and Zn preferentially bound to MT1 and Ag to MT2, and Cu evenly distributed between the two. The results presented here are the first of their kind for the important bioindicator species, the northern pike, as well as one of the rare comprehensive studies on the acidic/basic nature of metal-binding biomolecules in fish, which can contribute significantly to a better understanding of the behaviour and fate of metals in the fish organism, specifically in liver as main metabolic and detoxification organ.

The diiron active site is pivotal in catalyzing transformations in both biological and chemical systems. Recently, a range of biomimetic diiron catalysts have been synthesized, drawing inspiration from the active architecture of soluble methane monooxygenase (sMMO). These catalysts have been successfully deployed for the dehydrogenation of indolines, marking a significant advancement in the field. Using density functional theory (DFT) calculations, we have identified a novel mechanistic pathway that governs the dehydrogenation of indolines catalyzed by a biomimetic diiron complex. Specifically, this reaction is facilitated by the transfer of a hybrid atom from the C1 position of the substrate to the distal oxygen atom of the Fe(III)Fe(III)-1,1-μ-hydroperoxy active species. This transfer serves as the rate-limiting step for the heterolytic cleavage of the OO bond, ultimately generating the substrate cation. The mechanism we propose aligns well with mechanistic investigations incorporating both kinetic isotope effect (KIE) measurements and evaluations of stereochemical selectivity. This research contributes to the broader scientific understanding of catalysis involving biomimetic diiron complexes and offers valuable insights into the catalytic behaviors of non-heme diiron metalloenzymes.

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.

Xanthine oxidoreductase (XOR) is an enzyme found in various organisms. It converts hypoxanthine to xanthine and urate, which are crucial steps in purine elimination in humans. Elevated uric acid levels can lead to conditions like gout and hyperuricemia. Therefore, there is significant interest in developing drugs that target XOR for treating these conditions and other diseases. Oxipurinol, an analogue of xanthine, is a well-known inhibitor of XOR. Crystallographic studies have revealed that oxipurinol directly binds to the molybdenum cofactor (MoCo) in XOR. However, the precise details of the inhibition mechanism are still unclear, which would be valuable for designing more effective drugs with similar inhibitory functions. In this study, molecular dynamics and quantum mechanics/molecular mechanics calculations are employed to investigate the inhibition mechanism of XOR by oxipurinol. The study examines the structural and dynamic effects of oxipurinol on the pre-catalytic structure of the metabolite-bound system. Our results provide insights on the reaction mechanism catalyzed by the MoCo center in the active site, which aligns well with experimental findings. Furthermore, the results provide insights into the residues surrounding the active site and propose an alternative mechanism for developing alternative covalent inhibitors.

The COVID-19 pandemic has generated a major interest in designing inhibitors to prevent SARS-CoV-2 binding on host cells to protect against infection. One promising approach to such research utilizes molecular dynamics simulation to identify potential inhibitors that can prevent the interaction between spike (S) protein on the virus and angiotensin converting enzyme 2 (ACE2) receptor on the host cells. In these studies, many groups have chosen to exclude the ACE2-bound zinc (Zn) ion, which is critical for its enzymatic activity. While the relatively distant location of Zn ion from the S protein binding site (S1 domain), combined with the difficulties in modeling this ion has motivated the decision of exclusion, Zn can potentially contribute to the structural stability of the entire protein, and thus, may have implications on S protein-ACE2 interaction. In this study, the authors model both the ACE2-S1 and ACE2-inhibitor (mAb) system to investigate if there are variations in structure and the readouts due to the presence of Zn ion. Although distant from the S1 or inhibitor binding region, inclusion/exclusion of Zn has statistically significant effects on the structural stability and binding free energy in these systems. In particular, the binding free energy of the ACE2-S1 and ACE2-inhibitor structures is - 3.26 and - 14.8 kcal/mol stronger, respectively, in the Zn-bound structure than in the Zn-free structures. This finding suggests that including Zn may be important in screening potentially inhibitors and may be particularly important in modeling monoclonal antibodies, which may be more sensitive to changes in antigen structure.

We present a novel decomposition scheme for electronic interaction energies based on the flexible formulation of fragmentation schemes through fragment combination ranges (FCRs; J. Chem. Phys., 2021, 155, 164105). We devise a clear additive decomposition with contribution of nondisjoint fragments and correction terms for overlapping fragments and apply this scheme to the metalloenzyme-substrate complex of a lytic polysaccharide monooxygenase (LPMO) with an oligosaccharide. By this, we further illustrate the straightforward adaptability of the FCR-based schemes to novel systems. Our calculations suggest that the description of the electronic structure is a larger error source than the fragmentation scheme. In particular, we find a large impact of the basis set size on the interaction energies. Still, the introduction of three-body interaction terms in the fragmentation setup improves the agreement to the supermolecular reference. Yet, the qualitative results for the decomposition scheme with two-body terms only largely agree within the investigated electronic-structure approaches and basis sets, which are B97-3c, DFT (TPSS and B3LYP), and MP2 methods. The overlap contributions are found to be small, allowing analysis of the interaction energy into individual amino acid residues: We find a particularly strong interaction between the substrate and the LPMO copper active site.

OBJECTIVE: To investigate the potential correlation of transforming growth factor-β (TGF-β), matrix metalloprotein 9 (MMP-9), tissue inhibitor of metalloproteinases 1 (TIMP-1), Interleukin 1 (IL-1), IL-4, IL-6, IL-17, and tumor necrosis factor alpha (TNF-α) in refractory chronic rhinosinusitis.
METHODS: A total of 150 participants were retrospectively included in this study from August 2018 to February 2020. The people enrolled were equally allocated into refractory group (patients with refractory chronic rhinosinusitis), chronic group (patients with chronic rhinosinusitis), and control group (normal people). The level of TGF-β1, MMP-9, TIMP-1, IL-1, IL-4, IL-6, IL-17, and TNF-α were recorded. The unconditional multivariate binary logistic regression was used to analyze the factors affecting refractory chronic rhinosinusitis.
RESULTS: The Davos score, T&T olfactometer threshold test, and Lund-Mackay CT scores in refractory group were significantly higher than the chronic group (P<0.05). The level of TGF-β1, MMP-9, TIMP-1, IL-1, IL-4, IL-6, IL-17, and TNF-α in the refractory group were significantly higher than the chronic group and the control group (all P<0.05). Similarly, the level of the above mentioned indexes in the chronic group were significantly higher than the control group (P<0.05). The Davos score, T&T olfactometer threshold test score, Lund-Mackay CT score, and the level of TGF-β1, MMP-9, TIMP-1, IL-1, IL-4, IL-6, IL-17, and TNF-α positively correlated with refractory chronic rhinosinusitis. Moreover, the unconditional multivariate binary logistic regression showed that the influencing factors of refractory chronic rhinosinusitis included TGF-β1, MMP-9, TIMP-1, IL-1, IL-4, IL-6, IL-17, and TNF-α.
CONCLUSIONS: The findings of the present study provide evidence for TGF-β1, MMP-9, TIMP-1, IL-4, IL-6, IL-17, and TNF-α as the influencing factors of refractory chronic rhinosinusitis.

Molecular modelling applications in metalloenzyme design are still scarce due to a series of challenges. On top of that, the simulations of metal-mediated binding and the identification of catalytic competent geometries require both large conformational exploration and simulation of fine electronic properties. Here, we demonstrate how the incorporation of new tools in multiscale strategies, namely substrate diffusion exploration, allows taking a step further. As a showcase, the enantioselective profiles of the most outstanding variants of an artificial Rh2-based cyclopropanase (GSH, HFF and RFY) developed by Lewis and co-workers (Nat. Commun., 2015, 6, 7789 and Nat. Chem., 2018, 10, 318-324) have been rationalized. DFT calculations on the free-cofactor-mediated process identify the carbene insertion and the cyclopropanoid formation as crucial events, the latter being the enantiodetermining step, which displays up to 8 competitive orientations easily altered by the protein environment. The key intermediates of the reaction were docked into the protein scaffold showing that some mutated residues have direct interaction with the cofactor and/or the co-substrate. These interactions take the form of a direct coordination of Rh in GSH and HFF and a strong hydrophobic patch with the carbene moiety in RFY. Posterior molecular dynamics sustain that the cofactor induces global re-arrangements of the protein. Finally, massive exploration of substrate diffusion, based on the GPathFinder approach, defines this event as the origin of the enantioselectivity in GSH and RFY. For HFF, fine molecular dockings suggest that it is likely related to local interactions upon diffusion. This work shows how modelling of long-range mutations on the catalytic profiles of metalloenzymes may be unavoidable and software simulating substrate diffusion should be applied.

Oral iron supplementation constitutes the first line treatment for iron deficiency anemia (IDA), with daily doses between 80 mg and 200 mg of elemental iron. Ferrous salts, such as ferrous sulphate (FeSO4), while efficacious, frequently give rise to gastrointestinal side effects. In the present paper we attempted to directly compare the efficacy of an alternative to the FeSO4 formulation, which presents a better tolerability profile, iron protein succinylate (Ferplex®). In a diet-induced anemia model, rats were treated by oral gavage with vehicle, FeSO4, or Ferplex® at a human-dose equivalent of 80 mg and 200 mg of elemental iron. We evaluated the change in anemia-related hematological and biochemical parameters, conducting a histological examination of the intestine at sacrifice. Results indicate that both types of iron supplementation are equally effective in the treatment of IDA, restoring hemoglobin, hematocrit, erythrocytes, free iron and transferrin levels in 15 days, with no statistical differences between treated groups and control. The impact of anemia on body weight was also attenuated following treatment with both iron supplements. Thrombocyte and reticulocyte levels, altered by the anemic condition, returned to homeostasis after 15 days of either FeSO4 or Ferplex® treatment. Importantly, the lower and higher doses of iron were equally effective, thus supporting the current school of thought which states that lower therapeutic doses are sufficient for management of IDA. In addition, the study shows for the first time that oral treatment with Ferplex® does not increase serum hepcidin. Finally, Ferplex® induced minimal iron depositions in the intestinal tissue compared to FeSO4.