Intrinsic versus extrinsic stabilization of enzymes: the interaction of solutes and temperature on A4-lactate dehydrogenase orthologs from warm-adapted and cold-adapted marine fishes.

We examined the effects of temperature and stabilizing solutes on A4-lactate dehydrogenase (A4-LDH) from warm- and cold-adapted fishes, to determine how extrinsic stabilizers affect orthologs with different intrinsic stabilities. Conformational changes during substrate binding are rate-limiting for A4-LDH, thus stabilization due to intrinsic or extrinsic factors leads to decreased activity. A4-LDH from a warm-temperate goby (Gillichthys mirabilis), which has lower values for kcat and the Michaelis constant for pyruvate ( K m PYR), was intrinsically more stable than the orthologs of the cold-adapted Antarctic notothenioids Parachaenichthys charcoti and Chionodraco rastrospinosus, as shown by a higher apparent transition ('melting') temperature (Tm(APP)). We used four solutes, glycerol, sucrose, trimethylamine-N-oxide and poly(ethylene glycol) 8000, which stabilize proteins through different modes of preferential exclusion, to study temperature-solute interactions of the three orthologs. Changes in Tm(APP) were similar for all orthologs in each solute tested, but the catalytic rate of G. mirabilis A4-LDH was decreased most by solutes and increased most by temperature. In contrast, the K m PYR values of the Antarctic orthologs were more affected than that of the goby by both solutes and temperature. We conclude that (a) preferential exclusion of solutes functions within the native state of A4-LDH to favor conformational microstates with minimal surface area; (b) the varied effects of the different solutes on the kinetic properties are due to the interaction between this nonspecific stabilization and the differing intrinsic stabilities of the orthologs; (c) the catalytic rates of A4-LDH orthologs are equally affected by stabilizing solutes, if measurements are made at physiologically appropriate temperatures; and (d) global stability and localized flexibility of these A4-LDH orthologs may evolve independently.     

Evolutionary convergence in adaptation of proteins to temperature: A4-lactate dehydrogenases of Pacific damselfishes (Chromis spp.)

We have compared the kinetic properties (Michaelis-Menten constant [K(m)] and catalytic rate constant [k(cat)]) and amino acid sequences of orthologs of lactate dehydrogenase-A (A(4)-LDH) from congeners of Pacific damselfishes (genus Chromis) native to cold-temperate and tropical habitats to elucidate mechanisms of enzymatic adaptation to temperature. Specifically, we determined whether the sites of adaptive variation and the types of amino acids involved in substitutions at these sites were similar in the Chromis orthologs and other orthologs of warm-adapted and cold-adapted A(4)-LDH previously studied. We report striking evolutionary convergence in temperature adaptation of this protein and present further support for the hypothesis that enzyme adaptation to temperature involves subtle amino acid changes at a few sites that affect the mobility of the portions of the enzyme that are involved in rate-determining catalytic conformational changes. We tested the predicted effects of differences in sequence using site-directed mutagenesis. A single amino acid substitution in a key hinge region of the A(4)-LDH molecule is sufficient to change the kinetic characteristics of a temperate A(4)-LDH to that of a tropical ortholog. This substitution is at the same location that was identified in previous studies of adaptive variation in A(4)-LDH and was hypothesized to be important in adjusting K(m) and k(cat). Our results suggest that certain sites within an enzyme, notably those that establish the energy changes associated with rate-limiting movements of protein structure during catalysis, are "hot spots" of adaptation and that common types of amino acid substitutions occur at these sites to adapt structural "flexibility" and kinetic properties. Thus, despite the wide array of options that proteins have to adjust their structural stabilities in the face of thermal stress, the adaptive changes that couple "flexibility" to alterations of function may be limited in their diversity.     

Lactate dehydrogenase from the Antarctic eelpout, Lycodichthys dearborni

As part of our studies to examine the molecular basis of cold-adaptation, we have determined the kinetic properties, thermal stability and deduced amino acid sequence of the enzyme lactate dehydrogenase (LDH) from an Antarctic zoarcid fish, Lycodichthys dearborni. Unlike Antarctic notothenioid fish which are endemic to the Southern Ocean, zoarcid fish are cosmopolitan and have a substantially longer evolutionary history as a sub-order. The A₄-LDH isoform was isolated and purified from the white muscle of L. dearborni. The kinetic parameters K_m^PYR and k_cat were determined at temperatures from 0 to 25°C. K_m^PYR was substantially higher at low temperatures than those from Antarctic and temperate notothenioid fish, whereas k_cat at these temperatures was essentially the same as those of the other fish LDH in this study. The sequence of L. dearborni A₄-LDH was determined from cDNA derived from white muscle RNA and found to be similar to, but distinct from, the A₄-LDH sequences of Antarctic notothenioid fish. Molecular modelling based on the structure of the A₄-LDH from Pagothenia borchgrevinki suggested that three conservative amino acid changes within the core of the protein that are not directly part of the active site but which might nonetheless influence the active site, may be important in cold-adaptation in L. dearborni A₄-LDH, and that several other changes on the surface of the protein might also play a role in cold-adaptation.     

Solutes modify a conformational transition in a membrane transport protein

The bacterial outer-membrane vitamin B₁₂ transporter, BtuB, undergoes a dramatic order-to-disorder transition in its N-terminal energy-coupling motif (Ton box) upon substrate binding. Here, site-directed spin labeling (SDSL) is used to show that a range of solutes prevents this conformational change when ligand is bound to BtuB, resulting in a more ordered Ton box structure. For each solute examined, the data indicate that solutes effectively block this conformational transition through an osmotic mechanism. The molecular weight dependence of this solute effect has been examined for a series of polyethylene glycols, and a sharp molecular weight cutoff is observed. This cutoff indicates that solutes are preferentially excluded from a cavity within the protein as well as the protein surface. Furthermore, the sensitivity of the conformational change to solution osmolality is consistent with a structural model predicted by SDSL. When the Ton box is unfolded by detergents or mutations (rather than by ligand binding), solutes, such as polyethylene glycols and salts, also induce a more structured compacted conformation. These results suggest that conformational changes in this class of outer membrane transporters, which involve modest energy differences and changes in hydration, may be modulated by a range of solutes, including solutes typically used in protein crystallization.     

Escape of a Small Molecule from Inside T4 Lysozyme by Multiple Pathways

The T4 lysozyme L99A mutant is often used as a model system to study small-molecule binding to proteins, but pathways for ligand entry and exit from the buried binding site and the associated protein conformational changes have not been fully resolved. Here, molecular dynamics simulations were employed to model benzene exit from its binding cavity using the weighted ensemble (WE) approach to enhance sampling of low-probability unbinding trajectories. Independent WE simulations revealed four pathways for benzene exit, which correspond to transient tunnels spontaneously formed in previous simulations of apo T4 lysozyme. Thus, benzene unbinding occurs through multiple pathways partially created by intrinsic protein structural fluctuations. Motions of several α-helices and side chains were involved in ligand escape from metastable microstates. WE simulations also provided preliminary estimates of rate constants for each exit pathway. These results complement previous works and provide a semiquantitative characterization of pathway heterogeneity for binding of small molecules to proteins.     

Substrate adaptabilities of Thermotogae mannan binding proteins as a function of their evolutionary histories

The Thermotogae possess a large number of ATP-binding cassette (ABC) transporters, including two mannan binding proteins, ManD and CelE (previously called ManE). We show that a gene encoding an ancestor of these was acquired by the Thermotogae from the archaea followed by gene duplication. To address the functional evolution of these proteins as a consequence of their evolutionary histories, we measured the binding affinities of ManD and CelE orthologs from representative Thermotogae. Both proteins bind cellobiose, cellotriose, cellotetraose, β-1,4-mannotriose, and β-1,4-mannotetraose. The CelE orthologs additionally bind β-1,4-mannobiose, laminaribiose, laminaritriose and sophorose while the ManD orthologs additionally only weakly bind β-1,4-mannobiose. The CelE orthologs have higher unfolding temperatures than the ManD orthologs. An examination of codon sites under positive selection revealed that many of these encode residues located near or in the binding site, suggesting that the proteins experienced selective pressures in regions that might have changed their functions. The gene arrangement, phylogeny, binding properties, and putative regulatory networks suggest that the ancestral mannan binding protein was a CelE ortholog which gave rise to the ManD orthologs. This study provides a window on how one class of proteins adapted to new functions and temperatures to fit the physiologies of their new hosts.     

Exploring Early Stages of the Chemical Unfolding of Proteins at the Proteome Scale

After decades of using urea as denaturant, the kinetic role of this molecule in the unfolding process is still undefined: does urea actively induce protein unfolding or passively stabilize the unfolded state? By analyzing a set of 30 proteins (representative of all native folds) through extensive molecular dynamics simulations in denaturant (using a range of force-fields), we derived robust rules for urea unfolding that are valid at the proteome level. Irrespective of the protein fold, presence or absence of disulphide bridges, and secondary structure composition, urea concentrates in the first solvation shell of quasi-native proteins, but with a density lower than that of the fully unfolded state. The presence of urea does not alter the spontaneous vibration pattern of proteins. In fact, it reduces the magnitude of such vibrations, leading to a counterintuitive slow down of the atomic-motions that opposes unfolding. Urea stickiness and slow diffusion is, however, crucial for unfolding. Long residence urea molecules placed around the hydrophobic core are crucial to stabilize partially open structures generated by thermal fluctuations. Our simulations indicate that although urea does not favor the formation of partially open microstates, it is not a mere spectator of unfolding that simply displaces to the right of the folded←→unfolded equilibrium. On the contrary, urea actively favors unfolding: it selects and stabilizes partially unfolded microstates, slowly driving the protein conformational ensemble far from the native one and also from the conformations sampled during thermal unfolding.     

谷氨酸棒杆菌anti-σ因子CseE及其突变体与σ因子SigE的相互作用分析

[目的]谷氨酸棒杆菌是重要的氨基酸生产菌株,本研究针对SigE与ZAS家族蛋白CseE相互作用机制进行探索研究,重点分析CseE突变体影响与SigE结合能力的机制。[方法]本研究选择谷氨酸棒杆菌ATCC 13032来源的SigE和CseE蛋白为研究目标,利用遗传学方法获得过表达的重组谷氨酸棒杆菌,通过RT-qPCR研究SigE调控sigE和cseE的转录情况。同时,利用ITC和His pull-down实验验证ZAS家族的CseE蛋白与Zn²⁺及SigE的结合情况。之后对CseE蛋白进行功能域分析、多序列比对,研究功能域关键氨基酸位点对SigE结合能力的影响。其次对SigE和CseE蛋白进行分子对接和动力学模拟,分析关键氨基酸影响其结合的机制。[结果]谷氨酸棒杆菌SigE调控基因sigE和cseE的转录并且其活性受CseE蛋白控制。CseE蛋白为ZAS家族蛋白,具有Zn²⁺结合能力。CseE_His83A、CseE_cys87A和CseE_cys90A突变体不会影响与SigE的结合能力,而CseE_C87A-C90A和CseE_{His83A-C87A-C90A}突变体与SigE的结合能力略有下降。分子动力学模拟发现SigE-CseE_C87A-C90A和SigE-CseE_{His83A-C87A-C90A}之间的结合能量为-17.23 kcal/mol和-14.06 kcal/mol,分别比未突变体系结合能量降低22.8%及36.9%。[结论]谷氨酸棒杆菌SigE通过聚集RNA聚合酶来调控基因sigE和cseE的表达。CseE蛋白属于ZAS家族,具有Zn²⁺结合能力同时通过与SigE蛋白互作来抑制SigE活性。CseE_C87A-C90A及CseE_{His83A-C87A-C90A}突变体能影响与SigE结合的能力,减弱对SigE活性的控制。本研究产生的三维结构和确定的氨基酸关键位点为后续探索谷氨酸棒杆菌SigE和CseE响应环境压力机制提供了理论基础。     

Structural basis for the binding of compatible solutes by ProX from the hyperthermophilic archaeon Archaeoglobus fulgidus

Compatible solutes such as glycine betaine and proline betaine serve as protein stabilizers because of their preferential exclusion from protein surfaces. To use extracellular sources of this class of compounds as osmo-, cryo-, or thermoprotectants, Bacteria and Archaea have developed high affinity uptake systems of the ATP-binding cassette type. These transport systems require periplasmic- or extracellular-binding proteins that are able to bind the transported substance with high affinity. Therefore, binding proteins that bind compatible solutes have to avoid the exclusion of their ligands within the binding pocket. In the present study we addressed the question to how compatible solutes can be effectively bound by a protein at temperatures around 83 degrees C as this is done by the ligand-binding protein ProX from the hyperthermophilic archaeon Archaeoglobus fulgidus. We solved the structures of ProX without ligand and in complex with both of its natural ligands glycine betaine and proline betaine, as well as in complex with the artificial ligand trimethylammonium. Cation-pi interactions and non-classical hydrogen bonds between four tyrosine residues, a main chain carbonyl oxygen, and the ligand have been identified to be the key determinants in binding the quaternary amines of the three investigated ligands. The comparison of the ligand binding sites of ProX from A. fulgidus and the recently solved structure of ProX from Escherichia coli revealed a very similar solution for the problem of compatible solute binding, although both proteins share only a low degree of sequence identity. The residues involved in ligand binding are functionally equivalent but not conserved in the primary sequence.     

Decreases in activation energy and substrate affinity in cold-adapted A4-lactate dehydrogenase: evidence from the Antarctic notothenioid fish Chaenocephalus aceratus

Enzyme function is strongly affected by temperature, and orthologs from species adapted to different thermal environments often show temperature compensation in kinetic properties. Antarctic notothenioid fishes live in a habitat of constant, extreme cold (-1.86 +/- 2 degrees C), and orthologs of the enzyme A4-lactate dehydrogenase (A4-LDH) in these species have adapted to this environment through higher catalytic rates, lower Arrhenius activation energies (Ea), and increases in the apparent Michaelis constant for the substrate pyruvate (Km(PYR)). Here, site-directed mutagenesis was used to determine which amino acid substitutions found in A4-LDH of the notothenioid Chaenocephalus aceratus, with respect to orthologs from warm-adapted teleosts, are responsible for these adaptive changes in enzyme function. Km(PYR) was measured in eight single and two double mutants, and Ea was tested in five single and two double mutants in the temperature range 0 degrees C-20 degrees C. Of the four mutants that had an effect on these parameters, two increased Ea but did not affect Km(PYR) (Gly224Ser, Ala310Pro), and two increased both Ea and Km(PYR) (Glu233Met, Gln317Val). The double mutants Glu233Met/Ala310Pro and Glu233Met/Gln317Val increased Km(PYR) and Ea to levels not significantly different from the A4-LDH of a warm temperate fish (Gillichthys mirabilis, habitat temperature 10 degrees C-35 degrees C). The four single mutants are associated with two alpha-helices that move during the catalytic cycle; those that affect Ea but not Km(PYR) are further from the active site than those that affect both parameters. These results provide evidence that (1) cold adaptation in A4-LDH involves changes in mobility of catalytically important molecular structures; (2) these changes may alter activation energy alone or activation energy and substrate affinity together; and (3) the extent to which these parameters are affected may depend on the location of the substitutions within the mobile alpha-helices, perhaps due to differences in proximity to the active site.     

Binding thermodynamics of serotonin to rat-brain 5-HT1a, 5HT2a and 5-HT3 receptors

The thermodynamic parameters ΔG° , ΔH° and Δs° of the binding equilibrium of serotonin to 5-HT_1A, 5-HT_2A and 5-HT₃ rat-brain membrane receptors have been determined by means of affinity constant measurements at six temperatures in the range 0 –35 ° C and van't Hoff plots. At variance with 5-HT_1A and 5-HT₃, the binding at the 5-HT_2A receptors is strongly endothermic and entropy-driven. Comparison with the results obtained by other authors on 5-HT_2A receptors in rats and humans suggests that the observed differences can be explained by a single amino acid difference in the receptor sequence between these two species.     

Applying conformational selection theory to improve crossdocking efficiency in 3‐phosphoinositide dependent protein kinase‐1

The emerging picture of biomolecular recognition is that of conformational selection followed by induced‐fit. Conformational selection theory states that binding partners exist in various conformations in solution, with binding involving a “selection” between complementary conformers. In this study, we devise a docking protocol that mimics conformational selection in protein–ligand binding and demonstrate that it significantly enhances crossdocking accuracy over Glide's flexible docking protocol, which is widely used in the pharmaceutical industry. Our protocol uses a pregenerated conformational ensemble to simulate ligand flexibility. The ensemble was generated by thorough conformational sampling coupled with conformer minimization. The generated conformers were then rigidly docked in the active site of the protein along with a postdocking minimization step that allows limited induced fit effects to be modeled for the ligand. We illustrate the improved performance of our protocol through crossdocking of 31 ligands to cocomplexed proteins of the kinase 3‐phosphoinositide dependent protein kinase‐1 extracted from the crystal structures 1H1W (ATP bound), 1OKY (staurosporine bound) and 3QD0 (bound to a potent inhibitor). Consistent with conformational selection theory, the performance of our protocol was the best for crossdocking to the cognate protein bound to the natural ligand, ATP. Proteins 2014; 82:436–451. © 2013 Wiley Periodicals, Inc.     

Coupling Protein Side-Chain and Backbone Flexibility Improves the Re-design of Protein-Ligand Specificity

Interactions between small molecules and proteins play critical roles in regulating and facilitating diverse biological functions, yet our ability to accurately re-engineer the specificity of these interactions using computational approaches has been limited. One main difficulty, in addition to inaccuracies in energy functions, is the exquisite sensitivity of protein–ligand interactions to subtle conformational changes, coupled with the computational problem of sampling the large conformational search space of degrees of freedom of ligands, amino acid side chains, and the protein backbone. Here, we describe two benchmarks for evaluating the accuracy of computational approaches for re-engineering protein-ligand interactions: (i) prediction of enzyme specificity altering mutations and (ii) prediction of sequence tolerance in ligand binding sites. After finding that current state-of-the-art “fixed backbone” design methods perform poorly on these tests, we develop a new “coupled moves” design method in the program Rosetta that couples changes to protein sequence with alterations in both protein side-chain and protein backbone conformations, and allows for changes in ligand rigid-body and torsion degrees of freedom. We show significantly increased accuracy in both predicting ligand specificity altering mutations and binding site sequences. These methodological improvements should be useful for many applications of protein – ligand design. The approach also provides insights into the role of subtle conformational adjustments that enable functional changes not only in engineering applications but also in natural protein evolution.     

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.     

Function of muscle-type lactate dehydrogenase and citrate synthase of the Galápagos marine iguana, Amblyrhynchus cristatus, in relation to temperature

The Galápagos marine iguana, Amblyrhynchus cristatus, is unique among lizards in foraging subtidally, leading to activity across a broad range of ambient temperatures ( approximately 14-40 degrees C). To determine whether the marine iguana shows any biochemical changes consistent with maintaining enzyme function at both warm and cold body temperatures, we examined the function of the aerobic enzyme citrate synthase (CS) and the muscle isoform of the anaerobic enzyme lactate dehydrogenase (A(4)-LDH) in A. cristatus and a confamilial species, Iguana iguana, from 14 to 46 degrees C. We also deduced amino acid sequences from cDNA of each enzyme. In CS, despite two amino acid substitutions, we found no difference in the apparent Michaelis-Menten constant K(m) of oxaloacetate at any temperature, indicating that the substrate affinity of CS in A. cristatus has not adapted to changes in thermal environment. In A(4)-LDH, we used site-directed mutagenesis to show that the substitutions T9A and I283V (A. cristatus --> I. iguana) individually have no effect on kinetics, but together significantly decrease the K(m) of pyruvate and catalytic rate constant (k(cat)) of the A. cristatus ortholog. Thus, our data show that A. cristatus A(4)-LDH has not become cold adapted in response to this species' aquatic foraging behavior, and instead may be consistent with moderate warm adaptation with respect to the I. iguana ortholog.     

The effects of the deep-sea environment on transmembrane signaling

Membrane-associated processes may be particularly susceptible to perturbation by the high hydrostatic pressures and low temperatures of the deep ocean. Transmembrane signaling by guanyl nucleotide binding protein (G protein) coupled receptors (GPCRs) is affected at a number of steps: (1) agonist activation of the GPCR; (2) the interaction of the receptor with the heterotrimeric G protein; (3) the G protein GTPase cycle; and (4) the activation and function of the effector element, adenylyl cyclase. The effects of low temperature and high hydrostatic pressures on the A₁ adenosine receptor–inhibitory G protein (G_i)–adenylyl cyclase signaling complex were examined in teleost fishes from three families, Scorpaenidae, Macrouridae and Moridae. In a comparison of teleost fishes, rat and chicken, species with body temperatures from 1 to 40 °C, at atmospheric pressure, A₁ adenosine receptor agonist binding is conserved at the body temperature of the species. In the marine teleost fishes examined, increased pressure decreases agonist efficacy. There are differences among species in the effects of increased hydrostatic pressure on G protein interactions with receptors, GTP binding to G protein α subunits and the intrinsic GTPase activity of α subunits. Adenylyl cyclase activity and modulation are affected by increased pressure in all the species examined, except Antimora rostrata which was unaffected by pressure changes. At pressures approximating those which the species experience in situ adenylyl cyclase activity retains its sensitivity to modulators. To understand the physiological consequences of impaired cell signaling several prototypical human diseases are discussed.     

Communication over the Network of Binary Switches Regulates the Activation of A2A Adenosine Receptor

Dynamics and functions of G-protein coupled receptors (GPCRs) are accurately regulated by the type of ligands that bind to the orthosteric or allosteric binding sites. To glean the structural and dynamical origin of ligand-dependent modulation of GPCR activity, we performed total ~ 5 μsec molecular dynamics simulations of A_2A adenosine receptor (A_2AAR) in its apo, antagonist-bound, and agonist-bound forms in an explicit water and membrane environment, and examined the corresponding dynamics and correlation between the 10 key structural motifs that serve as the allosteric hotspots in intramolecular signaling network. We dubbed these 10 structural motifs “binary switches” as they display molecular interactions that switch between two distinct states. By projecting the receptor dynamics on these binary switches that yield 2¹⁰ microstates, we show that (i) the receptors in apo, antagonist-bound, and agonist-bound states explore vastly different conformational space; (ii) among the three receptor states the apo state explores the broadest range of microstates; (iii) in the presence of the agonist, the active conformation is maintained through coherent couplings among the binary switches; and (iv) to be most specific, our analysis shows that W246, located deep inside the binding cleft, can serve as both an agonist sensor and actuator of ensuing intramolecular signaling for the receptor activation. Finally, our analysis of multiple trajectories generated by inserting an agonist to the apo state underscores that the transition of the receptor from inactive to active form requires the disruption of ionic-lock in the DRY motif.