Kinetics of allosteric enzymes

Peptidase Activation by a Leader Peptide-Bound RiPP Recognition Element

The RiPP precursor recognition element (RRE) is a conserved domain found in many prokaryotic ribosomally synthesized and post-translationally modified peptide (RiPP) biosynthetic gene clusters (BGCs). RREs bind with high specificity and affinity to a recognition sequence within the N-terminal leader region of RiPP precursor peptides. Lasso peptide biosynthesis involves an RRE-dependent leader peptidase, which is discretely encoded or fused to the RRE as a di-domain protein. Here we leveraged thousands of predicted BGCs to define the RRE:leader peptidase interaction through evolutionary covariance analysis. Each interacting domain contributes a three-stranded β-sheet to form a hydrophobic β-sandwich-like interface. The bioinformatics-guided predictions were experimentally confirmed using proteins from discrete and fused lasso peptide BGC architectures. Support for the domain-domain interface derived from chemical shift perturbation, paramagnetic relaxation enhancement experiments, and rapid variant activity screening using cell-free biosynthesis. Further validation of selected variants was performed with purified proteins. We developed a p-nitroanilide-based leader peptidase assay to illuminate the role of RRE domains. Our data show that RRE domains play a dual function. RRE domains deliver the precursor peptide to the leader peptidase, and the rate is saturable as expected for a substrate. RRE domains also partially compose the elusive S2 proteolytic pocket that binds the penultimate threonine of lasso leader peptides. Because the RRE domain is required to form the active site, leader peptidase activity is greatly diminished when the RRE domain is supplied at substoichiometric levels. Full proteolytic activation requires RRE engagement with the recognition sequence-containing portion of the leader peptide. Together, our observations define a new mechanism for protease activity regulation.

Counterintuitive Electrostatics upon Metal Ion Coordination to a Receptor with Two Homotopic Binding Sites

The consecutive binding of two potassium ions to a bis(18-crown-6) analogue of Tröger’s base (BCETB) in water was studied by isothermal titration calorimetry using four different salts, KCl, KI, KSCN, and K₂SO₄. A counterintuitive result was observed: the enthalpy change associated with the binding of the second ion is more negative than that of the first (ΔH_bind,2^° < ΔH_bind,1^°). This remarkable finding is supported by continuum electrostatic theory as well as by atomic scale replica exchange molecular dynamics simulations, where the latter robustly reproduces experimental trends for all simulated salts, KCl, KI, and KSCN, using multiple force fields. While an enthalpic K⁺–K⁺attraction in water poses a small, but fundamentally important, contribution to the overall interaction, the probability of the collapsed conformation (COL) of BCETB, where both crown ether moieties (CEs) of BCETB are bent in toward the cavity, was found to increase successively upon binding of the first and second potassium ions. The promotion of the COL conformation reveals favorable intrinsic interactions between the potassium coordinated CEs, which further contribute to the observation that ΔH_bind,2^° < ΔH_bind,1^°. While the observed trend is independent of the counterion, the origin of the significantly larger magnitude of the difference ΔH_bind,2^° – ΔH_bind,1^° observed experimentally for KSCN was studied in light of the weaker hydration of the thiocyanate anion, resulting in an enrichment of thiocyanate ions close to BCETB compared to the other studied counterions.

Enzymatic Noncovalent Synthesis

Enzymatic reactions and noncovalent (i.e., supramolecular) interactions are two fundamental nongenetic attributes of life. Enzymatic noncovalent synthesis (ENS) refers to a process where enzymatic reactions control intermolecular noncovalent interactions for spatial organization of higher-order molecular assemblies that exhibit emergent properties and functions. Like enzymatic covalent synthesis (ECS), in which an enzyme catalyzes the formation of covalent bonds to generate individual molecules, ENS is a unifying theme for understanding the functions, morphologies, and locations of molecular ensembles in cellular environments. This review intends to provide a summary of the works of ENS within the past decade and emphasize ENS for functions. After comparing ECS and ENS, we describe a few representative examples where nature uses ENS, as a rule of life, to create the ensembles of biomacromolecules for emergent properties/functions in a myriad of cellular processes. Then, we focus on ENS of man-made (synthetic) molecules in cell-free conditions, classified by the types of enzymes. After that, we introduce the exploration of ENS of man-made molecules in the context of cells by discussing intercellular, peri/intracellular, and subcellular ENS for cell morphogenesis, molecular imaging, cancer therapy, and other applications. Finally, we provide a perspective on the promises of ENS for developing molecular assemblies/processes for functions. This review aims to be an updated introduction for researchers who are interested in exploring noncovalent synthesis for developing molecular science and technologies to address societal needs.

Intracellular Ser/Thr/Tyr phosphoproteome of the oral commensal Streptococcus gordonii DL1

Background

To respond and adapt to environmental challenges, prokaryotes regulate cellular processes rapidly and reversibly through protein phosphorylation and dephosphorylation. This study investigates the intracellular proteome and Ser/Thr/Tyr phosphoproteome of the oral commensal Streptococcus gordonii. Intracellular proteins from planktonic cells of S. gordonii DL1 were extracted and subjected to 2D-gel electrophoresis. Proteins in general were visualized using Coomassie Brilliant Blue and T-Rex staining. Phosphorylated proteins were visualized with Pro-Q Diamond Phosphoprotein Gel Stain. Proteins were identified by LC-MS/MS and sequence analysis.

Results

In total, sixty-one intracellular proteins were identified in S. gordonii DL1, many of which occurred at multiple isoelectric points. Nineteen of these proteins were present as one or more Ser/Thr/Tyr phosphorylated form. The identified phosphoproteins turned out to be involved in a variety of cellular processes.

Conclusion

Nineteen phosphoproteins involved in various cellular functions were identified in S. gordonii. This is the first time the global intracellular Ser/Thr/Tyr phosphorylation profile has been analysed in an oral streptococcus. Comparison with phosphoproteomes of other species from previous studies showed many similarities. Proteins that are consistently found in a phosphorylated state across several species and growth conditions may represent a core phosphoproteome profile shared by many bacteria.

The structure and reactivity of the HoxEFU complex from the cyanobacterium Synechocystis sp. PCC 6803

Cyanobacterial Hox is a [NiFe] hydrogenase that consists of the hydrogen (H₂)-activating subunits HoxYH, which form a complex with the HoxEFU assembly to mediate reactions with soluble electron carriers like NAD(P)H and ferredoxin (Fdx), thereby coupling photosynthetic electron transfer to energy-transforming catalytic reactions. Researchers studying the HoxEFUYH complex have observed that HoxEFU can be isolated independently of HoxYH, leading to the hypothesis that HoxEFU is a distinct functional subcomplex rather than an artifact of Hox complex isolation. Moreover, outstanding questions about the reactivity of Hox with natural substrates and the site(s) of substrate interactions and coupling of H₂, NAD(P)H, and Fdx remain to be resolved. To address these questions, here we analyzed recombinantly produced HoxEFU by electron paramagnetic resonance spectroscopy and kinetic assays with natural substrates. The purified HoxEFU subcomplex catalyzed electron transfer reactions among NAD(P)H, flavodoxin, and several ferredoxins, thus functioning in vitro as a shuttle among different cyanobacterial pools of reducing equivalents. Both Fdx1-dependent reductions of NAD⁺ and NADP⁺ were cooperative. HoxEFU also catalyzed the flavodoxin-dependent reduction of NAD(P)⁺, Fdx2-dependent oxidation of NADH and Fdx4- and Fdx11-dependent reduction of NAD⁺. MS-based mapping identified an Fdx1-binding site at the junction of HoxE and HoxF, adjacent to iron-sulfur (FeS) clusters in both subunits. Overall, the reactivity of HoxEFU observed here suggests that it functions in managing peripheral electron flow from photosynthetic electron transfer, findings that reveal detailed insights into how ubiquitous cellular components may be used to allocate energy flow into specific bioenergetic products.

Colorimetric Detection of Salivary α-Amylase Using Maltose as a Noncompetitive Inhibitor for Polysaccharide Cleavage

This paper describes an approach for colorimetric detection of salivary α-amylase, one of the potential biomarkers of autonomic nervous system (ANS) activity, for enabling assessment of fatigue. The ability of α-amylase to cleave α-bonds of polysaccharides is utilized for developing a colorimetric assay. In the proposed approach, 2-chloro-4-nitrophenyl-α-d-maltotrioside as substrate releases a colored byproduct upon cleavage by salivary α-amylase. Introduction of maltose as a noncompetitive inhibitor yields desirable linear responses in the physiologically relevant concentration range (20-500 μg/mL) with a limit of detection (LOD) of 8 μg/mL (in aqueous solution). The concentrations of substrate and noncompetitive inhibitor are subsequently optimized for colorimetric detection of salivary α-amylase. A facile paper-based "strip" assay is proposed for analysis of human saliva samples with marginal interference from saliva components. The proposed assay is rapid, specific, and easy-to-implement for colorimetric detection of salivary α-amylase between 20 and 500 μg/mL. Complementary RGB (red, green, blue components) analysis offers quantitative detection with a LOD of 11 μg/mL. The two assay formats are benchmarked against the Phadebas test, a state of the art method for spectrophotometric detection of α-amylase. The reported paper-based methodology possesses a high potential for estimation of altered ANS responses toward stressors that possibly could find applications in assessment of fatigue and for monitoring onset of fatigue.

A broader active site in Pyrococcus horikoshii CoA disulfide reductase accommodates larger substrates and reveals evidence of subunit asymmetry

Within the family of pyridine nucleotide disulfide oxidoreductase (PNDOR), enzymes are a group of single‐cysteine containing FAD‐dependent reductases that utilize a tightly bound coenzyme A to assist in the NAD(P)H‐dependent reduction of di‐, per‐, and polysulfide substrates in bacteria and archaea. For many of these homodimeric enzymes, it has proved difficult to determine the substrate specificity and metabolic function based on sequence and genome analysis alone. Coenzyme A‐disulfide reductase (CoADR) isolated from Pyrococcus horikoshii (phCoADR) reduces Co‐A per‐ and polysulfides, but, unlike other highly homologous members of this group, is a poor CoA disulfide reductase. The phCoADR structure has a narrower access channel for CoA substrates, which suggested that this restriction might be responsible for the enzyme's poor activity toward the bulky CoA disulfide substrate. To test this hypothesis, the substrate channel was widened by making four mutations along the channel wall (Y65A, Y66A, P67G, and H367G). The structure of the quadruple mutant shows a widened substrate channel, which is supported by a fourfold increase in k_cat for the NAD(P)H‐dependent reduction of CoA disulfide and enhanced activity toward the substrate at lower temperatures. Anaerobic titrations of the enzyme with NADH revealed a half‐site reactivity not observed with the wild‐type enzyme in which one subunit of the enzyme could be fully reduced to an EH₄ state, while the other remained in an EH₂ or EH₂·NADH state. These results suggest that for these closely related enzymes, substrate channel morphology is an important determinant of substrate specificity, and homology modeling will be the preferred technique for predicting function among PNDORs.

Mechanism of partial agonism in AMPA-type glutamate receptors

Neurotransmitters trigger synaptic currents by activating ligand-gated ion channel receptors. Whereas most neurotransmitters are efficacious agonists, molecules that activate receptors more weakly-partial agonists-also exist. Whether these partial agonists have weak activity because they stabilize less active forms, sustain active states for a lesser fraction of the time or both, remains an open question. Here we describe the crystal structure of an α-amino-3-hydroxy-5-methyl-4-isoxazolepropionate receptor (AMPAR) ligand binding domain (LBD) tetramer in complex with the partial agonist 5-fluorowillardiine (FW). We validate this structure, and others of different geometry, using engineered intersubunit bridges. We establish an inverse relation between the efficacy of an agonist and its promiscuity to drive the LBD layer into different conformations. These results suggest that partial agonists of the AMPAR are weak activators of the receptor because they stabilize multiple non-conducting conformations, indicating that agonism is a function of both the space and time domains.

Anaerobic Heme Degradation: ChuY Is an Anaerobilin Reductase That Exhibits Kinetic Cooperativity

Heme catabolism is an important biochemical process that many bacterial pathogens utilize to acquire iron. However, tetrapyrrole catabolites can be reactive and often require further processing for transport out of the cell or conversion to another useful cofactor. In previous work, we presented in vitro evidence of an anaerobic heme degradation pathway in Escherichia coli O157:H7. Consistent with reactions that have been reported for other radical S-adenosyl-l-methionine methyltransferases, ChuW transfers a methyl group to heme by a radical-mediated mechanism and catalyzes the β-scission of the porphyrin macrocycle. This facilitates iron release and the production of a new linear tetrapyrrole termed "anaerobilin". In this work, we describe the structure and function of ChuY, an enzyme expressed downstream from chuW within the same heme utilization operon. ChuY is structurally similar to biliverdin reductase and forms a dimeric complex in solution that reduces anaerobilin to the product we have termed anaerorubin. Steady state analysis of ChuY exhibits kinetic cooperativity that is best explained by a random addition mechanism with a kinetically preferred path for initial reduced nicotinamide adenine dinucleotide phosphate binding.

A Simple Model of Multivalent Adhesion and Its Application to Influenza Infection

Adhesion between biological surfaces, which is typically the result of molecular binding between receptors on one surface and ligands on another, plays a fundamental role in biology and is key to the infection mechanisms of certain viruses, including influenza. The physiological outcome of adhesion depends on both the number of bound cells (or viruses, or other biological particles) and the properties of the adhesion interface that is formed, including the equilibrium number of receptor-ligand connections. Here, we introduce a quantitative model for biological adhesion by adapting thermodynamic models developed for the related problem of multivalent molecular binding. In our model, adhesion affinity is approximated by a simple, analytical expression involving the numbers of ligands and receptors at the interface. Our model contains only two fitting parameters and is simple to interpret. When applied to the adhesion between the hemagglutinin ligands on influenza viruses and the sialic acid receptors on biosensors or on host cells, our model generates adhesion affinities consistent with experimental measurements performed over a range of numbers of receptors, and provides a semiquantitative estimate of the affinity range of the hemagglutinin-sialic acid interaction necessary for the influenza virus to successfully infect host cells. The model also provides a quantitative explanation for the experimental finding that a mutant avian virus gained transmissibility in mammals despite the mutations conferring only a less than twofold increase in the affinity of its hemagglutinin for mammalian receptors: the model predicts an order-of-magnitude improvement in adhesion to mammalian cells. We also extend our model to describe the competitive inhibition of adhesion: the model predicts that hemagglutinin inhibitors of relatively modest affinity can dramatically reduce influenza virus adhesion to host cells, suggesting that such inhibitors, if discovered, may be viable therapeutic agents against influenza.

Investigating the Reaction Mechanism of F420-Dependent Glucose-6-phosphate Dehydrogenase from Mycobacterium tuberculosis: Kinetic Analysis of the Wild-Type and Mutant Enzymes

F₄₂₀-dependent glucose-6-phosphate dehydrogenase (FGD) catalyzes the conversion of glucose-6-phosphate (G6P) to 6-phosphogluconolactone, using F₄₂₀ cofactor as the hydride transfer acceptor, within mycobacteria. A previous crystal structure of wild-type FGD led to a proposed mechanism suggesting that the active site residues His40, Trp44, and Glu109 could be involved in catalysis. We have characterized the wild-type FGD and five FGD variants (H40A, W44F, W44Y, W44A, and E109Q) by fluorescence binding assays and steady-state and pre-steady-state kinetic experiments. Compared to wild-type FGD, all the variants had lower binding affinities for F₄₂₀, thus suggesting that Trp44, His40, and Glu109 aid in F₄₂₀ binding. While all the variants had decreased catalytic efficiencies, FGD H40A and W44A were the least efficient, having lost ∼1000- and ∼2000-fold activity, respectively. This confirms a crucial catalytic role for His40 in the FGD reaction and suggests that aromaticity at residue 44 aids catalysis. To investigate the proposed roles of Glu109 and His40 in acid-base catalysis, the pH dependence of kinetic parameters has been determined for the E109Q and H40A mutants and compared to those of the wild-type enzyme. The log k_cat-pH profile of wild-type FGD and E109Q revealed two ionizable residues in the enzyme-substrate complex, while H40A displayed only one ionization event. The FGD E109Q variant displayed pH-dependent kinetic cooperativity with respect to the F₄₂₀ cofactor. The multiple-turnover pre-steady-state kinetics were biphasic for wild-type FGD, W44F, W44Y, and E109Q, while the H40A and W44A variants displayed only a single phase because of their reduced catalytic efficiency.

Hierarchical nanostructure and synergy of multimolecular signalling complexes

Signalling complexes are dynamic, multimolecular structures and sites for intracellular signal transduction. Although they play a crucial role in cellular activation, current research techniques fail to resolve their structure in intact cells. Here we present a multicolour, photoactivated localization microscopy approach for imaging multiple types of single molecules in fixed and live cells and statistical tools to determine the nanoscale organization, topology and synergy of molecular interactions in signalling complexes downstream of the T-cell antigen receptor. We observe that signalling complexes nucleated at the key adapter LAT show a hierarchical topology. The critical enzymes PLCγ1 and VAV1 localize to the centre of LAT-based complexes, and the adapter SLP-76 and actin molecules localize to the periphery. Conditional second-order statistics reveal a hierarchical network of synergic interactions between these molecules. Our results extend our understanding of the nanostructure of signalling complexes and are relevant to studying a wide range of multimolecular complexes.

Protein Allostery and Conformational Dynamics

The functions of many proteins are regulated through allostery, whereby effector binding at a distal site changes the functional activity (e.g., substrate binding affinity or catalytic efficiency) at the active site. Most allosteric studies have focused on thermodynamic properties, in particular, substrate binding affinity. Changes in substrate binding affinity by allosteric effectors have generally been thought to be mediated by conformational transitions of the proteins or, alternatively, by changes in the broadness of the free energy basin of the protein conformational state without shifting the basin minimum position. When effector binding changes the free energy landscape of a protein in conformational space, the change affects not only thermodynamic properties but also dynamic properties, including the amplitudes of motions on different time scales and rates of conformational transitions. Here we assess the roles of conformational dynamics in allosteric regulation. Two cases are highlighted where NMR spectroscopy and molecular dynamics simulation have been used as complementary approaches to identify residues possibly involved in allosteric communication. Perspectives on contentious issues, for example, the relationship between picosecond-nanosecond local and microsecond-millisecond conformational exchange dynamics, are presented.

Biochemical studies on Francisella tularensis RelA in (p)ppGpp biosynthesis

Francisella tularensis RelA shows significant sequence differences from other members of the RelA family of enzymes. In the present study, we describe the functional similarities and differences between F. tularensis RelA and the model RelA from Escherichia coli.

The bacterial stringent response is induced by nutrient deprivation and is mediated by enzymes of the RSH (RelA/SpoT homologue; RelA, (p)ppGpp synthetase I; SpoT, (p)ppGpp synthetase II) superfamily that control concentrations of the ‘alarmones’ (p)ppGpp (guanosine penta- or tetra-phosphate). This regulatory pathway is present in the vast majority of pathogens and has been proposed as a potential anti-bacterial target. Current understanding of RelA-mediated responses is based on biochemical studies using Escherichia coli as a model. In comparison, the Francisella tularensis RelA sequence contains a truncated regulatory C-terminal region and an unusual synthetase motif (EXSD). Biochemical analysis of F. tularensis RelA showed the similarities and differences of this enzyme compared with the model RelA from Escherichia coli. Purification of the enzyme yielded a stable dimer capable of reaching concentrations of 10 mg/ml. In contrast with other enzymes from the RelA/SpoT homologue superfamily, activity assays with F. tularensis RelA demonstrate a high degree of specificity for GTP as a pyrophosphate acceptor, with no measurable turnover for GDP. Steady state kinetic analysis of F. tularensis RelA gave saturation activity curves that best fitted a sigmoidal function. This kinetic profile can result from allosteric regulation and further measurements with potential allosteric regulators demonstrated activation by ppGpp (5′,3′-dibisphosphate guanosine) with an EC₅₀ of 60±1.9 μM. Activation of F. tularensis RelA by stalled ribosomal complexes formed with ribosomes purified from E. coli MRE600 was observed, but interestingly, significantly weaker activation with ribosomes isolated from Francisella philomiragia.

Evidence of Kinetic Cooperativity in Dimeric Ketopantoate Reductase from Staphylococcus aureus

Ketopantoate reductase (KPR) catalyzes the NADPH-dependent production of pantoate, an essential precursor in the biosynthesis of coenzyme A. Previous structural studies have been limited to Escherichia coli KPR, a monomeric enzyme that follows a sequential ordered mechanism. Here we report the crystal structure of the Staphylococcus aureus enzyme at 1.8 Å resolution, the first description of a dimeric KPR. Using sedimentation velocity analysis, we show that the S. aureus KPR dimer is stable in solution. In fact, our structural analysis shows that the dimeric assembly we identify is present in the majority of KPR crystal structures. Steady state analysis of S. aureus KPR reveals strong positive cooperativity with respect to NADPH (Hill coefficient of 2.5). In contrast, high concentrations of the substrate ketopantoate (KP) inhibit the activity of the enzyme. These observations are consistent with a random addition mechanism in which the initial binding of NADPH is the kinetically preferred path. In fact, Förster resonance energy transfer studies of the equilibrium binding of NADPH show only a small degree of cooperativity between subunits (Hill coefficient of 1.3). Thus, the apparently strong cooperativity observed in substrate saturation curves is due to a kinetic process that favors NADPH binding first. This interpretation is consistent with our analysis of the A181L substitution, which increases the Km of ketopantoate 844-fold, without affecting kcat. The crystal structure of KPRA181L shows that the substitution displaces Ser239, which is known to be important for the binding affinity of KP. The decrease in KP affinity would enhance the already kinetically preferred NADPH binding path, making the random mechanism appear to be sequentially ordered and reducing the kinetic cooperativity. Consistent with this interpretation, the NADPH saturation curve for KPRA181L is hyperbolic.

On the estimation of cooperativity in ion channel kinetics: activation free energy and kinetic mechanism of Shaker K+ channel

In this paper, we have explored generic criteria of cooperative behavior in ion channel kinetics treating it on the same footing with multistate receptor-ligand binding in a compact theoretical framework. We have shown that the characterization of cooperativity of ion channels in terms of the Hill coefficient violates the standard Hill criteria defined for allosteric cooperativity of ligand binding. To resolve the issue, an alternative measure of cooperativity is proposed here in terms of the cooperativity index that sets a unified criteria for both the systems. More importantly, for ion channel this index can be very useful to describe the cooperative kinetics as it can be readily determined from the experimentally measured ionic current combined with theoretical modelling. We have analyzed the correlation between the voltage value and slope of the voltage-activation curve at the half-activation point and consequently determined the standard free energy of activation of the ion channel using two well-established mechanisms of cooperativity, namely, Koshland-Nemethy-Filmer (KNF) and Monod-Wyman-Changeux (MWC) models. Comparison of the theoretical results for both the models with appropriate experimental data of mutational perturbation of Shaker K(+) channel supports the experimental fact that the KNF model is more suitable to describe the cooperative behavior of this class of ion channels, whereas the performance of the MWC model is unsatisfactory. We have also estimated the mechanistic performance through standard free energy of channel activation for both the models and proposed a possible functional disadvantage in the MWC scheme.

Organism-adapted specificity of the allosteric regulation of pyruvate kinase in lactic acid bacteria

Pyruvate kinase (PYK) is a critical allosterically regulated enzyme that links glycolysis, the primary energy metabolism, to cellular metabolism. Lactic acid bacteria rely almost exclusively on glycolysis for their energy production under anaerobic conditions, which reinforces the key role of PYK in their metabolism. These organisms are closely related, but have adapted to a huge variety of native environments. They include food-fermenting organisms, important symbionts in the human gut, and antibiotic-resistant pathogens. In contrast to the rather conserved inhibition of PYK by inorganic phosphate, the activation of PYK shows high variability in the type of activating compound between different lactic acid bacteria. System-wide comparative studies of the metabolism of lactic acid bacteria are required to understand the reasons for the diversity of these closely related microorganisms. These require knowledge of the identities of the enzyme modifiers. Here, we predict potential allosteric activators of PYKs from three lactic acid bacteria which are adapted to different native environments. We used protein structure-based molecular modeling and enzyme kinetic modeling to predict and validate potential activators of PYK. Specifically, we compared the electrostatic potential and the binding of phosphate moieties at the allosteric binding sites, and predicted potential allosteric activators by docking. We then made a kinetic model of Lactococcus lactis PYK to relate the activator predictions to the intracellular sugar-phosphate conditions in lactic acid bacteria. This strategy enabled us to predict fructose 1,6-bisphosphate as the sole activator of the Enterococcus faecalis PYK, and to predict that the PYKs from Streptococcus pyogenes and Lactobacillus plantarum show weaker specificity for their allosteric activators, while still having fructose 1,6-bisphosphate play the main activator role in vivo. These differences in the specificity of allosteric activation may reflect adaptation to different environments with different concentrations of activating compounds. The combined computational approach employed can readily be applied to other enzymes.

Computer simulation of assembly and co-operativity of hexameric AAA ATPases

AAA ATPases form a functionally diverse superfamily of proteins. Most members form homo-hexameric ring complexes, are catalytically active only in the fully assembled state, and show co-operativity among the six subunits. The mutual dependence among the subunits is clearly evidenced by the fact that incorporation of mutated, inactive subunits can decrease the activity of the remaining wild type subunits. For the first time, we develop here models to describe this form of allostery, evaluate them in a simulation study, and test them on experimental data. We show that it is important to consider the assembly reactions in the kinetic model, and to define a formal inhibition scheme. We simulate three inhibition scenarios explicitly, and demonstrate that they result in differing outcomes. Finally, we deduce fitting formulas, and test them on real and simulated data. A non-competitive inhibition formula fitted experimental and simulated data best. To our knowledge, our study is the first one that derives and tests formal allosteric schemes to explain the inhibitory effects of mutant subunits on oligomeric enzymes.

An enzymatic pathway for the biosynthesis of the formylhydroxyornithine required for rhodochelin iron coordination

Rhodochelin, a mixed catecholate-hydroxamate type siderophore isolated from Rhodococcus jostii RHA1, holds two L-δ-N-formyl-δ-N-hydroxyornithine (L-fhOrn) moieties essential for proper iron coordination. Previously, bioinformatic and genetic analysis proposed rmo and rft as the genes required for the tailoring of the L-ornithine (L-Orn) precursor [Bosello, M. (2011) J. Am. Chem. Soc.133, 4587-4595]. In order to investigate if both Rmo and Rft constitute a pathway for L-fhOrn biosynthesis, the enzymes were heterologously produced and assayed in vitro. In the presence of molecular oxygen, NADPH and FAD, Rmo monooxygenase was able to convert L-Orn into L-δ-N-hydroxyornithine (L-hOrn). As confirmed in a coupled reaction assay, this hydroxylated intermediate serves as a substrate for the subsequent N(10)-formyl-tetrahydrofolate-dependent (N(10)-fH(4)F) Rtf-catalyzed formylation reaction, establishing a route for the L-fhOrn biosynthesis, prior to its incorporation by the NRPS assembly line. It is of particular interest that a major improvement to this study has been reached with the use of an alternative approach to the chemoenzymatic FolD-dependent N(10)-fH(4)F conversion, also rescuing the previously inactive CchA, the Rft-homologue in coelichelin assembly line [Buchenau, B. (2004) Arch. Microbiol.182, 313-325; Pohlmann, V. (2008) Org. Biomol. Chem.6, 1843-1848].

Analytical expressions for the homotropic binding of ligand to protein dimers and trimers

Cooperative binding of a ligand to multiple subsites on a protein is a common theme among enzymes and receptors. The analysis of cooperative binding data (either positive or negative) often relies on the assumption that free ligand concentration, L, can be approximated by the total ligand concentration, L(T). When this approximation does not hold, such analyses result in inaccurate estimates of dissociation constants. Presented here are exact analytical expressions for equilibrium concentrations of all enzyme and ligand species (in terms of K(d) values and total concentrations of protein and ligand) for homotropic dimeric and trimeric protein-ligand systems. These equations circumvent the need to approximate L and are provided in Excel worksheets suitable for simulation and least-squares fitting. The equations and worksheets are expanded to treat cases where binding signals vary with distinct site occupancy.

Curcuminoid biosynthesis by two type III polyketide synthases in the herb Curcuma longa

Curcuminoids found in the rhizome of turmeric, Curcuma longa, possess various biological activities. Despite much attention regarding the biosynthesis of curcuminoids because of their pharmaceutically important properties and biosynthetically intriguing structures, no enzyme systems have been elucidated. Here we propose a pathway for curcuminoid biosynthesis in the herb C. longa, which includes two novel type III polyketide synthases. One of the type III polyketide synthases, named diketide-CoA synthase (DCS), catalyzed the formation of feruloyldiketide-CoA by condensing feruloyl-CoA and malonyl-CoA. The other, named curcumin synthase (CURS), catalyzed the in vitro formation of curcuminoids from cinnamoyldiketide-N-acetylcysteamine (a mimic of the CoA ester) and feruloyl-CoA. Co-incubation of DCS and CURS in the presence of feruloyl-CoA and malonyl-CoA yielded curcumin at high efficiency, although CURS itself possessed low activity for the synthesis of curcumin from feruloyl-CoA and malonyl-CoA. These findings thus revealed the curcumin biosynthetic route in turmeric, in which DCS synthesizes feruloyldiketide-CoA, and CURS then converts the diketide-CoA esters into a curcuminoid scaffold.

Substrate specificity and diastereoselectivity of strictosidine glucosidase, a key enzyme in monoterpene indole alkaloid biosynthesis

Strictosidine glucosidase (SGD) from Catharanthus roseus catalyzes the deglycosylation of strictosidine, an intermediate from which thousands of monoterpene indole alkaloids are derived. The steady-state kinetics of SGD with a variety of strictosidine analogs revealed the substrate preferences of this enzyme at two key positions of the strictosidine substrate. Additionally, SGD from C. roseus turns over both strictosidine and its stereoisomer vincoside, indicating that although this enzyme prefers the naturally occurring diastereomer, the enzyme is not completely diastereoselective. The implications of the substrate specificity of SGD in metabolic engineering efforts of C. roseus are highlighted.

Synthesis and Inclusion Properties of 1,3-alternate-bis[(ethoxycarbonyl) Methoxy]thiacalix[4]arene-mono(crown-4) Ether

A novel ditopic receptor possessing two complexation sites such as crown ether and (ethoxycarbonyl) methoxy groups bearing 1,3-alternate conformation based on thiacalix[4]arene was prepared. The binding behaviours with alkali metals have been examined by 1H NMR titration experiment. The exclusive formation of the heterogeneous dinuclear complexes of 1,3-alternate-3 with Li+ and K+ was observed.

A mutation in the S-switch region of the Runt domain alters the dynamics of an allosteric network responsible for CBFbeta regulation

The Runt domain is the DNA binding domain of the core binding factor (CBF) Runx subunits. The CBFs are transcription factors that play critical roles in hematopoiesis, bone, and neuron development in mammals. A common non-DNA binding CBFbeta subunit heterodimerizes with the Runt domain of the Runx proteins and allosterically regulates its affinity for DNA. Previous NMR dynamics studies suggested a model whereby CBFbeta allosterically regulates DNA binding by quenching conformational exchange in the Runt domain, particularly in the S-switch region and the betaE'-F loop. We sought to test this model, and to this end introduced all possible single amino acid substitutions into the S-switch region and the betaE'-F loop, and screened for mutations that enhanced DNA-binding. We demonstrate that one Runt domain mutant, R164N, binds both DNA and CBFbeta with higher affinity, but it is less sensitive to allosteric regulation by CBFbeta. Analysis of NMR relaxation data shows that the chemical exchange exhibited by the wild-type Runt domain is largely quenched by the R164N substitution. These data support a model in which the dynamic behavior of a network of residues connecting the CBFbeta and DNA binding sites on the Runt domain plays a critical role in the mechanism of allosteric regulation. This study provides an important functional link between dynamic behavior and protein allosteric function, consistent with results on other allosterically regulated proteins.

Slow transitions between two conformational states of band 3 (AE1) modulate divalent anion transport and DBDS binding to a second site on band 3 which is activated by lowering the pH (pK approximately 5.0)

Evidence is emerging which indicates that the anion transport activity of band 3 may be regulated. I review the molecular basis for regulation of the anion transport function of band 3 in terms of evidence showing that divalent anion transport involves a slow "hysteretic" transition between two functional states, mediated by interactions between subunits within band 3 oligomers. In addition, I briefly describe recent work from my laboratory where we have discovered a novel manifestation of slow conformational changes in band 3. This involves 4,4'-dibenzamido-2,2'-stilbenedisulfonate (DBDS) binding to a second, proton-activated site distinct from the primary stilbenedisulfonate site. This site is exposed on the outer aspect of band 3 when the pH is lowered (pK approximately 5.0). This is the same pK as the protonation site on band 3 involved in divalent anion-proton co-transport (APCT) [J. Gen. Physiol. 79 (1982) 87]. Significantly, we have found that DBDS binding to this proton-activated site has unusually slow kinetics, and increasing DBDS concentration causes a decrease in the apparent pseudo-first-order rate constant. These results suggest that a slow conformational pre-equilibrium is the rate limiting step in DBDS binding to the proton-activated site on band 3 observed at low pH. Our results support an allosteric two-state model for regulation of divalent anion transport by band 3 oligomers involving a slow conformational transition and interactions between subunits [Biochemistry 31 (1992) 7301].

CBFbeta allosterically regulates the Runx1 Runt domain via a dynamic conformational equilibrium

Core binding factors (CBFs) are heterodimeric transcription factors consisting of a DNA-binding CBFalpha subunit and non-DNA-binding CBFbeta subunit. The CBFbeta subunit increases the affinity of the DNA-binding Runt domain of CBFalpha for DNA while making no direct contacts to the DNA. We present evidence for conformational exchange in the S-switch region in a Runt domain-DNA complex that is quenched upon CBFbeta binding. Analysis of (15)N backbone relaxation parameters shows that binding of CBFbeta reduces the backbone dynamics in the microsecond-to-millisecond time frame for several regions of the Runt domain that make energetically important contacts with the DNA. The DNA also undergoes conformational exchange in the Runt domain-DNA complex that is quenched in the presence of CBFbeta. Our results indicate that allosteric regulation by the CBFbeta subunit is mediated by a shift in an existing dynamic conformational equilibrium of both the Runt domain and DNA.

The E. coli replication factor DnaC protein exists in two conformations with different nucleotide binding capabilities. I. Determination of the binding mechanism using ATP and ADP fluorescent analogues

The kinetic mechanism of binding of ATP and ADP fluorescent analogues to the E. coli replicative factor DnaC protein has been studied using the fluorescence stopped-flow technique. The experiments have been performed under pseudo-first-order conditions with respect to the nucleotide cofactor or the DnaC concentration. Three relaxation processes are observed at a large excess of the nucleotide, while only two relaxation processes are detected in the excess of the protein. Such behavior of the kinetic system is a diagnostic indication of the presence of the protein conformational equilibrium prior to the ligand binding. The obtained data indicate that the minimum mechanism that describes the observed kinetics includes the conformational transition of the DnaC protein, prior to nucleotide binding, followed by the two-step, sequential association of the cofactor to only one of the protein conformations, as defined by In the examined solution conditions, the conformation of the DnaC protein is shifted toward the state (DnaC)(2) that binds the nucleotide. The lack of any cofactor binding to the (DnaC)(1) state points to the existence of a stringent locking mechanism of the nucleotide binding-site in the protein. Binding of ATP and ADP analogues obeys the same mechanism, with similar rate constants, indicating that ATP and ADP analogues bind to the same protein conformation. The (C)(1) intermediate dominates the distribution of the DnaC protein population in the presence of cofactors. The formation of (C)(1) is accompanied by a low nucleotide fluorescence increase, indicating a hydrophilic environment around the ribose of bound cofactors. Transition to (C)(2) places the ribose region in a highly hydrophobic environment with relative molar fluorescence intensity approximately 8-fold higher than that of the free cofactor. The significance of these results for the functioning of the DnaC protein is discussed.

Multiple conformational changes in enzyme catalysis

Understanding the molecular mechanisms of enzyme catalysis and allosteric regulation has been a primary goal of biochemistry for many years. The dynamics of these processes, approached through a variety of kinetic methods, are discussed. The results obtained for many different enzymes suggest that multiple intermediates and conformations are general characteristics of the catalytic process and allosteric regulation. Ribonuclease, dihydrofolate reductase, chymotrypsin, aspartate aminotransferase, and aspartate transcarbamoylase are considered as specific examples. Typical and maximum rates of conformational changes and catalysis are also discussed, based on results obtained from model systems. The nature and rates of interconversion of the intermediates, along with structural information, can be used as the bases for understanding the incredible catalytic efficiency of enzymes. Potential roles of conformational changes in the catalytic process are discussed in terms of static and environmental effects, and in terms of dynamic coupling within the enzyme-substrate complex.

Physical and chemical properties of DMP 504, a polyalkylammonium-based bile acid sequestrant

The purpose of this study was to characterize the physicochemical properties of DMP 504 and lay the foundation for formulation development. Thermal properties were characterized by DSC and TGA and moisture sorption and desorption by TGA. The association rate and equilibrium binding capacity of the polymer for a prototype bile acid was evaluated using cholic acid, and solid state stability was examined at 25 degrees C, 40 degrees C (with and without 5% added water), 60 degrees C, and 600 foot candles/25 degrees C. The solid state excipient compatibility of binary mixtures of DMP 504 and several commonly used pharmaceutical excipients was also evaluated. Thermal analysis of the polymer showed a glass transition temperature at approximately 95 degrees C and no melting point, indicating a highly amorphous macromolecular structure with thermal stability up to 250 degrees C. Moisture sorption and desorption isotherms at controlled humidity ranging from 11% to 97% RH did not display hysteresis. Cholic acid associated with DMP 504 extremely rapidly so that binding was essentially complete within 5 min. Scatchard analysis of the equilibrium binding of cholic acid to DMP 504 was unconventional, and indicated that the system was exhibiting positive cooperative behavior. Modeling the binding curve for a system exhibiting cooperative behavior indicated a maximum binding capacity of DMP 504 for cholic acid in phosphate buffered saline (pH 7.0) of 4.9 mumol/mg, and a cooperativity value, P, of 2.2 implying that binding of one molecule promotes the binding of additional molecules. DMP 504, a hygroscopic, amorphous cross-linked polymer with a tendency to gain or lose moisture with ease, is stable in the solid state, either drug substance alone or in presence of excipients, at normal storage temperatures and light, and under controlled conditions of humidity.

Effect of a benzodiazepine (chlordiazepoxide) on a GABAA receptor from rat brain. Requirement of only one bound GABA molecule for channel opening

Chlordiazepoxide (CDPX) enhanced the rate of chloride exchange mediated by the major GABAA receptor found on sealed native membrane vesicles from rat cerebral cortex. The initial rate constant for chloride exchange for this receptor, (JA), a measure of open channel, was determined from the progress of GABA-mediated influx of 36Cl-. The dependence of JA on GABA concentration was hyperbolic in the presence of CDPX (150 microM, sufficient to give maximum enhancement of chloride exchange rate) but sigmoid in its absence. Enhancement of channel opening (10-fold at 0.3 microM GABA) decreased with increasing GABA concentration. The maximal response, above 1,000 microM GABA, was unaltered. The half-response concentration was reduced from 80 microM to 50 microM. CDPX alone caused no measurable 36Cl- exchange. In the presence of CDPX, channel opening occurred with only one bound GABA molecule, whereas in its absence, channel opening with two bound GABA molecules was much more favorable. This could not be direct allosteric modulation of the channel opening conformational change by binding of CDPX at effector sites, but could be explained by an additional change of the receptor on binding CDPX to give a closed state which gave channel opening mediated by a single GABA binding site. Another possibility is that CDPX could act at one of the channel opening binding sites without a postulated, second closed conformational state.

Molecular characterization of the Enterococcus faecalis cytolysin activator

The gene encoding component A (cylA), the activator protein of the Enterococcus faecalis cytolysin, has been localized on pAD1, and the nucleotide sequence was determined. cylA consists of a 1,236-bp open reading frame encoding a 412-amino-acid polypeptide. A search of the National Biomedical Research Foundation data base revealed significant homology between the inferred amino acid sequence of component A and subtilisin BPN'. Component A activation of the cytolysin precursor (component L) was observed to be inhibited by the serine protease inhibitor diisopropylfluorophosphate. Mature component A exhibits a molecular weight of approximately 30,000 and an isoelectric point of 4.5. Differences between the size of the primary translation product (45,625 daltons) and the mature enzyme suggest that, as for subtilisin, component A is secreted as a proenzyme. These results provide the basis for a model of component A activation of component L and a role for component A in protecting the cytolysin-producing cell from lysis.

Acetylcholine receptor (from Electrophorus electricus): a comparison of single-channel current recordings and chemical kinetic measurements

We report a direct comparison between two types of measurements of the dynamic properties of the acetylcholine receptor: single-channel currents recorded using the patch-clamp technique and chemical kinetic measurements. Electrophorus electricus electroplax cells, and membrane vesicles prepared from these cells, were used. Such a comparison, and single-channel currents recorded from these cells, have not previously been reported. We first give the theoretical basis for the comparison and define the conditions under which the comparisons are elegantly simple. We relate (i) measurements of currents through receptor channels in the cell membranes to measurements of the rates of ion translocation through the receptor channels in vesicles and (ii) measurements of the lifetimes of receptor states (for instance, the lifetime of the active state of the receptor--i.e., the state in which it can form open channels) to rate coefficients obtained in chemical kinetic measurements (for instance, those for the interconversions between different states of the receptor). In eel Ringer's solution we have found the single-channel conductance (gamma) of the receptor in E. electricus electroplax cells to be 53 pS. From this value, a specific reaction rate for ion translocation, J, of 5 X 10(7) M-1 X sec-1 was calculated. When membrane vesicles prepared from the electroplax cells and the same solution compositions were used, chemical kinetic measurements gave a J value of 3 X 10(7) M-1 X sec-1. The agreement between the two measurements is important because (i) they reflect different experimental conditions, which require different assumptions in interpreting the results, and (ii) it indicates that the two techniques can be used to obtain complementary information: the methods have different time resolutions and can be used in different ranges of acetylcholine concentrations.