An analytical method for determining relative specificities for sequential reactions catalyzed by the same enzyme: application to the hydrolysis of triacylglycerols by lipases

We propose a model for the sequential hydrolysis of ester bonds of triacylglycerols by lipases and use it as the basis for an analytical method for determining the relative specificity of the lipase for the various substrates with which it can react, when the substrates occur simultaneously in a single reaction system. We then apply the method to our own data and literature data involving the hydrolysis of triacylglycerols by lipases. Our model is able to fit well to most of the reaction profiles, enabling the estimation of relative specificities. We discuss the limitations and potential applications of our method.     

Phosphorylation-induced conformational changes in a mitogen-activated protein kinase substrate. Implications for tyrosine hydroxylase activation

Mitogen-activated protein (MAP) kinase-mediated phosphorylation of specific residues in tyrosine hydroxylase leads to an increase in enzyme activity. However, the mechanism whereby phosphorylation affects enzyme turnover is not well understood. We used a combination of fluorescence resonance energy transfer (FRET) measurements and molecular dynamics simulations to explore the conformational free energy landscape of a 10-residue MAP kinase substrate found near the N terminus of the enzyme. This region is believed to be part of an autoregulatory sequence that overlies the active site of the enzyme. FRET was used to measure the effect of phosphorylation on the ensemble of peptide conformations, and molecular dynamics simulations generated free energy profiles for both the unphosphorylated and phosphorylated peptides. We demonstrate how FRET transfer efficiencies can be calculated from molecular dynamics simulations. For both the unphosphorylated and phosphorylated peptides, the calculated FRET efficiencies are in excellent agreement with the experimentally determined values. Moreover, the FRET measurements and molecular simulations suggest that phosphorylation causes the peptide backbone to change direction and fold into a compact structure relative to the unphosphorylated state. These results are consistent with a model of enzyme activation where phosphorylation of the MAP kinase substrate causes the N-terminal region to adopt a compact structure away from the active site. The methods we employ provide a general framework for analyzing the accessible conformational states of peptides and small molecules. Therefore, they are expected to be applicable to a variety of different systems.     

Peptidases that terminate the action of enkephalins. Consideration of physiological importance for amino-, carboxy-, endo-, and pseudoenkephalinase

The term "enkephalinase" has been frequently applied to enzyme activity in a variety of tissue preparations. In some cases there has been the implication that cleavage of a specific peptide bond in the enkephalin molecule results from the action of a single enzyme with the major responsibility of inactivating synaptic enkephalin. It is not known to what extent diverse enkephalin-degrading enzymes, with differing peptide bond specificities, may act in concert at any given synapse. There do exist, however, enzymes having known characteristic specificities with respect both to peptide substrates, including enkephalins, and to identifiable peptide bonds. Thus, at any given site of enkephalin release there probably resides a characteristic assembly of peptidases concerned with inactivation of this neuromediator. We propose that the term "enkephalinase" be used to encompass the entire family of enkephalin-degrading enzymes, and that "aminoenkephalinase", "carboxyenkephalinase", "endoenkephalinase" and "pseudoenkephalinase" should designate enzymes of known specificities with respect to both peptide substrates and particular peptide bonds.     

Studies on the biochemical basis of spontaneous mutation. III. Rate model for DNA polymerase-effected nucleotide misincorporation

We propose a model to investigate the relation between insertion and excision activities of polymerases involved in DNA synthesis, and the frequency of errors resulting from substituting either mismatched bases or base analogues into a DNA molecule. An analytical equation is derived which expresses the error frequency as a function of nucleotide insertion and removal rates. For the general case, given arbitrary rates of insertion and removal, and allowing the enzyme to peel back by excising previously incorporated nucleotides, we have developed a computer simulation for the synthesis of a DNA molecule. In the special case, where insertion and removal frequencies are within the biologically interesting range for spontaneous mutations, the effect of “peelback” on error correction can be obtained analytically. Our results suggest that the magnitude of the removal frequency (3′-exonuclease activity) is the parameter that exerts the greatest influence on error correction capability; the frequency of errors is less sensitive to either the specificity for removal of mismatched relative to correctly matched bases, or to peelback.     

Identification of the sequence motif of glycoside hydrolase 13 family members

A bioinformatics analysis of sequences of enzymes of the glycoside hydrolase (GH) 13 family members such as α-amylase, cyclodextrin glycosyltransferase (CGTase), branching enzyme and cyclomaltodextrinase has been carried out in order to find out the sequence motifs that govern the reactions specificities of these enzymes by using hidden Markov model (HMM) profile. This analysis suggests the existence of such sequence motifs and residues of these motifs constituting the -1 to +3 catalytic subsites of the enzyme. Hence, by introducing mutations in the residues of these four subsites, one can change the reaction specificities of the enzymes. In general it has been observed that α -amylase sequence motif have low sequence conservation than rest of the motifs of the GH13 family members.     

The evolution of enzyme specificity in the metabolic replicator model of prebiotic evolution

The chemical machinery of life must have been catalytic from the outset. Models of the chemical origins have attempted to explain the ecological mechanisms maintaining a minimum necessary diversity of prebiotic replicator enzymes, but little attention has been paid so far to the evolutionary initiation of that diversity. We propose a possible first step in this direction: based on our previous model of a surface-bound metabolic replicator system we try to explain how the adaptive specialization of enzymatic replicator populations might have led to more diverse and more efficient communities of cooperating replicators with two different enzyme activities. The key assumptions of the model are that mutations in the replicator population can lead towards a) both of the two different enzyme specificities in separate replicators: efficient "specialists" or b) a "generalist" replicator type with both enzyme specificities working at less efficiency, or c) a fast-replicating, non-enzymatic "parasite". We show that under realistic trade-off constraints on the phenotypic effects of these mutations the evolved replicator community will be usually composed of both types of specialists and of a limited abundance of parasites, provided that the replicators can slowly migrate on the mineral surface. It is only at very weak trade-offs that generalists take over in a phase-transition-like manner. The parasites do not seriously harm the system but can freely mutate, therefore they can be considered as pre-adaptations to later, useful functions that the metabolic system can adopt to increase its own fitness.     

Amylo-1,6-glucosidase/1,4-α-glucan: 1,4-α-glucan 4-α-glycosyltransferase: Specificity toward polysaccharide substrates

The substrate specificities of the glucosidase-transferase debranching enzyme systems from yeast and rabbit muscle were examined by the use of polysaccharide substrates of defined outer chain lengths. The results were consistent with the specificities ascribed to the transferase portion of the debranching enzyme system by previous studies using maltosaccharide substrates. The specificities of the two enzyme systems were also examined in the reversion reaction. The results showed that both systems displayed an inverse specificity to that observed in the hydrolytic reaction. This suggested that the reversion reaction reflects the specificity of the glucosidase portion of the debranching system. The major differences between the specificities of the yeast and rabbit muscle systems were found to lie in the specificity of the transferase and in the ability of the yeast system to debranch native glycogen and amylopectin.     

Rational design of solution additives for the prevention of protein aggregation

We have developed a statistical-mechanical model of the effect of solution additives on protein association reactions. This model incorporates solvent radial distribution functions obtained from all-atom molecular dynamics simulations of particular proteins into simple models of protein interactions. In this way, the effects of additives can be computed along the entire association/dissociation reaction coordinate. We used the model to test our hypothesis that a class of large solution additives, which we term "neutral crowders," can slow protein association and dissociation by being preferentially excluded from protein-protein encounter complexes, in a manner analogous to osmotic stress. The magnitude of this proposed "gap effect" was probed for two simple model systems: the association of two spheres and the association of two planes. Our results suggest that for a protein of 20 A radius, an 8 A additive can increase the free energy barrier for association and dissociation by as much as 3-6 kcal/mol. Because the proposed gap effect is present only for reactions involving multiple molecules, it can be exploited to develop novel additives that affect protein association reactions although having little or no effect on unimolecular reactions such as protein folding. This idea has many potential applications in areas such as the stabilization of proteins against aggregation during folding and in pharmaceutical formulations.     

Sensitivity and control analysis of periodically forced reaction networks using the Green's function method

A general sensitivity and control analysis of periodically forced reaction networks with respect to small perturbations in arbitrary network parameters is presented. A well-known property of sensitivity coefficients for periodic processes in dynamical systems is that the coefficients generally become unbounded as time tends to infinity. To circumvent this conceptual obstacle, a relative time or phase variable is introduced so that the periodic sensitivity coefficients can be calculated. By employing the Green's function method, the sensitivity coefficients can be defined using integral control operators that relate small perturbations in the network's parameters and forcing frequency to variations in the metabolite concentrations and reaction fluxes. The properties of such operators do not depend on a particular parameter perturbation and are described by the summation and connectivity relationships within a control-matrix operator equation. The aim of this paper is to derive such a general control-matrix operator equation for periodically forced reaction networks, including metabolic pathways. To illustrate the general method, the two limiting cases of high and low forcing frequency are considered. We also discuss a practically important case where enzyme activities and forcing frequency are modulated simultaneously. We demonstrate the developed framework by calculating the sensitivity and control coefficients for a simple two reaction pathway where enzyme activities enter reaction rates linearly and specifically.     

Barrier Compression and Its Contribution to Both Classical and Quantum Mechanical Aspects of Enzyme Catalysis

It is generally accepted that enzymes catalyze reactions by lowering the apparent activation energy by transition state stabilization or through destabilization of ground states. A more controversial proposal is that enzymes can also accelerate reactions through barrier compression—an idea that has emerged from studies of H-tunneling reactions in enzyme systems. The effects of barrier compression on classical (over-the-barrier) reactions, and the partitioning between tunneling and classical reaction paths, have largely been ignored. We performed theoretical and computational studies on the effects of barrier compression on the shape of potential energy surfaces/reaction barriers for model (malonaldehyde and methane/methyl radical anion) and enzymatic (aromatic amine dehydrogenase) proton transfer systems. In all cases, we find that barrier compression is associated with an approximately linear decrease in the activation energy. For partially nonadiabatic proton transfers, we show that barrier compression enhances, to similar extents, the rate of classical and proton tunneling reactions. Our analysis suggests that barrier compression—through fast promoting vibrations, or other means—could be a general mechanism for enhancing the rate of not only tunneling, but also classical, proton transfers in enzyme catalysis.     

Combining Substrate Specificity Analysis with Support Vector Classifiers Reveals Feruloyl Esterase as a Phylogenetically Informative Protein Group

Background

Our understanding of how fungi evolved to develop a variety of ecological niches, is limited but of fundamental biological importance. Specifically, the evolution of enzymes affects how well species can adapt to new environmental conditions. Feruloyl esterases (FAEs) are enzymes able to hydrolyze the ester bonds linking ferulic acid to plant cell wall polysaccharides. The diversity of substrate specificities found in the FAE family shows that this family is old enough to have experienced the emergence and loss of many activities.

Methodology/Principal Findings

In this study we evaluate the relative activity of FAEs against a variety of model substrates as a novel predictive tool for Ascomycota taxonomic classification. Our approach consists of two analytical steps; (1) an initial unsupervised analysis to cluster the FAEs substrate specificity data which were generated by cultivation of 34 Ascomycota strains and then an analysis of the produced enzyme cocktail against 10 substituted cinnamate and phenylalkanoate methyl esters, (2) a second, supervised analysis for training a predictor built on these substrate activities. By applying both linear and non-linear models we were able to correctly predict the taxonomic Class (∼86% correct classification), Order (∼88% correct classification) and Family (∼88% correct classification) that the 34 Ascomycota belong to, using the activity profiles of the FAEs.

Conclusion/Significance

The good correlation with the FAEs substrate specificities that we have defined via our phylogenetic analysis not only suggests that FAEs are phylogenetically informative proteins but it is also a considerable step towards improved FAEs functional prediction.     

CRISPR-STAT: an easy and reliable PCR-based method to evaluate target-specific sgRNA activity

CRISPR/Cas9 has emerged as a versatile genome-engineering tool that relies on a single guide RNA (sgRNA) and the Cas9 enzyme for genome editing. Simple, fast and economical methods to generate sgRNAs have made targeted mutagenesis routine in cultured cells, mice, zebrafish and other model systems. Pre-screening of sgRNAs for target efficacy is desirable both for successful mutagenesis and minimizing wasted animal husbandry on targets with poor activity. Here, we describe an easy, quick and cost-effective fluorescent polymerase chain reaction (PCR)-based method, CRISPR Somatic Tissue Activity Test (CRISPR-STAT), to determine target-specific efficiency of sgRNA. As a proof of principle, we validated our method using 28 sgRNAs with known and varied levels of germline transmission efficiency in zebrafish by analysis of their somatic activity in injected embryos. Our data revealed a strong positive correlation between the fluorescent PCR profiles of the injected embryos and the germline transmission efficiency. Furthermore, the assay was sensitive enough to evaluate multiplex gene targeting. This method is easy to implement by laboratories with access to a capillary sequencer. Although we validated the method using CRISPR/Cas9 and zebrafish, it can be applied to other model systems and other genome targeting nucleases.     

Kinetic behaviour of proenzymes activation in the presence of different inhibitors for both activating and activated enzymes

In the present paper, a kinetic analysis of a general model for proenzyme activation, where the activating enzyme and also the activated one are reversibly inhibited in two steps by two different inhibitors, has been performed. The cases in which both inhibitors are the same, or in which the inhibition is irreversible (only one or the two inhibition routes) are treated as particular cases of the general model. In addition, the kinetic behaviour of many other proenzyme activation systems involving inhibition, particular cases of the reaction scheme under study, can be obtained. The total number of particular cases for the general model under study is 370, so this approach offers to the scientific community working in limited proteolysis regulation for the first time a method based on general solutions which only needs to be specified to their concrete problem of zymogen activation. Finally, new adimensional parameters are introduced, allowing the knowledgement, in the case that any of the inhibition routes is irreversible, the relative weight of both activation and irreversible inhibition routes.     

Anesthetic inhibition of firefly luciferase, a protein model for general anesthesia, does not exhibit pressure reversal.

The surprising observation that pressures of the order of 150 atmospheres can restore consciousness to an anesthetized animal has long been central to theories of the molecular mechanisms underlying general anesthesia. We have constructed a high-pressure gas chamber to test for "pressure reversal" of the best available protein model of general anesthetic target sites: the pure enzyme firefly luciferase, which accounts extremely well for animal potencies (over a 100,000-fold range). We found no significant pressure reversal for a variety of anesthetics of differing size and polarity. It thus appears that either firefly luciferase is not an adequate model for general anesthetic target sites or that pressure and anesthetics act at different molecular sites in the central nervous system.     

Transient model of thermal deactivation of enzymes

The kinetics of enzyme deactivation provide useful insights on processes that determine the level of biological function of any enzyme. Photinus pyralis (firefly) luciferase is a convenient enzyme system for studying mechanisms and kinetics of enzyme deactivation, refolding, and denaturation caused by various external factors, physical or chemical by nature. In this report we present a study of luciferase deactivation caused by increased temperature (i.e., thermal deactivation). We found that deactivation occurs through a reversible intermediate state and can be described by a Transient model that includes active and reversibly inactive states. The model can be used as a general framework for analysis of complex, multiexponential transient kinetics that can be observed for some enzymes by reaction progression assays. In this study the Transient model has been used to develop an analytical model for studying a time course of luciferase deactivation. The model might be applicable toward enzymes in general and can be used to determine if the enzyme exposed to external factors, physical or chemical by nature, undergoes structural transformation consistent with thermal mechanisms of deactivation.     

A new conceptual framework for enzyme catalysis. Hydrogen tunnelling coupled to enzyme dynamics in flavoprotein and quinoprotein enzymes.

Recent years have witnessed high levels of activity in identifying enzyme systems that catalyse H-transfer by quantum tunneling. Rather than being restricted to a small number of specific enzymes as perceived initially, it has now become an accepted mechanism for H-transfer in a growing number of enzymes. Furthermore, H-tunneling is driven by the thermally induced dynamics of the enzyme. In some of those enzymes that break stable C-H bonds the reaction proceeds purely by quantum tunneling, without the need to partially ascend the barrier. Enzymes studied that fall into this category include the flavoprotein and quinoprotein amine dehydrogenases, which have proved to be excellent model systems. These enzymes have enabled us to study the relationship between barrier shape and reaction kinetics. This has involved studies with "slow" and "fast" substrates and enzymes impaired by mutagenesis. A number of key questions now remain, including the nature of the coupling between protein dynamics and quantum tunneling. The wide-ranging implications of quantum tunneling introduce a paradigm shift in the conceptual framework for enzyme catalysis, inhibition and design.     

Conversion of a cyclodextrin glucanotransferase into an alpha-amylase: assessment of directed evolution strategies

Glycoside hydrolase family 13 (GH13) members have evolved to possess various distinct reaction specificities despite the overall structural similarity. In this study we investigated the evolutionary input required to effeciently interchange these specificities and also compared the effectiveness of laboratory evolution techniques applied, i.e., error-prone PCR and saturation mutagenesis. Conversion of our model enzyme, cyclodextrin glucanotransferase (CGTase), into an alpha-amylase like hydrolytic enzyme by saturation mutagenesis close to the catalytic core yielded a triple mutant (A231V/F260W/F184Q) with the highest hydrolytic rate ever recorded for a CGTase, similar to that of a highly active alpha-amylase, while cyclodextrin production was virtually abolished. Screening of a much larger, error-prone PCR generated library yielded far less effective mutants. Our results demonstrate that it requires only three mutations to change CGTase reaction specificity into that of another GH13 enzyme. This suggests that GH13 members may have diversified by introduction of a limited number of mutations to the common ancestor, and that interconversion of reaction specificites may prove easier than previously thought.     

Biapenem Inactivation by B2 Metallo β-Lactamases: Energy Landscape of the Hydrolysis Reaction

Background

A general mechanism has been proposed for metallo β-lactamases (MβLs), in which deprotonation of a water molecule near the Zn ion(s) results in the formation of a hydroxide ion that attacks the carbonyl oxygen of the β-lactam ring. However, because of the absence of X-ray structures that show the exact position of the antibiotic in the reactant state (RS) it has been difficult to obtain a definitive validation of this mechanism.

Methodology/Principal Findings

We have employed a strategy to identify the RS, which does not rely on substrate docking and/or molecular dynamics. Starting from the X-ray structure of the enzyme:product complex (the product state, PS), a QM/MM scan was used to drive the reaction uphill from product back to reactant. Since in this process also the enzyme changes from PS to RS, we actually generate the enzyme:substrate complex from product and avoid the uncertainties associated with models of the reactant state. We used this strategy to study the reaction of biapenem hydrolysis by B2 MβL CphA. QM/MM simulations were carried out under 14 different ionization states of the active site, in order to generate potential energy surfaces (PESs) corresponding to a variety of possible reaction paths.

Conclusions/Significance

The calculations support a model for biapenem hydrolysis by CphA, in which the nucleophile that attacks the β-lactam ring is not the water molecule located in proximity of the active site Zn, but a second water molecule, hydrogen bonded to the first one, which is used up in the reaction, and thus is not visible in the X-ray structure of the enzyme:product complex.     

Detecting DNA Modifications from SMRT Sequencing Data by Modeling Sequence Context Dependence of Polymerase Kinetic

DNA modifications such as methylation and DNA damage can play critical regulatory roles in biological systems. Single molecule, real time (SMRT) sequencing technology generates DNA sequences as well as DNA polymerase kinetic information that can be used for the direct detection of DNA modifications. We demonstrate that local sequence context has a strong impact on DNA polymerase kinetics in the neighborhood of the incorporation site during the DNA synthesis reaction, allowing for the possibility of estimating the expected kinetic rate of the enzyme at the incorporation site using kinetic rate information collected from existing SMRT sequencing data (historical data) covering the same local sequence contexts of interest. We develop an Empirical Bayesian hierarchical model for incorporating historical data. Our results show that the model could greatly increase DNA modification detection accuracy, and reduce requirement of control data coverage. For some DNA modifications that have a strong signal, a control sample is not even needed by using historical data as alternative to control. Thus, sequencing costs can be greatly reduced by using the model. We implemented the model in a R package named seqPatch, which is available at https://github.com/zhixingfeng/seqPatch.