Search

Hybrid quantum mechanics/molecular mechanics proton transfer calculations are performed using different protein structures and with different retinal geometries and active site water molecules in order to assess whether the transient protonation of Asp212 has participated in deprotonation of Asp85. The results have shown that a transient [O] conformer with protonated Asp212 is sampled during the long-distance proton transfer to the proton release group.

For the photocycle of the membrane protein bacteriorhodopsin to proceed efficiently, the thermal 13-cis to all-trans back-isomerization of the retinal chromophore must return the protein to its resting state on a time-scale of milliseconds. Here, we report on quantum mechanical/molecular mechanical energy calculations examining the structural and energetic determinants of the retinal cis-trans isomerization in the protein environment. The results suggest that a hydrogen-bonded network consisting of the retinal Schiff base, active site amino acid residues, and water molecules can stabilize the twisted retinal, thus reducing the intrinsic energy cost of the cis-trans thermal isomerization barrier. [ABSTRACT FROM AUTHOR]

The first proton transfer in the bacteriorhodopsin photocycle takes place during the L → M transition. Structural details of the pre proton transfer L intermediate have been investigated using experiments and computations. Here, we assess L-state structural models by performing hybrid quantum mechanical/molecular mechanical molecular dynamics and excitation energy calculations. The computations suggest that a water-bridged twisted retinal structure gives the closest agreement with the experimental L/bR shift in the excitation energy. [ABSTRACT FROM AUTHOR]

We present a computational methodology to sample rareevents in large biological enzymes that may involve electronicallypolarizing, reactive processes. The approach includes simultaneousdynamical treatment of electronic and nuclear degrees of freedom,where contributions from the electronic portion are computed usinghybrid density functional theory and the computational costs are reducedthrough a hybrid quantum mechanics/molecular mechanics (QM/MM) treatment.Thus, the paper involves a QM/MM dynamical treatment of rare events.The method is applied to probe the effect of the active site elementson the critical hydrogen transfer step in the soybean lipoxygenase-1(SLO-1) catalyzed oxidation of linoleic acid. It is found that thedynamical fluctuations and associated flexibility of the active siteare critical toward maintaining the electrostatics in the regime wherethe reactive process can occur smoothly. Physical constraints enforcedto limit the active site flexibility are akin to mutations and, inthe cases studied, have a detrimental effect on the electrostaticfluctuations, thus adversely affecting the hydrogen transfer process. [ABSTRACT FROM AUTHOR]

Using ab initio molecular dynamics (AIMD) simulations that facilitate the treatment of rare events, we probe the active site participation in the rate-determining hydrogen transfer step in the catalytic oxidation of linoleic acid by soybean lipoxygenase-1 (SLO-1). The role of two different active site components is probed. (a) On the hydrogen atom acceptor side of the active site, the hydrogen bonding propensity between the acceptor side hydroxyl group, which is bound to the iron cofactor, and the backbone carboxyl group of isoleucine (residue number 839) is studied toward its role in promoting the hydrogen transfer event. Primary and secondary (H/D) isotope effects are also probed and a definite correlation with subtle secondary H/D isotope effects is found. With increasing average nuclear kinetic energy, the increase in transfer probability is enhanced due to the presence of the hydrogen bond between the backbone carbonyl of I839 and the acceptor oxygen. Further increase in average nuclear kinetic energy reduces the strength of this secondary hydrogen bond which leads to a deterioration in hydrogen transfer rates and finally embrances an Arrhenius-like behavior. (b) On the hydrogen atom donor side, the coupling between vibrational modes predominantly localized on the donor-side linoleic acid group and the reactive mode is probed. There appears to be a qualitative difference in the coupling between modes that belong to linoleic acid and the hydrogen transfer mode, for hydrogen and deuterium transfer. For example, the donor side secondary hydrogen atom is much more labile (by nearly a factor of 5) during deuterium transfer as compared to the case for hydrogen transfer. This appears to indicate a greater coupling between the modes belonging to the linoleic acid scaffold and the deuterium transfer mode and also provides a new rationalization for the abnormal (nonclassical) secondary isotope effect results obtained by Knapp, Rickert, and Klinman in J. Am. Chem. Soc. , 2002 , 124 , 3865 . To substantiate our findings noted in point a above, we have suggested an I839 → A839 or I839 → V839 mutation. This will modify the bulkiness of hydrogen the bonding residue, allowing greater flexibility in the secondary hydrogen bond formation highlighted above and adversely affecting the reaction rate.

Identifying the group that acts as the proton storage/loading site is a challenging but important problem for understanding the mechanism of proton pumping in biomolecular proton pumps, such as bacteriorhodopsin (bR) and cytochrome c oxidase. Recent experimental studies of bR propelled the idea that the proton storage/release group (PRG) in bR is not an amino acid but a water cluster embedded in the protein. We argue that this idea is at odds with our knowledge of protein electrostatics, since invoking the water cluster as the PRG would require the protein to raise the pK(a) of a hydronium by almost 11 pK(a) units, which is difficult considering known cases of pK(a) shifts in proteins. Our recent quantum mechanics/molecular mechanics (QM/MM) simulations suggested an alternative "intermolecular proton bond" model in which the stored proton is shared between two conserved Glu residues (194 and 204). Here we show that this model leads to microscopic pK(a) values consistent with available experimental data and the functional requirement of a PRG. Extensive QM/MM simulations also show that, independent of a number of technical issues, such as the influence of QM region size, starting X-ray structure, and nuclear quantum effects, the "intermolecular proton bond" model is qualitatively consistent with available spectroscopic data. Potential of mean force calculations show explicitly that the stored proton strongly prefers the pair of Glu residues over the water cluster. The results and analyses help highlight the importance of considering protein electrostatics and provide arguments for why the "intermolecular proton bond" model is likely applicable to the PRG in biomolecular proton pumps in general.

Proton-transfer reactions in the bacteriorhodopsin light-driven proton pump are coupled with structural rearrangements of protein amino acids and internal water molecules. It is generally thought that the first proton-transfer step from retinal Schiff base to the nearby Asp85 is coupled with movement of the Arg82 side chain away from Asp85 and toward the extracellular proton release group. This movement of Arg82 likely triggers the release of the proton from the proton release group to the extracellular bulk. The exact timing of the movement of Arg82 and how this movement is coupled with proton transfer are still not understood in molecular detail. Here, we address these questions by computing the free energy for the movement of the Arg82 side chain. The calculations indicate that protonation of Asp85 leads to a fast reorientation of the Arg82 side chain toward the extracellular proton release group.