Impact of phosphorylation on structure and thermodynamics of the interaction between the N-terminal domain of enzyme I and the histidine phosphocarrier protein of the bacterial phosphotransferase system

Elucidation of the Mechanisms of Inter-domain Coupling in the Monomeric State of Enzyme I by High-pressure NMR

The catalytic cycle of Enzyme I (EI), a phosphotransferase enzyme responsible for converting phosphoenolpyruvate (PEP) into pyruvate, is characterized by a series of local and global conformational rearrangements. This multistep process includes a monomer-to-dimer transition, followed by an open-to-closed rearrangement of the dimeric complex upon PEP binding. In the present study, we investigate the thermodynamics of EI dimerization using a range of high-pressure solution NMR techniques complemented by SAXS experiments. 1H-15N TROSY and 1H-13C methyl TROSY NMR spectra combined with 15N relaxation measurements revealed that a native-like engineered variant of full-length EI fully dissociates into stable monomeric state above 1.5 kbar. Conformational ensembles of EI monomeric state were generated via a recently developed protocol combining coarse-grained molecular simulations with experimental backbone residual dipolar coupling measurements. Analysis of the structural ensembles provided detailed insights into the molecular mechanisms driving formation of the catalytically competent dimeric state, and reveals that each step of EI catalytical cycle is associated with a significant reduction in either inter- or intra-domain conformational entropy. Altogether, this study completes a large body work conducted by our group on EI and establishes a comprehensive structural and dynamical description of the catalytic cycle of this prototypical multidomain, oligomeric enzyme.

Kinetics of diffusion-influenced multisite phosphorylation with enzyme reactivation

The simplest way to account for the influence of diffusion on the kinetics of multisite phosphorylation is to modify the rate constants in the conventional rate equations of chemical kinetics. We have previously shown that this is not enough and new transitions between the reactants must also be introduced. Here we extend our results by considering enzymes that are inactive after modifying the substrate and need time to become active again. This generalization leads to a surprising result. The introduction of enzyme reactivation results in a diffusion-modified kinetic scheme with a new transition that has a negative rate constant. The reason for this is that mapping non-Markovian rate equations onto Markovian ones with time-independent rate constants is not a good approximation at short times. We then developed a non-Markovian theory that involves memory kernels instead of rate constants. This theory is now valid at short times, but is more challenging to use. As an example, the diffusion-modified kinetic scheme with new connections was used to calculate kinetics of double phosphorylation and steady-state response in a phosphorylation-dephosphorylation cycle. We have reproduced the loss of bistability in the phosphorylation-dephosphorylation cycle when the enzyme reactivation time decreases, which was obtained by particle-based computer simulations.

Protein Conformational Dynamics Underlie Selective Recognition of Thermophilic over Mesophilic Enzyme I by a Substrate Analogue

Substrate selectivity is an important preventive measure to decrease the possibility of cross interactions between enzymes and metabolites that share structural similarities. In addition, understanding the mechanisms that determine selectivity towards a particular substrate increases the knowledge base for designing specific inhibitors for target enzymes. Here, we combine NMR, molecular dynamics (MD) simulations, and protein engineering to investigate how two substrate analogues, allylicphosphonate (cPEP) and sulfoenolpyruvate (SEP), recognize the mesophilic (eEIC) and thermophilic (tEIC) homologues of the receptor domain of bacterial Enzyme I, which has been proposed as a target for antimicrobial research. Chemical Shift Perturbation (CSP) experiments show that cPEP and SEP recognize tEIC over the mesophilic homologue. Combined Principal Component Analysis of half-microsecond-long MD simulations reveals that incomplete quenching of a breathing motion in the eEIC-ligand complex destabilizes the interaction and makes the investigated substrate analogues selective toward the thermophilic enzyme. Our results indicate that residual protein motions need to be considered carefully when optimizing small molecule inhibitors of EI. In general, our work demonstrates that protein conformational dynamics can be exploited in the rational design and optimization of inhibitors with subfamily selectivity.

A Single Point Mutation Controls the Rate of Interconversion Between the g + and g - Rotamers of the Histidine 189 χ2 Angle That Activates Bacterial Enzyme I for Catalysis

Enzyme I (EI) of the bacterial phosphotransferase system (PTS) is a master regulator of bacterial metabolism and a promising target for development of a new class of broad-spectrum antibiotics. The catalytic activity of EI is mediated by several intradomain, interdomain, and intersubunit conformational equilibria. Therefore, in addition to its relevance as a drug target, EI is also a good model for investigating the dynamics/function relationship in multidomain, oligomeric proteins. Here, we use solution NMR and protein design to investigate how the conformational dynamics occurring within the N-terminal domain (EIN) affect the activity of EI. We show that the rotameric g⁺-to-g⁻ transition of the active site residue His¹⁸⁹ χ2 angle is decoupled from the state A-to-state B transition that describes a ∼90° rigid-body rearrangement of the EIN subdomains upon transition of the full-length enzyme to its catalytically competent closed form. In addition, we engineered EIN constructs with modulated conformational dynamics by hybridizing EIN from mesophilic and thermophilic species, and used these chimeras to assess the effect of increased or decreased active site flexibility on the enzymatic activity of EI. Our results indicate that the rate of the autophosphorylation reaction catalyzed by EI is independent from the kinetics of the g⁺-to-g⁻ rotameric transition that exposes the phosphorylation site on EIN to the incoming phosphoryl group. In addition, our work provides an example of how engineering of hybrid mesophilic/thermophilic chimeras can assist investigations of the dynamics/function relationship in proteins, therefore opening new possibilities in biophysics.

Structure elucidation of the elusive Enzyme I monomer reveals the molecular mechanisms linking oligomerization and enzymatic activity

Enzyme I (EI) is a phosphotransferase enzyme responsible for converting phosphoenolpyruvate (PEP) into pyruvate. This reaction initiates a five-step phosphorylation cascade in the bacterial phosphotransferase (PTS) transduction pathway. Under physiological conditions, EI exists in an equilibrium between a functional dimer and an inactive monomer. The monomer-dimer equilibrium is a crucial factor regulating EI activity and the phosphorylation state of the overall PTS. Experimental studies of EI's monomeric state have yet been hampered by the dimer's high thermodynamic stability, which prevents its characterization by standard structural techniques. In this study, we modified the dimerization domain of EI (EIC) by mutating three amino acids involved in the formation of intersubunit salt bridges. The engineered variant forms an active dimer in solution that can bind and hydrolyze PEP. Using hydrostatic pressure as an additional perturbation, we were then able to study the complete dissociation of the variant from 1 bar to 2.5 kbar in the absence and the presence of EI natural ligands. Backbone residual dipolar couplings collected under high-pressure conditions allowed us to determine the conformational ensemble of the isolated EIC monomeric state in solution. Our calculations reveal that three catalytic loops near the dimerization interface become unstructured upon monomerization, preventing the monomeric enzyme from binding its natural substrate. This study provides an atomic-level characterization of EI's monomeric state and highlights the role of the catalytic loops as allosteric connectors controlling both the activity and oligomerization of the enzyme.

An allosteric pocket for inhibition of bacterial Enzyme I identified by NMR-based fragment screening

Enzyme I (EI), which is the key enzyme to activate the bacterial phosphotransferase system, plays an important role in the regulation of several metabolic pathways and controls the biology of bacterial cells at multiple levels. The conservation and ubiquity of EI among different types of bacteria makes the enzyme a potential target for antimicrobial research. Here, we use NMR-based fragment screening to identify novel inhibitors of EI. We identify three molecular fragments that allosterically inhibit the phosphoryl transfer reaction catalyzed by EI by interacting with the enzyme at a surface pocket located more than 10 Å away from the substrate binding site. Interestingly, although the three molecules share the same binding pocket, we observe that two of the discovered EI ligands act as competitive inhibitors while the third ligand acts as a mixed inhibitor. Characterization of the EI-inhibitor complexes by NMR and Molecular Dynamics simulations reveals key interactions that perturb the fold of the active site and provides structural foundation for the different inhibitory activity of the identified molecular fragments. In particular, we show that contacts between the inhibitor and the side-chain of V292 are crucial to destabilize binding of the substrate to EI. In contrast, mixed inhibition is caused by additional contacts between the inhibitor and ⍺-helix 2 that perturb the active site structure and turnover in an allosteric manner. We expect our results to provide the basis for the development of second generation allosteric inhibitors of increased potency and to suggest novel molecular strategies to combat drug-resistant infections.

Hybrid Thermophilic/Mesophilic Enzymes Reveal a Role for Conformational Disorder in Regulation of Bacterial Enzyme I

Conformational disorder is emerging as an important feature of biopolymers, regulating a vast array of cellular functions, including signaling, phase separation, and enzyme catalysis. Here we combine NMR, crystallography, computer simulations, protein engineering, and functional assays to investigate the role played by conformational heterogeneity in determining the activity of the C-terminal domain of bacterial Enzyme I (EIC). In particular, we design chimeric proteins by hybridizing EIC from thermophilic and mesophilic organisms, and we characterize the resulting constructs for structure, dynamics, and biological function. We show that EIC exists as a mixture of active and inactive conformations and that functional regulation is achieved by tuning the thermodynamic balance between active and inactive states. Interestingly, we also present a hybrid thermophilic/mesophilic enzyme that is thermostable and more active than the wild-type thermophilic enzyme, suggesting that hybridizing thermophilic and mesophilic proteins is a valid strategy to engineer thermostable enzymes with significant low-temperature activity.

NMR Methods for Structural Characterization of Protein-Protein Complexes

Protein-protein interactions and the complexes thus formed are critical elements in a wide variety of cellular events that require an atomic-level description to understand them in detail. Such complexes typically constitute challenging systems to characterize and drive the development of innovative biophysical methods. NMR spectroscopy techniques can be applied to extract atomic resolution information on the binding interfaces, intermolecular affinity, and binding-induced conformational changes in protein-protein complexes formed in solution, in the cell membrane, and in large macromolecular assemblies. Here we discuss experimental techniques for the characterization of protein-protein complexes in both solution NMR and solid-state NMR spectroscopy. The approaches include solvent paramagnetic relaxation enhancement and chemical shift perturbations (CSPs) for the identification of binding interfaces, and the application of intermolecular nuclear Overhauser effect spectroscopy and residual dipolar couplings to obtain structural constraints of protein-protein complexes in solution. Complementary methods in solid-state NMR are described, with emphasis on the versatility provided by heteronuclear dipolar recoupling to extract intermolecular constraints in differentially labeled protein complexes. The methods described are of particular relevance to the analysis of membrane proteins, such as those involved in signal transduction pathways, since they can potentially be characterized by both solution and solid-state NMR techniques, and thus outline key developments in this frontier of structural biology.

The oligomerization state of bacterial enzyme I (EI) determines EI's allosteric stimulation or competitive inhibition by α-ketoglutarate

The bacterial phosphotransferase system (PTS) is a signal transduction pathway that couples phosphoryl transfer to active sugar transport across the cell membrane. The PTS is initiated by phosphorylation of enzyme I (EI) by phosphoenolpyruvate (PEP). The EI phosphorylation state determines the phosphorylation states of all other PTS components and is thought to play a central role in the regulation of several metabolic pathways and to control the biology of bacterial cells at multiple levels, for example, affecting virulence and biofilm formation. Given the pivotal role of EI in bacterial metabolism, an improved understanding of the mechanisms controlling its activity could inform future strategies for bioengineering and antimicrobial design. Here, we report an enzymatic assay, based on Selective Optimized Flip Angle Short Transient (SOFAST) NMR experiments, to investigate the effect of the small-molecule metabolite α-ketoglutarate (αKG) on the kinetics of the EI-catalyzed phosphoryl transfer reaction. We show that at experimental conditions favoring the monomeric form of EI, αKG promotes dimerization and acts as an allosteric stimulator of the enzyme. However, when the oligomerization state of EI is shifted toward the dimeric species, αKG functions as a competitive inhibitor of EI. We developed a kinetic model that fully accounted for the experimental data and indicated that bacterial cells might use the observed interplay between allosteric stimulation and competitive inhibition of EI by αKG to respond to physiological fluctuations in the intracellular environment. We expect that the mechanism for regulating EI activity revealed here is common to several other oligomeric enzymes.

The phosphocarrier protein HPr of the bacterial phosphotransferase system globally regulates energy metabolism by directly interacting with multiple enzymes in Escherichia coli

The histidine-phosphorylatable phosphocarrier protein (HPr) is an essential component of the sugar-transporting phosphotransferase system (PTS) in many bacteria. Recent interactome findings suggested that HPr interacts with several carbohydrate-metabolizing enzymes, but whether HPr plays a regulatory role was unclear. Here, we provide evidence that HPr interacts with a large number of proteins in Escherichia coli We demonstrate HPr-dependent allosteric regulation of the activities of pyruvate kinase (PykF, but not PykA), phosphofructokinase (PfkB, but not PfkA), glucosamine-6-phosphate deaminase (NagB), and adenylate kinase (Adk). HPr is either phosphorylated on a histidyl residue (HPr-P) or non-phosphorylated (HPr). PykF is activated only by non-phosphorylated HPr, which decreases the PykF K_half for phosphoenolpyruvate by 10-fold (from 3.5 to 0.36 mm), thus influencing glycolysis. PfkB activation by HPr, but not by HPr-P, resulted from a decrease in the K_half for fructose-6-P, which likely influences both gluconeogenesis and glycolysis. Moreover, NagB activation by HPr was important for the utilization of amino sugars, and allosteric inhibition of Adk activity by HPr-P, but not by HPr, allows HPr to regulate the cellular energy charge coordinately with glycolysis. These observations suggest that HPr serves as a directly interacting global regulator of carbon and energy metabolism and probably of other physiological processes in enteric bacteria.

Biophysical characterization of the domain association between cytosolic A and B domains of the mannitol transporter enzymes II(Mtl) in the presence and absence of a connecting linker

The mannitol transporter enzyme II(Mtl) of the bacterial phosphotransferase system is a multi-domain protein that catalyzes mannitol uptake and phosphorylation. Here we investigated the domain association between cytosolic A and B domains of enzyme II(Mtl) , which are natively connected in Escherichia coli, but separated in Thermoanaerobacter tengcongensis. NMR backbone assignment and residual dipolar couplings indicated that backbone folds were well conserved between the homologous domains. The equilibrium binding of separately expressed domains, however, exhibited ∼28-fold higher affinity compared to the natively linked ones. Phosphorylation of the active site loop significantly contributed to the binding by reducing conformational dynamics at the binding interface, and a few key mutations at the interface were critical to further stabilize the complex by hydrogen bonding and hydrophobic interactions. The affinity increase implicated that domain associations in cell could be maintained at an optimal level regardless of the linker.

Dynamic equilibrium between closed and partially closed states of the bacterial Enzyme I unveiled by solution NMR and X-ray scattering

Enzyme I (EI) is the first component in the bacterial phosphotransferase system, a signal transduction pathway in which phosphoryl transfer through a series of bimolecular protein-protein interactions is coupled to sugar transport across the membrane. EI is a multidomain, 128-kDa homodimer that has been shown to exist in two conformational states related to one another by two large (50-90°) rigid body domain reorientations. The open conformation of apo EI allows phosphoryl transfer from His189 located in the N-terminal domain α/β (EIN(α/β)) subdomain to the downstream protein partner bound to the EIN(α) subdomain. The closed conformation, observed in a trapped phosphoryl transfer intermediate, brings the EIN(α/β) subdomain into close proximity to the C-terminal dimerization domain (EIC), thereby permitting in-line phosphoryl transfer from phosphoenolpyruvate (PEP) bound to EIC to His189. Here, we investigate the solution conformation of a complex of an active site mutant of EI (H189A) with PEP. Simulated annealing refinement driven simultaneously by solution small angle X-ray scattering and NMR residual dipolar coupling data demonstrates unambiguously that the EI(H189A)-PEP complex exists in a dynamic equilibrium between two approximately equally populated conformational states, one corresponding to the closed structure and the other to a partially closed species. The latter likely represents an intermediate in the open-to-closed transition.

Development of a method for reconstruction of crowded NMR spectra from undersampled time-domain data

NMR is a unique methodology for obtaining information about the conformational dynamics of proteins in heterogeneous biomolecular systems. In various NMR methods, such as transferred cross-saturation, relaxation dispersion, and paramagnetic relaxation enhancement experiments, fast determination of the signal intensity ratios in the NMR spectra with high accuracy is required for analyses of targets with low yields and stabilities. However, conventional methods for the reconstruction of spectra from undersampled time-domain data, such as linear prediction, spectroscopy with integration of frequency and time domain, and analysis of Fourier, and compressed sensing were not effective for the accurate determination of the signal intensity ratios of the crowded two-dimensional spectra of proteins. Here, we developed an NMR spectra reconstruction method, "conservation of experimental data in analysis of Fourier" (Co-ANAFOR), to reconstruct the crowded spectra from the undersampled time-domain data. The number of sampling points required for the transferred cross-saturation experiments between membrane proteins, photosystem I and cytochrome b 6 f, and their ligand, plastocyanin, with Co-ANAFOR was half of that needed for linear prediction, and the peak height reduction ratios of the spectra reconstructed from truncated time-domain data by Co-ANAFOR were more accurate than those reconstructed from non-uniformly sampled data by compressed sensing.

Large interdomain rearrangement triggered by suppression of micro- to millisecond dynamics in bacterial Enzyme I

Enzyme I (EI), the first component of the bacterial phosphotransfer signal transduction system, undergoes one of the largest substrate-induced interdomain rearrangements documented to date. Here, we characterize the perturbations generated by two small molecules, the natural substrate phosphoenolpyruvate (PEP) and the inhibitor α-ketoglutarate (αKG), on the structure and dynamics of EI using NMR, small-angle X-ray scattering (SAXS) and biochemical techniques. The results indicate unambiguously that the open-to-closed conformational switch of EI is triggered by complete suppression of micro- to millisecond dynamics within the C-terminal domain of EI. Indeed, we show that a ligand-induced transition from a dynamic to a more rigid conformational state of the C-terminal domain stabilizes the interface between the N- and C-terminal domains observed in the structure of the closed state, thereby promoting the resulting conformational switch and autophosphorylation of EI. The mechanisms described here may be common to several other multidomain proteins and allosteric systems.

Encounter complexes and dimensionality reduction in protein-protein association

An outstanding challenge has been to understand the mechanism whereby proteins associate. We report here the results of exhaustively sampling the conformational space in protein–protein association using a physics-based energy function. The agreement between experimental intermolecular paramagnetic relaxation enhancement (PRE) data and the PRE profiles calculated from the docked structures shows that the method captures both specific and non-specific encounter complexes. To explore the energy landscape in the vicinity of the native structure, the nonlinear manifold describing the relative orientation of two solid bodies is projected onto a Euclidean space in which the shape of low energy regions is studied by principal component analysis. Results show that the energy surface is canyon-like, with a smooth funnel within a two dimensional subspace capturing over 75% of the total motion. Thus, proteins tend to associate along preferred pathways, similar to sliding of a protein along DNA in the process of protein-DNA recognition.

DOI:http://dx.doi.org/10.7554/eLife.01370.001

Proteins rarely act alone. Instead, they tend to bind to other proteins to form structures known as complexes. When two proteins come together to form a complex, they twist and turn through a series of intermediate states before they form the actual complex. These intermediate states are difficult to study because they don’t last for very long, which means that our knowledge of how complexes are formed remains incomplete.

One promising approach for studying the formation of complexes is called paramagnetic relaxation enhancement. In this technique certain areas in one of the proteins are labelled with magnetic particles, which produce signals when the two proteins are close to each other. Repeating the measurement several times with the magnetic particles in different positions provides information about the overall structure of the complex. Computational modelling can then be used to work out the fine details of the structure, including the shapes of the intermediate structures made by the proteins as they interact.

A computer method called docking can be used to predict the most favourable positions that the proteins can take, relative to one another, in a complex. This involves calculating the energy contained in the system, with the correct structure having the lowest energy. Docking methods also predict protein models with slightly higher energies, but with structures that are radically different. Modellers usually ignore these structures, but comparing the docking results to paramagnetic relaxation enhancement data, Kozakov et al. found that these structures actually represent the intermediate states.

Analysing the structure of the intermediate states revealed that the movement of the two proteins relative to one another is severely restricted as they form the final complex. Kozakov et al. found that proteins associate along preferred pathways, similar to the way a protein slides along DNA in the process of protein-DNA recognition. Knowing that the movement of the proteins is restricted in this way will enable researchers to improve the efficiency of docking calculations.

DOI:http://dx.doi.org/10.7554/eLife.01370.002

Thermodynamic dissection of large-scale domain motions coupled with ligand binding of enzyme I

Domain motions are central to the biological functions of many proteins. The energetics of the motions, however, is often difficult to characterize when motions are coupled with the ligand binding. Here, we determined the thermodynamic parameters of individual domain motions and ligand binding of enzyme I (EI) using strategic domain-deletion mutants that selectively removed particular motions. Upon ligand binding, EI employs two large-scale domain motions, the hinge motion and the swivel motion, to switch between conformational states of distinct domain-domain orientations. Calorimetric analysis of the EI mutants separated the free energy changes of the binding and motions, demonstrating that the unfavorable hinge motion (ΔG = 1.5 kcal mol(-1)) was driven by the favorable swivel motion (ΔG = -5.2 kcal mol(-1)). The large free energy differences could be explained by the physicochemical nature of the domain interfaces associated with the motions; the hinge motion employed much narrower interface than the swivel motion without any hydrogen bonds or salt bridges. The small heat capacity further suggested that the packing of the domain interfaces associated with the hinge motion was less compact than that commonly observed in proteins. Lastly, thermodynamic analysis of phosphorylated EI suggests that the domain motions are regulated by the ligand binding and the phosphorylation states. Taken together, the thermodynamic dissection approach illustrates how multiple motions and ligand binding are energetically connected during the functional cycle of EI.

Structural basis for enzyme I inhibition by α-ketoglutarate

Creating new bacterial strains in which carbon and nitrogen metabolism are uncoupled is potentially very useful for optimizing yields of microbial produced chemicals from renewable carbon sources. However, the mechanisms that balance carbon and nitrogen consumption in bacteria are poorly understood. Recently, α-ketoglutarate (αKG), the carbon substrate for ammonia assimilation, has been observed to inhibit Escherichia coli enzyme I (EI), the first component of the bacterial phosphotransferase system (PTS), thereby providing a direct biochemical link between central carbon and nitrogen metabolism. Here we investigate the EI-αKG interaction by NMR and enzymatic assays. We show that αKG binds with a KD of ∼2.2 mM at the active site of EI, acting as a competitive inhibitor. In addition, we use molecular docking simulations to derive a structural model of the enzyme-inhibitor complex that is fully consistent with NMR and analytical ultracentrifugation data. We expect that the EI-αKG structure presented here will provide a starting point for structure-based design of EI mutants resistant to αKG.

Cell signaling, post-translational protein modifications and NMR spectroscopy

Post-translationally modified proteins make up the majority of the proteome and establish, to a large part, the impressive level of functional diversity in higher, multi-cellular organisms. Most eukaryotic post-translational protein modifications (PTMs) denote reversible, covalent additions of small chemical entities such as phosphate-, acyl-, alkyl- and glycosyl-groups onto selected subsets of modifiable amino acids. In turn, these modifications induce highly specific changes in the chemical environments of individual protein residues, which are readily detected by high-resolution NMR spectroscopy. In the following, we provide a concise compendium of NMR characteristics of the main types of eukaryotic PTMs: serine, threonine, tyrosine and histidine phosphorylation, lysine acetylation, lysine and arginine methylation, and serine, threonine O-glycosylation. We further delineate the previously uncharacterized NMR properties of lysine propionylation, butyrylation, succinylation, malonylation and crotonylation, which, altogether, define an initial reference frame for comprehensive PTM studies by high-resolution NMR spectroscopy.

Stability and binding of the phosphorylated species of the N-terminal domain of enzyme I and the histidine phosphocarrier protein from the Streptomyces coelicolor phosphoenolpyruvate:sugar phosphotransferase system

The phosphotransferase system (PTS) is involved in the use of carbon sources in bacteria. It is formed by two general proteins: enzyme I (EI) and the histidine phosphocarrier (HPr), and various sugar-specific permeases. EI is formed by two domains, with the N-terminal domain (EIN) being responsible for the binding to HPr. In low-G+C Gram-positive bacteria, HPr becomes phosphorylated not only by phosphoenolpyruvate (PEP) at the active-site histidine, but also by ATP at a serine. In this work, we have characterized: (i) the stability and binding affinities between the active-site-histidine phosphorylated species of HPr and the EIN from Streptomyces coelicolor; and (ii) the stability and binding affinities of the species involving the phosphorylation at the regulatory serine of HPr(sc). Our results show that the phosphorylated active-site species of both proteins are less stable than the unphosphorylated counterparts. Conversely, the Hpr-S47D, which mimics phosphorylation at the regulatory serine, is more stable than wild-type HPr(sc) due to helical N-capping effects, as suggested by the modeled structure of the protein. Binding among the phosphorylated and unphosphorylated species is always entropically driven, but the affinity and the enthalpy vary widely.

Calorimetric and spectroscopic investigation of the interaction between the C-terminal domain of Enzyme I and its ligands

Enzyme I initiates a series of phosphotransfer reactions during sugar uptake in the bacterial phosphotransferase system. Here, we have isolated a stable recombinant C-terminal domain of Enzyme I (EIC) of Escherichia coli and characterized its interaction with the N-terminal domain of Enzyme I (EIN) and also with various ligands. EIC can phosphorylate EIN, but their binding is transient regardless of the presence of phosphoenolpyruvate (PEP). Circular dichroism and NMR indicate that ligand binding to EIC induces changes near aromatic groups but not in the secondary structure of EIC. Binding of PEP to EIC is an endothermic reaction with the equilibrium dissociation constant (K(D) ) of 0.28 mM, whereas binding of the inhibitor oxalate is an exothermic reaction with K(D) of 0.66 mM from calorimetry. The binding thermodynamics of EIC and PEP compared to that of Enzyme I (EI) and PEP reveals that domain-domain motion in EI can contribute as large as ∼-3.2 kcal/mol toward PEP binding.

Conformational selection and substrate binding regulate the monomer/dimer equilibrium of the C-terminal domain of Escherichia coli enzyme I

The bacterial phosphotransferase system (PTS) is a signal transduction pathway that couples phosphoryl transfer to active sugar transport across the cell membrane. The PTS is initiated by the binding of phosphoenolpyruvate (PEP) to the C-terminal domain (EIC) of enzyme I (EI), a highly conserved protein that is common to all sugar branches of the PTS. EIC exists in a dynamic monomer/dimer equilibrium that is modulated by ligand binding and is thought to regulate the overall PTS. Isolation of EIC has proven challenging, and conformational dynamics within the EIC domain during the catalytic cycle are still largely unknown. Here, we present a robust protocol for expression and purification of recombinant EIC from Escherichia coli and show that isolated EIC is capable of hydrolyzing PEP. NMR analysis and residual dipolar coupling measurements indicate that the isolated EIC domain in solution adopts a stable tertiary fold and quaternary structure that is consistent with previously reported crystallographic data. NMR relaxation dispersion measurements indicate that residues around the PEP binding site and in the β3α3 turn (residues 333-366), which is located at the dimer interface, undergo a rapid transition on the sub-millisecond time scale (with an exchange rate constant of ∼1500 s(-1)) between major open (∼97%) and minor closed (∼3%) conformations. Upon PEP binding, the β3α3 turn is effectively locked in the closed state by the formation of salt bridges between the phosphate group of PEP and the side chains of Lys(340) and Arg(358), thereby stabilizing the dimer.

Active site phosphoryl groups in the biphosphorylated phosphotransferase complex reveal dynamics in a millisecond time scale

The N-terminal domain of Enzyme I (EIN) and phosphocarrier HPr can form a biphosphorylated complex when they are both phosphorylated by excess cellular phosphoenolpyruvate. Here we show that the electrostatic repulsion between the phosphoryl groups in the biphosphorylated complex results in characteristic dynamics at the active site in a millisecond time scale. The dynamics is localized to phospho-His15 and the stabilizing backbone amide groups of HPr, and does not impact on the phospho-His189 of EIN. The dynamics occurs with the k(ex) of ~500 s(-1) which compares to the phosphoryl transfer rate of ~850 s(-1) between EIN and HPr. The conformational dynamics in HPr may be important for its phosphotransfer reactions with multiple partner proteins.

Solution structure of the IIAChitobiose-HPr complex of the N,N'-diacetylchitobiose branch of the Escherichia coli phosphotransferase system

The solution structure of the complex of enzyme IIA of the N,N'-diacetylchitobiose (Chb) transporter with the histidine phosphocarrier protein HPr has been solved by NMR. The IIA(Chb)-HPr complex completes the structure elucidation of representative cytoplasmic complexes for all four sugar branches of the bacterial phosphoryl transfer system (PTS). The active site His-89 of IIA(Chb) was mutated to Glu to mimic the phosphorylated state. IIA(Chb)(H89E) and HPr form a weak complex with a K(D) of ~0.7 mM. The interacting binding surfaces, concave for IIA(Chb) and convex for HPr, complement each other in terms of shape, residue type, and charge distribution, with predominantly hydrophobic residues, interspersed by some uncharged polar residues, located centrally, and polar and charged residues at the periphery. The active site histidine of HPr, His-15, is buried within the active site cleft of IIA(Chb) formed at the interface of two adjacent subunits of the IIA(Chb) trimer, thereby coming into close proximity with the active site residue, H89E, of IIA(Chb). A His89-P-His-15 pentacoordinate phosphoryl transition state can readily be modeled without necessitating any significant conformational changes, thereby facilitating rapid phosphoryl transfer. Comparison of the IIA(Chb)-HPr complex with the IIA(Chb)-IIB(Chb) complex, as well as with other cytoplasmic complexes of the PTS, highlights a unifying mechanism for recognition of structurally diverse partners. This involves generating similar binding surfaces from entirely different underlying structural elements, large interaction surfaces coupled with extensive redundancy, and side chain conformational plasticity to optimize diverse sets of intermolecular interactions.

Domain cooperativity in multidomain proteins: what can we learn from molecular alignment in anisotropic media?

Many proteins have modular design with multiple globular domains connected via flexible linkers. As a simple model of such system, we study a tandem construct consisting of two identical SH3 domains and a variable-length Gly/Ser linker. When the linker is short, this construct represents a dumbbell-shaped molecule with limited amount of domain-domain mobility. Due to its elongated shape, this molecule efficiently aligns in steric alignment media. As the length of the linker increases, the two domains become effectively uncoupled and begin to behave as independent entities. Consequently, their degree of alignment drops, approaching that found in the (near-spherical) isolated SH3 domains. To model the dependence of alignment parameters on the length of the interdomain linker, we have generated in silico a series of conformational ensembles representing SH3 tandems with different linker length. These ensembles were subsequently used as input for alignment prediction software PALES. The predicted alignment tensors were compared with the results of experimental measurements using a series of tandem-SH3 samples in PEG/hexanol alignment media. This comparison broadly confirmed the expected trends. At the same time, it has been found that the isolated SH3 domain aligns much stronger than expected. This finding can be attributed to complex morphology of the PEG/hexanol media and/or to weak site-specific interactions between the protein and the media. In the latter case, there are strong indications that electrostatic interactions may play a role. The fact that PEG/hexanol does not behave as a simple steric media should serve as a caution for studies that use PALES as a quantitative prediction tool (especially for disordered proteins). Further progress in this area depends on our ability to accurately model the anisotropic media and its site-specific interactions with protein molecules. Once this ability is improved, it should be possible to use the alignment parameters as a measure of domain-domain cooperativity, thus identifying the situations where two domains transiently interact with each other or become coupled through a partially structured linker.

Solution structure of the 128 kDa enzyme I dimer from Escherichia coli and its 146 kDa complex with HPr using residual dipolar couplings and small- and wide-angle X-ray scattering

The solution structures of free Enzyme I (EI, ∼128 kDa, 575 × 2 residues), the first enzyme in the bacterial phosphotransferase system, and its complex with HPr (∼146 kDa) have been solved using novel methodology that makes use of prior structural knowledge (namely, the structures of the dimeric EIC domain and the isolated EIN domain both free and complexed to HPr), combined with residual dipolar coupling (RDC), small- (SAXS) and wide- (WAXS) angle X-ray scattering and small-angle neutron scattering (SANS) data. The calculational strategy employs conjoined rigid body/torsion/Cartesian simulated annealing, and incorporates improvements in calculating and refining against SAXS/WAXS data that take into account complex molecular shapes in the description of the solvent layer resulting in a better representation of the SAXS/WAXS data. The RDC data orient the symmetrically related EIN domains relative to the C(2) symmetry axis of the EIC dimer, while translational, shape, and size information is provided by SAXS/WAXS. The resulting structures are independently validated by SANS. Comparison of the structures of the free EI and the EI-HPr complex with that of the crystal structure of a trapped phosphorylated EI intermediate reveals large (∼70-90°) hinge body rotations of the two subdomains comprising the EIN domain, as well as of the EIN domain relative to the dimeric EIC domain. These large-scale interdomain motions shed light on the structural transitions that accompany the catalytic cycle of EI.

The N-terminal domain of the enzyme I is a monomeric well-folded protein with a low conformational stability and residual structure in the unfolded state

The bacterial phosphoenolpyruvate-dependent sugar phosphotransferase system is a multiprotein complex that phosphorylates and, concomitantly, transports carbohydrates across the membrane into the cell. The first protein of the cascade is a multidomain protein so-called enzyme I (EI). The N-terminal domain of EI from Streptomyces coelicolor, EIN(sc), responsible for the binding to the second protein in the cascade (the histidine phosphocarrier, HPr), was cloned and successfully expressed and purified. We have previously shown that EI(sc) binds to HPr(sc) with smaller affinity than other members of the EI and HPr families [Hurtado-Gómez et al. (2008) Biophys. J., 95, 1336-1348]. We think that the study of the isolated binding HPr(sc) domain, that is EIN(sc), could shed light on the small affinity value measured. Therefore, in this work we present a detailed description of the structural features of the EIN domain, as a first step towards a complete characterization of the molecular recognition process between the two proteins. We show that EIN(sc) is a folded protein, with alpha-helix and beta-sheet structures and also random-coil conformations, as shown by circular dichroism (CD), FTIR and NMR spectroscopies. The acquisition of secondary and tertiary structures, and the burial of hydrophobic regions, occurred concomitantly at acidic pHs, but at very low pH, the domain acquired a molten-globule conformation. The EIN(sc) protein was not very stable, with an apparent conformational free energy change upon unfolding, DeltaG, of 4.1 +/- 0.4 kcal mol(-1), which was pH independent in the range explored (from pH 6.0 to 8.5). The thermal denaturation midpoint, which was also pH invariant, was similar to that measured in the isolated intact EI(sc). Although EIN(sc) shows thermal- and chemical denaturations that seems to follow a two-state mechanism, there is evidence of residual structure in the chemical and thermally unfolded states, as indicated by differential scanning calorimetry and CD measurements.

Solution structure of the IIAChitobiose-IIBChitobiose complex of the N,N'-diacetylchitobiose branch of the Escherichia coli phosphotransferase system

The solution structure of the IIA-IIB complex of the N,N'-diacetylchitobiose (Chb) transporter of the Escherichia coli phosphotransferase system has been solved by NMR. The active site His-89 of IIA(Chb) was mutated to Glu to mimic the phosphorylated state and the active site Cys-10 of IIB(Chb) was substituted by serine to prevent intermolecular disulfide bond formation. Binding is weak with a K(D) of approximately 1.3 mm. The two complementary interaction surfaces are largely hydrophobic, with the protruding active site loop (residues 9-16) of IIB(Chb) buried deep within the active site cleft formed at the interface of two adjacent subunits of the IIA(Chb) trimer. The central hydrophobic portion of the interface is surrounded by a ring of polar and charged residues that provide a relatively small number of electrostatic intermolecular interactions that serve to correctly align the two proteins. The conformation of the active site loop in unphosphorylated IIB(Chb) is inconsistent with the formation of a phosphoryl transition state intermediate because of steric hindrance, especially from the methyl group of Ala-12 of IIB(Chb). Phosphorylation of IIB(Chb) is accompanied by a conformational change within the active site loop such that its path from residues 11-13 follows a mirror-like image relative to that in the unphosphorylated state. This involves a transition of the phi/psi angles of Gly-13 from the right to left alpha-helical region, as well as smaller changes in the backbone torsion angles of Ala-12 and Met-14. The resulting active site conformation is fully compatible with the formation of the His-89-P-Cys-10 phosphoryl transition state without necessitating any change in relative translation or orientation of the two proteins within the complex.

Mechanistic details of a protein-protein association pathway revealed by paramagnetic relaxation enhancement titration measurements

Protein-protein association generally proceeds via the intermediary of a transient, lowly populated, encounter complex ensemble. The mechanism whereby the interacting molecules in this ensemble locate their final stereospecific structure is poorly understood. Further, a fundamental question is whether the encounter complex ensemble is an effectively homogeneous population of nonspecific complexes or whether it comprises a set of distinct structural and thermodynamic states. Here we use intermolecular paramagnetic relaxation enhancement (PRE), a technique that is exquisitely sensitive to lowly populated states in the fast exchange regime, to characterize the mechanistic details of the transient encounter complex interactions between the N-terminal domain of Enzyme I (EIN) and the histidine-containing phosphocarrier protein (HPr), two major bacterial signaling proteins. Experiments were conducted at an ionic strength of 150 mM NaCl to eliminate any spurious nonspecific associations not relevant under physiological conditions. By monitoring the dependence of the intermolecular transverse PRE (Gamma(2)) rates measured on (15)N-labeled EIN on the concentration of paramagnetically labeled HPr, two distinct types of encounter complex configurations along the association pathway are identified and dissected. The first class, which is in equilibrium with and sterically occluded by the specific complex, probably involves rigid body rotations and small translations near or at the active site. In contrast, the second class of encounter complex configurations can coexist with the specific complex to form a ternary complex ensemble, which may help EIN compete with other HPr binding partners in vivo by increasing the effective local concentration of HPr even when the active site of EIN is occupied.

Crystal structure of enzyme I of the phosphoenolpyruvate sugar phosphotransferase system in the dephosphorylated state

The bacterial phosphoenolpyruvate (PEP) sugar phosphotransferase system mediates sugar uptake and controls the carbon metabolism in response to carbohydrate availability. Enzyme I (EI), the first component of the phosphotransferase system, consists of an N-terminal protein binding domain (EIN) and a C-terminal PEP binding domain (EIC). EI transfers phosphate from PEP by double displacement via a histidine residue on EIN to the general phosphoryl carrier protein HPr. Here we report the 2.4 A crystal structure of the homodimeric EI from Staphylococcus aureus. EIN consists of the helical hairpin HPr binding subdomain and the phosphorylatable betaalpha phospho-histidine (P-His) domain. EIC folds into an (betaalpha)(8) barrel. The dimer interface of EIC buries 1833 A(2) of accessible surface per monomer and contains two Ca(2+) binding sites per dimer. The structures of the S. aureus and Escherichia coli EI domains (Teplyakov, A., Lim, K., Zhu, P. P., Kapadia, G., Chen, C. C., Schwartz, J., Howard, A., Reddy, P. T., Peterkofsky, A., and Herzberg, O. (2006) Proc. Natl. Acad. Sci. U.S.A. 103, 16218-16223) are very similar. The orientation of the domains relative to each other, however, is different. In the present structure the P-His domain is docked to the HPr binding domain in an orientation appropriate for in-line transfer of the phosphate to the active site histidine of the acceptor HPr. In the E. coli structure the phospho-His of the P-His domain projects into the PEP binding site of EIC. In the S. aureus structure the crystallographic temperature factors are lower for the HPr binding domain in contact with the P-His domain and higher for EIC. In the E. coli structure it is the reverse.

Using the experimentally determined components of the overall rotational diffusion tensor to restrain molecular shape and size in NMR structure determination of globular proteins and protein-protein complexes

This paper describes an approach for making use of the components of the experimentally determined rotational diffusion tensor derived from NMR relaxation measurements in macromolecular structure determination. The parameters of the rotational diffusion tensor describe the shape and size of the macromolecule or macromolecular complex, and are therefore complementary to traditional NMR restraints. The structural information contained in the rotational diffusion tensor is not dissimilar to that present in the small-angle region of solution X-ray scattering profiles. We demonstrate the utility of rotational diffusion tensor restraints for protein structure refinement using the N-terminal domain of enzyme I (EIN) as an example and validate the results by solution small-angle X-ray scattering. We also show how rotational diffusion tensor restraints can be used for docking complexes using the dimeric HIV-1 protease and the EIN-HPr complexes as examples. In the former case, the rotational diffusion tensor restraints are sufficient in their own right to determine the position of one subunit relative to another. In the latter case, rotational diffusion tensor restraints complemented by highly ambiguous distance restraints derived from chemical shift perturbation mapping and a hydrophobic contact potential are sufficient to correctly dock EIN to HPr. In each case, the cluster containing the lowest-energy structure corresponds to the correct solution.