An interpretation of fluctuations in enzyme catalysis rate, spectral diffusion, and radiative component of lifetimes in terms of electric field fluctuations

Time-dependent fluctuations in the catalysis rate ({delta}k(t)) observed in single-enzyme experiments were found in a particular study to have an autocorrelation function decaying on the same time scale as that of spectral diffusion {delta}{omega}0(t). To interpret this similarity, the present analysis focuses on a factor in enzyme catalysis, the local electrostatic interaction energy (E) at the active site and its effect on the activation free energy barrier. We consider the slow fluctuations of the electrostatic interaction energy ({delta}E(t)) as a contributor to {delta}k(t) and relate the latter to {delta}{omega}0(t). The resulting relation between {delta}k(t) and {delta}{omega}0(t) is a dynamic analog of the solvatochromism used in interpreting solvent effects on organic reaction rates. The effect of the postulated {delta}E(t) on fluctuations in the radiative component ({delta}{gamma}Formula(t)) of the fluorescence decay of chromophores in proteins also is examined, and a relation between {delta}{gamma}Formula(t) and {delta}{omega}0(t) is obtained. Experimental tests will determine whether the correlation functions for {delta}k(t), {delta}{omega}0(t), and {delta}{gamma}Formula are indeed similar for any enzyme. Measurements of dielectric dispersion, {varepsilon}({omega}), for the enzyme discussed elsewhere will provide further insight into the correlation function for {delta}E(t). They also will determine whether fluctuations in the nonradiative component {gamma}Formula of the lifetime decay has a different origin, fluctuations in distance for example

Enzyme Promiscuity in Enolase Superfamily. Theoretical Study of o-Succinylbenzoate Synthase Using QM/MM Methods

The promiscuous activity of the enzyme o-succinylbenzoate synthase (OSBS) from the actinobacteria Amycolatopsis is investigated by means of QM/MM methods, using both density functional theory and semiempirical Hamiltonians. This enzyme catalyzes not only the dehydration of 2-succinyl-6R-hydroxy-2,4-cyclohexadiene-1R-carboxylate but also catalyzes racemization of different acylamino acids, with N-succinyl-R-phenylglycine being the best substrate. We investigated the molecular mechanisms for both reactions exploring the potential energy surface. Then, molecular dynamics simulations were performed to obtain the free energy profiles and the averaged interaction energies of enzymatic residues with the reacting system. Our results confirm the plausibility of the reaction mechanisms proposed in the literature, with a good agreement between theoretical and experimentally derived activation free energies. Our simulations unravel the role played by the different residues in each of the two possible reactions. The presence of flexible loops in the active site and the selection of structural modifications in the substrate seem to be key elements to promote the promiscuity of this enzyme.This work was supported by the Spanish Ministerio de Economia y Competitividad project CTQ2012-36253-C03-03 ́ and FEDER funds. K.S. thanks the Polish National Science Center (NCN) for Grant 2011/02/A/ST4/00246. The authors acknowledge computational facilities of the Servei d’Informatica ̀ de la Universitat de Valencia in the ̀ “Tirant” supercomputer, which is part of the Spanish Supercomputing Network

Unraveling the role of protein dynamics in dihydrofolate reductase catalysis

Protein dynamics have controversially been proposed to be at the heart of enzyme catalysis, but identification and analysis of dynamical effects in enzyme-catalyzed reactions have proved very challenging. Here, we tackle this question by comparing an enzyme with its heavy (15N, 13C, 2H substituted) counterpart, providing a subtle probe of dynamics. The crucial hydride transfer step of the reaction (the chemical step) occurs more slowly in the heavy enzyme. A combination of experimental results, quantum mechanics/molecular mechanics simulations, and theoretical analyses identify the origins of the observed differences in reactivity. The generally slightly slower reaction in the heavy enzyme reflects differences in environmental coupling to the hydride transfer step. Importantly, the barrier and contribution of quantum tunneling are not affected, indicating no significant role for “promoting motions” in driving tunneling or modulating the barrier. The chemical step is slower in the heavy enzyme because protein motions coupled to the reaction coordinate are slower. The fact that the heavy enzyme is only slightly less active than its light counterpart shows that protein dynamics have a small, but measurable, effect on the chemical reaction rate

Evolutionarily Conserved Linkage between Enzyme Fold, Flexibility, and Catalysis

Proteins are intrinsically flexible molecules. The role of internal motions in a protein's designated function is widely debated. The role of protein structure in enzyme catalysis is well established, and conservation of structural features provides vital clues to their role in function. Recently, it has been proposed that the protein function may involve multiple conformations: the observed deviations are not random thermodynamic fluctuations; rather, flexibility may be closely linked to protein function, including enzyme catalysis. We hypothesize that the argument of conservation of important structural features can also be extended to identification of protein flexibility in interconnection with enzyme function. Three classes of enzymes (prolyl-peptidyl isomerase, oxidoreductase, and nuclease) that catalyze diverse chemical reactions have been examined using detailed computational modeling. For each class, the identification and characterization of the internal protein motions coupled to the chemical step in enzyme mechanisms in multiple species show identical enzyme conformational fluctuations. In addition to the active-site residues, motions of protein surface loop regions (>10 Å away) are observed to be identical across species, and networks of conserved interactions/residues connect these highly flexible surface regions to the active-site residues that make direct contact with substrates. More interestingly, examination of reaction-coupled motions in non-homologous enzyme systems (with no structural or sequence similarity) that catalyze the same biochemical reaction shows motions that induce remarkably similar changes in the enzyme–substrate interactions during catalysis. The results indicate that the reaction-coupled flexibility is a conserved aspect of the enzyme molecular architecture. Protein motions in distal areas of homologous and non-homologous enzyme systems mediate similar changes in the active-site enzyme–substrate interactions, thereby impacting the mechanism of catalyzed chemistry. These results have implications for understanding the mechanism of allostery, and for protein engineering and drug design

Let’s not forget tautomers

A compound exhibits tautomerism if it can be represented by two structures that are related by an intramolecular movement of hydrogen from one atom to another. The different tautomers of a molecule usually have different molecular fingerprints, hydrophobicities and pKa’s as well as different 3D shape and electrostatic properties; additionally, proteins frequently preferentially bind a tautomer that is present in low abundance in water. As a result, the proper treatment of molecules that can tautomerize, ~25% of a database, is a challenge for every aspect of computer-aided molecular design. Library design that focuses on molecular similarity or diversity might inadvertently include similar molecules that happen to be encoded as different tautomers. Physical property measurements might not establish the properties of individual tautomers with the result that algorithms based on these measurements may be less accurate for molecules that can tautomerize—this problem influences the accuracy of filtering for library design and also traditional QSAR. Any 2D or 3D QSAR analysis must involve the decision of if or how to adjust the observed Ki or IC50 for the tautomerization equilibria. QSARs and recursive partitioning methods also involve the decision as to which tautomer(s) to use to calculate the molecular descriptors. Docking virtual screening must involve the decision as to which tautomers to include in the docking and how to account for tautomerization in the scoring. All of these decisions are more difficult because there is no extensive database of measured tautomeric ratios in both water and non-aqueous solvents and there is no consensus as to the best computational method to calculate tautomeric ratios in different environments

Discovering Conformational Sub-States Relevant to Protein Function

Background: Internal motions enable proteins to explore a range of conformations, even in the vicinity of native state. The role of conformational fluctuations in the designated function of a protein is widely debated. Emerging evidence suggests that sub-groups within the range of conformations (or sub-states) contain properties that may be functionally relevant. However, low populations in these sub-states and the transient nature of conformational transitions between these substates present significant challenges for their identification and characterization. Methods and Findings: To overcome these challenges we have developed a new computational technique, quasianharmonic analysis (QAA). QAA utilizes higher-order statistics of protein motions to identify sub-states in the conformational landscape. Further, the focus on anharmonicity allows identification of conformational fluctuations that enable transitions between sub-states. QAA applied to equilibrium simulations of human ubiquitin and T4 lysozyme reveals functionally relevant sub-states and protein motions involved in molecular recognition. In combination with a reaction pathway sampling method, QAA characterizes conformational sub-states associated with cis/trans peptidyl-prolyl isomerization catalyzed by the enzyme cyclophilin A. In these three proteins, QAA allows identification of conformational sub-states, with critical structural and dynamical features relevant to protein function. Conclusions: Overall, QAA provides a novel framework to intuitively understand the biophysical basis of conformational diversity and its relevance to protein function. © 2011 Ramanathan et al

Temperature-dependent polymorphism of N-(4-fluorophenyl)-1,5-dimethyl-1H-imidazole-4-carboxamide 3-oxide: experimental and theoretical studies on intermolecular interactions in the crystal state

X-ray analysis of N-(4-fluorophenyl)-1,5-dimethyl-1H-imidazole-4-carboxamide 3-oxide reveals the temperature-dependent polymorphism associated with the crystallographic symmetry conversion. The observed crystal structure transformation corresponds to a symmetry reduction from I41 /a (I) to P43 (II) space groups. The phase transition mainly concerns the subtle but clearly noticeable reorganization of molecules in the crystal space, with the structure of individual molecules left almost unchanged. The Hirshfeld surface analysis shows that various intermolecular contacts play an important role in the crystal packing, revealing graphically the differences in spatial arrangements of the molecules in both polymorphs. The N-oxide oxygen atom acts as a formally negatively charged hydrogen bonding acceptor in intramolecular hydrogen bond of N–H…O− type. The combined crystallographic and theoretical DFT methods demonstrate that the observed intramolecular N-oxide N–H…O hydrogen bond should be classified as a very strong charge-assisted and closed-shell non-covalent interaction