Role of Bulk Water Environment in Regulation of Functional Hydrogen-Bond Network in Photoactive Yellow Protein

Ab Initio Evaluation of the Redox Potential of Cytochrome c

Various biochemical activities of metabolism and biosynthesis are fulfilled by redox processes with explicit electron exchange, which furnish redox enzymes with high chemical reactivity. However, theoretical investigation of a redox process, which simultaneously involves a complex electronic change at a redox metal center and conformational reorganization of the surrounding protein environment coupled to the electronic change, requires computationally conflicting approaches, highly accurate quantum chemical calculations, and long-time molecular dynamics (MD) simulations, limiting the physicochemical understanding of biological redox processes. Here, we theoretically examined a redox process of cytochrome c by means of a hybrid molecular simulation technique, which enables one to consistently treat the redox center at the ab initio quantum chemistry level of theory and the protein reorganization with long-time MD simulations on the microsecond timescale. The calculations successfully evaluated a large absolute redox potential, 4.34 eV, with errors of only 0.03 to 0.34 eV to the experimental ones without any problem-specific empirical parameters. Through the long-time MD sampling, large and nonlinear reorganization of the protein environment was unveiled and the molecular determinants for the redox potential were identified. The present ab initio approach significantly expands the applicability of theoretical investigation to biological redox systems with more electronically complicated redox centers such as polynuclear transition metal complexes.

Chemo-Mechanical Coupling in the Transport Cycle of a Heme ABC Transporter

The heme importer from pathogenic bacteria is a member of the ATP-binding cassette (ABC) transporter family, which uses the energy of ATP-binding and hydrolysis for extensive conformational changes. Previous studies have indicated that conformational changes after heme translocation are triggered by ATP-binding to nucleotide binding domains (NBDs) and then, in turn, induce conformational transitions of the transmembrane domains (TMDs). In this study, we applied a template-based iterative all-atom molecular dynamics (MD) simulation to predict the ATP-bound outward-facing conformation of the Burkholderia cenocepacia heme importer BhuUV-T. The resulting model showed a stable conformation of the TMD with the cytoplasmic gate in the closed state and the periplasmic gate in the open state. Furthermore, targeted MD simulation predicted the intermediate structure of an occluded form (Occ) with bound ATP, in which both ends of the heme translocation channel are closed. The MD simulation of the predicted Occ revealed that Ser147 on the ABC signature motifs (LSGG[Q/E]) of NBDs occasionally flips and loses the active conformation required for ATP-hydrolysis. The flipping motion was found to be coupled to the inter-NBD distance. Our results highlight the functional significance of the signature motif of ABC transporters in regulation of ATPase and chemo-mechanical coupling mechanism.

Theoretical approaches for dynamical ordering of biomolecular systems

BACKGROUND

Living systems are characterized by the dynamic assembly and disassembly of biomolecules. The dynamical ordering mechanism of these biomolecules has been investigated both experimentally and theoretically. The main theoretical approaches include quantum mechanical (QM) calculation, all-atom (AA) modeling, and coarse-grained (CG) modeling. The selected approach depends on the size of the target system (which differs among electrons, atoms, molecules, and molecular assemblies). These hierarchal approaches can be combined with molecular dynamics (MD) simulation and/or integral equation theories for liquids, which cover all size hierarchies.

SCOPE OF REVIEW

We review the framework of quantum mechanical/molecular mechanical (QM/MM) calculations, AA MD simulations, CG modeling, and integral equation theories. Applications of these methods to the dynamical ordering of biomolecular systems are also exemplified.

MAJOR CONCLUSIONS

The QM/MM calculation enables the study of chemical reactions. The AA MD simulation, which omits the QM calculation, can follow longer time-scale phenomena. By reducing the number of degrees of freedom and the computational cost, CG modeling can follow much longer time-scale phenomena than AA modeling. Integral equation theories for liquids elucidate the liquid structure, for example, whether the liquid follows a radial distribution function.

GENERAL SIGNIFICANCE

These theoretical approaches can analyze the dynamic behaviors of biomolecular systems. They also provide useful tools for exploring the dynamic ordering systems of biomolecules, such as self-assembly. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

Neutron crystallography of photoactive yellow protein reveals unusual protonation state of Arg52 in the crystal

Because of its high pK_a, arginine (Arg) is believed to be protonated even in the hydrophobic environment of the protein interior. However, our neutron crystallographic structure of photoactive yellow protein, a light sensor, demonstrated that Arg52 adopts an electrically neutral form. We also showed that the hydrogen bond between the chromophore and Glu46 is a so-called low barrier hydrogen bond (LBHB). Because both the neutral Arg and LBHB are unusual in proteins, these observations remain controversial. To validate our findings, we carried out neutron crystallographic analysis of the E46Q mutant of PYP. The resultant structure revealed that the proportion of the cationic form is higher in E46Q than in WT, although the cationic and neutral forms of Arg52 coexist in E46Q. These observations were confirmed by the occupancy of the deuterium atom bound to the N _η1 atom combined with an alternative conformation of the N_(η2)D₂ group comprising sp² hybridisation. Based on these results, we propose that the formation of the LBHB decreases the proton affinity of Arg52, stabilizing the neutral form in the crystal.

Photocycle of Photoactive Yellow Protein in Cell-Mimetic Environments: Molecular Volume Changes and Kinetics

Using various spectroscopic techniques such as UV-visible spectroscopy, circular dichroism spectroscopy, NMR spectroscopy, small-angle X-ray scattering, transient grating, and transient absorption techniques, we investigated how cell-mimetic environments made by crowding influence the photocycle of photoactive yellow protein (PYP) in terms of the molecular volume change and kinetics. Upon addition of molecular crowding agents, the ratio of the diffusion coefficient of the blue-shifted intermediate (pB) to that of the ground species (pG) significantly changes from 0.92 and approaches 1.0. This result indicates that the molecular volume change accompanied by the photocycle of PYP in molecularly crowded environments is much smaller than that which occurs in vitro and that the pB intermediate under crowded environments favors a compact conformation due to the excluded volume effect. The kinetics of the photocycle of PYP in cell-mimetic environments is greatly decelerated by the dehydration, owing to the interaction between the protein and small crowding agents, but is barely affected by the excluded volume effect. The results lead to the inference that the signaling transducer of PYP may not necessarily utilize the conformational change of PYP to sense the signaling state.