Relating Protein Motion to Catalysis

Cooperativity Principles in Self-Assembled Nanomedicine

Nanomedicine is a discipline that applies nanoscience and nanotechnology principles to the prevention, diagnosis, and treatment of human diseases. Self-assembly of molecular components is becoming a common strategy in the design and syntheses of nanomaterials for biomedical applications. In both natural and synthetic self-assembled nanostructures, molecular cooperativity is emerging as an important hallmark. In many cases, interplay of many types of noncovalent interactions leads to dynamic nanosystems with emergent properties where the whole is bigger than the sum of the parts. In this review, we provide a comprehensive analysis of the cooperativity principles in multiple self-assembled nanostructures. We discuss the molecular origin and quantitative modeling of cooperative behaviors. In selected systems, we describe the examples on how to leverage molecular cooperativity to design nanomedicine with improved diagnostic precision and therapeutic efficacy in medicine.

Direct observation of ultrafast large-scale dynamics of an enzyme under turnover conditions

The potential effect of conformational dynamics of enzymes on their chemical steps has been intensely debated recently. We use single-molecule FRET experiments on adenylate kinase (AK) to shed new light on this question. AK closes its domains to bring its two substrate close together for reaction. We show that domain closure takes only microseconds to complete, which is two orders of magnitude faster than the chemical reaction. Nevertheless, active-site mutants that reduce the rate of domain closure also reduce the reaction rate, suggesting a connection between the two phenomena. We propose that ultrafast domain closure is used by enzymes as a mechanism to optimize mutual orientation of substrates, a novel mode of coupling between conformational dynamics and catalysis.

The functional cycle of many proteins involves large-scale motions of domains and subunits. The relation between conformational dynamics and the chemical steps of enzymes remains under debate. Here we show that in the presence of substrates, domain motions of an enzyme can take place on the microsecond time scale, yet exert influence on the much-slower chemical step. We study the domain closure reaction of the enzyme adenylate kinase from Escherichia coli while in action (i.e., under turnover conditions), using single-molecule FRET spectroscopy. We find that substrate binding increases dramatically domain closing and opening times, making them as short as ∼15 and ∼45 µs, respectively. These large-scale conformational dynamics are likely the fastest measured to date, and are ∼100–200 times faster than the enzymatic turnover rate. Some active-site mutants are shown to fully or partially prevent the substrate-induced increase in domain closure times, while at the same time they also reduce enzymatic activity, establishing a clear connection between the two phenomena, despite their disparate time scales. Based on these surprising observations, we propose a paradigm for the mode of action of enzymes, in which numerous cycles of conformational rearrangement are required to find a mutual orientation of substrates that is optimal for the chemical reaction.

A reevaluation of the origin of the rate acceleration for enzyme-catalyzed hydride transfer

There is no consensus of opinion on the origin of the large rate accelerations observed for enzyme-catalyzed hydride transfer. The interpretation of recent results from studies on hydride transfer reactions catalyzed by alcohol dehydrogenase (ADH) focus on the proposal that the effective barrier height is reduced by quantum-mechanical tunneling through the energy barrier. This interpretation contrasts sharply with the notion that enzymatic rate accelerations are obtained through direct stabilization of the transition state for the nonenzymatic reaction in water. The binding energy of the dianion of substrate DHAP provides 11 kcal mol^-1 stabilization of the transition state for the hydride transfer reaction catalyzed by glycerol-3-phosphate dehydrogenase (GPDH). We summarize evidence that the binding interactions between (GPDH) and dianion activators are utilized directly for stabilization of the transition state for enzyme-catalyzed hydride transfer. The possibility is considered, and then discounted, that these dianion binding interactions are utilized for the stabilization of a tunnel ready state (TRS) that enables efficient tunneling of the transferred hydride through the energy barrier, and underneath the energy maximum for the transition state. It is noted that the evidence to support the existence of a tunnel-ready state for the hydride transfer reactions catalyzed by ADH is ambiguous. We propose that the rate acceleration for ADH is due to the utilization of the binding energy of the cofactor NAD+/NADH in the stabilization of the transition state for enzyme-catalyzed hydride transfer.

Genetically Fused T4L Acts as a Shield in Covalent Enzyme Immobilisation Enhancing the Rescued Activity

Enzyme immobilisation is a common strategy to increase enzymes resistance and reusability in a variety of excellent ‘green’ applications. However, the interaction with the solid support often leads to diminished specific activity, especially when non-specific covalent binding to the carrier takes place which affects the delicate architecture of the enzyme. Here we developed a broadly applicable strategy where the T4-lysozyme (T4L) is genetically fused at the N-terminus of different enzymes and used as inert protein spacer which directly attaches to the carrier preventing shape distortion of the catalyst. Halomonas elongata aminotransferase (HEWT), Bacillus subtilis engineered esterase (BS2m), and horse liver alcohol dehydrogenase (HLADH) were used as model enzymes to elucidate the benefits of the spacer. While HEWT and HLADH activity and expression were diminished by the fused T4L, both enzymes retained almost quantitative activity after immobilisation. In the case of BS2m, the protective effect of the T4L effectively was important and led to up to 10-fold improvement in the rescued activity.

Unified understanding of folding and binding mechanisms of globular and intrinsically disordered proteins

Extensive experimental and theoretical studies have advanced our understanding of the mechanisms of folding and binding of globular proteins, and coupled folding and binding of intrinsically disordered proteins (IDPs). The forces responsible for conformational changes and binding are common in both proteins; however, these mechanisms have been separately discussed. Here, we attempt to integrate the mechanisms of coupled folding and binding of IDPs, folding of small and multi-subdomain proteins, folding of multimeric proteins, and ligand binding of globular proteins in terms of conformational selection and induced-fit mechanisms as well as the nucleation–condensation mechanism that is intermediate between them. Accumulating evidence has shown that both the rate of conformational change and apparent rate of binding between interacting elements can determine reaction mechanisms. Coupled folding and binding of IDPs occurs mainly by induced-fit because of the slow folding in the free form, while ligand binding of globular proteins occurs mainly by conformational selection because of rapid conformational change. Protein folding can be regarded as the binding of intramolecular segments accompanied by secondary structure formation. Multi-subdomain proteins fold mainly by the induced-fit (hydrophobic collapse) mechanism, as the connection of interacting segments enhances the binding (compaction) rate. Fewer hydrophobic residues in small proteins reduce the intramolecular binding rate, resulting in the nucleation–condensation mechanism. Thus, the folding and binding of globular proteins and IDPs obey the same general principle, suggesting that the coarse-grained, statistical mechanical model of protein folding is promising for a unified theoretical description of all mechanisms.

Slow domain reconfiguration causes power-law kinetics in a two-state enzyme

Protein dynamics are typically captured well by rate equations that predict exponential decays for two-state reactions. Here, we describe a remarkable exception. The electron-transfer enzyme quiescin sulfhydryl oxidase (QSOX), a natural fusion of two functionally distinct domains, switches between open- and closed-domain arrangements with apparent power-law kinetics. Using single-molecule FRET experiments on time scales from nanoseconds to milliseconds, we show that the unusual open-close kinetics results from slow sampling of an ensemble of disordered domain orientations. While substrate accelerates the kinetics, thus suggesting a substrate-induced switch to an alternative free energy landscape of the enzyme, the power-law behavior is also preserved upon electron load. Our results show that the slow sampling of open conformers is caused by a variety of interdomain interactions that imply a rugged free energy landscape, thus providing a generic mechanism for dynamic disorder in multidomain enzymes.

Protein Mass Effects on Formate Dehydrogenase

Isotopically labeled enzymes (denoted as "heavy" or "Born-Oppenheimer" enzymes) have been used to test the role of protein dynamics in catalysis. The original idea was that the protein's higher mass would reduce the frequency of its normal-modes without altering its electrostatics. Heavy enzymes have been used to test if the vibrations in the native enzyme are coupled to the chemistry it catalyzes, and different studies have resulted in ambiguous findings. Here the temperature-dependence of intrinsic kinetic isotope effects of the enzyme formate dehydrogenase is used to examine the distribution of H-donor to H-acceptor distance as a function of the protein's mass. The protein dynamics are altered in the heavy enzyme to diminish motions that determine the transition state sampling in the native enzyme, in accordance with a Born-Oppenheimer-like effect on bond activation. Findings of this work suggest components related to fast frequencies that can be explained by Born-Oppenheimer enzyme hypothesis (vibrational) and also slower time scale events that are non-Born-Oppenheimer in nature (electrostatic), based on evaluations of protein mass dependence of donor-acceptor distance and forward commitment to catalysis along with steady state and single turnover measurements. Together, the findings suggest that the mass modulation affected both local, fast, protein vibrations associated with the catalyzed chemistry and the protein's macromolecular electrostatics at slower time scales; that is, both Born-Oppenheimer and non-Born-Oppenheimer effects are observed. Comparison to previous studies leads to the conclusion that isotopic labeling of the protein may have different effects on different systems, however, making heavy enzyme studies a very exciting technique for exploring the dynamics link to catalysis in proteins.

Patterns of coevolving amino acids unveil structural and dynamical domains

Patterns of interacting amino acids are so preserved within protein families that the sole analysis of evolutionary comutations can identify pairs of contacting residues. It is also known that evolution conserves functional dynamics, i.e., the concerted motion or displacement of large protein regions or domains. Is it, therefore, possible to use a pure sequence-based analysis to identify these dynamical domains? To address this question, we introduce here a general coevolutionary coupling analysis strategy and apply it to a curated sequence database of hundreds of protein families. For most families, the sequence-based method partitions amino acids into a few clusters. When viewed in the context of the native structure, these clusters have the signature characteristics of viable protein domains: They are spatially separated but individually compact. They have a direct functional bearing too, as shown for various reference cases. We conclude that even large-scale structural and functionally related properties can be recovered from inference methods applied to evolutionary-related sequences. The method introduced here is available as a software package and web server (spectrus.sissa.it/spectrus-evo_webserver).

Mechanism of salt-induced activity enhancement of a marine-derived laccase, Lac15

Laccase (benzenediol: oxygen oxidoreductases, EC1.10.3.2) is a multi-copper oxidase capable of oxidizing a variety of phenolic and other aromatic organic compounds. The catalytic power of laccase makes it an attractive candidate for potential applications in many areas of industry including biodegradation of organic pollutants and synthesis of novel drugs. Most laccases are vulnerable to high salt and have limited applications. However, some laccases are not only tolerant to but also activated by certain concentrations of salt and thus have great application potential. The mechanisms of salt-induced activity enhancement of laccases are unclear as yet. In this study, we used dynamic light scattering, size exclusion chromatography, analytical ultracentrifugation, intrinsic fluorescence emission, circular dichroism, ultraviolet-visible light absorption, and an enzymatic assay to investigate the potential correlation between the structure and activity of the marine-derived laccase, Lac15, whose activity is promoted by low concentrations of NaCl. The results showed that low concentrations of NaCl exert little influence on the protein structure, which was partially folded in the absence of the salt; moreover, the partially folded rather than the fully folded state seemed to be favorable for enzyme activity, and this partially folded state was distinctive from the so-called 'molten globule' occasionally observed in active enzymes. More data indicated that salt might promote laccase activity through mechanisms involving perturbation of specific local sites rather than a change in global structure. Potential binding sites for chloride ions and their roles in enzyme activity promotion are proposed.

Watching Proteins Wiggle: Mapping Structures with Two-Dimensional Infrared Spectroscopy

Proteins exhibit structural fluctuations over decades of time scales. From the picosecond side chain motions to aggregates that form over the course of minutes, characterizing protein structure over these vast lengths of time is important to understanding their function. In the past 15 years, two-dimensional infrared spectroscopy (2D IR) has been established as a versatile tool that can uniquely probe proteins structures on many time scales. In this review, we present some of the basic principles behind 2D IR and show how they have, and can, impact the field of protein biophysics. We highlight experiments in which 2D IR spectroscopy has provided structural and dynamical data that would be difficult to obtain with more standard structural biology techniques. We also highlight technological developments in 2D IR that continue to expand the scope of scientific problems that can be accessed in the biomedical sciences.

In silico screening, molecular docking, and molecular dynamics studies of SNP-derived human P5CR mutants

Pyrroline-5-carboxylate reductase (P5CR) encoded by PYCR1 gene is a housekeeping enzyme that catalyzes the reduction of P5C to proline using NAD(P)H as the cofactor. In this study, we used in silico approaches to examine the role of nonsynonymous single-nucleotide polymorphisms in the PYCR1 gene and their putative functions in the pathogenesis of Cutis Laxa. Among the 348 identified SNPs, 15 were predicted to be potentially damaging by both SIFT and PolyPhen tools; of them two SNP-derived mutations, R119G and G206W, have been previously reported to correlate with Cutis Laxa. These two mutations were therefore selected to be mapped to the wild-type (WT) P5CR structure for further structural and functional analyses. The results of comparative computational analyses using I-Mutant and Autodock reveal reductions in both stability and cofactor binding affinity of these two mutants. Comparative molecular dynamics (MD) simulations were performed to evaluate the changes in dynamic properties of P5CR upon mutations. The results reveal that the two mutations enhance the rigidity of P5CR structure, especially that of cofactor binding site, which could result in decreased kinetics of cofactor entrance and egress. Comparison between the structural properties of the WT and mutants during MD simulations shows that the enhanced rigidity of mutants results most likely from the increased number of inter-atomic interactions and the decreased number of dynamic hydrogen bonds. Our study provides novel insight into the deleterious effects of the R119G and G206W mutations on P5CR, and sheds light on the mechanisms by which these mutations mediate Cutis Laxa.

Two Dynamical Regimes of the Substrate Radical Rearrangement Reaction in B12-Dependent Ethanolamine Ammonia-Lyase Resolve Contributions of Native Protein Configurations and Collective Configurational Fluctuations to Catalysis

The kinetics of the substrate radical rearrangement reaction step in B₁₂-dependent ethanolamine ammonia-lyase (EAL) from Salmonella typhimurium are measured over a 92 K temperature range. The observed first-order rate constants display a piecewise-continuous Arrhenius dependence, with linear regions over 295 → 220 K (monoexponential) and 214 → 203 K (biexponential) that are delineated by a kinetic bifurcation and kinks at 219 and 217 K, respectively. The results are interpreted by using a free energy landscape model and derived microscopic kinetic mechanism. The bifurcation and kink transitions correspond to the effective quenching of two distinct sets of native collective protein configurational fluctuations that (1) reconfigure the protein within the substrate radical free energy minimum, in a reaction-enabling step, and (2) create the protein configurations associated with the chemical step. Below 217 K, the substrate radical decay reaction persists. Increases in activation enthalpy and entropy of both the microscopic enabling and reaction steps indicate that this non-native reaction coordinate is conducted by local, incremental fluctuations. Continuity in the Arrhenius relations indicates that the same sets of protein groups and interactions mediate the rearrangement over the 295 to 203 K range, but with a repertoire of configurations below 217 K that is restricted, relative to the native configurations accessible above 219 K. The experimental features of a culled reaction step, first-order kinetic measurements, and wide room-to-cryogenic temperature range, allow the direct demonstration and kinetic characterization of protein dynamical contributions to the core adiabatic, bond-making/bond-breaking reaction in EAL.

Better than Membranes at the Origin of Life?

Organelles without membranes are found in all types of cells and typically contain RNA and protein. RNA and protein are the constituents of ribosomes, one of the most ancient cellular structures. It is reasonable to propose that organelles without membranes preceded protocells and other membrane-bound structures at the origins of life. Such membraneless organelles would be well sheltered in the spaces between mica sheets, which have many advantages as a site for the origins of life.

Flexible Proteins at the Origin of Life

Almost all modern proteins possess well-defined, relatively rigid scaffolds that provide structural preorganization for desired functions. Such scaffolds require the sufficient length of a polypeptide chain and extensive evolutionary optimization. How ancestral proteins attained functionality, even though they were most likely markedly smaller than their contemporary descendants, remains a major, unresolved question in the origin of life. On the basis of evidence from experiments and computer simulations, we argue that at least some of the earliest water-soluble and membrane proteins were markedly more flexible than their modern counterparts. As an example, we consider a small, evolved in vitro ligase, based on a novel architecture that may be the archetype of primordial enzymes. The protein does not contain a hydrophobic core or conventional elements of the secondary structure characteristic of modern water-soluble proteins, but instead is built of a flexible, catalytic loop supported by a small hydrophilic core containing zinc atoms. It appears that disorder in the polypeptide chain imparts robustness to mutations in the protein core. Simple ion channels, likely the earliest membrane protein assemblies, could also be quite flexible, but still retain their functionality, again in contrast to their modern descendants. This is demonstrated in the example of antiamoebin, which can serve as a useful model of small peptides forming ancestral ion channels. Common features of the earliest, functional protein architectures discussed here include not only their flexibility, but also a low level of evolutionary optimization and heterogeneity in amino acid composition and, possibly, the type of peptide bonds in the protein backbone.

Applications of NMR and computational methodologies to study protein dynamics

Overwhelming evidence now illustrates the defining role of atomic-scale protein flexibility in biological events such as allostery, cell signaling, and enzyme catalysis. Over the years, spin relaxation nuclear magnetic resonance (NMR) has provided significant insights on the structural motions occurring on multiple time frames over the course of a protein life span. The present review article aims to illustrate to the broader community how this technique continues to shape many areas of protein science and engineering, in addition to being an indispensable tool for studying atomic-scale motions and functional characterization. Continuing developments in underlying NMR technology alongside software and hardware developments for complementary computational approaches now enable methodologies to routinely provide spatial directionality and structural representations traditionally harder to achieve solely using NMR spectroscopy. In addition to its well-established role in structural elucidation, we present recent examples that illustrate the combined power of selective isotope labeling, relaxation dispersion experiments, chemical shift analyses, and computational approaches for the characterization of conformational sub-states in proteins and enzymes.

Synthetic, Switchable Enzymes

The construction of switchable, radiation-controlled, aptameric enzymes - "swenzymes" - is, in principle, feasible. We propose a strategy to make such catalysts from 2 (or more) aptamers each selected to bind specifically to one of the substrates in, for example, a 2-substrate reaction. Construction of a combinatorial library of candidate swenzymes entails selecting a set of a million aptamers that bind one substrate and a second set of a million aptamers that bind the second substrate; the aptamers in these sets are then linked pairwise by a linker, thus bringing together the substrates. In the presence of the substrates, some linked aptamer pairs catalyze the reaction when exposed to external energy in the form of a specific frequency of low-intensity, nonionizing electromagnetic or acoustic radiation. Such swenzymes are detected via a separate product-capturing aptamer that changes conformation on capturing the product; this altered conformation allows it (1) to bind to every potential swenzyme in its vicinity (thereby giving a higher probability of capture to the swenzymes that generate the product) and (2) to bind to a sequence on a magnetic bead (thereby permitting purification of the swenzyme plus product-capturing aptamer by precipitation). Attempts to implement the swenzyme strategy may help elucidate fundamental problems in enzyme catalysis.

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.

Shedding light on the extra thermal stability of thermophilic proteins

An entropic stabilization mechanism has recently gained attention and credibility as the physical ground for the extra thermal stability of globular proteins from thermophilic microorganisms. An empirical result, obtained from the analysis of thermodynamic data for a large set of proteins, strengthens the general reliability of the theoretical approach originally devised to rationalize the occurrence of cold denaturation [Graziano, PCCP 2014, 16, 21755-21767]. It is shown that this theoretical approach can readily account for the entropic stabilization mechanism. On decreasing the conformational entropy gain associated with denaturation, the thermal stability of a model globular protein increases markedly.

On the indirect relationship between protein dynamics and enzyme activity

The behaviors of simple thermal systems have been well studied in physical chemistry and the principles obtained from such studies have been applied to complex thermal systems, such as proteins and enzymes. But the simple application of such principles is questionable and may lead to mistakes under some circumstances. In enzymology, the transition state theory of chemical reactions has been accepted as a fundamental theory, but the role of protein dynamics in enzyme catalysis is controversial in the context of transition state theory. By studying behaviors of complex thermal systems, we have revised the Arrhenius equation and transition state theory and our model is validated in enzymology. Formally speaking, the revised Arrhenius equation is apparently similar to a conventional Arrhenius equation, but the physical meanings of its parameters differ from that of traditional forms in principle. Within this model, the role of protein dynamics in enzyme catalysis is well defined and quantified.

Structure-based analysis of Bacilli and plasmid dihydrofolate reductase evolution

Dihydrofolate reductase (DHFR), a key enzyme in tetrahydrofolate-mediated biosynthetic pathways, has a structural motif known to be highly conserved over a wide range of organisms. Given its critical role in purine and amino acid synthesis, DHFR is a well established therapeutic target for treating a wide range of prokaryotic and eukaryotic infections as well as certain types of cancer. Here we present a structural-based computer analysis of bacterial (Bacilli) and plasmid DHFR evolution. We generated a structure-based sequence alignment using 7 wild-type DHFR x-ray crystal structures obtained from the RCSB Protein Data Bank and 350 chromosomal and plasmid homology models we generated from sequences obtained from the NCBI Protein Database. We used these alignments to compare active site and non-active site conservation in terms of amino acid residues, secondary structure and amino acid residue class. With respect to amino acid sequences and residue classes, active-site positions in both plasmid and chromosomal DHFR are significantly more conserved than non-active site positions. Secondary structure conservation was similar for active site and non-active site positions. Plasmid-encoded DHFR proteins have greater degree of sequence and residue class conservation, particularly in sequence positions associated with a network of concerted protein motions, than chromosomal-encoded DHFR proteins. These structure-based were used to build DHFR specific phylogenetic trees from which evidence for horizontal gene transfer was identified.

Putting Humpty-Dumpty Together: Clustering the Functional Dynamics of Single Biomolecular Machines Such as the Spliceosome

The spliceosome is a biomolecular machine that, in all eukaryotes, accomplishes site-specific splicing of introns from precursor messenger RNAs (pre-mRNAs) with high fidelity. Operating at the nanometer scale, where inertia and friction have lost the dominant role they play in the macroscopic realm, the spliceosome is highly dynamic and assembles its active site around each pre-mRNA anew. To understand the structural dynamics underlying the molecular motors, clocks, and ratchets that achieve functional accuracy in the yeast spliceosome (a long-standing model system), we have developed single-molecule fluorescence resonance energy transfer (smFRET) approaches that report changes in intra- and intermolecular interactions in real time. Building on our work using hidden Markov models (HMMs) to extract kinetic and conformational state information from smFRET time trajectories, we recognized that HMM analysis of individual state transitions as independent stochastic events is insufficient for a biomolecular machine as complex as the spliceosome. In this chapter, we elaborate on the recently developed smFRET-based Single-Molecule Cluster Analysis (SiMCAn) that dissects the intricate conformational dynamics of a pre-mRNA through the splicing cycle in a model-free fashion. By leveraging hierarchical clustering techniques developed for Bioinformatics, SiMCAn efficiently analyzes large datasets to first identify common molecular behaviors. Through a second level of clustering based on the abundance of dynamic behaviors exhibited by defined functional intermediates that have been stalled by biochemical or genetic tools, SiMCAn then efficiently assigns pre-mRNA FRET states and transitions to specific splicing complexes, with the potential to find heretofore undescribed conformations. SiMCAn thus arises as a general tool to analyze dynamic cellular machines more broadly.

The cooperativity effect in the reaction of soluble quinoprotein (PQQ-containing) glucose dehydrogenase is not due to subunit interaction but to substrate-assisted catalysis

Soluble quinoprotein (PQQ-containing) glucose dehydrogenase (sGDH, EC 1.1.99.35) catalyzes the oxidation of β-d-glucose to d-glucono-δ-lactone. Although sGDH has many analytical applications, the relationship between activity and substrate concentration is not well established. Previous steady-state kinetic studies revealed a negative cooperativity effect which has recently been ascribed to subunit interaction. To investigate this conclusion, stopped-flow kinetic experiments were carried out on the reaction in which oxidized enzyme (E_ox ) was reduced with substrates to E_red . The appearance of E_red is observed to be preceded by formation of an intermediate enzyme form, Int, which is mono-exponentially formed from E_ox . However, the rate of conversion of Int into E_red depends hyperbolically on the concentration of substrate (leading to a 35-fold stimulation in the case of glucose). Evidence is provided that substrate not only binds to E_ox but also to Int and E_red as well, and that the binding to Int causes the significant stimulation of Int decay. It is proposed that a proton shuffling step is involved in the decay, which is facilitated by binding of substrate to Int. Substituting the PQQ-activating Ca by a Ba ion lowered all reaction rates but did not change the stimulation factor. In summary, the previous proposal that the cooperativity effect of sGDH is due to interaction between its substrate-loaded subunits is incorrect; it is due to substrate-assisted catalysis of the enzyme.

ENZYMES

EC 1.1.99.35 - soluble quinoprotein glucose dehydrogenase.

Super RLuc8: A novel engineered Renilla luciferase with a red-shifted spectrum and stable light emission

Renilla luciferase is a bioluminescent enzyme which is broadly used as a reporter protein in molecular biosensors. In this study, a novel luciferase with desired light emission wavelength and thermostability is reported. The results indicated that the new luciferase, namely super RLuc8, had a red-shifted spectrum and showed stable light emission. Super RLuc8 showed a 10-fold (p-value=0.0084) increase in the thermostability at 37°C after 20min incubation, in comparison to the native enzyme. The optimum temperature of the mutant increased from 30 to 37°C. Molecular dynamics simulation analysis indicated that the increased thermostability was most probably caused by a better structural compactness and more local rigidity in the regions out of the emitter site.

Fold and flexibility: what can proteins' mechanical properties tell us about their folding nucleus?

The determination of a protein's folding nucleus, i.e. a set of native contacts playing an important role during its folding process, remains an elusive yet essential problem in biochemistry. In this work, we investigate the mechanical properties of 70 protein structures belonging to 14 protein families presenting various folds using coarse-grain Brownian dynamics simulations. The resulting rigidity profiles combined with multiple sequence alignments show that a limited set of rigid residues, which we call the consensus nucleus, occupy conserved positions along the protein sequence. These residues' side chains form a tight interaction network within the protein's core, thus making our consensus nuclei potential folding nuclei. A review of experimental and theoretical literature shows that most (above 80%) of these residues were indeed identified as folding nucleus member in earlier studies.

The Sampling of Conformational Dynamics in Ambient-Temperature Crystal Structures of Arginine Kinase

Arginine kinase provides a model for functional dynamics, studied through crystallography, enzymology, and nuclear magnetic resonance. Structures are now solved, at ambient temperature, for the transition state analog (TSA) complex. Analysis of quasi-rigid sub-domain displacements show that differences between the two TSA structures average about 5% of changes between substrate-free and TSA forms, and they are nearly co-linear. Small backbone hinge rotations map to sites that also flex on substrate binding. Anisotropic atomic displacement parameters (ADPs) are refined using rigid-body TLS constraints. Consistency between crystal forms shows that they reflect intrinsic molecular properties more than crystal lattice effects. In many regions, the favored directions of thermal/static displacement are appreciably correlated with movements on substrate binding. Correlation between ADPs and larger substrate-associated movements implies that the latter approximately follow paths of low-energy intrinsic motions.

Hydrostatic Pressure Studies Distinguish Global from Local Protein Motions in C-H Activation by Soybean Lipoxygenase-1

The proposed contributions of distinct classes of local versus global protein motions during enzymatic bond making/breaking processes has been difficult to verify. We employed soybean lipoxygenase-1 as a model system to investigate the impact of high pressure at variable temperatures on the hydrogen-tunneling properties of the wild-type protein and three single-site mutants. For all variants, pressure dramatically elevates the enthalpies of activation for the C-H activation. In contrast, the primary kinetic isotope effects (KIEs) for C-H activation and their corresponding temperature dependencies remain unchanged up to ca. 700 bar. The differential impact of elevated hydrostatic pressure on the temperature dependencies of rate constants versus substrate KIEs provides direct evidence for two distinct classes of protein motions: local, isotope-dependent donor-acceptor distance-sampling modes, and a more global, isotope-independent search for productive protein conformational sub-states.

Assessing Spectral Simulation Protocols for the Amide I Band of Proteins

We present a benchmark study of spectral simulation protocols for the amide I band of proteins. The amide I band is widely used in infrared spectroscopy of proteins due to the large signal intensity, high sensitivity to hydrogen bonding, and secondary structural motifs. This band has, thus, proven valuable in many studies of protein structure-function relationships. We benchmark spectral simulation protocols using two common force fields in combination with several electrostatic mappings and coupling models. The results are validated against experimental linear absorption and two-dimensional infrared spectroscopy for three well-studied proteins. We find two-dimensional infrared spectroscopy to be much more sensitive to the simulation protocol than linear absorption and report on the best simulation protocols. The findings demonstrate that there is still room for ideas to improve the existing models for the amide I band of proteins.

Examinations of the Chemical Step in Enzyme Catalysis

Advances in computational and experimental methods in enzymology have aided comprehension of enzyme-catalyzed chemical reactions. The main difficulty in comparing computational findings to rate measurements is that the first examines a single energy barrier, while the second frequently reflects a combination of many microscopic barriers. We present here intrinsic kinetic isotope effects and their temperature dependence as a useful experimental probe of a single chemical step in a complex kinetic cascade. Computational predictions are tested by this method for two model enzymes: dihydrofolate reductase and thymidylate synthase. The description highlights the significance of collaboration between experimentalists and theoreticians to develop a better understanding of enzyme-catalyzed chemical conversions.

Conformational flexibility related to enzyme activity: evidence for a dynamic active-site gatekeeper function of Tyr(215) in Aerococcus viridans lactate oxidase

L-Lactate oxidase (LOX) belongs to a large family of flavoenzymes that catalyze oxidation of α-hydroxy acids. How in these enzymes the protein structure controls reactivity presents an important but elusive problem. LOX contains a prominent tyrosine in the substrate binding pocket (Tyr(215) in Aerococcus viridans LOX) that is partially responsible for securing a flexible loop which sequesters the active site. To characterize the role of Tyr(215), effects of substitutions of the tyrosine (Y215F, Y215H) were analyzed kinetically, crystallographically and by molecular dynamics simulations. Enzyme variants showed slowed flavin reduction and oxidation by up to 33-fold. Pyruvate release was also decelerated and in Y215F, it was the slowest step overall. A 2.6-Å crystal structure of Y215F in complex with pyruvate shows the hydrogen bond between the phenolic hydroxyl and the keto oxygen in pyruvate is replaced with a potentially stronger hydrophobic interaction between the phenylalanine and the methyl group of pyruvate. Residues 200 through 215 or 216 appear to be disordered in two of the eight monomers in the asymmetric unit suggesting that they function as a lid controlling substrate entry and product exit from the active site. Substitutions of Tyr(215) can thus lead to a kinetic bottleneck in product release.

Solution NMR Spectroscopy for the Study of Enzyme Allostery

Allostery is a ubiquitous biological regulatory process in which distant binding sites within a protein or enzyme are functionally and thermodynamically coupled. Allosteric interactions play essential roles in many enzymological mechanisms, often facilitating formation of enzyme-substrate complexes and/or product release. Thus, elucidating the forces that drive allostery is critical to understanding the complex transformations of biomolecules. Currently, a number of models exist to describe allosteric behavior, taking into account energetics as well as conformational rearrangements and fluctuations. In the following Review, we discuss the use of solution NMR techniques designed to probe allosteric mechanisms in enzymes. NMR spectroscopy is unequaled in its ability to detect structural and dynamical changes in biomolecules, and the case studies presented herein demonstrate the range of insights to be gained from this valuable method. We also provide a detailed technical discussion of several specialized NMR experiments that are ideally suited for the study of enzymatic allostery.

Protein Conformational Changes Are Detected and Resolved Site Specifically by Second-Harmonic Generation

We present here a straightforward, broadly applicable technique for real-time detection and measurement of protein conformational changes in solution. This method is based on tethering proteins labeled with a second-harmonic generation (SHG) active dye to supported lipid bilayers. We demonstrate our method by measuring the conformational changes that occur upon ligand binding with three well-characterized proteins labeled at lysine residues: calmodulin (CaM), maltose-binding protein (MBP), and dihydrofolate reductase (DHFR). We also create a single-site cysteine mutant of DHFR engineered within the Met20 catalytic loop region and study the protein’s structural motion at this site. Using published x-ray crystal structures, we show that the changes in the SHG signals upon ligand binding are the result of structural motions that occur at the labeled sites between the apo and ligand-bound forms of the proteins, which are easily distinguished from each other. In addition, we demonstrate that different magnitudes of the SHG signal changes are due to different and specific ligand-induced conformational changes. Taken together, these data illustrate the potential of the SHG approach for detecting and measuring protein conformational changes for a wide range of biological applications.

Nonequilibrium Relaxation of Conformational Dynamics Facilitates Catalytic Reaction in an Elastic Network Model of T7 DNA Polymerase

Nucleotide-induced conformational closing of the finger domain of DNA polymerase is crucial for its catalytic action during DNA replication. Such large-amplitude molecular motion is often not fully accessible to either direct experimental monitoring or molecular dynamics simulations. However, a coarse-grained model can offer an informative alternative, especially for probing the relationship between conformational dynamics and catalysis. Here we investigate the dynamics of T7 DNA polymerase catalysis using a Langevin-type elastic network model incorporating detailed structural information on the open conformation without the substrate bound. Such a single-parameter model remarkably captures the induced conformational dynamics of DNA polymerase upon dNTP binding, and reveals its close coupling to the advancement toward transition state along the coordinate of the target reaction, which contributes to significant lowering of the activation energy barrier. Furthermore, analysis of stochastic catalytic rates suggests that when the activation energy barrier has already been significantly lowered and nonequilibrium relaxation toward the closed form dominates the catalytic rate, one must appeal to a picture of two-dimensional free energy surface in order to account for the full spectrum of catalytic modes. Our semiquantitative study illustrates the general role of conformational dynamics in achieving transition-state stabilization, and suggests that such an elastic network model, albeit simplified, possesses the potential to furnish significant mechanistic insights into the functioning of a variety of enzymatic systems.

Evolutionary Conserved Positions Define Protein Conformational Diversity

Conformational diversity of the native state plays a central role in modulating protein function. The selection paradigm sustains that different ligands shift the conformational equilibrium through their binding to highest-affinity conformers. Intramolecular vibrational dynamics associated to each conformation should guarantee conformational transitions, which due to its importance, could possibly be associated with evolutionary conserved traits. Normal mode analysis, based on a coarse-grained model of the protein, can provide the required information to explore these features. Herein, we present a novel procedure to identify key positions sustaining the conformational diversity associated to ligand binding. The method is applied to an adequate refined dataset of 188 paired protein structures in their bound and unbound forms. Firstly, normal modes most involved in the conformational change are selected according to their corresponding overlap with structural distortions introduced by ligand binding. The subspace defined by these modes is used to analyze the effect of simulated point mutations on preserving the conformational diversity of the protein. We find a negative correlation between the effects of mutations on these normal mode subspaces associated to ligand-binding and position-specific evolutionary conservations obtained from multiple sequence-structure alignments. Positions whose mutations are found to alter the most these subspaces are defined as key positions, that is, dynamically important residues that mediate the ligand-binding conformational change. These positions are shown to be evolutionary conserved, mostly buried aliphatic residues localized in regular structural regions of the protein like β-sheets and α-helix.

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.

Tuned by metals: the TET peptidase activity is controlled by 3 metal binding sites

TET aminopeptidases are dodecameric particles shared in the three life domains involved in various biological processes, from carbon source provider in archaea to eye-pressure regulation in humans. Each subunit contains a dinuclear metal site (M1 and M2) responsible for the enzyme catalytic activity. However, the role of each metal ion is still uncharacterized. Noteworthy, while mesophilic TETs are activated by Mn(2+), hyperthermophilic TETs prefers Co(2+). Here, by means of anomalous x-ray crystallography and enzyme kinetics measurements of the TET3 aminopeptidase from the hyperthermophilic organism Pyrococcus furiosus (PfTET3), we show that M2 hosts the catalytic activity of the enzyme, while M1 stabilizes the TET3 quaternary structure and controls the active site flexibility in a temperature dependent manner. A new third metal site (M3) was found in the substrate binding pocket, modulating the PfTET3 substrate preferences. These data show that TET activity is tuned by the molecular interplay among three metal sites.