Semiquantitative calculations of catalytic free energies in genetically modified enzymes

The catalytic free energy and binding free energies of the native and the Asn-155----Thr, Asn-155----Leu, and Asn-155----Ala mutants of subtilisin are calculated by the empirical valence bond method and a free energy perturbation method. Two simple procedures are used; one "mutates" the substrate, and the other "mutates" the enzyme. The calculated changes in free energies (delta delta G not equal to cat and delta delta Gbind) between the mutant and native enzymes are within 1 kcal/mol of the corresponding observed values. This indicates that we are approaching a quantitative structure-function correlation. The calculated changes in catalytic free energies are almost entirely due to the electrostatic interaction between the enzyme-water system and the charges of the reacting system. This supports the idea that the electrostatic free energy associated with the changes of charges of the reacting system is the key factor in enzyme catalysis.     

Exploring the origin of the ion selectivity of the KcsA potassium channel

The availability of structural information about biological ion channels provides an opportunity to gain a detailed understanding of the control of ion selectivity by biological systems. However, accomplishing this task by computer simulation approaches is very challenging. First, although the activation barriers for ion transport can be evaluated by microscopic simulations, it is hard to obtain accurate results by such approaches. Second, the selectivity is related to the actual ion current and not directly to the individual activation barriers. Thus, it is essential to simulate the ion currents and this cannot be accomplished at present by microscopic MD approaches. In order to address this challenge, we developed and refined an approach capable of evaluating ion current while still reflecting the realistic features of the given channel. Our method involves generation of semimacroscopic free energy surfaces for the channel/ions system and Brownian dynamics (BD) simulations of the corresponding ion current. In contrast to most alternative macroscopic models, our approach is able to reproduce the difference between the free energy surfaces of different ions and thus to address the selectivity problem. Our method is used in a study of the selectivity of the KcsA channel toward the K+ and Na+ ions. The BD simulations with the calculated free energy profiles produce an appreciable selectivity. To the best of our knowledge, this is the first time that the trend in the selectivity in the ion current is produced by a computer simulation approach. Nevertheless, the calculated selectivity is still smaller than its experimental estimate. Recognizing that the calculated profiles are not perfect, we examine how changes in these profiles can account for the observed selectivity. It is found that the origin of the selectivity is more complex than generally assumed. The observed selectivity can be reproduced by increasing the barrier at the exit and the entrance of the selectivity filter, but the necessary changes in the barrier approach the limit of the error in the PDLD/S-LRA calculations. Other options that can increase the selectivity are also considered, including the difference between the Na+...Na+ and K+...K+ interaction. However, this interesting effect does not appear to lead to a major difference in selectivity since the Na+ ions at the limit of strong interaction tend to move in a less concerted way than the K+ ions. Changes in the relative binding energies at the different binding sites are also not so effective in changing the selectivity. Finally, it is pointed out that using the calculated profiles as a starting point and forcing the model to satisfy different experimentally based constraints, should eventually provide more detailed understanding of the different complex factors involved in ion selectivity of biological channels.     

Computer simulation studies of the fidelity of DNA polymerases

Computer simulations can provide in principle quantitative correlation between the structures of DNA polymerases and the replication fidelity. This paper describes our progress in this direction. Using several theoretical approaches, including the free energy perturbation (FEP), linear response approximation (LRA), and the empirical valence bond (EVB) methods, we examined the stability of several mismatched base pairs in DNA duplex in aqueous solution, the contribution of binding energy to the fidelity of DNA polymerases beta and T7, and the mechanism and energetics of the polymerization reaction catalyzed by T7 DNA polymerase.     

How does GAP catalyze the GTPase reaction of Ras? A computer simulation study

The formation of a complex between p21(ras) and GAP accelerates the GTPase reaction of p21(ras) and terminates the signal for cell proliferation. The understanding of this rate acceleration is important for the elucidation of the role of Ras mutants in tumor formation. In principle there are two main options for the origin of the effect of GAP. One is a direct electrostatic interaction between the residues of GAP and the transition state of the Ras-GAP complex and the other is a GAP-induced shift of the structure of Ras to a configuration that increases the stabilization of the transition state. This work examines the relative importance of these options by computer simulations of the catalytic effect of Ras. The simulations use the empirical valence bond (EVB) method to study the GTPase reaction along the alternative associative and dissociative paths. This approach reproduces the trend in the overall experimentally observed catalytic effect of GAP: the calculated effect is 7 +/- 3 kcal/mol as compared to the observed effect of approximately 6.6 kcal/mol. Furthermore, the calculated effect of mutating Arg789 to a nonpolar residue is 3-4 kcal/mol as compared to the observed effect of 4.5 kcal/mol for the Arg789Ala mutation. It is concluded, in agreement with previous proposals, that the effect of Arg789 is associated with its direct interaction with the transition state charge distribution. However, calculations that use the coordinates of Ras from the Ras-GAP complex (referred to here as Ras') reproduce a significant catalytic effect relative to the Ras coordinates. This indicates that part of the effect of GAP involves a stabilization of a catalytic configuration of Ras. This configuration increases the positive electrostatic potential on the beta-phosphate (relative to the corresponding situation in the free Ras). In other words, GAP stabilizes the GDP bound configuration of Ras relative to that of the GTP-bound conformation. The elusive oncogenic effect of mutating Gln61 is also explored. The calculated effect of such mutations in the Ras-GAP complex are found to be small, while the observed effect is very large (8.7 kcal/mol). Since the Ras is locked in its Ras-GAP configuration in our simulations, we conclude that the oncogenic effect of mutation of Gln61 is indirect and is associated most probably with the structural changes of Ras upon forming the Ras-GAP complex. In view of these and the results for the Ras' we conclude that GAP activates Ras by both direct electrostatic stabilization of the transition state and an indirect allosteric effect that stabilizes the GDP-bound form. The present study also explored the feasibility of the associative and dissociative mechanism in the GTPase reaction of Ras. It is concluded that the reaction is most likely to involve an associative mechanism.     

Reorganization energy of the initial electron-transfer step in photosynthetic bacterial reaction centers.

The reorganization energy (lambda) for electron transfer from the primary electron donor (P*) to the adjacent bacteriochlorophyll (B) in photosynthetic bacterial reaction centers is explored by molecular-dynamics simulations. Relatively long (40 ps) molecular-dynamics trajectories are used, rather than free energy perturbation techniques. When the surroundings of the reaction center are modeled as a membrane, lambda for P* B --> P+ B- is found to be approximately 1.6 kcal/mol. The results are not sensitive to the treatment of the protein's ionizable groups, but surrounding the reaction center with water gives higher values of lambda (approximately 6.5 kcal/mol). In light of the evidence that P+ B- lies slightly below P* in energy, the small lambda obtained with the membrane model is consistent with the speed and temperature independence of photochemical charge separation. The calculated reorganization energy is smaller than would be expected if the molecular dynamics trajectories had sampled the full conformational space of the system. Because the system does not relax completely on the time scale of electron transfer, the lambda obtained here probably is more pertinent than the larger value that would be obtained for a fully equilibrated system.     

Oscillations of the energy gap for the initial electron-transfer step in bacterial reaction centers

Oscillations in the electrostatic energy gap [V_elec(t)] for electron transfer from the primary electron donor (P) to the adjacent bacteriochlorophyll (B) in photosynthetic bacterial reaction centers are examined by molecular-dynamics simulations. Autocorrelation functions of V_elec in the reactant state (PB) include prominent oscillations with an energy of 17 cm^–1. This feature is much weaker if the trajectory is propagated in the product state P⁺B^–. The autocorrelation functions also include oscillations in the regions of 5, 80 and 390 cm^–1 in both states, and near 25 and 48 cm^–1 in P⁺B^–. The strong 17-cm^–1 oscillation could involve motions that modulate the distance between P and B, because a similar oscillation occurs in the direct electrostatic interactions between the electron carriers.