Computational Analysis of the Proton Translocation from Asp96 to Schiff Base in Bacteriorhodopsin

Quantifying Bacteriorhodopsin Activity as a Function of its Local Environment with a Raman-Based Assay

Bacteriorhodopsin (bR) is a transmembrane protein that functions as a light-driven proton pump in halophilic archaea. The bR photocycle has been well-characterized; however, these measurements almost exclusively measured purified bR, outside of its native membrane. To investigate what effect the cellular environment has on the bR photocycle, we have developed a Raman-based assay that can monitor the activity of the bR in a variety of conditions, including in its native membrane. The assay uses two continuous-wave lasers, one to initiate photochemistry and one to monitor bR activity. The excitation leads to the steady-state depletion of ground-state bR, which directly relates to the population of photocycle intermediate states. We have used this assay to monitor bR activity both in vitro and in vivo. Our in vitro measurements confirm that our assay is sensitive to bulk environmental changes reported in the literature. Our in vivo measurements show a decrease in bR activity with increasing extracellular pH for bR in its native membrane. The difference in activity with increasing pH indicates that the native membrane environment affects the function of bR. This assay opens the door to future measurements into understanding how the local environment of this transmembrane protein affects function.

Hydroxide Ion Mechanism for Long-Range Proton Pumping in the Third Proton Transfer of Bacteriorhodopsin

In bacteriorhodopsin, representative light-driven proton pump, five proton transfers yield vectorial active proton translocation, resulting in a proton gradient in microbes. Third proton transfer occurs from Asp96 to the Schiff base on the photocycle, which is expected to be a long-range proton transfer via the Grotthuss mechanism through internal water molecules. Here, large-scale quantum molecular dynamics simulations are performed for the third proton transfer, where all the atoms (~50000 atoms) are treated quantum-mechanically. The simulations demonstrate that two reaction paths exist along the water wire, namely, via hydronium and via hydroxide ions. The free energy analysis confirms that the path via hydroxide ions is considerably favorable and consistent with the observed lifetime of the transient water wire. Therefore, the proposed hydroxide ion mechanism, as in the first proton transfer, is responsible for the third long-range proton transfer.

Prediction of the Structure of Complexes Comprised of Proteins and Glycosaminoglycans Using Docking Simulation and Cluster Analysis

A typical docking simulation provides information on the structure of ligand-receptor complexes and their binding affinity in terms of a docking energy. We have developed a potent method combining a docking simulation with cluster analysis to extract adequate docking structures from the many possible output structures of the simulation. First, we tried to predict the structure of basic fibroblast growth factor (bFGF) bound to heparin, using the docking simulation program AutoDock 3.0. Two X-ray crystal structures had already been obtained for bFGF. One was a complex of the protein and heparin, a kind of glycosaminoglycan, and the other, only the protein itself, hereafter called a simplex. We docked a heparin molecule onto the protein simplex and generated many trial structures for the bFGF-heparin complex. The structures of those docked complexes were optimized through energy minimization by AMBER8. Although neither the docking energy calculated by AMBER8 nor that calculated by AutoDock 3.0 could be used satisfactorily by themselves to select a proper heparin-binding complex from the output structures, the majority of the structures generated by AutoDock 3.0 were fairly close to each other in atom geometry, and the averaged geometry over these structures was also close to that of the crystal. Hence, we utilized only the atom geometry for evaluation and carried out cluster analysis with the collection of geometries. This procedure enabled selection of a structure considerably close to the crystal's. We applied this approach to two other heparin-binding proteins:  antithrombin and annexin V. Two crystal structures, a complex and a simplex, had been elucidated for these proteins as well as for bFGF. Our trials gave an exact prediction of the heparin-binding structures of these proteins, showing the approach in this study is effective in studying the docking of ligands that have a variety of docking conformations due to the presence of multiple rotatable bonds and charged chemical groups.

Deprotonation of D96 in bacteriorhodopsin opens the proton uptake pathway

Despite extensive investigation, the precise mechanism controlling the opening of the cytoplasmic proton uptake pathway in bacteriorhodopsin (bR) has remained a mystery. From an analysis of the X-ray structure of the D96G/F171C/F219L triple mutant of bR and 60 independent molecular dynamics simulations of bR photointermediates, we report that the deprotonation of D96, a key residue in proton transfer reactions, serves two roles that occur sequentially. First, D96 donates a proton to the Schiff base. Subsequently, the deprotonation of D96 serves to "unlatch" the cytoplasmic side. The latching function of D96 appears to be remarkably robust, functioning to open hydration channels in all photointermediate structures. These results suggest that the protonation state of D96 may be the critical biophysical cue controlling the opening and closing of the cytoplasmic half-channel in bR. We suspect that this protonation-switch mechanism could also be utilized in other proton pumps to minimize backflow and reinforce directionality.

Mechanistic and conformational flexibility of the covalent linkage formed during β-lyase activity on an AP-site: application to hOgg1

The β/δ-lyase activity of bifunctional glycosylases on damaged nucleotides in DNA involves the formation of a covalent linkage between the protein (lysine or N-terminal proline) and DNA (C1' of the damaged nucleotide). In the present study, the conformational and mechanistic flexibility of the cross-link is examined. Repair of 8-oxoguanine damage by hOgg1 is considered as a representative system, and the glycosylase through β-lyase steps are investigated using density functional theory. (PCM/SMD)-M06-2X/6-311+G(2df,2p)//PCM-B3LYP/6-31G(d) energetics were determined for eight unique mechanisms differing in the conformation of the imine linkage (E/Z), the proton (pro-S/R) abstracted during elimination, and whether the ring-opening step is base catalyzed. This initial study used a model system limited to the damaged nucleoside 3'-monophosphate and a model nucleophile to investigate this series of complex reaction steps. The great flexibility exhibited by the linkage and clustered β-elimination energetics indicate sterics will play a large role in predicting the preferred lyase mechanism for a given enzyme. The stationary points identified herein can be overlaid into a protein structure to assist in generating initial guesses for large model systems. By comparing the characterized geometries and enzyme active sites, methods for catalysis of the various chemical steps can be identified, and these possibilities are discussed in detail for hOgg1. Interestingly, the most stable structure on the potential energy surface occurs before elimination of the 3'-phosphate. Hydrolysis of the protein-DNA cross-link at this point would yield an AP-site, which provides support for the recently observed monofunctional activity of hOgg1.

A QM/MM-based computational investigation on the catalytic mechanism of saccharopine reductase

Saccharopine reductase from Magnaporthe grisea, an NADPH-containing enzyme in the α-aminoadipate pathway, catalyses the formation of saccharopine, a precursor to L-lysine, from the substrates glutamate and α-aminoadipate-δ-semialdehyde. Its catalytic mechanism has been investigated using quantum mechanics/molecular mechanics (QM/MM) ONIOM-based approaches. In particular, the overall catalytic pathway has been elucidated and the effects of electron correlation and the anisotropic polar protein environment have been examined via the use of the ONIOM(HF/6-31G(d):AMBER94) and ONIOM(MP2/6-31G(d)//HF/6-31G(d):AMBER94) methods within the mechanical embedding formulism and ONIOM(MP2/6-31G(d)//HF/6-31G(d):AMBER94) and ONIOM(MP2/6-311G(d,p)//HF/6-31G(d):AMBER94) within the electronic embedding formulism. The results of the present study suggest that saccharopine reductase utilises a substrate-assisted catalytic pathway in which acid/base groups within the cosubstrates themselves facilitate the mechanistically required proton transfers. Thus, the enzyme appears to act most likely by binding the three required reactant molecules glutamate, α-aminoadipate-δ-semialdehyde and NADPH in a manner and polar environment conducive to reaction.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.