Electron Transfer Capability in Atomic Hydrogen Reactions for Imidazole Groups Bound to the Insulating Alkanethiolate Layer on Au(111).

The charge transfer capability associated with chemical reactions at metal-organic interfaces was studied via the atomic H addition reaction for an imidazole-terminated alkanethiolate self-assembled monolayer (Im-SAM) film on Au(111) at room temperature, using near-edge X-ray absorption fine structure spectroscopy, infrared reflection absorption spectroscopy, work function measurements, and density functional theory calculations. The imidazolium cation is a stable species in liquids; therefore, it is pertinent to determine whether the hydrogenation reactions of the imidazole groups produce imidazolium cations accompanied by electron transfer to the Au substrate, even in the absence of solvate and/or counterions on the insulating alkanethiolate layer. The experiments made it clear that the imidazolium moieties were formed during the irradiation of Im-SAM with atomic H. Theoretical model calculations also revealed that the total energies and molecular orbital levels satisfied the imidazolium cation formation associated with electron transfer. In a detailed analysis of the work function change depending on H irradiation, we confirmed that some of the imidazolium radicals became cations in Im-SAM.

Phosphorescence Properties of Boron/Bismuth Hybrid Conjugated Materials.

By introducing main-group elements such as boron and bismuth to π-conjugated systems, it is possible to modify the optical properties of π-conjugated materials through orbital interactions between the orbital on the elements and π/π*-orbitals, and the heavy atom effect. Moreover, bismuth, which is the heaviest stable element, induces a significant heavy atom effect, making organobismuth compounds promising for applications as phosphorescent materials. In this study, we synthesized new room-temperature phosphorescent materials by incorporating bismuth into thiophene units. The phosphorescence properties of these materials, such as emission lifetime and wavelength, could be further controlled by combining tricoordinate boron with the thienylbismuth structures. The synthesized bismuth- and boron-containing thiophene compounds exhibited phosphorescence at room temperature in both solution and solid states. Furthermore, the introduction of boron raised the energy of the triplet state in the π-conjugated system, resulting in a blue shift of the phosphorescence wavelength. The analysis of photoluminescence properties and TD-DFT calculations revealed that the introduction of bismuth enhances phosphorescence properties, whereas the introduction of boron further promotes intersystem crossing.

Furanyl bis(indolyl)methane as a palladium ion-selective chromogenic agent.

The colorless solution of furan-2-yl bis(indolyl)methane (BIM) is newly revealed to work as a palladium (Pd2+) ion-selective chromogenic agent by turning orange. 5-(N-Methyl-N-phenyl-aminomethyl)-furan-2-yl BIM could be synthesized from 5-chloromethylfurfural as a biorenewable feedstock via one-pot and double functionalization, and a mixture of its solution and Pd2+ ions showed the highest absorbance at 465 nm in UV-Vis analysis. On the other hand, other metal ions (Cu2+, Cr2+, Cr3+, Fe2+, Fe3+, Ni2+, Zn2+, In2+, Pt2+, or Ce3+) exhibited no response.

Computational Exploration of Minimum Energy Reaction Pathway of N2O Formation from Intermediate I of P450nor Using an Active Center Model.

P450nor is a heme-containing enzyme that catalyzes the conversion of nitric oxide (NO) to nitrous oxide (N2O). Its catalytic mechanism has attracted attention in chemistry, biology, and environmental engineering. The catalytic cycle of P450nor is proposed to consist of three major steps. The reaction mechanism for the last step, N2O generation, remains unknown. In this study, the reaction pathway of the N2O generation from the intermediate I was explored with the B3LYP calculations using an active center model after the examination of the validity of the model. In the validation, we compared the heme distortions between P450nor and other oxidoreductases, suggesting a small effect of protein environment on the N2O generation reaction in P450nor. We then evaluated the electrostatic environment effect of P450nor on the hydride affinity to the active site with quantum mechanics/molecular mechanics (QM/MM) calculations, confirming that the affinity was unchanged with or without the protein environment. The active center model for P450nor showed that the N2O generation process in the enzymatic reaction undergoes a reasonable barrier height without protein environment. Consequently, our findings strongly suggest that the N2O generation reaction from the intermediate I depends sorely on the intrinsic reactivity of the heme cofactor bound on cysteine residue.

Quantum Effects of Hydrogen Nuclei on the Nuclear Magnetic Shielding Tensor of Ice Ih.

Ice is the most fundamental hydrogen-bonded system in which the hydrogen nuclear quantum effect significantly impacts the structure and relevant thermochemical and spectroscopic properties. While ice was experimentally investigated using proton nuclear magnetic resonance spectroscopy more than 40 years ago, the corresponding theoretical investigations have been rarely reported due to the difficulty in evaluating how the proton nuclear quantum effect influences the spectral characteristics of such a condensed material. In this study, we applied a combination of the ONIOM and multicomponent molecular orbital (MC_MO) methods for calculating the anisotropic and isotropic components of the nuclear magnetic shielding tensor of the hexagonal ice crystal to quantify the effects of nuclear quantum fluctuations on the spectroscopic properties of ice. The nuclear magnetic shielding values computed by incorporating the hydrogen nuclear quantum effect reasonably agree with the experimental values. The nuclear quantum effects were found to increase the anisotropic component of the magnetic shielding tensor while decreasing the isotropic component. Such a difference can be explained by their distinct dependence on the electrostatic field and hydrogen-bonding structural parameters.

Prediction of Protein Function from Tertiary Structure of the Active Site in Heme Proteins by Convolutional Neural Network.

Structure-function relationships in proteins have been one of the crucial scientific topics in recent research. Heme proteins have diverse and pivotal biological functions. Therefore, clarifying their structure-function correlation is significant to understand their functional mechanism and is informative for various fields of science. In this study, we constructed convolutional neural network models for predicting protein functions from the tertiary structures of heme-binding sites (active sites) of heme proteins to examine the structure-function correlation. As a result, we succeeded in the classification of oxygen-binding protein (OB), oxidoreductase (OR), proteins with both functions (OB-OR), and electron transport protein (ET) with high accuracy. Although the misclassification rate for OR and ET was high, the rates between OB and ET and between OB and OR were almost zero, indicating that the prediction model works well between protein groups with quite different functions. However, predicting the function of proteins modified with amino acid mutation(s) remains a challenge. Our findings indicate a structure-function correlation in the active site of heme proteins. This study is expected to be applied to the prediction of more detailed protein functions such as catalytic reactions.

Structures and mechanisms of actin ATP hydrolysis.

The major cytoskeleton protein actin undergoes cyclic transitions between the monomeric G-form and the filamentous F-form, which drive organelle transport and cell motility. This mechanical work is driven by the ATPase activity at the catalytic site in the F-form. For deeper understanding of the actin cellular functions, the reaction mechanism must be elucidated. Here, we show that a single actin molecule is trapped in the F-form by fragmin domain-1 binding and present their crystal structures in the ATP analog-, ADP-Pi-, and ADP-bound forms, at 1.15-Å resolutions. The G-to-F conformational transition shifts the side chains of Gln137 and His161, which relocate four water molecules including W1 (attacking water) and W2 (helping water) to facilitate the hydrolysis. By applying quantum mechanics/molecular mechanics calculations to the structures, we have revealed a consistent and comprehensive reaction path of ATP hydrolysis by the F-form actin. The reaction path consists of four steps: 1) W1 and W2 rotations; 2) PG-O3B bond cleavage; 3) four concomitant events: W1-PO3- formation, OH- and proton cleavage, nucleophilic attack by the OH- against PG, and the abstracted proton transfer; and 4) proton relocation that stabilizes the ADP-Pi-bound F-form actin. The mechanism explains the slow rate of ATP hydrolysis by actin and the irreversibility of the hydrolysis reaction. While the catalytic strategy of actin ATP hydrolysis is essentially the same as those of motor proteins like myosin, the process after the hydrolysis is distinct and discussed in terms of Pi release, F-form destabilization, and global conformational changes.

Elucidation of the Correlation between Heme Distortion and Tertiary Structure of the Heme-Binding Pocket Using a Convolutional Neural Network.

Heme proteins serve diverse and pivotal biological functions. Therefore, clarifying the mechanisms of these diverse functions of heme is a crucial scientific topic. Distortion of heme porphyrin is one of the key factors regulating the chemical properties of heme. Here, we constructed convolutional neural network models for predicting heme distortion from the tertiary structure of the heme-binding pocket to examine their correlation. For saddling, ruffling, doming, and waving distortions, the experimental structure and predicted values were closely correlated. Furthermore, we assessed the correlation between the cavity shape and molecular structure of heme and demonstrated that hemes in protein pockets with similar structures exhibit near-identical structures, indicating the regulation of heme distortion through the protein environment. These findings indicate that the tertiary structure of the heme-binding pocket is one of the factors regulating the distortion of heme porphyrin, thereby controlling the chemical properties of heme relevant to the protein function; this implies a structure-function correlation in heme proteins.

Global Analysis of Heme Proteins Elucidates the Correlation between Heme Distortion and the Heme-Binding Pocket.

Heme proteins play diverse and important biological roles, from electron transfer and chemical catalysis to oxygen transport and/or storage. Although the distortion of heme porphyrin correlates with the physical properties of heme, such as the redox potential and oxygen affinity, the relationship between heme distortion and the heme protein environment is unclear. Here, we tested the hypothesis that the protein environment of the heme-binding pocket determines heme distortion (conformation). We analyzed the correlations between the amino acid composition of the heme-binding pocket and the magnitude of heme distortion along 12 vibrational modes using machine learning. A correlation was detected in the three lowest vibrational modes. Analysis of heme distortions in nearly the same environments of the heme-binding pocket supported this notion. Our analyses indicate that the heme-binding pocket environment is a major factor impacting the distortion of heme porphyrin along the three lowest vibrational modes. In addition, statistical analysis of the distortion of heme porphyrin revealed that the peaks of distributions of the ruffling and breathing distortions are shifted from 0 (the equilibrium structure). Both the ruffling and breathing distortions are correlated with the redox potential of heme, so that heme molecules with these distortions have a lower redox potential than planar molecules. These findings explain the structure-function relationship of heme.

Hydrogen/Deuterium Transfer from Anisole to Methoxy Radicals: A Theoretical Study of a Deuterium-Labeled Drug Model.

Recently, deuterium-labeled drugs, such as deutetrabenazine, have attracted considerable attention. Consequently, understanding the reaction mechanisms of deuterium-labeled drugs is crucial, both fundamentally and for real applications. To understand the mechanisms of H- and D-transfer reactions, in this study, we used deuterated anisole as a deutetrabenazine model and computationally considered the nuclear quantum effects of protons, deuterons, and electrons. We demonstrated that geometrical differences exist in the partially and fully deuterated methoxy groups and hydrogen-bonded structures of intermediates and transition states due to the H/D isotope effect. The observed geometrical features and electronic structures are ascribable to the different nuclear quantum effects of protons and deuterons. Primary and secondary kinetic isotope effects (KIEs) were calculated for H- and D-transfer reactions from deuterated and undeuterated anisole, with the calculated primary KIEs in good agreement with the corresponding experimental data. These results reveal that the nuclear quantum effects of protons and deuterons need to be considered when analyzing the reaction mechanisms of H- and D-transfer reactions and that a theoretical approach that directly includes nuclear quantum effects is a powerful tool for the analysis of H/D isotope effects in H- and D-transfer reactions.

Short-lived intermediate in N2O generation by P450 NO reductase captured by time-resolved IR spectroscopy and XFEL crystallography.

Nitric oxide (NO) reductase from the fungus Fusarium oxysporum is a P450-type enzyme (P450nor) that catalyzes the reduction of NO to nitrous oxide (N2O) in the global nitrogen cycle. In this enzymatic reaction, the heme-bound NO is activated by the direct hydride transfer from NADH to generate a short-lived intermediate ( I ), a key state to promote N-N bond formation and N-O bond cleavage. This study applied time-resolved (TR) techniques in conjunction with photolabile-caged NO to gain direct experimental results for the characterization of the coordination and electronic structures of I TR freeze-trap crystallography using an X-ray free electron laser (XFEL) reveals highly bent Fe-NO coordination in I , with an elongated Fe-NO bond length (Fe-NO = 1.91 Å, Fe-N-O = 138°) in the absence of NAD+ TR-infrared (IR) spectroscopy detects the formation of I with an N-O stretching frequency of 1,290 cm-1 upon hydride transfer from NADH to the Fe3+-NO enzyme via the dissociation of NAD+ from a transient state, with an N-O stretching of 1,330 cm-1 and a lifetime of ca. 16 ms. Quantum mechanics/molecular mechanics calculations, based on these crystallographic and IR spectroscopic results, demonstrate that the electronic structure of I is characterized by a singly protonated Fe3+-NHO•- radical. The current findings provide conclusive evidence for the N2O generation mechanism via a radical-radical coupling of the heme nitroxyl complex with the second NO molecule.

PyDISH: database and analysis tools for heme porphyrin distortion in heme proteins.

Heme participates in a wide range of biological functions such as oxygen transport, electron transport, oxygen reduction, transcriptional regulation and so on. While the mechanism of each function has been investigated for many heme proteins, the origin of the diversity of the heme functions is still unclear and a crucial scientific issue. We have constructed a database of heme proteins, named Python-based database and analyzer for DIStortion of Heme porphyrin (PyDISH), which also contains some analysis tools. The aim of PyDISH is to integrate the information on the structures of hemes and heme proteins and the functions of heme proteins. This database will provide the structure-function relationships focusing on heme porphyrin distortion and lead to the elucidation of the origin of the functional diversity of heme proteins. In addition, the insights obtained from the database can be used for the design of protein function. PyDISH contains the structural data of more than 13 000 hemes extracted from the Protein Data Bank, including heme porphyrin distortion, axial ligands coordinating to the heme and the orientation of the propionate sidechains of heme. PyDISH also has information about the protein domains, including Uniprot ID, protein fold by CATH ID, organism, coordination distance and so on. The analytical tools implemented in PyDISH allow users to not only browse and download the data but also analyze the structures of heme porphyrin by using the analytical tools implemented in PyDISH. PyDISH users will be able to utilize the obtained results for the design of protein function. Database URL: http://pydish.bio.info.hiroshima-cu.ac.jp/.

Strong Hydrogen Bonds at the Interface between Proton-Donating and -Accepting Self-Assembled Monolayers on Au(111).

Hydrogen-bonding heterogeneous bilayers on substrates have been studied as a base for new functions of molecular adlayers by means of atomic force microscopy (AFM), X-ray photoelectron spectroscopy (XPS), infrared reflection absorption spectroscopy (IRAS), and density functional theory (DFT) calculations. Here, we report the formation of the catechol-fused bis(methylthio)tetrathiafulvalene (H2Cat-BMT-TTF) adlayer hydrogen bonding with an imidazole-terminated alkanethiolate self-assembled monolayer (Im-SAM) on Au(111). The heterogeneous bilayer is realized by sequential two-step immersions in solutions for the individual Im-SAM and H2Cat-BMT-TTF adlayer formations. In the measurements by AFM, a grained H2Cat-BMT-TTF adlayer on Im-SAM is revealed. The coverage and the chemical states of H2Cat-BMT-TTF on Im-SAM are specified by XPS. On the vibrational spectrum measured by IRAS, the strong hydrogen bonds between H2Cat-BMT-TTF and Im-SAM are characterized by the remarkably red-shifted OH stretching mode at 3140 cm-1, which is much lower than that for hydrogen-bonding water (typically ∼3300 cm-1). The OH stretching mode frequency and the adsorption strength for the H2Cat-BMT-TTF molecule hydrogen bonding with imidazole groups are quantitatively examined on the basis of DFT calculations.

Capturing an initial intermediate during the P450nor enzymatic reaction using time-resolved XFEL crystallography and caged-substrate.

Time-resolved serial femtosecond crystallography using an X-ray free electron laser (XFEL) in conjunction with a photosensitive caged-compound offers a crystallographic method to track enzymatic reactions. Here we demonstrate the application of this method using fungal NO reductase, a heme-containing enzyme, at room temperature. Twenty milliseconds after caged-NO photolysis, we identify a NO-bound form of the enzyme, which is an initial intermediate with a slightly bent Fe-N-O coordination geometry at a resolution of 2.1 Å. The NO geometry is compatible with those analyzed by XFEL-based cryo-crystallography and QM/MM calculations, indicating that we obtain an intact Fe3+-NO coordination structure that is free of X-ray radiation damage. The slightly bent NO geometry is appropriate to prevent immediate NO dissociation and thus accept H- from NADH. The combination of using XFEL and a caged-compound is a powerful tool for determining functional enzyme structures during catalytic reactions at the atomic level.

Theoretical study of the H/D isotope effect on phase transition of hydrogen-bonded organic conductor κ-H3(Cat-EDT-TTF)2.

κ-H3(Cat-EDT-TTF)2 (H-TTF) is a hydrogen-bonded π-electron system which was found to reveal C2/c symmetry at 50-293 K, while its isotopologue, κ-D3(Cat-EDT-TTF)2 (D-TTF), showed the phase transition at 185 K from C2/c to P1[combining macron]. To elucidate the origin of such a difference, we calculated the potential energy curves (PECs) for the hydrogen transfer along the H-bonds in these conductors. We found that both the π-stacking and the hydrogen nuclear quantum effect drastically affected the hydrogen transfer energy. By taking account of both effects, we obtained a symmetric single-well effective PEC for H-TTF, which indicated that the hydrogen was always located at the center of the H-bond. By contrast, the effective PEC of D-TTF was a low-barrier double-well, indicating that the position of the H-bonded deuterium would change according to the temperature. We concluded that the π-stacking and the nuclear quantum effect were the key factors for the appearance of phase transition only in D-TTF.

Quantum Mechanics/Molecular Mechanics Study of the Sialyltransferase Reaction Mechanism.

The sialyltransferase is an enzyme that transfers the sialic acid moiety from cytidine 5'-monophospho-N-acetyl-neuraminic acid (CMP-NeuAc) to the terminal position of glycans. To elucidate the catalytic mechanism of sialyltransferase, we explored the potential energy surface along the sialic acid transfer reaction coordinates by the hybrid quantum mechanics/molecular mechanics method on the basis of the crystal structure of sialyltransferase CstII. Our calculation demonstrated that CstII employed an SN1-like reaction mechanism via the formation of a short-lived oxocarbenium ion intermediate. The computational barrier height was 19.5 kcal/mol, which reasonably corresponded with the experimental reaction rate. We also found that two tyrosine residues (Tyr156 and Tyr162) played a vital role in stabilizing the intermediate and the transition states by quantum mechanical interaction with CMP.

Inverse Ubbelohde effect in the short hydrogen bond of photosystem II: Relation between H/D isotope effect and symmetry in potential energy profile.

The short hydrogen bond between tyrosine Yz and D1-His190 of photosystem II (PSII) was investigated using multicomponent quantum mechanics, where the quantum fluctuation of a hydrogen nucleus was incorporated into electronic structure calculation. Our computation demonstrated that the deuteration for hydrogen in the short hydrogen bond of PSII led to the reduction of the O…N distance. It indicated an inverse Ubbelohde effect typically recognized in strong and symmetric hydrogen-bonding clusters such as FHF(-) and H3O2-. We confirmed that the relation between the geometric isotope effect and the symmetry of the potential energy profile of FHF(-) was reasonably agreed with that of PSII. According to this agreement, the short hydrogen bond in PSII can be regarded as a short strong hydrogen bond. © 2016 Wiley Periodicals, Inc.

Vibrational analysis on the revised potential energy curve of the low-barrier hydrogen bond in photoactive yellow protein.

Photoactive yellow protein (PYP) has a characteristic hydrogen bond (H bond) between p-coumaric acid chromophore and Glu46, whose OH bond length has been observed to be 1.21 Å by the neutron diffraction technique [Proc. Natl. Acad. Sci. 106, 440-4]. Although it has been expected that such a drastic elongation of the OH bond could be caused by the quantum effect of the hydrogen nucleus, previous theoretical computations including the nuclear quantum effect have so far underestimated the bond length by more than 0.07 Å. To elucidate the origin of the difference, we performed a vibrational analysis of the H bond on potential energy curve with O…O distance of 2.47 Å on the equilibrium structure, and that with O…O distance of 2.56 Å on the experimental crystal structure. While the vibrationally averaged OH bond length for equilibrium structure was underestimated, the corresponding value for crystal structure was in reasonable agreement with the corresponding experimental values. The elongation of the O…O distance by the quantum mechanical or thermal fluctuation would be indispensable for the formation of a low-barrier hydrogen bond in PYP.

Isotope effect on the circular dichroism spectrum of methyl α-D-glucopyranoside in aqueous solution.

H/D isotope effect on the circular dichroism spectrum of methyl α-D-glucopyranoside in aqueous solution has been analyzed by multicomponent density functional theory calculations using the polarizable continuum model. By comparing the computational spectra with the corresponding experimental spectrum obtained with a vacuum-ultraviolet circular dichroism spectrophotometer, it was demonstrated that the isotope effect provides insights not only into the isotopic difference of the intramolecular interactions of the solutes, but also into that of the solute-solvent intermolecular interaction.

Performance test of multicomponent quantum mechanical calculation with polarizable continuum model for proton chemical shift.

Multicomponent quantum mechanical (MC_QM) calculations with polarizable continuum model (PCM) have been tested against liquid (1)H NMR chemical shifts for a test set of 80 molecules. Improvement from conventional quantum mechanical calculations was achieved for MC_QM calculations. The advantage of the multicomponent scheme could be attributed to the geometrical change from the equilibrium geometry by the incorporation of the hydrogen nuclear quantum effect, while that of PCM can be attributed to the change of the electronic structure according to the polarization by solvent effects.