CONVERSION OF LIGHT ENERGY TO ELECTROSTATIC ENERGY IN THE PROTON PUMP OF HALOBACTERIUM HALOBIUM

Abstract— The energetics of the proton pump in Halobacterium halobium is considered on a molecular level. It is proposed that absorption of light by bacteriorhodopsin leads to charge separation and sequential changes in the proton affinities some of the protein groups resulting in proton pumping against the pH gradient. Two alternative models of 'through-space' and 'through-bond' charge separation are proposed. Novel diagrams that describe the energetics of the pumping process are presented. These diagrams can be used both to examine different models of bacteriorhodopsin and for considering other photobiological processes that involve conversion of light energy to electrostatic energy.     

Calculations of the electrostatic free energy contributions to the binding free energy of sulfonamides to carbonic anhydrase

The interactions between biologically important enzymes and drugs are of great interest. In order to address some aspects of these interactions we have initiated a program to investigate enzymedrug interactions. Specifically, the interactions between one of the isozymes of carbonic anhydrase and a family of drugs known as sulfonamides have been studied using computational methods. In particular the electrostatic free energy of binding of carbonic anhydrase II with acetazolamide, methazolamide,p-chlorobenzenesulfonamide,p-aminobenzenesulfonamide and three new compounds (MK1, MK2, and MK3) has been computed using finite-difference Poisson-Boltzmann (FDPB) [1] method and the semimacroscopic version [2, 3] of the protein dipole Langevin dipole (PDLD) method [4]. Both methods, FDPB and PDLD, give similar results for the electrostatic free energy of binding even though different charges and different treatments were used for the protein. The calculated electrostatic binding free energies are in reasonable agreement with the experimental data. The potential and the limitation of electrostatic models for studies of binding energies are discussed.     

A microscopic model for calculations of chemical processes in aqueous solutions

A new microscopic model for calculations of chemical processes in aqueous solutions is presented. The model, referred to here as the “surface constrained soft sphere dipoles” (SCSSD) model avoids the problems of the continuum models by explicitly including the solvent molecules. Each solvent molecule is represented as a point dipole attached to the center of a soft sphere. The solvation energy is evaluated by minimizing the solute-solvent energy with respect to the orientations and positions of those dipoles while constraining the surface dipoles to have the orientations and positions of the bulk solvent. The model is demonstrated by calculating the energetic of charge separation in aqueous solution and evaluating the corresponding dielectric constant. The SCSSD model can be used for quantitative studies of ionic reactions in solutions. This is demonstrated by calculation of the potential surface for the dissociation of formic acid in aqueous solution.     

Associative Versus Dissociative Mechanisms of Phosphate Monoester Hydrolysis: On the Interpretation of Activation Entropies

Phosphate monoester and anhydride hydrolysis is ubiquitous in biology, being involved in, amongst other things, signal transduction, energy production, and the regulation of protein function. Therefore, this reaction has understandably been the focus of intensive research. Nevertheless, the precise mechanism by which phosphate monoester hydrolysis proceeds remains controversial. Traditionally, it has been assumed and frequently implied that a near‐zero activation entropy is indicative of a dissociative pathway. Herein, we examine free‐energy surfaces for the hydrolysis of the methyl phosphate dianion and the methyl pyrophosphate trianion in aqueous solution. In both cases, the reaction can proceed through either compact or expansive concerted (A_ND_N) transition states, with fairly similar barriers. We have evaluated the activation entropies for each transition state and demonstrate that both associative and dissociative transition states have near‐zero entropies of activation that are in good agreement with experimental values. Therefore, we believe that the activation entropy alone is not a useful diagnostic tool, as it depends not only on bond orders at the transition state, but also on other issues that include (but are not limited to) steric factors determining the configurational volumes available to reactants during the reaction, solvation and desolvation effects that may be associated with charge redistribution upon approaching the transition state and entropy changes associated with intramolecular degrees of freedom as the transition state is approached.     

On the mechanism of hydrolysis of phosphate monoesters dianions in solutions and proteins

The nature of the hydrolysis of phosphate monoester dianions in solutions and in proteins is a problem of significant current interest. The present work explores this problem by systematic calculations of the potential surfaces of the reactions of a series of phosphate monoesters with different leaving groups. These calculations involve computational studies ranging from ab initio calculations with implicit solvent models to ab initio QM/MM free energy calculations. The calculations reproduce the observed linear free energy relationship (LFER) for the solution reaction and thus are consistent with the overall experimental trend and can be used to explore the nature of the transition state (TS) region, which is not accessible to direct experimental studies. It is found that the potential surface for the associative and dissociative paths is very flat and that the relative height of the associative and dissociative TS is different in different systems. In general, the character of the TS changes from associative to dissociative upon decrease in the pKa of the leaving group. It is also demonstrated that traditional experimental markers such as isotope effects and the LFER slope cannot be used in a conclusive way to distinguish between the two classes of transition states. In addition it is found that the effective charges of the TS do not follow the previously assumed simple rule. Armed with that experience we explore the free energy surface for the GTPase reaction of the RasGap system. In this case it is found that the surface is flat but that the lowest TS is associative. The present study indicates that the nature of the potential surfaces for the phosphoryl transfer reactions in solution and proteins is quite complicated and cannot be determined in a conclusive way without the use of careful theoretical studies that should, of course, reproduce the available experimental information.     

Simulations of the large kinetic isotope effect and the temperature dependence of the hydrogen atom transfer in lipoxygenase

Elucidating the role of nuclear quantum mechanical (NQM) effects in enzyme catalysis is a topic of significant current interest. Despite the great experimental progress in this field it is important to have theoretical approaches capable of evaluating and analyzing nuclear quantum mechanical contributions to catalysis. In this study, we use the catalytic reaction of lipoxygenase, which is characterized by an extremely large kinetic isotope effect, as a challenging test case for our simulation approach. This is done by applying the quantum classical path (QCP) method with an empirical valence bond potential energy surface. Our computational strategy evaluates the relevant NQM corrections and reproduces the large observed kinetic isotope effect and the temperature dependence of the H atom transfer reaction while being less successful with the D atom transfer reaction. However, the main point of our study is not so much to explore the temperature dependence of the isotope effect but rather to develop and validate an approach for calculations of nuclear quantum mechanical contributions to activation free energies. Here, we find that the deviation between the calculated and observed activation free energies is small for both H and D at all investigated temperatures. The present study also explores the nature of the reorganization energy in the enzyme and solution reactions. It is found that the outer-sphere reorganization energy is extremely small. This reflects the fact that the considered reaction involves a very small charge transfer. The implication of this finding is discussed in the framework of the qualitative vibronic model. The main point of the present study is, however, that the rigorous QCP approach provides a reliable computational tool for evaluating NQM contributions to catalysis even when the given reaction includes large tunneling contributions. Interestingly, our results indicate that the NQM effects in the lipoxygenase reaction are similar in the enzyme and in the reference solution reactions, and thus do not contribute to catalysis. We also reached similar conclusions in studies of other enzymes.