Observation of a promethium complex in solution.

Lanthanide rare-earth metals are ubiquitous in modern technologies1-5, but we know little about chemistry of the 61st element, promethium (Pm)6, a lanthanide that is highly radioactive and inaccessible. Despite its importance7,8, Pm has been conspicuously absent from the experimental studies of lanthanides, impeding our full comprehension of the so-called lanthanide contraction phenomenon: a fundamental aspect of the periodic table that is quoted in general chemistry textbooks. Here we demonstrate a stable chelation of the 147Pm radionuclide (half-life of 2.62 years) in aqueous solution by the newly synthesized organic diglycolamide ligand. The resulting homoleptic PmIII complex is studied using synchrotron X-ray absorption spectroscopy and quantum chemical calculations to establish the coordination structure and a bond distance of promethium. These fundamental insights allow a complete structural investigation of a full set of isostructural lanthanide complexes, ultimately capturing the lanthanide contraction in solution solely on the basis of experimental observations. Our results show accelerated shortening of bonds at the beginning of the lanthanide series, which can be correlated to the separation trends shown by diglycolamides9-11. The characterization of the radioactive PmIII complex in an aqueous environment deepens our understanding of intra-lanthanide behaviour12-15 and the chemistry and separation of the f-block elements16.

Massively parallel GPU enabled third-order cluster perturbation excitation energies for cost-effective large scale excitation energy calculations.

We present here a massively parallel implementation of the recently developed CPS(D-3) excitation energy model that is based on cluster perturbation theory. The new algorithm extends the one developed in Baudin et al. [J. Chem. Phys., 150, 134110 (2019)] to leverage multiple nodes and utilize graphical processing units for the acceleration of heavy tensor contractions. Furthermore, we show that the extended algorithm scales efficiently with increasing amounts of computational resources and that the developed code enables CPS(D-3) excitation energy calculations on large molecular systems with a low time-to-solution. More specifically, calculations on systems with over 100 atoms and 1000 basis functions are possible in a few hours of wall clock time. This establishes CPS(D-3) excitation energies as a computationally efficient alternative to those obtained from the coupled-cluster singles and doubles model.

[Fe4S4] cubane in sulfite reductases: new insights into bonding properties and reactivity.

The dissimilatory sulfite reductase enzyme has very characteristic active site where the substrate binds to an iron site, ligated by a siroheme macrocycle and a thiol directly connected to a [Fe4S4] cluster. This arrangement gives the enzyme remarkable efficiency in reducing sulfite and nitrite all the way to hydrogen sulfide and ammonia. For the first time we present a theoretical study where substrate binding modalities and activation are elucidated using active site models containing proton supply side chains and the [Fe4S4] cluster. Density functional theory (DFT) was deployed in conjunction with the energy decomposition scheme (as implemented in AMS), the quantum theory of atoms in molecules (QTAIM), and conceptual DFT (cDFT) descriptors. We quantified the role of the electrostatic interactions inside the active site created by the side chains as well as the influence of the [Fe4S4] cluster on the substrate binding. Furthermore, using conceptual DFT results we shed light of the activation process, thus, laying foundation for further mechanistic studies. We found that the bonding of the ligands to the iron complex is dominated by electrostatic interactions, but the presence of the [Fe4S4] cubane leads to substantial changes in electronic interaction. The spin state of the cubane, however, affects the binding energy only marginally. The conceptual DFT results show that the presence of the [Fe4S4] cubane affects the reactivity of the active site as it is involved in electron transfer. This is corroborated by an increase in the electrophilicity index, thus making the active site more prone to react with the ligands. The interaction energies between the ligand and the siroheme group are also increased upon the presence of the cubane group, thus, suggesting that the siroheme group is not an innocent spectator but plays an active role in the reactivity of the dSIR active site.

Multilayer Divide-Expand-Consolidate Coupled-Cluster Method: Demonstrative Calculations of the Adsorption Energy of Carbon Dioxide in the Mg-MOF-74 Metal-Organic Framework.

The implementation and evaluation of a multilayer extension of the divide-expand-consolidate (DEC) scheme within the LSDalton program is presented. The DEC scheme is a linear-scaling, fragmentation-based local coupled-cluster (CC) method that provides a means of overcoming the scaling wall associated with canonical CC electronic structure calculations on large molecular systems. Taking advantage of the local nature of correlation effects, the correlation energy for the full molecule is calculated from a set of independent fragments using localized molecular orbitals. However, when only a small subsystem of a larger system is of interest, for example, adsorption sites or catalytically active sites, the majority of the computational time may be spent evaluating the correlation energy of fragments which have little effect on the properties in the area of interest (AOI). The multilayer DEC (ML-DEC) scheme addresses this by taking advantage of the independent nature of the fragments in order to evaluate the correlation energy of various regions of the system at different levels of theory. Regions far from the AOI are evaluated at lower (cheaper) levels of theory such as Hartree-Fock (HF) or Møller-Plesset second-order perturbation theory (MP2), while the area immediately surrounding the AOI is treated with a higher level CC model. Through the ML-DEC scheme, the computational cost of CC calculations on these types of systems can be significantly reduced while maintaining the accuracy of higher-level calculations. Results from HF/RI-MP2 and RI-MP2/CCSD ML-DEC calculations of the binding energy of a fatty acid dimer are presented. We find that the ML-DEC scheme is capable of reproducing DEC energy differences at a target level of theory, provided that the region treated at the target level of theory is chosen to be sufficiently large. Time-to-solution is found to be significantly reduced, particularly in the RI-MP2/CCSD calculations. Finally, the ML-DEC scheme is applied to the calculation of CO2 adsorption in a Mg-MOF-74 channel.

Cluster perturbation theory. III. Perturbation series for coupled cluster singles and doubles excitation energies.

The cluster perturbation series, CPS(D), for coupled cluster singles and doubles excitation energies is considered. It is demonstrated that the second-order model CPS(D-2) is identical to the configuration interaction singles with perturbative doubles, CIS(D) model. The third-order model, CPS(D-3), provides excitation energies of coupled cluster singles and doubles (CCSD) quality in the sense that the difference between CPS(D-3) and CCSD excitation energies is of the same size or smaller than the effect of adding triples corrections to CCSD excitation energies. We further show that the third-order corrections can be efficiently implemented, in particular, when the resolution of the identity approximation is used for integrals. We also show that the CPS(D-3) excitation energies can be determined for system sizes that are far beyond what can be considered in conventional CCSD excitation energy calculations.

The GPU-enabled divide-expand-consolidate RI-MP2 method (DEC-RI-MP2).

We report porting of the Divide-Expand-Consolidate Resolution of the Identity second-order Møller-Plesset perturbation (DEC-RI-MP2) method to the graphic processing units (GPUs) using OpenACC compiler directives. It is shown that the OpenACC compiler directives implementation efficiently accelerates the rate-determining step of the DEC-RI-MP2 method with minor implementation effort. Moreover, the GPU acceleration results in a better load balance and thus in an overall scaling improvement of the DEC algorithm. The resulting cross-platform hybrid MPI/OpenMP/OpenACC implementation has scalable and portable performance on heterogeneous HPC architectures. The GPU-enabled code was benchmarked using a reduced version of the S12L test set of Stefan Grimme (Grimme, Chem. Eur. J. 2012, 18, 9955) consisting of supramolecular complexes up to 158 atoms and 4292 contracted basis functions (cc-pVTZ). The test set results demonstrate the general applicability of the DEC-RI-MP2 method showing results consistent with the DEC-RI-MP2 introductory paper (Baudin et al., J. Chem. Phys. 2016, 144, 054102) on molecules of complicated electronic structures. © 2016 Wiley Periodicals, Inc.

The molecular gradient using the divide-expand-consolidate resolution of the identity second-order Møller-Plesset perturbation theory: The DEC-RI-MP2 gradient.

We report an implementation of the molecular gradient using the divide-expand-consolidate resolution of the identity second-order Møller-Plesset perturbation theory (DEC-RI-MP2). The new DEC-RI-MP2 gradient method combines the precision control as well as the linear-scaling and massively parallel features of the DEC scheme with efficient evaluations of the gradient contributions using the RI approximation. We further demonstrate that the DEC-RI-MP2 gradient method is capable of calculating molecular gradients for very large molecular systems. A test set of supramolecular complexes containing up to 158 atoms and 1960 contracted basis functions has been employed to demonstrate the general applicability of the DEC-RI-MP2 method and to analyze the errors of the DEC approximation. Moreover, the test set contains molecules of complicated electronic structures and is thus deliberately chosen to stress test the DEC-RI-MP2 gradient implementation. Additionally, as a showcase example the full molecular gradient for insulin (787 atoms and 7604 contracted basis functions) has been evaluated.

Electronic structure of the unique [4Fe-3S] cluster in O2-tolerant hydrogenases characterized by 57Fe Mossbauer and EPR spectroscopy.

Iron-sulfur clusters are ubiquitous electron transfer cofactors in hydrogenases. Their types and redox properties are important for H(2) catalysis, but, recently, their role in a protection mechanism against oxidative inactivation has also been recognized for a [4Fe-3S] cluster in O(2)-tolerant group 1 [NiFe] hydrogenases. This cluster, which is uniquely coordinated by six cysteines, is situated in the proximity of the catalytic [NiFe] site and exhibits unusual redox versatility. The [4Fe-3S] cluster in hydrogenase (Hase) I from Aquifex aeolicus performs two redox transitions within a very small potential range, forming a superoxidized state above +200 mV vs. standard hydrogen electrode (SHE). Crystallographic data has revealed that this state is stabilized by the coordination of one of the iron atoms to a backbone nitrogen. Thus, the proximal [4Fe-3S] cluster undergoes redox-dependent changes to serve multiple purposes beyond classical electron transfer. In this paper, we present field-dependent (57)Fe-Mössbauer and EPR data for Hase I, which, in conjunction with spectroscopically calibrated density functional theory (DFT) calculations, reveal the distribution of Fe valences and spin-coupling schemes for the iron-sulfur clusters. The data demonstrate that the electronic structure of the [4Fe-3S] core in its three oxidation states closely resembles that of corresponding conventional [4Fe-4S] cubanes, albeit with distinct differences for some individual iron sites. The medial and distal iron-sulfur clusters have similar electronic properties as the corresponding cofactors in standard hydrogenases, although their redox potentials are higher.

Six-electron reduction of nitrite to ammonia by cytochrome c nitrite reductase: insights from density functional theory studies.

In this Forum Article, an extensive discussion of the mechanism of six-electron, seven-proton nitrite reduction by the cytochrome c nitrite reductase enzyme is presented. On the basis of previous studies, the entire mechanism is summarized and a unified picture of the most plausible sequence of elementary steps is presented. According to this scheme, the mechanism can be divided into five functional stages. The first phase of the reaction consists of substrate binding and N-O bond cleavage. Here His277 plays a crucial role as a proton donor. In this step, the N-O bond is cleaved heterolytically through double protonation of the substrate. The second phase of the mechanism consists of two proton-coupled electron-transfer events, leading to an HNO intermediate. The third phase involves the formation of hydroxylamine, where Arg114 provides the necessary proton for the reaction. The second N-O bond is cleaved in the fourth phase of the mechanism, again triggered by proton transfer from His277. The Tyr218 side chain governs the fifth and last phase of the mechanism. It consists of radical transfer and ammonia formation. Thus, this mechanism implies that all conserved active-site side chains work in a concerted way in order to achieve this complex chemical transformation from nitrite to ammonia. The Forum Article also provides a detailed discussion of the density functional theory based cluster model approach to bioinorganic reactivity. A variety of questions are considered: the resting state of enzyme and substrate binding modes, interaction with the metal site and with active-site side chains, electron- and proton-transfer events, substrate dissociation, etc.

Heme-bound nitroxyl, hydroxylamine, and ammonia ligands as intermediates in the reaction cycle of cytochrome c nitrite reductase: a theoretical study.

In this article, we consider, in detail, the second half-cycle of the six-electron nitrite reduction mechanism catalyzed by cytochrome c nitrite reductase. In total, three electrons and four protons must be provided to reach the final product, ammonia, starting from the HNO intermediate. According to our results, the first event in this half-cycle is the reduction of the HNO intermediate, which is accomplished by two PCET reactions. Two isomeric radical intermediates, HNOH(•) and H2NO(•), are formed. Both intermediates are readily transformed into hydroxylamine, most likely through intramolecular proton transfer from either Arg114 or His277. An extra proton must enter the active site of the enzyme to initiate heterolytic cleavage of the N-O bond. As a result of N-O bond cleavage, the H2N(+) intermediate is formed. The latter readily picks up an electron, forming H2N(+•), which in turn reacts with Tyr218. Interestingly, evidence for Tyr218 activity was provided by the mutational studies of Lukat (Biochemistry 47:2080, 2008), but this has never been observed in the initial stages of the overall reduction process. According to our results, an intramolecular reaction with Tyr218 in the final step of the nitrite reduction process leads directly to the final product, ammonia. Dissociation of the final product proceeds concomitantly with a change in spin state, which was also observed in the resonance Raman investigations of Martins et al. (J Phys Chem B 114:5563, 2010).

Coupled cluster theory on modern heterogeneous supercomputers.

This study examines the computational challenges in elucidating intricate chemical systems, particularly through ab-initio methodologies. This work highlights the Divide-Expand-Consolidate (DEC) approach for coupled cluster (CC) theory-a linear-scaling, massively parallel framework-as a viable solution. Detailed scrutiny of the DEC framework reveals its extensive applicability for large chemical systems, yet it also acknowledges inherent limitations. To mitigate these constraints, the cluster perturbation theory is presented as an effective remedy. Attention is then directed towards the CPS (D-3) model, explicitly derived from a CC singles parent and a doubles auxiliary excitation space, for computing excitation energies. The reviewed new algorithms for the CPS (D-3) method efficiently capitalize on multiple nodes and graphical processing units, expediting heavy tensor contractions. As a result, CPS (D-3) emerges as a scalable, rapid, and precise solution for computing molecular properties in large molecular systems, marking it an efficient contender to conventional CC models.

Corrigendum: Coupled cluster theory on modern heterogeneous supercomputers.

[This corrects the article DOI: 10.3389/fchem.2023.1154526.].

Reductive activation of the heme iron-nitrosyl intermediate in the reaction mechanism of cytochrome c nitrite reductase: a theoretical study.

Cytochrome c nitrite reductase catalyzes the six-electron, seven-proton reduction of nitrite to ammonia without release of any detectable reaction intermediate. This implies a unique flexibility of the active site combined with a finely tuned proton and electron delivery system. In the present work, we employed density functional theory to study the recharging of the active site with protons and electrons through the series of reaction intermediates based on nitrogen monoxide [Fe(II)-NO(+), Fe(II)-NO·, Fe(II)-NO(-), and Fe(II)-HNO]. The activation barriers for the various proton and electron transfer steps were estimated in the framework of Marcus theory. Using the barriers obtained, we simulated the kinetics of the reduction process. We found that the complex recharging process can be accomplished in two possible ways: either through two consecutive proton-coupled electron transfers (PCETs) or in the form of three consecutive elementary steps involving reduction, PCET, and protonation. Kinetic simulations revealed the recharging through two PCETs to be a means of overcoming the predicted deep energetic minimum that is calculated to occur at the stage of the Fe(II)-NO· intermediate. The radical transfer role for the active-site Tyr(218), as proposed in the literature, cannot be confirmed on the basis of our calculations. The role of the highly conserved calcium located in the direct proximity of the active site in proton delivery has also been studied. It was found to play an important role in the substrate conversion through the facilitation of the proton transfer steps.

Hit Expansion of a Noncovalent SARS-CoV-2 Main Protease Inhibitor.

Inhibition of the SARS-CoV-2 main protease (Mpro) is a major focus of drug discovery efforts against COVID-19. Here we report a hit expansion of non-covalent inhibitors of Mpro. Starting from a recently discovered scaffold (The COVID Moonshot Consortium. Open Science Discovery of Oral Non-Covalent SARS-CoV-2 Main Protease Inhibitor Therapeutics. bioRxiv 2020.10.29.339317) represented by an isoquinoline series, we searched a database of over a billion compounds using a cheminformatics molecular fingerprinting approach. We identified and tested 48 compounds in enzyme inhibition assays, of which 21 exhibited inhibitory activity above 50% at 20 µM. Among these, four compounds with IC50 values around 1 µM were found. Interestingly, despite the large search space, the isoquinolone motif was conserved in each of these four strongest binders. Room-temperature X-ray structures of co-crystallized protein-inhibitor complexes were determined up to 1.9 Å resolution for two of these compounds as well as one of the stronger inhibitors in the original isoquinoline series, revealing essential interactions with the binding site and water molecules. Molecular dynamics simulations and quantum chemical calculations further elucidate the binding interactions as well as electrostatic effects on ligand binding. The results help explain the strength of this new non-covalent scaffold for Mpro inhibition and inform lead optimization efforts for this series, while demonstrating the effectiveness of a high-throughput computational approach to expanding a pharmacophore library.

Dynamic hydrogen-bonding network in the distal pocket of the nitrosyl complex of Pseudomonas aeruginosa cd1 nitrite reductase.

cd(1) nitrite reductase (NIR) is a key enzyme in the denitrification process that reduces nitrite to nitric oxide (NO). It contains a specialized d(1)-heme cofactor, found only in this class of enzymes, where the substrate, nitrite, binds and is converted to NO. For a long time, it was believed that NO must be released from the ferric d(1)-heme to avoid enzyme inhibition by the formation of ferrous-nitroso complex, which was considered as a dead-end product. However, recently an enhanced rate of NO dissociation from the ferrous form, not observed in standard b-type hemes, has been reported and attributed to the unique d(1)-heme structure (Rinaldo, S.; Arcovito, A.; Brunori, M.; Cutruzzolà, F. J. Biol. Chem. 2007, 282, 14761-14767). Here, we report on a detailed study of the spatial and electronic structure of the ferrous d(1)-heme NO complex from Pseudomonas aeruginosa cd(1) NIR and two mutants Y10F and H369A/H327A in solution, searching for the unique properties that are responsible for the relatively fast release. There are three residues at the "distal" side of the heme (Tyr(10), His(327), and His(369)), and in this work we focus on the identification and characterization of possible H-bonds they can form with the NO, thereby affecting the stability of the complex. For this purpose, we have used high field pulse electron-nuclear double resonance (ENDOR) combined with density functional theory (DFT) calculations. The DFT calculations were essential for assigning and interpreting the ENDOR spectra in terms of geometric structure. We have shown that the NO in the nitrosyl d(1)-heme complex of cd(1) NIR forms H-bonds with Tyr(10) and His(369), whereas the second conserved histidine, His(327), appears to be less involved in NO H-bonding. This is in contrast to the crystal structure that shows that Tyr(10) is removed from the NO. We have also observed a larger solvent accessibility to the distal pocket in the mutants as compared to the wild-type. Moreover, it was shown that the H-bonding network within the active site is dynamic and that a change in the protonation state of one of the residues does affect the strength and position of the H-bonds formed by the others. In the Y10F mutant, His(369) is closer to the NO, whereas mutation of both distal histidines displaces Tyr(10), removing its H-bond. The implications of the H-bonding network found in terms of the complex stability and catalysis are discussed.

Substrate binding and activation in the active site of cytochrome c nitrite reductase: a density functional study.

Cytochrome c nitrite reductase is a homodimeric enzyme, containing five covalently attached c-type hemes per subunit. Four of the heme irons are bishistidine-ligated, whereas the fifth, the active site of the protein, has an unusual lysine coordination and calcium site nearby. A fascinating feature of this enzyme is that the full six-electron reduction of the nitrite is achieved without release of any detectable reaction intermediate. Moreover, the enzyme is known to work over a wide pH range. Both findings suggest a unique flexibility of the active site in the complicated six-electron, seven-proton reduction process. In the present work, we employed density functional theory to study the energetics and kinetics of the initial stages of nitrite reduction. The possible role of second-sphere active-site amino acids as proton donors was investigated by taking different possible protonation states and geometrical conformations into account. It was found that the most probable proton donor is His(277), whose spatial orientation and fine-tuned acidity lead to energetically feasible, low-barrier protonation reactions. However, substrate protonation may also be accomplished by Arg(114). The calculated barriers for this pathway are only slightly higher than the experimentally determined value of 15.2 kcal/mol for the rate-limiting step. Hence, having proton-donating side chains of different acidity within the active site may increase the operational pH range of the enzyme. Interestingly, Tyr(218), which was proposed to play an important role in the overall mechanism, appears not to take part in the reaction during the initial stage.

Expanding manganese(IV) aqueous chemistry: unusually stable water-soluble hexahydrazide clathrochelate complexes.

Mn cage complexes are rare, and the ones successfully isolated in the solid state are not stable in water and organic solvents. Herein, we present the first report of mononuclear Mn clathrochelates, in which the encapsulated metal exists in the oxidation state +4. The complexes are extremely stable in the crystalline state and in solutions and show rich redox chemistry.

Activation of olefins via asymmetric Brønsted acid catalysis.

The activation of olefins for asymmetric chemical synthesis traditionally relies on transition metal catalysts. In contrast, biological enzymes with Brønsted acidic sites of appropriate strength can protonate olefins and thereby generate carbocations that ultimately react to form natural products. Although chemists have recently designed chiral Brønsted acid catalysts to activate imines and carbonyl compounds, mimicking these enzymes to protonate simple olefins that then engage in asymmetric catalytic reactions has remained a substantial synthetic challenge. Here, we show that a class of confined and strong chiral Brønsted acids enables the catalytic asymmetric intramolecular hydroalkoxylation of unbiased olefins. The methodology gives rapid access to biologically active 1,1-disubstituted tetrahydrofurans, including (-)-Boivinianin A.

Indefinitely stable iron(IV) cage complexes formed in water by air oxidation.

In nature, iron, the fourth most abundant element of the Earth's crust, occurs in its stable forms either as the native metal or in its compounds in the +2 or +3 (low-valent) oxidation states. High-valent iron (+4, +5, +6) compounds are not formed spontaneously at ambient conditions, and the ones obtained synthetically appear to be unstable in polar organic solvents, especially aqueous solutions, and this is what limits their studies and use. Here we describe unprecedented iron(IV) hexahydrazide clathrochelate complexes that are assembled in alkaline aqueous media from iron(III) salts, oxalodihydrazide and formaldehyde in the course of a metal-templated reaction accompanied by air oxidation. The complexes can exist indefinitely at ambient conditions without any sign of decomposition in water, nonaqueous solutions and in the solid state. We anticipate that our findings may open a way to aqueous solution and polynuclear high-valent iron chemistry that remains underexplored and presents an important challenge.