Neutron crystallography and quantum chemical analysis of bilin reductase PcyA mutants reveal substrate and catalytic residue protonation states.

PcyA, a ferredoxin-dependent bilin pigment reductase, catalyzes the site-specific reduction of the two vinyl groups of biliverdin (BV), producing phycocyanobilin. Previous neutron crystallography detected both the neutral BV and its protonated form (BVH+) in the wildtype (WT) PcyA-BV complex, and a nearby catalytic residue Asp105 was found to have two conformations (protonated and deprotonated). Semiempirical calculations have suggested that the protonation states of BV are reflected in the absorption spectrum of the WT PcyA-BV complex. In the previously determined absorption spectra of the PcyA D105N and I86D mutants, complexed with BV, a peak at 730 nm, observed in the WT, disappeared and increased, respectively. Here, we performed neutron crystallography and quantum chemical analysis of the D105N-BV and I86D-BV complexes to determine the protonation states of BV and the surrounding residues and study the correlation between the absorption spectra and protonation states around BV. Neutron structures elucidated that BV in the D105N mutant is in a neutral state, whereas that in the I86D mutant is dominantly in a protonated state. Glu76 and His88 showed different hydrogen bonding with surrounding residues compared with WT PcyA, further explaining why D105N and I86D have much lower activities for phycocyanobilin synthesis than the WT PcyA. Our quantum mechanics/molecular mechanics calculations of the absorption spectra showed that the spectral change in D105N arises from Glu76 deprotonation, consistent with the neutron structure. Collectively, our findings reveal more mechanistic details of bilin pigment biosynthesis.

Difference in reaction mechanism between ZnZrOx and InZrOx for CO2 hydrogenation.

Oxide solid-solution catalysts, such as Zn-doped ZrO2 (ZnZrOx) and In-doped ZrO2 (InZrOx), exhibit distinctive catalytic capabilities for CH3OH synthesis via CO2 hydrogenation. We investigated the active site structures of these catalysts and their associated reaction mechanisms using both experimental and computational approaches. Electron microscopy and X-ray absorption spectroscopy reveal that the primary active sites are isolated cations, such as Zn2+ and In3+, dissolved in tetragonal ZrO2. Notably, for Zn2+, decomposition of the methoxy group, which is an essential intermediate in CH4 synthesis, is partially suppressed because of the relatively high stability of the methoxy group. Conversely, the methyl group strongly adsorbs on In3+, facilitating the conversion of the methoxy species into methyl groups. The decomposition of CH3OH is also suggested to contribute to CH4 synthesis. These results highlight the generation of CH4 as a byproduct of the InZrOx catalyst. Understanding the active site structure and elucidating the reaction mechanism at the atomic level are anticipated to contribute significantly to the future development of oxide solid-solution catalysts.

Effects of third-order susceptibility in sum frequency generation spectra: a molecular dynamics study in liquid water.

When sum frequency generation (SFG) spectroscopy is applied to charged solid/liquid interfaces, the observed SFG signals include both the second-order and third-order polarizations. The latter is called the χ(3) effect, which mainly includes induced molecular orientation by electric fields at charged interfaces. We theoretically evaluate the χ(3) effect on the SFG spectroscopy of liquid water using molecular dynamics (MD) simulations. The MD simulations enable us to definitely calculate the χ(3) susceptibility as a bulk property, and thereby separating it from the usual χ(2) effect shown in the SFG spectra. The calculated results of χ(3) for liquid water are fairly consistent with the experimental estimates. The present finding is utilized to analyze the spectral change of SFG at the air/water interface under electric fields and at the charged silica/water interface. The present analysis of the spectral changes allows for distinguishing the intrinsic change of the interface structure and the χ(3) effect from bulk liquid.

Theoretical and experimental examination of SFG polarization analysis at acetonitrile-water solution surfaces.

Sum frequency generation (SFG) spectroscopy is widely used to observe molecular orientation at interfaces through a combination of various types of polarization. The present work thoroughly examines the relation between the polarization dependence of SFG signals and the molecular orientation, by comparing SFG measurements and molecular dynamics (MD) simulations of acetonitrile/water solutions. The present SFG experiment and MD simulations yield quite consistent results on the ratios of χ(2) elements, supporting the reliability of both means. However, the subsequent polarization analysis tends to derive more upright tilt angles of acetonitrile than the direct MD calculations. The reasons for discrepancy are examined in terms of three issues; (i) anisotropy of the Raman tensor, (ii) cross-correlation, and (iii) orientational distribution. The analysis revealed that the issues (i) and (iii) are the main causes of errors in the conventional polarization analysis of SFG spectra. In methyl CH stretching, the anisotropy of Raman tensor cannot be estimated from the simple bond polarizability model. The neglect of the orientational distribution is shown to systematically underestimate the tilt angle of acetonitrile. Further refined use of polarization analysis in collaboration with MD simulations should be proposed.

Molecular Mechanism of Autodissociation in Liquid Water: Ab Initio Molecular Dynamics Simulations.

Autodissociation in liquid water is one of the most important processes in various topics of physical chemistry, such as acid-base chemistry. Molecular simulations have elucidated most of the molecular mechanisms at the atomic level, yet quantitative analysis to compare with experiments using the potential of mean force (PMF) remains a hurdle, including the definition of reaction coordinates and the accuracy of liquid structures by ab initio molecular dynamics (AIMD) simulations with density functional theory (DFT) methods. Here, we perform AIMD simulations with the revPBE-D3 exchange-correlation functional to compute the PMF profiles of autoionization, or proton transfer (PT), in liquid water. For the quantitative analysis with physically meaningful reaction coordinates, we employ a PT coordinate, donor-acceptor (OH--H3O+) distance, and hydrogen (H)-bond number. The one-dimensional (1D) PMF profile along the PT coordinate shows no local minimum in the product state of PT (OH- and H3O+), which is necessary to accurately compute the acid dissociation constant (or pKa). On the other hand, the 2D PMF profiles along the PT coordinate and donor-acceptor distance show local minima in the product state and reaction barriers, and the computed pKw is comparable to the experiment. In addition, the 2D PMF profiles along the PT coordinate and the H-bond number reveal the molecular mechanism of the H-bond rearrangement concomitant with PT, in which the H-bond breaking before PT is slightly preferable. These findings indicate that an accurate evaluation of pKa by MD simulations requires the donor-acceptor distance in addition to the conventional PT coordinate.

Dynamics of proton transfer and vibrational relaxation in dilute hydrofluoric acid.

The molecular mechanisms in both vibrational relaxation and proton transfer (PT) associated with infrared (IR)-induced PT in a dilute hydrofluoric acid solution at ambient temperature are studied by molecular dynamics (MD) simulations with the multistate empirical valence bond model. To investigate the solvation dynamics, a collective solvent coordinate and its perpendicular bath modes are defined from the diabatic energy gap and their motions are examined by the generalized Langevin equation (GLE) formalism. The GLE analysis using the equilibrium MD simulation shows that the major solvent reorganizations in the PT are represented by the libration and hindered translation. In particular, the libration gives the stronger coupling to the solvent reorganization and the faster relaxation. The nonequilibrium MD simulation demonstrated that both the HF stretching vibration and the solvent reorganization relax on a similar time scale and thus compete in the PT. It also supported the "presolvation mechanism" for the PT in this system.

Hydration structure in dilute hydrofluoric acid.

We have performed the multistate empirical valence bond (MS-EVB) molecular dynamics simulations of a dilute hydrofluoric acid solution at ambient temperature to study the hydration structure associated with its weak acidity. The developed MS-EVB model showed reasonable agreement with experimental and previous ab initio molecular dynamics and reference interaction site model self-consistent field simulations for the free energy and structural properties. The local tetrahedral and translational order parameters around the fluorine atom significantly increase in the transition and product states of the HF dissociation reaction. This indicates that the angular and translational rearrangements of the hydrogen-bond topology are necessary especially around the fluorine atom. At the transition state of the proton transfer, the tetrahedral order parameters are very large, whereas the translational order parameters are not. This suggests that for the proton transfer to occur the large angular rearrangements of the hydrogen-bond topology are more necessary than the translational ones.

Vibrational spectroscopy and relaxation of an anharmonic oscillator coupled to harmonic bath.

The vibrational spectroscopy and relaxation of an anharmonic oscillator coupled to a harmonic bath are examined to assess the applicability of the time correlation function (TCF), the response function, and the semiclassical frequency modulation (SFM) model to the calculation of infrared (IR) spectra. These three approaches are often used in connection with the molecular dynamics simulations but have not been compared in detail. We also analyze the vibrational energy relaxation (VER), which determines the line shape and is itself a pivotal process in energy transport. The IR spectra and VER are calculated using the generalized Langevin equation (GLE), the Gaussian wavepacket (GWP) method, and the quantum master equation (QME). By calculating the vibrational frequency TCF, a detailed analysis of the frequency fluctuation and correlation time of the model is provided. The peak amplitude and width in the IR spectra calculated by the GLE with the harmonic quantum correction are shown to agree well with those by the QME though the vibrational frequency is generally overestimated. The GWP method improves the peak position by considering the zero-point energy and the anharmonicity although the red-shift slightly overshoots the QME reference. The GWP also yields an extra peak in the higher-frequency region than the fundamental transition arising from the difference frequency of the center and width oscillations of a wavepacket. The SFM approach underestimates the peak amplitude of the IR spectra but well reproduces the peak width. Further, the dependence of the VER rate on the strength of an excitation pulse is discussed.

Constrained Density Functional Theory Molecular Dynamics Simulation of Deprotonation in Aqueous Silicic Acid.

The solubility of silicic acid in water is so low that the molecular mechanism behind the physical properties such as pKa remains poorly understood, despite the importance in fields such as chemistry and geology. Theoretical calculations provide the molecular-level information on such a rare chemical species yet face difficulties in the selection of reaction coordinates and the rare-event sampling of proton transfers (PTs) in condensed phases. We propose a straightforward calculation scheme of pKa with ab initio molecular dynamics (MD) simulation and the constrained density functional theory (CDFT), which provides structural and dynamical properties such as radial distribution functions, vibrational spectra, and reaction paths. The calculated reaction free energies of deprotonations agreed with experiments within a few kcal/mol. Analysis of the solvation structure shows that, after deprotonation, the hydronium ion D3O+ repels away from the deprotonated silicic acid SiO(OD)3- without forming a stable contact-ion pair. The calculated vibrational spectra are consistent with the spectroscopic measurements, and the dynamical analysis of the reaction path quantifies the couplings of the OD stretch vibrations of silicic acid and water to the PT reaction in terms of the vertical energy-gap coordinate that includes both the solute and the solvent components.

Efficient Computation of Difference Vibrational Spectra in Isothermal-Isobaric Ensemble.

Difference spectroscopy between two close systems is widely used to augment its selectivity to the different parts of the observed system, though the molecular dynamics calculation of tiny difference spectra would be computationally extraordinary demanding by subtraction of two spectra. Therefore, we have proposed an efficient computational algorithm of difference spectra without resorting to the subtraction. The present paper reports our extension of the theoretical method in the isothermal-isobaric (NPT) ensemble. The present theory expands our applications of analysis including pressure dependence of the spectra. We verified that the present theory yields accurate difference spectra in the NPT condition as well, with remarkable computational efficiency over the straightforward subtraction by several orders of magnitude. This method is further applied to vibrational spectra of liquid water with varying pressure and succeeded in reproducing tiny difference spectra by pressure change. The anomalous pressure dependence is elucidated in relation to other properties of liquid water.

Improved Theory of Difference Vibrational Spectroscopy and Application to Water.

Measuring difference spectrum of two closely related systems is a versatile experimental method to highlight the difference in the two systems. However, the computation of a minuscule difference between two spectra by molecular dynamics (MD) simulation is far more challenging than that of each spectrum in terms of statistical convergence. Therefore, we have proposed a theory of difference spectra, which accelerates the calculation of difference spectra by several orders of magnitude [ Sakaguchi, S. et al. J. Chem. Phys. 2014 , 140 , 144109 ]. The present paper reports our subsequent advances in the computational method to greatly improve its accuracy, numerical stability, and applicability. The present method of computation based on the nonequilibrium MD simulation allows for general molecular models including polarizable ones without sacrificing its computational efficiency. The improved method was applied to polarizable liquid water and yielded difference spectra of infrared, Raman, and vibrational sum frequency generation spectroscopies. The present method enables us to analyze various difference spectra of large molecular systems using MD simulation.

Vibrational Quantum Decoherence in Liquid Water.

Traditional descriptions of vibrational energy transfer consider a quantum oscillator interacting with a classical environment. However, a major limitation of this simplified description is the neglect of quantum decoherence induced by the different interactions between two distinct quantum states and their environment, which can strongly affect the predicted energy-transfer rate and vibrational spectra. Here, we use quantum-classical molecular dynamics simulations to determine the vibrational quantum decoherence time for an OH stretch vibration in liquid heavy water. We show that coherence is lost on a sub-100 fs time scale due to the different responses of the first shell neighbors to the ground and excited OH vibrational states. This ultrafast decoherence induces a strong homogeneous contribution to the linear infrared spectrum and suggests that resonant vibrational energy transfer in H2O may be more incoherent than previously thought.