Regulation of tryptophan synthase by temperature, monovalent cations, and an allosteric ligand. Evidence from Arrhenius plots, absorption spectra, and primary kinetic isotope effects

To investigate the linkage between enzyme conformation and catalysis, we have determined the effects of temperature on catalytic properties of the tryptophan synthase alpha(2)beta(2) complex and beta(2) subunit in the absence or presence of different monovalent cations (Cs(+), Na(+), and GuH(+)) and of an allosteric ligand, alpha-glycerol 3-phosphate. Arrhenius plots of the activity data between 5 and 50 degrees C are nonlinear in the presence of certain ligands but not others. The conditions that yield nonlinear Arrhenius plots also yield temperature-dependent changes in the equilibrium distribution of enzyme-substrate intermediates and in primary kinetic isotope effects. The results provide evidence that the nonlinear Arrhenius plots are caused by a temperature-dependent conformational change that precedes the rate-limiting step in catalysis. Thermodynamic analysis of the data associated with the conformational change shows that the activation energies are much higher at low temperatures than at high temperatures. We correlate the results with a model in which the enzyme is converted by increased temperature under certain conditions from a low-activity "open" conformation to a high-activity "closed" conformation. The allosteric ligand and different monovalent cations, including GuH(+), which also acts as a chaotropic agent, affect the equilibrium between the open and closed forms. The large positive entropy changes in the conformational conversion suggest that the closed conformation results from tightened hydrophobic interactions that exclude water from the active site of the beta subunit.     

Packing interface energetics in different crystal forms of the λ Cro dimer

Variation among crystal structures of the λ Cro dimer highlights conformational flexibility. The structures range from a wild type closed to a mutant fully open conformation, but it is unclear if each represents a stable solution state or if one may be the result of crystal packing. Here we use molecular dynamics (MD) simulation to investigate the energetics of crystal packing interfaces and the influence of site‐directed mutagenesis on them in order to examine the effect of crystal packing on wild type and mutant Cro dimer conformation. Replica exchange MD of mutant Cro in solution shows that the observed conformational differences between the wild type and mutant protein are not the direct consequence of mutation. Instead, simulation of Cro in different crystal environments reveals that mutation affects the stability of crystal forms. Molecular Mechanics Poisson‐Boltzmann Surface Area binding energy calculations reveal the detailed energetics of packing interfaces. Packing interfaces can have diverse properties in strength, energetic components, and some are stronger than the biological dimer interface. Further analysis shows that mutation can strengthen packing interfaces by as much as ～5 kcal/mol in either crystal environment. Thus, in the case of Cro, mutation provides an additional energetic contribution during crystal formation that may stabilize a fully open higher energy state. Moreover, the effect of mutation in the lattice can extend to packing interfaces not involving mutation sites. Our results provide insight into possible models for the effect of crystallization on Cro conformational dynamics and emphasize careful consideration of protein crystal structures. Proteins 2014; 82:1128–1141. © 2013 Wiley Periodicals, Inc.     

Conformational Changes of Three Periplasmic Receptors for Bacterial Chemotaxis and Transport: The Maltose-, Glucose/Galactose- and Ribose-binding Proteins.

Small-angle X-ray scattering experiments were carried out for the maltose-, glucose/galactose- and ribose-binding proteins of Gram negative bacteria. All were shown to be monomers that decrease in radius of gyration on ligand binding.The results obtained for the maltose-binding protein agree well with crystal structures of the closed, ligand-bound, and open, ligand-free protein, suggesting that these are indeed the primary forms in solution. The closed form is stabilized by protein – sugar interactions, while the open conformation is stabilized by close contacts between the two domains. Since it is the proper spacial relationship of the domains in the closed form that is most important for interaction with chemotaxis and transport partners, the stabilization of the open form would help keep ligand-free molecules from interfering in function.The scattering results also provide evidence that a large conformational change takes place in association with ligand binding to the glucose/galactose- and ribose-binding proteins, and that the two changes are similar. Modeling suggests that the open forms resemble those found in the related leucine and leucine/isoleucine/valine-binding proteins, but are different from those observed for the maltose-binding protein and the related purine repressor.     

Conformational equilibria and free energy profiles for the allosteric transition of the ribose-binding protein

The ribose-binding protein (RBP) is a sugar-binding bacterial periplasmic protein whose function is associated with a large allosteric conformational change from an open to a closed conformation upon binding to ribose. The crystal structures of RBP in open and closed conformations have been solved. It has been hypothesized that the open and closed conformations exist in a dynamic equilibrium in solution, and that sugar binding shifts the population from open conformations to closed conformations. Here, we study by computer simulations the thermodynamic changes that accompany this conformational change, and model the structural changes that accompany the allosteric transition, using umbrella sampling molecular dynamics and the weighted histogram analysis method. The open state is comprised of a diverse ensemble of conformations; the open ribose-free X-ray crystal conformations being representative of this ensemble. The unligated open form of RBP is stabilized by conformational entropy. The simulations predict detectable populations of closed ribose-free conformations in solution. Additional interdomain hydrogen bonds stabilize this state. The predicted shift in equilibrium from the open to the closed state on binding to ribose is in agreement with experiments. This is driven by the energetic stabilization of the closed conformation due to ribose-protein interactions. We also observe a significant population of a hitherto unobserved ribose-bound partially open state. We believe that this state is the one that has been suggested to play a role in the transfer of ribose to the membrane-bound permease complex.     

The influence of 1-alkanols and external pressure on the lateral pressure profiles of lipid bilayers

The suggestion by Robert Cantor, that drug-induced pressure changes in lipid bilayers can change the conformational equilibrium between open and closed states of membrane proteins and thereby cause anesthesia, attracted much attention lately. Here, we studied the effect of both large external pressure and of 1-alkanols of different chain lengths—some of them anesthetics, others not—on the lateral pressure profiles across dimyristoylphosphatidylcholine (DMPC) bilayers by molecular dynamics simulations. For a pure DMPC bilayer, high pressure both reduced and broadened the tension at the interface hydrophobic/hydrophilic and diminished the repulsion between the phospholipid headgroups. Whereas the effect of ethanol on the lateral pressure profile was similar to the effect of a large external pressure on a DMPC bilayer, long-chain 1-alkanols significantly amplified local maxima and minima in the lateral pressure profile. For most 1-alkanols, external pressure had moderate effects and did not reverse the changes 1-alkanols exerted on the pressure profile. Nevertheless, assuming the bent helix model as a simple geometric model for the transmembrane region of a membrane protein, protein conformational equilibria were shifted in opposite directions by addition of 1-alkanols and additional application of external pressure.     

Protein conformer selection by ligand binding observed with crystallography.

A large-scale movement between "closed" and "open" conformations of a protein loop was observed directly with protein crystallography by trapping individual conformers through binding of an exogenous ligand and characterization with solution kinetics. The buried indole ring of Trp191 in cytochrome c peroxidase (CCP) was displaced by exogenous ligands, causing a conformational change of loop Pro190-Asn195 and exposing Trp191 to the protein surface. Kinetic measurements are consistent with a two-step binding mechanism in which the rate-limiting step is a transition of the protein to the open state, which then binds the ligand. This large-scale conformational change of a functionally important region of CCP is independent of ligand and indicates that about 4% of the wild-type protein is in the open form in solution at any given time.     

Computational design of an integrin I domain stabilized in the open high affinity conformation

We have taken a computational approach to design mutations that stabilize a large protein domain of approximately 200 residues in two alternative conformations. Mutations in the hydrophobic core of the alphaMbeta2 integrin I domain were designed to stabilize the crystallographically defined open or closed conformers. When expressed on the cell surface as part of the intact heterodimeric receptor, binding of the designed open and closed I domains to the ligand iC3b, a form of the complement component C3, was either increased or decreased, respectively, compared to wild type. Moreover, when expressed in isolation from other integrin domains using an artificial transmembrane domain, designed open I domains were active in ligand binding, whereas designed closed and wild type I domains were inactive. Comparison to a human expert designed open mutant showed that the computationally designed mutants are far more active. Thus, computational design can be used to stabilize a molecule in a desired conformation, and conformational change in the I domain is physiologically relevant to regulation of ligand binding.     

Structural studies of an engineered zinc biosensor reveal an unanticipated mode of zinc binding

Protein engineering was used previously to convert maltose-binding protein (MBP) into a zinc biosensor. Zn(2+) binding by the engineered MBP was thought to require a large conformational change from "open" to "closed", similar to that observed when maltose is bound by the wild-type protein. We show that although this re-designed MBP molecule binds Zn(2+) with high affinity as previously reported, it does not adopt a closed conformation in solution as assessed by small-angle X-ray scattering. High-resolution crystallographic studies of the engineered Zn(2+)-binding MBP molecule demonstrate that Zn(2+) is coordinated by residues on the N-terminal lobe only, and therefore Zn(2+) binding does not require the protein to adopt a fully closed conformation. Additional crystallographic studies indicate that this unexpected Zn(2+) binding site can also coordinate Cu(2+) and Ni(2+) with only subtle changes in the overall conformation of the protein. This work illustrates that the energetic barrier to domain closure, which normally functions to maintain MBP in an open concentration in the absence of ligand, is not easily overcome by protein design. A comparison to the mechanism of maltose-induced domain rearrangement is discussed.     

Calcium-dependent open/closed conformations and interfacial energy maps of reconstituted hemichannels

Using an atomic force microscope, we have studied three-dimensional molecular topography and calcium-sensitive conformational changes of individual hemichannels. Full-length (non-truncated) Cx43 hemichannels (connexons), when reconstituted in lipid bilayer, appear as randomly distributed individual particles and clusters. They show a lack of preferential orientation of insertion into lipid membrane; in a single bilayer, connexons with protrusion of either the extracellular face or the large non-truncated cytoplasmic face are observed. Extracellular domains of these undocked hemichannels are structurally different from hemichannels in the docked gap junctional plaques examined after their exposure by force dissection or chemical dissection. Calcium induced a reversible change in the extracellular pore diameter. Hemichannels imaged in a physiological buffer with 1.8 mm Ca(+2) had the pore diameter of approximately 1.8 nm, consistent with the closed channel conformation. Reducing Ca(+2) concentration to approximately 1.4, 1, and 0 mm, which changes hemichannels from the closed to open conformation, increased the pore diameter to approximately 2.5 nm for approximately 27, 74, and 100% of hemichannels, respectively. Thus, open/close probability of the hemichannel appears to be [Ca(2+)]-dependent. Computational analysis of the atomic force microscopy phase mode imaging reveals a significantly higher interfacial energy for open hemichannels that results from the interactions between the atomic force microscope probe and the hydrophobic domains. Thus, hydrophobic extracellular domains of connexins regulate calcium-dependent conformational changes.     

Osmolyte‐induced conformational changes in the Hsp90 molecular chaperone

Osmolytes are small molecules that play a central role in cellular homeostasis and the stress response by maintaining protein thermodynamic stability at controlled levels. The underlying physical chemistry that describes how different osmolytes impact folding free energy is well understood, however little is known about their influence on other crucial aspects of protein behavior, such as native‐state conformational changes. Here we investigate this issue with the Hsp90 molecular chaperone, a large dimeric protein that populates a complex conformational equilibrium. Using small angle X‐ray scattering we observe dramatic osmolyte‐dependent structural changes within the native ensemble. The degree to which different osmolytes affect the Hsp90 conformation strongly correlates with thermodynamic metrics of their influence on stability. This observation suggests that the well‐established osmolyte principles that govern stability also apply to large‐scale conformational changes, a proposition that is corroborated by structure‐based fitting of the scattering data, surface area comparisons and m‐value analysis. This approach shows how osmolytes affect a highly cooperative open/closed structural transition between two conformations that differ by a domain‐domain interaction. Hsp90 adopts an additional ligand‐specific conformation in the presence of ATP and we find that osmolytes do not significantly affect this conformational change. Together, these results extend the scope of osmolytes by suggesting that they can maintain protein conformational heterogeneity at controlled levels using similar underlying principles that allow them to maintain protein stability; however the relative impact of osmolytes on different structural states can vary significantly.     

Quantifying allosteric effects in proteins

In allosteric regulation, protein activity is altered when ligand binding causes changes in the protein conformational distribution. Little is known about which aspects of protein design lead to effective allosteric regulation, however. To increase understanding of the relation between protein structure and allosteric effects, we have developed theoretical tools to quantify the influence of protein-ligand interactions on probability distributions of reaction rates and protein conformations. We define the rate divergence, Dk, and the allosteric potential, Dx, as the Kullback-Leibler divergence between either the reaction-rate distributions or protein conformational distributions with and without the ligand bound. We then define Dx as the change in the conformational distribution of the combined protein/ligand system, derive Dx in the harmonic approximation, and identify contributions from 3 separate terms: the first term, D[stackxomega], results from changes in the eigenvalue spectrum; the second term, D[stackxDeltax], results from changes in the mean conformation; and the third term, Dxv, corresponds to changes in the eigenvectors. Using normal modes analysis, we have calculated these terms for a natural interaction between lysozyme and the ligand tri-N-acetyl-D-glucosamine, and compared them with calculations for a large number of simulated random interactions. The comparison shows that interactions in the known binding-site are associated with large values of Dxv. The results motivate using allosteric potential calculations to predict functional binding sites on proteins, and suggest the possibility that, in Nature, effective ligand interactions occur at intrinsic control points at which binding induces a relatively large change in the protein conformational distribution.     

6-Aminohexanoic acid-plasminogen interactions studied by fluorescence

Circular polarization of luminescence spectra of human plasminogen and of its derivatives were measured in solutions of ligand-free proteins and with saturating amounts of 6-aminohexanoic acid. Spectroscopic changes induced by the ligand reveal similar perturbations of the binding sites in all the protein derivatives. It is concluded that the gross conformational change induced by 6-aminohexanoic acid binding to the native plasminogen involves changes of sterical relations of entire protein domains.     

Energetic and spatial parameters for gating of the bacterial large conductance mechanosensitive channel, MscL

MscL is multimeric protein that forms a large conductance mechanosensitive channel in the inner membrane of Escherichia coli. Since MscL is gated by tension transmitted through the lipid bilayer, we have been able to measure its gating parameters as a function of absolute tension. Using purified MscL reconstituted in liposomes, we recorded single channel currents and varied the pressure gradient (P) to vary the tension (T). The tension was calculated from P and the radius of curvature was obtained using video microscopy of the patch. The probability of being open (Po) has a steep sigmoidal dependence on T, with a midpoint (T1/2) of 11.8 dyn/cm. The maximal slope sensitivity of Po/Pc was 0.63 dyn/cm per e-fold. Assuming a Boltzmann distribution, the energy difference between the closed and fully open states in the unstressed membrane was DeltaE = 18.6 kBT. If the mechanosensitivity arises from tension acting on a change of in-plane area (DeltaA), the free energy, TDeltaA, would correspond to DeltaA = 6.5 nm2. MscL is not a binary channel, but has four conducting states and a closed state. Most transition rates are independent of tension, but the rate-limiting step to opening is the transition between the closed state and the lowest conductance substate. This transition thus involves the greatest DeltaA. When summed over all transitions, the in-plane area change from closed to fully open was 6 nm2, agreeing with the value obtained in the two-state analysis. Assuming a cylindrical channel, the dimensions of the (fully open) pore were comparable to DeltaA. Thus, the tension dependence of channel gating is primarily one of increasing the external channel area to accommodate the pore of the smallest conducting state. The higher conducting states appear to involve conformational changes internal to the channel that don't involve changes in area.     

Reaction of canine plasminogen with 6-aminohexanoate: a thermodynamic study combining fluorescence, circular dichroism, and isothermal titration calorimetry

The thermodynamics of the binding of 6-aminohexanoate (6-AH) to dog glu-plasminogen has been studied. Fluorescence titrations revealed four binding sites. Three yielded positive fluorescence changes on ligand binding; one yielded a negative fluorescence change. The fluorescence data gave no indication of cooperative interactions. Binding was studied using circular dichroism (CD). Near 295 nm there were small changes associated with binding ligand. These were magnified at 235 nm, a wavelength that is mainly associated with tryptophan bands. The dissociation constants obtained from the fluorescence were applied to the CD data and fit quite well. Below 220 nm, there were no significant differences between samples with or without 6-AH and, therefore, no substantial change in the secondary structure of the protein. Isothermal titration calorimetry was used in combination with the binding constants from fluorescence to study the enthalpy and entropy contributions to 6-AH binding. The enthalpies of association for the four sites are all negative. Their absolute values are small for the tight sites and large for the weakest. -TDeltaS is negative for the tight sites and positive for the weakest. The binding of 6-AH to plasminogen is entropically driven for the two tightest sites and enthalpically driven for the weakest site. The binding of 6-AH to lys-plasminogen has been studied and differs slightly from binding to glu-plasminogen. Most importantly, the binding of 6-AH for the weak site goes from enthalpy- to entropy-driven as is found with the other sites.     

Conformational dynamics of a lipid-interacting protein: MD simulations of saposin B

Saposin B is a water soluble alpha-helical protein which can bind to membranes and extract selected lipids, especially cerebroside sulfates. The X-ray structure of saposin B is homodimeric. There are two conformations of the dimer in the crystal-one with a closed central cavity (the AB dimer) and one (the CD dimer) with a more open cavity. We have conducted a series of short (5 ns) molecular dynamics simulations of saposin B, starting from both the AB and CD conformations and with/without bound lipid and/or water molecules within the central hydrophobic cavity. The more open (CD) dimer showed greater conformational drift than the AB dimer. The conformational drift was also somewhat higher in the absence of bound lipid. Two more extended (30 ns) simulations of AB and CD dimers were performed and analyzed in terms of changes in intersubunit packing within the dimers. The AB dimer remained largely unchanged in conformation over the duration of the extended simulation. In contrast, the CD dimer underwent a substantial conformational change corresponding to a 'scissor' motion of the two monomers so as to compress the central cavity to a more closed conformation than that seen in the AB dimer structure. A H-bond between the Q53 and Y54 side chains of the alpha3 helices of the two opposing monomers seems to hold the dimer in this 'scissor-closed' conformation. We suggest that a cycle of conformational changes, expanding and compressing the central cavity of the saposin B dimer, may play a key role in facilitating lipid extraction from bilayers.     

Effect of ionic strength and oxidation on the P-loop conformation of the protein tyrosine phosphatase-like phytase, PhyAsr

The protein tyrosine phosphatase (PTP)-like phytase, PhyAsr, from Selenomonas ruminantium is a novel member of the PTP superfamily, and the only described member that hydrolyzes myo-inositol-1,2,3,4,5,6-hexakisphosphate. In addition to the unique substrate specificity of PhyAsr, the phosphate-binding loop (P-loop) has been reported to undergo a conformational change from an open (inactive) to a closed (active) conformation upon ligand binding at low ionic strength. At high ionic strengths, the P-loop was observed in the closed, active conformation in both the presence and absence of ligand. To test whether the P-loop movement can be induced by changes in ionic strength, we examined the effect that ionic strength has on the catalytic efficiency of PhyAsr, and determined the structure of the enzyme at several ionic strengths. The catalytic efficiency of PhyAsr is highly sensitive to ionic strength, with a seven-fold increase in k(cat)/K(m) and a ninefold decrease in K(m) when the ionic strength is increased from 100 to 500 mm. Surprisingly, the P-loop is observed in the catalytically competent conformation at all ionic strengths, despite the absence of a ligand. Here we provide structural evidence that the ionic strength dependence of PhyAsr and the conformational change in the P-loop are not linked. Furthermore, we demonstrate that the previously reported P-loop conformational change is a result of irreversible oxidation of the active site thiolate. Finally, we rationalize the observed P-loop conformational changes observed in all oxidized PTP structures.     

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.