Enzymes control chemical reactions that are key to life processes, and allow them to take place on the time scale needed for synchronization between the relevant reaction cycles. In addition to general interest in their biological roles, these proteins present a fundamental scientific puzzle, since the origin of their tremendous catalytic power is still unclear. While many different hypotheses have been put forward to rationalize this, one of the proposals that has become particularly popular in recent years is the idea that dynamical effects contribute to catalysis. Here, we present a critical review of the dynamical idea, considering all reasonable definitions of what does and does not qualify as a dynamical effect. We demonstrate that no dynamical effect (according to these definitions) has ever been experimentally shown to contribute to catalysis. Furthermore, the existence of non-negligible dynamical contributions to catalysis is not supported by consistent theoretical studies. Our review is aimed, in part, at readers with a background in chemical physics and biophysics, and illustrates that despite a substantial body of experimental effort, there has not yet been any study that consistently established a connection between an enzyme’s conformational dynamics and a significant increase in the catalytic contribution of the chemical step. We also make the point that the dynamical proposal is not a semantic issue but a well-defined scientific hypothesis with well-defined conclusions.

The enzyme catechol O-methyltransferase (COMT) catalyzes the transfer of a methyl group from S-adenosylmethionine to dopamine and related catechols. The search for the origin of COMT catalysis has led to different proposals and hypothesis, including the entropic, the NAC and the compression proposals as well as the more reasonable electrostatic idea. Thus it is important to understand the catalytic power of this enzyme and to examine the validity of different proposals and in particularly the repeated recent implication of the compression idea. The corresponding analysis should be done by well-defined physically based analysis that involves computations rather than circular interpretations of experimental results. Thus we explore here the origin of the catalytic efficiency of COMT by using the empirical valence bond and the linear response approximation approaches. The results demonstrate that the catalytic effect of COMT is mainly due to electrostatic preorganization effects. It is also shown that the compression, NAC and entropic proposals do not account for the catalytic effect.

In this study, the mechanism of dimerization of the full-length Alzheimer amyloid beta (Aβ42) peptide and structural properties of the three most stable dimers have been elucidated through 0.8 μs classical molecular dynamics (MD) simulations. The Aβ42 dimer has been reported to be the smallest neurotoxic species that adversely affects both memory and synaptic plasticity. On the basis of interactions between the distinct regions of the Aβ42 monomer, 10 different starting configurations were developed from their native folded structures. However, only six of them were found to form dimers and among them the three most stable (X(P), C-C(AP), and N-N(P)) were chosen for the detailed analysis. The structural properties of these dimers were compared with the available experimental and theoretical data. The MD simulations show that hydrophobic regions of both monomers play critical roles in the dimerization process. The high content of the α-helical structure in all the dimers is in line with its experimentally proposed role in the oligomerization. The formation of a zipper-like structure in X(P) is also in accordance with its existence in the aggregates of several short amyloidogenic peptides. The computed values of translational (D(T)) and rotational (D(R)) diffusion constants of 0.63 × 10(-6) cm(2)/s and 0.035 ns(-1), respectively, for this dimer are supported by the corresponding values of the Aβ42 monomer. These simulations have also elucidated several other key structural properties of these peptides. This information will be very useful to design small molecules for the inhibition and disruption of the critical Aβ42 dimers.

Computer-aided enzyme design presents a major challenge since in most cases it has not resulted in an impressive catalytic power. The reasons for the problems with computational design include the use of nonquantitative approaches, but they may also reflect other difficulties that are not completely obvious. Thus, it is very useful to try to learn from the trend in directed evolution experiments. Here we explore the nature of the refinement of Kemp eliminases by directed evolution, trying to gain an understanding of related requirements from computational design. The observed trend in the directed evolution refinement of KE07 and HG3 are reproduced, showing that in the case of KE07 the directed evolution leads to ground-state destabilization, whereas in the case of HG3 the directed evolution leads to transition-state stabilization. The nature of the different paths of the directed evolution is examined and discussed. The present study seems to indicate that computer-aided enzyme design may require more than calculations of the effect of single mutations and should be extended to calculations of the effect of simultaneous multiple mutations (that make a few residues preorganized effectively). However, the analysis of two known evolution paths can still be accomplished using the relevant sequences and structures. Thus, by comparing two directed evolution paths of Kemp eliminases we reached the important conclusion that the more effective path leads to transition-state stabilization.

In this B3LYP study, the catalytic mechanisms for the hydrolysis of the three different peptide bonds (Lys28-Gly29, Phe19-Phe20, and His14-Gln15) of Alzheimer amyloid beta (Abeta) peptide by insulin-degrading enzyme (IDE) have been elucidated. For all these peptides, the nature of the substrate was found to influence the structure of the active enzyme-substrate complex. The catalytic mechanism is proposed to proceed through the following three steps: (1) activation of the metal-bound water molecule, (2) formation of the gem-diol intermediate, and (3) cleavage of the peptide bond. With the computed barrier of 14.3, 18.8, and 22.3 kcal/mol for the Lys28-Gly29, Phe19-Phe20, and His14-Gln15 substrates, respectively, the process of water activation was found to be the rate-determining step for all three substrates. The computed energetics show that IDE is the most efficient in hydrolyzing the Lys28-Gly29 (basic polar-neutral nonpolar) peptide bond followed by the Phe19-Phe20 (neutral nonpolar-neutral nonpolar) and His14-Gln15 (basic polar-neutral polar) bonds of the Abeta substrate.