Novel Approach for Identifying Key Residues in Enzymatic Reactions: Proton Abstraction in Ketosteroid Isomerase

Theoretical Improvements in Enzyme Efficiency Associated with Noisy Rate Constants and Increased Dissipation

Previous studies have revealed the extraordinarily large catalytic efficiency of some enzymes. High catalytic proficiency is an essential accomplishment of biological evolution. Natural selection led to the increased turnover number, kcat, and enzyme efficiency, kcat/KM, of uni-uni enzymes, which convert a single substrate into a single product. We added or multiplied random noise with chosen rate constants to explore the correlation between dissipation and catalytic efficiency for ten enzymes: beta-galactosidase, glucose isomerase, β-lactamases from three bacterial strains, ketosteroid isomerase, triosephosphate isomerase, and carbonic anhydrase I, II, and T200H. Our results highlight the role of biological evolution in accelerating thermodynamic evolution. The catalytic performance of these enzymes is proportional to overall entropy production-the main parameter from irreversible thermodynamics. That parameter is also proportional to the evolutionary distance of β-lactamases PC1, RTEM, and Lac-1 when natural or artificial evolution produces the optimal or maximal possible catalytic efficiency. De novo enzyme design and attempts to speed up the rate-limiting catalytic steps may profit from the described connection between kinetics and thermodynamics.

Computational Methods for Biochemical Simulations Implemented in GAMESS

Computational methods for modeling biochemical processes implemented in GAMESS package are reviewed; in particular, quantum mechanics combined with molecular mechanics (QM/MM), semi-empirical, and fragmentation approaches. A detailed summary of capabilities is provided for the QM/MM implementation in QuanPol program and the fragment molecular orbital (FMO) method. Molecular modeling and visualization packages useful for biochemical simulations with GAMESS are described. GAMESS capabilities with corresponding references are tabulated for reader's convenience.

Fragment molecular orbital study of the cAMP-dependent protein kinase catalyzed phosphoryl transfer: a comparison with the differential transition state stabilization method

The importance of key residues to the activity of the cAMP-dependent protein kinase catalyzed phosphoryl transfer and to the stabilization of the transition state of the reaction has been investigated by means of the fragment molecular orbital (FMO) method. To evaluate the accuracy of the method and its capability of fragmenting covalent bonds, we have compared stabilization energies due to the interactions between individual residues and the reaction center to results obtained with the differential transition state stabilization method (Szarek, et al., J. Phys. Chem. B, 2008, 112, 11819-11826) and observe, despite a size difference in the fragment describing the reaction center, near-quantitative agreement. We have also computed deletion energies to investigate the effect of virtual deletion of key residues on the activation energy. These results are consistent with the stabilization energies and yield additional information as they clearly capture the effect of secondary interactions, i.e. interactions in the second coordination layer of the reaction center. We find that using FMO to calculate deletion energies is a powerful and time efficient approach to analyze the importance of key residues to the activity of an enzyme catalyzed reaction.

Electric Fields and Fast Protein Dynamics in Enzymes

In recent years, there has been much discussion regarding the origin of enzymatic catalysis and whether including protein dynamics is necessary for understanding catalytic enhancement. An important contribution in this debate was made with the application of the vibrational Stark effect spectroscopy to measure electric fields in the active site. This provided a window on electric fields at the transition state in enzymatic reactions. We performed computational studies on two enzymes where we have shown that fast dynamics is part of the reaction mechanism and calculated the electric field near the bond-breaking event. We found that the fast motions that we had identified lead to an increase of the electric field, thus preparing an enzymatic configuration that is electrostatically favorable for the catalytic chemical step. We also studied the enzyme that has been the subject of Stark spectroscopy, ketosteroid isomerase, and found electric fields of a similar magnitude to the two previous examples.

Rapid Estimation of Catalytic Efficiency by Cumulative Atomic Multipole Moments: Application to Ketosteroid Isomerase Mutants

We propose a simple atomic multipole electrostatic model to rapidly evaluate the effects of mutation on enzyme activity and test its performance on wild-type and mutant ketosteroid isomerase. The predictions of our atomic multipole model are similar to those obtained with symmetry-adapted perturbation theory at a fraction of the computational cost. We further show that this approach is relatively insensitive to the precise amino acid side chain conformation in mutants and may thus be useful in computational enzyme (re)design.

[Applications of the Fragment Molecular Orbital Method in Drug Discovery]

  Recently, ab initio quantum mechanical calculations have been applied to large molecules, including biomolecular systems. The fragment molecular orbital (FMO) method is one of the most efficient approaches for the quantum mechanical investigation of such molecules. In the FMO method, dividing a target molecule into small fragments reduces computational effort. The clear definition of inter-fragment interaction energy (IFIE) as an expression of total energy is another valuable feature of the FMO method because it provides the ability to analyze interactions in biomolecules. Thus, the FMO method is expected to be useful for drug discovery. This study demonstrates applications of the FMO method related to drug discovery. First, IFIE, according to FMO calculations, was used in the optimization of drug candidates for the development of anti-prion compounds. The second example involved interaction analysis of the human immunodeficiency virus type 1 (HIV-1) protease and a drug compound that used a novel analytical method for dispersion interaction, i.e., fragment interaction analysis based on LMP2 (FILM).