Real-time observation of a metal complex-driven reaction intermediate using a porous protein crystal and serial femtosecond crystallography.

Determining short-lived intermediate structures in chemical reactions is challenging. Although ultrafast spectroscopic methods can detect the formation of transient intermediates, real-space structures cannot be determined directly from such studies. Time-resolved serial femtosecond crystallography (TR-SFX) has recently proven to be a powerful method for capturing molecular changes in proteins on femtosecond timescales. However, the methodology has been mostly applied to natural proteins/enzymes and limited to reactions promoted by synthetic molecules due to structure determination challenges. This work demonstrates the applicability of TR-SFX for investigations of chemical reaction mechanisms of synthetic metal complexes. We fix a light-induced CO-releasing Mn(CO)3 reaction center in porous hen egg white lysozyme (HEWL) microcrystals. By controlling light exposure and time, we capture the real-time formation of Mn-carbonyl intermediates during the CO release reaction. The asymmetric protein environment is found to influence the order of CO release. The experimentally-observed reaction path agrees with quantum mechanical calculations. Therefore, our demonstration offers a new approach to visualize atomic-level reactions of small molecules using TR-SFX with real-space structure determination. This advance holds the potential to facilitate design of artificial metalloenzymes with precise mechanisms, empowering design, control and development of innovative reactions.

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.

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

Origin of Homochirality in Amino Acids Induced by Lyman-α Irradiation in the Early Stage of the Milky Way.

The enantiomeric excess (ee) of l-form amino acids found in the Murchison meteorite poses some issues about the cosmic origin of their chirality. Circular dichroism (CD) spectra of amino acids in the far-ultraviolet (FUV) at around 6.8 eV (182 nm) indicate that the circularly polarized light can induce ee through photochemical reactions. Here, we resort to ab initio calculations to extract the CD spectra up to the vacuum-ultraviolet (VUV) region (∼11 eV), and we propose a novel equation to compute the ee applicable to a wider range of light frequency than what is available to date. This allows us to show that the strength of the induced ee (|ee|) in the 10 eV VUV region is comparable to the one in the 6.8 eV FUV region. This feature is common for some key amino acids (alanine, 2-aminobutyric acid, and valine). In space, intense Lyman-α (Lyα) light of 10.2 eV is emitted from star forming regions. This study provides a theoretical basis that Lyα emitter from an early starburst in the Milky Way plays a crucial role in initiating the ee of amino acids.

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.

Enantiomeric Excesses of Aminonitrile Precursors Determine the Homochirality of Amino Acids.

High enantiomeric excesses (ee's) of l-amino acids, including non-proteinogenic amino acid isovaline (Iva), were discovered in the Murchison meteorite, but the detailed molecular mechanism responsible for the observed ee of amino acids remains elusive and inconsistent, because Iva has an inverted circular dichroism (CD) spectrum with respect to α-H amino acids, e.g., alanine. To address this issue, we resort to accurate ab initio calculations for amino acids and their precursors in the Strecker synthesis. We evaluated their photolysis-induced ee in the range 5-11 eV including the Lyman alpha emission line (Lyα), the typical intensive 10.2 eV radiation ascribed to the early phase of galactic evolution. We show that only the aminonitrile precursors are characterized by positive ee in the Lyα region, explaining why right-handed circularly polarized Lyα (R-CP-Lyα) induces homologous l-amino acids. This study shows that the homochirality of amino acids is produced at the aminonitrile precursors stage.

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.

Molecular mechanism of a large conformational change of the quinone cofactor in the semiquinone intermediate of bacterial copper amine oxidase.

Copper amine oxidase from Arthrobacter globiformis (AGAO) catalyses the oxidative deamination of primary amines via a large conformational change of a topaquinone (TPQ) cofactor during the semiquinone formation step. This conformational change of TPQ occurs in the presence of strong hydrogen bonds and neighboring bulky amino acids, especially the conserved Asn381, which restricts TPQ conformational changes over the catalytic cycle. Whether such a semiquinone intermediate is catalytically active or inert has been a matter of debate in copper amine oxidases. Here, we show that the reaction rate of the Asn381Ala mutant decreases 160-fold, and the X-ray crystal structures of the mutant reveals a TPQ-flipped conformation in both the oxidized and reduced states, preceding semiquinone formation. Our hybrid quantum mechanics/molecular mechanics (QM/MM) simulations show that the TPQ conformational change is realized through the sequential steps of the TPQ ring-rotation and slide. We determine that the bulky side chain of Asn381 hinders the undesired TPQ ring-rotation in the oxidized form, favoring the TPQ ring-rotation in reduced TPQ by a further stabilization leading to the TPQ semiquinone form. The acquired conformational flexibility of TPQ semiquinone promotes a high reactivity of Cu(i) to O2, suggesting that the semiquinone form is catalytically active for the subsequent oxidative half-reaction in AGAO. The ingenious molecular mechanism exerted by TPQ to achieve the "state-specific" reaction sheds new light on a drastic environmental transformation around the catalytic center.

Roles of the Flexible Primary Coordination Sphere of the Mn4CaOx Cluster: What Are the Immediate Decay Products of the S3 State?

The primary coordination sphere of the multinuclear cofactor (Mn4CaOx) in the oxygen-evolving complex (OEC) of photosystem II is absolutely conserved to maintain its structure and function. Recent time-resolved serial femtosecond crystallography identified large reorganization of the primary coordination sphere in the S2 to S3 transition, which elicits a cascade of events involving Mn oxidation and water molecule binding to a putative catalytic Mn site. We examined how the crystallographic fields, created by transient conformational states of the OEC at various time points, affect the thermodynamics of various isomers of the Mn cluster using DFT calculations, with an aim of comprehending the functional roles of the flexible primary coordination sphere in the S2 to S3 transition and in the recovery of the S2 state. The results show that the relative movements of surrounding residues change the size and shape of the cavity of the cluster and thereby affect the thermodynamics of various catalytic intermediates as well as the ability to capture a new water molecule at a coordinatively unsaturated site. The implication of these findings is that the protein dynamics may serve to gate the catalytic reaction efficiently by controlling the sequence of Mn oxidation/reduction and water binding/release. This interpretation is consistent with EPR experiments; g ∼ 5 and g ∼ 3 signals obtained after near-infrared (NIR) excitation of the S3 state at 4 K and a g ∼ 5 only signal produced after prolonged incubation of the S3 state at 77 K can be best explained as originating from water-bound S2 clusters (Stotal = 7/2) under a S3 ligand field, i.e., the immediate one-electron reduction products of the oxyl-oxo (Stotal = 6) and hydroxo-oxo (Stotal = 3) species in the S3 state.

Theoretical Investigation into a Possibility of Formation of Propylene Oxide Homochirality in Space.

The preferential synthesis or destruction of a single enantiomer by ultraviolet circularly polarized light (UV-CPL) has been proposed as a possible triggering mechanism for the extraterrestrial origin of homochirality. Herein, we investigate the photoabsorption property of propylene oxide (c-C3H6O) for UV-CPL in the Lyman-α region. Our calculations show that c-C3H6O was produced by CH3+ and CH3CH(OH)CH3 or C3H7• and O (triplet). The computed electronic circular dichroism spectra show that c-C3H6O and the intermediate (CH3CH(OH)CH2+) could absorb the UV-CPL originating from the Lyman-α emitter spectrum, suggesting that the photolysis of c-C3H6O or CH3CH(OH)CH2+ upon irradiation could induce chiral symmetry breakage.

Comprehensive Search of Stable Isomers of Alanine and Alanine Precursors in Prebiotic Syntheses.

Enantiomeric excesses of l-amino acids have been detected in meteorites; however, their molecular mechanism and prebiotic syntheses are still a matter of debate. To elucidate the origin of homochirality, alanine and the chiral precursors formed in prebiotic processes were investigated with regard to their stabilities among their isomers by employing the minimum energy principle, namely, the abundancy of a molecule in the interstellar medium is directly correlated to the stability among isomers. To facilitate the search for possible isomers, we developed a new isomer search algorithm, the random connection method, and performed a thorough search for all the stable isomers within a given chemical formula. We found that alanine and most of its precursors are located at higher energy by more than 5.7 kcal mol-1, with respect to the most stable isomer that consists of a linear-chain structure, whereas only the 2-aminopropanenitrile is the most stable isomer among all others possible. The inherent stability of the α-amino nitrile suggests that the 2-aminopropanenitrile is the dominant contribution in the formation of the common enantiomeric excess over α-amino acids.

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.

Estimation of the relative contributions to the electronic energy transfer rates based on Förster theory: The case of C-phycocyanin chromophores.

In the present study, we provide a reformulation of the theory originally proposed by Förster which allows for simple and convenient formulas useful to estimate the relative contributions of transition dipole moments of a donor and acceptor (chemical factors), their orientation factors (intermolecular structural factors), intermolecular center-to-center distances (intermolecular structural factors), spectral overlaps of absorption and emission spectra (photophysical factors), and refractive index (material factor) to the excitation energy transfer (EET) rate constant. To benchmark their validity, we focused on the EET occurring in C-phycocyanin (C-PC) chromophores. To this aim, we resorted to quantum chemistry calculations to get optimized molecular structures of the C-PC chromophores within the density functional theory (DFT) framework. The absorption and emission spectra, as well as transition dipole moments, were computed by using the time-dependent DFT (TDDFT). Our method was applied to several types of C-PCs showing that the EET rates are determined by an interplay of their specific physical, chemical, and geometrical features. These results show that our formulas can become a useful tool for a reliable estimation of the relative contributions of the factors regulating the EET transfer rate.

pyProGA-A PyMOL plugin for protein residue network analysis.

The field of protein residue network (PRN) research has brought several useful methods and techniques for structural analysis of proteins and protein complexes. Many of these are ripe and ready to be used by the proteomics community outside of the PRN specialists. In this paper we present software which collects an ensemble of (network) methods tailored towards the analysis of protein-protein interactions (PPI) and/or interactions of proteins with ligands of other type, e.g. nucleic acids, oligosaccharides etc. In parallel, we propose the use of the network differential analysis as a method to identify residues mediating key interactions between proteins. We use a model system, to show that in combination with other, already published methods, also included in pyProGA, it can be used to make such predictions. Such extended repertoire of methods allows to cross-check predictions with other methods as well, as we show here. In addition, the possibility to construct PRN models from various kinds of input is so far a unique asset of our code. One can use structural data as defined in PDB files and/or from data on residue pair interaction energies, either from force-field parameters or fragment molecular orbital (FMO) calculations. pyProGA is a free open-source software available from https://gitlab.com/Vlado_S/pyproga.

Weak O2 binding and strong H2O2 binding at the non-heme diiron center of trypanosome alternative oxidase.

Alternative oxidase (AOX) catalyzes the four-electron reduction of dioxygen to water as an additional terminal oxidase, and the catalytic reaction is critical for the parasite to survive in its bloodstream form. Recently, the X-ray crystal structure of trypanosome alternative oxidase (TAO) complexed with ferulenol was reported and the molecular structure of the non-heme diiron center was determined. The binding of O2 was a unique side-on type compared to other iron proteins. In order to characterize the O2 binding state of TAO, the O2 binding states were searched at a quantum mechanics/molecular mechanics (QM/MM) theoretical level in the present study. We found that the most stable O2 binding state is the end-on type, and the binding states of the side-on type are higher in energy. Based on the binding energies and electronic structure analyses, O2 binds very weakly to the TAO iron center (ΔE =6.7 kcal mol-1) in the electronic state of Fe(II)…OO, not in the suggested charge transferred state such as the superoxide state (Fe(III)OO· -) as seen in hemerythrin. Coordination of other ligands such as water, Cl-, CN-, CO, N3- and H2O2 was also examined, and H2O2 was found to bind most strongly to the Fe(II) site by ΔE = 14.0 kcal mol-1. This was confirmed experimentally through the measurement of ubiquinol oxidase activity of TAO and Cryptosporidium parvum AOX which was found to be inhibited by H2O2 in a dose-dependent and reversible manner.

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.

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.

Unique protonation states of aspartate and topaquinone in the active site of copper amine oxidase.

The oxidative deamination of biogenic amines, crucial in the metabolism of a wealth of living organisms, is catalyzed by copper amine oxidases (CAOs). In this work, on the ground of accurate molecular modeling, we provide a clear insight into the unique protonation states of the key catalytic aspartate residue Asp298 and the prosthetic group of topaquinone (TPQ) in the CAO of Arthrobacter globiformis (AGAO). This provides both extensions and complementary information to the crystal structure determined by our recent neutron diffraction (ND) experiment. The hybrid quantum mechanics/molecular mechanics (QM/MM) simulations suggest that the ND structure closely resembles a state in which Asp298 is protonated and the TPQ takes an enolate form. The TPQ keto form can coexist in the fully protonated state. The energetic and structural analyses indicate that the active site structure of the AGAO crystal is not a single state but rather a mixture of the different protonation and conformational states identified in this work.