Preorganization and protein dynamics in enzyme catalysis

Protein Dynamics and Enzymatic Catalysis

This Perspective presents a review of our work and that of others in the highly controversial topic of the coupling of protein dynamics to reaction in enzymes. We have been involved in studying this topic for many years. Thus, this perspective will naturally present our own views, but it also is designed to present an overview of the variety of viewpoints of this topic, both experimental and theoretical. This is obviously a large and contentious topic.

Complex Loop Dynamics Underpin Activity, Specificity, and Evolvability in the (βα)8 Barrel Enzymes of Histidine and Tryptophan Biosynthesis

Enzymes are conformationally dynamic, and their dynamical properties play an important role in regulating their specificity and evolvability. In this context, substantial attention has been paid to the role of ligand-gated conformational changes in enzyme catalysis; however, such studies have focused on tremendously proficient enzymes such as triosephosphate isomerase and orotidine 5'-monophosphate decarboxylase, where the rapid (μs timescale) motion of a single loop dominates the transition between catalytically inactive and active conformations. In contrast, the (βα)8-barrels of tryptophan and histidine biosynthesis, such as the specialist isomerase enzymes HisA and TrpF, and the bifunctional isomerase PriA, are decorated by multiple long loops that undergo conformational transitions on the ms (or slower) timescale. Studying the interdependent motions of multiple slow loops, and their role in catalysis, poses a significant computational challenge. This work combines conventional and enhanced molecular dynamics simulations with empirical valence bond simulations to provide rich details of the conformational behavior of the catalytic loops in HisA, PriA, and TrpF, and the role of their plasticity in facilitating bifunctionality in PriA and evolved HisA variants. In addition, we demonstrate that, similar to other enzymes activated by ligand-gated conformational changes, loops 3 and 4 of HisA and PriA act as gripper loops, facilitating the isomerization of the large bulky substrate ProFAR, albeit now on much slower timescales. This hints at convergent evolution on these different (βα)8-barrel scaffolds. Finally, our work reemphasizes the potential of engineering loop dynamics as a tool to artificially manipulate the catalytic repertoire of TIM-barrel proteins.

The role of conformational dynamics on the activity of polymer-conjugated CalB in organic solvents

A perennial interest in enzyme catalysis has been expanding its applicability from aqueous phase where enzymes are naturally evolved to organic solvents in which the majority of industrial chemical synthesis...

Networks of Networks: An Essay on Multi-Level Biological Organization

The multi-level organization of nature is self-evident: proteins do interact among them to give rise to an organized metabolism, while in the same time each protein (a single node of such interaction network) is itself a network of interacting amino-acid residues allowing coordinated motion of the macromolecule and systemic effect as allosteric behavior. Similar pictures can be drawn for structure and function of cells, organs, tissues, and ecological systems. The majority of biologists are used to think that causally relevant events originate from the lower level (the molecular one) in the form of perturbations, that “climb up” the hierarchy reaching the ultimate layer of macroscopic behavior (e.g., causing a specific disease). Such causative model, stemming from the usual genotype-phenotype distinction, is not the only one. As a matter of fact, one can observe top-down, bottom-up, as well as middle-out perturbation/control trajectories. The recent complex network studies allow to go further the pure qualitative observation of the existence of both non-linear and non-bottom-up processes and to uncover the deep nature of multi-level organization. Here, taking as paradigm protein structural and interaction networks, we review some of the most relevant results dealing with between networks communication shedding light on the basic principles of complex system control and dynamics and offering a more realistic frame of causation in biology.

Metal Cation-Binding Mechanisms of Q-Proline Peptoid Macrocycles in Solution

The rational design of foldable and functionalizable peptidomimetic scaffolds requires the concerted application of both computational and experimental methods. Recently, a new class of designed peptoid macrocycle incorporating spiroligomer proline mimics (Q-prolines) has been found to preorganize when bound by monovalent metal cations. To determine the solution-state structure of these cation-bound macrocycles, we employ a Bayesian inference method (BICePs) to reconcile enhanced-sampling molecular simulations with sparse ROESY correlations from experimental NMR studies to predict and design conformational and binding properties of macrocycles as functional scaffolds for peptidomimetics. Conformations predicted to be most populated in solution were then simulated in the presence of explicit cations to yield trajectories with observed binding events, revealing a highly preorganized all-trans amide conformation, whose formation is likely limited by the slow rate of cis/trans isomerization. Interestingly, this conformation differs from a racemic crystal structure solved in the absence of cation. Free energies of cation binding computed from distance-dependent potentials of mean force suggest Na⁺ has a higher affinity to the macrocycle than K⁺, with both cations binding much more strongly in acetonitrile than water. The simulated affinities are able to correctly rank the extent to which different macrocycle sequences exhibit preorganization in the presence of different metal cations and solvents, suggesting our approach is suitable for solution-state computational design.

Editing Domain Motions Preorganize the Synthetic Active Site of Prolyl-tRNA Synthetase

Prolyl-tRNA synthetases (ProRSs) catalyze the covalent attachment of proline onto cognate tRNAs, an indispensable step for protein synthesis in all living organisms. ProRSs are modular enzymes and the "prokaryotic-like" ProRSs are distinguished from "eukaryotic-like" ProRSs by the presence of an editing domain (INS) inserted between motifs 2 and 3 of the main catalytic domain. Earlier studies suggested the presence of coupled-domain dynamics could contribute to catalysis; however, the role that the distal, highly mobile INS domain plays in catalysis at the synthetic active site is not completely understood. In the present study, a combination of theoretical and experimental approaches has been used to elucidate the precise role of INS domain dynamics. Quantum mechanical/molecular mechanical simulations were carried out to model catalytic Pro-AMP formation by Enterococcus faecalis ProRS. The energetics of the adenylate formation by the wild-type enzyme was computed and contrasted with variants containing active site mutations, as well as a deletion mutant lacking the INS domain. The combined results revealed that two distinct types of dynamics contribute to the enzyme's catalytic power. One set of motions is intrinsic to the INS domain and leads to conformational preorganization that is essential for catalysis. A second type of motion, stemming from the electrostatic reorganization of active site residues, impacts the height and width of the energy profile and has a critical role in fine tuning the substrate orientation to facilitate reactive collisions. Thus, motions in a distal domain can preorganize the active site of an enzyme to optimize catalysis.

Driving Protein Conformational Cycles in Physiology and Disease: "Frustrated" Amino Acid Interaction Networks Define Dynamic Energy Landscapes: Amino Acid Interaction Networks Change Progressively Along Alpha Tryptophan Synthase's Catalytic Cycle

A general framework by which dynamic interactions within a protein will promote the necessary series of structural changes, or "conformational cycle," required for function is proposed. It is suggested that the free-energy landscape of a protein is biased toward this conformational cycle. Fluctuations into higher energy, although thermally accessible, conformations drive the conformational cycle forward. The amino acid interaction network is defined as those intraprotein interactions that contribute most to the free-energy landscape. Some network connections are consistent in every structural state, while others periodically change their interaction strength according to the conformational cycle. It is reviewed here that structural transitions change these periodic network connections, which then predisposes the protein toward the next set of network changes, and hence the next structural change. These concepts are illustrated by recent work on tryptophan synthase. Disruption of these dynamic connections may lead to aberrant protein function and disease states.

Harnessing Conformational Plasticity to Generate Designer Enzymes

Recent years have witnessed an explosion of interest in understanding the role of conformational dynamics both in the evolution of new enzymatic activities from existing enzymes and in facilitating the emergence of enzymatic activity de novo on scaffolds that were previously non-catalytic. There are also an increasing number of examples in the literature of targeted engineering of conformational dynamics being successfully used to alter enzyme selectivity and activity. Despite the obvious importance of conformational dynamics to both enzyme function and evolvability, many (although not all) computational design approaches still focus either on pure sequence-based approaches or on using structures with limited flexibility to guide the design. However, there exist a wide variety of computational approaches that can be (re)purposed to introduce conformational dynamics as a key consideration in the design process. Coupled with laboratory evolution and more conventional existing sequence- and structure-based approaches, these techniques provide powerful tools for greatly expanding the protein engineering toolkit. This Perspective provides an overview of evolutionary studies that have dissected the role of conformational dynamics in facilitating the emergence of novel enzymes, as well as advances in computational approaches that allow one to target conformational dynamics as part of enzyme design. Harnessing conformational dynamics in engineering studies is a powerful paradigm with which to engineer the next generation of designer biocatalysts.

Crowders Steal Dihydrofolate Reductase Ligands through Quinary Interactions

Dihydrofolate reductase (DHFR) reduces dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. Due to its role in one carbon metabolism, chromosomal DHFR is the target of the antibacterial drug, trimethoprim. Resistance to trimethoprim has resulted in a type II DHFR that is not structurally related to the chromosomal enzyme target. Because of its metabolic significance, understanding DHFR kinetics and ligand binding behavior in more cell-like conditions, where the total macromolecule concentration can be as great as 300 mg/mL, is important. The progress-curve kinetics and ligand binding properties of the drug target (chromosomal E. coli DHFR) and the drug resistant (R67 DHFR) enzymes were studied in the presence of macromolecular cosolutes. There were varied effects on NADPH oxidation and binding to the two DHFRs, with some cosolutes increasing affinity and others weakening binding. However, DHF binding and reduction in both DHFRs decreased in the presence of all cosolutes. The decreased binding of ligands is mostly attributed to weak associations with the macromolecules, as opposed to crowder effects on the DHFRs. Computer simulations found weak, transient interactions for both ligands with several proteins. The net charge of protein cosolutes correlated with effects on NADP⁺ binding, with near neutral and positively charged proteins having more detrimental effects on binding. For DHF binding, effects correlated more with the size of binding pockets on the protein crowders. These nonspecific interactions between DHFR ligands and proteins predict that the in vivo efficiency of DHFRs may be much lower than expected from their in vitro rates.

Exploring PfDHFR reaction surface: A combined molecular dynamics and QM/MM analysis

The substrate to the enzyme PfDHFR (Plasmodium falciparum Dihydrofolate Reductase) is a small molecule dihydrofolate (DHF), it gets converted to tetrahydrofolate (THF) in the active site of the enzyme. The PfDHFR reaction surface involves the protonation of DHF to DHFP as an initial step before the catalytic conversion. The binding affinities of all these species (DHF, DHFP and THF) contribute to the mechanism of DHFR catalytic action. Molecular dynamics (MD) simulations and Quantum Mechanics/Molecular Mechanics (QM/MM) analysis were performed to evaluate the binding affinity and molecular recognition interactions of the substrate DHF/DHFP and the product THF, in the active site of wild-type PfDHFR (wtPfDHFR). The binding affinities of the cofactor NADPH/NADP⁺ were also estimated in all the three complexes. The molecular dynamics (MD) simulations of the substrate, product and cofactor in the cavities of wtPfDHFR revealed the variation of the atomic level interactions during the course of the catalytic conversion. It was found that the DHFP binds very strongly to the PfDHFR active site and pulls the cofactor NADPH closer to itself. The QM/MM analysis revealed that the binding energy of DHFP (-59.82 kcal/mol) and NADPH (-100.24 kcal/mol) in DHFP-wtPfDHFR complex, is higher in comparison to the binding energy of DHF (-38.67 kcal/mol) and NADPH (-77.53 kcal/mol) in DHF-wtPfDHFR complex and the binding energy of THF (-30.72 kcal/mol) and NADP⁺ (-73.72 kcal/mol) in THF-wtPfDHFR complex. The hydride ion donor-acceptor distance (DAD) analysis was also carried out. This combined MD and QM/MM analysis revealed that the protonation of DHF increases the proximity between the substrate and the cofactor, thus facilitates the reaction profile of PfDHFR.

Flexibility of the "rigid" classics or rugged bottom of the folding funnels of myoglobin, lysozyme, RNase A, chymotrypsin, cytochrome c, and carboxypeptidase A1

The abilities to crystalize of a globular protein and to solve its crystal structure seem to represent triumph of the lock-and-key model of protein functionality, where the presence of unique 3D structure resembling aperiodic crystal is considered as a prerequisite for a given protein to possess specific biologic activity. The history of protein crystallography has its roots in first crystal structures of myoglobin, lysozyme, RNase A, chymotrypsin, cytochrome c, and carboxypeptidase A1 solved more than 50 y ago. This article briefly considers extensive structural information currently available for these proteins and shows that the bottoms of their folding funnels (i.e., the lowest parts of their potential energy landscapes) are not smoothed but rugged. In other words, these crystallization classics are characterized by significant conformational flexibility and are not rigid (immobile) crystal-like entities.

The intramolecular hydrogen bonded-halogen bond: a new strategy for preorganization and enhanced binding

Natural and synthetic molecules use weak noncovalent forces to preorganize structure and enable remarkable function. Herein, we introduce the intramolecular hydrogen bonded-halogen bond (HB-XB) as a novel method to preorganize halogen bonding (XBing) molecules, while generating a polarization-enhanced XB. Positioning a fluoroaniline between two iodopyridinium XB donors engendered intramolecular hydrogen bonding (HBing) to the electron-rich belt of both XB donors. NMR solution studies established the efficacy of the HB-XB. The receptor with HB-XBs (G2XB) displayed a nearly 9-fold increase in halide binding over control receptors. Gas-phase density functional theory conformational analysis indicated that the amine stabilizes the bidentate conformation. Furthermore, gas-phase interaction energies showed that the bidentate HB-XBs of G2XBme2+ are more than 3.2 kcal mol-1 stronger than the XBs in a control without the intramolecular HB. Additionally, crystal structures confirm that HB-XBs form tighter contacts with I- and Br- and produce receptors that are more planar. Collectively the results establish the intramolecular HB-XB as a tractable strategy to preorganize XB molecules and regulate XB strength.

High-resolution studies of hydride transfer in the ferredoxin:NADP+ reductase superfamily

Ferredoxin: NADP⁺ reductase (FNR) is an FAD-containing enzyme best known for catalysing the transfer of electrons from ferredoxin (Fd) to NADP⁺ to make NADPH during photosynthesis. It is also the prototype for a broad enzyme superfamily, including the NADPH oxidases (NOXs) that all catalyse similar FAD-enabled electron transfers between NAD(P)H and one-electron carriers. Here, we define further mechanistic details of the NAD(P)H ⇌ FAD hydride-transfer step of the reaction based on spectroscopic studies and high-resolution (~ 1.5 Å) crystallographic views of the nicotinamide-flavin interaction in crystals of corn root FNR Tyr316Ser and Tyr316Ala variants soaked with either nicotinamide, NADP⁺ , or NADPH. The spectra obtained from FNR crystal complexes match those seen in solution and the complexes reveal active site packing interactions and patterns of covalent distortion of the FAD that imply significant active site compression that would favour catalysis. Furthermore, anisotropic B-factors show that the mobility of the C4 atom of the nicotinamide in the FNR:NADP⁺ complex has a directionality matching that expected for boat-like excursions of the nicotinamide ring thought to enhance hydride transfer. Arguments are made for the relevance of this binding mode to catalysis, and specific consideration is given to how the results extrapolate to provide insight to structure-function relations for the membrane-bound NOX enzymes for which little structural information has been available.

DATABASES

Structural data are available in the PDB database under the accession numbers 3LO8 (wild-type), 5VW4 [Y316S:nicotinamide (P3₂ 21)], 5VW9 [Y316S:nicotinamide (P3₁ 21)], 5VW3 [Y316S:NADP⁺ (P3₂ 21)], 5VW8 [Y316S:NADP⁺ (P3₁ 21)], 5VW2 [Y316S:NADPH (P3₂ 21)], 5VW5 [Y316A:nicotinamide (P3₂ 21)], 5VW6 [Y316A:NADP⁺ (P3₂ 21)], 5VW7 [Y316A:NADPH (P3₂ 21)], 5VWA [Y316F (P3₂ 21)], and 5VWB [Y316F:NADP⁺ (P3₁ 21)]. Enzyme Commission number: ferredoxin:NADP⁺ reductase - E C1.18.1.2.

Perspective: Defining and quantifying the role of dynamics in enzyme catalysis

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.

p53 Proteoforms and Intrinsic Disorder: An Illustration of the Protein Structure-Function Continuum Concept

Although it is one of the most studied proteins, p53 continues to be an enigma. This protein has numerous biological functions, possesses intrinsically disordered regions crucial for its functionality, can form both homo-tetramers and isoform-based hetero-tetramers, and is able to interact with many binding partners. It contains numerous posttranslational modifications, has several isoforms generated by alternative splicing, alternative promoter usage or alternative initiation of translation, and is commonly mutated in different cancers. Therefore, p53 serves as an important illustration of the protein structure–function continuum concept, where the generation of multiple proteoforms by various mechanisms defines the ability of this protein to have a multitude of structurally and functionally different states. Considering p53 in the light of a proteoform-based structure–function continuum represents a non-canonical and conceptually new contemplation of structure, regulation, and functionality of this important protein.

Ligand-Dependent Conformational Dynamics of Dihydrofolate Reductase

Enzymes are known to change among several conformational states during turnover. The role of such dynamic structural changes in catalysis is not fully understood. The influence of dynamics in catalysis can be inferred, but not proven, by comparison of equilibrium structures of protein variants and protein-ligand complexes. A more direct way to establish connections between protein dynamics and the catalytic cycle is to probe the kinetics of specific protein motions in comparison to progress along the reaction coordinate. We have examined the enzyme model system dihydrofolate reductase (DHFR) from Escherichia coli with tryptophan fluorescence-probed temperature-jump spectroscopy. We aimed to observe the kinetics of the ligand binding and ligand-induced conformational changes of three DHFR complexes to establish the relationship among these catalytic steps. Surprisingly, in all three complexes, the observed kinetics do not match a simple sequential two-step process. Through analysis of the relationship between ligand concentration and observed rate, we conclude that the observed kinetics correspond to the ligand binding step of the reaction and a noncoupled enzyme conformational change. The kinetics of the conformational change vary with the ligand's identity and presence but do not appear to be directly related to progress along the reaction coordinate. These results emphasize the need for kinetic studies of DHFR with highly specific spectroscopic probes to determine which dynamic events are coupled to the catalytic cycle and which are not.

Mg2+ binds to the surface of thymidylate synthase and affects hydride transfer at the interior active site

Thymidylate synthase (TSase) produces the sole intracellular de novo source of thymidine (i.e., the DNA base T) and thus is a common target for antibiotic and anticancer drugs. Mg(2+) has been reported to affect TSase activity, but the mechanism of this interaction has not been investigated. Here we show that Mg(2+) binds to the surface of Escherichia coli TSase and affects the kinetics of hydride transfer at the interior active site (16 Å away). Examination of the crystal structures identifies a Mg(2+) near the glutamyl moiety of the folate cofactor, providing the first structural evidence for Mg(2+) binding to TSase. The kinetics and NMR relaxation experiments suggest that the weak binding of Mg(2+) to the protein surface stabilizes the closed conformation of the ternary enzyme complex and reduces the entropy of activation on the hydride transfer step. Mg(2+) accelerates the hydride transfer by ~7-fold but does not affect the magnitude or temperature dependence of the intrinsic kinetic isotope effect. These results suggest that Mg(2+) facilitates the protein motions that bring the hydride donor and acceptor together, but it does not change the tunneling ready state of the hydride transfer. These findings highlight how variations in cellular Mg(2+) concentration can modulate enzyme activity through long-range interactions in the protein, rather than binding at the active site. The interaction of Mg(2+) with the glutamyl tail of the folate cofactor and nonconserved residues of bacterial TSase may assist in designing antifolates with polyglutamyl substitutes as species-specific antibiotic drugs.

Protein conformational disorder and enzyme catalysis

Though lacking a well-defined three-dimensional structure, intrinsically unstructured proteins are ubiquitous in nature. These molecules play crucial roles in many cellular processes, especially signaling and regulation. Surprisingly, even enzyme catalysis can tolerate substantial disorder. This observation contravenes conventional wisdom but is relevant to an understanding of how protein dynamics modulates enzyme function. This chapter reviews properties and characteristics of disordered proteins, emphasizing examples of enzymes that lack defined structures, and considers implications of structural disorder for catalytic efficiency and evolution.

Unusual biophysics of intrinsically disordered proteins

Research of a past decade and a half leaves no doubt that complete understanding of protein functionality requires close consideration of the fact that many functional proteins do not have well-folded structures. These intrinsically disordered proteins (IDPs) and proteins with intrinsically disordered protein regions (IDPRs) are highly abundant in nature and play a number of crucial roles in a living cell. Their functions, which are typically associated with a wide range of intermolecular interactions where IDPs possess remarkable binding promiscuity, complement functional repertoire of ordered proteins. All this requires a close attention to the peculiarities of biophysics of these proteins. In this review, some key biophysical features of IDPs are covered. In addition to the peculiar sequence characteristics of IDPs these biophysical features include sequential, structural, and spatiotemporal heterogeneity of IDPs; their rough and relatively flat energy landscapes; their ability to undergo both induced folding and induced unfolding; the ability to interact specifically with structurally unrelated partners; the ability to gain different structures at binding to different partners; and the ability to keep essential amount of disorder even in the bound form. IDPs are also characterized by the "turned-out" response to the changes in their environment, where they gain some structure under conditions resulting in denaturation or even unfolding of ordered proteins. It is proposed that the heterogeneous spatiotemporal structure of IDPs/IDPRs can be described as a set of foldons, inducible foldons, semi-foldons, non-foldons, and unfoldons. They may lose their function when folded, and activation of some IDPs is associated with the awaking of the dormant disorder. It is possible that IDPs represent the "edge of chaos" systems which operate in a region between order and complete randomness or chaos, where the complexity is maximal. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.

A remote mutation affects the hydride transfer by disrupting concerted protein motions in thymidylate synthase

The role of protein flexibility in enzyme-catalyzed activation of chemical bonds is an evolving perspective in enzymology. Here we examine the role of protein motions in the hydride transfer reaction catalyzed by thymidylate synthase (TSase). Being remote from the chemical reaction site, the Y209W mutation of Escherichia coli TSase significantly reduces the protein activity, despite the remarkable similarity between the crystal structures of the wild-type and mutant enzymes with ligands representing their Michaelis complexes. The most conspicuous difference between these two crystal structures is in the anisotropic B-factors, which indicate disruption of the correlated atomic vibrations of protein residues in the mutant. This dynamically altered mutant allows a variety of small thiols to compete for the reaction intermediate that precedes the hydride transfer, indicating disruption of motions that preorganize the protein environment for this chemical step. Although the mutation causes higher enthalpy of activation of the hydride transfer, it only shows a small effect on the temperature dependence of the intrinsic KIE, suggesting marginal changes in the geometry and dynamics of the H-donor and -acceptor at the tunneling ready state. These observations suggest that the mutation disrupts the concerted motions that bring the H-donor and -acceptor together during the pre- and re-organization of the protein environment. The integrated structural and kinetic data allow us to probe the impact of protein motions on different time scales of the hydride transfer reaction within a complex enzymatic mechanism.

Hot spots for allosteric regulation on protein surfaces

Recent work indicates a general architecture for proteins in which sparse networks of physically contiguous and coevolving amino acids underlie basic aspects of structure and function. These networks, termed sectors, are spatially organized such that active sites are linked to many surface sites distributed throughout the structure. Using the metabolic enzyme dihydrofolate reductase as a model system, we show that: (1) the sector is strongly correlated to a network of residues undergoing millisecond conformational fluctuations associated with enzyme catalysis, and (2) sector-connected surface sites are statistically preferred locations for the emergence of allosteric control in vivo. Thus, sectors represent an evolutionarily conserved "wiring" mechanism that can enable perturbations at specific surface positions to rapidly initiate conformational control over protein function. These findings suggest that sectors enable the evolution of intermolecular communication and regulation.

Kinetic and thermodynamic study of cloned thermostable endo-1,4-β-xylanase from Thermotoga petrophila in mesophilic host

The 1,044 bp endo-1,4-β-xylanase gene of a hyperthermophilic Eubacterium, "Thermotoga petrophila RKU 1" (T. petrophila) was amplified, from the genomic DNA of donor bacterium, cloned and expressed in mesophilic host E. coli strain BL21 Codon plus. The extracellular target protein was purified by heat treatment followed by anion and cation exchange column chromatography. The purified enzyme appeared as a single band, corresponding to molecular mass of 40 kDa, upon SDS-PAGE. The pH and temperature profile showed that enzyme was maximally active at 6.0 and 95 °C, respectively against birchwood xylan as a substrate (2,600 U/mg). The enzyme also exhibited marked activity towards beech wood xylan (1,655 U/mg). However minor activity against CMC (61 U/mg) and β-Glucan barley (21 U/mg) was observed. No activity against Avicel, Starch, Laminarin and Whatman filter paper 42 was observed. The K(m), V(max) and K (cat) of the recombinant enzyme were found to be 3.5 mg ml(-1), 2778 μmol mg(-1)min(-1) and 2,137,346.15 s(-1), respectively against birchwood xylan as a substrate. The recombinant enzyme was found very stable and exhibited half life (t(½)) of 54.5 min even at temperature as high as 96 °C, with enthalpy of denaturation (ΔH*(D)), free energy of denaturation (ΔG*(D)) and entropy of denaturation (ΔS*(D)) of 513.23 kJ mol(-1), 104.42 kJ mol(-1) and 1.10 kJ mol(-1)K(-1), respectively at 96 °C. Further the enthalpy (ΔH*), Gibbs free energy (ΔG*) and entropy (ΔS*) for birchwood xylan hydrolysis by recombinant endo-1,4-β-xylanase were calculated at 95 °C as 62.45 kJ mol(-1), 46.18 kJ mol(-1) and 44.2 J mol(-1) K(-1), respectively.

Thermodynamics and solvent effects on substrate and cofactor binding in Escherichia coli chromosomal dihydrofolate reductase

Chromosomal dihydrofolate reductase from Escherichia coli catalyzes the reduction of dihydrofolate to tetrahydrofolate using NADPH as a cofactor. The thermodynamics of ligand binding were examined using an isothermal titration calorimetry approach. Using buffers with different heats of ionization, zero to a small, fractional proton release was observed for dihydrofolate binding, while a proton was released upon NADP(+) binding. The role of water in binding was additionally monitored using a number of different osmolytes. Binding of NADP(+) is accompanied by the net release of ∼5-24 water molecules, with a dependence on the identity of the osmolyte. In contrast, binding of dihydrofolate is weakened in the presence of osmolytes, consistent with "water uptake". Different effects are observed depending on the identity of the osmolyte. The net uptake of water upon dihydrofolate binding was previously observed in the nonhomologous R67-encoded dihydrofolate reductase (dfrB or type II enzyme) [Chopra, S., et al. (2008) J. Biol. Chem. 283, 4690-4698]. As R67 dihydrofolate reductase possesses a nonhomologous sequence and forms a tetrameric structure with a single active site pore, the observation of weaker DHF binding in the presence of osmolytes in both enzymes implicates cosolvent effects on free dihydrofolate. Consistent with this analysis, stopped flow experiments find betaine mostly affects DHF binding via changes in k(on), while betaine mostly affects NADPH binding via changes in k(off). Finally, nonadditive enthalpy terms when binary and ternary cofactor binding events are compared suggest the presence of long-lived conformational transitions that are not included in a simple thermodynamic cycle.

Impact of mutation on proton transfer reactions in ketosteroid isomerase: insights from molecular dynamics simulations

The two proton transfer reactions catalyzed by ketosteroid isomerase (KSI) involve a dienolate intermediate stabilized by hydrogen bonds with Tyr14 and Asp99. Molecular dynamics simulations based on an empirical valence bond model are used to examine the impact of mutating these residues on the hydrogen-bonding patterns, conformational changes, and van der Waals and electrostatic interactions during the proton transfer reactions. While the rate constants for the two proton transfer steps are similar for wild-type (WT) KSI, the simulations suggest that the rate constant for the first proton transfer step is smaller in the mutants due to the significantly higher free energy of the dienolate intermediate relative to the reactant. The calculated rate constants for the mutants D99L, Y14F, and Y14F/D99L relative to WT KSI are qualitatively consistent with the kinetic experiments indicating a significant reduction in the catalytic rates along the series of mutants. In the simulations, WT KSI retained two hydrogen-bonding interactions between the substrate and the active site, while the mutants typically retained only one hydrogen-bonding interaction. A new hydrogen-bonding interaction between the substrate and Tyr55 was observed in the double mutant, leading to the prediction that mutation of Tyr55 will have a greater impact on the proton transfer rate constants for the double mutant than for WT KSI. The electrostatic stabilization of the dienolate intermediate relative to the reactant was greater for WT KSI than for the mutants, providing a qualitative explanation for the significantly reduced rates of the mutants. The active site exhibited restricted motion during the proton transfer reactions, but small conformational changes occurred to facilitate the proton transfer reactions by strengthening the hydrogen-bonding interactions and by bringing the proton donor and acceptor closer to each other with the proper orientation for proton transfer. Thus, these calculations suggest that KSI forms a preorganized active site but that the structure of this preorganized active site is altered upon mutation. Moreover, small conformational changes due to stochastic thermal motions are required within this preorganized active site to facilitate the proton transfer reactions.

At the dawn of the 21st century: Is dynamics the missing link for understanding enzyme catalysis?

Enzymes play a key role in almost all biological processes, accelerating a variety of metabolic reactions as well as controlling energy transduction, the transcription, and translation of genetic information, and signaling. They possess the remarkable capacity to accelerate reactions by many orders of magnitude compared to their uncatalyzed counterparts, making feasible crucial processes that would otherwise not occur on biologically relevant timescales. Thus, there is broad interest in understanding the catalytic power of enzymes on a molecular level. Several proposals have been put forward to try to explain this phenomenon, and one that has rapidly gained momentum in recent years is the idea that enzyme dynamics somehow contributes to catalysis. This review examines the dynamical proposal in a critical way, considering basically all reasonable definitions, including (but not limited to) such proposed effects as "coupling between conformational and chemical motions," "landscape searches" and "entropy funnels." It is shown that none of these proposed effects have been experimentally demonstrated to contribute to catalysis, nor are they supported by consistent theoretical studies. On the other hand, it is clarified that careful simulation studies have excluded most (if not all) dynamical proposals. This review places significant emphasis on clarifying the role of logical definitions of different catalytic proposals, and on the need for a clear formulation in terms of the assumed potential surface and reaction coordinate. Finally, it is pointed out that electrostatic preorganization actually accounts for the observed catalytic effects of enzymes, through the corresponding changes in the activation free energies.

Conformational differences between the wild type and V30M mutant transthyretin modulate its binding to genistein: implications to tetramer stability and ligand-binding

Transthyretin (TTR) is a tetrameric beta-sheet-rich transporter protein directly involved in human amyloid diseases. It was recently found that the isoflavone genistein (GEN) potently inhibits TTR amyloid fibril formation (Green et al., 2005) and is therefore a promising candidate for TTR amyloidosis treatment. Here we used structural and biophysical approaches to characterize genistein binding to the wild type (TTRwt) and to its most frequent amyloidogenic variant, the V30M mutant. In a dose-dependent manner, genistein elicited considerable increases in both mutant and TTRwt stability as demonstrated by high hydrostatic pressure (HHP) and acid-mediated dissociation/denaturation assays. TTR:GEN crystal complexes and isothermal titration calorimetry (ITC) experiments showed that the binding mechanisms of genistein to the TTRwt and to V30M are different and are dependent on apoTTR structure conformations. Furthermore, we could also identify potential allosteric movements caused by genistein binding to the wild type TTR that explains, at least in part, the frequently observed negatively cooperative process between the two sites of TTRwt when binding ligands. These findings show that TTR mutants may present different ligand recognition and therefore are of value in ligand design for inhibiting TTR amyloidosis.

Coupling effects of distal loops on structural stability and enzymatic activity of Escherichia coli dihydrofolate reductase revealed by deletion mutants

Residues distal from the active site in dihydrofolate reductase (DHFR) have regulatory roles in catalytic reaction and also folding stability. The couplings of the distal residues to the ones in the active site have been analyzed using site-directed mutants. To expand our understanding of the structural and functional influences of distal residue mutation, we explored the structural stability and enzymatic activity of deletion mutants. Deletion has greater structural and dynamical impacts on the corresponding part than site-directed mutation does. Thus, deletion amplifies the effects caused by distal mutations, which should make the mutual couplings among the distant residues more apparent. We focused on residues 52, 67, 121, and 145 in the four distinct loops of DHFR. All the single-residue deletion mutants showed marked reduction in stability, except for Delta52 in an alphaC-betaC loop. Double deletion mutants showed that the loop alphaC-betaC has nonadditive couplings with the betaF-betaG and betaG-betaH loops regarding stability. Single deletion to the loops alphaC-betaC or betaC-betaD resulted in considerable activity reduction, demonstrating that the loops couple to the residues near the active site. The four loops were shown to be functionally interdependent from the double deletion experiments.

Surface sites for engineering allosteric control in proteins

Statistical analyses of protein families reveal networks of coevolving amino acids that functionally link distantly positioned functional surfaces. Such linkages suggest a concept for engineering allosteric control into proteins: The intramolecular networks of two proteins could be joined across their surface sites such that the activity of one protein might control the activity of the other. We tested this idea by creating PAS-DHFR, a designed chimeric protein that connects a light-sensing signaling domain from a plant member of the Per/Arnt/Sim (PAS) family of proteins with Escherichia coli dihydrofolate reductase (DHFR). With no optimization, PAS-DHFR exhibited light-dependent catalytic activity that depended on the site of connection and on known signaling mechanisms in both proteins. PAS-DHFR serves as a proof of concept for engineering regulatory activities into proteins through interface design at conserved allosteric sites.

Catalytic cycle of human glutathione reductase near 1 A resolution

Efficient enzyme catalysis depends on exquisite details of structure beyond those resolvable in typical medium- and high-resolution crystallographic analyses. Here we report synchrotron-based cryocrystallographic studies of natural substrate complexes of the flavoenzyme human glutathione reductase (GR) at nominal resolutions between 1.1 and 0.95 A that reveal new aspects of its mechanism. Compression in the active site causes overlapping van der Waals radii and distortion in the nicotinamide ring of the NADPH substrate, which enhances catalysis via stereoelectronic effects. The bound NADPH and redox-active disulfide are positioned optimally on opposite sides of the flavin for a 1,2-addition across a flavin double bond. The new structures extend earlier observations to reveal that the redox-active disulfide loop in GR is an extreme case of sequential peptide bonds systematically deviating from planarity--a net deviation of 53 degrees across five residues. But this apparent strain is not a factor in catalysis, as it is present in both oxidized and reduced structures. Intriguingly, the flavin bond lengths in oxidized GR are intermediate between those expected for oxidized and reduced flavin, but we present evidence that this may not be due to the protein environment but instead due to partial synchrotron reduction of the flavin by the synchrotron beam. Finally, of more general relevance, we present evidence that the structures of synchrotron-reduced disulfide bonds cannot generally be used as reliable models for naturally reduced disulfide bonds.

R67, the other dihydrofolate reductase: rational design of an alternate active site configuration

R67 dihydrofolate reductase (DHFR) bears no sequence or structural homologies with chromosomal DHFRs. The gene for this enzyme produces subunits that are 78 amino acids long, which assemble into a homotetramer possessing 222 symmetry. More recently, a tandem array of four gene copies linked in-frame was constructed, which produces a monomer containing 312 amino acids named Quad3. Asymmetric mutations in Quad3 have also been constructed to probe the role of Q67 and K32 residues in catalysis. This present study mixes and matches mutations to determine if the Q67H mutation, which tightens binding approximately 100-fold to both dihydrofolate (DHF) and NADPH, can help rescue the K32M mutation. While the latter mutation weakens DHF binding over 60-fold, it concurrently increases kcat by a factor of 5. Two Q67H mutations were added to gene copies 1 and 4 in conjunction with the K32M mutation in gene copies 1 and 3. Addition of these Q67H mutations tightens binding 40-fold, and the catalytic efficiency (kcat/Km(DHF)) of the resulting protein is similar to that of Quad3. Since these Q67H mutations can mostly compensate for the K32M lesion, K32 must not be necessary for DHF binding. Another multimutant combines the K32M mutation in gene copies 1 and 3 with the Q67H mutation in all gene copies. This mutant is inhibited by DHF but not NADPH, indicating that NADPH binds only to the wild type half of the pore, while DHF can bind to either the wild type or mutant half of the pore. This inhibition pattern contrasts with the mutant containing only the Q67H substitution in all four gene copies, which is severely inhibited by both NADPH and substrate. Since gene duplication and divergence are evolutionary tools for gaining function, these constructs are a first step toward building preferences for NADPH and DHF in each half of the active site pore of this primitive enzyme.

A Trimetal Site and Substrate Distortion in a Family II Inorganic Pyrophosphatase

We report the first crystal structures of a family II pyrophosphatase complexed with a substrate analogue, imidodiphosphate (PNP). These provide new insights into the catalytic reaction mechanism of this enzyme family. We were able to capture the substrate complex both by fluoride inhibition and by site-directed mutagenesis providing complementary snapshots of the Michaelis complex. Structures of both the fluoride-inhibited wild type and the H98Q variant of the PNP-Bacillus subtilis pyrophosphatase complex show a unique trinuclear metal center. Each metal ion coordinates a terminal oxygen on the electrophilic phosphate and a lone pair on the putative nucleophile, thus placing it in line with the scissile bond without any coordination by protein. The nucleophile moves further away from the electrophilic phosphorus site, to the opposite side of the trimetal plane, upon binding of substrate. In comparison with earlier product complexes, the side chain of Lys296 has swung in and so three positively charged side chains, His98, Lys205 and Lys296, now surround the bridging nitrogen in PNP. Finally, one of the active sites in the wild-type structure appears to show evidence of substrate distortion. Binding to the enzyme may thus strain the substrate and thus enhance the catalytic rate.

Neutron diffraction studies of Escherichia coli dihydrofolate reductase complexed with methotrexate

Hydrogen atoms play a central role in many biochemical processes yet are difficult to visualize by x-ray crystallography. Spallation neutron sources provide a new arena for protein crystallography with TOF measurements enhancing data collection efficiency and allowing hydrogen atoms to be located in smaller crystals of larger biological macromolecules. Here we report a 2.2-A resolution neutron structure of Escherichia coli dihydrofolate reductase (DHFR) in complex with methotrexate (MTX). Neutron data were collected on a 0.3-mm(3) D(2)O-soaked crystal at the Los Alamos Neutron Scattering Center. This study provides an example of using spallation neutrons to study protein dynamics, to identify protonation states directly from nuclear density maps, and to analyze solvent structure. Our structure reveals that the occluded loop conformation [monomer (mon.) A] of the DHFR.MTX complex undergoes greater H/D exchange compared with the closed-loop conformer (mon. B), partly because the Met-20 and beta(F-G) loops readily exchange in mon. A. The eight-stranded beta sheet of both DHFR molecules resists H/D exchange more than the helices and loops. However, the C-terminal strand, betaH, in mon. A is almost fully exchanged. Several D(2)Os form hydrogen bonds with exchanged amides. At the active site, the N1 atom of MTX is protonated and thus charged when bound to DHFR. Several D(2)Os are observed at hydrophobic surfaces, including two pockets near the MTX-binding site. A previously unidentified D(2)O hydrogen bonds with the catalytic D27 in mon. B, stabilizing its negative charge.

Relating Protein Motion to Catalysis

This review examines the linkage between protein conformational motions and enzyme catalysis. The fundamental issues related to this linkage are probed in the context of two enzymes that catalyze hydride transfer, namely dihydrofolate reductase and liver alcohol dehydrogenase. The extensive experimental and theoretical studies addressing the role of protein conformational changes in these enzyme reactions are summarized. Evidence is presented for a network of coupled motions throughout the protein fold that facilitate the chemical reaction. This network is comprised of fast thermal motions that are in equilibrium as the reaction progresses along the reaction coordinate and that lead to slower equilibrium conformational changes conducive to the chemical reaction.

Ensemble-based signatures of energy propagation in proteins: a new view of an old phenomenon

The ability of a protein to transmit the energetic effects of binding from one site to another constitutes the underlying basis for allosterism and signal transduction. Despite clear experimental evidence indicating the ability of proteins to transmit the effects of binding, the means by which this propagation is facilitated is not well understood. Using our previously developed ensemble-based description of the equilibrium, we investigated the physical basis of energy propagation and identified several fundamental and general aspects of energetic coupling between residues in a protein. First, partitioning of a conformational ensemble into four distinct sub-ensembles allows for explanation of the range of experimentally observed coupling behaviors (i.e., positive, neutral, and negative coupling between various regions of the protein structure). Second, the relative thermodynamic properties of these four sub-ensembles define the energetic coupling between residues as either positive, neutral, or negative. Third, analysis of the structural and thermodynamic features of the states within each sub-ensemble reveals significant variability. This third result suggests that a quantitative description of energy propagation in proteins requires an understanding of the structural and energetic features of more than just one or a few low-energy states, but also of many high-energy states. Such findings illuminate the difficulty in interpreting energy propagation in proteins in terms of a structural pathway that physically links coupled sites.

Structure and hydride transfer mechanism of a moderate thermophilic dihydrofolate reductase from Bacillus stearothermophilus and comparison to its mesophilic and hyperthermophilic homologues

Dihydrofolate reductase (DHFR) from a moderate thermophilic organism, Bacillus stearothermophilus, has been cloned and expressed. Physical characterization of the protein (BsDHFR) indicates that it is a monomeric protein with a molecular mass of 18,694.6 Da (0.8), coincident with the mass of 18 694.67 Da calculated from the primary sequence. Determination of the X-ray structure of BsDHFR provides the first structure for a monomeric DHFR from a thermophilic organism, indicating a high degree of conservation of structure in relation to all chromosomal DHFRs. Structurally based sequence alignment of DHFRs indicates the following levels of sequence identity and similarity for BsDHFR: 38 and 58% with Escherichia coli, 35 and 56% with Lactobacillus casei, and 23 and 40% with Thermotoga maritima, respectively. Steady state kinetic isotope effect studies indicate an ordered kinetic mechanism at elevated temperatures, with NADPH binding first to the enzyme. This converts to a more random mechanism at reduced temperatures, reflected in a greatly reduced K(m) for dihydrofolate at 20 degrees C in relation to that at 60 degrees C. A reduction in either temperature or pH reduces the degree to which the hydride transfer step is rate-determining for the second-order reaction of DHF with the enzyme-NADPH binary complex. Transient state kinetics have been used to study the temperature dependence of the isotope effect on hydride transfer at pH 9 between 10 and 50 degrees C. The data support rate-limiting hydride transfer with a moderate enthalpy of activation (E(a) = 5.5 kcal/mol) and a somewhat greater temperature dependence for the kinetic isotope effect than predicted from classical behavior [A(H)/A(D) = 0.57 (0.15)]. Comparison of kinetic parameters for BsDHFR to published data for DHFR from E. coli and T. maritima shows a decreasing trend in efficiency of hydride transfer with increasing thermophilicity of the protein. These results are discussed in the context of the capacity of each enzyme to optimize H-tunneling from donor (NADPH) to acceptor (DHF) substrates.

Dynamics of immobilized and native Escherichia coli dihydrofolate reductase by quasielastic neutron scattering

The internal dynamics of native and immobilized Escherichia coli dihydrofolate reductase (DHFR) have been examined using incoherent quasielastic neutron scattering. These results reveal no difference between the high frequency vibration mean-square displacement of the native and the immobilized E. coli DHFR. However, length-scale-dependent, picosecond dynamical changes are found. On longer length scales, the dynamics are comparable for both DHFR samples. On shorter length scales, the dynamics is dominated by local jump motions over potential barriers. The residence time for the protons to stay in a potential well is tau = 7.95 +/- 1.02 ps for the native DHFR and tau = 20.36 +/- 1.80 ps for the immobilized DHFR. The average height of the potential barrier to the local motions is increased in the immobilized DHFR, and may increase the activation energy for the activity reaction, decreasing the rate as observed experimentally. These results suggest that the local motions on the picosecond timescale may act as a lubricant for those associated with DHFR activity occurring on a slower millisecond timescale. Experiments indicate a significantly slower catalytic reaction rate for the immobilized E. coli DHFR. However, the immobilization of the DHFR is on the exterior of the enzyme and essentially distal to the active site, thus this phenomenon has broad implications for the action of drugs distal to the active site.

Searching sequence space: two different approaches to dihydrofolate reductase catalysis

There are numerous examples of proteins that catalyze the same reaction while possessing different structures. This review focuses on two dihydrofolate reductases (DHFRs) that have disparate structures and discusses how the catalytic strategies of these two DHFRs are driven by their respective scaffolds. The two proteins are E. coli chromosomal DHFR (Ec DHFR) and a type II R-plasmid-encoded DHFR, typified by R67 DHFR. The former has been described as a very well evolved enzyme with an efficiency of 0.15, while the latter has been suggested to be a model for a "primitive" enzyme that has not yet been optimized by evolution. This comparison underlines what is important to catalysis in these two enzymes and concurrently highlights fundamental issues in enzyme catalysis.

Defining the role of active-site loop fluctuations in dihydrofolate reductase catalysis

Dynamic processes are implicit in the catalytic function of all enzymes. To obtain insights into the relationship between the dynamics and thermodynamics of protein fluctuations and catalysis, we have measured millisecond time scale motions in the enzyme dihydrofolate reductase using NMR relaxation methods. Studies of a ternary complex formed from the substrate analog folate and oxidized NADP+ cofactor revealed conformational exchange between a ground state, in which the active site loops adopt a closed conformation, and a weakly populated (4.2% at 30 degrees C) excited state with the loops in the occluded conformation. Fluctuations between these states, which involve motions of the nicotinamide ring of the cofactor into and out of the active site, occur on a time scale that is directly relevant to the structural transitions involved in progression through the catalytic cycle.

Defining the binding site of homotetrameric R67 dihydrofolate reductase and correlating binding enthalpy with catalysis

R67 dihydrofolate reductase (DHFR) is a novel protein that possesses 222 symmetry. A single active site pore traverses the length of the homotetramer. Although the 222 symmetry implies that four symmetry-related binding sites should exist for each substrate as well as each cofactor, isothermal titration calorimetry (ITC) studies indicate only two molecules bind. Three possible combinations include two dihydrofolate molecules, two NADPH molecules, or one substrate with one cofactor. The latter is the productive ternary complex. To evaluate the roles of A36, Y46, T51, G64, and V66 residues in binding and catalysis, a site-directed mutagenesis approach was employed. One mutation per gene produces four mutations per active site pore, which often result in large cumulative effects. Conservative mutations at these positions either eliminate the ability of the gene to confer trimethoprim resistance or have no effect on catalysis. This result, in conjunction with previous mutagenesis studies on K32, K33, S65, Q67, I68, and Y69 [Strader, M. B., et al. (2001) Biochemistry 40, 11344-11352; Hicks, S. N., et al. (2003) Biochemistry 42, 10569-10578; Park, H., et al. (1997) Protein Eng. 10, 1415-1424], allows mapping of the active site surface. Residues for which conservative mutations have large effects on binding and catalysis include K32, Q67, I68, and Y69. These residues form a stripe that establishes the ligand binding surface. Residues that accommodate conservative mutations that do not greatly affect catalysis include K33, Y46, T51, S65, and V66. Isothermal titration calorimetry studies were also conducted on many of the mutants described above to determine the enthalpy of folate binding to the R67 DHFR.NADPH complex. A linear correlation between this DeltaH value and log k(cat)/K(m) is observed. Since structural tightness appears to be correlated with the exothermicity of the binding interaction, this leads to the hypothesis that enthalpy-driven formation of the ternary complex in these R67 DHFR variants plays a strong role in catalysis. Use of the alternate cofactor, NADH, extends this correlation, indicating preorganization of the ternary complex determines the efficiency of the reaction. This hypothesis is consistent with data suggesting R67 DHFR uses an endo transition state (where the nicotinamide ring of cofactor overlaps the more bulky side of the substrate's pteridine ring).

Microsecond freeze-hyperquenching: development of a new ultrafast micro-mixing and sampling technology and application to enzyme catalysis

A novel freeze-quench instrument with a characteristic <<dead-time>> of 137 +/- 18 micros is reported. The prototype has several key features that distinguish it from conventional freeze-quench devices and provide a significant improvement in time resolution: (a) high operating pressures (up to 400 bar) result in a sample flow with high linear rates (up to 200 m s(-1)); (b) tangential micro-mixer with an operating volume of approximately 1 nl yields short mixing times (up to 20 micros); (c) fast transport between the mixer and the cryomedium results in short reaction times: the ageing solution exits the mixer as a free-flowing jet, and the chemical reaction occurs "in-flight" on the way to the cryomedium; (d) a small jet diameter (approximately 20 microm) and a high jet velocity (approximately 200 m s(-1)) provide high sample-cooling rates, resulting in a short cryofixation time (up to 30 micros). The dynamic range of the freeze-quench device is between 130 micros and 15 ms. The novel tangential micro-mixer efficiently mixes viscous aqueous solutions, showing more than 95% mixing at eta < or = 4 (equivalent to protein concentrations up to 250 mg ml(-1)), which makes it an excellent tool for the preparation of pre-steady state samples of concentrated protein solutions for spectroscopic structure analysis. The novel freeze-quench device is characterized using the reaction of binding of azide to metmyoglobin from horse heart. Reaction samples are analyzed using 77 K optical absorbance spectroscopy, and X-band EPR spectroscopy. A simple procedure of spectral analysis is reported that allows (a) to perform a quantitative analysis of the reaction kinetics and (b) to identify and characterize novel reaction intermediates. The reduction of dioxygen by the bo3-type quinol oxidase from Escherichia coli is assayed using the MHQ technique. In these pilot experiments, low-temperature optical absorbance measurements show the rapid oxidation of heme o3 in the first 137 micros of the reaction, accompanied by the formation of an oxo-ferryl species. X-band EPR spectroscopy shows that a short-living radical intermediate is formed during the oxidation of heme o3. The radical decays within approximately 1 ms concomitant with the oxidation of heme b, and can be attributed to the PM reaction intermediate converting to the oxoferryl intermediate F. The general field of application of the freeze-quench methodology is discussed.

Structure, dynamics, and catalytic function of dihydrofolate reductase

Molecular motions are widely regarded as contributing factors in many aspects of protein function. The enzyme dihydrofolate reductase (DHFR), and particularly that from Escherichia coli, has become an important system for investigating the linkage between protein dynamics and catalytic function, both because of the location and timescales of the motions observed and because of the availability of a large amount of structural and mechanistic data that provides a detailed context within which the motions can be interpreted. Changes in protein dynamics in response to ligand binding, conformational change, and mutagenesis have been probed using numerous experimental and theoretical approaches, including X-ray crystallography, fluorescence, nuclear magnetic resonance (NMR), molecular dynamics simulations, and hybrid quantum/classical dynamics methods. These studies provide a detailed map of changes in conformation and dynamics throughout the catalytic cycle of DHFR and give new insights into the role of protein motions in the catalytic activity of this enzyme.

The active site of the ribosome is composed of two layers of conserved nucleotides with distinct roles in peptide bond formation and peptide release

Peptide bond formation and peptide release are catalyzed in the active site of the large subunit of the ribosome where universally conserved nucleotides surround the CCA ends of the peptidyl- and aminoacyl-tRNA substrates. Here, we describe the use of an affinity-tagging system for the purification of mutant ribosomes and analysis of four universally conserved nucleotides in the innermost layer of the active site: A2451, U2506, U2585, and A2602. While pre-steady-state kinetic analysis of the peptidyl transferase activity of the mutant ribosomes reveals substantially reduced rates of peptide bond formation using the minimal substrate puromycin, their rates of peptide bond formation are unaffected when the substrates are intact aminoacyl-tRNAs. These mutant ribosomes do, however, display substantial defects in peptide release. These results reveal a view of the catalytic center in which an inner shell of conserved nucleotides is pivotal for peptide release, while an outer shell is responsible for promoting peptide bond formation.

Reaction-path energetics and kinetics of the hydride transfer reaction catalyzed by dihydrofolate reductase

We have studied the hydride transfer reaction catalyzed by the enzyme dihydrofolate reductase (DHFR) and the coenzyme nicotinamide adenine dinucleotide phosphate (NADPH); the substrate is 5-protonated 7,8-dihydrofolate, and the product is tetrahydrofolate. The potential energy surface is modeled by a combined quantum mechanical-molecular mechanical (QM/MM) method employing Austin model 1 (AM1) and a simple valence bond potential for 69 QM atoms and employing the CHARMM22 and TIP3P molecular mechanics force fields for the other 21 399 atoms; the QM and MM regions are joined by two boundary atoms treated by the generalized hybrid orbital (GHO) method. All simulations are carried out using periodic boundary conditions at neutral pH and 298 K. In stage 1, a reaction coordinate is defined as the difference between the breaking and forming bond distances to the hydride ion, and a quasithermodynamic free energy profile is calculated along this reaction coordinate. This calculation includes quantization effects on bound vibrations but not on the reaction coordinate, and it is used to locate the variational transition state that defines a transition state ensemble. Then, the key interactions at the reactant, variational transition state, and product are analyzed in terms of both bond distances and electrostatic energies. The results of both analyses support the conclusion derived from previous mutational studies that the M20 loop of DHFR makes an important contribution to the electrostatic stabilization of the hydride transfer transition state. Third, transmission coefficients (including recrossing factors and multidimensional tunneling) are calculated and averaged over the transition state ensemble. These averaged transmission coefficients, combined with the quasithermodynamic free energy profile determined in stage 1, allow us to calculate rate constants, phenomenological free energies of activation, and primary and secondary kinetic isotope effects. A primary kinetic isotope effect (KIE) of 2.8 has been obtained, in good agreement with the experimentally determined value of 3.0 and with the value 3.2 calculated previously. The primary KIE is mainly a consequence of the quantization of bound vibrations. In contrast, the secondary KIE, with a value of 1.13, is almost entirely due to dynamical effects on the reaction coordinate, especially tunneling.