Search

THz absorption spectra have been recorded for the weakly bound molecular complexes of H2O with C2H4 and C2H2 embedded in cryogenic neon matrices at 2.8 K. The observation and assignment of a large-amplitude acceptor OH librational mode of the C2H2—H2O complex at 145.5 cm−1 confirms an intermolecular CH⋯O hydrogen-bonded configuration of C2v symmetry with the H2O subunit acting as the hydrogen bond acceptor. The observation and assignment of two large-amplitude donor OH librational modes of the C2H4—H2O complex at 255.0 and 187.5 cm−1, respectively, confirms an intermolecular OH⋯π hydrogen-bonded configuration with the H2O subunit acting as the hydrogen bond donor to the π-cloud of C2H4. A (semi)-empirical value for the change of vibrational zero-point energy of 4.0–4.1 kJ mol−1 is proposed and the combination with quantum chemical calculations at the CCSD(T)-F12b/aug-cc-pVQZ level provides a reliable estimate of 7.1 ± 0.3 kJ mol−1 for the dissociation energy D0 of the C2H4—H2O complex. In addition, tentative assignments for the two strongly infrared active OH librational modes of the ternary C2H4—HOH—C2H4 complex having H2O as a doubly OH⋯π hydrogen bond donor are proposed at 213.6 and 222.3 cm−1. The present findings demonstrate that the relative stability of the weak hydrogen bond motifs is not entirely rooted in differences of electronic energy but also to a large extent by differences in the vibrational zero-point energy contributions arising from the class of large-amplitude intermolecular modes.

Gene sequences form the primary basis for understanding the relationships among extant plant groups, but genetic data are unavailable from fossils to evaluate the affinities of extinct taxa. Here we show that geothermally resistant fossil cuticles of seed-bearing plants, analysed with Fourier transform infrared (FTIR) spectroscopy and hierarchical cluster analysis (HCA), retain biomolecular suites that consistently distinguish major taxa even after experiencing different diagenetic histories. Our results reveal that similarities between the cuticular biochemical signatures of major plant groups (extant and fossil) are mostly consistent with recent phylogenetic hypotheses based on molecular and morphological data. Our novel chemotaxonomic data also support the hypothesis that the extinct Nilssoniales and Bennettitales are closely allied, but only distantly related to Cycadales. The chemical signature of the cuticle of Czekanowskia (Leptostrobales) is strongly similar to that of Ginkgo leaves and supports a close evolutionary relationship between these groups. Finally, our results also reveal that the extinct putative araucariacean, Allocladus, when analysed through HCA, is grouped closer to Ginkgoales than to conifers. Thus, in the absence of modern relatives yielding molecular information, FTIR spectroscopy provides valuable proxy biochemical data complementing morphological characters to distinguish fossil taxa and to help elucidate extinct plant relationships., Additional funding from:Swedish Research Council (VR) under grant LUCCI (Lund University Carbon Cycle Centre); and the Utrecht Network Young researchers’ grant.

The true global potential energy minimum configuration of the formaldehyde dimer (CH2O)2, including the presence of a single or a double weak intermolecular CH⋯O hydrogen bond motif, has been a long-standing subject among both experimentalists and theoreticians as two different energy minima conformations of Cs and C2h symmetry have almost identical energies. The present work demonstrates how the class of large-amplitude hydrogen bond vibrational motion probed in the THz region provides excellent direct spectroscopic observables for these weak intermolecular CH⋯O hydrogen bond motifs. The combination of concentration dependency measurements, observed isotopic spectral shifts associated with H/D substitutions and dedicated annealing procedures, enables the unambiguous assignment of three large-amplitude infrared active hydrogen bond vibrational modes for the non-planar Cs configuration of (CH2O)2 embedded in cryogenic neon and enriched para-hydrogen matrices. A (semi)-empirical value for the change of vibrational zero-point energy of 5.5 ± 0.3 kJ mol−1 is proposed for the dimerization process. These THz spectroscopic observations are complemented by CCSD(T)-F12/aug-cc-pV5Z (electronic energies) and MP2/aug-cc-pVQZ (force fields) electronic structure calculations yielding a (semi)-empirical value of 13.7 ± 0.3 kJ mol−1 for the dissociation energy D0 of this global potential energy minimum.

Feathers are amongst the most complex epidermal structures known, and they have a well-documented evolutionary trajectory across non-avian dinosaurs and basal birds. Moreover, melanosome-like microbodies preserved in association with fossil plumage have been used to reconstruct original colour, behaviour and physiology. However, these putative ancient melanosomes might alternatively represent microorganismal residues, a conflicting interpretation compounded by a lack of unambiguous chemical data. We therefore used sensitive molecular imaging, supported by multiple independent analytical tests, to demonstrate that the filamentous epidermal appendages in a new specimen of the Jurassic paravian Anchiornis comprise remnant eumelanosomes and fibril-like microstructures, preserved as endogenous eumelanin and authigenic calcium phosphate. These results provide novel insights into the early evolution of feathers at the sub-cellular level, and unequivocally determine that melanosomes can be preserved in fossil feathers.

The far-infrared absorption spectra have been recorded for hydrogen-bonded complexes of water with ethanol embedded in cryogenic neon matrices at 2.8 K. The partial isotopic H/D-substitution of the ethanol subunit enabled by a dual inlet deposition procedure enables the observation and unambiguous assignment of the intermolecular high-frequency out-of-plane and the low-frequency in-plane donor OH librational modes for two different conformations of the mixed binary ethanol/water complex. The resolved donor OH librational bands confirm directly previous experimental evidence that ethanol acts as the O⋯HO hydrogen bond acceptor in the two most stable conformations. In the most stable conformation, the water subunit forces the ethanol molecule into its less stable gauche configuration upon dimerization owing to a cooperative secondary weak O⋯HC hydrogen bond interaction evidenced by a significantly blue-shift of the low-frequency in-plane donor OH librational band origin. The strong correlation between the low-frequency in-plane donor OH librational motion and the secondary intermolecular O⋯HC hydrogen bond is demonstrated by electronic structure calculations. The experimental findings are further supported by CCSD(T)-F12/aug-cc-pVQZ calculations of the conformational energy differences together with second-order vibrational perturbation theory calculations of the large-amplitude donor OH librational band origins.

Terahertz absorption spectra have been recorded for the weakly bound CO2-H2O complex embedded in cryogenic neon matrices at 2.8 K. The three high-frequency van der Waals vibrational transitions associated with out-of-plane wagging, in-plane rocking, and torsional motion of the isotopic H2O subunit have been assigned and provide crucial observables for benchmark theoretical descriptions of this systems' flat intermolecular potential energy surface. A (semi)-empirical value for the zero-point energy of 273 ± 15 cm(-1) from the class of intermolecular van der Waals vibrations is proposed and the combination with high-level quantum chemical calculations provides a value of 726 ± 15 cm(-1) for the dissociation energy D0.

Lately, there has been great interest in performing free-energy perturbation (FEP) at the combined quantum mechanics and molecular mechanics (QM/MM) level, e.g. for enzyme reactions. Such calculations require extensive sampling of phase space, which typically is prohibitive with density-functional theory or ab initio methods. Therefore, such calculations have mostly been performed with semiempirical QM (SQM) methods, or by using a thermodynamic cycle involving sampling at the MM level and perturbations between the MM and QM/MM levels of theory. However, the latter perturbations typically have convergence problems, unless the QM system is kept fixed during the simulations, because the MM and QM/MM descriptions of the internal degrees of freedom inside the QM system are too dissimilar. We have studied whether the convergence of the MM → QM/MM perturbation can be improved by using a thoroughly parameterised force field or by using SQM/MM methods. As a test case we use the first half-reaction of haloalkane dehalogenase and the QM calculations are performed with the PBE, B3LYP, and TPSSH density-functional methods. We show that the convergence can be improved with a tailored force field, but only locally around the parameterised state. Simulations based on SQM/MM methods using the MNDO, AM1, PM3, RM1, PDDG-MNDO, and PDDG-PM3 Hamiltonians have slightly better convergence properties, but very long simulations are still needed (∼10 ns) and convergence is obtained only if electrostatic interactions between the QM system and the surroundings are ignored. This casts some doubts on the common practice to base QM/MM FEPs on semiempirical simulations without any reweighting of the trajectories.

We have calculated the reorganization energy for the intramolecular electron transfer between the reduced type 1 copper site and the peroxy intermediate of the trinuclear cluster in the multicopper oxidase CueO. The calculations are performed at the combined quantum mechanics and molecular mechanics (QM/MM) level, based on molecular dynamics simulations with tailored potentials for the two copper sites. We obtain a reorganization energy of 91-133 kJ/mol, depending on the theoretical treatment. The two Cu sites contribute by 12 and 22 kJ/mol to this energy, whereas the solvent contribution is 34 kJ/mol. The rest comes from the protein, involving small contributions from many residues. We have also estimated the energy difference between the two electron-transfer states and show that the reduction of the peroxy intermediate is exergonic by 43-87 kJ/mol, depending on the theoretical method. Both the solvent and the protein contribute to this energy difference, especially charged residues close to the two Cu sites. We compare these estimates with energies obtained from QM/MM optimizations and QM calculations in a vacuum and discuss differences between the results obtained at various levels of theory.

A common approach for the computational modeling of enzyme reactions is to study a rather small model of the active site (20-200 atoms) with quantum mechanical (QM) methods, modeling the rest of the surroundings by a featureless continuum with a dielectric constant of approximately 4. In this paper, we discuss how the residues included in the QM model should be selected and how many residues need to be included before reaction energies converge. As a test case, we use a proton-transfer reaction between a first-sphere cysteine ligand and a second-sphere histidine group in the active site of [Ni,Fe] hydrogenase. We show that it is not a good approach to add groups according to their distance to the active site. A better approach is to add groups according to their contributions to the QM/MM energy difference. However, the energies can still vary by up to 50 kJ/mol for QM systems of sizes up to 230 atoms. In fact, the QM-only approach is based on the hope that a large number of sizable contributions will cancel. Interactions with neutral groups are, in general, short-ranged, with net energy contributions of less than 4 kJ/mol at distances above 5 A from the active site. Interactions with charged groups are much more long-ranged, and interactions with buried charges 20 A from the active site can still contribute by 5 kJ/mol to the reaction energy. Thus, to accurately model the influence of the surroundings on enzyme reaction energies, a detailed and unbiased atomistic account of the surroundings needs to be included.

Catalytic metal sites in enzymes frequently have second-sphere carboxylate groups that neutralize the charge of the site and share protons with first-sphere ligands. This gives rise to an ambiguity concerning the position of this proton, which has turned out to be hard to settle with experimental, as well as theoretical, methods. We study three such proton-transfer reactions in two proteins and show that, in [Ni,Fe] hydrogenase, the bridging Cys-546 ligand is deprotonated by His-79, whereas in oxidized copper nitrite reductase, the His-100 ligand is neutral and the copper-bound water molecule is deprotonated by Asp-98. We show that these reactions strongly depend on the electrostatic interactions with the surrounding protein and solvent, because there is a large change in the dipole moment of the active site (2-6 D). Neither vacuum quantum mechanical (QM) calculations with large models, a continuum solvent, or a Poisson-Boltzmann treatment of the surroundings, nor combined QM and molecular mechanics (QM/MM) optimizations give reliable estimates of the proton-transfer energies (mean absolute deviations of over 20 kJ/mol). Instead, QM/MM free-energy perturbations are needed to obtain reliable estimates of the reaction energies. These calculations also indicate what interactions and residues are important for the energy, showing how the quantum system may be systematically enlarged. With such a procedure, results with an uncertainty of similar to 10 kJ/mol can be obtained, provided that a proper QM method is used.

This chapter illustrates how quantum chemical calculations can be used to elucidate structural and functional aspects of tetrapyrrole cofactors, focusing on porphyrins, cobalamins, coenzyme F430, and chlorophyll. A particular emphasis is put on the biochemical significance of axial ligands, which can tune the function of the tetrapyrroles. With the use of quantum chemical calculations, it is possible to draw important conclusions regarding aspects of tetrapyrroles that could not otherwise be accessed. The results show that the general reactivity is mainly determined by the metal and the tetrapyrrole ring system, whereas the electronic structure and reactivity are tuned by the choice of axial ligands, providing a unique insight into the design of cofactors in nature.