DszA Catalyzes C-S Bond Cleavage through N5-Hydroperoxyl Formation.

Due to its detrimental impact on human health and the environment, regulations demand ultralow sulfur levels on fossil fuels, in particular in diesel. However, current desulfurization techniques are expensive and cannot efficiently remove heteroaromatic sulfur compounds, which are abundant in crude oil and concentrate in the diesel fraction after distillation. Biodesulfurization via the four enzymes of the metabolic 4S pathway of the bacterium Rhodococcus erythropolis (DszA-D) is a possible solution. However, the 4S pathway needs to operate at least 500 times faster for industrial applicability, a goal currently pursued through enzyme engineering. In this work, we unveil the catalytic mechanism of the flavin monooxygenase DszA. Surprisingly, we found that this enzyme follows a recently proposed atypical mechanism that passes through the formation of an N5OOH intermediate at the re side of the cofactor, aided by a well-defined, predominantly hydrophobic O2 pocket. Besides clarifying the unusual chemical mechanism of the complex DszA enzyme, with obvious implications for understanding the puzzling chemistry of flavin-mediated catalysis, the result is crucial for the rational engineering of DszA, contributing to making biodesulfurization attractive for the oil refining industry.

The Challenging World of Simple Inorganic Rings: Revisiting Roesky's Ketone and Roesky's Sulfoxide.

The surprising differences between the experimental solid-state and calculated gas-phase structures of 5-oxo-1,3,2,4-dithiadiazole (Roesky's ketone, 1) and 1-oxo-1,2,4,3,5-trithiadiazole (Roesky's sulfoxide, 2), identified and studied in a series of papers published between 2004 and 2010 but then never satisfactorily explained, have been revisited, making use of the more advanced computational possibilities currently available. The previous calculations' considerable overestimations of the C-S and S-S bond lengths in 1 and 2, respectively, have been partly explained based on the results of periodic calculations and the application of Valence Bond (VB) Theory. In the case of 1, the crystal environment appears to stabilize a structure with a highly polarized C=O bond, which features a C-S bond with considerable double-bond character - an effect which does not exist for the isolated molecule - explaining the much shorter bond in the solid state. For 2, a similar conclusion can be drawn for the S-S distance. For both compounds, though, packing effects are not the sole source of the differences: the inability of Density Functional Theory (DFT) to properly deal with the electronic structures of these apparently simple main-group systems remains a contributing factor.

The Halogen Bond in Weakly Bonded Complexes and the Consequences for Aromaticity and Spin-Orbit Coupling.

The halogen bond complexes CF3X⋯Y and C2F3X⋯Y, with Y = furan, thiophene, selenophene and X = Cl, Br, I, have been studied by using DFT and CCSD(T) in order to understand which factors govern the interaction between the halogen atom X and the aromatic ring. We found that PBE0-dDsC/QZ4P gives an adequate description of the interaction energies in these complexes, compared to CCSD(T) and experimental results. The interaction between the halogen atom X and the π-bonds in perpendicular orientation is stronger than the interaction with the in-plane lone pairs of the heteroatom of the aromatic cycle. The strength of the interaction follows the trend Cl < Br < I; the chalcogenide in the aromatic ring nor the hybridization of the C−X bond play a decisive role. The energy decomposition analysis shows that the interaction energy is dominated by all three contributions, viz., the electrostatic, orbital, and dispersion interactions: not one factor dominates the interaction energy. The aromaticity of the ring is undisturbed upon halogen bond formation: the π-ring current remains equally strong and diatropic in the complex as it is for the free aromatic ring. However, the spin-orbit coupling between the singlet and triplet π→π* states is increased upon halogen bond formation and a faster intersystem crossing between these states is therefore expected.

Different Enzyme Conformations Induce Different Mechanistic Traits in HIV-1 Protease.

The influence of the dynamical flexibility of enzymes on reaction mechanisms is a cornerstone in biological sciences. In this study, we aim to 1) study the convergence of the activation free energy by using the first step of the reaction catalysed by HIV-1 protease as a case study, and 2) provide further evidence for a mechanistic divergence in this enzyme, as two different reaction pathways were seen to contribute to this step. We used quantum mechanics/molecular mechanics molecular dynamics simulations, on four different initial conformations that led to different barriers in a previous study. Despite the sampling, the four activation free energies still spanned a range of 5.0 kcal ⋅ mol-1 . Furthermore, the new simulations did confirm the occurrence of an unusual mechanistic divergence, with two different mechanistic pathways displaying equivalent barriers. An active-site water molecule is proposed to influence the mechanistic pathway.

Towards the Accurate Thermodynamic Characterization of Enzyme Reaction Mechanisms.

We employed QM/MM molecular dynamics (MD) simulations to characterize the rate-limiting step of the glycosylation reaction of pancreatic α-amylase with combined DFT/molecular dynamics methods (PBE/def2-SVP : AMBER). Upon careful choice of four starting active site conformations based on thorough reactivity criteria, Gibbs energy profiles were calculated with umbrella sampling simulations within a statistical convergence of 1-2 kcal ⋅ mol-1 . Nevertheless, Gibbs activation barriers and reaction energies still varied from 11.0 to 16.8 kcal ⋅ mol-1 and -6.3 to +3.8 kcal ⋅ mol-1 depending on the starting conformations, showing that despite significant state-of-the-art QM/MM MD sampling (0.5 ns/profile) the result still depends on the starting structure. The results supported the one step dissociative mechanism of Asp197 glycosylation preceded by an acid-base reaction by the Glu233, which are qualitatively similar to those from multi-PES QM/MM studies, and thus support the use of the latter to determine enzyme reaction mechanisms.

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

Electronic couplings for singlet fission: Orbital choice and extrapolation to the complete basis set limit.

For the search for promising singlet fission candidates, the calculation of the effective electronic coupling, which is required to estimate the singlet fission rate between the initially excited state (S0S1) and the multiexcitonic state (1TT, two triplets on neighboring molecules, coupled into a singlet), should be sufficiently reliable and fast enough to explore the configuration space. We propose here to modify the calculation of the effective electronic coupling using a nonorthogonal configuration interaction approach by: (a) using only one set of orbitals, optimized for the triplet state of the molecules, to describe all molecular electronic states, and (b) only taking the leading configurations into consideration. Furthermore, we also studied the basis set convergence of the electronic coupling, and we found, by comparison to the complete basis set limit obtained using the cc-pVnZ series of basis sets, that both the aug-cc-pVDZ and 6-311++G** basis sets are a good compromise between accuracy and computational feasibility. The proposed approach enables future work on larger clusters of molecules than dimers.

The electronic structure of carbones revealed: insights from valence bond theory.

In this contribution, we studied the OC-C bond in carbon suboxide and related allene compounds using the valence bond method. The nature of this bond has been the subject of debate, whether it is a regular, electron sharing bond or a dative bond. We compared the nature of this bond in carbon suboxide with the gold-CO bond in Au(CO)2+, which is a typical dative bond, and we studied its charge-shift bond character. We found that the C-CO bond in carbon suboxide is unique in the sense that it cannot be assigned as either a dative or electron sharing bond, but it is an admixture of electron sharing and dative components, together with a high contribution of ionic character. These findings provide a clear basis for distinguishing the commonly found dative bonds between ligands and transition metals and the present case of what may be described as coordinative bonding to carbon.

Lessons from combined experimental and theoretical examination of the FTIR and 2D-IR spectroelectrochemistry of the amide I region of cytochrome c.

Amide I difference spectroscopy is widely used to investigate protein function and structure changes. In this article, we show that the common approach of assigning features in amide I difference signals to distinct secondary structure elements in many cases may not be justified. Evidence comes from Fourier transform infrared (FTIR) and 2D-IR spectroelectrochemistry of the protein cytochrome c in the amide I range, in combination with computational spectroscopy based on molecular dynamics (MD) simulations. This combination reveals that each secondary structure unit, such as an alpha-helix or a beta-sheet, exhibits broad overlapping contributions, usually spanning a large part of the amide I region, which in the case of difference absorption experiments (such as in FTIR spectroelectrochemistry) may lead to intensity-compensating and even sign-changing contributions. We use cytochrome c as the test case, as this small electron-transferring redox-active protein contains different kinds of secondary structure units. Upon switching its redox-state, the protein exhibits a different charge distribution while largely retaining its structural scaffold. Our theoretical analysis suggests that the change in charge distribution contributes to the spectral changes and that structural changes are small. However, in order to confidently interpret FTIR amide I difference signals in cytochrome c and proteins in general, MD simulations in combination with additional experimental approaches such as isotope labeling, the insertion of infrared labels to selectively probe local structural elements will be required. In case these data are not available, a critical assessment of previous interpretations of protein amide I 1D- and 2D-IR difference spectroscopy data is warranted.

Total Synthesis of the Alleged Structure of Crenarchaeol Enables Structure Revision*.

Crenarchaeol is a glycerol dialkyl glycerol tetraether lipid produced exclusively in Archaea of the phylum Thaumarchaeota. This membrane-spanning lipid is undoubtedly the structurally most sophisticated of all known archaeal lipids and an iconic molecule in organic geochemistry. The 66-membered macrocycle possesses a unique chemical structure featuring 22 mostly remote stereocenters, and a cyclohexane ring connected by a single bond to a cyclopentane ring. Herein we report the first total synthesis of the proposed structure of crenarchaeol. Comparison with natural crenarchaeol allowed us to propose a revised structure of crenarchaeol, wherein one of the 22 stereocenters is inverted.

Hydrophobic Collapse in N-Methylacetamide-Water Mixtures.

Aqueous N-methylacetamide solutions were investigated by polarization-resolved pump-probe and 2D infrared spectroscopy (2D IR), using the amide I mode as a reporter. The 2D IR results are compared with molecular dynamics simulations and spectral calculations to gain insight into the molecular structures in the mixture. N-Methylacetamide and water molecules tend to form clusters with "frozen" amide I dynamics. This is driven by a hydrophobic collapse as the methyl groups of the N-methylacetamide molecules cluster in the presence of water. Since the studied system can be considered as a simplified model for the backbone of proteins, the present study forms a convenient basis for understanding the structural and vibrational dynamics in proteins. It is particularly interesting to find out that a hydrophobic collapse as the one driving protein folding is observed in such a simple system.

Hydrogen bond and lifetime dynamics in diluted alcohols.

Hydrogen-bonding plays a crucial role in many chemical and biochemical reactions. Alcohols, with their hydrophilic and hydrophobic groups, constitute an important class of hydrogen-bonding molecules with functional tuning possibilities through changes in the hydrophobic tails. Recent studies demonstrated that for solutions of alcohols changes in the hydrophobic tail significantly affect a broad range of dynamics properties of the liquid. Still, the understanding is lacking on the origin of such differences in terms of a solvent- versus a solute-dominated effect. Here we reveal this origin by studying hydrogen-bond dynamics in a number of alcohol molecules - from methanol to butanol - diluted in a hydrogen-bond accepting environment, acetonitrile. The dynamics were investigated by pump-probe and 2D infrared spectroscopy combined with molecular dynamics-spectral simulations, using the OH stretching mode as a reporter. For all the considered alcohols, the vibrational lifetime of the OH stretching mode was found to be ∼3 ps. The hydrogen-bond dynamics exhibit similar behavior with a fast (∼200 fs) initial relaxation dominated by librational motion and a slow (∼4 ps) relaxation due to hydrogen-bond exchange dynamics. The similar dynamics over such a broad range of alcohols led us to conclude that the previously observed differences in dynamics in bulk alcohols originate from the dependence of the solvent properties on the hydrophobic tail, while the solute properties as found herein are essentially independent of the hydrophobic tail.

2D IR spectroscopy of high-pressure phases of ice.

We present experimental and simulated 2D IR spectra of some high-pressure forms of isotope-pure D2O ice and compare the results to those of ice Ih published previously [F. Perakis and P. Hamm, Phys. Chem. Chem. Phys. 14, 6250 (2012); L. Shi et al., ibid. 18, 3772 (2016)]. Ice II, ice V, and ice XIII have been chosen for this study, since this selection covers many aspects of the polymorphism of ice. That is, ice II is a hydrogen-ordered phase of ice, in contrast to ice Ih, while ice V and ice XIII are a hydrogen-disordered/ordered couple that shares essentially the same oxygen structure and hydrogen-bonded network. For the transmission 2D IR spectroscopy, a novel method had to be developed for the preparation of ultrathin films (1-2 µm) of high-pressure ices with good optical quality. We also simulated 2D IR spectra based on molecular dynamics simulations connected to a vibrational exciton picture. These simulations agree with the experimental results in a semi-quantitative manner for ice II, while the same approach failed for ice V and ice XIII. From the perspective of 2D IR spectroscopy, ice II appears to be more inhomogeneously broadened than ice Ih, despite its hydrogen-order, which we attribute to the fact that ice II is structurally more complex with four distinguishable hydrogen bonds that mix due to exciton coupling. Ice V and ice XIII, on the other hand, behave as expected with the hydrogen-disordered case (ice V) being more inhomogenously broadened. Furthermore, in all hydrogen-ordered forms (ice II and ice XIII), cross peaks could be identified in the anisotropic 2D IR spectrum, whose signs reveal the relative direction of the corresponding excitonic states.

Hydrogen bond dynamics in bulk alcohols.

Hydrogen-bonded liquids play a significant role in numerous chemical and biological phenomena. In the past decade, impressive developments in multidimensional vibrational spectroscopy and combined molecular dynamics-quantum mechanical simulation have established many intriguing features of hydrogen bond dynamics in one of the fundamental solvents in nature, water. The next class of a hydrogen-bonded liquid--alcohols--has attracted much less attention. This is surprising given such important differences between water and alcohols as the imbalance between the number of hydrogen bonds, each molecule can accept (two) and donate (one) and the very presence of the hydrophobic group in alcohols. Here, we use polarization-resolved pump-probe and 2D infrared spectroscopy supported by extensive theoretical modeling to investigate hydrogen bond dynamics in methanol, ethanol, and isopropanol employing the OH stretching mode as a reporter. The sub-ps dynamics in alcohols are similar to those in water as they are determined by similar librational and hydrogen-bond stretch motions. However, lower density of hydrogen bond acceptors and donors in alcohols leads to the appearance of slow diffusion-controlled hydrogen bond exchange dynamics, which are essentially absent in water. We anticipate that the findings herein would have a potential impact on fundamental chemistry and biology as many processes in nature involve the interplay of hydrophobic and hydrophilic groups.

Changes in aromaticity of spin-crossover complexes: a signature for non-innocent ligands.

The influence of the spin state of the metal centre in spin crossover compounds on the aromaticity of the ligands has been investigated for iron(II)tris-bipyridine (Fe(bpy)32+), and Fe(II)(formazanate)2 (as a truncated model and the full phenyl substituted compound). It was found that the aromaticity of the bipyridine ligands is unaffected by changing the spin state of the central iron atom, but that of the formazanate ligands is reduced upon transition to the high-spin state. This change in aromaticity is rationalized using the symmetry selection rules for aromaticity in terms of virtual excitations from occupied to empty orbitals. A further consequence of this loss in aromaticity is a shift to higher energy in the ring vibrations of the formazanate compounds that can be observed in either its IR or Raman spectrum; this prediction has been confirmed here. This change in aromaticity as a consequence of change in spin state can be regarded as an indication for non-innocent ligands.

Revisiting the reaction pathways for phospholipid hydrolysis catalyzed by phospholipase A2 with QM/MM methods.

Secreted phospholipase A2 (sPLA2) is a Ca2+-dependent, widely distributed enzyme superfamily in almost all mammalian tissues and bacteria. It is also a critical component of the venom of nearly all snakes, as well as many invertebrate species. In non-venomous contexts, sPLA2 assumes significance in cellular signaling pathways by binding cell membranes and catalyzing ester bond hydrolysis at the sn-2 position of phospholipids. Elevated levels of GIIA sPLA2 have been detected in the synovial fluid of arthritis patients, where it exhibits a pro-inflammatory function. Consequently, identifying sPLA2 inhibitors holds promise for creating highly effective pharmaceutical treatments. Beyond arthritis, the similarities among GIIA sPLA2s offer an opportunity for developing treatments against snakebite envenoming, the deadliest neglected tropical disease. Despite decades of study, the details of PLA2 membrane-binding, substrate-binding, and reaction mechanism remain elusive, demanding a comprehensive understanding of the sPLA2 catalytic mechanism. This study explores two reaction mechanism hypotheses, involving one or two water molecules, and distinct roles for the Ca2+ cofactor. Our research focuses on the human synovial sPLA2 enzyme bound to lipid bilayers of varying phospholipid compositions, and employing adiabatic QM/MM and QM/MM MD umbrella sampling methods to energetically and geometrically characterize the structures found along both reaction pathways. Our studies demonstrate the higher frequency of productive conformations within the single-water pathway, also revealing a lower free energy barrier for hydrolyzing POPC. Furthermore, we observe that the TS of this concerted one-step reaction closely resembles transition state geometries observed in X-ray crystallography complexes featuring high-affinity transition state analogue inhibitors.

Photoactivation of titanium-oxo cluster [Ti6O6(OR)6(O2C t Bu)6]: mechanism, photoactivated structures, and onward reactivity with O2 to a peroxide complex.

The molecular titanium-oxo cluster [Ti6O6(OiPr)6(O2C t Bu)6] (1) can be photoactivated by UV light, resulting in a deeply coloured mixed valent (photoreduced) Ti (iii/iv) cluster, alongside alcohol and ketone (photooxidised) organic products. Mechanistic studies indicate that a two-electron (not free-radical) mechanism occurs in this process, which utilises the cluster structure to facilitate multielectron reactions. The photoreduced products [Ti6O6(OiPr)4(O2C t Bu)6(sol)2], sol = iPrOH (2) or pyridine (3), can be isolated in good yield and are structurally characterized, each with two, uniquely arranged, antiferromagnetically coupled d-electrons. 2 and 3 undergo onward oxidation under air, with 3 cleanly transforming into peroxide complex, [Ti6O6(OiPr)4(O2C t Bu)6(py)(O2)] (5). 5 reacts with isopropanol to regenerate the initial cluster (1) completing a closed cycle, and suggesting opportunities for the deployment of these easily made and tuneable clusters for sustainable photocatalytic processes using air and light. The redox reactivity described here is only possible in a cluster with multiple Ti sites, which can perform multi-electron processes and can adjust its shape to accommodate changes in electron density.

Temperature-Induced Collapse of Elastin-like Peptides Studied by 2DIR Spectroscopy.

Elastin-like peptides are hydrophobic biopolymers that exhibit a reversible coacervation transition when the temperature is raised above a critical point. Here, we use a combination of linear infrared spectroscopy, two-dimensional infrared spectroscopy, and molecular dynamics simulations to study the structural dynamics of two elastin-like peptides. Specifically, we investigate the effect of the solvent environment and temperature on the structural dynamics of a short (5-residue) elastin-like peptide and of a long (450-residue) elastin-like peptide. We identify two vibrational energy transfer processes that take place within the amide I' band of both peptides. We observe that the rate constant of one of the exchange processes is strongly dependent on the solvent environment and argue that the coacervation transition is accompanied by a desolvation of the peptide backbone where up to 75% of the water molecules are displaced. We also study the spectral diffusion dynamics of the valine(1) residue that is present in both peptides. We find that these dynamics are relatively slow and indicative of an amide group that is shielded from the solvent. We conclude that the coacervation transition of elastin-like peptides is probably not associated with a conformational change involving this residue.

Interplay between Hydrogen Bonding and Vibrational Coupling in Liquid N-Methylacetamide.

Intrinsically disordered proteins play an important role in biology, and unraveling their labile structure presents a vital challenge. However, the dynamical structure of such proteins thwarts their study by standard techniques such as X-ray diffraction and NMR spectroscopy. Here, we use a neat liquid composed of N-methylacetamide molecules as a model system to elucidate dynamical and structural properties similar to those one can expect to see in intrinsically disordered proteins. To examine the structural dynamics in the neat liquid, we combine molecular dynamics, response-function-based spectral simulations, and two-dimensional polarization-resolved infrared spectroscopy in the amide I (CO stretch) region. The two-dimensional spectra reveal a delicate interplay between hydrogen bonding and intermolecular vibrational coupling effects, observed through a fast anisotropy decay. The present study constitutes a general platform for understanding the structure and dynamics of highly disordered proteins.

Assessing Spectral Simulation Protocols for the Amide I Band of Proteins.

We present a benchmark study of spectral simulation protocols for the amide I band of proteins. The amide I band is widely used in infrared spectroscopy of proteins due to the large signal intensity, high sensitivity to hydrogen bonding, and secondary structural motifs. This band has, thus, proven valuable in many studies of protein structure-function relationships. We benchmark spectral simulation protocols using two common force fields in combination with several electrostatic mappings and coupling models. The results are validated against experimental linear absorption and two-dimensional infrared spectroscopy for three well-studied proteins. We find two-dimensional infrared spectroscopy to be much more sensitive to the simulation protocol than linear absorption and report on the best simulation protocols. The findings demonstrate that there is still room for ideas to improve the existing models for the amide I band of proteins.