Ligand binding to a high‐energy partially unfolded protein

The conformational energy landscape of a protein determines populations of all possible conformations of the protein and also determines the kinetics of the conversion between the conformations. Interaction with ligands influences the conformational energy landscapes of proteins and shifts populations of proteins in different conformational states. To investigate the effect of ligand binding on partial unfolding of a protein, we use Escherichia coli dihydrofolate reductase (DHFR) and its functional ligand NADP⁺ as a model system. We previously identified a partially unfolded form of DHFR that is populated under native conditions. In this report, we determined the free energy for partial unfolding of DHFR at varying concentrations of NADP⁺ and found that NADP⁺ binds to the partially unfolded form as well as the native form. DHFR unfolds partially without releasing the ligand, though the binding affinity for NADP⁺ is diminished upon partial unfolding. Based on known crystallographic structures of NADP⁺‐bound DHFR and the model of the partially unfolded protein we previously determined, we propose that the adenosine‐binding domain of DHFR remains folded in the partially unfolded form and interacts with the adenosine moiety of NADP⁺. Our result demonstrates that ligand binding may affect the conformational free energy of not only native forms but also high‐energy non‐native forms.     

Pre-existing soft modes of motion uniquely defined by native contact topology facilitate ligand binding to proteins

Modeling protein flexibility constitutes a major challenge in accurate prediction of protein-ligand and protein-protein interactions in docking simulations. The lack of a reliable method for predicting the conformational changes relevant to substrate binding prevents the productive application of computational docking to proteins that undergo large structural rearrangements. Here, we examine how coarse-grained normal mode analysis has been advantageously applied to modeling protein flexibility associated with ligand binding. First, we highlight recent studies that have shown that there is a close agreement between the large-scale collective motions of proteins predicted by elastic network models and the structural changes experimentally observed upon ligand binding. Then, we discuss studies that have exploited the predicted soft modes in docking simulations. Two general strategies are noted: pregeneration of conformational ensembles that are then utilized as input for standard fixed-backbone docking and protein structure deformation along normal modes concurrent to docking. These studies show that the structural changes apparently "induced" upon ligand binding occur selectively along the soft modes accessible to the protein prior to ligand binding. They further suggest that proteins offer suitable means of accommodating/facilitating the recognition and binding of their ligand, presumably acquired by evolutionary selection of the suitable three-dimensional structure.     

How different are structurally flexible and rigid binding sites? Sequence and structural features discriminating proteins that do and do not undergo conformational change upon ligand binding

Proteins are dynamic molecules and often undergo conformational change upon ligand binding. It is widely accepted that flexible loop regions have a critical functional role in enzymes. Lack of consideration of binding site flexibility has led to failures in predicting protein functions and in successfully docking ligands with protein receptors. Here we address the question: which sequence and structural features distinguish the structurally flexible and rigid binding sites? We analyze high-resolution crystal structures of ligand bound (holo) and free (apo) forms of 41 proteins where no conformational change takes place upon ligand binding, 35 examples with moderate conformational change, and 22 cases where a large conformational change has been observed. We find that the number of residue-residue contacts observed per-residue (contact density) does not distinguish flexible and rigid binding sites, suggesting a role for specific interactions and amino acids in modulating the conformational changes. Examination of hydrogen bonding and hydrophobic interactions reveals that cases that do not undergo conformational change have high polar interactions constituting the binding pockets. Intriguingly, the large, aromatic amino acid tryptophan has a high propensity to occur at the binding sites of examples where a large conformational change has been noted. Further, in large conformational change examples, hydrophobic-hydrophobic, aromatic-aromatic, and hydrophobic-polar residue pair interactions are dominant. Further analysis of the Ramachandran dihedral angles (phi, psi) reveals that the residues adopting disallowed conformations are found in both rigid and flexible cases. More importantly, the binding site residues adopting disallowed conformations clustered narrowly into two specific regions of the L-Ala Ramachandran map. Examination of the dihedral angles changes upon ligand binding shows that the magnitude of phi, psi changes are in general minimal, although some large changes particularly between right-handed alpha-helical and extended conformations are seen. Our work further provides an account of conformational changes in the dihedral angles space. The findings reported here are expected to assist in providing a framework for predicting protein-ligand complexes and for template-based prediction of protein function.     

Small-angle X-ray scattering studies of enzymes

Enzyme function requires conformational changes to achieve substrate binding, domain rearrangements, and interactions with partner proteins, but these movements are difficult to observe. Small-angle X-ray scattering (SAXS) is a versatile structural technique that can probe such conformational changes under solution conditions that are physiologically relevant. Although it is generally considered a low-resolution structural technique, when used to study conformational changes as a function of time, ligand binding, or protein interactions, SAXS can provide rich insight into enzyme behavior, including subtle domain movements. In this perspective, we highlight recent uses of SAXS to probe structural enzyme changes upon ligand and partner-protein binding and discuss tools for signal deconvolution of complex protein solutions.     

Development of a fluorescently labeled thermostable DHFR for studying conformational changes associated with inhibitor binding

A fluorescently-labeled, conformationally-sensitive Bacillus stearothermophilus (Bs) dihydrofolate reductase (DHFR) (C73A/S131C_MDCC DHFR) was developed and used to investigate kinetics and protein conformational motions associated with methotrexate (MTX) binding. This construct bears a covalently-attached fluorophore, N-[2-(1-maleimidyl)ethyl]-7-(diethylamino)coumarin-3-carboxamide (MDCC) attached at a distal cysteine, introduced by mutagenesis. The probe is sensitive to the local molecular environment, reporting on changes in the protein structure associated with ligand binding. Intrinsic tryptophan fluorescence of the unlabeled Bs DHFR construct (C73A/S131C DHFR) also showed changes upon MTX association. Stopped-flow analysis of all data can be understood by invoking the presence of two native state DHFR conformers that bind to MTX at different rates (20.2 and 0.067 μM⁻¹ s⁻¹), similar to previously published findings for Escherichia coli DHFR. Probe fluorescence of C73A/S131C_MDCC DHFR predominantly reports on MTX binding to one of the conformers while intrinsic tryptophan fluorescence of C73A/S131C DHFR reports on binding to the other conformer. This study demonstrates the use of an extrinsic fluorophore attached to a distal region to investigate ligand binding interactions that are not experimentally accessible via intrinsic tryptophan fluorescence alone. The thermostability of C73A/S131C_MDCC DHFR provides an important new tool with applications for investigating the temperature dependence of DHFR conformational changes associated with binding and catalysis.     

What is the relationship between the global structures of apo and holo proteins?

It is well known that ligand binding and release may induce a wide range of structural changes in a receptor protein, varying from small movements of loops or side chains in the binding pocket to large‐scale domain hinge‐bending and shear motions or even partial unfolding that facilitates the capture and release of a ligand. An interesting question is what in general are the conformational changes triggered by ligand binding? The aim of this work is analyze the magnitude of structural changes in a protein resulting from ligand binding to assess if the state of ligand binding needs to be included in template‐based protein structure prediction algorithms. To address this issue, a nonredundant dataset of 521 paired protein structures in the ligand‐free and ligand‐bound form was created and used to estimate the degree of both local and global structure similarity between the apo and holo forms. In most cases, the proteins undergo relatively small conformational rearrangements of their tertiary structure upon ligand binding/release (most root‐mean‐square‐deviations from native, RMSD, are <1 Å). However, a clear difference was observed between single‐ and multiple‐domain proteins. For the latter, RMSD changes greater than 1 Å and sometimes larger were found for almost 1/3 of the cases; these are mainly associated with large‐scale hinge‐bending movements of entire domains. The changes in the mutual orientation of individual domains in multiple‐domain proteins upon ligand binding were investigated using a mechanistic model based on mass‐weighted principal axes as well as interface buried surface calculations. Some preferences toward the anticipated mechanism of protein domain movements are predictable based on the examination of just the ligand‐free structural form. These results have applications to protein structure prediction, particularly in the context of protein domain assembly, if additional information concerning ligand binding is exploited. Proteins 2008. © 2007 Wiley‐Liss, Inc.     

The coupling of structural fluctuations to hydride transfer in dihydrofolate reductase

The energy barrier for hydride transfer in wild-type G121V and G121S variants of Escherichia coli dihydrofolate reductase (DHFR) fluctuates in a time-dependent manner. This fluctuation may be attributed to structural changes in the protein that modulate the site of chemistry. Despite being far from the active site, mutations at position 121 of DHFR reduce the hydride transfer rate of the enzyme. This occurrence has been suggested to arise from modifications to the conformational ensemble of the protein. We elucidate the effects of the G121S and G121V mutations on the hydride transfer barrier by identifying structural changes in the protein that correlate with lowered barriers. The effect of these structural parameters on the hydride transfer barrier may be rationalized by simple considerations of the geometric constraints of the hydride transfer reaction. Fluctuations of these properties are associated with specific backbone dihedral angles of residues within the Methione-20 (M20) loop. The dihedral angle preferences are mediated by interactions with the region of the enzyme in the vicinity of residue 121 and are translated into distinct ligand conformations. We predict mutations within the M20 loop that may alter the conformational space explored by DHFR. Such mutational changes are anticipated to adjust the hydride transfer efficacy of DHFR by modifying equilibrium distributions of hydride transfer barriers found in the enzyme.     

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.     

NMR structures of apo L. casei dihydrofolate reductase and its complexes with trimethoprim and NADPH: contributions to positive cooperative binding from ligand-induced refolding, conformational changes, and interligand hydrophobic interactions

In order to examine the origins of the large positive cooperativity (ΔG(0)(coop) = -2.9 kcal mol(-1)) of trimethoprim (TMP) binding to a bacterial dihydrofolate reductase (DHFR) in the presence of NADPH, we have determined and compared NMR solution structures of L. casei apo DHFR and its binary and ternary complexes with TMP and NADPH and made complementary thermodynamic measurements. The DHFR structures are generally very similar except for the A-B loop region and part of helix B (residues 15-31) which could not be directly detected for L. casei apo DHFR because of line broadening from exchange between folded and unfolded forms. Thermodynamic and NMR measurements suggested that a significant contribution to the cooperativity comes from refolding of apo DHFR on binding the first ligand (up to -0.95 kcals mol(-1) if 80% of A-B loop requires refolding). Comparisons of Cα-Cα distance differences and domain rotation angles between apo DHFR and its complexes indicated that generally similar conformational changes involving domain movements accompany formation of the binary complexes with either TMP or NADPH and that the binary structures are approaching that of the ternary complex as would be expected for positive cooperativity. These favorable ligand-induced structural changes upon binding the first ligand will also contribute significantly to the cooperative binding. A further substantial contribution to cooperative binding results from the proximity of the bound ligands in the ternary complex: this reduces the solvent accessible area of the ligand and provides a favorable entropic hydrophobic contribution (up to -1.4 kcal mol(-1)).     

Effects of ligand binding on the internal dynamics of maltose-binding protein.

Ligand binding to proteins often causes large conformational changes. A typical example is maltose-binding protein (MBP), a member of the family of periplasmic binding proteins of Gram-negative bacteria. Upon binding of maltose, MBP undergoes a large structural change that closes the binding cleft, i.e. the distance between its two domains decreases. In contrast, binding of the larger, nonphysiological ligand beta-cyclodextrin does not result in closure of the binding cleft. We have investigated the dynamic properties of MBP in its different states using time-resolved tryptophan fluorescence anisotropy. We found that the 'empty' protein exhibits strong internal fluctuations that almost vanish upon ligand binding. The measured relaxation times corresponding to internal fluctuations can be interpreted as originating from two types of motion: wobbling of tryptophan side-chains relative to the protein backbone, and orientational fluctuations of entire domains. After binding of a ligand, domain motions are no longer detectable and the fluctuations of some of the tryptophan side-chains become rather restricted. This transformation into a more rigid state is observed upon binding of both ligands, maltose and the larger beta-cyclodextrin. The fluctuations of tryptophan side-chains in direct contact with the ligand, however, are affected in a slightly different way by the two ligands.     

Relative solvent accessible surface area predicts protein conformational changes upon binding

Protein interactions are often accompanied by significant changes in conformation. We have analyzed the relationships between protein structures and the conformational changes they undergo upon binding. Based upon this, we introduce a simple measure, the relative solvent accessible surface area, which can be used to predict the magnitude of binding-induced conformational changes from the structures of either monomeric proteins or bound subunits. Applying this to a large set of protein complexes suggests that large conformational changes upon binding are common. In addition, we observe considerable enrichment of intrinsically disordered sequences in proteins predicted to undergo large conformational changes. Finally, we demonstrate that the relative solvent accessible surface area of monomeric proteins can be used as a simple proxy for protein flexibility. This reveals a powerful connection between the flexibility of unbound proteins and their binding-induced conformational changes, consistent with the conformational selection model of molecular recognition.     

A systematic study of the energetics involved in structural changes upon association and connectivity in protein interaction networks

The study of protein binding mechanisms is a major topic of research in structural biology. Here, we implement a combination of metrics to systematically assess the cost of backbone conformational changes that protein domains undergo upon association. Through the analyses of 2090 unique unbound → bound transitions, from over 12,000 structures, we show that two-thirds of these proteins do not suffer significant structural changes upon binding, and could thus fit the lock-and-key model well. Among the remaining proteins, one-third explores the bound conformation in the unbound state (conformational selection model) and, while most transitions are possible from an energetic perspective, a few do require external help to break the thermodynamic barrier (induced fit model). We also analyze the relationship between conformational transitions and protein connectivity, finding that, in general, domains interacting with many partners undergo smaller changes upon association, and are less likely to freely explore larger conformational changes.     

19F NMR study of the leucine-specific binding protein of Escherichia coli: mutagenesis and assignment of the 5-fluorotryptophan-labeled residues

The Escherichia coli L-leucine receptor is an aqueous protein and the first component in the distinct transport pathway for hydrophobic amino acids. L-leucine binding induces a conformational change, which enables the receptor to dock to the membrane components. To investigate the ligand-induced conformational change and binding properties of this protein, we used (19)F NMR to probe the four tryptophan residues located in the two lobes of the protein. The four tryptophan residues were labeled with 5-fluorotryptophan and assigned by site-directed mutagenesis. The (19)F NMR spectra of the partially ligand free proteins show broadened peaks which sharpen when L-leucine is bound, showing that the labeled wild-type protein and mutants are functional. The titration of L-phenylalanine into the 5-fluorotryptophan labeled wild-type protein shows the presence of closed and open conformers. Urea-induced denaturation studies support the NMR results that the wild-type protein binds L-phenylalanine in a different manner to L-leucine. Our studies showed that the tryptophan to phenylalanine mutations on structural units linked to the binding pocket produce subtle changes in the environment of Trp18 located directly in the binding cleft.     

Electron paramagnetic resonance probing of macromolecules: a comparison of structure-function relationships in chymotrypsinogen, -chymotrypsin and anhydrochymotrypsin

Chymotrypsinogen, chymotrypsin and anhydrochymotrypsin have been covalently spin-labeled by an analog of bromoacetamide, and the latter two proteins have been labeled by an analog of 1-chloro-3-tosylamido-4-phenyl butanone. The electron paramagnetic resonance spectra of the labeled proteins indicate protein conformational changes accompanying (1) activation of the zymogen and (2) the binding of protons and substrates by the native and anhydro enzymes, and tertiary structural differences between these protein forms which are at once informative and predictable. A spin-label linked to the thioether side-chain of methionine 192 in Chymotrypsinogen may be in contact with a hydrophobic surface. This interaction is lost upon zymogen activation with little change in the isotropic rotational freedom of the nitroxide group. The rotational freedom of the group increases sigmoidally with pH; a spectral dependence upon an ionizing group (pK_a = 8.9) is demonstrated. The binding of indole to the labeled enzyme raises the pK_a of the ionizing group to 10.2. A spin-label linked to histidine 57 in chymotrypsin senses both indole binding and pH changes directly; the same label in anhydrochymotrypsin responds directly only to changes in pH. Neither histidine-labeled derivative exhibits enzymic activity. The electron paramagnetic resonance spectra of these two labeled proteins at high pH indicate a decrease in the motional freedom of the spin label. The spectral data show that the conformational state of the labeled zymogen is not similar to the high-pH conformational state of the labeled enzyme. Furthermore, the pH-dependent conformational transition of labeled chymotrypsin requires neither the serine 195 hydroxyl nor the histidine 57 imidazole, since the transition occurs normally in derivatized and chemically modified protein forms. The chemical reactivity of histidine 57 in anhydrochymotrypsin is evaluated and the catalytic activities of two histidine alkylated enzymes are compared.     

Protein substrate binding induces conformational changes in the chaperonin GroEL. A suggested mechanism for unfoldase activity

Chaperonins are molecules that assist proteins during folding and protect them from irreversible aggregation. We studied the chaperonin GroEL and its interaction with the enzyme human carbonic anhydrase II (HCA II), which induces unfolding of the enzyme. We focused on conformational changes that occur in GroEL during formation of the GroEL-HCA II complex. We measured the rate of GroEL cysteine reactivity toward iodo[2-(14)C]acetic acid and found that the cysteines become more accessible during binding of a cysteine free mutant of HCA II. Spin labeling of GroEL with N-(1-oxyl-2,2,5, 5-tetramethyl-3-pyrrolidinyl)iodoacetamide revealed that this additional binding occurred because buried cysteine residues become accessible during HCA II binding. In addition, a GroEL variant labeled with 6-iodoacetamidofluorescein exhibited decreased fluorescence anisotropy upon HCA II binding, which resembles the effect of GroES/ATP binding. Furthermore, by producing cysteine-modified GroEL with the spin label N-(1-oxyl-2,2,5, 5-tetramethyl-3-pyrrolidinyl)iodoacetamide and the fluorescent label 5-((((2-iodoacetyl)amino)ethyl)amino)naphthalene-1-sulfonic acid, we detected increases in spin-label mobility and fluorescence intensity in GroEL upon HCA II binding. Together, these results show that conformational changes occur in the chaperonin as a consequence of protein substrate binding. Together with previous results on the unfoldase activity of GroEL, we suggest that the chaperonin opens up as the substrate protein binds. This opening mechanism may induce stretching of the protein, which would account for reported unfoldase activity of GroEL and might explain how GroEL can actively chaperone proteins larger than HCA II.     

Differentiation of the binding of two ligands to a tetrameric protein with a single symmetric active site by 19F NMR

R67 dihydrofolate reductase (R67 DHFR) is a plasmid‐encoded enzyme that confers resistance to the antibacterial drug trimethoprim. R67 DHFR is a tetramer with a single active site that is unusual as both cofactor and substrate are recognized by symmetry‐related residues. Such promiscuity has limited our previous efforts to differentiate ligand binding by NMR. To address this problem, we incorporated fluorine at positions 4, 5, 6, or 7 of the indole rings of tryptophans 38 and 45 and characterized the spectra to determine which probe was optimal for studying ligand binding. Two resonances were observed for all apo proteins. Unexpectedly, the W45 resonance appeared broad, and truncation of the disordered N‐termini resulted in the appearance of one sharp W45 resonance. These results are consistent with interaction of the N‐terminus with W45. Binding of the cofactor broadened W38 for all fluorine probes, whereas substrate, dihydrofolate, binding resulted in the appearance of three new resonances for 4‐ and 5‐fluoroindole labeled protein and severe line broadening for 6‐ and 7‐fluoroindole R67 DHFR. W45 became slightly broader upon ligand binding. With only two peaks in the ¹⁹F NMR spectra, our data were able to differentiate cofactor and substrate binding to the single, symmetric active site of R67 DHFR and yield binding affinities.     

Characterization of ligand-induced conformational states in the beta 2 adrenergic receptor

Drugs acting at G protein coupled receptors can be classified in biological assays as either agonists, partial agonists, neutral antagonists, or as inverse agonists. Very little is known about the actual molecular events and structural changes that occur in the receptor following ligand binding and during transmission of a signal across the membrane. Therefore, the structural basis for the biological classification of drug action remains unknown. To date, the conformational state of G protein coupled receptors has been inferred from the activity of the effector enzyme modulated by the G protein. We have used two different approaches to monitor conformational changes in beta 2 adrenergic receptor. Fluorescence spectroscopy can be used to directly monitor structural changes in purified beta 2 adrenergic receptor in real-time. The emission from many fluorescent molecules is strongly dependent on the polarity of the environment in which they are located. Thus, fluorescent probes covalently bound to proteins can be used as sensitive indicators of conformational changes and protein-protein interactions. In addition, we examined functional differences between agonists and partial agonists using fusion proteins between wild-type beta 2 receptor or a constitutively active beta 2 receptor mutant and Gs alpha. These receptor-G protein fusion proteins guarantee highly efficient coupling with a defined stoichiometry. The results of these experiments will be discussed in the context of current models of G protein coupled receptor activation.     

Stabilization of the shikimate pathway enzyme dehydroquinase by covalently bound ligand

Reversible binding of a ligand to an enzyme active site can elicit a variety of changes in the protein, such as conformational changes (close to the site of binding or communicated over long distances), changes in the ionization state of surrounding amino acid side chains, changes in the interaction of the target protein with other subunits (or other proteins), or even changes in the thermodynamic stability of the protein. Relatively little attention has been given to studying these effects in proteins to which the ligand has been irreversibly bound, yet this can be a convenient way of studying the effects of ligand binding in the absence of association/dissociation equilibria. We report the dramatic changes which occur to the shikimate pathway enzyme dehydroquinase when ligand is attached to its active site after borohydride reduction of the mechanistically important Schiff's base intermediates. The effects of this modification have been characterized by limited proteolysis, circular dichroism, guanidine hydrochloride denaturation, and differential scanning calorimetry. The conclusions from these studies are that although anchoring the ligand at the active site does not cause a gross change in conformation, it does increase markedly the conformational stability of the protein. This is conclusively established by three separate experiments: 1) the modified protein is completely resistant to proteases, whereas the unmodified protein is very susceptible to proteolysis; 2) the concentration of guanidine hydrochloride required to unfold the ligand-linked dehydroquinase is 3-4-fold greater than that of the unmodified protein; 3) the melting temperature (Tm) of the modified protein is 40 degrees C higher than that of the unmodified protein. These results are a very clear example of the thermodynamic link between ligand binding, conformational stability, and proteolytic susceptibility in vitro and will be a useful system for dissecting the contributions of individual protein-ligand interactions to these parameters.