The molten globule state of a chimera of human alpha-lactalbumin and equine lysozyme.

The molten globule state of equine lysozyme is more stable than that of alpha-lactalbumin and is stabilized by non-specific hydrophobic interactions and native-like hydrophobic interactions. We constructed a chimeric protein which is produced by replacing the flexible loop (residues 105-110) in human alpha-lactalbumin with the helix D (residues 109-114) in equine lysozyme to investigate the possible role of the helix D for the high stability and native-like packing interaction in the molten globule state of equine lysozyme. The stability of the molten globule state formed by the chimeric protein to guanidine hydrochloride-induced unfolding is the same as that of equine lysozyme and is substantially greater than that of human alpha-lactalbumin, although only six residues come from equine lysozyme. Our results also suggest that the non-native interaction in the molten globule state of alpha-lactalbumin changes to the native-like packing interaction due to helix substitution. The solvent-accessibility of the Trp residues in the molten globule state of the chimeric protein is similar to that in the molten globule state of equine lysozyme in which packing interaction around the Trp residues in the native state is partially preserved. Therefore, the helix D in equine lysozyme is one of the contributing factors to the high stability and native-like packing interaction in the molten globule state of equine lysozyme. Our results indicate that the native-like packing interaction can stabilize the rudimentary intermediate which is stabilized by the non-specific hydrophobic interactions.     

Poisson-Boltzmann analysis of the lambda repressor-operator interaction.

A theoretical study of the ion atmosphere contribution to the binding free energy of the lambda repressor-operator complex is presented. The finite-difference form of the Poisson-Boltzmann equation was solved to calculate the electrostatic interaction energy of the amino-terminal domain of the lambda repressor with a 9 or 45 base pair oligonucleotide. Calculations were performed at various distances between repressor and operator as well as at different salt concentrations to determine ion atmosphere contributions to the total electrostatic interaction. Details in the distribution of charges on DNA and protein atoms had a strong influence on the calculated total interaction energies. In contrast, the calculated salt contributions are relatively insensitive to changes in the details of the charge distribution. The results indicate that the ion atmosphere contribution favors association at all protein-DNA distances studied. The theoretical number of ions released upon repressor-operator binding appears to be in reasonable agreement with experimental data.     

Mechanism of polyethylene glycol interaction with the molten globule folding intermediate of bovine carbonic anhydrase B.

Polyethylene glycol has been shown to bind to the molten globule intermediate on the bovine carbonic anhydrase B folding pathway. The mechanism of this interaction has been extensively probed. Polyethylene glycol (PEG) binds weakly to the molten globule first intermediate as measured by hydrophobic interaction chromatography, but PEG does not bind to either the native state or the second intermediate. The binding of PEG to the molten globule has been confirmed with both intrinsic fluorescence and fluorescence quenching experiments which indicate a single PEG-binding site on the molten globule. Electron paramagnetic resonance spectroscopic studies with nitroxide-labeled PEG also indicate a single binding site. Additional electron paramagnetic resonance studies with spin-labeled carbonic anhydrase B suggest that a conformational change occurs in the molten globule intermediate after PEG binds to the surface. The formation of a PEG-molten globule complex results in a reduction in self-association of this compact hydrophobic structure. PEG-molten globule complex formation is analogous to the observed interaction between chaperonins and a molten globule intermediate (Martin, J., Langer, T., Boteva, R., Schramel, A., Horwich, A.L., and Hartl, F.U. (1991) Nature 352, 36-42).     

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.     

Medium reorganization energy and enzymatic reaction activation energy

Reorganization and activation energies for charge transfer reactions occurring inside a dielectric sphere have been calculated by solving the problem of polar medium reorganization within and outside a dielectric sphere placed in another infinite dielectric. The dielectric sphere is assumed to simulate a protein globule, i.e. an enzyme molecule. It has been shown that for some reaction types the activation energy tends to decrease as the globule radius increases and that for each of the reaction types considered there is an optimal globule radius an increase of which does not bring about any tangible activation energy reduction. The calculated optimal radii for different processes are in good agreement with the increasing molecular sizes in the series: ribonuclease less than or equal to lysozyme less than serine proteinases approximately equal to cysteine proteinases less than NAD-dependent dehydrogenases. The calculated radii are usually about 1.5 to 1.7 times (and molecular masses about 4-5 times) smaller than the experimental ones. The reasons for this discrepancy are discussed and it has been suggested that the approximate nature of the treatment of a protein globule as a structureless dielectric is the main reason. It is shown that charge transfer at an acute angle to the globule surface is the optimum process. For endoergonic reaction stages it is the net charge transfer towards the periphery and for exoergonic ones that in the reverse direction which are advantageous. These conclusions are consistent with the data about the structure of the above-mentioned enzymes.     

Cofactor effects on the protein folding reaction: acceleration of alpha-lactalbumin refolding by metal ions

About 30% of proteins require cofactors for their proper folding. The effects of cofactors on the folding reaction have been investigated with alpha-lactalbumin as a model protein and metal ions as cofactors. Metal ions accelerate the refolding of alpha-lactalbumin by lessening the energy barrier between the molten globule state and the transition state, mainly by decreasing the difference of entropy between the two states. These effects are linked to metal ion binding to the protein in the native state. Hence, relationships between the metal affinities for the intermediate states and those for the native state are observed. Some residual specificity for the calcium ion is still observed in the molten globule state, this specificity getting closer in the transition state to that of the native state. The comparison between kinetic and steady-state data in association with the Phi value method indicates the binding of the metal ions on the unfolded state of alpha-lactalbumin. Altogether, these results provide insight into cofactor effects on protein folding. They also suggest new possibilities to investigate the presence of residual native structures in the unfolded state of protein and the effects of such structures on the protein folding reaction and on protein stability.     

Polyol-induced molten globule of cytochrome c: an evidence for stabilization by hydrophobic interaction.

To address the contribution of hydrophobic interaction to the stability of molten globule (MG) of proteins, the effects of various polyols (ethylene glycol, glycerol, erythritol, xylitol, sorbitol, and inositol) on the structure of acid-unfolded horse cytochrome c were examined at pH 2, by means of circular dichroism (CD), partial specific volume, adiabatic compressibility, and differential scanning calorimetry (DSC). Addition of polyols induced the characteristic CD spectra of MG, the effect being enhanced with an increase in their concentration and chain length (the number of OH groups) of polyols except for ethylene glycol. The free energy change of MG formation by sorbitol was comparable with those for the salt-induced MG formation but the heat capacity change was negligibly small. The partial specific volume did not change within the experimental error but the adiabatic compressibility largely increased by MG formation. The sorbitol-induced MG showed a highly cooperative DSC thermogram with a large heat capacity change in comparison with the salt-induced one. These results demonstrate that polyols can stabilize the MG state of this protein through the enhanced hydrophobic interaction overcoming the electrostatic repulsion between charged residues. The stabilizing mechanism and structure of MG state induced by polyols were discussed in terms of the preferential solvent interactions and osmotic pressure of the medium, in comparison with the salt-induced one.     

Cofactor-induced refolding: refolding of molten globule carbonic anhydrase induced by Zn(II) and Co(II)

The stability versus unfolding to the molten globule intermediate of bovine carbonic anhydrase II (BCA II) in guanidine hydrochloride (GuHCl) was found to depend on the metal ion cofactor [Zn(II) or Co(II)], and the apoenzyme was observed to be least stable. Therefore, it was possible to find a denaturant concentration (1.2 M GuHCl) at which refolding from the molten globule to the native state could be initiated merely by adding the metal ion to the apo molten globule. Thus, refolding could be performed without changing the concentration of the denaturant. The molten globule intermediate of BCA II could still bind the metal cofactor. Cofactor-effected refolding from the molten globule to the native state can be summarized as follows: (1) initially, the metal ion binds to the molten globule; (2) compaction of the metal-binding site region is then induced by the metal ion binding; (3) a functioning active center is formed; and (4) finally, the native tertiary structure is generated in the outer parts of the protein.     

Folding thermodynamics of model four-strand antiparallel beta-sheet proteins.

The thermodynamic properties for three different types of off-lattice four-strand antiparallel beta-strand protein models interacting via a hybrid Go-type potential have been investigated. Discontinuous molecular dynamic simulations have been performed for different sizes of the bias gap g, an artificial measure of a model protein's preference for its native state. The thermodynamic transition temperatures are obtained by calculating the squared radius of gyration R(g)(2), the root-mean-squared pair separation fluctuation Delta(B), the specific heat C(v), the internal energy of the system E, and the Lindemann disorder parameter Delta(L). Despite these models' simplicity, they exhibit a complex set of protein transitions, consistent with those observed in experimental studies on real proteins. Starting from high temperature, these transitions include a collapse transition, a disordered-to-ordered globule transition, a folding transition, and a liquid-to-solid transition. The high temperature transitions, i.e., the collapse transition and the disordered-to-ordered globule transition, exist for all three beta-strand proteins, although the native-state geometry of the three model proteins is different. However the low temperature transitions, i.e., the folding transition and the liquid-to-solid transition, strongly depend on the native-state geometry of the model proteins and the size of the bias gap.     

Dimensions,energetics, and denaturant effects of the protein unstructured state

Determining the energetics of the unfolded state of a protein is essential for understanding the folding mechanics of ordered proteins and the structure–function relation of intrinsically disordered proteins. Here, we adopt a coil‐globule transition theory to develop a general scheme to extract interaction and free energy information from single‐molecule fluorescence resonance energy transfer spectroscopy. By combining protein stability data, we have determined the free energy difference between the native state and the maximally collapsed denatured state in a number of systems, providing insight on the specific/nonspecific interactions in protein folding. Both the transfer and binding models of the denaturant effects are demonstrated to account for the revealed linear dependence of inter‐residue interactions on the denaturant concentration, and are thus compatible under the coil‐globule transition theory to further determine the dimension and free energy of the conformational ensemble of the unfolded state. The scaling behaviors and the effective θ‐state are also discussed.     

The small heat-shock chaperone protein, alpha-crystallin, does not recognize stable molten globule states of cytosolic proteins

The small heat-shock protein (sHsp), alpha-crystallin, acts as a molecular chaperone by interacting with destabilized 'substrate' proteins to prevent their precipitation from solution under conditions of stress. alpha-Crystallin and all sHsps are intracellular proteins. Similarly to other chaperones, the 'substrate' protein is in an intermediately folded, partly structured molten globule state when it interacts and complexes with alpha-crystallin. In this study, stable molten globule states of the cytosolic proteins, gamma-crystallin and myoglobin, have been prepared. Within the lens, gamma-crystallin naturally interacts with alpha-crystallin and myoglobin and alpha-crystallin are present together in muscle tissue. The molten globule states of gamma-crystallin and myoglobin were prepared by reacting gamma-crystallin with glucose 6-phosphate and by removing the haem group of myoglobin. Following spectroscopic characterisation of these modified proteins, their interaction with alpha-crystallin was examined by a variety of spectroscopic and protein chemical techniques. In both cases, there was no interaction with alpha-crystallin that led to complexation. It is concluded that alpha-crystallin does not recognise stable molten globule states of cytosolic 'substrate' proteins and only interacts with molten globule states of proteins that are on the irreversible pathway towards an aggregated and precipitated form.     

Effects of medium polarization and pre-existing field on activation energy of enzymatic charge-transfer reactions

The highly organized spatial structure of proteins' polar groups results in the existence of a permanent intraprotein electric field and in protein's weak dielectric response, i.e. its low dielectric constant. The first factor affects equilibrium free energy gap of a charge-transfer reaction, the second (medium polarization effect) influences both equilibrium and non-equilibrium (reorganization) energies, decreasing the latter substantially. In the framework of the rigorous 'fixed-charge-density' formalism, the medium polarization component of the reaction activation energy has been calculated, both for the activation energy of the elementary act proper, and the effective activation energy accounting for the charges' transfer from water into a low-dielectric structureless medium. In all typical cases of reactions, the energy spent for charge transfer from water into structureless 'protein' is larger than the gain in activation energy due to the protein's low reorganization energy. Therefore, the low dielectric constant of proteins is not sufficient to ensure their high catalytic activity, and an additional effect of the pre-existing intraprotein electric field, compensating for an excessive charging energy, is necessary. Only a combined action of low reorganization energy and pre-existing electric field provides proteins with their high catalytic activity. The dependence of activation energy on the globule geometry has been analyzed. It is shown that, for each reaction, an optimum set of geometric parameters exists. For five hydrolytic enzymes, the optimum globule radii have been calculated using the experimental geometry of their active sites. The calculated radii agree satisfactorily with the real sizes of these macromolecules, both by absolute and by relative values.     

Direct energy transfer to study the 3D structure of non-native proteins: AGH complex in molten globule state of apomyoglobin.

The direct energy transfer technique was modified and applied to probe the relative localization of apomyoglobin A-, G- and H-helixes, which are partly protected from deuterium exchange in the equilibrium molten globule state and in the molten globule-like kinetic intermediate. The non-radiative transfer of tryptophan electronic energy to 3-nitrotyrosine was studied in different conformational states of apomyoglobin (native, molten globule, unfolded) and interpreted in terms of average distances between groups of the protein chain. The experimental data show that the distance between the middle of A-helix and the N-terminus of G-helix as well as the distance between the middle of the A-helix and the C-terminus of the H-helix in the molten globule state are close to those in the native state. This is a strong argument in favor of similarity of the overall architecture of the molten globule and native states.     

Toroidal globular state of circular DNA

The influence of torsional elasticity of the double helix on compactization and structure of circular DNA in a compact form is studied in the case when the compact (globular) particle has a torus shape. For closed circular DNA the topological invariant, linking number of two strains, yields strict connection between conformation of double helix, considered as a unifilar homopolymer, and elastic energy of torsional twisting. The contribution of torsional elasticity to free energy of the toruslike globule is calculated. This contribution is shown to be proportional to the square of superturn's density. Torsional elasticity decreases the equilibrium radius of the toruslike globule formed by circular DNA in comparison with the case of linear DNA. Closure of linear DNA into a ring widens the stability range of the relatively short DNA compact form and tightens it for long DNA.     

Molten globule formation in apomyoglobin monitored by the fluorescent probe Nile Red

The interaction of nile red (NR) with apomyoglobin (ApoMb) in the native (pH 7) and molten globule (pH 4) states was investigated using experimental and computational methods. NR binds to hydrophobic locations in ApoMb with higher affinity (K(d) = 25 +/- 5 microM) in the native state than in the molten globule state (K(d) = 52 +/- 5 microM). In the molten globule state, NR is located in a more hydrophobic environment. The dye does not bind to the holoprotein, suggesting that the binding site is located at the heme pocket. In addition to monitoring steady-state properties, the fluorescence emission of NR is capable of tracking submillisecond, time-resolved structural rearrangements of the protein, induced by a nanosecond pH jump. Molecular dynamics simulations were run on ApoMb at neutral pH and at pH 4. The structure obtained for the molten globule state is consistent with the experimentally available structural data. The docking of NR with the crystal structure shows that the ligand binds into the binding pocket of the heme group, with an orientation bringing the planar ring system of NR to overlap with the position of two of the heme porphyrin rings in Mb. The docking of NR with the ApoMb structure at pH 4 shows that the dye binds to the heme pocket with a slightly less favorable binding energy, in keeping with the experimental K(d) value. Under these conditions, NR is positioned in a different orientation, reaching a more hydrophobic environment in agreement with the spectroscopic data.     

Conformational mobility of proteins and activation energy of enzymatic charge transfer reactions: an analysis in dielectric formalism

In calculating the medium reorganization energy and the activation energy of charge transfer enzymatic reactions, an allowance is made for the enhanced conformational mobility of the protein external region. The two-layer model is proposed, the outer layer having a higher static dielectric constant. The calculations show that the higher mobility in the outer layer causes some quantitative rather than qualitative changes. The main result obtained earlier is confirmed: the reorganization energy for charge transfer reaction in protein globule is much lower than in water and for this reason the activation energy also decreases. The higher dielectric constant of the outer layer somewhat favours the introduction of charge into active site and hence favours the natural selection of proteins as enzymes. This effect cannot exclude the necessity of other factors stabilizing ionic forms inside the protein globule. Freezing of conformational mobility (say, at low temperatures) hinders the charge transfer process as a consequence of the difficulty in equalizing the initial and final energy levels.     

Hierarchy of the interaction energy distribution in the spatial structure of globular proteins and the problem of domain definition.

An algorithm for determining of protein domain structure is proposed. Domain structures resulted from the algorithm application have been obtained and compared with available data. The method is based on entirely physical model of van der Waals interactions that reflects as illustrated in this work the distribution of electron density. Various levels of hierarchy in the protein spatial structure are discerned by analysis of the energy interaction between structural units of different scales. Thus the level of energy hierarchy plays role of sole parameter, and the method obviates the use of complicated geometrical criteria with numerous fitting parameters. The algorithm readily and accurately locates domains formed by continuous segments of the protein chain as well as those comprising non-sequential segments, sets no limit to the number of segments in a domain. We have analyzed 309 protein structures. Among 277 structures for which our results could be compared with the domain definitions made in other works, 243 showed complete or partial coincidence, and only in 34 cases the domain structures proved substantially different. The domains delineated with our approach may coincide with reference definition at different levels of the globule hierarchy. Along with defining the domain structure, our approach allows one to consider the protein spatial structure in terms of the spatial distribution of the interaction energy in order to establish the correspondence between the hierarchy of energy distribution and the hierarchy of structural elements.     

On the possibility for a protein to exist in a manifold of partially folded states

With the help of the methods of tryptophan fluorescence and room-temperature phosphorescence and using Escherichia coli alkaline phosphatase as an example, the ability of a protein to exist in a manifold of partially folded thermodynamically stable states differing in conformation, the internal dynamics, and functional activity was shown. Such intermediate (between native and unfolded) structures may form during unfolding or folding of the protein. It was shown that the degree of destruction of the native structural organization of the globule depends on both the nature and the mode of action of the destroying agent and the structure of the protein. Conformational transitions of the globule can change the kind of the internal dynamics (fast, slow), and shifts of dynamics can initiate conformation changes of the protein and precede them. A scheme of the structural and functional transformations of the protein during denaturation is presented, which takes into account the possibility of globule transitions into a manifold of functional active and inactive partially folded states. The role of partially folded forms of cell proteins in the development of pathology is discussed.     

Theory of excitable membranes. I. A simple model for a three-state artificial membrane.

The effect of the thermal fluctuation on the orientation distribution pattern of globular protein molecules in a two-dimensional lattice was investigated by the method of computer simulation. A set of interaction parameters was assigned to interaction sites on each molecule and the interaction energy between two molecules was given by the product of the parameters of facing sites. Orientation fluctuation was assumed to take place with the probability proportional to the Boltzmann factor. Patterns having different degrees of order appeared with the change of temperature. The entropy and other thermodynamic quantities of these patterns were calculated, and gradual and transitional changes of the pattern were discussed in comparison with the case of simple atoms or molecules.