On the link between conformational changes,ligand binding and heat capacity

BackgroundConformational changes coupled to ligand binding constitute the structural and energetics basis underlying cooperativity, allostery and, in general, protein regulation. These conformational rearrangements are associated with heat capacity changes. ITC is a unique technique for studying binding interactions because of the simultaneous determination of the binding affinity and enthalpy, and for providing the best estimates of binding heat capacity changes.Scope of reviewStill controversial issues in ligand binding are the discrimination between the “conformational selection model” and the “induced fit model”, and whether or not conformational changes lead to temperature dependent apparent binding heat capacities. The assessment of conformational changes associated with ligand binding by ITC is discussed. In addition, the “conformational selection” and “induced fit” models are reconciled, and discussed within the context of intrinsically (partially) unstructured proteins.Major conclusionsConformational equilibrium is a major contribution to binding heat capacity changes. A simple model may explain both conformational selection and induced fit scenarios. A temperature-independent binding heat capacity does not necessarily indicate absence of conformational changes upon ligand binding. ITC provides information on the energetics of conformational changes associated with ligand binding (and other possible additional coupled equilibria).General significancePreferential ligand binding to certain protein states leads to an equilibrium shift that is reflected in the coupling between ligand binding and additional equilibria. This represents the structural/energetic basis of the widespread dependence of ligand binding parameters on temperature, as well as pH, ionic strength and the concentration of other chemical species. This article is part of a Special Issue entitled Microcalorimetry in the BioSciences — Principles and Applications, edited by Fadi Bou-Abdallah.     

Conformational selection and induced changes along the catalytic cycle of Escherichia coli dihydrofolate reductase

Protein function often involves changes between different conformations. Central questions are how these conformational changes are coupled to the binding or catalytic processes during which they occur, and how they affect the catalytic rates of enzymes. An important model system is the enzyme dihydrofolate reductase (DHFR) from Escherichia coli, which exhibits characteristic conformational changes of the active‐site loop during the catalytic step and during unbinding of the product. In this article, we present a general kinetic framework that can be used (1) to identify the ordering of events in the coupling of conformational changes, binding, and catalysis and (2) to determine the rates of the substeps of coupled processes from a combined analysis of nuclear magnetic resonance R₂ relaxation dispersion experiments and traditional enzyme kinetics measurements. We apply this framework to E. coli DHFR and find that the conformational change during product unbinding follows a conformational‐selection mechanism, that is, the conformational change occurs predominantly prior to unbinding. The conformational change during the catalytic step, in contrast, is an induced change, that is, the change occurs after the chemical reaction. We propose that the reason for these conformational changes, which are absent in human and other vertebrate DHFRs, is robustness of the catalytic rate against large pH variations and changes to substrate/product concentrations in E. coli. Proteins 2012;. © 2012 Wiley Periodicals, Inc.     

Hydrolases on silica surfaces: Coverage-activity–molecular property relationships revealed

A set of recommendations to maintain high activity of immobilized enzymes is developed based on direct observation via AFM. This helps to close knowledge gaps that often lead to poor performance of nanobiocatalysts for chemical synthesis. Molecule-level height and volume distribution analyses from high-resolution AFM images were applied to Candida antarctica Lipase B (CALB), subtilisin Carlsberg, and the Lipase from Thermomyces lanuginosus (TLL) deposited on model silica surfaces. Ensembles of flexible or “soft” enzymes appear separated when interactions with the surface are considerable at low surface coverage but form highly entangled structures of increased conformational stability at high surface coverage. By contrast, ensembles of rigid or “hard” enzymes appear to maintain stable aggregates even under strong interaction with the surface. The more rigid the enzyme the higher its tendency to remain in a densely packed state that is able to withstand surface-induced conformational transitions detrimental to catalysis. Weakening of surface-protein interactions for “soft” enzymes will prevent single-molecule immobilization, which reduces catalytic competency through structural changes. Multi-layer coverage in enzyme immobilization should generally be avoided due to mass transfer limitations.     

Conformational selection in protein binding and function

Protein binding and function often involves conformational changes. Advanced nuclear magnetic resonance (NMR) experiments indicate that these conformational changes can occur in the absence of ligand molecules (or with bound ligands), and that the ligands may “select” protein conformations for binding (or unbinding). In this review, we argue that this conformational selection requires transition times for ligand binding and unbinding that are small compared to the dwell times of proteins in different conformations, which is plausible for small ligand molecules. Such a separation of timescales leads to a decoupling and temporal ordering of binding/unbinding events and conformational changes. We propose that conformational‐selection and induced‐change processes (such as induced fit) are two sides of the same coin, because the temporal ordering is reversed in binding and unbinding direction. Conformational‐selection processes can be characterized by a conformational excitation that occurs prior to a binding or unbinding event, while induced‐change processes exhibit a characteristic conformational relaxation that occurs after a binding or unbinding event. We discuss how the ordering of events can be determined from relaxation rates and effective on‐ and off‐rates determined in mixing experiments, and from the conformational exchange rates measured in advanced NMR or single‐molecule fluorescence resonance energy transfer experiments. For larger ligand molecules such as peptides, conformational changes and binding events can be intricately coupled and exhibit aspects of conformational‐selection and induced‐change processes in both binding and unbinding direction.     

How conformational changes can affect catalysis,inhibition and drug resistance of enzymes with induced-fit binding mechanism such as the HIV-1 protease

A central question is how the conformational changes of proteins affect their function and the inhibition of this function by drug molecules. Many enzymes change from an open to a closed conformation upon binding of substrate or inhibitor molecules. These conformational changes have been suggested to follow an induced-fit mechanism in which the molecules first bind in the open conformation in those cases where binding in the closed conformation appears to be sterically obstructed such as for the HIV-1 protease. In this article, we present a general model for the catalysis and inhibition of enzymes with induced-fit binding mechanism. We derive general expressions that specify how the overall catalytic rate of the enzymes depends on the rates for binding, for the conformational changes, and for the chemical reaction. Based on these expressions, we analyze the effect of mutations that mainly shift the conformational equilibrium on catalysis and inhibition. If the overall catalytic rate is limited by product unbinding, we find that mutations that destabilize the closed conformation relative to the open conformation increase the catalytic rate in the presence of inhibitors by a factor exp(ΔΔG_C/RT) where ΔΔG_C is the mutation-induced shift of the free-energy difference between the conformations. This increase in the catalytic rate due to changes in the conformational equilibrium is independent of the inhibitor molecule and, thus, may help to understand how non-active-site mutations can contribute to the multi-drug-resistance that has been observed for the HIV-1 protease. A comparison to experimental data for the non-active-site mutation L90M of the HIV-1 protease indicates that the mutation slightly destabilizes the closed conformation of the enzyme. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.     

The association reaction of collagen model polypeptides (Pro-Pro-Gly)n

The characterization of recently synthesized (Pro-Pro-Gly)_n, n = 7, 8 is described, along with melting profile studies of its association equilibrium, and thermal quenching studies of the kinetics of its association reaction. The order of the kinetic reaction is about 3, implying that three peptide chains are involved in the activated state of the rate-limiting step. The reaction rate was found to exhibit a negative temperature coefficient. With the (Pro-Pro-Gly)₇ peptide, the concentration dependence of the (Pro-Pro-Gly)_n association equilibrium was observed for the first time. Detailed thermodynamic analysis for these n = 7, 8 data, together with literature data for n = 10, 15, 20 were carried out for both the simple “all-or-none” binding model and for a series of complex equilibrium models. For the latter, all of the (Pro-Pro-Gly)_n data (in 10% acetic acid) are fit best with a maximally cooperative near-neighbor model with a standard enthalpy change ΔH = ?650 cal/mole of residues, and a standard entropy change ΔS = ?14.63 ?10/n cal/deg-mole of residues, wherein the ?10 eu represents an end-effect contribution to the binding free energy. With regard to optical rotatory properties and thermodynamic parameters, the data for the new n = 7, 8 peptides match rather well with the literature data for the n = 10, 15, and 20 peptides. The enthalpic stabilization per residue of the triple-helical form of (Pro-Pro-Gly)_n was nearly an order of magnitude smaller than the enthalpic stabilization per additional proline obtained from direct calorimetric measurements on native collagens of different (and much lower) proline contents by Privalov and Tiktopulo. [Biopolymers (1970) 9 , 127–139.] Possible explanations for this phenomenon are discussed.     

Building-block architecture of botulinum toxin complex: Conformational changes provide insights into the hemagglutination ability of the complex

Clostridium botulinum produces the botulinum neurotoxin (BoNT). Previously, we provided evidence for the “building-block” model of botulinum toxin complex (TC). In this model, a single BoNT is associated with a single nontoxic nonhemagglutinin (NTNHA), yielding M-TC; three HA-70 molecules are attached and form M-TC/HA-70, and one to three “arms” of the HA-33/HA-17 trimer (two HA-33 and one HA-17) further bind to M-TC/HA-70 via HA-17 and HA-70 binding, yielding one-, two-, and three-arm L-TC. Of all TCs, only the three-arm L-TC caused hemagglutination. In this study, we determined the solution structures for the botulinum TCs using small-angle X-ray scattering (SAXS). The mature three-arm L-TC exhibited the shape of a “bird spreading its wings”, in contrast to the model having three “arms”, as revealed by transmission electron microscopy. SAXS images indicated that one of the three arms of the HA-33/HA-17 trimer bound to both HA-70 and BoNT. Taken together, these findings regarding the conformational changes in the building-block architecture of TC may explain why only three-arm L-TC exhibited hemagglutination.     

Binding of alcohols to the peptide CO group of poly-L-proline in the I and II conformation. II. Binding constants from infrared measurements in solution

Trifluoroethanol, benzyl alcohol, and n-butanol bind to the peptide and acelyl CO groups of poly-O-acetyl-L -hydroxyproline in dichloromethane via hydrogen bonds. The binding aflinity decreases from trifhioroelhanol to n-buitanol. For the acelyl CO groups the binding does not depend on the conformation of the polymer but for the peptide CO groups the binding constants are larger by a factor of two to five time when it is in the helix II conformation (all peptide bonds trans) than when it assumes the helix I conformation (all peptide bonds cis). This preference is explained by the higher accessibility of the peptide CO groups in the II helix. The small additional energy which results from the preferential binding is sufficient, to induce a complete I → II transition because of the very high cooperativily of the system. The quantitative dependence of the equilibrium constant s for the propagation step of the transition on solvent composition (ratio of trifluoroethanol or benzyl alcohol to n-butanol) is derived from the binding data. It agrees satisfactorily with the empirical relation obtained from a best fit to transition curves of Ganseret al. The I ? II conversion of poly-L -proline is therefore an example of a conformational transition whose solvent dependence can be explained by a binding mechanism.     

Specificity in the interaction between poly(L-glutamate) and iron (III) complex ions in aqueous solution

The binding process between sodium poly(L -glutamate) and trans-2,2′,2″,2?-tetrapyridyl-Fe(III) complex ions in aqueous solution at pH around 7 has been studied by means of equilibrium dialysis and optical measurements. The binding isotherm indicates the occurrence of a cooperative process, whereby bound molecules facilitate the association of additional molecules. According to circular dichroism (CD) data, this effect is coupled with that which sees a conformational change in the charged polypeptide upon progessive binding of complex counterions. All these features are discussed in the light of the structural characteristics of the interacting species. A stereochemical model of the association “complex” is proposed.     

Conformational energy analysis of the leucine repeat regions of C/EBP,GCN4, and the proteins of themyc,jun, andfos oncogenes

It has been recently proposed that certain DNA binding proteins (including C/EBP, GCN4 and themyc, jun, andfos oncogene proteins) share a common structural motif based on helix-promoting regions containing heptad repeat sequences of leucines. It has been suggested that this structure is critical to the biological activity of these proteins, since it facilitates the formation of functional dimers held together by interdigitating leucine side-chains along the hydrophobic interfaces between long α-helical regions of the polypeptide chains in a configuration termed the “leucine zipper.” In this paper, conformational energy analysis is used to determine the preferred three-dimensional structures of the leucine repeat regions of these proteins. The results indicate that, in all cases, the global minimum energy conformation for these regions is an amphipathic α-helix with the leucine side-chains arrayed on one side in such a way to favor “leucine zipper” dimerization. Furthermore, amino acid substitutions in these regions (such as Pro for Leu), that are known to inhibit dimer formation and prevent DNA binding, are found to produce significant conformational changes that disrupt the amphipathic helical structure. Thus, these results provide support for the proposed “leucine zipper” configuration as a critical structural feature of this class of DNA binding proteins.

     

Selected-fit versus induced-fit protein binding: kinetic differences and mutational analysis

The binding of a ligand molecule to a protein is often accompanied by conformational changes of the protein. A central question is whether the ligand induces the conformational change (induced-fit), or rather selects and stabilizes a complementary conformation from a pre-existing equilibrium of ground and excited states of the protein (selected-fit). We consider here the binding kinetics in a simple four-state model of ligand-protein binding. In this model, the protein has two conformations, which can both bind the ligand. The first conformation is the ground state of the protein when the ligand is off, and the second conformation is the ground state when the ligand is bound. The induced-fit mechanism corresponds to ligand binding in the unbound ground state, and the selected-fit mechanism to ligand binding in the excited state. We find a simple, characteristic difference between the on- and off-rates in the two mechanisms if the conformational relaxation into the ground states is fast. In the case of selected-fit binding, the on-rate depends on the conformational equilibrium constant, whereas the off-rate is independent. In the case of induced-fit binding, in contrast, the off-rate depends on the conformational equilibrium, while the on-rate is independent. Whether a protein binds a ligand via selected-fit or induced-fit thus may be revealed by mutations far from the protein's binding pocket, or other "perturbations" that only affect the conformational equilibrium. In the case of selected-fit, such mutations will only change the on-rate, and in the case of induced-fit, only the off-rate.     

Evidence for allosteric effects on p53 oligomerization induced by phosphorylation

p53 is a tetrameric protein with a thermodynamically unstable deoxyribonucleic acid (DNA)‐binding domain flanked by intrinsically disordered regulatory domains that control its activity. The unstable and disordered segments of p53 allow high flexibility as it interacts with binding partners and permits a rapid on/off switch to control its function. The p53 tetramer can exist in multiple conformational states, any of which can be stabilized by a particular modification. Here, we apply the allostery model to p53 to ask whether evidence can be found that the “activating” C‐terminal phosphorylation of p53 stabilizes a specific conformation of the protein in the absence of DNA. We take advantage of monoclonal antibodies for p53 that measure indirectly the following conformations: unfolded, folded, and tetrameric. A double antibody capture enzyme linked‐immunosorbent assay was used to observe evidence of conformational changes of human p53 upon phosphorylation by casein kinase 2 in vitro. It was demonstrated that oligomerization and stabilization of p53 wild‐type conformation results in differential exposure of conformational epitopes PAb1620, PAb240, and DO12 that indicates a reduction in the “unfolded” conformation and increases in the folded conformation coincide with increases in its oligomerization state. These data highlight that the oligomeric conformation of p53 can be stabilized by an activating enzyme and further highlight the utility of the allostery model when applied to understanding the regulation of unstable and intrinsically disordered proteins.     

The role of localization in the operation of the mitotic spindle assembly checkpoint

The mitotic spindle assembly checkpoint (^MSAC) is an important regulatory mechanism of the cell cycle, ensuring proper chromosome segregation in mitosis. It delays the transition to anaphase until all chromosomes are properly attached to the mitotic spindle by emitting a diffusible “wait anaphase”-signal from unattached kinetochores. Current models of the checkpoint disregard important spatial properties like localization, diffusion and realistic numbers of kinetochores. To allow for in silico studies of the dynamics of these models in a more realistic environment, we introduce a mathematical framework for quasi-spatial simulation of localized biochemical processes that are typically observed during checkpoint activation and maintenance. The “emitted inhibition” model of the ^MSAC by Doncic et al. (Proc Natl Acad Sci USA 2005; 102:6332–7) assumes instantaneous activation of the diffusible “wait anaphase”-signal upon kinetochore encounter. We modify this model to account for binding kinetics with finite rates and use the developed framework to determine the feasible range of the binding parameters. We find that for proper activation, the binding rate constant has to be fast and above a critical value. Furthermore, this critical value depends significantly on the amount of local binding sites at each kinetochore. The critical values lie in a physiological realistic regime (10⁴–10⁶ M^-1s^-1). We also determine the feasible parameter range for fast checkpoint activation of the “Mad2 template” model, for which the kinetic parameters have recently been studied in vitro by Simonetta et al. (PLoS Biology 2009; 7:1000010). We find critical values for binding and catalysis rate constants, both significantly higher than the measured values. Our results suggest that yet unknown mechanisms at the kinetochores facilitate binding and catalysis in vivo. We conclude that quantitative models of the ^MSAC have to account for the limited availability of binding sites at kinetochores.     

Solid angle as steric parameter of acyclic saturated groups

With the early aim of quantifying steric consequences of chirality, efforts to define a nonempirical steric parameter of chemical groups are reported. Steric hindrance of a reacting center by any acyclic saturated R group has been characterized by a geometric “axial steric parameter”: the solid angle of R. When the group is a “symmetric top substituent” (i.e., when all the terminal atoms are equivalent), the solid angle matches the solid angle of a cone envelope of R. The definition of this cone is compared with Tolman's definition of a ligand cone in organometallic complexes. The chemical significance of this parameter is shown by an excellent correlation with the Dubois' experimental steric parameter E′_s. Modeling steric repulsion by the cone of R, and correcting solid angles for conformational effects, only 3 empirical coefficients are needed to calculate 33 values of E′_s with less than 10% error. The cone model is suggested to be relevant within the limits of random and independent free rotations about all the bonds in the C–R group. A separation between “axial cone steric hindrance” and other steric effects is proposed. The basic model and the corrections proposed allow the conformational features of esters to be discussed.     

Crystal structure of the coxsackievirus A16 RNA-dependent RNA polymerase elongation complex reveals novel features in motif A dynamics

<正>Dear Editor,Coxsackievirus A16(CV A16)and enterovirus 71(EV71)are currently the two primary causative agents of handfoot-and-mouth disease(HFMD)(Solomon et al.,2010;Mao et al.,2014),threatening health of children worldwide.They both belong to the Enterovirus genus of the     

Relationship between proton–proton nmr coupling constants and substituent electronegativities. III. Conformational analysis of proline rings in solution using a generalized Karplus equation

The relationship between published vicinal proton–proton coupling constants and the pseudorotation properties of the pyrrolidine ring in L -proline, 4-hydroxy-L -proline, 4-fluoro-L -proline, and several linear and cyclic model proline peptides is investigated. Compared to earlier studies, several important improvements are incorporated: (1) a new empirical generalization of the classical Karplus equation is utilized, which allows a valid correction for the effects of electronegativity and orientation of substitutents on ³J_HH; (2) an empirical correlation between proton–proton torsion angles and the pseudorotational parameters P and τ_m is derived; and (3) the best fit of the conformational parameters to the experimental coupling constants is obtained by means of a computerized iterative least-squares procedure. Two pseudorotation ranges were considered, classified as type N (χ₂ positive sign) and type S (χ₂ negative sign). The conformational equilibrium is fully described in terms of four geometrical parameters (P_N, τ_N, P_S, τ_S) and the equilibrium constant K. The present results indicate that, in general, the geometrical properties found in x-ray studies of proline and hydroxyproline residues are well preserved in solution. Several novel features are encountered, however. It is demonstrated that the proline ring occurs in a practically 1:1 conformational equilibrium between well-defined N- and S-type forms. Introduction of an amide group at the C-terminal end has no observable effect on this equilibrium, but the formation of a peptide bond at the imino nitrogen site results in a pronounced, but not exclusive, preference for an S-type form which is roughly 1.1 kcal/mol more stable than its N-type counterpart. The hydroxyproline ring system in neutral or acidic medium displays a pure N-type state, but N-acetylation results in the appearance of a minor (S-type) conformation. Cyclic proline dipeptides similarly exist in a biased conformational equilibrium. The major form (77–88%) corresponds to the N-type conformer observed in the solid state; the minor S-form has not been observed before. In contrast, cyclic hydroxyproline dipeptides display complete conformational purity. Ranges of endocyclic torsion angles deduced for the various classes of pyrrolidine derivatives in solution are presented. Each torsion appears confined to a surprisingly narrow range, comprising about 4°–8° in most cases. In all, the proline ring is far less “floppy” than hitherto assumed.     

Co-Factor Binding Confers Substrate Specificity to Xylose Reductase from Debaryomyces hansenii

Binding of substrates into the active site, often through complementarity of shapes and charges, is central to the specificity of an enzyme. In many cases, substrate binding induces conformational changes in the active site, promoting specific interactions between them. In contrast, non-substrates either fail to bind or do not induce the requisite conformational changes upon binding and thus no catalysis occurs. In principle, both lock and key and induced-fit binding can provide specific interactions between the substrate and the enzyme. In this study, we present an interesting case where cofactor binding pre-tunes the active site geometry to recognize only the cognate substrates. We illustrate this principle by studying the substrate binding and kinetic properties of Xylose Reductase from Debaryomyces hansenii (DhXR), an AKR family enzyme which catalyzes the reduction of carbonyl substrates using NADPH as co-factor. DhXR reduces D-xylose with increased specificity and shows no activity towards “non-substrate” sugars like L-rhamnose. Interestingly, apo-DhXR binds to D-xylose and L-rhamnose with similar affinity (K_d∼5.0–10.0 mM). Crystal structure of apo-DhXR-rhamnose complex shows that L-rhamnose is bound to the active site cavity. L-rhamnose does not bind to holo-DhXR complex and thus, it cannot competitively inhibit D-xylose binding and catalysis even at 4–5 fold molar excess. Comparison of K_d values with K_m values reveals that increased specificity for D-xylose is achieved at the cost of moderately reduced affinity. The present work reveals a latent regulatory role for cofactor binding which was previously unknown and suggests that cofactor induced conformational changes may increase the complimentarity between D-xylose and active site similar to specificity achieved through induced-fit mechanism.     

Structural basis for binding and transfer of heme in bacterial heme‐acquisition systems

Periplasmic heme‐binding proteins (PBPs) in Gram‐negative bacteria are components of the heme acquisition system. These proteins shuttle heme across the periplasmic space from outer membrane receptors to ATP‐binding cassette (ABC) heme importers located in the inner‐membrane. In the present study, we characterized the structures of PBPs found in the pathogen Burkholderia cenocepacia (BhuT) and in the thermophile Roseiflexus sp. RS‐1 (RhuT) in the heme‐free and heme‐bound forms. The conserved motif, in which a well‐conserved Tyr interacts with the nearby Arg coordinates on heme iron, was observed in both PBPs. The heme was recognized by its surroundings in a variety of manners including hydrophobic interactions and hydrogen bonds, which was confirmed by isothermal titration calorimetry. Furthermore, this study of 3 forms of BhuT allowed the first structural comparison and showed that the heme‐binding cleft of BhuT adopts an “open” state in the heme‐free and 2‐heme‐bound forms, and a “closed” state in the one‐heme‐bound form with unique conformational changes. Such a conformational change might adjust the interaction of the heme(s) with the residues in PBP and facilitate the transfer of the heme into the translocation channel of the importer.     

Mechanism of allosteric propagation across a β‐sheet structure investigated by molecular dynamics simulations

The bacterial adhesin FimH consists of an allosterically regulated mannose‐binding lectin domain and a covalently linked inhibitory pilin domain. Under normal conditions, the two domains are bound to each other, and FimH interacts weakly with mannose. However, under tensile force, the domains separate and the lectin domain undergoes conformational changes that strengthen its bond with mannose. Comparison of the crystallographic structures of the low and the high affinity state of the lectin domain reveals conformational changes mainly in the regulatory inter‐domain region, the mannose binding site and a large β sheet that connects the two distally located regions. Here, molecular dynamics simulations investigated how conformational changes are propagated within and between different regions of the lectin domain. It was found that the inter‐domain region moves towards the high affinity conformation as it becomes more compact and buries exposed hydrophobic surface after separation of the pilin domain. The mannose binding site was more rigid in the high affinity state, which prevented water penetration into the pocket. The large central β sheet demonstrated a soft spring‐like twisting. Its twisting motion was moderately correlated to fluctuations in both the regulatory and the binding region, whereas a weak correlation was seen in a direct comparison of these two distal sites. The results suggest a so called “population shift” model whereby binding of the lectin domain to either the pilin domain or mannose locks the β sheet in a rather twisted or flat conformation, stabilizing the low or the high affinity state, respectively. Proteins 2016; 84:990–1008. © 2016 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.