High helicity of peptide fragments corresponding to beta-strand regions of beta-lactoglobulin observed by 2D-NMR spectroscopy

Polarizable and non-polarizable force fields: Protein folding, unfolding, and misfolding

Molecular dynamics simulations are an invaluable tool to characterize the dynamic motions of proteins in atomistic detail. However, the accuracy of models derived from simulations inevitably relies on the quality of the underlying force field. Here, we present an evaluation of current non-polarizable and polarizable force fields (AMBER ff14SB, CHARMM 36m, GROMOS 54A7, and Drude 2013) based on the long-standing biophysical challenge of protein folding. We quantify the thermodynamics and kinetics of the β-hairpin formation using Markov state models of the fast-folding mini-protein CLN025. Furthermore, we study the (partial) folding dynamics of two more complex systems, a villin headpiece variant and a WW domain. Surprisingly, the polarizable force field in our set, Drude 2013, consistently leads to destabilization of the native state, regardless of the secondary structure element present. All non-polarizable force fields, on the other hand, stably characterize the native state ensembles in most cases even when starting from a partially unfolded conformation. Focusing on CLN025, we find that the conformational space captured with AMBER ff14SB and CHARMM 36m is comparable, but the ensembles from CHARMM 36m simulations are clearly shifted toward disordered conformations. While the AMBER ff14SB ensemble overstabilizes the native fold, CHARMM 36m and GROMOS 54A7 ensembles both agree remarkably well with experimental state populations. In addition, GROMOS 54A7 also reproduces experimental folding times most accurately. Our results further indicate an over-stabilization of helical structures with AMBER ff14SB. Nevertheless, the presented investigations strongly imply that reliable (un)folding dynamics of small proteins can be captured in feasible computational time with current additive force fields.

Effect of Antioxidants on Heavy Metals Induced Conformational Alteration of Cytochrome C and Myoglobin

BACKGROUND

The exposure to heavy metals due to unrestrained industrialization, pollution and non-degradability imposes a significant risk to human health. Proteins are prime targets of heavy metal stress, however, the underlying mechanisms and its impact on heme proteins is still not entirely clear.

OBJECTIVE

To analyze the deleterious effect of heavy metals such as cadmium, chromium and mercury on conformation of two proteins namely, cytochrome c and myoglobin. The protective effect of glycine and ascorbic acid (animal origin), gallic acid and sesamol (plant origin) on heavy metal exposure was studied.

METHODS

Far- and near-UV Circular Dichroism (CD) measurements monitored the changes in secondary and tertiary structure. Absorption Soret spectroscopy study revealed changes in heme-protein interaction. Peroxidase activity has been assayed to measure the absorption of tetraguaiacol. The interaction of heme proteins with different heavy metals was done using docking study.

RESULTS

Far- and near-UV CD measurements reveal that heavy metals disrupt the secondary and tertiary structure of heme proteins. Antioxidants counteract the deleterious effect of heavy metals. Absorption spectroscopy revealed changes in the Soret region of these heme proteins. Changes in peroxidase activity was observed on addition of heavy metals and antioxidants. Molecular docking validated interaction of the heavy metals with proteins with a significant binding affinity (-2.3 kcal/- mol).

CONCLUSION

Heavy metals interfered and disrupted both the heme proteins and mercury showed the maximum deleterious effect, further, chromium showed detrimental effect at very small concentration. The antioxidants from animal origin exhibited better protective response than those from plant source.

Toward Ab Initio Protein Folding: Inherent Secondary Structure Propensity of Short Peptides from the Bioinformatics and Quantum-Chemical Perspective

By combining bioinformatics with quantum-chemical calculations, we attempt to address quantitatively some of the physical principles underlying protein folding. The former allowed us to identify tripeptide sequences in existing protein three-dimensional structures with a strong preference for either helical or extended structure. The selected representatives of pro-helical and pro-extended sequences were converted into "isolated" tripeptides-capped at N- and C-termini-and these were subjected to an extensive conformational sampling and geometry optimization (typically thousands to tens of thousands of conformers for each tripeptide). For each conformer, the QM(DFT-D3)/COSMO-RS free-energy value was then calculated, G_conf(solv). The Δ G_conf(solv) is expected to provide an objective, unbiased, and quantitatively accurate measure of the conformational preference of the particular tripeptide sequence. It has been shown that irrespective of the helical vs extended preferences of the selected tripeptide sequences in context of the protein, most of the low-energy conformers of isolated tripeptides prefer the R-helical structure. Nevertheless, pro-helical tripeptides show slightly stronger helix preference than their pro-extended counterparts. Furthermore, when the sampling is repeated in the presence of a partner tripeptide to mimic the situation in a β-sheet, pro-extended tripeptides (exemplified by the VIV) show a larger free-energy benefit than pro-helical tripeptides (exemplified by the EAM). This effect is even more pronounced in a hydrophobic solvent, which mimics the less polar parts of a protein. This is in line with our bioinformatic results showing that the majority of pro-extended tripeptides are hydrophobic. The preference for a specific secondary structure by the studied tripeptides is thus governed by the plasticity to adopt to its environment. In addition, we show that most of the "naturally occurring" conformations of tripeptide sequences, i.e., those found in existing three-dimensional protein structures, are within ∼10 kcal·mol^-1 from their global minima. In summary, our "ab initio" data suggest that complex protein structures may start to emerge already at the level of their small oligopeptidic units, which is in line with a hierarchical nature of protein folding.

Non-Native α-Helices in the Initial Folding Intermediate Facilitate the Ordered Assembly of the β-Barrel in β-Lactoglobulin

The roles of non-native α-helices frequently observed in the initial folding stage of β-sheet proteins have been examined for many years. We herein investigated the residue-level structures of several mutants of bovine β-lactoglobulin (βLG) in quenched-flow pH-pulse labeling experiments. βLG assumes a collapsed intermediate with a non-native α-helical structure (I₀) in the early stage of folding, although its native form is predominantly composed of β-structures. The protection profile in I₀ of pseudo-wild type (WT*) βLG was found to deviate from the pattern of the "average area buried upon folding" (AABUF). In particular, the level of protection at the region of strand A, at which non-native α-helices form in the I₀ state, was significantly low compared to AABUF. G17E, the mutant with an increased helical propensity, showed a similar protection pattern. In contrast, the protection pattern for I₀ of E44L, the mutant with an increased β-sheet propensity, was distinct from that of WT* and resembled the AABUF pattern. Transverse relaxation measurements demonstrated that the positions of the residual structures in the unfolded states of these mutants were consistent with those of the protected residues in the respective I₀ states. On the basis of the slower conversion of I₀ to the native state for E44L to that for WT*, non-native α-helices facilitate the ordered assembly of the β-barrel by preventing interactions that trap folding.

Higher-order structure of the Rous sarcoma virus SP assembly domain

Purified retroviral Gag proteins can assemble in vitro to form immature virus-like particles (VLPs). By cryoelectron tomography, Rous sarcoma virus VLPs show an organized hexameric lattice consisting chiefly of the capsid (CA) domain, with periodic stalk-like densities below the lattice. We hypothesize that the structure represented by these densities is formed by amino acid residues immediately downstream of the folded CA, namely, the short spacer peptide SP, along with a dozen flanking residues. These 24 residues comprise the SP assembly (SPA) domain, and we propose that neighboring SPA units in a Gag hexamer coalesce to form a six-helix bundle. Using in vitro assembly, alanine scanning mutagenesis, and biophysical analyses, we have further characterized the structure and function of SPA. Most of the amino acid residues in SPA could not be mutated individually without abrogating assembly, with the exception of a few residues near the N and C termini, as well as three hydrophilic residues within SPA. We interpret these results to mean that the amino acids that do not tolerate mutations contribute to higher-order structures in VLPs. Hydrogen-deuterium exchange analyses of unassembled Gag compared that of assembled VLPs showed strong protection at the SPA region, consistent with a higher-order structure. Circular dichroism revealed that a 29mer SPA peptide shifts from a random coil to a helix in a concentration-dependent manner. Analytical ultracentrifugation showed concentration-dependent self-association of the peptide into a hexamer. Taken together, these results provide strong evidence for the formation of a critical six-helix bundle in Gag assembly.

IMPORTANCE

The structure of a retrovirus like HIV is created by several thousand molecules of the viral Gag protein, which assemble to form the known hexagonal protein lattice in the virus particle. How the Gag proteins pack together in the lattice is incompletely understood. A short segment of Gag known to be critical for proper assembly has been hypothesized to form a six-helix bundle, which may be the nucleating event that leads to lattice formation. The experiments reported here, using the avian Rous sarcoma virus as a model system, further define the nature of this segment of Gag, show that it is in a higher-order structure in the virus particle, and provide the first direct evidence that it forms a six-helix bundle in retrovirus assembly. Such knowledge may provide underpinnings for the development of antiretroviral drugs that interfere with virus assembly.

Transient non-native helix formation during the folding of β-lactoglobulin

In ideal proteins, only native interactions are stabilized step-by-step in a smooth funnel-like energy landscape. In real proteins, however, the transient formation of non-native structures is frequently observed. In this review, the transient formation of non-native structures is described using the non-native helix formation during the folding of β-lactoglobulin as a prominent example. Although β-lactoglobulin is a predominantly β-sheet protein, it has been shown to form non-native helices during the early stage of folding. The location of non-native helices, their stabilization mechanism, and their role in the folding reaction are discussed.

Spectroscopic and theoretical investigation of oxali-palladium interactions with β-lactoglobulin

The possibility of using a small cheap dairy protein, β-lactoglobulin (β-LG), as a carrier for oxali-palladium for drug delivery was studied. Their binding in an aqueous solution at two temperatures of 25 and 37°C was investigated using spectroscopic techniques in combination with a molecular docking study. Fluorescence intensity changes showed combined static and dynamic quenching during β-LG oxali-palladium binding, with the static mode being predominant in the quenching mechanism. The binding and thermodynamic parameters were determined by analyzing the results of quenching and those of the van't Hoff equation. According to obtained results the binding constants at two temperatures of 25 and 37°C are 3.3×10(9) M(-1) and 18.4×10(6) M(-1) respectively. Fluorescence resonance energy transfer (FRET) showed that the experimental results and the molecular docking results were coherent. An absence change of β-LG secondary structure was confirmed by the CD results. Molecular docking results agreed fully with the experimental results since the fluorescence studies also revealed the presence of two binding sites with a negative value for the Gibbs free energy of binding of oxali-palladium to β-LG. Furthermore, molecular docking and experimental results suggest that the hydrophobic effect plays a critical role in the formation of the oxali-palladium complex with β-LG. This agreement between molecular docking and experimental results implies that docking studies may be a suitable method for predicting and confirming experimental results, as shown in this study. Hence, the combination of molecular docking and spectroscopy methods is an effective innovative approach for binding studies, particularly for pharmacophores.

Transient helical structure during PI3K and Fyn SH3 domain folding

A growing list of proteins, including the β-sheet-rich SH3 domain, is known to transiently populate a compact α-helical intermediate before settling into the native structure. Examples have been discovered in cryogenic solvent as well as by pressure jumps. Earlier studies of λ repressor mutants showed that transient states with excess helix are robust in an all-α protein. Here we extend a previous study of src SH3 domain to two new SH3 sequences, phosphatidylinositol 3-kinase (PI3K) and a Fyn mutant, to see how robust such helix-rich transients are to sequence variations in this β-sheet fold. We quantify helical structure by circular dichroism (CD), protein compactness by small-angle X-ray scattering (SAXS), and transient helical populations by cryo-stopped-flow CD. Our results show that transient compact helix-rich intermediates are easily accessible on the folding landscape of different SH3 domains. In molecular dynamics simulations, force field errors are often blamed for transient non-native structure. We suggest that experimental examples of very fast α-rich transient misfolding could become a more subtle test for further force field improvements than observation of the native state alone.

A circumventing role for the non-native intermediate in the folding of β-lactoglobulin

Folding experiments have suggested that some proteins have kinetic intermediates with a non-native structure. A simple G ̅o model does not explain such non-native intermediates. Therefore, the folding energy landscape of proteins with non-native intermediates should have characteristic properties. To identify such properties, we investigated the folding of bovine β-lactoglobulin (βLG). This protein has an intermediate with a non-native α-helical structure, although its native form is predominantly composed of β-structure. In this study, we prepared mutants whose α-helical and β-sheet propensities are modified and observed their folding using a stopped-flow circular dichroism apparatus. One interesting finding was that E44L, whose β-sheet propensity was increased, showed a folding intermediate with an amount of β-structure similar to that of the wild type, though its folding took longer. Thus, the intermediate seems to be a trapped intermediate. The high α-helical propensity of the wild-type sequence likely causes the folding pathway to circumvent such time-consuming intermediates. We propose that the role of the non-native intermediate is to control the pathway at the beginning of the folding reaction.

Biosynthesis, purification, and characterization of a cannabinoid receptor 2 fragment (CB2(271-326))

Obtaining sufficient amount of purified G-protein coupled receptors (GPCRs) is almost always one of the major challenges for their structural studies. CB2(271-326), a human cannabinoid receptor 2 (CB2) fragment comprising part of the third extracellular loop (EL3), the seventh transmembrane domain (TM7) and C-terminal juxtamembrane region of the receptor, was over-expressed as a fusion protein into inclusion body (IB) of Escherichia coli. The fusion protein was purified by histidine-selected nickel affinity chromatography under denaturing conditions. Then, the fusion protein IBs were solubilized in detergent (Brij58) and the expression fusion leader sequence (TrpLE) was specifically cleaved with tobacco etch virus (TEV) protease. The target fragment, CB2(271-326), was subsequently purified by reverse-phase HPLC and confirmed by SDS-PAGE and mass spectrometry. This hydrophobic fragment can refold in mild detergents digitonin and Brij58. Circular dichroism (CD) spectroscopy of CB2(271-326) in digitonin and Brij58 micelles showed that the fragment adopts a more than 75% alpha-helical structure, with the remainder having beta-strand structure. Fluorescence spectroscopy and quenching studies suggested that the C-terminal region lies near the surface of the digitonin micelles and the TM7 region is folded relatively close to the center of the micelles. This study may provide an alternative strategy for the production and structure/functional studies of GPCRs such as CB2 receptor protein produced in the form of IBs.

Recognition of conformational changes in beta-lactoglobulin by molecularly imprinted thin films

Pathogenesis in protein conformational diseases is initiated by changes in protein secondary structure. This molecular restructuring presents an opportunity for novel shape-based detection approaches, as protein molecular weight and chemistry are otherwise unaltered. Here we apply molecular imprinting to discriminate between distinct conformations of the model protein beta-lactoglobulin (BLG). Thermal- and fluoro-alcohol-induced BLG isoforms were imprinted in thin films of 3-aminophenylboronic acid on quartz crystal microbalance chips. Enhanced rebinding of the template isoform was observed in all cases when compared to the binding of nontemplate isoforms over the concentration range of 1-100 microg mL(-1). Furthermore, it was observed that the greater the changes in the secondary structure of the template protein the lower the binding of native BLG challenges to the imprint, suggesting a strong steric influence in the recognition system. This feasibility study is a first demonstration of molecular imprints for recognition of distinct conformations of the same protein.

Conformation of a double-membrane-spanning fragment of a G protein-coupled receptor: Effects of hydrophobic environment and pH

Overcoming the problems associated with the expression, purification and in vitro handling of membrane proteins requires an understanding of the factors governing the folding and stability of such proteins in detergent solutions. As a sequel to our earlier report (Biochim. Biophys. Acta 1747(2005), 133-140), we describe an improved purification procedure and a detailed structural analysis of a fragment of the mu-opioid receptor ('TM2-3') that comprises the second and third transmembrane segments and the extracellular loop that connects them. Circular dichroism (CD) spectroscopy of TM2-3 in 2,2,2-trifluoroethanol gave a helical content similar to that predicted by published homology models, while spectra acquired in several detergents showed significantly lower helical contents. This indicates that this part of the mu-opioid receptor has an intrinsic propensity to be highly helical in membrane-like environments, but that in detergent solutions, this helical structure is not fully formed. Proteolysis of TM2-3 with trypsin showed that the helical portions of TM2 and TM3 are both shorter than their predicted lengths, indicating that helix-helix interactions in the full-length receptor are apparently important for stabilizing their conformation. Lengthening the alkyl chain of the detergent led to a small but significant increase in the helicity of TM2-3, suggesting that hydrophobic mismatch could play an important role in the stabilization of transmembrane helices by detergents. Protonation of aspartic acid residues in detergent-solubilized TM2-3 also caused a significant increase in helicity. Our results thus suggest that detergent alkyl chain-length and pH may influence membrane protein stability by modulating the stability of individual transmembrane segments.

The sequence TGAAKAVALVL from glyceraldehyde-3-phosphate dehydrogenase displays structural ambivalence and interconverts between α-helical and β-hairpin conformations mediated by collapsed conformational states

The peptide TGAAKAVALVL from glyceraldehyde-3-phosphate dehydrogenase adopts a helical conformation in the crystal structure and is a site for two hydrated helical segments, which are thought to be helical folding intermediates. Overlapping sequences of four to five residues from the peptide, sample both helical and strand conformations in known protein structures, which are dissimilar to glyceraldehyde-3-phosphate dehydrogenase suggesting that the peptide may have a structural ambivalence. Molecular dynamics simulations of the peptide sequence performed for a total simulation time of 1.2 micros, starting from the various initial conformations using GROMOS96 force field under NVT conditions, show that the peptide samples a large number of conformational forms with transitions from alpha-helix to beta-hairpin and vice versa. The peptide, therefore, displays a structural ambivalence. The mechanism from alpha-helix to beta-hairpin transition and vice versa reveals that the compact bends and turns conformational forms mediate such conformational transitions. These compact structures including helices and hairpins have similar hydrophobic radius of gyration (Rgh) values suggesting that similar hydrophobic interactions govern these conformational forms. The distribution of conformational energies is Gaussian with helix sampling lowest energy followed by the hairpins and coil. The lowest potential energy of the full helix may enable the peptide to take up helical conformation in the crystal structure of the glyceraldehyde-3-phosphate dehydrogenase, even though the peptide has a preference for hairpin too. The relevance of folding and unfolding events observed in our simulations to hydrophobic collapse model of protein folding are discussed.

Comparison of the conformational stability of the non-native α-helical intermediate of thiol-modified β-lactoglobulin upon interaction with sodium n-alkyl sulfates at two different pH

Bovine beta-lactoglobulin assumes a dimeric native conformation at neutral pH, while the conformation at pH 2 is monomeric but still native. beta-lactoglobulin has a free thiol at Cys121, which is buried between the beta-barrel and the C-terminal major or alpha-helix. This thiol group was specifically reacted with DTNB (5,5'-dithiobis(2-nitrobenzoic acid)) at pH 7.5 and 2, producing a modified beta-lactoglobulin containing a mix disulfide bond with 5-thio-2-nitrobenzoic acid (TNB). beta-Lactoglobulin is a predominantly beta-sheet protein, although it has a markedly high intrinsic preference for alpha-helical structure. The formation of non-native alpha-helical intermediate of thiol modified beta-lactoglobulin (TNB-beta-LG) was induced by n-alkyl sulfates including sodium octyl sulfate, SOS; sodium decyl sulfate, SDeS; sodium dodecyl sulfate, SDS; and sodium tetradecyl sulfate, STS at pH 7.5 and 2. The conformation and stability of non-native alpha-helical intermediate (alphaI) state of TNB-beta-LG were studied by circular dichroism (CD), fluorescence and differential scanning calorimetry (DSC) techniques. The effect of n-alkyl sulfates on the structure of alphaI state at both pH was utilized to investigate the contribution of hydrophobic interactions to the stability of alphaI intermediate. The present results suggest that the folding reaction of beta-LG follows a non-hierarchical mechanism and hydrophobic interactions play important roles in stabilizing the native state of beta-LG at pH 2 with more positive charges repulsion than at pH 7.5. Then TNB-beta-LG will become a useful model to analyze the conformation and stability of the intermediate of protein folding.

Alpha-helix formation in melittin and beta-lactoglobulin A induced by fluorinated dialcohols

Extensive study of the effect of fluorinated alcohols on protein conformations, notably the induction of alpha-helix formation, is important because of its wide range of applications. Circular dichroism (CD) was used to show that the enhancement of helix induction in beta-lactoglobulin A and melittin by the fluorinated diols 2,2,3,3-tetrafluoro-1,4-butanediol (TFBD), 2,2,3,3,4,4-hexafluoro-1,6-pentanediol (HFPD), and 2,2,3,3,4,4,5,5-octafluoro-1,6-hexanediol (OFHD) increases in the order TFBD < HFPD < OFHD. For fluorinated diols and monoalcohols the effectiveness of helix induction was found to increase exponentially with increasing number of fluorine atoms per alcohol molecule, and OFHD was found to be more effective than any previously reported fluorinated alcohol. Formation of standard micelles was ruled out as the cause of the enhanced helix induction by the fluorinated diols. The negligible red-edge excitation shift in the fluorescence of melittin indicated that the fluorinated diol/water solvent shell surrounding the tryptophan chromophore is less immobilized than are molecules in a lamellar vesicle.

Cooperative α-helix formation of β-lactoglobulin induced by sodium n-alkyl sulfates

It is generally assumed that folding intermediates contain partially formed native-like secondary structures. However, if we consider the fact that the conformational stability of the intermediate state is simpler than that of the native state, it would be expected that the secondary structures in a folding intermediate would not necessarily be similar to those of the native state. beta-Lactoglobulin is a predominantly beta-sheet protein, although it has a markedly high intrinsic preference for alpha-helical structure. The formation of non-native alpha-helical intermediate of beta-lactoglobulin was induced by n-alkyl sulfates including sodium octyl sulfate, SOS; sodium decyl sulfate, SDeS; sodium dodecyl sulfate, SDS; and sodium tetradecyl sulfate, STS at special condition. The effect of n-alkyl sulfates on the structure of native beta-lactoglobulin at pH 2 was utilized to investigate the contribution of hydrophobic interactions to the stability of non-native alpha-helical intermediate. The addition of various concentrations of n-alkyl sulfates to the native state of beta-lactoglobulin (pH 2) appears to support the stabilized form of non-native alpha-helical intermediate at pH 2. The m values of the intermediate state of beta-lactoglobulin by SOS, SDeS, SDS and STS showed substantial variation. The enhancement of m values as the stability criterion of non-native alpha-helical intermediate state corresponded with increasing chain length of the cited n-alkyl sulfates. The present results suggest that the folding reaction of beta-lactoglobulin follows a non-hierarchical mechanism and hydrophobic interactions play important roles in stabilizing the non-native alpha-helical intermediate state.

Helical and Expanded Conformation of Equine β-Lactoglobulin in the Cold-denatured State

The thermal unfolding transition of equine beta-lactoglobulin (ELG) was investigated by circular dichroism (CD) over a temperature range of -15 degrees C to 85 degrees C. In the presence of 2 M urea, a cooperative unfolding transition was observed both with increasing and decreasing temperature. The CD spectrum indicated that the heat and cold-denatured states of ELG have substantial secondary structures but lack persistent tertiary packing of the side-chains. In order to clarify the relation between the heat or cold-denatured state and the acid-denatured (A) state characterized previously, we have attempted to observe the temperature dependence of the CD spectrum at pH 1.5. The CD spectrum in the heat-denatured state is similar to that in the A state. The CD spectrum in the A state does not change cooperatively with increasing temperature. These results indicate that the heat-denatured state and the A state are the same structural state. On the other hand, the CD intensity at acid pH cooperatively increased with decreasing temperature. The CD spectrum at low temperature and acid pH is consistent with that in the cold-denatured state. Therefore, the cold-denatured state is distinguished from the heat-denatured state or the A state, and ELG assumes a larger amount of non-native alpha-helices in the cold-denatured state. Small angle X-ray scattering and analytical ultracentrifugation have indicated that ELG assumes an expanded chain-like conformation in the cold-denatured state in contrast to the compact globular conformation in the A state. The relation between the molecular size and the helical content in the partially folded states is discussed.

Diffusion of water and sodium counter-ions in nanopores of a β-lactoglobulin crystal: a molecular dynamics study

The dynamics of water and sodium counter-ions (Na(+)) in a C222(1) orthorhombic β-lactoglobulin crystal is investigated by means of 5 ns molecular dynamics simulations. The effect of the fluctuation of the protein atoms on the motion of water and sodium ions is studied by comparing simulations in a rigid and in a flexible lattice. The electrostatic interactions of sodium ions with the positively charged LYS residues inside the crystal channels significantly influence the ionic motion. According to our results, water molecules close to the protein surface undergo an anomalous diffusive motion. On the other hand, the motion of water molecules further away from the protein surface is normal diffusive. Protein fluctuations affect the diffusion constant of water, which increases from 0.646 ± 0.108 to 0.887 ± 0.41 nm(2) ns(-1), when protein fluctuations are taken into account. The pore size (0.63-1.05 nm) and the water diffusivities are in good agreement with previous experimental results. The dynamics of sodium ions is disordered. LYS residues inside the pore are the main obstacles to the motion of sodium ions. However, the simulation time is still too short for providing a precise description of anomalous diffusion of sodium ions. The results are not only of interest for studying ion and water transport through biological nanopores, but may also elucidate water-protein and ion-protein interactions in protein crystals.

Cooperative effect of hydrophobic and electrostatic forces on alcohol-induced α-helix formation of α1-acid glycoprotein

Alpha1-acid glycoprotein (AGP) is a serum glycoprotein that mainly binds basic drugs. Previous reports have shown that AGP converts from a beta-sheet to an alpha-helix upon interaction with biomembranes. In the current studies, we found that alkanols, diols, and halogenols all induce this conformational change. Increased length and bulkiness of the hydrocarbon group and the presence of a halogen atom promoted this conversion, whereas the presence of a hydroxyl group inhibited it. Moreover, the effect was dependent on the hydrophobic and electrostatic properties of the alcohols. These results indicate that, in a membrane environment, hydrophobic and electrostatic factors cooperatively induce the transition of AGP from a beta-sheet to an alpha-helix.

Expression and spectroscopic characterization of a large fragment of the μ-opioid receptor

We report here a procedure for the production in Escherichia coli and subsequent purification and characterization of an 80-residue fragment of the human mu-opioid receptor. The fragment ('TM2-3'), which comprises the second and third transmembrane segments as well as the first extracellular loop of the receptor, was expressed as a fusion with glutathione-S-transferase. The fusion protein, which accumulated in insoluble inclusion bodies, was solubilized with N-lauroylsarcosine, and TM2-3 was obtained by thrombin cleavage of the fusion protein followed by reversed-phase HPLC purification. CD spectroscopy of TM2-3 in lysophosphatidylcholine micelles showed that TM2-3 adopts approximately 50% alpha-helical structure in this environment, with the remainder consisting of disordered and/or beta-structure. This is consistent with the assumption of an alpha-helical structure by the two membrane-spanning regions and a nonhelical structure in the loop region of TM2-3. Fluorescence spectroscopy and fluorescence quenching experiments suggested that the extracellular loop lies near the surface of the lysophosphatidylcholine micelle. Our work shows that the study of large receptor fragments is a technically accessible approach to the study of the structural properties of the mu-opioid receptor and, possibly, other G-protein-coupled receptors as well.

Conformational changes of beta-lactoglobulin in sodium bis(2-ethylhexyl) sulfosuccinate reverse micelles. A fluorescence and CD study

The effect of beta-lactoglobulin encapsulation in sodium bis(2-ethylhexyl) sulfosuccinate reverse micelles on the environment of protein and on Trp was analysed at different water contents (omega0). CD data underlined the distortion of the beta-sheet and a less constrained tertiary structure as the omega0 increased, in agreement with a concomitant red shift and a decrease in the signal intensity obtained in steady-state fluorescence measurements. Fluorescence lifetimes, evaluated by biexponential analysis, were tau1 = 1.28 ns and tau2 = 3.36 ns in neutral water. In reverse micelles, decay-associated spectra indicated the occurrence of important environmental changes associated with omega0. Bimolecular fluorescence quenching by CCl4 and acrylamide was employed to analyse alterations in the accessibility of the two Trp residues in beta-lactoglobulin, induced by changes in omega0. The average bimolecular quenching constant <kq(CCl4)> was found not to depend on omega0, confirming the insolubility of this quencher in the aqueous interface, while <kq(arcylamide)> increases with omega0. The drastic decrease with omega0 of kq, associated with the longest lifetime kq2(CCl4), comparatively to the increase of kq2(acrylamide), emphasizes the location of beta-lactoglobulin in the aqueous interfacial region especially at omega0> or = 10. The fact that (omega0 = 30) >> kq2(acrylamide) (water) also confirms the important conformational changes of encapsulated beta-lactoglobulin.

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.

Coupled prediction of protein secondary and tertiary structure

The strong coupling between secondary and tertiary structure formation in protein folding is neglected in most structure prediction methods. In this work we investigate the extent to which nonlocal interactions in predicted tertiary structures can be used to improve secondary structure prediction. The architecture of a neural network for secondary structure prediction that utilizes multiple sequence alignments was extended to accept low-resolution nonlocal tertiary structure information as an additional input. By using this modified network, together with tertiary structure information from native structures, the Q3-prediction accuracy is increased by 7-10% on average and by up to 35% in individual cases for independent test data. By using tertiary structure information from models generated with the ROSETTA de novo tertiary structure prediction method, the Q3-prediction accuracy is improved by 4-5% on average for small and medium-sized single-domain proteins. Analysis of proteins with particularly large improvements in secondary structure prediction using tertiary structure information provides insight into the feedback from tertiary to secondary structure.

Mechanism of fast protein folding

An explosion of in vitro experimental data on the folding of proteins has revealed many examples of folding in the millisecond or faster timescale, often occurring in the absence of stable intermediate states. We review experimental methods for measuring fast protein folding kinetics, and then discuss various analytical models used to interpret these data. Finally, we classify general mechanisms that have been proposed to explain fast protein folding into two catagories, heterogeneous and homogeneous, reflecting the nature of the transition state. One heterogeneous mechanism, the diffusion-collision mechanism, can be used to interpret experimental data for a number of proteins.

Refolding of beta-lactoglobulin studied by stopped-flow circular dichroism at subzero temperatures

Refolding of bovine beta-lactoglobulin was studied by stopped-flow circular dichroism at subzero temperatures. In ethylene glycol 45%-buffer 55% at -15 degrees C, the isomerization rate from the kinetic intermediate rich in alpha-helix to the native state is approximately 300-fold slower than that at 4 degrees C in the absence of ethylene glycol, whereas the initial folding is completed within the dead time of the stopped-flow apparatus (10 ms). At -28 degrees C, we observed at least three phases; the fastest process, accompanied by an increase of alpha-helix content, is completed within the dead time of the stopped-flow apparatus (10 ms), the second phase, accompanied by an increase of alpha-helix content with the rate of 2 s(-1), and the third phase, accompanied by a decrease of alpha-helix content. This last phase, corresponding to the isomerization process at -15 degrees C described above, was so slow that we could not monitor any changes within 4 h. Based on the findings above, we propose that rapid alpha-helix formation and their concurrent collapse are common even in proteins rich in beta-structure in their native forms.

Salt-dependent monomer-dimer equilibrium of bovine beta-lactoglobulin at pH 3

Although bovine beta-lactoglobulin assumes a monomeric native structure at pH 3 in the absence of salt, the addition of salts stabilizes the dimer. Thermodynamics of the monomer-dimer equilibrium dependent on the salt concentration were studied by sedimentation equilibrium. The addition of NaCl, KCl, or guanidine hydrochloride below 1 M stabilized the dimer in a similar manner. On the other hand, NaClO(4) was more effective than other salts by about 20-fold, suggesting that anion binding is responsible for the salt-induced dimer formation, as observed for acid-unfolded proteins. The addition of guanidine hydrochloride at 5 M dissociated the dimer into monomers because of the denaturation of protein structure. In the presence of either NaCl or NaClO(4), the dimerization constant decreased with an increase in temperature, indicating that the enthalpy change (DeltaH(D)) of dimer formation is negative. The heat effect of the dimer formation was directly measured with an isothermal titration calorimeter by titrating the monomeric beta-lactoglobulin at pH 3.0 with NaClO(4). The net heat effects after subtraction of the heat of salt dilution, corresponding to DeltaH(D), were negative, and were consistent with those obtained by the sedimentation equilibrium. From the dependence of dimerization constant on temperature measured by sedimentation equilibrium, we estimated the DeltaH(D) value at 20 degrees C and the heat capacity change (DeltaC(p)) of dimer formation. In both NaCl and NaClO(4), the obtained DeltaC(p) value was negative, indicating the dominant role of burial of the hydrophobic surfaces upon dimer formation. The observed DeltaC(p) values were consistent with the calculated value from the X-ray dimeric structure using a method of accessible surface area. These results indicated that monomer-dimer equilibrium of beta-lactoglobulin at pH 3 is determined by a subtle balance of hydrophobic and electrostatic effects, which are modulated by the addition of salts or by changes in temperature.

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.

Initial denaturing conditions influence the slow folding phase of acylphosphatase associated with proline isomerization

The folding kinetics of human common-type acylphosphatase (cAcP) from its urea- and TFE-denatured states have been determined by stopped-flow fluorescence techniques. The refolding reaction from the highly unfolded state formed in urea is characterized by double exponential behavior that includes a slow phase associated with isomerism of the Gly53-Pro54 peptide bond. However, this slow phase is absent when refolding is initiated by dilution of the highly a-helical denatured state formed in the presence of 40% trifluoroethanol (TFE). NMR studies of a peptide fragment corresponding to residues Gly53-Gly69 of cAcP indicate that only the native-like trans isomer of the Gly-Pro peptide bond is significantly populated in the presence of TFE, whereas both the cis and trans isomers are found in an approximately 1:9 ratio for the peptide bond in aqueous solution. Molecular modeling studies in conjunction with NMR experiments suggest that the trans isomer of the Gly53-Pro54 peptide bond is stabilized in TFE by the formation of a nonnative-like hydrogen bond between the CO group of Gly53 and the NH group of Lys57. These results therefore reveal that a specific nonnative interaction in the denatured state can increase significantly the overall efficiency of refolding.

Structural cassette mutagenesis in a de novo designed protein: proof of a novel concept for examining protein folding and stability

The solution to the protein folding problem lies in defining the relative energetic contributions of short-range and long-range interactions. In other words, the tendency of a stretch of amino acids to adopt a final secondary structural fold is context dependent. Our approach to this problem is to address whether an amino acid sequence, a "cassette," with a defined secondary structure in the three-dimensional structure of a native protein, can adopt a different conformation when placed into a different protein environment. Thus, we designed de novo a disulfide-bridged two-stranded alpha-helical parallel coiled coil, where each polypeptide chain consisted of 39 residues, as a "cassette holder." The 11-residue cassette would be inserted into the center of each polypeptide chain between the two nucleating alpha-helices to replace the control sequence. This Structural Cassette Mutagenesis model permits the analysis of short-range interactions within the inserted cassette as well as long-range interactions between the nucleating helices and the cassette region. The cassette holder, with a control sequence as the cassette, had a GdnHCl transition midpoint during denaturation of 5.6M. To demonstrate the feasibility of our model, an 11-residue beta-strand cassette from an immunoglobulin fold was inserted. The cassette was fully induced into the alpha-helical conformation with a [GdnHCl]1/2 value of 3.2M. To demonstrate the importance of short-range interactions (beta-sheet/alpha-helical propensities of amino acid side chains) in modulating structure and stability, a series of 1-5 threonine residues (highest beta-sheet propensity) were substituted into the solvent-exposed portions of the cassette in the alpha-helical conformation. Each successive substitution systematically decreased the stability of the coiled coil with peptide T4b (4 Thr residues) having a [GdnHCl]1/2 value of 2.2M. The single substitution of Ile in the hydrophobic core of the cassette with Ala or Thr had the most dramatic effect on protein stability (peptide 120T, [GdnHCl]1/2 value of 1.4M). Though these substitutions were able to modulate stability, they were not able to disrupt the alpha-helical conformation of the cassette, showing the importance of the nucleating alpha-helices on either side of the cassette in controlling conformation of the cassette. We have demonstrated the feasibility of our model protein to accept a beta-strand cassette. The effect of cassettes containing other beta-strands, beta-turns, loops, regions of undefined structure, and helical segments on conformation and stability of our model protein will also be determined.

Design of peptides undergoing self-catalytic alpha-to-beta transition and amyloidogenesis

Improved understanding of amyloidogenic peptides and proteins such as prion proteins and Alzheimer's beta peptides has attracted much attention to the elucidation of the molecular mechanisms of such amyloidogenesis. As a representative, in the prion protein, the conformational transitions from alpha-helix to beta-structure undergo along with the amyloidogenesis in a self-catalytic manner. Moreover, recent studies by the de novo design of peptides and proteins as well as the amyloidogenesis of peptides and proteins including pathogenic protein mutants have provided insight into the conformational changes essential to amyloidogenesis and correct folding.

Molten globule structure of equine beta-lactoglobulin probed by hydrogen exchange

The molten globule structure of equine beta-lactoglobulin has been inferred from the hydrogen exchange protection of the backbone amide protons. In order to make it possible to measure the hydrogen exchange kinetics of the individual backbone amide protons, the uniformly (15)N-labeled recombinant protein was expressed in Escherichia coli and the NMR peak assignment was obtained for most of the backbone protons. The chemical shift and NOE results obtained under the condition where the protein assumes the native structure are fully consistent with the known secondary structure of bovine beta-lactoglobulin, indicating that the equine protein has a similar native conformation to that of the bovine protein. The hydrogen exchange rate of the individual backbone amide protons was measured under the conditions where the protein assumes the native and molten globule states. In the native state, strong protection was observed for the residues located in the eight (A to H) strands, which form a barrel structure, and residues of a major helix. In the molten globule state at acidic pH conditions, significant protection from the exchange has been observed for residues located in the A, F, G and H strands in the native structure. The pattern of protection is consistent with a native-like beta-sheet formation by these strands. The residues located in a major helix of the native structure are also protected, suggesting that this helix is formed in the molten globule and is packed against the sheet as in the native structure. These results indicate that a native-like subdomain is formed in the molten globule state of equine beta-lactoglobulin.

Identification of peptides in aggregates formed during hydrolysis of beta-lactoglobulin B with a Glu and Asp specific microbial protease

The purpose of the present study was to identify the peptides responsible for aggregate formation during hydrolysis of beta-lactoglobulin by BLP at neutral pH. Hydrolysates taken at various stages of aggregate formation were separated into a precipitate and a soluble phase and each was analyzed by CE and mass spectrometry. The aggregates consisted of six to seven major peptides of which four were tentatively identified. The peptides were positively charged at neutral pH and had a high charge-to-mass ratio at low pH. The fragment f135-158 seemed to be the initiator of aggregation, since it was present at high concentration in the aggregates at all stages, and the concentration of this peptide remained low in the supernatant. F135-158 contains several basic and acid amino acids alternating with hydrophobic amino acids, which is in accordance with formation of noncovalently linked aggregates, as previously shown.

Is folding of beta-lactoglobulin non-hierarchic? Intermediate with native-like beta-sheet and non-native alpha-helix

The refolding of beta-lactoglobulin, a beta-barrel protein consisting of beta strands betaA-betaI and one major helix, is unusual because non-native alpha-helices are formed at the beginning of the process. We studied the refolding kinetics of bovine beta-lactoglobulin A at pH 3 using the stopped-flow circular dichroism and manual H/(2)H exchange pulse labeling coupled with heteronuclear NMR. The protection pattern from the H/(2)H exchange of the native state indicated the presence of a stable hydrophobic core consisting of betaF, betaG and betaH strands. The protection pattern of the kinetic intermediate obtained about one second after initiating the reaction was compared with that of the native state. In this relatively late kinetic intermediate, which still contains some non-native helical structure, the disulfide-bonded beta-hairpin made up of betaG and betaH strands was formed, but the rest of the molecule was fluctuating, where the non-native alpha-helices may reside. Subsequently, the core beta-sheet extends, accompanied by a further alpha-helix to beta-sheet transition. Thus, the refolding of beta-lactoglobulin exhibits two elements: the critical role of the core beta-sheet is consistent with the hierarchic mechanism, whereas the alpha-helix to beta-sheet transition suggests the non-hierarchic mechanism.

New insight into the pH-dependent conformational changes in bovine beta-lactoglobulin from Raman optical activity

We have studied the conformation of beta-lactoglobulin in aqueous solution at room temperature over the pH range approximately 2.0-9.0 using vibrational Raman optical activity (ROA). The ROA spectra clearly show that the basic up and down beta-barrel core is preserved over the entire pH range, in agreement with other studies. However, from the shift of a sharp positive ROA band at approximately 1268 to approximately 1294 cm(-1) on going from pH values below that of the Tanford transition, which is centered at pH approximately 7.5, to values above, the Tanford transition appears to be associated with changes in the local conformations of residues in loop sequences possibly corresponding to a migration into the alpha-helical region of the Ramachandran surface from a nearby region. These changes may be related to those detected in X-ray crystal structures which revealed that the Tanford transition is associated with conformational changes in loops which form a doorway to the interior of the protein. The results illustrate how the ability of ROA to detect loop and turn structure separately from secondary structure is useful for studying conformational plasticity in proteins.

Polymer principles and protein folding

This paper surveys the emerging role of statistical mechanics and polymer theory in protein folding. In the polymer perspective, the folding code is more a solvation code than a code of local phipsi propensities. The polymer perspective resolves two classic puzzles: (1) the Blind Watchmaker's Paradox that biological proteins could not have originated from random sequences, and (2) Levinthal's Paradox that the folded state of a protein cannot be found by random search. Both paradoxes are traditionally framed in terms of random unguided searches through vast spaces, and vastness is equated with impossibility. But both processes are partly guided. The searches are more akin to balls rolling down funnels than balls rolling aimlessly on flat surfaces. In both cases, the vastness of the search is largely irrelevant to the search time and success. These ideas are captured by energy and fitness landscapes. Energy landscapes give a language for bridging between microscopics and macroscopics, for relating folding kinetics to equilibrium fluctuations, and for developing new and faster computational search strategies.

Equilibrium unfolding CD studies of bovine beta-lactoglobulin and its 14-52 fragment at acidic pH

Bovine beta-lactoglobulin represents an interesting example of context-dependent secondary structure induction. In fact, secondary structure predictions indicated that this beta-barrel protein has a surprisingly high alpha-helical preference, which was retained for short fragments. Cooperative transitions from the native beta-sheet to alpha-helical structures were additionally induced by organic solvents, in particular trifluoroethanol. As a result of this high alpha-helical preference, it has been proposed that non-native alpha-helical intermediates could be formed in the unfolding pathway of this protein. In order to provide a better understanding of the processes that underlie conformational plasticity in this protein, CD measurements in the presence of increasing amounts of urea and in the presence of organic solvents were performed. Urea unfolding studies, performed at pH 2.1 and 37 degrees C, revealed an apparent two-state transition, and afforded no evidence of non native alpha-helical intermediates. The protein treated with up to 6M urea, refolded to the native structure, while treatment with higher molar concentration urea, lead to partial misfolding. A 29-mer peptide covering the region of strands a and b of the intact protein, characterized by the presence of 4/3 heptad repeats, was synthesized and studied by CD in the presence of different solvents. On the basis of the obtained results, a mechanism was proposed to explain the structural transition from the beta to alpha structure, provoked by organic solvents in the intact protein.

Optimization of hydrophobic domains in peptides that undergo transformation from alpha-helix to beta-fibril

Recent studies on peptide fibrillogenesis by the de novo method as well as amyloidogenic proteins including prion proteins and Alzheimer's beta-peptides have provided insights into the conformational changes, such as alpha-helix to beta-structure, involved in folding and misfolding processes. We have found that an exposed hydrophobic nucleation domain at N-terminal causes a structural transition of a peptide from alpha-helix to beta-fibril. It became clear that N-terminal acyl groups of particular lengths in a 2alpha-helix peptide caused the peptide to undergo an alpha-to-beta transition. The peptide with the octanoyl group (C8-2alpha) showed the highest rate of transformation. The study of the designed peptides revealed that these alpha-to-beta transitions were closely related to the initial alpha-helix conformation and its stability. Engineering peptides that undergo alpha-to-beta transitions are attractive not only to the study of pathogenic proteins such as prion proteins, but also to the control of self-assembly of peptides, which will lead to the development of peptidyl self-assembling materials.

Prion protein fragments spanning helix 1 and both strands of beta sheet (residues 125-170) show evidence for predominantly helical propensity by CD and NMR

BACKGROUND

Transmissible spongiform encephalopathies are a group of neurodegenerative disorders of man and animals that are believed to be caused by an alpha-helical to beta-sheet conformational change in the prion protein, PrP. Recently determined NMR structures of recombinant PrP (residues 121-231 and 90-231) have identified a short two-stranded anti-parallel beta sheet in the normal cellular form of the protein (PrPC). This beta sheet has been suggested to be involved in seeding the conformational transition to the disease-associated form (PrPSc) via a partially unfolded intermediate state.

RESULTS

We describe CD and NMR studies of three peptides (125-170, 142-170 and 156-170) that span the beta-sheet and helix 1 region of PrP, forming a large part of the putative PrPSc-PrPC binding site that has been proposed to be important for self-seeding replication of PrPSc. The data suggest that all three peptides in water have predominantly helical propensities, which are enhanced in aqueous methanol (as judged by deviations from random-coil Halpha chemical shifts and 3JHalpha-NH values). Although the helical propensity is most marked in the region corresponding to helix 1 (144-154), it is also apparent for residues spanning the two beta-strand sequences.

CONCLUSIONS

We have attempted to model the conformational properties of a partially unfolded state of PrP using peptide fragments spanning the region 125-170. We find no evidence in the sequence for any intrinsic conformational preference for the formation of extended beta-like structure that might be involved in promoting the PrPC-PrPSc conformational transition.