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.

CONSPECTUS: The selective hydrolysis of a peptide or amide bond (-(O═)C-NH-) by a synthetic metallopeptidase is required in a wide range of biological, biotechnological, and industrial applications. In nature, highly specialized enzymes known as proteases and peptidases are used to accomplish this daunting task. Currently, many peptide bond cleaving enzymes and synthetic reagents have been utilized to achieve efficient peptide hydrolysis. However, they possess some serious limitations. To overcome these inadequacies, a variety of metal complexes have been developed that mimic the activities of natural enzymes (metallopeptidases). However, in comparison to metallopeptidases, the hydrolytic reactions facilitated by their existing synthetic analogues are considerably slower and occur with lower catalytic turnover. This could be due to the following reasons: (1) they lack chemical properties of amino acid residues found within enzyme active sites; (2) they contain a higher metal coordination number compared with naturally occurring enzymes; and (3) they do not have access to second coordination shell residues that provide substantial rate enhancements in enzymes. Additionally, the critical structural and mechanistic information required for the development of the next generation of synthetic metallopeptidases cannot be readily obtained through existing experimental techniques. This is because most experimental techniques cannot follow the individual chemical steps in the catalytic cycle due to the fast rate of enzymes. They are also limited by the fact that the diamagnetic d(10) Zn(II) center is silent to electronic, electron spin resonance, and (67)Zn NMR spectroscopies. Therefore, we have employed molecular dynamics (MD), quantum mechanics (QM), and hybrid quantum mechanics/molecular mechanics (QM/MM) techniques to derive this information. In particular, the role of the metal ions, ligands, and microenvironment in the functioning of mono- and binuclear metal center containing enzymes such as insulin degrading enzyme (IDE) and bovine lens leucine aminopeptidase (BILAP), respectively, and their synthetic analogues have been investigated. Our results suggested that in the functioning of IDE, the chemical nature of the peptide bond played a role in the energetics of the reaction and the peptide bond cleavage occurred in the rate-limiting step of the mechanism. In the cocatalytic mechanism used by BILAP, one metal center polarized the scissile peptide bond through the formation of a bond between the metal and the carbonyl group of the substrate, while the second metal center delivered the hydroxyl nucleophile. The Zn(N3) [Zn(His, His, His)] core of matrix metalloproteinase was better than the Zn(N2O) [Zn(His, His, Glu)] core of IDE for peptide hydrolysis. Due to the synergistic interaction between the two metal centers, the binuclear metal center containing Pd2(μ-OH)([18]aneN6)](4+) complex was found to be ∼100 times faster than the mononuclear [Pd(H2O)4](2+) complex. A successful small-molecule synthetic analogue of a ...

Detailed understanding of the action of biological molecular machines must overcome the challenge of gaining a clear knowledge of the corresponding free-energy landscape. An example for this is the elucidation of the nature of converting chemical energy to torque and work in the rotary molecular motor of F1-ATPase. A major part of the challenge involves understanding the rotary–chemical coupling from a non-phenomenological structure/energy description. Here we focused on using a coarse-grained model of F1-ATPase to generate a structure-based free-energy landscape of the rotary–chemical process of the whole system. In particular, we concentrated on exploring the possible impact of the position of the catalytic dwell on the efficiency and torque generation of the molecular machine. It was found that the experimentally observed torque can be reproduced with landscapes that have different positions for the catalytic dwell on the rotary–chemical surface. Thus, although the catalysis is undeniably required for torque generation, the experimentally observed position of the catalytic dwell at 80° might not have a clear advantage for the force generation by F1-ATPase. This further implies that the rotary–chemical couplings in these biological motors are quite robust and their efficiencies do not depend explicitly on the position of the catalytic dwells. Rather, the specific positioning of the dwells with respect to the rotational angle is a characteristic arising due to the structural construct of the molecular machine and might not bear any clear connection to the thermodynamic efficiency for the system.

In this study, interactions of the two full-length Alzheimer amyloid beta peptides (Abeta40 and Abeta42) with the fully active form of insulin degrading enzyme (IDE) through unrestrained, all-atom MD simulations have been investigated. This enzyme is a Zn-containing metallopeptidase that catalyzes the degradation of the monomeric forms of these peptides, and this process is critical for preventing the progression of Alzheimer's disease (AD). The available X-ray structures of the free and small fragment-bound (Asp1-Glu3 and Lys16-Asp23 of Abeta40 and Asp1-Glu3 and Lys16-Glu22 of Abeta42) mutated forms of IDE and NMR structures of the full-length Abeta40 and Abeta42 have been used to build the starting structures for these simulations. The most representative structures derived from the Abeta40-IDE and Abeta42-IDE simulations accurately reproduced the locations of the active site Zn(2+) metal and small fragments of the substrates and their interactions with the enzyme from the X-ray structures. The remaining fragments of both the substrates were found to interact with IDE through several hydrogen bonding, pi-pi, CH-pi, and NH-pi interactions. In comparison to Abeta40, Abeta42 is more flexible and interacts through a smaller number (17-22) of hydrogen bonds in the catalytic chamber of IDE. Both the substrates adopted more beta-sheet character in the IDE environment, an observation that is in line with experiments. Their structural characteristics inside IDE are significantly different than the ones observed in aqueous solution. The atomistic level details provided by these simulations can help in the elucidation of binding and degrading mechanisms of the Abeta peptides by IDE.