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.

The high phenotypic plasticity developed by plants includes rapid responses and adaptations to aggressive or changing environments. To achieve this, they evolved extremely efficient mechanisms of signaling mediated by a wide range of molecules, including small signal molecules. Among them, hydrogen cyanide (HCN) has been largely ignored due to its toxic characteristics. However, not only is it present in living organisms, but it has been shown that it serves several functions in all kingdoms of life. Research using model plants has changed the traditional point of view, and it has been demonstrated that HCN plays a positive role in the plant response to pathogens independently of its toxicity. Indeed, HCN induces a response aimed at protecting the plant from pathogen attack, and the HCN is provided either exogenously (in vitro or by some cyanogenic bacteria species present in the rhizosphere) or endogenously (in reactions involving ethylene, camalexin, or other cyanide-containing compounds). The contribution of different mechanisms to HCN function, including a new post-translational modification of cysteines in proteins, namely S-cyanylation, is discussed here. This work opens up an expanding \'HCN field\' of research related to plants and other organisms.

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.

In this study we propose to coin the term Wolfium bond (WfB) to refer to a net attractive force (noncovalent interaction) between any element of group 6 and electron donor atoms (neutral molecules or anions) and to differentiate it from a coordination bond (metal-ligand interaction). We provide evidence of the existence of this interaction by inspecting the X-ray crystal structure of proteins containing Molybdopterin and Tungstopterin cofactors from the Protein Data Bank (PDB). The plausible biological role of the interaction as well as its physical nature (antibonding Wf-Ligand orbital involved) are also analyzed by means of ab initio calculations (RI-MP2/def2-TZVP level of theory), Atoms in Molecules (AIM), Natural Bond Orbital (NBO) and Noncovalent Interactions plot (NCIplot) analyses.

The Mo/W-containing metalloenzyme formate dehydrogenase (FDH) is an efficient and selective natural catalyst that reversibly converts CO2 to formate under ambient conditions. In this study, we investigate the impact of the greater protein environment on the electrostatic potential (ESP) of the active site. To model the enzyme environment, we used a combination of classical molecular dynamics and multiscale quantum-mechanical (QM)/molecular-mechanical (MM) simulations. We leverage charge shift analysis to systematically construct QM regions and analyze the electronic environment of the active site by evaluating the degree of charge transfer between the core active site and the protein environment. The contribution of the terminal chalcogen ligand to the ESP of the metal center is substantial and dependent on the chalcogen identity, with similar, less negative ESPs for Se and S terminal chalcogens in comparison to O regardless of whether the metal is Mo or W. The orientation of the side chains and conformations of the cofactor also affect the ESP, highlighting the importance of sampling dynamic fluctuations in the protein. Overall, our observations suggest that the terminal chalcogen ligand identity plays an important role in the enzymatic activity of FDH, suggesting opportunities for a rational bioinspired catalyst design.

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.

Mycobacterium tuberculosis (Mtb) has an old history as a human pathogen and still kills over one million people every year. One key feature of this bacterium is its dormancy: a phenomenon responsible for major changes in its metabolism and replication that have been associated with the need for a lengthy therapy for Mtb. This process is regulated by key heme-based sensors, particularly DosT and DevS (DosS), among other co-regulators, and also linked to nitrogen utilization (nitrate/nitrite) and stringent responses. In face of the current threat of tuberculosis, there is an urgent need to develop new therapeutic agents capable of targeting the dormant state, associated with the need for a lengthy therapy. Interestingly, many of those key proteins are indeed metallo-containing or metallo-dependent biomolecules, opening exciting bioinorganic opportunities. Here, we critically reviewed a series of small molecules targeting key proteins involved in these processes, including DosT/DevS/DevR, RegX3, MprA, MtrA, NarL, PknB, Rel, PPK, nitrate and nitrite reductases, GlnA1, aiming for new opportunities and alternative therapies. In the battle against Mycobacterium tuberculosis, new drug targets must be searched, in particular  those involved in dormancy. A series of exciting cases for drug development involving metallo-containing or metallo-dependent biomolecules are reviewed, opening great opportunities for the bioinorganic chemistry community.

We experimentally studied the pressure dependence of the zero-field splitting (ZFS) parameter of hemin (iron(III) protoporphyrin IX chloride), which is a model complex of hemoproteins, via high-frequency and high-field electron paramagnetic resonance (HFEPR) under pressure. Owing to the large ZFS, the pressure effect on the electronic structure of iron-porphyrin complexes has not yet been explored using EPR. Therefore, we systematically studied this effect using our newly developed sub-terahertz EPR spectroscopy system in the frequency range of 80-515 GHz, under magnetic fields up to 10 T and pressure up to 2 GPa. We observed a systematic shift of the resonance fields of hemin upon pressure application, from which the axial component of the ZFS parameter was found to increase from D = 6.9 to 7.9 cm-1 at 2 GPa. In contrast to the previous methods used to study proteins under pressure, which mainly focused on conformational changes, our HFEPR technique can obtain more microscopic insights into the electronic structures of metal ions under pressure. In this sense, our technique provides novel opportunities to study the pressure effects on biofunctional active centers of versatile metalloproteins.

Molybdenum cofactor deficiency is an autosomal recessive inborn error of metabolism, which results from mutations in genes involved in Moco biosynthesis. Moco serves as a cofactor of several enzymes, including sulfite oxidase. MoCD is clinically characterized by intractable seizures and severe, rapidly progressing neurodegeneration leading to death in early childhood in the majority of known cases. Here we report a patient with an unusual late disease onset and mild phenotype, characterized by a lack of seizures, normal early development, a decline triggered by febrile illness and a subsequent dystonic movement disorder. Genetic analysis revealed a homozygous c.1338delG MOCS1 mutation causing a frameshift (p.S442fs) with a premature termination of the MOCS1AB translation product at position 477 lacking the entire MOCS1B domain. Surprisingly, urine analysis detected trace amounts (1% of control) of the Moco degradation product urothione, suggesting a residual Moco synthesis in the patient, which was consistent with the mild clinical presentation. Therefore, we performed bioinformatic analysis of the patient\'s mutated MOCS1 transcript and found a potential Kozak-sequence downstream of the mutation site providing the possibility of an independent expression of a MOCS1B protein. Following the expression of the patient\'s MOCS1 cDNA in HEK293 cells we detected two proteins: a truncated MOCS1AB protein and a 22.4 kDa protein representing MOCS1B. Functional studies of both proteins confirmed activity of MOCS1B, but not of the truncated MOCS1AB. This finding demonstrates an unusual mechanism of translation re-initiation in the MOCS1 transcript, which results in trace amounts of functional MOCS1B protein being sufficient to partially protect the patient from the most severe symptoms of MoCD.

Zn-metalloproteins are a major class of targets for drug design. They constitute a demanding testing ground for polarizable molecular mechanics/dynamics aimed at extending the realm of quantum chemistry (QC) to very long-duration molecular dynamics (MD). The reliability of such procedures needs to be demonstrated upon comparing the relative stabilities of competing candidate complexes of inhibitors with the recognition site stabilized in the course of MD. This could be necessary when no information is available regarding the experimental structure of the inhibitor-protein complex. Thus, this study bears on the phosphomannose isomerase (PMI) enzyme, considered as a potential therapeutic target for the treatment of several bacterial and parasitic diseases. We consider its complexes with 5-phospho-d-arabinonohydroxamate and three analog ligands differing by the number and location of their hydroxyl groups. We evaluate the energy accuracy expectable from a polarizable molecular mechanics procedure, SIBFA. This is done by comparisons with ab initio quantum-chemistry (QC) calculations in the following cases: (a) the complexes of the four ligands in three distinct structures extracted from the entire PMI-ligand energy-minimized structures, and totaling up to 264 atoms; (b) the solvation energies of several energy-minimized complexes of each ligand with a shell of 64 water molecules; (c) the conformational energy differences of each ligand in different conformations characterized in the course of energy-minimizations; and (d) the continuum solvation energies of the ligands in different conformations. The agreements with the QC results appear convincing. On these bases, we discuss the prospects of applying the procedure to ligand-macromolecule recognition problems. © 2016 Wiley Periodicals, Inc.