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.

Development of serial X-ray fluorescence holography for radiation-sensitive protein crystals.

X-ray fluorescence holography (XFH) is a powerful atomic resolution technique capable of directly imaging the local atomic structure around atoms of a target element within a material. Although it is theoretically possible to use XFH to study the local structures of metal clusters in large protein crystals, the experiment has proven difficult to perform, especially on radiation-sensitive proteins. Here, the development of serial X-ray fluorescence holography to allow the direct recording of hologram patterns before the onset of radiation damage is reported. By combining a 2D hybrid detector and the serial data collection used in serial protein crystallography, the X-ray fluorescence hologram can be directly recorded in a fraction of the measurement time needed for conventional XFH measurements. This approach was demonstrated by obtaining the Mn Kα hologram pattern from the protein crystal Photosystem II without any X-ray-induced reduction of the Mn clusters. Furthermore, a method to interpret the fluorescence patterns as real-space projections of the atoms surrounding the Mn emitters has been developed, where the surrounding atoms produce large dark dips along the emitter-scatterer bond directions. This new technique paves the way for future experiments on protein crystals that aim to clarify the local atomic structures of their functional metal clusters, and for other related XFH experiments such as valence-selective XFH or time-resolved XFH.

Unraveling Binding Mechanism and Stability of Urease Inhibitors: A QM/MM MD Study.

Soil bacteria can produce urease, which catalyzes the hydrolysis of urea to ammonia (NH3) and carbamate. A variety of urease inhibitors have been proposed to reduce NH3 volatilization by interfering with the urease activity. We report a quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) study on the mechanism employed for the inhibition of urease by three representative competitive inhibitors; namely, acetohydroxamic acid (AHA), hydroxyurea (HU), and N-(n-butyl)phosphorictriamide (NBPTO). The possible connections between the structural and thermodynamical properties and the experimentally observed inhibition efficiency were evaluated and characterized. We demonstrate that the binding affinity decreases in the order NBPTO >> AHA > HU in terms of the computed activation and reaction free energies. This trend also indicates that NBPTO shows the highest inhibitory activity and the lowest IC50 value of 2.1 nM, followed by AHA (42 µM) and HU (100 µM). It was also found that the X=O moiety (X = carbon or phosphorous) plays a crucial role in the inhibitor binding process. These findings not only elucidate why the potent urease inhibitors are effective but also have implications for the design of new inhibitors.

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.

Effects of Active-Center Reduction of Plant-Type Ferredoxin on Its Structure and Dynamics: Computational Analysis Using Molecular Dynamics Simulations.

"Plant-type" ferredoxins (Fds) in the thylakoid membranes of plants, algae, and cyanobacteria possess a single [2Fe-2S] cluster in active sites and mediate light-induced electron transfer from Photosystem I reaction centers to various Fd-dependent enzymes. Structural knowledge of plant-type Fds is relatively limited to static structures, and the detailed behavior of oxidized and reduced Fds has not been fully elucidated. It is important that the investigations of the effects of active-center reduction on the structures and dynamics for elucidating electron-transfer mechanisms. In this study, model systems of oxidized and reduced Fds were constructed from the high-resolution crystal structure of Chlamydomonas reinhardtii Fd1, and three 200 ns molecular dynamics simulations were performed for each system. The force field parameters of the oxidized and reduced active centers were independently obtained using quantum chemical calculations. There were no substantial differences in the global conformations of the oxidized and reduced forms. In contrast, active-center reduction affected the hydrogen-bond network and compactness of the surrounding residues, leading to the increased flexibility of the side chain of Phe61, which is essential for the interaction between Fd and the target protein. These computational results will provide insight into the electron-transfer mechanisms in the Fds.

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.

Depolarizing Effects in Hydrogen Bond Energy in 310-Helices Revealed by Quantum Chemical Analysis.

Hydrogen-bond (H-bond) energies in 310-helices of short alanine peptides were systematically examined by precise DFT calculations with the negative fragmentation approach (NFA), a modified method based on the molecular tailoring approach. The contribution of each H-bond was evaluated in detail from the 310-helical conformation of total energies (whole helical model, WH3-10 model), and the results were compared with the property of H-bond in α-helix from our previous study. The H-bond energies of the WH3-10 model exhibited tendencies different from those exhibited by the α-helix in that they depended on the helical position of the relevant H-bond pair. H-bond pairs adjacent to the terminal H-bond pairs were observed to be strongly destabilized. The analysis of electronic structures indicated that structural characteristics cause the destabilization of the H-bond in 310-helices. We also found that the longer the helix length, the more stable the H-bond in the terminal pairs of the WH3-10 model, suggesting the action of H-bond cooperativity.

Predicting Reaction Mechanisms for the Threonine-Residue Stereoinversion Catalyzed by a Dihydrogen Phosphate Ion.

The stereoinversion of amino acid residues in proteins is considered to trigger various age-related diseases. Serine (Ser) residues are relatively prone to stereoinversion. It is assumed that threonine (Thr) residues also undergo stereoinversion, which results in the formation of the d-allo-Thr residue, by the same mechanisms as those for Ser-residue stereoinversion; however, d-allo-Thr residues have not been detected in vivo. To date, although Ser-residue stereoinversion has been suggested to progress via enolization, plausible reaction mechanisms for Thr-residue stereoinversion have not been proposed. In this study, we investigated the pathway of Thr-residue enolization and successfully identified the three types of plausible reaction pathways of Thr-residue stereoinversion catalyzed by a dihydrogen phosphate ion. The geometries of reactant complexes, transition states, and enolized product complexes were optimized using B3LYP density functional methods, and single-point calculations were performed for all optimized geometries using Møller-Plesset perturbation theory to obtain reliable energies. As a result, the calculated activation energies of Thr-residue stereoinversion were 105-106 kJ mol-1, which were comparable with those of Ser-residue stereoinversion reported previously. The infrequency of Thr-residue stereoinversion may be due to other factors, such as the hydrophobicity and/or the steric hindrance of the γ-methyl group, rather than the high activation energies.

QM/MM Molecular Dynamics Simulations Revealed Catalytic Mechanism of Urease.

Urease catalyzes the hydrolysis of urea to form ammonia and carbamate, inducing an overall pH increase that affects both human health and agriculture. Inhibition, mutagenesis, and kinetic studies have provided insights into its enzymatic role, but there have been debates on the substrate binding mode as well as the reaction mechanism. In the present study, we report quatum mechanics-only (QM-only) and quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) calculations on urease that mainly investigate the binding mode of urea and the mechanism of the urease-catalyzed hydrolysis reaction. Comparison between the experimental data and our QM(GFN2-xTB)/MM metadynamics results demonstrates that urea hydrolysis via a complex with bidentate-bound urea is much more favorable than via that with monodentate-bound urea for both nucleophilic attack and the subsequent proton transfer steps. We also indicate that the bidentate coordination of urea fits the active site with a closed conformation of the mobile flap and can facilitate the stabilization of transition states and intermediates by forming multiple hydrogen bonds with certain active site residues.

Development of a one-step analysis method for several amino acids using a microfluidic paper-based analytical device.

A one-step analysis method was developed for four types of amino acids using a microfluidic paper-based analytical device fabricated from chromatography filtration paper and laminate films. Aminoacyl-tRNA synthetase was used to detect each amino acid. The obtained laminated paper-based analytical device (LPAD) contained four enzymatic reaction areas. Colorimetric detection was performed based on the molybdenum blue reaction. A model method for the simple, easy, and simultaneous detection of several amino acid concentrations was suggested, in contrast to the conventional methods such as HPLC or LC-MS. The method provided a selective quantification at the ranges of 3.6-100 µM for tryptophan, 10.1-100 µM for glycine, 5.9-100 µM for histidine and 5.6-100 µM for lysine with a detection limit of 1.1 µM, 3.3 µM, 1.9 µM and 1.8 µM, respectively. LPAD fabrication was considerably simple, and the subsequent detection process was easy and required a short period of time (within 15 min).

Analysis of Fluctuation in the Heme-Binding Pocket and Heme Distortion in Hemoglobin and Myoglobin.

Heme is located in the active site of proteins and has diverse and important biological functions, such as electron transfer and oxygen transport and/or storage. The distortion of heme porphyrin is considered an important factor for the diverse functions of heme because it correlates with the physical properties of heme, such as oxygen affinity and redox potential. Therefore, clarification of the relationship between heme distortion and the protein environment is crucial in protein science. Here, we analyzed the fluctuation in heme distortion in the protein environment for hemoglobin and myoglobin using molecular dynamics (MD) simulations and quantum mechanical (QM) calculations as well as statistical analysis of the protein structures of hemoglobin and myoglobin stored in Protein Data Bank. Our computation and statistical analysis showed that the protein environment for hemoglobin and myoglobin prominently affects the doming distortion of heme porphyrin, which correlates with its oxygen affinity, and that the magnitude of distortion is different between hemoglobin and myoglobin. These results suggest that heme distortion is affected by its protein environment and fluctuates around its fitted conformation, leading to physical properties that are appropriate for protein functions.

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.

Quantum chemical studies on hydrogen bonds in helical secondary structures.

We present a brief review of our recent computational studies of hydrogen bonds (H-bonds) in helical secondary structures of proteins, α-helix and 310-helix, using a Negative Fragmentation Approach with density functional theory. We found that the depolarized electronic structures of the carbonyl oxygen of the ith residue and the amide hydrogen of the (i + 4)th residue cause weaker H-bond in an α-helix than in an isolated H-bond. Our calculations showed that the H-bond energies in the 310-helix were also weaker than those of the isolated H-bonds. In the 310-helices, the adjacent N-H group at the (i + 1)th residue was closer to the C=O group of the H-bond pair than the adjacent C=O group in the 310-helices, whereas the adjacent C=O group at the (i + 1)th residue was close to the H-bond acceptor in α-helices. Therefore, the destabilization of the H-bond is attributed to the depolarization caused by the adjacent residue of the helical backbone connecting the H-bond donor and acceptor. The differences in the change in electron density revealed that such depolarizations were caused by the local electronic interactions in their neighborhood inside the helical structure and redistributed the electron density. We also present the improvements in the force field of classical molecular simulation, based on our findings. Supplementary Information: The online version contains supplementary material available at 10.1007/s12551-022-01034-5.

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.

Spin-projected QM/MM Free Energy Simulations for Oxidation Reaction of Guanine in B-DNA by Singlet Oxygen.

Guanine is the most susceptible base to oxidation damage induced by reactive oxygen species including singlet oxygen (1 O2 , 1 Δg ). We clarify whether the first step of guanine oxidation in B-DNA proceeds via either a zwitterionic or a diradical intermediate. The free energy profiles are calculated by means of a combined quantum mechanical and molecular mechanical (QM/MM) method coupled with the adaptive biasing force (ABF) method. To describe the open-shell electronic structure of 1 O2 correctly, the broken-symmetry spin-unrestricted density functional theory (BS-UDFT) with an approximate spin projection (AP) correction is applied to the QM region. We find that the effect of spin contamination on the activation and reaction free energies is up to ∼8 kcal mol-1 , which is too large to be neglected. The QM(AP-ULC-BLYP)/MM-based free energy calculations also reveal that the reaction proceeds through a diradical transition state, followed by a conversion to a zwitterionic intermediate. Our computed activation energy of 5.2 kcal mol-1 matches experimentally observed range (0∼6 kcal mol-1 ).

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

Microfluidic Paper-Based Analytical Device for Histidine Determination.

A laminated paper-based analytical device (LPAD) for histidine detection was fabricated from a chromatography filtration paper and laminate films. Histidine recognition was effected by histidyl-tRNA synthetase (HisRS), and its detection was signaled colorimetrically based on the molybdenum blue reaction. The analytical conditions and detectable concentration range of histidine were examined. The method provided selective quantification from 1 to 100 µM histidine. LPAD fabrication is considerably simple, involving only the craft-cutting of the chromatography filtration paper and laminate film, and is cost-effective.

Effects of a remote mutation from the contact paratope on the structure of CDR-H3 in the anti-HIV neutralizing antibody PG16.

PG16 is a broadly neutralizing antibody to the human immunodeficiency virus (HIV). A crystal structure of PG16 revealed that the unusually long 28-residue complementarity determining region (CDR) H3 forms a unique subdomain, referred to as a "hammerhead", that directly contacts the antigen. The hammerhead apparently governs the function of PG16 while a previous experimental assay showed that the mutation of TyrH100Q to Ala, which does not directly contact the antigen, decreased the neutralization ability of PG16. However, the molecular mechanism by which a remote mutation from the hammerhead or contact paratope affects the neutralization potency has remained unclear. Here, we performed molecular dynamics simulations of the wild-type and variants (TyrH100Q to Ala, and TyrH100Q to Phe) of PG16, to clarify the effects of these mutations on the dynamics of CDR-H3. Our simulations revealed that the structural rigidity of the CDR-H3 in PG16 is attributable to the hydrogen bond interaction between TyrH100Q and ProH99, as well as the steric support by TyrH100Q. The loss of both interactions increases the intrinsic fluctuations of the CDR-H3 in PG16, leading to a conformational transition of CDR-H3 toward an inactive state.

Hydrogen bond donors and acceptors are generally depolarized in α-helices as revealed by a molecular tailoring approach.

Hydrogen-bond (H-bond) interaction energies in α-helices of short alanine peptides were systematically examined by precise density functional theory calculations, followed by a molecular tailoring approach. The contribution of each H-bond interaction in α-helices was estimated in detail from the entire conformation energies, and the results were compared with those in the minimal H-bond models, in which only H-bond donors and acceptors exist with the capping methyl groups. The former interaction energies were always significantly weaker than the latter energies, when the same geometries of the H-bond donors and acceptors were applied. The chemical origin of this phenomenon was investigated by analyzing the differences among the electronic structures of the local peptide backbones of the α-helices and those of the minimal H-bond models. Consequently, we found that the reduced H-bond energy originated from the depolarizations of both the H-bond donor and acceptor groups, due to the repulsive interactions with the neighboring polar peptide groups in the α-helix backbone. The classical force fields provide similar H-bond energies to those in the minimal H-bond models, which ignore the current depolarization effect, and thus they overestimate the actual H-bond energies in α-helices. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.