A machine learning potential construction based on radial distribution function sampling.

Sampling reference data is crucial in machine learning potential (MLP) construction. Inadequate coverage of local configurations in reference data may lead to unphysical behaviors in MLP-based molecular dynamics (MLP-MD) simulations. To address this problem, this study proposes a new on-the-fly reference data sampling method called radial distribution function (RDF)-based data sampling for MLP construction. This method detects and extracts anomalous structures from the trajectories of MLP-MD simulations by focusing on the shapes of RDFs. The detected structures are added to the reference data to improve the accuracy of the MLP. This method allows us to realize a reasonable MLP construction for liquid water with minimal additional data. We prepare data from an H2O molecular cluster system and verify whether the constructed MLPs are practical for bulk water systems. MLP-MD simulations without RDF-based data sampling show unphysical behaviors, such as atomic collisions. In contrast, after applying this method, we obtain MLP-MD trajectories with features, such as RDF shapes and angle distributions, that are comparable to those of ab initio MD simulations. Our simulation results demonstrate that the RDF-based data sampling approach is useful for constructing MLPs that are robust to extrapolations from molecular cluster systems to bulk systems without any specialized know-how.

A new small molecule DoNA binding to CAG repeat RNA.

We here report a new molecule DoNA binding to a CAG repeat RNA. DoNA is a dimer of the NA molecule that we previously reported. NA binds with high affinity to a CAG repeat DNA but not significantly to a CAG repeat RNA. Binding analyses using SPR and CSI-TOF MS indicated a significant increase in the affinity of DoNA to a single stranded CAG repeat RNA compared to NA. Systematic investigation of the RNA motifs bound by DoNA using hairpin RNA models revealed that DoNA binds to the CAG units at overhang and terminal positions, and notably, it binds to the structurally flexible internal and hairpin loop region.

Theoretical elucidation of the structure, bonding, and reactivity of the CaMn4Ox clusters in the whole Kok cycle for water oxidation embedded in the oxygen evolving center of photosystem II. New molecular and quantum insights into the mechanism of the O-O bond formation.

This paper reviews our historical developments of broken-symmetry (BS) and beyond BS methods that are applicable for theoretical investigations of metalloenzymes such as OEC in PSII. The BS hybrid DFT (HDFT) calculations starting from high-resolution (HR) XRD structure in the most stable S1 state have been performed to elucidate structure and bonding of whole possible intermediates of the CaMn4Ox cluster (1) in the Si (i = 0 ~ 4) states of the Kok cycle. The large-scale HDFT/MM computations starting from HR XRD have been performed to elucidate biomolecular system structures which are crucial for examination of possible water inlet and proton release pathways for water oxidation in OEC of PSII. DLPNO CCSD(T0) computations have been performed for elucidation of scope and reliability of relative energies among the intermediates by HDFT. These computations combined with EXAFS, XRD, XFEL, and EPR experimental results have elucidated the structure, bonding, and reactivity of the key intermediates, which are indispensable for understanding and explanation of the mechanism of water oxidation in OEC of PSII. Interplay between theory and experiments have elucidated important roles of four degrees of freedom, spin, charge, orbital, and nuclear motion for understanding and explanation of the chemical reactivity of 1 embedded in protein matrix, indicating the participations of the Ca(H2O)n ion and tyrosine(Yz)-O radical as a one-electron acceptor for the O-O bond formation. The Ca-assisted Yz-coupled O-O bond formation mechanisms for water oxidation are consistent with recent XES and very recent time-resolved SFX XFEL and FTIR results.

Specific zinc binding to heliorhodopsin.

Heliorhodopsins (HeRs), a recently discovered family of rhodopsins, have an inverted membrane topology compared to animal and microbial rhodopsins. The slow photocycle of HeRs suggests a light-sensor function, although the actual function remains unknown. Although HeRs exhibit no specific binding of monovalent cations or anions, recent ATR-FTIR spectroscopy studies have demonstrated the binding of Zn2+ to HeR from Thermoplasmatales archaeon (TaHeR) and 48C12. Even though ion-specific FTIR spectra were observed for many divalent cations, only helical structural perturbations were observed for Zn2+-binding, suggesting a possible modification of the HeR function by Zn2+. The present study shows that Zn2+-binding lowers the thermal stability of TaHeR, and slows back proton transfer to the retinal Schiff base (M decay) during its photocycle. Zn2+-binding was similarly observed for a TaHeR opsin that lacks the retinal chromophore. We then studied the Zn2+-binding site by means of the ATR-FTIR spectroscopy of site-directed mutants. Among five and four mutants of His and Asp/Glu, respectively, only E150Q exhibited a completely different spectral feature of the α-helix (amide-I) in ATR-FTIR spectroscopy, suggesting that E150 is responsible for Zn2+-binding. Molecular dynamics (MD) simulations built a coordination structure of Zn2+-bound TaHeR, where E150 and protein bound water molecules participate in direct coordination. It was concluded that the specific binding site of Zn2+ is located at the cytoplasmic side of TaHeR, and that Zn2+-binding affects the structure and structural dynamics, possibly modifying the unknown function of TaHeR.

Enantioselective amino acid interactions in solution.

Enantiomeric excesses (ee) of L-amino acids in meteorites are higher than 10%, especially for isovaline (Iva). This suggests the existence of some kind of triggering mechanism responsible for the amplification of the ee from an initial small value. Here, we investigate the dimeric molecular interactions of alanine (Ala) and Iva in solution as an initial nucleation step of crystals at an accurate first-principles level. We find that the dimeric interaction of Iva is more chirality-dependent than that of Ala, thus providing a clear molecular-level insight into the enantioselectivity of amino acids in solution.

Organocatalytic-racemization reaction elucidation of aspartic acid by density functional theory.

Aldehydes and carboxylic acids are widely used as catalysts for efficient racemization process of amino acids. However, the detailed reaction mechanism remains unclear. This work aims to clarify the racemization mechanism of aspartic acid (Asp) catalyzed by salicylaldehyde and acetic acid by using computational approaches. Density functional theory was used to obtain the structures and relative energies of 10 intermediates and five transition states, thus characterizing the main stages of the reaction. The calculated energy diagram shows that the dehydration step has the highest energy barrier, followed by the reaction step to change the chirality of Asp, which is a crucial process for racemization. In the dehydration reaction, water molecules can induce a remarkable decrease in the required energy.

Neutron crystallography of copper amine oxidase reveals keto/enolate interconversion of the quinone cofactor and unusual proton sharing.

Recent advances in neutron crystallographic studies have provided structural bases for quantum behaviors of protons observed in enzymatic reactions. Thus, we resolved the neutron crystal structure of a bacterial copper (Cu) amine oxidase (CAO), which contains a prosthetic Cu ion and a protein-derived redox cofactor, topa quinone (TPQ). We solved hitherto unknown structures of the active site, including a keto/enolate equilibrium of the cofactor with a nonplanar quinone ring, unusual proton sharing between the cofactor and the catalytic base, and metal-induced deprotonation of a histidine residue that coordinates to the Cu. Our findings show a refined active-site structure that gives detailed information on the protonation state of dissociable groups, such as the quinone cofactor, which are critical for catalytic reactions.

The Nature of the Chemical Bonds of High-Valent Transition-Metal Oxo (M=O) and Peroxo (MOO) Compounds: A Historical Perspective of the Metal Oxyl-Radical Character by the Classical to Quantum Computations.

This review article describes a historical perspective of elucidation of the nature of the chemical bonds of the high-valent transition metal oxo (M=O) and peroxo (M-O-O) compounds in chemistry and biology. The basic concepts and theoretical backgrounds of the broken-symmetry (BS) method are revisited to explain orbital symmetry conservation and orbital symmetry breaking for the theoretical characterization of four different mechanisms of chemical reactions. Beyond BS methods using the natural orbitals (UNO) of the BS solutions, such as UNO CI (CC), are also revisited for the elucidation of the scope and applicability of the BS methods. Several chemical indices have been derived as the conceptual bridges between the BS and beyond BS methods. The BS molecular orbital models have been employed to explain the metal oxyl-radical character of the M=O and M-O-O bonds, which respond to their radical reactivity. The isolobal and isospin analogy between carbonyl oxide R2C-O-O and metal peroxide LFe-O-O has been applied to understand and explain the chameleonic chemical reactivity of these compounds. The isolobal and isospin analogy among Fe=O, O=O, and O have also provided the triplet atomic oxygen (3O) model for non-heme Fe(IV)=O species with strong radical reactivity. The chameleonic reactivity of the compounds I (Cpd I) and II (Cpd II) is also explained by this analogy. The early proposals obtained by these theoretical models have been examined based on recent computational results by hybrid DFT (UHDFT), DLPNO CCSD(T0), CASPT2, and UNO CI (CC) methods and quantum computing (QC).

Heme-bound tyrosine vibrations in hemoglobin M: Resonance Raman, crystallography, and DFT calculation.

Hemoglobins M (Hbs M) are human hemoglobin variants in which either the α or ß subunit contains a ferric heme in the α2ß2 tetramer. Though the ferric subunit cannot bind O2, it regulates O2 affinity of its counterpart ferrous subunit. We have investigated resonance Raman spectra of two Hbs, M Iwate (α87His → tyrosine [Tyr]) and M Boston (α58His → Tyr), having tyrosine as a heme axial ligand at proximal and distal positions, respectively, that exhibit unassigned resonance Raman bands arising from ferric (not ferrous) hemes at 899 and 876 cm-1. Our quantum chemical calculations using density functional theory on Fe-porphyrin models with p-cresol and/or 4-methylimidazole showed that the unassigned bands correspond to the breathing-like modes of Fe3+-bound Tyr and are sensitive to the Fe-O-C(Tyr) angle. Based on the frequencies of the Raman bands, the Fe-O-C(Tyr) angles of Hbs M Iwate and M Boston were predicted to be 153.5° and 129.2°, respectively. Consistent with this prediction, x-ray crystallographic analysis showed that the Fe-O-C(Tyr) angles of Hbs M Iwate and M Boston in the T quaternary structure were 153.6° and 134.6°, respectively. It also showed a similar Fe-O bond length (1.96 and 1.97 Å) and different tilting angles.

Reaction mechanism of N-cyclopropylglycine oxidation by monomeric sarcosine oxidase.

Monomeric sarcosine oxidase (MSOX) is a fundamental - yet one of the simplest - member of a family of flavoenzymes able to catalyze the oxidation of sarcosine (N-methylglycine) and other secondary amines. MSOX is one of the best characterized members of the amine oxidoreductases (AOs), however, its reaction mechanism is still controversial. A single electron transfer (SET) process was suggested on the basis of studies with N-cyclopropylglycine (CPG), although a hydride transfer mechanism would be more consistent in general for AOs. To shed some light on the detailed reaction mechanisms of CPG in MSOX, we performed hybrid quantum mechanical/molecular mechanical (QM/MM) simulations. We found that the polar mechanism is energetically the most favorable. The free energy profile indicates that the first rate-limiting step is the CPG binding to the flavin ring which simultaneously proceeds with the ring-opening of the CPG cyclopropyl group. This reaction step of the CPG adduct formation corresponds to the nucleophilic attack of the cyclopropyl group (C3 atom) to the flavin ring (C4a atom), whereas the expected radical species formation in the SET mechanism was not observed. The following inactivated species, which accumulates during the CPG oxidation in MSOX, can be ascribed to an imine state, and not an enamine state, on the basis of the computed UV/Vis spectra. The conformation of CPG was found to be crucial for reactions following the CPG adduct formation.

UNO DMRG CAS CI calculations of binuclear manganese complex Mn(IV)2 O2 (NHCHCO2 )4 : Scope and applicability of Heisenberg model.

Both direct exchange and super-exchange interactions cooperate to realize inter-spin magnetic interaction in binuclear manganese complex Mn(IV)2 O2 (NHCHCO2 )4 with a di-µ-oxo path. We revisited this spin system using DMRG CAS methods and CAS selection procedures. Our results indicate that our previous "dynamically extended spin polarization" (DE-SP) procedure for organic polyradicals and so forth does not work well. Thus, we have examined another selection procedure, the "dynamically extended super-exchange" (DE-SE) procedure. DMRG CASCI [18,18] by UB3LYP(HS)-UNO(DE-SE) can realize antiferromagnetic J values similar to experimental ones (-87 cm-1 ). In addition, all J values between all spin states (HS[septet],IS[quintet],IS[triplet],LS[singlet])were also shown to be correct under sufficiently large M values. © 2018 Wiley Periodicals, Inc.

Theoretical and computational investigations of geometrical, electronic and spin structures of the CaMn4 OX (X = 5, 6) cluster in the Kok cycle Si (i = 0-3) of oxygen evolving complex of photosystem II.

The optimized geometries of the CaMn4 OX (X = 5, 6) cluster in the oxygen evolving complex (OEC) of photosystem II (PSII) by large-scale quantum mechanics (QM) and molecular mechanics (MM) calculations are compared with recent serial femtosecond crystallography (SFX) results for the Si (i = 0-3) states. The valence states of four Mn ions by the QM/MM calculations are also examined in relation to the experimental results by the X-ray emission spectroscopy (XES) for the Si intermediates. Geometrical and valence structures of right-opened Mn-hydroxide, Mn-oxo and Mn-peroxide intermediates in the S3 state are investigated in detail in relation to recent SFX and XES experiments for the S3 state. Interplay between theory and experiment indicates that the Mn-oxo intermediate is a new possible candidate for the S3 state. Implications of the computational results are discussed in relation to possible mechanisms of the oxygenoxygen bond formation for water oxidation in OEC of PSII.

Photosubstitution Reaction of cis-[Ru(bpy)2(CH3CN)2]2+ and cis-[Ru(bpy)2(NH3)2]2+ in Aqueous Solution via Monoaqua Intermediate.

The photoinduced ligand exchange reaction of Ru(II) complexes in aqueous solution was studied using density functional theory (DFT). The optimized structures of the lowest triplet state of cis-[Ru(bpy)2(CH3CN)2]2+ (bpy = bipyridine), cis-[Ru(bpy)2(NH3)2]2+, and their monoaqua complexes were analyzed. The metal-centered (3MC) structure was lower than the metal-to-ligand charge transfer (3MLCT) structure for cis-[Ru(bpy)2(CH3CN)2]2+, whereas the 3MLCT structure was lower than the 3MC structure for cis-[Ru(bpy)2(NH3)2]2+. Such a difference would correlate with the higher quantum yield of the former complex. For the monoaqua complexes, the most stable local minimum structure was the 3MC structure, in which the Ru-OH2O and Ru-Nbpy ( trans to the oxygen) bonds were elongated. Therefore, the dissociation of the H2O ligand would be preferred to that of the CH3CN (or NH3) ligand from the monoaqua intermediate, which might result in the reformation of the monoaqua intermediate, and thus, the formation of the bis-aqua product would take a longer time than that of the monoaqua intermediate.

First-Principles Study of the Reaction Mechanism of CHO + H on Graphene Surface.

Many organic molecules observed in the interstellar medium are considered to be formed on dust grains and populated into the gas phase. We analyzed the reaction of HCO + H on a graphene surface using ab initio molecular dynamics simulations as a case study of the formation and desorption of organic molecules on interstellar dust particles. During the reactions of chemisorbed CHO (chemisorbed at the C atom) with free H, CO was generated and efficiently desorbed from the surface. These results suggest that the reactions, of which the reactant forms a covalent bond with the surface while the product does not, cause efficient desorption of the product upon reaction. In such reactions a repulsive force between the product and the surface would be generated and accelerate translation of the product in a specific direction. In addition, it was also shown that the branching ratio of the reactions between radical species on the surface would be affected by the form of the adsorption on the surface, e.g., when a free H reacted with the CHO chemisorbed at the C atom, CH2O was not generated.

Molecular association model of PPARα and its new specific and efficient ligand, pemafibrate: Structural basis for SPPARMα.

Peroxisome proliferator-activated receptor-α (PPARα) is a ligand-activated transcription factor involved in the regulation of lipid homeostasis and improves hypertriglyceridemia. Pemafibrate is a novel selective PPARα modulator (SPPARMα) that activates PPARα transcriptional activity. Here, we computationally constructed the structure of the human PPARα in a complex with pemafibrate, along with that of hPPARα complexed with the classical fenofibrate, and studied their interactions quantitatively by using the first-principles calculations-based fragment molecular orbital (FMO) method. Comprehensive structural and protein-ligand binding elucidation along with the in vitro luciferase analysis let us to identify pemafibrate as a novel SPPARMα. Unlike known fibrate ligands, which bind only with the arm I of the Y-shaped ligand binding pocket, the Y-shaped pemafibrate binds to the entire cavity region. This lock and key nature causes enhanced induced fit in pemafibrate-ligated PPARα. Importantly, this selective modulator allosterically changes PPARα conformation to form a brand-new interface, which in turn binds to PPARα co-activator, PGC-1α, resulting in the full activation of PPARα. The structural basis for the potent effects of pemafibrate on PPARα transcriptional activity predicted by the in silico FMO methods was confirmed by in vitro luciferase assay for mutants. The unique binding mode of pemafibrate reveals a new pattern of nuclear receptor ligand recognition and suggests a novel basis for ligand design, offering cues for improving the binding affinity and selectivity of ligand for better clinical consequences. The findings explain the high affinity and efficacy of pemafibrate, which is expected to be in the clinical use soon.

Redox Potential-Dependent Formation of an Unusual His-Trp Bond in Bilirubin Oxidase.

Bilirubin oxidase (BOD) belongs to the family of blue multicopper oxidases, and catalyzes the concomitant oxidation of bilirubin to biliverdin and the reduction of molecular oxygen to water via a four-electron reduction system. The active sites of BOD comprise four copper atoms; type I copper (T1Cu) forms a mononuclear site, and a cluster of three copper atoms forms a trinuclear center. In the present study, we determined the high-resolution crystal structures of BOD from the fungus Myrothecium verrucaria. We investigated wild-type (WT) BOD and a BOD mutant called Met467Gln, which is inactive against bilirubin. The structures revealed that a novel post-translational crosslink between Trp396 and His398 is formed in the vicinity of the T1Cu site in WT BOD, whereas it is absent in the Met467Gln mutant. Our structural and computational studies suggest that the His-Trp crosslink adjusts the redox potential of T1Cu to that of bilirubin to efficiently abstract electrons from the substrate.

Structural Changes of the Trinuclear Copper Center in Bilirubin Oxidase upon Reduction.

Geometric and electronic structure changes in the copper (Cu) centers in bilirubin oxidase (BOD) upon a four-electron reduction were investigated by quantum mechanics/molecular mechanics (QM/MM) calculations. For the QM region, the unrestricted density functional theory (UDFT) method was adopted for the open-shell system. We found new candidates of the native intermediate (NI, intermediate II) and the resting oxidized (RO) states, i.e., NIH+ and RO0. Elongations of the Cu-Cu atomic distances for the trinuclear Cu center (TNC) and very small structural changes around the type I Cu (T1Cu) were calculated as the results of a four-electron reduction. The QM/MM optimized structures are in good agreement with recent high-resolution X-ray structures. As the structural change in the TNC upon reduction was revealed to be the change in the size of the triangle spanned by the three Cu atoms of TNC, we introduced a new index (l) to characterize the specific structural change. Not only the wild-type, but also the M467Q, which mutates the amino acid residue coordinating T1Cu, were precisely analyzed in terms of their molecular orbital levels, and the optimized redox potential of T1Cu was theoretically reconfirmed.

Large-scale QM/MM calculations of the CaMn4O5 cluster in the S3 state of the oxygen evolving complex of photosystem II. Comparison between water-inserted and no water-inserted structures.

Large-scale QM/MM calculations were performed to elucidate an optimized geometrical structure of a CaMn4O5 cluster with and without water insertion in the S3 state of the oxygen evolving complex (OEC) of photosystem II (PSII). The left (L)-opened structure was found to be stable under the assumption of no hydroxide anion insertion in the S3 state, whereas the right (R)-opened structure became more stable if one water molecule is inserted to the Mn4Ca cluster. The optimized Mna(4)-Mnd(1) distance determined by QM/MM was about 5.0 Å for the S3 structure without an inserted hydroxide anion, but this is elongated by 0.2-0.3 Å after insertion. These computational results are discussed in relation to the possible mechanisms of O-O bond formation in water oxidation by the OEC of PSII.

The reaction mechanism of sarcosine oxidase elucidated using FMO and QM/MM methods.

Monomeric sarcosine oxidase (MSOX) is a flavoprotein that oxidizes sarcosine to the corresponding imine product and is widely used in clinical diagnostics to test renal function. In the past decade, several experimental studies have been performed to elucidate the underlying mechanism of this oxidation reaction. However, the details of the molecular mechanism remain unknown. In this study, we theoretically examined three possible reaction mechanisms, namely, the single-electron transfer, hydride-transfer, and polar mechanisms, using the fragment molecular orbital (FMO) and mixed quantum mechanics/molecular mechanics (QM/MM) methods. We found that, of the three possible reaction pathways, hydride-transfer is the most energetically favorable mechanism. Significantly, hydrogen is not transferred in the hydride state (H-) but in a hydrogen atom state (H˙). Furthermore, a single electron is simultaneously transferred from sarcosine to flavin through their overlapping orbitals. Therefore, based on a detailed theoretical analysis of the calculated reaction pathway, the reaction mechanism of MSOX can be labeled the "hydrogen-atom-coupled electron-transfer" (HACET) mechanism instead of being categorized as the classical hydride-transfer mechanism. QM/MM and FMO calculations revealed that sarcosine is moved close to the flavin ring because of a small charge transfer (about 0.2 electrons in state 1 (MSOX-sarcosine complex)) and that the positively charged residues (Arg49, Arg52, and Lys348) near the active site play a prominent role in stabilizing the sarcosine-flavin complex. These results indicate that strong Coulombic interactions primarily control amine oxidation in the case of MSOX. The new reaction mechanism, HACET, will be important for all the flavoprotein-catalyzed oxidation reactions.

Chemical Equilibrium Models for the S3 State of the Oxygen-Evolving Complex of Photosystem II.

We have performed hybrid density functional theory (DFT) calculations to investigate how chemical equilibria can be described in the S3 state of the oxygen-evolving complex in photosystem II. For a chosen 340-atom model, 1 stable and 11 metastable intermediates have been identified within the range of 13 kcal mol(-1) that differ in protonation, charge, spin, and conformational states. The results imply that reversible interconversion of these intermediates gives rise to dynamic equilibria that involve processes with relocations of protons and electrons residing in the Mn4CaO5 cluster, as well as bound water ligands, with concomitant large changes in the cluster geometry. Such proton tautomerism and redox isomerism are responsible for reversible activation/deactivation processes of substrate oxygen species, through which Mn-O and O-O bonds are transiently ruptured and formed. These results may allow for a tentative interpretation of kinetic data on substrate water exchange on the order of seconds at room temperature, as measured by time-resolved mass spectrometry. The reliability of the hybrid DFT method for the multielectron redox reaction in such an intricate system is also addressed.