Temperature Dependence of Intrinsically Disordered Proteins in Simulations: What are We Missing?

The temperature dependence of the conformational properties in simulations of the intrinsically disordered model protein histatin 5 has been investigated using different combinations of force fields, water models, and atomistic and coarse-grained methods. The results have been compared to experimental data obtained from NMR, SAXS, and CD experiments to assess the accuracy and validity of the simulations. The results showed that neither simulations completely agreed with the experimental data, nor did they agree with each other. It was however possible to conclude that the observed conformational changes upon variations in temperature were not at all driven by electrostatic interactions. The final conclusion was that none of the simulations that were investigated in this study was able to accurately capture the temperature induced conformational changes of our model IDP.

A growth hormone receptor SNP promotes lung cancer by impairment of SOCS2-mediated degradation.

Both humans and mice lacking functional growth hormone (GH) receptors are known to be resistant to cancer. Further, autocrine GH has been reported to act as a cancer promoter. Here we present the first example of a variant of the GH receptor (GHR) associated with cancer promotion, in this case lung cancer. We show that the GHRP495T variant located in the receptor intracellular domain is able to prolong the GH signal in vitro using stably expressing mouse pro-B-cell and human lung cell lines. This is relevant because GH secretion is pulsatile, and extending the signal duration makes it resemble autocrine GH action. Signal duration for the activated GHR is primarily controlled by suppressor of cytokine signalling 2 (SOCS2), the substrate recognition component of the E3 protein ligase responsible for ubiquitinylation and degradation of the GHR. SOCS2 is induced by a GH pulse and we show that SOCS2 binding to the GHR is impaired by a threonine substitution at Pro 495. This results in decreased internalisation and degradation of the receptor evident in TIRF microscopy and by measurement of mature (surface) receptor expression. Mutational analysis showed that the residue at position 495 impairs SOCS2 binding only when a threonine is present, consistent with interference with the adjacent Thr494. The latter is key for SOCS2 binding, together with nearby Tyr487, which must be phosphorylated for SOCS2 binding. We also undertook nuclear magnetic resonance spectroscopy approach for structural comparison of the SOCS2 binding scaffold Ile455-Ser588, and concluded that this single substitution has altered the structure of the SOCS2 binding site. Importantly, we find that lung BEAS-2B cells expressing GHRP495T display increased expression of transcripts associated with tumour proliferation, epithelial-mesenchymal transition and metastases (TWIST1, SNAI2, EGFR, MYC and CCND1) at 2 h after a GH pulse. This is consistent with prolonged GH signalling acting to promote cancer progression in lung cancer.

Formation of hydrogen bonds precedes the rate-limiting formation of persistent structure in the folding of ACBP.

A burst phase in the early folding of the four-helix two-state folder protein acyl-coenzyme A binding protein (ACBP) has been detected using quenched-flow in combination with site-specific NMR-detected hydrogen exchange. Several of the burst phase structures coincide with a structure consisting of eight conserved hydrophobic residues at the interface between the two N and C-terminal helices. Previous mutation studies have shown that the formation of this structure is rate limiting for the final folding of ACBP. The burst phase structures observed in ACBP are different from the previously reported collapsed types of burst phase intermediates observed in the folding of other proteins.

Acyl-coenzyme A binding protein (ACBP).

Acyl-coenzyme A binding proteins are known from a large group of eukaryote species and to bind a long chain length acyl-CoA ester with very high affinity. Detailed biochemical mapping of ligand binding properties has been obtained as well as in-depth structural studies on the bovine apo-protein and of the complex with palmitoyl-CoA using NMR spectroscopy. In the four alpha-helix bundle structure, a set of 21 highly conserved residues present in more that 90% of all known sequences of acyl-coenzyme A binding proteins constitutes three separate mini-cores. These residues are predominantly located at the helix-helix interfaces. From studies of a large set of mutant proteins the role of the conserved residues has been related to structure, function, folding and stability.

The formation of a native-like structure containing eight conserved hydrophobic residues is rate limiting in two-state protein folding of ACBP.

The acyl-coenzyme A-binding proteins (ACBPs) contain 26 highly conserved sequence positions. The majority of these have been mutated in the bovine protein, and their influence on the rate of two-state folding and unfolding has been measured. The results identify eight sequence positions, out of 24 probed, that are critical for fast productive folding. The residues are all hydrophobic and located in the interface between the N- and C-terminal helices. The results suggest that one specific site dominated by conserved hydrophobic residues forms the structure of the productive rate-determining folding step and that a sequential framework model can describe the protein folding reaction.

Conserved residues and their role in the structure, function, and stability of acyl-coenzyme A binding protein.

In the family of acyl-coenzyme A binding proteins, a subset of 26 sequence sites are identical in all eukaryotes and conserved throughout evolution of the eukaryotic kingdoms. In the context of the bovine protein, the importance of these 26 sequence positions for structure, function, stability, and folding has been analyzed using single-site mutations. A total of 28 mutant proteins were analyzed which covered 17 conserved sequence positions and three nonconserved positions. As a first step, the influence of the mutations on the protein folding reaction has been probed, revealing a folding nucleus of eight hydrophobic residues formed between the N- and C-terminal helices [Kragelund, B. B., et al. (1999) Nat. Struct. Biol. (In press)]. To fully analyze the role of the conserved residues, the function and the stability have been measured for the same set of mutant proteins. Effects on function were measured by the extent of binding of the ligand dodecanoyl-CoA using isothermal titration calorimetry, and effects on protein stability were measured with chemical denaturation followed by intrinsic tryptophan and tyrosine fluorescence. The sequence sites that have been conserved for direct functional purposes have been identified. These are Phe5, Tyr28, Tyr31, Lys32, Lys54, and Tyr73. Binding site residues are mainly polar or charged residues, and together, four of these contribute approximately 8 kcal mol-1 of the total free energy of binding of 11 kcal mol-1. The sequence sites conserved for stability of the structure have likewise been identified and are Phe5, Ala9, Val12, Leu15, Leu25, Tyr28, Lys32, Gln33, Tyr73, Val77, and Leu80. Essentially, all of the conserved residues that maintain the stability are hydrophobic residues at the interface of the helices. Only one conserved polar residue, Gln33, is involved in stability. The results indicate that conservation of residues in homologous proteins may result from a summed optimization of an effective folding reaction, a stable native protein, and a fully active binding site. This is important in protein design strategies, where optimization of one of these parameters, typically function or stability, may influence any of the others markedly.

Mapping the lifetimes of local opening events in a native state protein.

The rate constants for the processes that lead to local opening and closing of the structures around hydrogen bonds in native proteins have been determined for most of the secondary structure hydrogen bonds in the four-helix protein acyl coenzyme A binding protein. In an analysis that combines these results with the energies of activation of the opening processes and the stability of the local structures, three groups of residues in the protein structure have been identified. In one group, the structures around the hydrogen bonds have frequent openings, every 600 to 1,500 s, and long lifetimes in the open state, around 1 s. In another group of local structures, the local opening is a very rare event that takes place only every 15 to 60 h. For these the lifetime in the open state is also around 1 s. The majority of local structures have lifetimes between 2,000 and 20,000 s and relatively short lifetimes of the open state in the range between 30 and 400 ms. Mapping of these groups of amides to the tertiary structure shows that the openings of the local structures are not cooperative at native conditions, and they rarely if ever lead to global unfolding. The results suggest a mechanism of hydrogen exchange by progressive local openings.

Hydrophobic core substitutions in calbindin D9k: effects on Ca2+ binding and dissociation.

Hydrophobic core residues have a marked influence on the Ca2+-binding properties of calbindin D9k, even though there are no direct contacts between these residues and the bound Ca2+ ions. Eleven different mutants with substitutions in the hydrophobic core were produced, and their equilibrium Ca2+-binding constants measured from Ca2+ titrations in the presence of chromophoric chelators. The Ca2+-dissociation rate constants were estimated from Ca2+ titrations followed by 1H NMR1 and were measured more accurately using stopped-flow fluorescence. The parameters were measured at four KCl concentrations to assess the salt dependence of the perturbations. The high similarity between the NMR spectra of mutants and wild-type calbindin D9k suggests that the structure is largely unperturbed by the substitutions. More detailed NMR investigations of the mutant in which Val61 is substituted by Ala showed that the mutation causes only very minimal perturbations in the immediate vicinity of residue 61. Substitutions of alanines or glycines for bulky residues in the center of the core were found to have significant effects on both Ca2+ affinity and dissociation rates. These substitutions caused a reduction in affinity and an increase in off-rate. Small effects, both increases and decreases, were observed for substitutions involving residues far from the Ca2+ sites and toward the outer part of the hydrophobic core. The mutant with the substitution Phe66 --> Trp behaved differently from all other mutants, and displayed a 25-fold increase in overall affinity of binding two Ca2+ ions and a 6-fold reduction in calcium dissociation rate. A strong correlation (R = 0.94) was found between the observed Ca2+-dissociation rates and affinities, as well as between the salt dependence of the off-rate and the distance to the nearest Ca2+-coordinating atom. There was also a strong correlation (R = 0.95) between the Ca2+ affinity and stability of the Ca2+ state and a correlation (R = 0. 69) between the Ca2+ affinity and stability of the apo state, as calculated from the results in the present and preceding paper in this issue [Julenius, K., Thulin, E., Linse, S., and Finn, B. E. (1998) Biochemistry 37, 8915-8925]. The change in salt dependencies of koff and cooperativity were most pronounced for residues completely buried in the core of the protein (solvent accessible surface area approximately 0). Altogether, the results suggest that the hydrophobic core residues promote Ca2+ binding both by contributing to the preformation of the Ca2+ sites in the apo state and by preferentially stabilizing the Ca2+-bound state.

Barley lipid-transfer protein complexed with palmitoyl CoA: the structure reveals a hydrophobic binding site that can expand to fit both large and small lipid-like ligands.

BACKGROUND: . Plant nonspecific lipid-transfer proteins (nsLTPs) bind a variety of very different lipids in vitro, including phospholipids, glycolipids, fatty acids and acyl coenzyme As. In this study we have determined the structure of a nsLTP complexed with palmitoyl coenzyme A (PCoA) in order to further our understanding of the structural mechanism of the broad specificity of these proteins and its relation to the function of nsLTPs in vivo. RESULTS: . 1H and 13C nuclear magnetic resonance spectroscopy (NMR) have been used to study the complex between a nsLTP isolated from barley seeds (bLTP) and the ligand PCoA. The resonances of 97% of the 1H atoms were assigned for the complexed bLTP and nearly all of the resonances were assigned in the bound PCoA ligand. The palmitoyl chain of the ligand was uniformly 13C-labelled allowing the two ends of the hydrocarbon chain to be assigned. The comparison of a subset of 20 calculated structures to an average structure showed root mean square deviations of 1.89 +/- 0.19 for all C, N, O, P and S atoms of the entire complex and of 0.57 +/- 0.09 for the peptide backbone atoms of the four alpha helices of the complexed bLTP. The four-helix topology of the uncomplexed bLTP is maintained in the complexed form of the protein. The bLTP only binds the hydrophobic parts of PCoA with the rest of the ligand remaining exposed to the solvent. The palmitoyl chain moiety of the ligand is placed in the interior of the protein and bent in a U-shape. This part of the ligand is completely buried within a hydrophobic pocket of the protein. CONCLUSIONS: . A comparison of the structures of bLTP in the free and bound forms suggests that bLTP can accommodate long olefinic ligands by expansion of the hydrophobic binding site. This expansion is achieved by a bend of one helix, HA, and by conformational changes in both the C terminus and helix HC. This mode of binding is different from that seen in the structure of maize nsLTP in complex with palmitic acid, where binding of the ligand is not associated with structural changes.

Thermodynamics of ligand binding to acyl-coenzyme A binding protein studied by titration calorimetry.

Ligand binding to recombinant bovine acyl-CoA binding protein (ACBP) was examined using isothermal microcalorimetry. Microcalorimetric measurements confirm that the binding affinity of acyl-CoA esters for ACBP is strongly dependent on the length of the acyl chain with a clear preference for acyl-CoA esters containing more than eight carbon atoms and that the 3'-phosphate of the ribose accounts for almost half of the binding energy. Binding of acyl-CoA esters, with increasing chain length, to ACBP was clearly enthalpically driven with a slightly unfavorable entropic contribution. Accessible surface areas derived from the measured enthalpies were compared to those calculated from sets of three-dimensional solution structures and showed reasonable correlation, confirming the enthalphically driven binding. Binding of dodecanoyl-CoA to ACBP was studied at various temperatures and was characterized by a weak temperature dependence on delta G zero and a strong enthalpy-entropy compensation. This was a direct consequence of a large heat capacity delta Cp caused by the presence of strong hydrophobic interactions. Furthermore, the binding of dodecanoyl-CoA was studied at various pH values and ionic strengths. The data presented here state that ACBP binds long-chain acyl-CoA esters with very high affinity and suggest that ACBP acts as a housekeeping protein with no pronounced built-in specificity.

Fast and one-step folding of closely and distantly related homologous proteins of a four-helix bundle family.

Bovine acyl-coenzyme A binding protein is a four-helix bundle protein belonging to a group of homologous eukaryote proteins that binds medium and long-chain acyl-coenzyme A esters with a very high affinity. The three-dimensional structure of both the free and the ligated protein together with the folding kinetics have been described in detail for the bovine protein and with four new sequences reported here, a total of 16 closely related sequences ranging from yeasts and plants to human are known. The kinetics of folding and unfolding in different concentrations of guanidine hydrochloride together with equilibrium unfolding have been measured for bovine, rat and yeast acyl-coenzyme A binding protein. The bovine and rat sequences are closely related whereas the yeast is more distantly related to these. In addition to the three natural variants, kinetics of a bovine mutant protein, Tyr31 --> Asn, have been studied. Both the folding and unfolding rates in water of the yeast protein are 15 times faster than those of bovine. The folding rates in water of the two mammalian forms, rat and bovine, are similar, though still significantly different. A faster unfolding rate both for rat and the bovine mutant protein results from a lower stability of the native states of these. These hydrophobic regions, mini cores, have been identified in the three-dimensional structure of the bovine protein and found to be formed primarily by residues that have been conserved throughout the entire eukaryote evolution from yeasts to both plants and mammals as seen in the sample of 16 sequences. The conserved residues are found to stabilize helix-helix interactions and serve specific functional purposes for ligand binding. The fast one-step folding mechanism of ACBP has been shown to be a feature that seems to be maintained throughout evolution despite numerous differences in sequence and even dramatic differences in folding kinetics and protein stability. The protein study raises the question to what extent does the conserved hydrophobic residues provide a scaffold for an efficient one-step folding mechanism.

Local perturbations by ligand binding of hydrogen deuterium exchange kinetics in a four-helix bundle protein, acyl coenzyme A binding protein (ACBP).

Amide hydrogen exchange kinetics of the individual amides in a four-helix bundle protein, acyl-coenzyme A binding protein, have been studied by nuclear magnetic resonance spectroscopy. The kinetics of amides with exchange rate constants in the range of 10(-25) to 10(-6.5) S-1 at pH 6.65 in free protein and the ligand-protein complex have been measured, and the effect of binding the ligand, palmitoyl-coenzyme A, on individual exchange rates has been analysed. Specific correlations between exchange kinetics and the structural properties of the individual amides known from the three-dimensional structure of acyl-coenzyme A binding protein have been examined. Furthermore, an analysis has been performed comparing the structural perturbations of the protein-ligand interactions known from the three dimensional structure of the complex of palmitoyl-coenzyme A and acyl-coenzyme A binding protein with the ligand-induced perturbations on the amides exchange kinetics. Finally, the ligand-induced perturbations on hydrogen exchange have been compared with those on 15N relaxation. The results suggest that hydrogen exchange kinetics in the individual sites of acyl-coenzyme A binding protein are primarily determined by local structural features; they show that ligand binding gives rise mainly to changes localized at the sites of interaction between protein and ligand; they imply that the perturbation of exchange kinetics caused by ligation can be either, as in one example a local stabilisation of the pre-exchange equilibrium induced by formation of a hydrogen bond, or as seen here in several examples a reduction of the dynamic processes that lead to the opening and closing processes of the pre-exchange equilibrium. The results seem not to indicate changes in the rate of the final chemical exchange step.

Folding of a four-helix bundle: studies of acyl-coenzyme A binding protein.

The refolding from denaturing conditions of a small four-helix bundle, the acyl-coenzyme A binding protein, has been investigated by utilizing an array of fast-reaction techniques. Stopped-flow tryptophan fluorescence for measuring the incorporation of aromatic residues into the protein core and far- and near-ultraviolet circular dichroism to measure the formation of secondary and tertiary structure, respectively, together with the formation of persistent structure measured by hydrogen exchange pulse labeling experiments analyzed by electrospray ionisation mass spectrometry all show that 90% of the acyl-coenzyme A binding protein molecules achieve their fully folded and active, native state with a time constant of less than 5 ms at 25 degrees C and of ca. 30 ms at 5 degrees C. The kinetic parameters measured by the different techniques are closely similar, indicating that the different elements of structure form effectively concomitantly. There is no evidence for a significant population of any partially structured intermediate states, and the kinetics are identical whether refolding occurs from an unfolded state generated either by low pH or by addition of guanidine hydrochloride. The kinetics of both refolding and unfolding are monophasic processes for practically 90% of the molecules, and can be described by a two-state model. The results add to our knowledge of the folding scheme of different structural motifs and are discussed in terms of current views of the mechanisms of protein folding.

Three-dimensional structure of the complex between acyl-coenzyme A binding protein and palmitoyl-coenzyme A.

Multidimensional 1H, 13C and 15N nuclear magnetic resonance spectroscopy has been used to study the complex between palmitoyl-coenzyme A and acyl-coenzyme A binding protein. The 1H and the 15N spectra of the holo-protein have been almost completely assigned and so has most of the 1H spectrum of the coenzyme A part of the protein-bound ligand. The palmitoyl part of the ligand has been uniformly labelled with 13C and the nuclear magnetic resonance signals of the carbon atoms and their protons have been assigned at the two ends of the hydrocarbon chain. A total of 1251 distance restraints from nuclear Overhauser effects and 131 dihedral angle restraints from three-bond coupling constants provided the basis for the structure calculation. A comparison of 20 structures calculated from these data to the average structure showed that they could be aligned with an atomic root-mean-square deviation of 1.3(+/- 0.2) A for all C, N, O, P and S atoms in protein and ligand. The apo-protein is a four-helix protein and this structure is maintained in the holo-protein. The four alpha-helices are Ac1 of residues 3 to 15, Ac2 from residue 20 to 36, Ac3 from 51 to 62, and Ac4 from 65 to 84. For the four alpha-helices of the peptide backbone of the holo-protein the root-mean-square deviation for the C, C alpha and N atoms was 0.42(+/- 0.08) A. The binding site for the palmitoyl-chain stretches between the N-terminal end of Ac3 where the carboxyl part binds, to the N-terminal of Ac3 where the omega-end of the palmitoyl part binds. The adenosine-3'-phosphate is bound near residues of each of the four helices in an arrangement where it can form salt bridges and/or hydrogen bonds to either backbone or side-chain atoms of Ala9, Tyr28, Lys32, Lys54 and Tyr73. The polar parts of the pantetheine and the pyrophosphate are structured in the bound ligand to form an interface with the solvent. Also the ligand forms a set of non-polar intramolecular interactions where the adenine, the pantetheine, and the palmitoyl-chain are associated, so overall the structure of the bound ligand seems to be organized to protect the lipophilic palmitoyl part from the polar solvent.