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.

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.

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.

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

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.

Molecular Mechanism of Depolarization-Dependent Inactivation in W366F Mutant of Kv1.2

Voltage-gated potassium channels play crucial roles in regulating membrane potential. They are activated by membrane depolarization, allowing the selective permeation of K⁺ ions across the plasma membrane, and enter a nonconducting state after lasting depolarization, a process known as inactivation. Inactivation in voltage-activated potassium channels occurs through two distinct mechanisms, N-type and C-type inactivation. C-type inactivation is caused by conformational changes in the extracellular mouth of the channel, whereas N-type inactivation is elicited by changes in the cytoplasmic mouth of the protein. The W434F-mutated Shaker channel is known as a nonconducting mutant and is in a C-type inactivation state at a depolarizing membrane potential. To clarify the structural properties of C-type inactivated protein, we performed molecular dynamics simulations of the wild-type and W366F (corresponding to W434F in Shaker) mutant of the Kv1.2-2.1 chimera channel. The W366F mutant was in a nearly nonconducting state with a depolarizing voltage and recovered from inactivation with a reverse voltage. Our simulations and three-dimensional reference interaction site model analysis suggested that structural changes in the selectivity filter upon membrane depolarization trap K⁺ ions around the inner mouth of the selectivity filter and prevent ion permeation. This pore restriction is involved in the molecular mechanism of C-type inactivation.

The evolutionary process of mammalian sex determination genes focusing on marsupial SRYs

BACKGROUND

Maleness in mammals is genetically determined by the Y chromosome. On the Y chromosome SRY is known as the mammalian male-determining gene. Both placental mammals (Eutheria) and marsupial mammals (Metatheria) have SRY genes. However, only eutherian SRY genes have been empirically examined by functional analyses, and the involvement of marsupial SRY in male gonad development remains speculative.

RESULTS

In order to demonstrate that the marsupial SRY gene is similar to the eutherian SRY gene in function, we first examined the sequence differences between marsupial and eutherian SRY genes. Then, using a parsimony method, we identify 7 marsupial-specific ancestral substitutions, 13 eutherian-specific ancestral substitutions, and 4 substitutions that occurred at the stem lineage of therian SRY genes. A literature search and molecular dynamics computational simulations support that the lineage-specific ancestral substitutions might be involved with the functional differentiation between marsupial and eutherian SRY genes. To address the function of the marsupial SRY gene in male determination, we performed luciferase assays on the testis enhancer of Sox9 core (TESCO) using the marsupial SRY. The functional assay shows that marsupial SRY gene can weakly up-regulate the luciferase expression via TESCO.

CONCLUSIONS

Despite the sequence differences between the marsupial and eutherian SRY genes, our functional assay indicates that the marsupial SRY gene regulates SOX9 as a transcription factor in a similar way to the eutherian SRY gene. Our results suggest that SRY genes obtained the function of male determination in the common ancestor of Theria (placental mammals and marsupials). This suggests that the marsupial SRY gene has a function in male determination, but additional experiments are needed to be conclusive.

Structural basis for the membrane association of ankyrinG via palmitoylation

By clustering various ion channels and transporters, ankyrin-G (AnkG) configures the membrane-excitation platforms in neurons and cardiomyocytes. AnkG itself localizes to specific areas on the plasma membrane via s-palmitoylation of Cys. However, the structural mechanism by which AnkG anchors to the membrane is not understood. In this study, we solved the crystal structures of the reduced and oxidized forms of the AnkG s-palmitoylation domain and used multiple long-term coarse-grained molecular dynamics simulations to analyze their membrane association. Here we report that the membrane anchoring of AnkG was facilitated by s-palmitoylation, defining a stable binding interface on the lipid membrane, and that AnkG without s-palmitoylation also preferred to stay near the membrane but did not have a unique binding interface. This suggests that AnkG in the juxtamembrane region is primed to accept lipid modification at Cys, and once that happens AnkG constitutes a rigid structural base upon which a membrane-excitation platform can be assembled.

Enhanced exchange algorithm without detailed balance condition for replica exchange method

The replica exchange method (REM) is a powerful tool for the conformational sampling of biomolecules. In this study, we propose an enhanced exchange algorithm for REM not meeting the detailed balance condition (DBC), but satisfying the balance condition in all considered exchanges between two replicas. Breaking the DBC can minimize the rejection rate and make an exchange process rejection-free as the number of replicas increases. To enhance the efficiency of REM, all possible pairs--not only the nearest neighbor--were considered in the exchange process. The test simulations of the alanine dipeptide confirmed the correctness of our method. The average traveling distance of each replica in the temperature distribution was also increased in proportion to an increase in the exchange rate. Furthermore, we applied our algorithm to the conformational sampling of the 10-residue miniprotein, chignolin, with an implicit solvent model. The results showed a faster convergence in the calculation of its free energy landscape, compared to that achieved using the normal exchange method of adjacent pairs. This algorithm can also be applied to the conventional near neighbor method and is expected to reduce the required number of replicas.