Global changes in local protein dynamics reduce the entropic cost of carbohydrate binding in the arabinose-binding protein

Protein dynamics make important but poorly understood contributions to molecular recognition phenomena. To address this, we measure changes in fast protein dynamics that accompany the interaction of the arabinose-binding protein (ABP) with its ligand, d-galactose, using NMR relaxation and molecular dynamics simulation. These two approaches present an entirely consistent view of the dynamic changes that occur in the protein backbone upon ligand binding. Increases in the amplitude of motions are observed throughout the protein, with the exception of a few residues in the binding site, which show restriction of dynamics. These counter-intuitive results imply that a localised binding event causes a global increase in the extent of protein dynamics on the pico- to nanosecond timescale. This global dynamic change constitutes a substantial favourable entropic contribution to the free energy of ligand binding. These results suggest that the structure and dynamics of ABP may be adapted to exploit dynamic changes to reduce the entropic costs of binding.     

Water in protein hydration and ligand recognition

This review describes selected basics of water in biomolecular recognition. We focus on a qualitative understanding of the most important physical aspects, how these change in magnitude between bulk water and protein environment, and how the roles that water plays for proteins arise from them. These roles include mechanical support, thermal coupling, dielectric screening, mass and charge transport, and the competition with a ligand for the occupation of a binding site. The presence or absence of water has ramifications that range from the thermodynamic binding signature of a single ligand up to cellular survival. The large inhomogeneity in water density, polarity and mobility around a solute is hard to assess in experiment. This is a source of many difficulties in the solvation of protein models and computational studies that attempt to elucidate or predict ligand recognition. The influence of water in a protein binding site on the experimental enthalpic and entropic signature of ligand binding is still a point of much debate. The strong water‐water interaction in enthalpic terms is counteracted by a water molecule's high mobility in entropic terms. The complete arrest of a water molecule's mobility sets a limit on the entropic contribution of a water displacement process, while the solvent environment sets limits on ligand reactivity.     

Effect of prototypic drugs ibuprofen and warfarin on global chaotropic unfolding of human serum heme-albumin: a fast-field-cycling 1H-NMR relaxometric study

Human serum albumin (HSA) is the most prominent protein in plasma, but it is also found in tissues and secretions throughout the body. The three-domain design of HSA provides a variety of binding sites for many ligands, including heme and drugs. HSA has been used as a model multidomain protein to investigate how interdomain interactions affect the global folding/unfolding process. Here, we report on the reversible chemical denaturation of heme-HSA involving three different conformational states (F, N, and B, occurring at pH 4.0, 7.0, and 9.0, respectively) and on the effect of prototypic drugs ibuprofen and warfarin on thermodynamics of the reversible unfolding process. Chaotropic unfolding of heme-HSA in the F, N, and B conformations is governed by different thermodynamic regimes, with the B form showing an entropic stabilization of the structure that compensates an enthalpic destabilization, and the F form easily unfolding under entropic control. Warfarin and ibuprofen binding stabilizes heme-HSA in both N and B states.     

Substrate-Modulated Thermal Fluctuations Affect Long-Range Allosteric Signaling in Protein Homodimers: Exemplified in CAP

The role of conformational dynamics in allosteric signaling of proteins is increasingly recognized as an important and subtle aspect of this ubiquitous phenomenon. Cooperative binding is commonly observed in proteins with twofold symmetry that bind two identical ligands. We construct a coarse-grained model of an allosteric coupled dimer and show how the signal can be propagated between the distant binding sites via change in slow global vibrational modes alone. We demonstrate that modulation on substrate binding of as few as 5-10 slow modes can give rise to cooperativity observed in biological systems and that the type of cooperativity is given by change of interaction between the two monomers upon ligand binding. To illustrate the application of the model, we apply it to a challenging test case: the catabolite activator protein (CAP). CAP displays negative cooperativity upon association with two identical ligands. The conformation of CAP is not affected by the binding, but its vibrational spectrum undergoes a strong modification. Intriguingly, the first binding enhances thermal fluctuations, yet the second quenches them. We show that this counterintuitive behavior is, in fact, necessary for an optimal anticooperative system, and captured within a well-defined region of the model's parameter space. From analyzing the experimental results, we conclude that fast local modes take an active part in the allostery of CAP, coupled to the more-global slow modes. By including them into the model, we elucidate the role of the modes on different timescales. We conclude that such dynamic control of allostery in homodimers may be a general phenomenon and that our model framework can be used for extended interpretation of thermodynamic parameters in other systems.     

Allosteric Dynamic Control of Binding

Proteins have a highly dynamic nature and there is a complex interrelation between their structural dynamics and binding behavior. By assuming various conformational ensembles, they perform both local and global fluctuations to interact with other proteins in a dynamic infrastructure adapted to functional motion. Here, we show that there is a significant association between allosteric mutations, which lead to high-binding-affinity changes, and the hinge positions of global modes, as revealed by a large-scale statistical analysis of data in the Structural Kinetic and Energetic Database of Mutant Protein Interactions (SKEMPI). We further examined the mechanism of allosteric dynamics by conducting studies on human growth hormone (hGH) and pyrin domain (PYD), and the results show how mutations at the hinge regions could allosterically affect the binding-site dynamics or induce alternative binding modes by modifying the ensemble of accessible conformations. The long-range dissemination of perturbations in local chemistry or physical interactions through an impact on global dynamics can restore the allosteric dynamics. Our findings suggest a mechanism for the coupling of structural dynamics to the modulation of protein interactions, which remains a critical phenomenon in understanding the effect of mutations that lead to functional changes in proteins.     

Streptavidin-biotin binding energetics

The high affinity energetics in the streptavidin-biotin system provide an excellent model system for studying how proteins balance enthalpic and entropic components to generate an impressive overall free energy for ligand binding. We review here concerted site-directed mutagenesis, biophysical, and computational studies of aromatic and hydrogen bonding interaction energetics between streptavidin and biotin. These results also have provided insight into how streptavidin builds a large activation barrier to dissociation by managing the enthalpic and entropic activation components. Finally, we review recent studies of the biotin dissociation pathway that address the fundamental question of how ligands exit protein binding pockets.     

A Direct Coupling between Global and Internal Motions in a Single Domain Protein? MD Investigation of Extreme Scenarios

Proteins are not rigid molecules, but exhibit internal motions on timescales ranging from femto- to milliseconds and beyond. In solution, proteins also experience global translational and rotational motions, sometimes on timescales comparable to those of the internal fluctuations. The possibility that internal and global motions may be directly coupled has intriguing implications, given that enzymes and cell signaling proteins typically associate with binding partners and cellular scaffolds. Such processes alter their global motion and may affect protein function. Here, we present molecular dynamics simulations of extreme case scenarios to examine whether a possible relationship exists. In our model protein, a ubiquitin-like RhoGTPase binding domain of plexin-B1, we removed either internal or global motions. Comparisons with unrestrained simulations show that internal and global motions are not appreciably coupled in this single-domain protein. This lack of coupling is consistent with the observation that the dynamics of water around the protein, which is thought to permit, if not stimulate, internal dynamics, is also largely independent of global motion. We discuss implications of these results for the structure and function of proteins.     

Hsp104 and Potentiated Variants Can Operate as Distinct Nonprocessive Translocases

Heat shock protein (Hsp) 104 is a hexameric ATPases associated with diverse cellular activities motor protein that enables cells to survive extreme stress. Hsp104 couples the energy of ATP binding and hydrolysis to solubilize proteins trapped in aggregated structures. The mechanism by which Hsp104 disaggregates proteins is not completely understood but may require Hsp104 to partially or completely translocate polypeptides across its central channel. Here, we apply transient state, single turnover kinetics to investigate the ATP-dependent translocation of soluble polypeptides by Hsp104 and Hsp104^A503S, a potentiated variant developed to resolve misfolded conformers implicated in neurodegenerative disease. We establish that Hsp104 and Hsp104^A503S can operate as nonprocessive translocases for soluble substrates, indicating a “partial threading” model of translocation. Remarkably, Hsp104^A503S exhibits altered coupling of ATP binding to translocation and decelerated dissociation from polypeptide substrate compared to Hsp104. This altered coupling and prolonged substrate interaction likely increases entropic pulling forces, thereby enabling more effective aggregate dissolution by Hsp104^A503S.     

An entropic mechanism to generate highly cooperative and specific binding from protein phosphorylations

Cooperative interactions are essential to the operation of many biochemical networks. Such networks then respond ultrasensitively in a nonlinear manner to linear changes in network input, and network output, for example, levels of a phosphorylated protein or of gene expression, becomes a sigmoidal function of concentrations of input molecules. We present a novel, entropic ultrasensitivity mechanism that generates highly cooperative and specific binding between two proteins. We consider a disordered protein with multiple phosphorylation sites that binds to a single binding site on an interacting protein. We assume that each phosphorylation locally orders the protein. Such local order affects protein conformational entropy nonlinearly and generates binding that is a highly cooperative function of the number of protein phosphorylations (with Hill coefficients well above 10). Substantial binding may only occur once the disordered protein is phosphorylated a critical number of times or more. Cooperativity is determined by the size of the disordered region of the protein, the binding affinity, and unusually the concentration of the interacting protein. Given the widespread occurrence of disordered, multiply phosphorylated proteins, its highly ultrasensitive character, and the ease of its control, entropic, phosphorylation-driven cooperativity may be extensively exploited intracellularly.     

Dynamic Allostery of the Catabolite Activator Protein Revealed by Interatomic Forces

The Catabolite Activator Protein (CAP) is a showcase example for entropic allostery. For full activation and DNA binding, the homodimeric protein requires the binding of two cyclic AMP (cAMP) molecules in an anti-cooperative manner, the source of which appears to be largely of entropic nature according to previous experimental studies. We here study at atomic detail the allosteric regulation of CAP with Molecular dynamics (MD) simulations. We recover the experimentally observed entropic penalty for the second cAMP binding event with our recently developed force covariance entropy estimator and reveal allosteric communication pathways with Force Distribution Analyses (FDA). Our observations show that CAP binding results in characteristic changes in the interaction pathways connecting the two cAMP allosteric binding sites with each other, as well as with the DNA binding domains. We identified crucial relays in the mostly symmetric allosteric activation network, and suggest point mutants to test this mechanism. Our study suggests inter-residue forces, as opposed to coordinates, as a highly sensitive measure for structural adaptations that, even though minute, can very effectively propagate allosteric signals.     

Immunoglobulin structure exhibits control over CDR motion

Motions of the IgG structure are evaluated using normal mode analysis of an elastic network model to detect hinges, the dominance of low frequency modes, and the most important internal motions. One question we seek to answer is whether or not IgG hinge motions facilitate antigen binding. We also evaluate the protein crystal and packing effects on the experimental temperature factors and disorder predictions. We find that the effects of the protein environment on the crystallographic temperature factors may be misleading for evaluating specific functional motions of IgG. The extent of motion of the antigen binding domains is computed to show their large spatial sampling. We conclude that the IgG structure is specifically designed to facilitate large excursions of the antigen binding domains. Normal modes are shown as capable of computationally evaluating the hinge motions and the spatial sampling by the structure. The antigen binding loops and the major hinge appear to behave similarly to the rest of the structure when we consider the dominance of the low frequency modes and the extent of internal motion. The full IgG structure has a lower spectral dimension than individual F(ab) domains, pointing to more efficient information transfer through the antibody than through each domain. This supports the claim that the IgG structure is specifically constructed to facilitate antigen binding by coupling motion of the antigen binding loops with the large scale hinge motions.     

Wrapping of flanking non-operator DNA in lac repressor-operator complexes: implications for DNA looping

In our studies of lac repressor tetramer (T)-lac operator (O) interactions, we observed that the presence of extended regions of non-operator DNA flanking a single lac operator sequence embedded in plasmid DNA produced large and unusual cooperative and anticooperative effects on binding constants (Kobs) and their salt concentration dependences for the formation of 1:1 (TO) and especially 1:2 (TO2) complexes. To explore the origin of this striking behavior we report and analyze binding data on 1:1 (TO) and 1:2 (TO2) complexes between repressor and a single O(sym) operator embedded in 40 bp, 101 bp, and 2514 bp DNA, over very wide ranges of [salt]. We find large interrelated effects of flanking DNA length and [salt] on binding constants (K(TO)obs, K(TO2)obs) and on their [salt]-derivatives, and quantify these effects in terms of the free energy contributions of two wrapping modes, designated local and global. Both local and global wrapping of flanking DNA occur to an increasing extent as [salt] decreases. Global wrapping of plasmid-length DNA is extraordinarily dependent on [salt]. We propose that global wrapping is driven at low salt concentration by the polyelectrolyte effect, and involves a very large number (>/similar 20) of coulombic interactions between DNA phosphates and positively charged groups on lac repressor. Coulombic interactions in the global wrap must involve both the core and the second DNA-binding domain of lac repressor, and result in a complex which is looped by DNA wrapping. The non-coulombic contribution to the free energy of global wrapping is highly unfavorable ( approximately +30-50 kcal mol(-1)), which presumably results from a significant extent of DNA distortion and/or entropic constraints. We propose a structural model for global wrapping, and consider its implications for looping of intervening non-operator DNA in forming a complex between a tetrameric repressor (LacI) and one multi-operator DNA molecule in vivo and in vitro. The existence of DNA wrapping in LacI-DNA interactions motivates the proposal that most if not all DNA binding proteins may have evolved the capability to wrap and thereby organize flanking regions of DNA.     

The molecular mechanism of toxin-induced conformational changes in a potassium channel: relation to C-type inactivation

Recently, a solid-state NMR study revealed that scorpion toxin binding leads to conformational changes in the selectivity filter of potassium channels. The exact nature of the conformational changes, however, remained elusive. We carried out all-atom molecular dynamics simulations that enabled us to cover the complete pathway of toxin approach and binding, and we validated our simulation results by using solid-state NMR data and electrophysiological measurements. Our structural model revealed a mechanism of cooperative toxin-induced conformational changes that accounts both for the signal changes observed in solid-state NMR and for the tight interaction between KcsA-Kv1.3 and Kaliotoxin. We show that this mechanism is structurally and functionally closely related to recovery from C-type inactivation. Furthermore, our simulations indicate heterogeneity in the binding modes of Kaliotoxin, which might serve to enhance its affinity for KcsA-Kv1.3 further by entropic stabilization.     

Entropic elastic processes in protein mechanisms. II. Simple (passive) and coupled (active) development of elastic forces

The first part of this review on entropic elastic processes in protein mechanisms (Urry, 1988) demonstrated with the polypentapeptide of elastin (Val¹-Pro²-Gly³-Val⁴-Gly⁵)_n that elastic structure develops as the result of an inverse temperature transition and that entropic elasticity is due to internal chain dynamics in a regular nonrandom structure. This demonstration is contrary to the pervasive perspective of entropic protein elasticity of the past three decades wherein a network of random chains has been considered the necessary structural consequence of the occurrence of dominantly entropic elastomeric force. That this is not the case provides a new opportunity for understanding the occurrence and role of entropic elastic processes in protein mechanisms. Entropic elastic processes are considered in two classes: passive and active. The development of elastomeric force on deformation is class I (passive) and the development of elastomeric force as the result of a chemical process shifting the temperature of a transition is class II (active). Examples of class I are elastin, the elastic filament of muscle, elastic force changes in enzyme catalysis resulting from binding processes and resulting in the straining of a scissile bond, and in the turning on and off of channels due to changes in transmembrane potential. Demonstration of the consequences of elastomeric force developing as the result of an inverse temperature transition are seen in elastin, where elastic recoil is lost on oxidation, i.e., on decreasing the hydrophobicity of the chain and shifting the temperature for the development of elastomeric force to temperatures greater than physiological. This is relevant in general to loss of elasticity on aging and more specifically to the development of pulmonary emphysema. Since random chain networks are not the products of inverse temperature transitions and the temperature at which an inverse temperature transition occurs depends on the hydrophobicity of the polypeptide chain, it now becomes possible to consider chemical processes for turning elastomeric force on and off by reversibly changing the hydrophobicity of the polypeptide chain. This is herein called mechanochemical coupling of the first kind; this is the chemical modulation of the temperature for the transition from a less-ordered less elastic state to a more-ordered more elastic state. In the usual considerations to date, development of elastomeric force is the result of a standard transition from a more-ordered less elastic state to a less-ordered more elastic state. When this is chemically modulated, it is herein called mechanochemical coupling of the second kind. For elastin and the polypentapeptide of elastin, since entropic elastomeric force results on formation of a regular nonrandom structure and thermal randomization of chains results in loss of elastic modulus to levels of limited use in protein mechanisms, consideration of regular spiral-like structures rather than ramdom chain networks or random coils are proposed for mechanochemical coupling of the second kind. Chemical processes to effect mechanochemical coupling in biological systems are most obviously phosphorylation-dephosphorylation and changes in calcium ion activity but also changes in pH. These issues are considered in the events attending parturition in muscle contraction and in cell motility.     

Ligand binding to proteins: the binding landscape model.

Models of ligand binding are often based on four assumptions: (1) steric fit: that binding is determined mainly by shape complementarity; (2) native binding: that ligands mainly bind to native states; (3) locality: that ligands perturb protein structures mainly at the binding site; and (4) continuity: that small changes in ligand or protein structure lead to small changes in binding affinity. Using a generalization of the 2D HP lattice model, we study ligand binding and explore these assumptions. We first validate the model by showing that it reproduces typical binding behaviors. We observe ligand-induced denaturation, ANS and heme-like binding, and "lock-and-key" and "induced-fit" specific binding behaviors characterized by Michaelis-Menten or more cooperative types of binding isotherms. We then explore cases where the model predicts violations of the standard assumptions. For example, very different binding modes can result from two ligands of identical shape. Ligands can sometimes bind highly denatured states more tightly than native states and yet have Michaelis-Menten isotherms. Even low-population binding to denatured states can cause changes in global stability, hydrogen-exchange rates, and thermal B-factors, contrary to expectations, but in agreement with experiments. We conclude that ligand binding, similar to protein folding, may be better described in terms of energy landscapes than in terms of simpler mass-action models.     

Molecular basis for co-operativity in Ca2+ binding to calbindin D9k. 1H nuclear magnetic resonance studies of (Cd2+)1-bovine calbindin D9k

The molecular basis for the co-operativity in binding of calcium ions by bovine calbindin D9k has been addressed by carrying out a comparative analysis of the solution conformation and dynamics of the apo, half saturated and fully saturated species using two-dimensional 1H nuclear magnetic resonance spectroscopy. Since the half saturated calcium form of the protein is not significantly populated under equilibrium conditions due to the co-operativity in binding of calcium ions, the half saturated cadmium form of the protein has been substituted for the calcium form. To verify that cadmium forms of calbindin D9k represent viable models for the calcium-bound species, the fully saturated cadmium form has been prepared and compared to the calcium-saturated protein. Virtually complete 1H resonance assignments have been obtained for both the (Cd2+)1 and the (Cd2+)2 states. Secondary structure elements and the global folding pattern were determined from nuclear Overhauser effects, backbone spin-spin coupling constants and slowly exchanging amide protons. Comparisons of the half saturated protein with the apo and calcium-saturated forms of calbindin D9k show that all three structures are highly similar. However, a change in the structural and dynamic properties of the protein does occur upon binding of the first ion; the half saturated form is found to be more similar to the calcium-saturated form than to the apo form. These results have important implications concerning the molecular basis for the co-operativity, and suggest that entropic effects associated with the protein dynamics play an important role.     

Characterization of protein matrix motions in the Rb. sphaeroides photosynthetic reaction center

We use Normal Mode Analysis to investigate motions in the photosynthetic reaction center (RC) protein. We identify the regions involved in concerted fluctuations of the protein matrix and analyze the normalized amplitudes and the directionality of the first few dominant modes. We also seek to quantify the coupling of normal modes to long-range electron transfer (ET). We find that a quasi-continuous spectrum of protein motions rather than one individual mode contributes to light-driven electron transfer. This is consistent with existing theoretical models (e.g. the spin-boson/dispersed polaron model) for the coupling of the protein and solvent "bath" to charge separation events. [Figure: see text].