The core lipocalin, bovine beta-lactoglobulin

The lipocalin family became established shortly after the structural similarity was noted between plasma retinol binding protein and the bovine milk protein, beta-lactoglobulin. During the past 60 years, beta-lactoglobulin has been studied by essentially every biochemical technique available and so there is a huge literature upon its properties. Despite all of these studies, no specific biological function has been ascribed definitively to the protein, although several possibilities have been suggested. During the processing of milk on an industrial scale, the unpredictable nature of the process has been put down to the presence of beta-lactoglobulin and certainly the whey protein has been implicated in the initiation of aggregation that leads to the fouling of heat exchangers. This short review of the properties of the protein will concentrate mainly on studies carried out under essentially physiological conditions and will review briefly some of the possible functions for the protein that have been described.     

Characterization of beta-lactotensin,a bioactive peptide derived from bovine beta-lactoglobulin,as a neurotensin agonist

beta-Lactotensin (beta-LT: His-Ile-Arg-Leu) is an ileum-contracting peptide derived from residues No. 146-149 of bovine beta-lactoglobulin. The ileum-contracting activity of beta-LT was blocked by the NT1 antagonist SR48692. beta-LT was selective for the neurotensin NT2 receptor while neurotensin was selective for the NT1 receptor. beta-LT is the first natural ligand showing selectivity for the NT2 receptor. beta-LT showed hypertensive activity after intravenous administration at a dose of 30 mg/kg in conscious rats, while neurotensin showed hypotensive activity. The hypertensive activity of beta-LT was blocked by levocabastine (1 mg/kg, i.v.), an NT2 antagonist. SR48692, which blocked the hypotensive activity of neurotensin, had no effect on the hypertensive activity of beta-LT. These results suggest that the hypertensive activity of beta-LT is mediated by the NT2 receptor. It was concluded that the NT1 and NT2 receptors mediate the opposite effect on blood pressure.     

Isolation of lactoperoxidase, lactoferrin, alpha-lactalbumin, beta-lactoglobulin B and beta-lactoglobulin A from bovine rennet whey using ion exchange chromatography

A mild and rapid method is described for isolating various milk proteins from bovine rennet whey. beta-Lactoglobulin from bovine rennet whey was easily adsorbed on and desorbed from a weak anion exchanger, diethylaminoethyl-Toyopearl. However, alpha-lactalbumin could not be adsorbed onto the resin. alpha-Lactalbumin and beta-lactoglobulin from rennet whey could also be adsorbed and separated using a strong anion exchanger, quaternary aminoethyl-Toyopearl. The rennet whey was passed through a strong cation exchanger, sulphopropyl-Toyopearl, to separate lactoperoxidase and lactoferrin. alpha-Lactalbumin and beta-lactoglobulin were adsorbed onto quaternary aminoethyl-Toyopearl. alpha-Lactalbumin was eluted using a linear (0-0.15 M) concentration gradient of NaCl in 0.05 M Tris-HCl buffer (pH 8.5). Subsequently, beta-lactoglobulin B and beta-lactoglobulin A were eluted from the column with 0.05 M Tris-HCl (pH 6.8), using a linear (0.1-0.25 M) concentration gradient of NaCl. The yields were 1260 mg alpha-lactalbumin, 1290 mg beta-lactoglobulin B and 2280 mg beta-lactoglobulin A from 1 l rennet whey.     

beta-Lactotensin, a neurotensin agonist peptide derived from bovine beta-lactoglobulin, enhances memory consolidation in mice

beta-Lactotensin (His-Ile-Arg-Leu) is an ileum-contracting tetrapeptide isolated from bovine beta-lactoglobulin. We previously reported that a neurotensin agonist beta-lactotensin shows antinociceptive effect through neurotensin NT(2) receptor. We found that centrally or orally administered beta-lactotensin at a dose of 60nmol/mouse or 300-500mg/kg, respectively, increased memory consolidation in the step-through-type inhibitory avoidance test in mice. The memory-enhancing activity of beta-lactotensin was inhibited by the dopamine D(2) receptor antagonist raclopride but not the D(1) receptor antagonist SCH23390. Taken together, beta-lactotensin might improve memory consolidation through activating the dopamine D(2) receptor.     

Reversible unfolding of bovine beta-lactoglobulin mutants without a free thiol group

Bovine beta-lactoglobulin (beta-lg) has been used extensively as a model for studying protein folding. One of the problems preventing clarification of the folding mechanism is the incomplete reversibility from the unfolded state, probably caused by the thiol-disulfide exchange between a free thiol at Cys-121 and two disulfide bonds. We constructed and expressed three beta-lg subtype A mutants in which Cys-121 was replaced by Ala, Ser, or Val (i.e. C121A, C121S, and C121V). We studied the reversibilities of these mutants from urea denaturation using circular dichroism, tryptophan fluorescence, reversed-phase and gel-filtration high performance liquid chromatographies, and SDS-PAGE. The folded structure of each mutant was similar to that of wild-type beta-lg. Urea-induced unfolding at pH 7.0 and 3.0 showed that although the C121S mutation notably decreases the stability, the destabilizing effects of the C121A and C121V mutations are less severe. For all of the mutants, complete refolding from the unfolded state in 8 M urea at both pH 7.0 and 3.0 was observed. Kinetics of the formation of the irreversibly unfolded species of wild-type beta-lg in 8 M urea at pH 7.0 indicated that, first, an intramolecular thiol-disulfide exchange occurs to produce a mixture of species with non-native disulfide bonds followed by the intermolecular thiol-disulfide exchange producing the oligomers. These results indicate that intramolecular and intermolecular thiol-disulfide exchange reactions cause the low reversibility of wild-type beta-lg especially at neutral pH and that the mutation of Cys-121 improves the reversibility, enabling us to study the folding of beta-lg more exactly under various conditions.     

Measuring the folding landscape of a harmonically constrained biopolymer

Pioneering studies have shown that the probability distribution of opening length for a DNA hairpin, recorded under constant force using an optical trap, can be used to reconstruct the energy landscape of the transition. However, measurements made under constant force are subject to some limitations. Under constant force a system with a sufficiently high energy barrier spends most of its time in the closed or open conformation, with relatively few statistics collected in the transition state region. We describe a measurement scheme in which the system is driven progressively through the transition by an optical trap and an algorithm is used to extract the energy landscape of the transition from the fluctuations recorded during this process. We illustrate this technique in simulations and demonstrate its effectiveness in experiments on a DNA hairpin. We find that the combination of this technique with the use of short DNA handles facilitates a high-resolution measurement of the hairpin's folding landscape with a very short measurement time.     

Peptide models of folding initiation sites of bovine beta-lactoglobulin: identification of nativelike hydrophobic interactions involving G and H strands

In an attempt to characterize the early folding events in bovine beta-lactoglobulin (BLG), a set of peptides, covering the flexible N-terminal region and the stable C-terminus beta-core, was synthesized and analyzed by circular dichroism and by nuclear magnetic resonance in water, trifluoroethanol (TFE), and sodium dodecyl sulfate (SDS) below and above the critical micellar concentration. The role of local and long-range hydrophobic interactions in guiding the folding has been investigated. For the peptide fragment covering the more flexible N-terminal region of BLG (beta-strands A, B), where both theoretical predictions and kinetic refolding experiments suggested the formation of non-native alpha-helix, no native long-range contacts were identified, and a helical secondary structure was stabilized only in the presence of 25 mM SDS. At variance, in 50% (v/v) TFE, native, long-range hydrophobic interactions were observed in the peptide covering the core region comprising G and H beta-strands. The side chains involved in these interactions form a nativelike hydrophobic cluster, thus suggesting that the GH region may act as the folding initiation site for BLG. This result is reinforced by the identification, in the urea denaturated BLG, of residual structure located at the level of the GH interface, as evidenced by NMR analysis. These results, in excellent agreement with kinetic, thermodynamic, and cold denaturation folding data, once more underline the utmost importance of the GH region for the stability and folding of BLG. Severe aggregation effects prevented the structural analysis of the peptide covering the EFGH region, indicating that this larger segment does not represent an independent folding domain and that the terminal alpha-helix is necessary for stabilizing the BLG folding core.     

Morphogenesis of a protein: folding pathways and the energy landscape

Current knowledge on the reaction whereby a protein acquires its native three-dimensional structure was obtained by and large through characterization of the folding mechanism of simple systems. Given the multiplicity of amino acid sequences and unique folds, it is not so easy, however, to draw general rules by comparing folding pathways of different proteins. In fact, quantitative comparison may be jeopardized not only because of the vast repertoire of sequences but also in view of a multiplicity of structures of the native and denatured states. We have tackled the problem of the relationships between the sequence information and the folding pathway of a protein, using a combination of kinetics, protein engineering and computational methods, applied to relatively simple systems. Our strategy has been to investigate the folding mechanism determinants using two complementary approaches, i.e. (i) the study of members of the same family characterized by a common fold, but substantial differences in amino acid sequence, or (ii) heteromorphic pairs characterized by largely identical sequences but with different folds. We discuss some recent data on protein-folding mechanisms by presenting experiments on different members of the PDZ domain family and their circularly permuted variants. Characterization of the energetics and structures of intermediates and TSs (transition states), obtained by Φ-value analysis and restrained MD (molecular dynamics) simulations, provides a glimpse of the malleability of the dynamic states and of the role of the topology of the native states and of the denatured states in dictating folding and misfolding pathways.     

A conserved interdomain interaction is a determinant of folding cooperativity in the GST fold

The canonical glutathione transferase (GST) fold found in many monomeric and dimeric proteins consists of two domains that differ in structure and conformational dynamics. However, no evidence exists that the two domains unfold/fold independently at equilibrium, indicating the significance of interdomain interactions in governing cooperativity between domains. Bioinformatics analyses indicate the interdomain interface of the GST fold is large, predominantly hydrophobic with a high packing density explaining cooperative interdomain behavior. Structural alignments reveal a topologically conserved lock-and-key interaction across the domain interface in which a bulky hydrophobic residue ("key") protrudes from the surface of the N-domain and inserts into a pocket ("lock") in the C-domain. To better understand the molecular basis for the contribution of interdomain interactions toward cooperativity within the GST fold in the absence of any influence from quaternary interactions, studies were done with two monomeric GST proteins: Escherichia coli Grx2 (EcGrx2) and human CLIC1 (hCLIC1). Replacing the methionine "key" residue with alanine is structurally nondisruptive, whereas it significantly diminishes the folding cooperativity of both proteins. The loss in cooperativity between domains in the mutants is reflected by a change in the equilibrium folding mechanism from a wild-type two-state process to a three-state process, populating a stable folding intermediate.     

Alternate pathways for folding in the flavodoxin fold family revealed by a nucleation-growth model

A recent study of experimental results for flavodoxin-like folds suggests that proteins from this family may exhibit a similar, signature pattern of folding intermediates. We study the folding landscapes of three proteins from the flavodoxin family (CheY, apoflavodoxin, and cutinase) using a simple nucleation and growth model that accurately describes both experimental and simulation results for the transition state structure, and the structure of on-pathway and misfolded intermediates for CheY. Although the landscape features of these proteins agree in basic ways with the results of the study, the simulations exhibit a range of folding behaviours consistent with two alternate folding routes corresponding to nucleation and growth from either side of the central beta-strand.     

Equilibrium and kinetic folding pathways of a TIM barrel with a funneled energy landscape

The role of native contact topology in the folding of a TIM barrel model based on the alpha-subunit of tryptophan synthase (alphaTS) from Salmonella typhimurium (Protein Data Bank structure 1BKS) was studied using both equilibrium and kinetic simulations. Equilibrium simulations of alphaTS reveal the population of two intermediate ensembles, I1 and I2, during unfolding/refolding at the folding temperature, Tf = 335 K. Equilibrium intermediate I1 demonstrates discrete structure in regions alpha0-beta6 whereas intermediate I2 is a loose ensemble of states with N-terminal structure varying from at least beta1-beta3 (denoted I2A) to alpha0-beta4 at most (denoted I2B). The structures of I1 and I2 match well with the two intermediate states detected in equilibrium folding experiments of Escherichia coli alphaTS. Kinetic folding simulations of alphaTS reveal the sequential population of four intermediate ensembles, I120Q, I200Q, I300Q, and I360Q, during refolding. Kinetic intermediates I120Q, I200Q, and I300Q are highly similar to equilibrium alphaTS intermediates I2A, I2B, and I1, respectively, consistent with kinetic experiments on alphaTS from E. coli. A small population (approximately 10%) of kinetic trajectories are trapped in the I120Q intermediate ensemble and require a slow and complete unfolding step to properly refold. Both the on-pathway and off-pathway I120Q intermediates show structure in beta1-beta3, which is also strikingly consistent with kinetic folding experiments of alphaTS. In the off-pathway intermediate I(120Q), helix alpha2 is wrapped in a nonnative chiral arrangement around strand beta3, sterically preventing the subsequent folding step between beta3 and beta4. These results demonstrate the success of combining kinetic and equilibrium simulations of minimalist protein models to explore TIM barrel folding and the folding of other large proteins.     

Spatial confinement is a major determinant of the folding landscape of human chromosomes

The global architecture of the cell nucleus and the spatial organization of chromatin play important roles in gene expression and nuclear function. Single-cell imaging and chromosome conformation capture-based techniques provide a wealth of information on the spatial organization of chromosomes. However, a mechanistic model that can account for all observed scaling behaviors governing long-range chromatin interactions is missing. Here we describe a model called constrained self-avoiding chromatin (C-SAC) for studying spatial structures of chromosomes, as the available space is a key determinant of chromosome folding. We studied large ensembles of model chromatin chains with appropriate fiber diameter, persistence length and excluded volume under spatial confinement. We show that the equilibrium ensemble of randomly folded chromosomes in the confined nuclear volume gives rise to the experimentally observed higher-order architecture of human chromosomes, including average scaling properties of mean-square spatial distance, end-to-end distance, contact probability and their chromosome-to-chromosome variabilities. Our results indicate that the overall structure of a human chromosome is dictated by the spatial confinement of the nuclear space, which may undergo significant tissue- and developmental stage-specific size changes.     

Protein topology affects the appearance of intermediates during the folding of proteins with a flavodoxin-like fold

The topology of a native protein influences the rate with which it is formed, but does topology affect the appearance of folding intermediates and their specific role in kinetic folding as well? This question is addressed by comparing the folding data recently obtained on apoflavodoxin from Azotobacter vinelandii with those available on all three other alpha-beta parallel proteins the kinetic folding mechanism of which has been studied, i.e. Anabaena apoflavodoxin, Fusarium solani pisi cutinase and CheY. Two kinetic folding intermediates, one on-pathway and the other off-pathway, seem to be present during the folding of proteins with an alpha-beta parallel, also called flavodoxin-like, topology. The on-pathway intermediate lies on a direct route from the unfolded to the native state of the protein involved. The off-pathway intermediate needs to unfold to allow the production of native protein. Available simulation data of the folding of CheY show the involvement of two intermediates with characteristics that resemble those of the two intermediates experimentally observed. Apparently, protein topology governs the appearance and kinetic roles of protein folding intermediates during the folding of proteins that have a flavodoxin-like fold.     

Effects of a helix substitution on the folding mechanism of bovine alpha-lactalbumin

The structure, stability, and unfolding-refolding kinetics of a chimeric protein, in which the amino acid sequence of the flexible loop region (residues 105-110) comes from equine lysozyme and the remainder of the sequence comes from bovine alpha-lactalbumin were studied by circular dichroism spectroscopy and stopped-flow measurements, and the results were compared with those of bovine alpha-lactalbumin. The substitution of the flexible loop in bovine alpha-lactalbumin with the helix D of equine lysozyme destabilizes the molten globule state, although the native state is significantly stabilized by substitution of the flexible loop region. The kinetic refolding and unfolding experiments showed that the chimeric protein refolds significantly faster and unfolds substantially slower than bovine alpha-lactalbumin. To characterize the transition state between the molten globule and the native states, we investigated the guanidine hydrochloride concentration dependence of the rate constants of refolding and unfolding. Despite the significant differences in the stabilities of both the molten globule and native states between the chimeric protein and bovine alpha-lactalbumin, the free energy level of the transition state is not affected by the amino acid substitution in the flexible loop region. Our results suggest that the destabilization in the molten globule state of the chimeric protein is caused by the disruption of the non-native interaction in the flexible loop region and that the disruption of the non-native interaction reduces the free energy barrier of refolding. We conclude that the non-native interaction in the molten globule state may act as a kinetic trap for the folding of alpha-lactalbumin.     

Folding-unfolding equilibrium and kinetics of equine beta-lactoglobulin: equivalence between the equilibrium molten globule state and a burst-phase folding intermediate

The denaturant-induced equilibrium unfolding transition of equine beta-lactoglobulin was investigated by ultraviolet absorption, fluorescence, and circular dichroism (CD) spectra. An equilibrium intermediate populates at moderate denaturant concentrations, and its CD spectrum is similar to that of the molten globule state previously observed for this protein at acid pH [Ikeguchi, M., Kato, S., Shimizu, A., and Sugai, S. (1997) Proteins: Struct., Funct., Genet. 27, 567-575]. The unfolding and refolding kinetics were also investigated by the stopped-flow CD and fluorescence. A significant change in the CD intensity was observed within the dead time of measurements (25 ms) when the refolding reaction was initiated by diluting the urea-unfolded protein solution, indicating the transient accumulation of the folding intermediate. The CD spectrum of this burst-phase intermediate agrees well with that of the molten globule state at acid pH. The stability of the burst-phase intermediate was also estimated from the urea-concentration dependence of the burst-phase amplitude, and it shows a fair agreement with that of the equilibrium intermediate. These results indicate that the molten globule state of equine beta-lactoglobulin populates at moderate urea concentration as well as at acid pH and it is equivalent with the kinetic folding intermediate.     

The kinetics of bovine growth hormone folding are consistent with a framework model

The framework model of protein folding requires the hydrogen-bonded secondary structure to be formed early in folding (i.e. the formation of secondary structure precedes the tertiary structure) (Kim, P. S., and Baldwin, R. L. (1982) Annu. Rev. Biochem. 51, 459-489). To test the framework model directly the kinetics of bovine growth hormone (bGH) folding were compared utilizing two methods of detection, one that measures the secondary structure (far ultraviolet circular dichroism) and another that measures the tertiary structure (near ultraviolet absorbance). The results demonstrate that, under identical folding conditions, the kinetics observed by far ultraviolet circular dichroism are faster than those observed by ultraviolet absorption. The faster kinetics observed by circular dichroism indicate the existence of a helix-containing intermediate which is consistent with the framework model. The effect of protein concentration and denaturant concentration on the kinetics of refolding were studied. The rate of refolding measured by absorbance and circular dichroism was dependent on protein concentration. The protein concentration dependence on refolding is due to the transient formation of an associated intermediate. The concentration dependence of folding is taken as evidence that folding is a sequential process with partially folded monomers responsible for the observed association effect. At dilute protein concentrations the refolding can be studied independent of the association phenomena. The growth hormones utilized in this study were derived from Escherichia coli through recombinant DNA technology and from bovine pituitaries. The pituitary-derived bGH has been shown to be heterogeneous at the NH2 terminus (Lorenson, M. F., and Ellis, S. (1975) Endocrinology 96, 833-838), whereas the recombinant DNA-derived bGH contains a single NH2 terminus. No differences in the folding kinetics between the recombinant DNA and pituitary derived-bGH were observed. It is concluded that the heterogeneity of the NH2 terminus of growth hormone obtained from bovine pituitaries does not affect the observed in vitro folding kinetics.