Molecular basis of cooperativity in protein folding IV. Core: A general cooperative folding model

The cooperative nature of the protein folding process is independent of the characteristic fold and the specific secondary structure attributes of a globular protein. A general folding/unfolding model should, therefore, be based upon structural features that transcend the peculiarities of α-helices, β-sheets, and other structural motifs found in proteins. The studies presented in this paper suggest that a single structural characteristic common to all globular proteins is essential for cooperative folding. The formation of a partly folded state from the native state results in the exposure to solvent of two distinct regions: (1) the portions of the protein that are unfolded; and (2) the “complementary surfaces,” located in the regions of the protein that remain folded. The cooperative character of the folding/unfolding transition is determined largely by the energetics of exposing complementary surface regions to the solvent. By definition, complementary regions are present only in partly folded states; they are absent from the native and unfolded states. An unfavorable free energy lowers the probability of partly folded states and increases the cooperativity of the transition. In this paper we present a mathematical formulation of this behavior and develop a general cooperative folding/unfolding model, termed the “complementary region” (CORE) model. This model successfully reproduces the main properties of folding/unfolding transitions without limiting the number of partly folded states accessible to the protein, thereby permitting a systematic examination of the structural and solvent conditions under which intermediates become populated. It is shown that the CORE model predicts two-state folding/unfolding behavior, even though the two-state character is not assumed in the model. © 1993 Wiley-Liss, Inc.     

Molecular basis of cooperativity in protein folding. V. Thermodynamic and structural conditions for the stabilization of compact denatured states

The heat-denatured state of proteins has been usually assumed to be a fully hydrated random coil. It is now evident that under certain solvent conditions or after chemical or genetic modifications, the protein molecule may exhibit a hydrophobic core and residual secondary structure after thermal denaturation. This state of the protein has been called the “compact denatured” or “molten globule” state. Recently is has been shown that α-lactalbumin at pH < 5 denatures into a molten globule state upon increasing the temperature (Griko, Y., Freire, E., Privalov, P. L. Biochemistry 33:1889–1899, 1994). This state has a lower heat capacity and a higher enthalpy at low temperatures than the unfolded state. At those temperatures the stabilization of the molten globule state is of an entropic origin since the enthalpy contributes unfavorably to the Gibbs free energy. Since the molten globule is more structured than the unfolded state and, therefore, is expected to have a lower configurational entropy, the net entropic gain must originate primarily from solvent related entropy arising from the hydrophobic effect, and to a lesser extent from protonation or electrostatic effects. In this work, we have examined a large ensemble of partly folded states derived from the native structure of α-lactalbumin in order to identify those states that satisfy the energetic criteria of the molten globule. It was found that only few states satisfied the experimental constraints and that, furthermore, those states were part of the same structural family. In particular, the regions corresponding to the A, B, and C helices were found to be folded, while the β sheet and the D helix were found to be unfolded. At temperatures below 45°C the states exhibiting those structural characteristics are enthalpically higher than the unfolded state in agreement with the experimental data. Interestingly, those states have a heat capacity close to that observed for the acid pH compact denatured state of α-lactalbumin [980 cal (mol.K)^?l]. In addition, the folded regions of these states include those residues found to be highly protected by NMR hydrogen exchange experiments. This work represents an initial attempt to model the structural origin of the thermodynamic properties of partly folded states. The results suggest a number of structural features that are consistent with experimental data. © 1994 Wiley-Liss, Inc.     

Crystal structure of AlgQ2, a macromolecule (alginate)-binding protein of Sphingomonas sp. A1 at 2.0A resolution

Sphingomonas sp. A1 possesses a high molecular mass (average 25,700 Da) alginate uptake system mediated by a novel pit-dependent ABC transporter. The X-ray crystallographic structure of AlgQ2 (57,200 Da), an alginate-binding protein in the system, was determined by the multiple isomorphous replacement method and refined at 2.0 A resolution with a final R-factor of 18.3% for 15 to 2.0 A resolution data. The refined structure of AlgQ2 was comprised of 492 amino acid residues, 172 water molecules, and one calcium ion. AlgQ2 was composed of two globular domains with a deep cleft between them, which is expected to be the alginate-binding site. The overall structure is basically similar to that of maltose/maltodextrin-binding protein, except for the presence of an N2-subdomain. The entire calcium ion-binding site is similar to the site in the EF-hand motif, but comprises a ten residue loop. This calcium ion-binding site is about 40 A away from the alginate-binding site.     

Conformational gating of dimannose binding to the antiviral protein cyanovirin revealed from the crystal structure at 1.35 A resolution

Cyanovirin (CV-N) is a small lectin with potent HIV neutralization activity, which could be exploited for a mucosal defense against HIV infection. The wild-type (wt) protein binds with high affinity to mannose-rich oligosaccharides on the surface of gp120 through two quasi-symmetric sites, located in domains A and B. We recently reported on a mutant of CV-N that contained a single functional mannose-binding site, domain B, showing that multivalent binding to oligomannosides is necessary for antiviral activity. The structure of the complex with dimannose determined at 1.8 A resolution revealed a different conformation of the binding site than previously observed in the NMR structure of wt CV-N. Here, we present the 1.35 A resolution structure of the complex, which traps three different binding conformations of the site and provides experimental support for a locking and gating mechanism in the nanoscale time regime observed by molecular dynamics simulations.     

Role of electrostatic interactions for the stability and folding behavior of cold shock protein

The impacts of three charged‐residue‐involved mutations, E46A, R3E, and R3E/L66E, on the thermostability and folding behavior of the cold shock protein from the themophile Bacillus caldolyticus (Bc‐Csp) were investigated by using a modified Gō‐like model, in which the nonspecific electrostatic interactions of charged residues were taken into account. Our simulation results show that the wild‐type Bc‐Csp and its three mutants are all two‐sate folders, which is consistent with the experimental observations. It is found that these three mutations all lead to a decrease of protein thermodynamical stability, and the effect of R3E mutation is the strongest. The lower stability of these three mutants is due to the increase of the enthalpy of the folded state and the entropy of the unfolded state. Using this model, we also studied the folding kinetics and the folding/unfolding pathway of the wild‐type Bc‐Csp as well as its three mutants and then discussed the effects of electrostatic interactions on the folding kinetics. The results indicate that the substitutions at positions 3 and 46 largely decrease the folding kinetics, whereas the mutation of residue 66 only slightly decreases the folding rate. This result agrees well with the experimental observations. It is also found that these mutations have little effects on the folding transition state and the folding pathway, in which the N‐terminal β sheet folds earlier than the C‐terminal region. We also investigated the detailed unfolding pathway and found that it is really the reverse of the folding pathway, providing the validity of our simulation results. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Crystal structure of human coactosin-like protein at 1.9 A resolution

Human coactosin-like protein (CLP) shares high homology with coactosin, a filamentous (F)-actin binding protein, and interacts with 5LO and F-actin. As a tumor antigen, CLP is overexpressed in tumor tissue cells or cell lines, and the encoded epitopes can be recognized by cellular and humoral immune systems. To gain a better understanding of its various functions and interactions with related proteins, the crystal structure of CLP expressed in Escherichia coli has been determined to 1.9 A resolution. The structure features a central beta-sheet surrounded by helices, with two very tight hydrophobic cores on each side of the sheet. CLP belongs to the actin depolymerizing protein superfamily, and is similar to yeast cofilin and actophilin. Based on our structural analysis, we observed that CLP forms a polymer along the crystallographic b axis with the exact same repeat distance as F-actin. A model for the CLP polymer and F-actin binding has therefore been proposed.     

The crystal structure of staphylococcal nuclease refined at 1.7 A resolution

The crystal structure of staphylococcal nuclease has been determined to 1.7 A resolution with a final R-factor of 16.2% using stereochemically restrained Hendrickson-Konnert least-squares refinement. The structure reveals a number of conformational changes relative to the structure of the ternary complex of staphylococcal nuclease 1,2 bound with deoxythymidine-3',5'-diphosphate and Ca2+. Tyr-113 and Tyr-115, which pack against the nucleotide base in the nuclease complex, are rotated outward creating a more open binding pocket in the absence of nucleotide. The side chains of Ca2+ ligands Asp-21 and Asp-40 shift as does Glu-43, the proposed general base in the hydrolysis of the 5'-phosphodiester bond. The significance of some changes in the catalytic site is uncertain due to the intrusion of a symmetry related Lys-70 side chain which hydrogen bonds to both Asp-21 and Glu-43. The position of a flexible loop centered around residue 50 is altered, most likely due to conformational changes propagated from the Ca2+ site. The side chains of Arg-35, Lys-84, Tyr-85, and Arg-87, which hydrogen bond to the 3'- and 5'-phosphates of the nucleotide in the nuclease complex, are unchanged in conformation, with packing interactions with adjacent protein side chains sufficient to fix the geometry in the absence of ligand. The nuclease structure presented here, in combination with the stereochemically restrained refinement of the nuclease complex structure at 1.65 A, provides a wealth of structural information for the increasing number of studies using staphylococcal nuclease as a model system of protein structure and function.     

Apoflavodoxin (un)folding followed at the residue level by NMR

The denaturant-induced (un)folding of apoflavodoxin from Azotobacter vinelandii has been followed at the residue level by NMR spectroscopy. NH groups of 21 residues of the protein could be followed in a series of 1H-15N heteronuclear single-quantum coherence spectra recorded at increasing concentrations of guanidinium hydrochloride despite the formation of protein aggregate. These NH groups are distributed throughout the whole apoflavodoxin structure. The midpoints of unfolding determined by NMR coincide with the one obtained by fluorescence emission spectroscopy. Both techniques give rise to unfolding curves with transition zones at significantly lower denaturant concentrations than the one obtained by circular dichroism spectroscopy. The NMR (un)folding data support a mechanism for apoflavodoxin folding in which a relatively stable intermediate is involved. Native apoflavodoxin is shown to cooperatively unfold to a molten globule-like state with extremely broadened NMR resonances. This initial unfolding step is slow on the NMR chemical shift timescale. The subsequent unfolding of the molten globule is faster on the NMR chemical shift timescale and the limited appearance of 1H-15N HSQC cross peaks of unfolded apoflavodoxin in the denaturant range studied indicates that it is noncooperative.     

The 2.0 A resolution crystal structure of prostaglandin H2 synthase-1: structural insights into an unusual peroxidase

Prostaglandin H2 synthase (EC 1.14.99.1) is an integral membrane enzyme containing a cyclooxygenase site, which is the target for the non-steroidal anti-inflammatory drugs, and a spatially distinct peroxidase site. Previous crystallographic studies of this clinically important drug target have been hindered by low resolution. We present here the 2.0 A resolution X-ray crystal structure of ovine prostaglandin H2 synthase-1 in complex with alpha-methyl-4-biphenylacetic acid, a defluorinated analog of the non-steroidal anti-inflammatory drug flurbiprofen. Detergent molecules are seen to bind to the protein's membrane-binding domain, and their positions suggest the depth to which this domain is likely to penetrate into the lipid bilayer. The relation of the enzyme's proximal heme ligand His388 to the heme iron is atypical for a peroxidase; the iron-histidine bond is unusually long and a substantial tilt angle is observed between the heme and imidazole planes. A molecule of glycerol, used as a cryoprotectant during diffraction experiments, is seen to bind in the peroxidase site, offering the first view of any ligand in this active site. Insights gained from glycerol binding may prove useful in the design of a peroxidase-specific ligand.     

The structure of the major transition state for folding of an FF domain from experiment and simulation

We have analysed the transition state of folding of the four-helix FF domain from HYPA/FBP11 by high-resolution experiment and simulation as part of a continuing effort to understand the principles of folding and the refinement of predictive methods. The major transition state for folding was subjected to a Phi-value analysis utilising 50 mutants. The transition state contained a nucleus for folding centred around the end of helix 1 (H1) and the beginning of helix 2 (H2). Secondary structure in this region was fully formed (PhiF=0.9-1) and tertiary interactions were well developed. Interactions in the distal part of the native structure were weak (PhiF=0-0.2). The hydrophobic core and other parts of the protein displayed intermediate Phi-values, suggesting that interactions coalesce as the end of H1 and beginning of H2 are in the process of being formed. The distribution of Phi-values resembled that of barnase, which folds via an intermediate, rather than that of CI2 which folds by a concerted nucleation-condensation mechanism. The overall picture of the transition state structure identified in molecular dynamics simulations is in qualitative agreement, with the turn connecting H1 and H2 being formed, a loosened core, and H4 partially unfolded and detached from the core. There are some differences in the details and interpretation of specific Phi-values.     

The crystal structure of the olfactory marker protein at 2.3 A resolution

Olfactory marker protein (OMP) is a highly expressed and phylogenetically conserved cytoplasmic protein of unknown function found almost exclusively in mature olfactory sensory neurons. Electrophysiological studies of olfactory epithelia in OMP knock-out mice show strongly retarded recovery following odorant stimulation leading to an impaired response to pulsed odor stimulation. Although these studies show that OMP is a modulator of the olfactory signal-transduction cascade, its biochemical role is not established. In order to facilitate further studies on the molecular function of OMP, its crystal structure has been determined at 2.3 A resolution using multiwavelength anomalous diffraction experiments on selenium-labeled protein. OMP is observed to form a modified beta-clamshell structure with eight antiparallel beta-strands. While OMP has no significant sequence homology to proteins of known structure, it has a similar fold to a domain found in a variety of existing structures, including in a large family of viral capsid proteins. The surface of OMP is mostly convex and lacking obvious small molecule binding sites, suggesting that it is more likely to be involved in modulating protein-protein interaction than in interacting with small molecule ligands. Three highly conserved regions have been identified as leading candidates for protein-protein interaction sites in OMP. One of these sites represents a loop known to mediate ligand interactions in the structurally homologous EphB2 receptor ligand-binding domain. This site is partially buried in the crystal structure but fully exposed in the NMR solution structure of OMP due to a change in the orientation of an alpha-helix that projects outward from the structurally invariant beta-clamshell core. Gating of this conformational change by molecular interactions in the signal-transduction cascade could be used to control access to OMP's equivalent of the EphB2 ligand-interaction loop, thereby allowing OMP to function as a molecular switch.     

The molecular structure of UDP-galactose 4-epimerase from Escherichia coli determined at 2.5 A resolution.

UDP-galactose 4-epimerase catalyzes the conversion of UDP-galactose to UDP-glucose during normal galactose metabolism. The molecular structure of UDP-galactose 4-epimerase from Escherichia coli has now been solved to a nominal resolution of 2.5 A. As isolated from E. coli, the molecule is a dimer of chemically identical subunits with a total molecular weight of 79,000. Crystals of the enzyme used for this investigation were grown as a complex with the substrate analogue, UDP-benzene, and belonged to the space group P2(1)2(1)2(1) with unit cell dimensions of a = 76.3 A, b = 83.1 A, c = 132.1 A, and one dimer per asymmetric unit. An interpretable electron density map calculated to 2.5 A resolution was obtained by a combination of multiple isomorphous replacement with six heavy atom derivatives, molecular averaging, and solvent flattening. Each subunit of epimerase is divided into two domains. The larger N-terminal domain, composed of amino acid residues 1-180, shows a classic NAD+ binding motif with seven strands of parallel beta-pleated sheet flanked on either side of alpha-helices. The seventh strand of the beta-pleated sheet is contributed by amino acid residues from the smaller domain. In addition, this smaller C-terminal domain, consisting of amino acid residues 181-338, contains three strands of beta-pleated sheet, two major alpha-helices and one helical turn. The substrate analogue, UDP-benzene, binds in the cleft located between the two domains with its phenyl ring in close proximity to the nicotinamide ring of NAD+. Contrary to the extensive biochemical literature suggesting that epimerase binds only one NAD+ per functional dimer, the map clearly shows electron density for two nicotinamide cofactors binding in symmetry-related positions in the dimer. Likewise, each subunit in the dimer also binds one substrate analogue.     

Crystal structure of an Escherichia coli tRNA(Gly) microhelix at 2.0 A resolution

tRNA identity elements determine the correct aminoacylation by the cognate aminoacyl-tRNA synthetase. In class II aminoacyl tRNA synthetase systems, tRNA specificity is assured by rather few and simple recognition elements, mostly located in the acceptor stem of the tRNA. Here we present the crystal structure of an Escherichia coli tRNA(Gly) aminoacyl stem microhelix at 2.0 A resolution. The tRNA(Gly) microhelix crystallizes in the space group P3(2)21 with the cell constants a=b=35.35 A, c=130.82 A, gamma=120 degrees . The helical parameters, solvent molecules and a potential magnesium binding site are discussed.     

A thermodynamic and kinetic analysis of the folding pathway of an SH3 domain entropically stabilised by a redesigned hydrophobic core

The folding thermodynamics and kinetics of the alpha-spectrin SH3 domain with a redesigned hydrophobic core have been studied. The introduction of five replacements, A11V, V23L, M25V, V44I and V58L, resulted in an increase of 16% in the overall volume of the side-chains forming the hydrophobic core but caused no remarkable changes to the positions of the backbone atoms. Judging by the scanning calorimetry data, the increased stability of the folded structure of the new SH3-variant is caused by entropic factors, since the changes in heat capacity and enthalpy upon the unfolding of the wild-type and mutant proteins were identical at 298 K. It appears that the design process resulted in an increase in burying both the hydrophobic and hydrophilic surfaces, which resulted in a compensatory effect upon the changes in heat capacity and enthalpy. Kinetic analysis shows that both the folding and unfolding rate constants are higher for the new variant, suggesting that its transition state becomes more stable compared to the folded and unfolded states. The phi(double dagger-U) values found for a number of side-chains are slightly lower than those of the wild-type protein, indicating that although the transition state ensemble (TSE) did not change overall, it has moved towards a more denatured conformation, in accordance with Hammond's postulate. Thus, the acceleration of the folding-unfolding reactions is caused mainly by an improvement in the specific and/or non-specific hydrophobic interactions within the TSE rather than by changes in the contact order. Experimental evidence showing that the TSE changes globally according to its hydrophobic content suggests that hydrophobicity may modulate the kinetic behaviour and also the folding pathway of a protein.     

Structural motifs for pyridoxal-5'-phosphate binding in decarboxylases: an analysis based on the crystal structure of the Lactobacillus 30a ornithine decarboxylase.

Two of the five domains in the structure of the ornithine decarboxylase (OrnDC) from Lactobacillus 30a share similar structural folds around the pyridoxal-5''-phosphate (PLP)-binding pocket with the aspartate aminotransferases (AspATs). Sequence comparisons focusing on conserved residues of the aligned structures reveal that this structural motif is also present in a number of other PLP-dependent enzymes including the histidine, dopa, tryptophan, glutamate, and glycine decarboxylases as well as tryptophanase and serine-hydroxymethyl transferase. However, this motif is not present in eukaryotic OrnDCs, the diaminopimelate decarboxylases, nor the Escherichia coli or oat arginine decarboxylases. The identification and comparison of residues involved in defining the different classes are discussed.     

Molecular structure of the bilin binding protein (BBP) from Pieris brassicae after refinement at 2.0 A resolution

The bilin binding protein (BBP) from the insect Pieris brassicae has been analysed for amino acid sequence, spectral properties and three-dimensional structure. The crystal structure that had been determined by isomorphous replacement has been refined at 2.0 A (1 A = 0.1 nm) resolution to an R-value of 0.20. The asymmetric unit contains four independent subunits of BBP. The co-ordinate differences are 0.25 A, in accord with the estimated error in co-ordinates. The polypeptide chain fold is characterized by an eight-stranded barrel. The connecting loops splay out at the upper end of the barrel and open it, whilst the lower end is closed. The overall shape resembles a calyx. The biliverdin IX gamma chromophore is located in a central cleft at the upper end of the barrel. The bilatriene moiety is in cyclic helical geometry with configuration Z,Z,Z and conformation syn,syn,syn. The geometry is in accord with the spectral properties and permits a correlation between sign of the circular dichroism bands and sense of the bilatriene helices. The fold of BBP is related to retinol binding protein (RBP), as had been recognized in the preliminary analysis, although the amino acid sequences of RBP and BBP show only 10% homology. There are large differences in the loops at the upper end of the barrel, whilst the segments of the centre and the lower end of the barrel superimpose closely. The ligands of BBP and RBP, biliverdin and retinol, respectively, are also similarly located.     

The structure of the N-terminal domain of riboflavin synthase in complex with riboflavin at 2.6A resolution

Riboflavin synthase of Escherichia coli is a homotrimer with a molecular mass of 70 kDa. The enzyme catalyzes the dismutation of 6,7-dimethyl-8-(1'-D-ribityl)-lumazine, affording riboflavin and 5-amino-6-ribitylamino-2,4(1H,3H)-pyrimidinedione. The N-terminal segment (residues 1-87) and the C-terminal segment (residues 98-187) form beta-barrels with similar fold and a high degree of sequence similarity. A recombinant peptide comprising amino acid residues 1-97 forms a dimer, which binds riboflavin with high affinity. Here, we report the structure of this construct in complex with riboflavin at 2.6A resolution. It is demonstrated that the complex can serve as a model for ligand-binding in the native enzyme. The structure and riboflavin-binding mode is in excellent agreement with structural information obtained from the native enzyme from Escherichia coli and riboflavin synthase from Schizosaccharomyces pombe. The implications for the binding specificity and the regiospecificity of the catalyzed reaction are discussed.     

Prediction of folding transition-state position (betaT) of small, two-state proteins from local secondary structure content

Folding kinetics of proteins is governed by the free energy and position of transition states. But attempts to predict the position of folding transition state on reaction pathway from protein structure have been met with only limited success, unlike the folding-rate prediction. Here, we find that the folding transition-state position is related to the secondary structure content of native two-state proteins. We present a simple method for predicting the transition-state position from their alpha-helix, turn and polyproline secondary structures. The method achieves 81% correlation with experiment over 24 small, two-state proteins, suggesting that the local secondary structure content, especially for content of alpha-helix, is a determinant of the solvent accessibility of the transition state ensemble and size of folding nucleus.     

Conformational plasticity in folding of the split beta-alpha-beta protein S6: evidence for burst-phase disruption of the native state

An increasing number of folding studies of two-state proteins shows that point mutations sometimes change the kinetic m-values, leading to kinks and curves in the chevron plots. The molecular origin of these changes is yet unclear although it is speculated that they are linked to structural rearrangement of the transition state or to accumulation of meta-stable intermediates. To shed more light on this issue, we present here a combined m and phi-value analysis of the split beta-alpha-beta protein S6. Wild-type S6 displays classical two-state kinetics with v-shaped chevron plot, but a majority of its mutants display distinct m-value changes or curved chevrons. We observe that this kinetic aberration of S6 is linked to mutations that are clustered in distinct regions of the native structure. The most pronounced changes, i.e. decrease in the m-value for the unfolding rate constant, are seen upon truncation of interactions between the N and C termini, whereas mutations in the centre of the hydrophobic core show smaller or even opposed effects. As a consequence, the calculated phi-values display a systematic increase upon addition of denaturant. In the case of S6, the phenomenon seems to arise from a general plasticity of the different species on the folding pathway. That is, the structure of the denatured ensemble, the transition state, and the native ground-state for unfolding seem to change upon mutation. From these changes, it is concluded that interactions spanning the centre of the hydrophobic core form early in folding, whereas the entropically disfavoured interactions linking the N and C termini consolidate very late, mainly on the down-hill-side of the folding barrier.