Allosteric signaling in the biotin repressor occurs via local folding coupled to global dampening of protein dynamics

The biotin repressor is an allosterically regulated, site-specific DNA-binding protein. Binding of the small ligand bio-5′-AMP activates repressor dimerization, which is a prerequisite to DNA binding. Multiple disorder-to-order transitions, some of which are known to be important for the functional allosteric response, occur in the vicinity of the ligand-binding site concomitant with effector binding to the repressor monomer. In this work, the extent to which these local changes are coupled to additional changes in the structure/dynamics of the repressor was investigated using hydrogen/deuterium exchange coupled to mass spectrometry. Measurements were performed on the apo-protein and on complexes of the protein bound to four different effectors that elicit a range of thermodynamic responses in the repressor. Global exchange measurements indicate that binding of any effector to the intact protein is accompanied by protection from exchange. Mass spectrometric analysis of pepsin-cleavage products generated from the exchanged complexes reveals that the protection is distributed throughout the protein. Furthermore, the magnitude of the level of protection in each peptide from hydrogen/deuterium exchange correlates with the magnitude of the functional allosteric response elicited by a ligand. These results indicate that local structural changes in the binding site that occur concomitant with effector binding nucleate global dampening of dynamics. Moreover, the magnitude of dampening of repressor dynamics tracks with the magnitude of the functional response to effector binding.     

Evaluation of energetic and dynamic coupling networks in a PDZ domain protein

A number of computational and experimental studies have identified intramolecular communication "pathways" or "networks" important for transmitting allostery. Here, we have used mutagenesis and NMR relaxation methods to investigate the scope and nature of the communication networks found in the second post-synaptic density-95/discs large/zonula occludens-1 (PDZ) domain of the human protein tyrosine phosphatase 1E protein (hPTP1E) (PDZ2). It was found that most mutations do not have a significant energetic contribution to peptide ligand binding. Three mutants that showed significant changes in binding also displayed context-dependent dynamic effects. Both a mutation at a partially exposed site (H71Y) and a buried core position (I35V) had a limited response in side-chain (2)H-based dynamics when compared to wild-type PDZ2. In contrast, a change at a second core position (I20F) that had previously been shown to be part of an energetic and dynamic network, resulted in extensive changes in side-chain dynamics. This response is reminiscent to that seen previously upon peptide ligand binding. These results shed light on the nature of the PDZ2 dynamic network and suggest that position 20 in PDZ2 acts as a "hub" that is energetically and dynamically critical for transmitting changes in dynamics throughout the PDZ domain.     

The role of slow and fast protein motions in allosteric interactions

Allostery is fundamentally thermodynamic in nature. Long-range communication in proteins may be mediated not only by changes in the mean conformation with enthalpic contribution but also by changes in dynamic fluctuations with entropic contribution. The important role of protein motions in mediating allosteric interactions has been established by NMR spectroscopy. By using CAP as a model system, we have shown how changes in protein structure and internal dynamics can allosterically regulate protein function and activity. The results indicate that changes in conformational entropy can give rise to binding enhancement, binding inhibition, or have no effect in the expected affinity, depending on the magnitude and sign of enthalpy–entropy compensation. Moreover, allosteric interactions can be regulated by the modulation a low-populated conformation states that serve as on-pathway intermediates for ligand binding. Taken together, the interplay between fast internal motions, which are intimately related to conformational entropy, and slow internal motions, which are related to poorly populated conformational states, can regulate protein activity in a way that cannot be predicted on the basis of the protein’s ground-state structure.     

Ligand-release pathways in the pheromone-binding protein of Bombyx mori

Pheromone-binding proteins (PBP) supply olfactory neuron cells with pheromones by binding the ligands they are tailored for and carrying them to their receptor. The function of a PBP as an efficient carrier requires fast ligand uptake and release. The molecular basis of the ligand-binding mechanism was addressed here for the intriguing case of the PBP of the silk moth Bombyx mori. This PBP completely encapsulates its ligand bombykol without displaying any obvious ligand entrance/exit sites. Here, two opposite dissociation routes were identified as the most likely entrance/exit paths by replica-exchange molecular dynamics, essential dynamics, and force-probe molecular dynamics simulations. One of the paths runs along a flexible front lid; the other along the termini at the back. Calculated forces and energies suggest that both routes are physiologically relevant. The multiplicity of pathways may reduce or tune the entropic barrier for ligand binding.     

Dynamically driven protein allostery

Allosteric interactions are typically considered to proceed through a series of discrete changes in bonding interactions that alter the protein conformation. Here we show that allostery can be mediated exclusively by transmitted changes in protein motions. We have characterized the negatively cooperative binding of cAMP to the dimeric catabolite activator protein (CAP) at discrete conformational states. Binding of the first cAMP to one subunit of a CAP dimer has no effect on the conformation of the other subunit. The dynamics of the system, however, are modulated in a distinct way by the sequential ligand binding process, with the first cAMP partially enhancing and the second cAMP completely quenching protein motions. As a result, the second cAMP binding incurs a pronounced conformational entropic penalty that is entirely responsible for the observed cooperativity. The results provide strong support for the existence of purely dynamics-driven allostery.     

Mechanism of auxin interaction with Auxin Binding Protein (ABP1): a molecular dynamics simulation study

Auxin Binding Protein 1 (ABP1) is ubiquitous in green plants. It binds the phytohormone auxin with high specificity and affinity, but its role in auxin-induced processes is unknown. To understand the proposed receptor function of ABP1 we carried out a detailed molecular modeling study. Molecular dynamics simulations showed that ABP1 can adopt two conformations differing primarily in the position of the C-terminus and that one of them is stabilized by auxin binding. This is in agreement with experimental evidence that auxin induces changes at the ABP1 C-terminus. Simulations of ligand egress from ABP1 revealed three main routes by which an auxin molecule can enter or leave the ABP1 binding site. Assuming the previously proposed orientation of ABP1 to plant cell membranes, one of the routes leads to the membrane and the other two to ABP1's aqueous surroundings. A network of hydrogen-bonded water molecules leading from the bulk water to the zinc-coordinated ligands in the ABP1 binding site was formed in all simulations. Water entrance into the zinc coordination sphere occurred simultaneously with auxin egress. These results suggest that the hydrogen-bonded water molecules may assist in protonation and deprotonation of auxin molecules and their egress from the ABP1 binding site.     

Probing protein dynamics by nuclear magnetic resonance

Proteins are dynamic molecules that often undergo conformational changes while performing their specific functions, such as target recognition, ligand binding and catalysis. NMR spectroscopy is uniquely suited to study protein dynamics, because site-specific information can be obtained for motions that span a broad range of time scales. The information obtained from NMR dynamics experiments has provided insights into specific structural changes or conformational energetics associated with molecular function. In the last decade, a number of new advancements in NMR methodologies have further extended our ability to characterize protein dynamics. Here, we present an overview of current NMR technology that is used to monitor the dynamic properties of proteins.     

Long-Timescale Molecular-Dynamics Simulations of the Major Urinary Protein Provide Atomistic Interpretations of the Unusual Thermodynamics of Ligand Binding

The mouse major urinary protein (MUP) has proved to be an intriguing test bed for detailed studies on protein-ligand recognition. NMR, calorimetric, and modeling investigations have revealed that the thermodynamics of ligand binding involve a complex interplay between competing enthalpic and entropic terms. We performed six independent, 1.2 μs molecular-dynamics simulations on MUP—three replicates on the apo-protein, and three on the complex with the pheromone isobutylmethoxypyrazine. Our findings provide the most comprehensive picture to date of the structure and dynamics of MUP, and how they are modulated by ligand binding. The mechanical pathways by which amino acid side chains can transmit information regarding ligand binding to surface loops and either increase or decrease their flexibility (entropy-entropy compensation) are identified. Dewetting of the highly hydrophobic binding cavity is confirmed, and the results reveal an aspect of ligand binding that was not observed in earlier, shorter simulations: bound ligand retains extensive rotational freedom. Both of these features have significant implications for interpretations of the entropic component of binding. More generally, these simulations test the ability of current molecular simulation methods to produce a reliable and reproducible picture of protein dynamics on the microsecond timescale.     

Microcalorimetric indications for ligand binding as a function of the protein for galactoside-specific plant and avian lectins

The process cascade leading to the final accommodation of the carbohydrate ligand in the lectin's binding site comprises enthalpic and entropic contributions of the binding partners and solvent molecules. With emphasis on lactose, N-acetyllactosamine, and thiodigalactoside as potent inhibitors of binding of galactoside-specific lectins, the question was addressed to what extent these parameters are affected as a function of the protein. The microcalorimetric study of carbohydrate association to the galectin from chicken liver (CG-16) and the agglutinin from Viscum album (VAA) revealed enthalpy-entropy compensation with evident protein type-dependent changes for N-acetyllactosamine. Reduction of the entropic penalty by differential flexibility of loops or side chains and/or solvation properties of the protein will have to be reckoned with to assign a molecular cause to protein type-dependent changes in thermodynamic parameters for lectins sharing the same monosaccharide specificity.     

Conformational and dynamics changes induced by bile acids binding to chicken liver bile acid binding protein

The correlation between protein motions and function is a central problem in protein science. Several studies have demonstrated that ligand binding and protein dynamics are strongly correlated in intracellular lipid binding proteins (iLBPs), in which the high degree of flexibility, principally occurring at the level of helix-II, CD, and EF loops (the so-called portal area), is significantly reduced upon ligand binding. We have recently investigated by NMR the dynamic properties of a member of the iLBP family, chicken liver bile acid binding protein (cL-BABP), in its apo and holo form, as a complex with two bile salts molecules. Binding was found to be regulated by a dynamic process and a conformational rearrangement was associated with this event. We report here the results of molecular dynamics (MD) simulations performed on apo and holo cL-BABP with the aim of further characterizing the protein regions involved in motion propagation and of evaluating the main molecular interactions stabilizing bound ligands. Upon binding, the root mean square fluctuation values substantially decrease for CD and EF loops while increase for the helix-loop-helix region, thus indicating that the portal area is the region mostly affected by complex formation. These results nicely correlate with backbone dynamics data derived from NMR experiments. Essential dynamics analysis of the MD trajectories indicates that the major concerted motions involve the three contiguous structural elements of the portal area, which however are dynamically coupled in different ways whether in the presence or in the absence of the ligands. Motions of the EF loop and of the helical region are part of the essential space of both apo and holo-BABP and sample a much wider conformational space in the apo form. Together with NMR results, these data support the view that, in the apo protein, the flexible EF loop visits many conformational states including those typical of the holo state and that the ligand acts stabilizing one of these pre-existing conformations. The present results, in agreement with data reported for other iLBPs, sharpen our knowledge on the binding mechanism for this protein family.     

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.     

Accounting for global protein deformability during protein-protein and protein-ligand docking

Computational docking methods are valuable tools aimed to simplify the costly process of drug development and improvement. Most current approaches assume a rigid receptor structure to allow virtual screening of large numbers of possible ligands and putative binding sites on a receptor molecule. However, inclusion of receptor flexibility can be of critical importance since binding of a ligand can lead to changes in the receptor protein conformation that are sterically necessary to accommodate a ligand. Recent approaches to efficiently account for receptor flexibility during docking simulations are reviewed. In particular, accounting efficiently for global conformational changes of the protein backbone during docking is a still challenging unsolved problem. An approximate method has recently been suggested that is based on relaxing the receptor conformation during docking in pre-calculated soft collective degrees of freedom (M. Zacharias, Rapid protein-ligand docking using soft modes from molecular dynamics simulations to account for protein deformability: binding of FK506 to FKBP, Proteins: Struct., Funct., Genet. 54 (2004) 759-767). Test applications on protein-protein docking and on docking the inhibitor staurosporine to the apo-form of cAMP-dependent protein kinase A catalytic domain indicate significant improvement of docking results compared to rigid docking at a very modest computational demand. Accounting for receptor conformational changes in pre-calculated global degrees of freedom might offer a promising route to improve systematic docking screening simulations.     

Propagation of dynamic changes in barnase upon binding of barstar: an NMR and computational study

NMR spectroscopy and computer simulations were used to examine changes in chemical shifts and in dynamics of the ribonuclease barnase that result upon binding to its natural inhibitor barstar. Although the spatial structures of free and bound barnase are very similar, binding results in changes of the dynamics of both fast side-chains, as revealed by (2)H relaxation measurements, and NMR chemical shifts in an extended beta-sheet that is located far from the binding interface. Both side-chain dynamics and chemical shifts are sensitive to variations in the ensemble populations of the inter-converting molecular states, which can escape direct structural observation. Molecular dynamics simulations of free barnase and barnase in complex with barstar, as well as a normal mode analysis of barnase using a Gaussian network model, reveal relatively rigid domains that are separated by the extended beta-sheet mentioned above. The observed changes in NMR parameters upon ligation can thus be rationalized in terms of changes in inter-domain dynamics and in populations of exchanging states, without measurable structural changes. This provides an alternative model for the propagation of a molecular response to ligand binding across a protein that is based exclusively on changes in dynamics.     

Conformational dynamics of a multidomain protein by neutron scattering and computational analysis

The flexible conformations of a multidomain protein are responsible for its biological functions. Although MurD, a 47-kDa protein that consists of three domains, sequentially changes its domain conformation from an open form to a closed form through a semiclosed form in its enzymatic reaction, the domain dynamics in each conformation remains unclear. In this study, we verify the conformational dynamics of MurD in the corresponding three states (apo and ATP- and inhibitor-bound states) with a combination of small-angle x-ray and neutron scattering (SAXS and SANS), dynamic light scattering (DLS), neutron backscattering (NBS), neutron spin echo (NSE) spectroscopy, and molecular dynamics (MD) simulations. Applying principal component analysis of the MD trajectories, twisting and open-closed domain modes are identified as the major collective coordinates. The deviations of the experimental SAXS profiles from the theoretical calculations based on the known crystal structures become smaller in the ATP-bound state than in the apo state, and a further decrease is evident upon inhibitor binding. These results suggest that domain motions of the protein are suppressed step by step of each ligand binding. The DLS and NBS data yield collective and self-translational diffusion constants, respectively, and we used them to extract collective domain motions in nanometer and nanosecond scales from the NSE data. In the apo state, MurD shows both twisting and open-closed domain modes, whereas an ATP binding suppresses twisting domain motions, and a further reduction of open-closed mode is seen in the inhibitor-binding state. These observations are consistent with the structure modifications measured by the small-angle scattering as well as the MD simulations. Such changes in the domain dynamics associated with the sequential enzymatic reactions should be related to the affinity and reaction efficiency with a ligand that binds specifically to each reaction state.     

Differences in steroid specificity for rat androgen binding protein and the cytoplasmic receptor.

Two proteins in the rat, androgen binding protein (ABP) and the cytoplasmic receptor (CR), have high affinity and limited capacity for binding androgens. To determine the structural requirements for binding with high affinity, each protein was partially purified and the ability of over 100 steroids to compete with [3H]dihydrotestosterone (17 beta-hydroxy-5 alpha-androstan-3-one) for binding sites was assessed. The results indicate marked differences in the steroid specificities of the two proteins. Some alterations of dihydrotestosterone at C-2 or C-2 and C-3 increase binding to ABP two to four-fold. Similarly, the affinity of 17 beta-hydroxy-7 alpha-methyl-4-estren-3-one for ABP increases two-fold when a double bond is created at C-14. Addition of a methyl group in the alpha position at C-7 or C-17, or an ethinyl group at C-17 cause little change in affinity; however, modifications at C-11 and C-17 beta, and deletion of the methyl group at C-10 significantly impair binding to ABP. Binding to the CR is maintained or increased by deletion of the methyl group at C-10. Binding is lessened by modifications at C-3 and C-17 beta. Most alterations at C-2, C-7, C-11, and C-17 alpha have only minor effects on binding to the CR. These studies should provide a molecular basis for predicting the effects of specific structural modifications. When some modifications at C-2 or C-2 and C-3 are combined with changes at C-17 beta, the resulting steroids retain very high affinity for ABP and very limited binding to the CR. Such steroids may provide a means for assessing the function of ABP.     

Transporter-to-trap conversion: a disulfide bond formation in cellular retinoic acid binding protein I mutant triggered by retinoic acid binding irreversibly locks the ligand inside the protein

Transport proteins must bind their ligands reversibly to enable release at the point of delivery, while irreversible binding is usually associated with the extreme cases of ligand sequestration. Protein conformational dynamics is an important modulator of binding kinetics, as increased flexibility in the regions adjacent to the binding site may facilitate both association and dissociation processes. Ligand entry to, and exit from, the internal binding site of the cellular retinoic acid binding protein I (CRABP I) occurs via a flexible portal region, which functions as a dynamic aperture. We designed and expressed a CRABP I mutant (A35C/T57C), in which a small-scale conformational switch caused by the ligand binding event triggers formation of a disulfide bond in the portal region, thereby arresting structural fluctuations and effectively locking the ligand inside the binding cavity. At the same time, no formation of the disulfide bond is observed in the apo form of the mutant, and most characteristics of the mutant, including protein stability, are very similar to those of the wild-type protein in the absence of retinoic acid. The mutation does not alter the kinetics of retinoic acid binding to the protein, although the disulfide formation makes the binding effectively irreversible, as suggested by the absence of retinoic acid transfer from the holo form of the mutant to lipid vesicles in the absence of a reducing agent. Taken together, these data suggest that the disulfide bond formation in the portal region arrests large-scale structural fluctuations, which are required for retinoic acid release from the protein. The unique properties of the CRABP I mutant described in this work can be used to inspire and guide a design of nanodevices for multiple tasks ranging from sequestering small-molecule toxins in both tissue and circulation to nutrient deprivation of pathogens.     

Differences between apo and three holo forms of the intestinal fatty acid binding protein seen by molecular dynamics computer calculations

It is commonly believed that binding affinity can be estimated by consideration of local changes of ligand and protein. This paper discusses a set of molecular dynamics simulations of intestinal fatty acid binding protein addressing the protein's response to presence or absence of different ligands. A 5-ns simulation was performed of the protein without a ligand, and three simulations (one 5-ns and two 2-ns) were performed with different fatty acids bound. The results indicate that, although the basic protein structure is unchanged by the presence of the ligand, other properties are significantly affected by ligand binding. For example, zero-time covariance patterns between protein, bound waters, and ligand vary between the different simulations. Moreover, the interaction energies between ligand and specific residues indicate that different ligands are stabilized in different ways. In sum, the results suggest that binding thermodynamics within this system will need to be calculated not from a subset of nearby protein:ligand interactions, but will depend on a knowledge of the motions coupling together water, protein, and ligand.     

Epididymal androgen binding protein in protein-calorie malnourished rats

Androgen binding protein (ABP) was measured at five-day intervals from 20 to 50 days of age in the epididymides of rats fed a low protein diet. The results were compared to findings made in age-matched controls fed a standard laboratory diet. ABP was consistently lower in the caputs of the malnourished rats throughout the entire study. Moreover, no significant changes were noted in caput ABP concentrations as these rats aged. In contrast, this protein was observed to increase during the course of the month-long experiment in those animals fed the standard laboratory diet. These findings indicate that dietary protein deficiency has a marked debilitative effect on Sertoli cell function.     

Role of the Adenine Ligand on the Stabilization of the Secondary and Tertiary Interactions in the Adenine Riboswitch

Riboswitches are RNA-based genetic control elements that function via a conformational transition mechanism when a specific target molecule binds to its binding pocket. To facilitate an atomic detail interpretation of experimental investigations on the role of the adenine ligand on the conformational properties and kinetics of folding of the add adenine riboswitch, we performed molecular dynamics simulations in both the presence and the absence of the ligand. In the absence of ligand, structural deviations were observed in the J23 junction and the P1 stem. Destabilization of the P1 stem in the absence of ligand involves the loss of direct stabilizing interactions with the ligand, with additional contributions from the J23 junction region. The J23 junction of the riboswitch is found to be more flexible, and the tertiary contacts among the junction regions are altered in the absence of the adenine ligand; results suggest that the adenine ligand associates and dissociates from the riboswitch in the vicinity of J23. Good agreement was obtained with the experimental data with the results indicating dynamic behavior of the adenine ligand on the nanosecond time scale to be associated with the dynamic behavior of hydrogen bonding with the riboswitch. Results also predict that direct interactions of the adenine ligand with U74 of the riboswitch are not essential for stable binding although it is crucial for its recognition. The possibility of methodological artifacts and force-field inaccuracies impacting the present observations was checked by additional molecular dynamics simulations in the presence of 2,6-diaminopurine and in the crystal environment.