Hydration numbers of biologically relevant divalent metal cations from ab initio molecular dynamics and continuum solvation methods.

Hydration and, in particular, the coordination number of a metal ion is of paramount importance as it defines many of its (bio)physicochemical properties. It is not only essential for understanding its behavior in aqueous solutions but also determines the metal ion reference state and its binding energy to (bio)molecules. In this paper, for divalent metal cations Ca2+, Cd2+, Cu2+, Fe2+, Hg2+, Mg2+, Ni2+, Pb2+, and Zn2+, we compare two approaches for predicting hydration numbers: (1) a mixed explicit/continuum DFT-D3//COSMO-RS solvation model and (2) density functional theory based ab initio molecular dynamics. The former approach is employed to calculate the Gibbs free energy change for the sequential hydration reactions, starting from [M(H2O)2]2+ aqua complexes to [M(H2O)9]2+, allowing explicit water molecules to bind in the first or second coordination sphere and determining the most stable [M(H2O)n]2+ structure. In the latter approach, the hydration number is obtained by integrating the ion-water radial distribution function. With a couple of exceptions, the metal ion hydration numbers predicted by the two approaches are in mutual agreement, as well as in agreement with the experimental data.

What are the minimal folding seeds in proteins? Experimental and theoretical assessment of secondary structure propensities of small peptide fragments.

Certain peptide sequences, some of them as short as amino acid triplets, are significantly overpopulated in specific secondary structure motifs in folded protein structures. For example, 74% of the EAM triplet is found in α-helices, and only 3% occurs in the extended parts of proteins (typically ß-sheets). In contrast, other triplets (such as VIV and IYI) appear almost exclusively in extended parts (79% and 69%, respectively). In order to determine whether such preferences are structurally encoded in a particular peptide fragment or appear only at the level of a complex protein structure, NMR, VCD, and ECD experiments were carried out on selected tripeptides: EAM (denoted as pro-'α-helical' in proteins), KAM(α), ALA(α), DIC(α), EKF(α), IYI(pro-ß-sheet or more generally, pro-extended), and VIV(ß), and the reference α-helical CATWEAMEKCK undecapeptide. The experimental data were in very good agreement with extensive quantum mechanical conformational sampling. Altogether, we clearly showed that the pro-helical vs. pro-extended propensities start to emerge already at the level of tripeptides and can be fully developed at longer sequences. We postulate that certain short peptide sequences can be considered minimal "folding seeds". Admittedly, the inherent secondary structure propensity can be overruled by the large intramolecular interaction energies within the folded and compact protein structures. Still, the correlation of experimental and computational data presented herein suggests that the secondary structure propensity should be considered as one of the key factors that may lead to understanding the underlying physico-chemical principles of protein structure and folding from the first principles.

Acidity constants and protonation sites of cyclic dinucleotides determined by capillary electrophoresis, quantum chemical calculations, and NMR spectroscopy.

Cyclic dinucleotides (CDNs) are important second messengers in bacteria and eukaryotes. Detailed characterization of their physicochemical properties is a prerequisite for understanding their biological functions. Herein, we examine acid-base and electromigration properties of selected CDNs employing capillary electrophoresis (CE), density functional theory (DFT), and nuclear magnetic resonance (NMR) spectroscopy to provide benchmark pKa values, as well as to unambiguously determine the protonation sites. Acidity constants (pKa ) of the NH+ moieties of adenine and guanine bases and actual and limiting ionic mobilities of CDNs were determined by nonlinear regression analysis of the pH dependence of their effective electrophoretic mobilities measured by CE in aqueous background electrolytes in a wide pH range (0.98-11.48), at constant temperature (25°C), and constant ionic strength (25 mM). The thermodynamic pKa values were found to be in the range 3.31-4.56 for adenine and 2.28-3.61 for guanine bases, whereas the pKa of enol group of guanine base was in the range 10.21-10.40. Except for systematic shifts of ∼2 pKa , the pKa values calculated by the DFT-D3//COSMO-RS composite protocol that included large-scale conformational sampling and "cross-morphing" were in a relatively good agreement with the pKa s determined by CE and predict N1 atom of adenine and N7 atom of guanine as the protonation sites. The protonation of the N1 atom of adenine and N7 atom of guanine in acidic background electrolytes (BGEs) and the dissociation of the enol group of guanine in alkaline BGEs was confirmed also by NMR spectroscopy.

Experimental Evidence and Mechanistic Description of the Phenolic H-Transfer to the Cu2O2 Active Site of oxy-Tyrosinase.

Tyrosinase is a ubiquitous coupled binuclear copper enzyme that activates O2 toward the regioselective monooxygenation of monophenols to catechols via a mechanism that remains only partially defined. Here, we present new mechanistic insights into the initial steps of this monooxygenation reaction by employing a pre-steady-state, stopped-flow kinetics approach that allows for the direct measurement of the monooxygenation rates for a series of para-substituted monophenols by oxy-tyrosinase. The obtained biphasic Hammett plot and the associated solvent kinetic isotope effect values provide direct evidence for an initial H-transfer from the protonated phenolic substrate to the Cu2O2 core of oxy-tyrosinase. The correlation of these experimental results to quantum mechanics/molecular mechanics calculations provides a detailed mechanistic description of this H-transfer step. These new mechanistic insights revise and expand our fundamental understanding of Cu2O2 active sites in biology.

Discovery of a Druggable, Cryptic Pocket in SARS-CoV-2 nsp16 Using Allosteric Inhibitors.

A collaborative, open-science team undertook discovery of novel small molecule inhibitors of the SARS-CoV-2 nsp16-nsp10 2'-O-methyltransferase using a high throughput screening approach with the potential to reveal new inhibition strategies. This screen yielded compound 5a, a ligand possessing an electron-deficient double bond, as an inhibitor of SARS-CoV-2 nsp16 activity. Surprisingly, X-ray crystal structures revealed that 5a covalently binds within a previously unrecognized cryptic pocket near the S-adenosylmethionine binding cleft in a manner that prevents occupation by S-adenosylmethionine. Using a multidisciplinary approach, we examined the mechanism of binding of compound 5a to the nsp16 cryptic pocket and developed 5a derivatives that inhibited nsp16 activity and murine hepatitis virus replication in rat lung epithelial cells but proved cytotoxic to cell lines canonically used to examine SARS-CoV-2 infection. Our study reveals the druggability of this newly discovered SARS-CoV-2 nsp16 cryptic pocket, provides novel tool compounds to explore the site, and suggests a new approach for discovery of nsp16 inhibition-based pan-coronavirus therapeutics through structure-guided drug design.

Comprehensive Mechanistic View of the Hydrolysis of Oxadiazole-Based Inhibitors by Histone Deacetylase 6 (HDAC6).

Histone deacetylase (HDAC) inhibitors used in the clinic typically contain a hydroxamate zinc-binding group (ZBG). However, more recent work has shown that the use of alternative ZBGs, and, in particular, the heterocyclic oxadiazoles, can confer higher isoenzyme selectivity and more favorable ADMET profiles. Herein, we report on the synthesis and biochemical, crystallographic, and computational characterization of a series of oxadiazole-based inhibitors selectively targeting the HDAC6 isoform. Surprisingly, but in line with a very recent finding reported in the literature, a crystal structure of the HDAC6/inhibitor complex revealed that hydrolysis of the oxadiazole ring transforms the parent oxadiazole into an acylhydrazide through a sequence of two hydrolytic steps. An identical cleavage pattern was also observed both in vitro using the purified HDAC6 enzyme as well as in cellular systems. By employing advanced quantum and molecular mechanics (QM/MM) and QM calculations, we elucidated the mechanistic details of the two hydrolytic steps to obtain a comprehensive mechanistic view of the double hydrolysis of the oxadiazole ring. This was achieved by fully characterizing the reaction coordinate, including identification of the structures of all intermediates and transition states, together with calculations of their respective activation (free) energies. In addition, we ruled out several (intuitively) competing pathways. The computed data (ΔG‡ ≈ 21 kcal·mol-1 for the rate-determining step of the overall dual hydrolysis) are in very good agreement with the experimentally determined rate constants, which a posteriori supports the proposed reaction mechanism. We also clearly (and quantitatively) explain the role of the -CF3 or -CHF2 substituent on the oxadiazole ring, which is a prerequisite for hydrolysis to occur. Overall, our data provide compelling evidence that the oxadiazole warheads can be efficiently transformed within the active sites of target metallohydrolases to afford reaction products possessing distinct selectivity and inhibition profiles.

Can Copper(I) and Silver(I) be Hydrogen Bond Acceptors?

Gold(I) centers can form moderately strong (Au⋅⋅⋅H) hydrogen bonds with tertiary ammonium groups, as has been demonstrated in the 3AuCl+ (3+ =1-(tert-butyl)-3-phenyl-4-(2-((dimethylammonio)methyl)phenyl)-1,2,4-triazol-5-ylidene) complex. However, similar hydrogen bonding interactions with isoelectronic silver(I) or copper(I) centers are unknown. Herein, we first explored whether the Au⋅⋅⋅H bond originally observed in 3AuCl+ can be strengthened by replacing Cl with Br or I. Experimental gas-phase IR spectra in the ∼3000 cm-1 region showed only a small effect of the halogen on the Au⋅⋅⋅H bond. Next, we measured the spectra of 3AgCl+ , which exhibited significant differences compared to its 3AuX+ counterparts. The difference has been explained by DFT calculations which indicated that the Ag⋅⋅⋅H interaction is only marginal in this complex, and a Cl⋅⋅⋅H hydrogen bond is formed instead. Calculations predicted the same for the 3CuCl+ complex. However, we noticed that for Ag and Cu complexes with less flexible ligands, such as dimethyl(2-(dimethylammonio)phenyl)phosphine (L7 H+ ), the computations predict the presence of the respective Ag⋅⋅⋅H and Cu⋅⋅⋅H hydrogen bonds, with a strength similar to the Au⋅⋅⋅H bond in 3AuCl+ . We, therefore, propose possible complexes where the presence of (Ag/Cu)⋅⋅⋅H bonds could be experimentally verified to broaden our understanding of these unusual interactions.

Redox Potentials with COSMO-RS: Systematic Benchmarking with Different Databases.

Recent techniques of computational electrochemistry can yield redox potentials with accuracy as good as 0.1 V. Yet, many such methods are not universal, easy to use, or computationally efficient. Herein, we provide a systematic benchmarking of a relatively cheap and straightforward computational approach for fairly accurate computations of redox potentials. It is based on a combination of the conductor-like screening model for real solvents (COSMO-RS) and the density functional theory (DFT). The benchmarking is done with databases covering diverse redox systems, including transition-metal aquacomplexes and various organic and inorganic compounds. In addition, we also present our own test set aiming at maximum chemical diversity and maximum range of redox potential values. We assess the performance of the fairly efficient computational protocol combining the COSMO-RS with the BP86 DFT functional. This is done by calibrating it against different high-level state-of-the-art techniques, in particular, polarizable continuum model (PCM) connected to composite G3(MP2,CC)(+) method, domain-based pair natural orbital implementation of coupled cluster theory, or complete basis set CBS-QB3 method. We also elaborate on the absolute reduction potential value of standard hydrogen electrode to be used with COSMO-RS, and we propose the value of approx. 4.4 V. The COSMO-RS/BP86-D3(BJ) combination outperforms other methods on a wide range of redox systems. However, we show that its accuracy is based on a balanced error cancelation and, therefore, it cannot be further systematically improved. As a result, the proposed procedure represents a pragmatic choice for large-scale screening, yet it could be combined with more advanced methods for detailed studies.

Quantum Bio-Inorganic Chemistry (QBIC) Society Special Collection.

This Editorial by Guest Editors L. Rulísek, M. Gruden, M. Orio and QBIC Society President R. J. Deeth summarizes the Chemistry Europe Special Collection on quantum bioinorganic chemistry, published in collaboration with the Quantum Bioinorganic Chemistry (QBIC) Society. The Collection, which recognizes of the strides made in this area of research, spans several of our portfolio journals Chemistry-A European Journal, the European Journal of Inorganic Chemistry, ChemPlusChem, ChemBioChem, ChemPhysChem, and ChemistryOpen, includes outstanding work in the field by authors across the world and we hope you enjoy reading them!

Elucidation of the tyrosinase/O2/monophenol ternary intermediate that dictates the monooxygenation mechanism in melanin biosynthesis.

Melanins are highly conjugated biopolymer pigments that provide photoprotection in a wide array of organisms, from bacteria to humans. The rate-limiting step in melanin biosynthesis, which is the ortho-hydroxylation of the amino acid L-tyrosine to L-DOPA, is catalyzed by the ubiquitous enzyme tyrosinase (Ty). Ty contains a coupled binuclear copper active site that binds O2 to form a µ:η2:η2-peroxide dicopper(II) intermediate (oxy-Ty), capable of performing the regioselective monooxygenation of para-substituted monophenols to catechols. The mechanism of this critical monooxygenation reaction remains poorly understood despite extensive efforts. In this study, we have employed a combination of spectroscopic, kinetic, and computational methods to trap and characterize the elusive catalytic ternary intermediate (Ty/O2/monophenol) under single-turnover conditions and obtain molecular-level mechanistic insights into its monooxygenation reactivity. Our experimental results, coupled with quantum-mechanics/molecular-mechanics calculations, reveal that the monophenol substrate docks in the active-site pocket of oxy-Ty fully protonated, without coordination to a copper or cleavage of the µ:η2:η2-peroxide O-O bond. Formation of this ternary intermediate involves the displacement of active-site water molecules by the substrate and replacement of their H bonds to the µ:η2:η2-peroxide by a single H bond from the substrate hydroxyl group. This H-bonding interaction in the ternary intermediate enables the unprecedented monooxygenation mechanism, where the µ-η2:η2-peroxide O-O bond is cleaved to accept the phenolic proton, followed by substrate phenolate coordination to a copper site concomitant with its aromatic ortho-hydroxylation by the nonprotonated µ-oxo. This study provides insights into O2 activation and reactivity by coupled binuclear copper active sites with fundamental implications in biocatalysis.

Cationic Gold(II) Complexes: Experimental and Theoretical Study.

Gold(II) complexes are rare, and their application to the catalysis of chemical transformations is underexplored. The reason is their easy oxidation or reduction to more stable gold(III) or gold(I) complexes, respectively. We explored the thermodynamics of the formation of [AuII (L)(X)]+ complexes (L=ligand, X=halogen) from the corresponding gold(III) precursors and investigated their stability and spectral properties in the IR and visible range in the gas phase. The results show that the best ancillary ligands L for stabilizing gaseous [AuII (L)(X)]+ complexes are bidentate and tridentate ligands with nitrogen donor atoms. The electronic structure and spectral properties of the investigated gold(II) complexes were correlated with quantum chemical calculations. The results show that the molecular and electronic structure of the gold(II) complexes as well as their spectroscopic properties are very similar to those of analogous stable copper(II) complexes.

Exhaustive Mapping of the Conformational Space of Natural Dipeptides by the DFT-D3//COSMO-RS Method.

We extensively mapped energy landscapes and conformations of 22 (including three His protonation states) proteinogenic α-amino acids in trans configuration and the corresponding 484 (222) dipeptides. To mimic the environment in a protein chain, the N- and C-termini of the studied systems were capped with acetyl and N-methylamide groups, respectively. We systematically varied the main chain dihedral angles (ϕ, ψ) by 40° steps and all side chain angles by 90° or 120° steps. We optimized the molecular geometries with the GFN2-xTB semiempirical (SQM) method and performed single point density functional theory calculations at the BP86-D3/DGauss-DZVP//COSMO-RS level in water, 1-octanol, N,N-dimethylformamide, and n-hexane. For each restrained (nonequilibrium) structure, we also calculated energy gradients (in water) and natural atomic charges. The exhaustive and unprecedented QM-based sampling enabled us to construct Ramachandran plots of quantum mechanical (QM(BP86-D3)//COSMO-RS) energies calculated on SQM structures, for all 506 (484 dipeptides and 22 amino acids) studied systems. We showed how the character of an amino acid side chain influences the conformational space of single amino acids and dipeptides. With clustering techniques, we were able to identify unique minima of amino acids and dipeptides (i.e., minima on the GFN2-xTB potential energy surfaces) and analyze the distribution of their BP86-D3//COSMO-RS conformational energies in all four solvents. We also derived an empirical formula for the number of unique minima based on the overall number of rotatable bonds within each peptide. The final peptide conformer data set (PeptideCs) comprises over 400 million structures, all of them annotated with QM(BP86-D3)//COSMO-RS energies. Thanks to its completeness and unbiased nature, the PeptideCs can serve, inter alia, as a data set for the validation of new methods for predicting the energy landscapes of protein structures. This data set may also prove to be useful in the development and reparameterization of biomolecular force fields. The data set is deposited at Figshare (10.25452/figshare.plus.19607172) and can be accessed using a simple web interface at http://peptidecs.uochb.cas.cz.

Cyclic Octapeptides Composed of Two Glutathione Units Outperform the Monomer in Lead Detoxification.

A rationally-designed scaffold of cyclic octapeptides composed of two units of the natural tripeptide glutathione (GSH) was optimized to strongly and selectively capture toxic lead ions (Pb(II)). Using state-of-the-art computational tools, a list of eleven plausible peptides was shortened to five analogs based on their calculated affinity to Pb(II) ions. We then synthesized and investigated them for their abilities to recover Pb-poisoned human cells. A clear pattern was observed from the in vitro detoxification results, indicating the importance of cavity size and polar moieties to enhance metal capturing. These, together with the apparent benefit of cyclizing the peptides, improved the detoxification of the two lead peptides by approximately two folds compared to GSH and the benchmark chelating agents against Pb poisoning. Moreover, the two peptides did not show any toxicity and, therefore, were thoroughly investigated to determine their potential as next-generation remedies for Pb poisoning.

Comprehensive Theoretical View of the [Cu2 O2 ] Side-on-Peroxo-/Bis-µ-Oxo Equilibria.

Coupled binuclear copper (CBC) sites are employed by many metalloenzymes to catalyze a broad set of biochemical transformations. Typically, the CBC catalytic sites are activated by the O2 molecule to form various [Cu2 O2 ] reactive species. This has also inspired synthesis and development of various biomimetic inorganic complexes featuring the CBC core. From theoretical perspective, the [Cu2 O2 ] reactivity often hinges on the side-on-peroxo-dicopper(II) (P) vs. bis-µ-oxo-dicopper(III) (O) isomerism - an equilibrium that has become almost iconic in theoretical bioinorganic chemistry. Herein, we present a comprehensive calibration and evaluation of the performance of various composite computational protocols available in contemporary computational chemistry, involving coupled-cluster and multireference (relativistic) wave function methods, popular density functionals and solvation models. Starting with the well-studied reference [Cu2 O2 (NH3 )6 ]2+ system, we compared the performance of electronic structure methods and discussed the relativistic effects. This allowed us to select several 'calibrated' DFT functionals that can be conveniently employed to study ten experimentally well-characterized [Cu2 O2 ] inorganic systems. We mostly predicted the lowest-energy structures (P vs. O) of the studied systems correctly. In addition, we present calibration of the used electronic structure methods for prediction of the spectroscopic features of the [Cu2 O2 ] core, mostly provided by the resonance Raman (rR) spectroscopy.

Understanding desaturation/hydroxylation activity of castor stearoyl Δ9-Desaturase through rational mutagenesis.

A recently proposed reaction mechanism of soluble Δ9 desaturase (Δ9D) allowed us to identify auxiliary residues His203, Asp101, Thr206 and Cys222 localized near the di-iron active site that are supposedly involved in the proton transfer (PT) to and from the active site. The PT, along with the electron transfer (ET), seems to be crucial for efficient desaturation. Thus, perturbing the major PT chains is expected to impair the native reaction and (potentially) amplify minor reaction channels, such as the substrate hydroxylation. To verify this hypothesis, we mutated the four residues mentioned above into their counterparts present in a soluble methane monooxygenase (sMMO), and determined the reaction products of mutants. We found that the mutations significantly promote residual monohydroxylation activities on stearoyl-CoA, often at the expense of native desaturation activity. The favored hydroxylation positions are C9, followed by C10 and C11. Reactions with unsaturated substrate, oleoyl-CoA, yield erythro-9,10-diol, cis-9,10-epoxide and a mixture of allylic alcohols. Additionally, using 9- and 11-hydroxystearoyl-CoA, we showed that the desaturation reaction can proceed only with the hydroxyl group at position C11, whereas the hydroxylation reaction is possible in both cases, i.e. with hydroxyl at position C9 or C11. Despite the fact that the overall outcome of hydroxylation is rather modest and that it is mostly the desaturation/hydroxylation ratio that is affected, our results broaden understanding of the origin of chemo- and stereoselectivity of the Δ9D and provide further insight into the catalytic action of the NHFe2 enzymes.

Evidence for H-bonding interactions to the µ-η2:η2-peroxide of oxy-tyrosinase that activate its coupled binuclear copper site.

The factors that control the diverse reactivity of the µ-η2:η2-peroxo dicopper(II) oxy-intermediates in the coupled binuclear copper proteins remain elusive. Here, spectroscopic and computational methods reveal H-bonding interactions between active-site waters and the µ-η2:η2-peroxide of oxy-tyrosinase, and define their effects on the Cu(II)2O2 electronic structure and O2 activation.

Predicting Effects of Site-Directed Mutagenesis on Enzyme Kinetics by QM/MM and QM Calculations: A Case of Glutamate Carboxypeptidase II.

Quantum and molecular mechanics (QM/MM) and QM-only (cluster model) modeling techniques represent the two workhorses in mechanistic understanding of enzyme catalysis. One of the stringent tests for QM/MM and/or QM approaches is to provide quantitative answers to real-world biochemical questions, such as the effect of single-point mutations on enzyme kinetics. This translates into predicting the relative activation energies to 1-2 kcal·mol-1 accuracy; such predictions can be used for the rational design of novel enzyme variants with desired/improved characteristics. Herein, we employ glutamate carboxypeptidase II (GCPII), a dizinc metallopeptidase, also known as the prostate specific membrane antigen, as a model system. The structure and activity of this major cancer antigen have been thoroughly studied, both experimentally and computationally, which makes it an ideal model system for method development. Its reaction mechanism is quite well understood: the reaction coordinate comprises a "tetrahedral intermediate" and two transition states and experimental activation Gibbs free energy of ∼17.5 kcal·mol-1 can be inferred for the known kcat ≈ 1 s-1. We correlate experimental kinetic data (including the E424H variant, newly characterized in this work) for various GCPII mutants (kcat = 8.6 × 10-5 s-1 to 2.7 s-1) with the energy profiles calculated by QM/MM and QM-only (cluster model) approaches. We show that the near-quantitative agreement between the experimental values and the calculated activation energies (ΔH⧧) can be obtained and recommend the combination of the two protocols: QM/MM optimized structures and cluster model (QM) energetics. The trend in relative activation energies is mostly independent of the QM method (DFT functional) used. Last but not least, a satisfactory correlation between experimental and theoretical data allows us to provide qualitative and fairly simple explanations of the observed kinetic effects which are thus based on a rigorous footing.