Overlooked Hydrogen Bond in a Blue Copper Protein Uncovered by Neutron and Sub-Ångström Resolution X-ray Crystallography.

Metalloproteins play fundamental roles in organisms and are utilized as starting points for the directed evolution of artificial enzymes. Knowing the strategies of metalloproteins, by which they exquisitely tune their activities, will not only lead to an understanding of biochemical phenomena but also contribute to various applications. The blue copper protein (BCP) has been a renowned model system to understand the biology, chemistry, and physics of metalloproteins. Pseudoazurin (Paz), a blue copper protein, mediates electron transfer in the bacterial anaerobic respiratory chain. Its redox potential is finely tuned by hydrogen (H) bond networks; however, difficulty in visualizing H atom positions in the protein hinders the detailed understanding of the protein's structure-function relationship. We here used neutron and sub-ångström resolution X-ray crystallography to directly observe H atoms in Paz. The 0.86-Å-resolution X-ray structure shows that the peptide bond between Pro80 and the His81 Cu ligand deviates from the ideal planar structure. The 1.9-Å-resolution neutron structure confirms a long-overlooked H bond formed by the amide of His81 and the S atom of another Cu ligand Cys78. Quantum mechanics/molecular mechanics calculations show that this H bond increases the redox potential of the Cu site and explains the experimental results well. Our study demonstrates the potential of neutron and sub-ångström resolution X-ray crystallography to understand the chemistry of metalloproteins at atomic and quantum levels.

Investigation on substrate specificity and catalytic activity of serine protease neuropsin.

Neuropsin is one of serine proteases mainly found at the hippocampus and the amygdala, where it contributes to the long-term potentiation and memory acquisition by rebuilding of synaptic connections. Despite of the importance of neuropsin, the substrate specificity and regulation mechanisms of neuropsin have been unclear. Thus, we investigated the substrate specificity and the catalytic activity of neuropsin by the protein-ligand docking and molecular dynamics (MD) simulations and succeeded to reproduce the trend of the experimental results. Our study revealed that the substrate specificity and the activity of neuropsin depended on multiple factors: the substrate charge, the substrate orientation, the hydrogen bond network within the catalytic triad and the substrate, and the formation of the oxyanion hole. The apo neuropsin was not reactive without proper alignment of catalytic triad. The substrate binding induced the reactive alignment of catalytic triad. Then the substrate-neuropsin interaction forms the oxyanion hole that stabilizes the transition state and reduces the free-energy barrier of the following scission reaction.

Molecular Dynamics Simulation Study on Allosteric Regulation of CD44-Hyaluronan Binding as a Force Sensing Mechanism.

CD44 protein exists on surfaces of a variety of human cells, acts as a receptor for the hyaluronan (HA) molecule, and mediates cell adhesion via the HA binding in leukocyte trafficking, cell rolling, and so on. The molecular structures of both CD44 and HA are well known, and the previous work shows that the external-mechanical force induces the partially disordered (PD) conformation from the ordered (O) conformation of CD44. The PD conformation has the higher HA affinity compared to the O conformation. However, the details of force-sensing mechanics have remained unclear. This study provides new insights into allosteric regulation of HA binding by conformational shift from the O to the PD conformation of the CD44 HA binding domain by using the classical molecular dynamics simulations. The O conformation was more favorable than the PD conformation under the equilibrium state, and the O conformation showed weak HA-binding affinity. Our simulation suggests that the PD conformation induced by the external force can refold to a compact structure similar to the O conformation keeping the bound HA. This new conformation showed a higher affinity than the O and PD conformations. Our results show that the unfolding of a remote disordered region from the ligand binding site by the external force allosterically regulates the HA affinity. This study promotes understanding not only the mechanism of CD44-mediated cell rolling but also the allosteric regulation induced by the external mechanical force.

Intra-electron transfer induced by protonation in copper-containing nitrite reductase.

The inter- and intra-electron and proton transfers in the nitrite reduction of copper-containing nitrite reductase (CuNiR) were investigated by using the QM/MM method with the calculational models containing type 1 (T1) and type 2 (T2) Cu sites. The electron transfer from the outer electron donor protein to the T1 Cu site occurred both before and after nitrite binding, and nitrite binding lowered the reduction potential of the Cu T1 site. The protonation of catalytic His244 subsequent to nitrite binding and T1 Cu reduction induced partial intra-electron transfer from T1 to T2 Cu sites. The proton transfer from His244 to nitrite bound on the T2 Cu site via the hydrogen bond network induced intra-electron transfer from the T1 to T2 Cu site. The interaction of the T1 Cu ligand with the second sphere amino acid residues and water molecules affected the reduction potential of the T1 Cu site. The water molecules in the so-called proton pool have an important role in the regulation of the basicity of His244. The conformation of the sensor loop did not change along the reaction, but the water molecule network extending along the sensor loop was changed by nitrite binding.

QM/MM Calculation of the Enzyme Catalytic Cycle Mechanism for Copper- and Zinc-Containing Superoxide Dismutase.

The entire enzyme catalytic mechanism including the electron and the proton transfers of the copper- and zinc-containing extracellular superoxide dismutase (SOD3) was investigated by using QM/MM method. In the first step, the electron transfer from O2·- to SOD3 occurred without the bond formation between the donor and the acceptor and formed the triplet oxygen molecule and reduced SOD3. In the reduced SOD3, the distorted tetrahedral structure of Cu(I) atom was maintained. The reduction of Cu(II) atom induced the protonation of His113, which bridges between the Cu(II) and Zn(II) atoms in the resting state. Since the protonation of His113 broke the bond between Cu(I) and His113, three-coordinated Cu(I) was formed. Further, we suggest the binding of O2·- formed hydrogen peroxide and the resting state after both the Cu reduction and the protonation of His113. The protonation of His113 caused the conformational change of Arg186 located at the entrance of the reactive site. The electrostatic potential surface around the reactive site showed that Arg186 plays an important role as electrostatic guidance for the negatively charged substrates only after the protonation of His113. The rotation of Arg186 switched the proton supply routes via Glu108 or Glu179 for transferring two protons from the bulk solvent.

DFT Study on Enzyme Turnover Including Proton and Electron Transfers of Copper-Containing Nitrite Reductase.

The reaction mechanism of copper-containing nitrite reductase (CuNiR) has been proposed to include two important events, an intramolecular electron transfer and a proton transfer. The two events have been suggested to be coupled, but the order of these events is currently under debate. We investigated the entire enzyme reaction mechanism of nitrite reduction at the T2 Cu site in thermophilic Geobacillus CuNiR from Geobacillus thermodenitrificans NG80-2 (GtNiR) using density functional theory calculations. We found significant conformational changes of His ligands coordinated to the T2 Cu site upon nitrite binding during the catalytic reaction. The reduction potentials and pKa values calculated for the relevant protonation and reduction states show two possible routes, A and B. Reduction of the T2 Cu site in the resting state is followed by endothermic nitrite binding in route A, while exothermic nitrite binding occurs prior to reduction of the T2 Cu site in route B. We concluded that our results support the random-sequential mechanism rather than the ordered mechanism.

DFT Study on Nitrite Reduction Mechanism in Copper-Containing Nitrite Reductase.

Dissimilatory reduction of nitrite by copper-containing nitrite reductase (CuNiR) is an important step in the geobiochemical nitrogen cycle. The proposed mechanisms for the reduction of nitrite by CuNiRs include intramolecular electron and proton transfers, and these two events are understood to couple. Proton-coupled electron transfer is one of the key processes in enzyme reactions. We investigated the geometric structure of bound nitrite and the mechanism of nitrite reduction on CuNiR using density functional theory calculations. Also, the proton transfer pathway, the key residues, and their roles in the reaction mechanism were clarified in this study. In our results, the reduction of T2 Cu site promotes the proton transfer, and the hydrogen bond network around the binding site has an important role not only to stabilize the nitrite binding but also to promote the proton transfer to nitrite.

Structural insights into the function of a thermostable copper-containing nitrite reductase.

Copper-containing nitrite reductase (CuNIR) catalyzes the reduction of nitrite (NO(-)2) to nitric oxide (NO) during denitrification. We determined the crystal structures of CuNIR from thermophilic gram-positive bacterium, Geobacillus thermodenitrificans (GtNIR) in chloride- and formate-bound forms of wild type at 1.15 Šresolution and the nitrite-bound form of the C135A mutant at 1.90 Šresolution. The structure of C135A with nitrite displays a unique η(1)-O coordination mode of nitrite at the catalytic copper site (T2Cu), which has never been observed at the T2Cu site in known wild-type CuNIRs, because the mobility of two residues essential to catalytic activity, Asp98 and His244, are sterically restricted in GtNIR by Phe109 on a characteristic loop structure that is found above Asp98 and by an unusually short CH-O hydrogen bond observed between His244 and water, respectively. A detailed comparison of the WT structure with the nitrite-bound C135A structure implies the replacement of hydrogen-bond networks around His244 and predicts the flow path of protons consumed by nitrite reduction. On the basis of these observations, the reaction mechanism of GtNIR through the η(1)-O coordination manner is proposed.

Free Energy Profile of APOBEC3G Protein Calculated by a Molecular Dynamics Simulation.

The human APOBEC3G protein (A3G) is a single-stranded DNA deaminase that inhibits the replication of retrotransposons and retroviruses, including HIV-1. Atomic details of A3G's catalytic mechanism have started to emerge, as the structure of its catalytic domain (A3Gctd) has been revealed by NMR and X-ray crystallography. The NMR and crystal structures are similar overall; however, differences are apparent for ß2 strand (ß2) and loops close to the catalytic site. To add some insight into these differences and to better characterize A3Gctd dynamics, we calculated its free energy profile by using the Generalized-Born surface area (GBSA) method accompanied with a molecular dynamics simulation. The GBSA method yielded an enthalpy term for A3Gctd's free energy, and we developed a new method that takes into account the distribution of the protein's dihedral angles to calculate its entropy term. The structure solved by NMR was found to have a lower energy than that of the crystal structure, suggesting that this conformation is dominant in solution. In addition, ß2-loop-ß2' configuration was stable throughout a 20-ns molecular dynamics (MD) simulation. This finding suggests that in solution A3Gctd is not likely to adopt the continuous ß2 strand configuration present in the APOBEC2 crystal structure. In the NMR structure, the solvent water accessibility of the catalytic Zn2+ was limited throughout the 20-ns MD simulation. This result explains previous observations in which A3G did not bind or catalyze single cytosine nucleotide, even when at excessive concentrations.

Hydration and temperature dependence of 13C and 1H NMR spectra of the DMPC phospholipid membrane and complete resonance assignment of its crystalline state.

Inhomogeneous line broadening due to conformational distributions of molecules is one of the troublesome problems in solid-state NMR spectroscopy. The best possible way to avoid it is to crystallize the sample. Here, we present a highly resolved (13)C cross-polarization (CP) magic angle spinning (MAS) NMR spectrum of the highly ordered crystalline 1,2-dimyrystoyl-sn-glycero-3-phosphocholine (DMPC) and completely assigned it using two-dimensional (2D) solid-state NMR spectra, dipolar heteronuclear correlation (HETCOR) spectra, scalar heteronuclear J coupling based chemical shift correlation (MAS-J-HMQC) spectra, and Dipolar Assisted Rotational Resonance (DARR) spectra. A comparison between assigned chemical shift values by solid-state NMR in this study and the calculated chemical shift values for X-ray crystal DMPC structures shows good agreement, indicating that the two isomers in the crystalline DMPC take the same conformation as the X-ray crystal structure. The phase diagram of the low hydration level of DMPC (3 ≤ n(W) ≤ 12) determined by (1)H and (13)C NMR spectra indicates that DMPC takes a crystalline state only in a very narrow region around n(W) = 4 and T < 313 K. These findings provide us with conformational information on crystalline DMPC and the physical properties of DMPC at a low hydration level and can possibly help us obtain a highly resolved solid-state NMR spectrum of microcrystalline membrane-associated protein samples.

Development of Chemical Compound Libraries for In Silico Drug Screening.

Chemical compound libraries are the basic database for virtual (in silico) drug screening, and the number of entries has reached 20 million. Many drug-like compound libraries for virtual drug screening have been developed and released. In this review, the process of constructing a database for virtual screening is reviewed, and several popular databases are introduced. Several kinds of focused libraries have been developed. The author has developed databases for metalloproteases, and the details of the libraries are described. The library for metalloproteases was developed by improving the generation of the dominant-ion forms. For instance, the SH group is treated as S- in this library while all SH groups are protonated in the conventional libraries. In addition, metal complexes were examined as new candidates of drug-like compounds. Finally, a method for generating chemical space is introduced, and the diversity of compound libraries is discussed.

Adsorption of small molecules with the hydroxyl group on sodium halide cluster ions.

We have investigated adsorption of molecules with hydroxyl group, ROH, on sodium halide cluster ions, Na(n)X(n-1)(+) (X = F and I, n = 10-17) by mass spectrometry and by theoretical calculations. From analysis of the cluster ion intensities, the adsorption of one water molecule (R = H) is most efficient for Na(13)X(12)(+), whose structure has a NaX defect from a 3 x 3 x 3 cubic structure of n = 14. This result suggests that the defect has an important role in the adsorption reaction. However, it is also found that the reactivity diminishes with increasing bulk size of the R group from H to CH(3), (CH(3))(2)CH, and (CH(3))(3)C. These results imply that the adsorption reactivity is dominated by steric hindrance; the smaller molecules are adsorbed inside the basket structures of Na(13)X(12)(+). Reactivity dependence on the basket size is also discussed by comparing the results of Na(n)F(n-1)(+) and Na(n)I(n-1)(+).

Size-dependent structures of NanI+n-1 cluster ions with a methanol adsorbate: a combined study by photodissociation spectroscopy and density-functional theory calculation.

Methanol adsorption sites on NanI+n-1 ions were investigated. Photoexcitation to charge-transfer states of NanI+n-1 (methanol) predominantly produces two fragment ions: Nan-1I+n-2 (methanol) (neutral NaI loss) and Nan-1I+n-2(neutral NaI and methanol loss), without forming NanI+n-1 (methanol loss). The relative intensities of these fragments are correlated with the geometries and binding energies.