Temperature dependence of the thermodynamics of helix-coil transition

Mechanisms of helix induction by the closed loop

Although various factors affecting the α-helix stability have been known, they are limited to short-range interactions that occur within the helix and its flanking regions. For better understanding of protein folding coupled with disulfide bond formation, the effect of disulfide bonds on the α-helix stability must be clarified. Using a protein fragment in which two helical regions are linked by a disulfide bond, we investigated the influence of the number of residues within the loop closed by the disulfide bond on the helix stability. We modified the number of residues within the loop by inserting glycine residues in the nonhelical region of the protein fragment. Circular dichroism and nuclear magnetic resonance spectra showed that increasing the number of inserted glycine residues led to a decreased helix content, while the helical regions remained unchanged. The helical fractions of individual residues were derived from chemical shift values, and their dependences on the number of inserted glycine residues were investigated. The results of this study and previous studies support the hypothesis that the helices are nucleated in the loop and propagated to both N- and C-termini. In addition, the fact that the helical fractions decrease with the number of inserted glycine residues suggests the mechanism in which the loop reduces the number of possible conformations, thereby promoting helix nucleation within the loop.

A new design and computational survey on RGD biofunctionalized RADA16-I self-assembling peptide for tissue engineering applications

Tissue engineering constantly needs innovative and biocompatible materials, and peptide-based materials seem very inspiring. Here we developed two new self-assembling peptides based on RADA16-I and RGD peptides and studied their potential in forming nanofibers under various conditions using all-atom and coarse-grained molecular dynamics simulation methods. First, a double-tailed RGD (dtRGD) peptide was designed by attaching two RADA16-I tails to an RGD-containing loop in which two disulfide bonds stabilized the loop integrity. In the second design, we bonded one side of the loop to the DA16-I tail (otRGD). The dtRGD peptides exhibited a remarkable propensity to form beta-sheet structures during all-atom MD simulations, starting from the initial random coil structure. The most promising outcomes in nanofiber formation were observed when simulating these peptides in a salt concentration that mimics the extracellular matrix. The representation of the RGD epitope was also significantly evident under these conditions. In the otRGD design, the final structure displayed a globular-like morphology, predominantly possessing coils and alpha-helices secondary structures, while maintaining effective RGD peptide exposure. This investigation signified the possibility of a new RGD representing biomaterial for tissue engineering purposes, however, further theoretical and experimental investigations are imperative to unlock their capabilities and applications.

Estimation of Peptide Helicity from Circular Dichroism Using the Ensemble Model

An established method for the quantitation of the helix content in peptides using circular dichroism (CD) relies on the linear spectroscopic model. This model assumes an average value of the helix-length correction for all peptide conformers, irrespective of the length of the helical segment. Here we assess the validity of this approximation and introduce a more physically realistic ensemble-based analysis of the CD signal in which the length correction is assigned specifically to each ensemble conformer. We demonstrate that the linear model underestimates peptide helicity, with the difference depending on the ensemble composition. We developed a computer program that implements the ensemble model to estimate the peptide helicity. Using this model and the CD data set covering a broad range of helicities, we recalibrate CD baseline parameters and redetermine helix-coil parameters for the alanine-rich peptide. We show that the ensemble model leverages small differences in signal between conformers to extract more information from the experimental data, enabling the determination of several poorly defined quantities, such as the nucleation constant and heat capacity change associated with helix folding. Overall, the presented ensemble-based treatment of the CD signal, together with the recalibrated values of the spectroscopic baseline parameters, provides a coherent framework for the analysis of the peptide helix content.

Development of a new heparan sulfate proteoglycan (HSPG) chromolith LC column to study the pH dependence binding of peptide vaccines to HSPG and role of human serum albumin on its binding

Heparan sulfate proteoglycan (HSPG) expressed on immune cell surface participate in antitumor T-cell responses generated in the acidic lymph node (LN) microenvironment. In this work, HSPG was immobilized for the first time on a HPLC chromolith support for studying the effect of extra cellular acidosis in LNs on the binding to HSPG of two peptide vaccines (universal cancer peptide UCP2 and UCP4). This home-made HSPG column enabling to work at high flow-rates, was resistance to change in pH, had a long - life time, an excellent repeatability and negligible non-specific binding sites. The performance of this affinity HSPG column was confirmed by the evaluation of recognition assay for a series of known ligand of HSPG. It was shown that at 37 °C, the UCP2 binding to HSPG versus pH described a sigmoidal shape while UCP4 remained relatively constant in the pH range 5.0-7.5 and lower than the one of UCP2. By the use of an HSA HPLC column, it was shown at 37 °C and in acidic conditions a loss of affinity of UCP2 and UCP4 to HSA. It was demonstrated that upon UCP2/HSA binding, the protonation of the histidine residue in the cluster R(arg) Q(Gln) Hist (H) of the UCP2 peptide allowed to expose more favorably than UCP4 its polar and cationic groups to the negative net charge of HSPG on immune cells. Acidic pHs led to the protonation of the UCP2 residue histidine by flipping the His switch to the on position with a concomitant increase in affinity for the negative net charge of HSPG confirming that UCP2 was more immunogenic than UCP4. As well this HSPG chromolith LC column developed in this work could be used in the feature for other protein - HSPG binding studies or for a separative mode.

Pseudo-Isolated α-Helix Platform for the Recognition of Deep and Narrow Targets

Although interest in stabilized α-helical peptides as next-generation therapeutics for modulating biomolecular interfaces is increasing, peptides have limited functionality and stability due to their small size. In comparison, α-helical ligands based on proteins can make steric clash with targets due to their large size. Here, we report the design of a monomeric pseudo-isolated α-helix (mPIH) system in which proteins behave as if they are peptides. The designed proteins contain α-helix ligands that do not require any covalent chemical modification, do not have frayed ends, and importantly can make sterically favorable interactions similar to isolated peptides. An optimal mPIH showed a more than 100-fold increase in target selectivity, which might be related to the advantages in conformational selection due to the absence of frayed ends. The α-helical ligand in the mPIH displayed high thermal stability well above human body temperature and showed reversible and rapid folding/unfolding transitions. Thus, mPIH can become a promising protein-based platform for developing stabilized α-helix pharmaceuticals.

From Protein Design to the Energy Landscape of a Cold Unfolding Protein

Understanding protein folding is crucial for protein sciences. The conformational spaces and energy landscapes of cold (unfolded) protein states, as well as the associated transitions, are hardly explored. Furthermore, it is not known how structure relates to the cooperativity of cold transitions, if cold and heat unfolded states are thermodynamically similar, and if cold states play important roles for protein function. We created the cold unfolding 4-helix bundle DCUB1 with a de novo designed bipartite hydrophilic/hydrophobic core featuring a hydrogen bond network which extends across the bundle in order to study the relative importance of hydrophobic versus hydrophilic protein-water interactions for cold unfolding. Structural and thermodynamic characterization resulted in the discovery of a complex energy landscape for cold transitions, while the heat unfolded state is a random coil. Below ∼0 °C, the core of DCUB1 disintegrates in a largely cooperative manner, while a near-native helical content is retained. The resulting cold core-unfolded state is compact and features extensive internal dynamics. Below -5 °C, two additional cold transitions are seen, that is, (i) the formation of a water-mediated, compact, and highly dynamic dimer, and (ii) the onset of cold helix unfolding decoupled from cold core unfolding. Our results suggest that cold unfolding is initiated by the intrusion of water into the hydrophilic core network and that cooperativity can be tuned by varying the number of core hydrogen bond networks. Protein design has proven to be invaluable to explore the energy landscapes of cold states and to robustly test related theories.

The crystallization enthalpy and entropy of protein solutions: microcalorimetry, van't Hoff determination and linearized Poisson–Boltzmann model of tetragonal lysozyme crystals

Microcalorimetric and van't Hoff determinations as well as a theoretical description provide a consistent picture of the crystallization enthalpy and entropy of protein solutions and their dependence on physicochemical solution parameters.

Gallic acid influence on bovine serum albumin thermal stability

A thermoanalytical approach reveals the dual action of GA on BSA thermal stability.

Analysis of the motion of vacuolar volutin granules in Saccharomyces cerevisiae

The moving volutin (polyphosphate) granules known as "dancing bodies" can be observed in the vacuoles of the yeast cells. The aim of work was to study the effects of cultivation conditions and influences of physico-chemical factors on the motion of vacuolar volutin granules in Saccharomyces cerevisiae cells. The motion of granules is a non-Markovian process. It does not depend on the cell cycle phase, but depends on the growth stage. The maximal number of cells with "dancing bodies" was observed under cultivation of yeast at 25-28 °C and pH 5.4-5.8. Irradiation by non-ionizing electromagnetic radiation (EMR) of extremely high frequency (61.22 GHz, 100 μW, 30 min) had no effect on granule motion. After irradiation by non-ionizing EMR of very high frequency (40.68 MHz, 30 W, 30 min) the number of cells with "dancing bodies" decreased significantly and in 2 h restored almost to the control value. The possible nature of the moving volutin granules phenomenon due to metabolic processes is discussed.

Influence of Glu/Arg, Asp/Arg, and Glu/Lys Salt Bridges on α-Helical Stability and Folding Kinetics

Using a combination of ultraviolet circular dichroism, temperature-jump transient-infrared spectroscopy, and molecular dynamics simulations, we investigate the effect of salt bridges between different types of charged amino-acid residue pairs on α-helix folding. We determine the stability and the folding and unfolding rates of 12 alanine-based α-helical peptides, each of which has a nearly identical composition containing three pairs of positively and negatively charged residues (either Glu(-)/Arg(+), Asp(-)/Arg(+), or Glu(-)/Lys(+)). Within each set of peptides, the distance and order of the oppositely charged residues in the peptide sequence differ, such that they have different capabilities of forming salt bridges. Our results indicate that stabilizing salt bridges (in which the interacting residues are spaced and ordered such that they favor helix formation) speed up α-helix formation by up to 50% and slow down the unfolding of the α-helix, whereas salt bridges with an unfavorable geometry have the opposite effect. Comparing the peptides with different types of charge pairs, we observe that salt bridges between side chains of Glu(-) and Arg(+) are most favorable for the speed of folding, probably because of the larger conformational space of the salt-bridging Glu(-)/Arg(+) rotamer pairs compared to Asp(-)/Arg(+) and Glu(-)/Lys(+). We speculate that the observed impact of salt bridges on the folding kinetics might explain why some proteins contain salt bridges that do not stabilize the final, folded conformation.

An Ancestral Tryptophanyl-tRNA Synthetase Precursor Achieves High Catalytic Rate Enhancement without Ordered Ground-State Tertiary Structures

Urzymes-short, active core modules derived from enzyme superfamilies-prepared from the two aminoacyl-tRNA synthetase (aaRS) classes contain only the modules shared by all related family members. They have been described as models for ancestral forms. Understanding them currently depends on inferences drawn from the crystal structures of the full-length enzymes. As aaRS Urzymes lack much of the mass of modern aaRS's, retaining only a small portion of the hydrophobic cores of the full-length enzymes, it is desirable to characterize their structures. We report preliminary characterization of (15)N tryptophanyl-tRNA synthetase Urzyme by heteronuclear single quantum coherence (HSQC) NMR spectroscopy supplemented by circular dichroism, thermal melting, and induced fluorescence of bound dye. The limited dispersion of (1)H chemical shifts (0.5 ppm) is inconsistent with a narrow ensemble of well-packed structures in either free or substrate-bound forms, although the number of resonances from the bound state increases, indicating a modest, ligand-dependent gain in structure. Circular dichroism spectroscopy shows the presence of helices and evidence of cold denaturation, and all ligation states induce Sypro Orange fluorescence at ambient temperatures. Although the term "molten globule" is difficult to define precisely, these characteristics are consistent with most such definitions. Active-site titration shows that a majority of molecules retain ∼60% of the transition state stabilization free energy observed in modern synthetases. In contrast to the conventional view that enzymes require stable tertiary structures, we conclude that a highly flexible ground-state ensemble can nevertheless bind tightly to the transition state for amino acid activation.

A comparative spectroscopic and calorimetric investigation of the interaction of amsacrine with heme proteins, hemoglobin and myoglobin

The binding of the anilido aminoacridine derivative amsacrine with the heme proteins, hemoglobin, and myoglobin, was characterized by various spectroscopic and calorimetric methods. The binding affinity to hemoglobin was (1.21 ± .05) × 10⁵ M^-1, while that to myoglobin was three times higher (3.59 ± .15) × 10⁵ M^-1. The temperature-dependent fluorescence study confirmed the formation of ground-state complexes with both the proteins. The stronger binding to myoglobin was confirmed from both spectroscopic and calorimetric studies. The binding was exothermic in both cases at the three temperatures studied, and was favored by both enthalpy and entropy changes. Circular dichroism results, three-dimensional (3D) and synchronous fluorescence studies confirmed that the binding of amsacrine significantly changed the secondary structure of hemoglobin, while the change in the secondary structure of myoglobin was much less. New insights, in terms of structural and energetic aspects of the interaction of amsacrine with the heme proteins, presented here may help in understanding the structure-activity relationship, therapeutic efficacy, and drug design aspects of acridines.

Importance of Hydrophilic Hydration and Intramolecular Interactions in the Thermodynamics of Helix-Coil Transition and Helix-Helix Assembly in a Deca-Alanine Peptide

For a model deca-alanine peptide the cavity (ideal hydrophobic) contribution to hydration favors the helix state over extended states and the paired helix bundle in the assembly of two helices. The energetic contributions of attractive protein-solvent interactions are separated into quasi-chemical components consisting of a short-range part arising from interactions with solvent in the first hydration shell and the remaining long-range part that is well described by a Gaussian. In the helix-coil transition, short-range attractive protein-solvent interactions outweigh hydrophobic hydration and favor the extended coil states. Analysis of enthalpic effects shows that it is the favorable hydration of the peptide backbone that favors the unfolded state. Protein intramolecular interactions favor the helix state and are decisive in favoring folding. In the pairing of two helices, the cavity contribution outweighs the short-range attractive protein-water interactions. However, long-range, protein-solvent attractive interactions can either enhance or reverse this trend depending on the mutual orientation of the helices. In helix-helix assembly, change in enthalpy arising from change in attractive protein-solvent interactions favors disassembly. In helix pairing as well, favorable protein intramolecular interactions are found to be as important as hydration effects. Overall, hydrophilic protein-solvent interactions and protein intramolecular interactions are found to play a significant role in the thermodynamics of folding and assembly in the system studied.

Thermal expansivities of peptides, polypeptides and proteins as measured by pressure perturbation calorimetry

The main goal of this work was to provide direct experimental evidence that the expansivity of peptides, polypeptides and proteins as measured by pressure perturbation calorimetry (PPC), can serve as a proxy to characterize relative compactness of proteins, especially the denatured state ensemble. This is very important as currently only small angle X-ray scattering (SAXS), intrinsic viscosity and, to a lesser degree, fluorescence resonance transfer (FRET) experiments are capable of reporting on the compactness of denatured state ensembles. We combined the expansivity measurements with other biophysical methods (far-UV circular dichroism spectroscopy, differential scanning calorimetry, and small angle X-ray scattering). Three case studies of the effects of conformational changes on the expansivity of polypeptides in solution are presented. We have shown that expansivity appears to be insensitive to the helix-coil transition, and appears to reflect the changes in hydration of the side-chains. We also observed that the expansivity is sensitive to the global conformation of the polypeptide chain and thus can be potentially used to probe hydration of different collapsed states of denatured or even intrinsically disordered proteins.

Thermodynamic characterization of the unfolding of the prion protein

The prion protein appears to be unusually susceptible to conformational change, and unlike nearly all other proteins, it can easily be made to convert to alternative misfolded conformations. To understand the basis of this structural plasticity, a detailed thermodynamic characterization of two variants of the mouse prion protein (moPrP), the full-length moPrP (23-231) and the structured C-terminal domain, moPrP (121-231), has been carried out. All thermodynamic parameters governing unfolding, including the changes in enthalpy, entropy, free energy, and heat capacity, were found to be identical for the two protein variants. The N-terminal domain remains unstructured and does not interact with the C-terminal domain in the full-length protein at pH 4. Moreover, the enthalpy and entropy of unfolding of moPrP (121-231) are similar in magnitude to values reported for other proteins of similar size. However, the protein has an unusually high native-state heat capacity, and consequently, the change in heat capacity upon unfolding is much lower than that expected for a protein of similar size. It appears, therefore, that the native state of the prion protein undergoes substantial fluctuations in enthalpy and hence, in structure.

Efficient Parameter Estimation of Generalizable Coarse-Grained Protein Force Fields Using Contrastive Divergence: A Maximum Likelihood Approach

Maximum Likelihood (ML) optimization schemes are widely used for parameter inference. They maximize the likelihood of some experimentally observed data, with respect to the model parameters iteratively, following the gradient of the logarithm of the likelihood. Here, we employ a ML inference scheme to infer a generalizable, physics-based coarse-grained protein model (which includes Go̅-like biasing terms to stabilize secondary structure elements in room-temperature simulations), using native conformations of a training set of proteins as the observed data. Contrastive divergence, a novel statistical machine learning technique, is used to efficiently approximate the direction of the gradient ascent, which enables the use of a large training set of proteins. Unlike previous work, the generalizability of the protein model allows the folding of peptides and a protein (protein G) which are not part of the training set. We compare the same force field with different van der Waals (vdW) potential forms: a hard cutoff model, and a Lennard-Jones (LJ) potential with vdW parameters inferred or adopted from the CHARMM or AMBER force fields. Simulations of peptides and protein G show that the LJ model with inferred parameters outperforms the hard cutoff potential, which is consistent with previous observations. Simulations using the LJ potential with inferred vdW parameters also outperforms the protein models with adopted vdW parameter values, demonstrating that model parameters generally cannot be used with force fields with different energy functions. The software is available at https://sites.google.com/site/crankite/.

Properties of hydrophobic free energy found by gas-liquid transfer

The hydrophobic free energy in current use is based on transfer of alkane solutes from liquid alkanes to water, and it has been argued recently that these values are incorrect and should be based instead on gas-liquid transfer data. Hydrophobic free energy is measured here by gas-liquid transfer of hydrocarbon gases from vapor to water. The new definition reduces more than twofold the values of the apparent hydrophobic free energy. Nevertheless, the newly defined hydrophobic free energy is still the dominant factor that drives protein folding as judged by ΔCp, the change in heat capacity, found from the free energy change for heat-induced protein unfolding. The ΔCp for protein unfolding agrees with ΔCp values for solvating hydrocarbon gases and disagrees with ΔCp for breaking peptide hydrogen bonds, which has the opposite sign. The ΔCp values for the enthalpy of liquid-liquid and gas-liquid transfer are similar. The plot of free energy against the apparent solvent-exposed surface area is given for linear alkanes, but only for a single conformation, the extended conformation, of these flexible-chain molecules. The ability of the gas-liquid hydrophobic factor to predict protein stability is tested and reasonable agreement is found, using published data for the dependences on temperature of the unfolding enthalpy of ribonuclease T1 and the solvation enthalpies of the nonpolar and polar groups.

Bacterial expression and purification of the amyloidogenic peptide PAPf39 for multidimensional NMR spectroscopy

PAPf39 is a 39 residue peptide fragment from human prostatic acidic phosphatase that forms amyloid fibrils in semen. These fibrils have been implicated in facilitating HIV transmission. To enable structural studies of PAPf39 by NMR spectroscopy, efficient methods allowing the production of milligram quantities of isotopically labeled peptide are essential. Here, we report the high-yield expression and purification of uniformly (13)C- and (15)N-labeled PAPf39 peptide, through expression as a fusion to ubiquitin at the N-terminus and an intein at the C-terminus. This allows the study of the PAPf39 monomer conformational ensemble by NMR spectroscopy. To this end, we performed the NMR chemical shift assignment of the PAPf39 peptide in the monomeric state at low pH.

Natural phenylalanine hydroxylase variants that confer a mild phenotype affect the enzyme's conformational stability and oligomerization equilibrium

Hyperphenylalaninemias are genetic diseases prevalently caused by mutations in the phenylalanine hydroxylase (PAH) gene. The wild-type PAH enzyme is a homotetramer regulated by its substrate, cofactor and phosphorylation. We reproduced a full-length wild-type protein and seven natural full-length PAH variants, p.I65M, p.N223Y, p.R297L, p.F382L, p.K398N, p.A403V, and p.Q419R, and analyzed their biochemical and biophysical behavior. All mutants exhibited reduced enzymatic activity, namely from 38% to 69% of wild-type activity. Biophysical characterization was performed by size-exclusion chromatography, light scattering and circular dichroism. In the purified wild-type PAH, we identified the monomer in equilibrium with the dimer and tetramer. In most mutants, the equilibrium shifted toward the dimer and most tended to form aggregates. All PAH variants displayed different biophysical behaviors due to loss of secondary structure and thermal destabilization. Specifically, p.F382L was highly unstable at physiological temperature. Moreover, using confocal microscopy with the number and brightness technique, we studied the effect of BH4 addition directly in living human cells expressing wild-type PAH or p.A403V, a mild mutant associated with BH4 responsiveness in vivo. Our results demonstrate that BH4 addition promotes re-establishment of the oligomerization equilibrium, thus indicating that the dimer-to-tetramer shift in pA403V plays a key role in BH4 responsiveness. In conclusion, we show that the oligomerization process and conformational stability are altered by mutations that could affect the physiological behavior of the enzyme. This endorses the hypothesis that oligomerization and folding defects of PAH variants are the most common causes of HPAs, particularly as regards mild human phenotypes.

The effect of the placement and total charge of the basic amino acid clusters on antibacterial organism selectivity and potency

Extensive circular dichroism, isothermal titration calorimetry and induced calcein leakage studies were conducted on a series of antimicrobial peptides (AMPs), with a varying number of Lys residues located at either the C-terminus or the N-terminus to gain insight into their effect on the mechanisms of binding with zwitterionic and anionic membrane model systems. Different CD spectra were observed for these AMPs in the presence of zwitterionic DPC and anionic SDS micelles indicating that they adopt different conformations on binding to the surfaces of zwitterionic and anionic membrane models. Different CD spectra were observed for these AMPs in the presence of zwitterionic POPC and anionic mixed 4:1 POPC/POPG LUVs and SUVs, indicating that they adopt very different conformations on interaction with these two types of LUVs and SUVs. In addition, ITC and calcein leakage data indicated that all the AMPs studied interact via very different mechanisms with anionic and zwitterionic LUVs. ITC data suggest these peptides interact primarily with the surface of zwitterionic LUVs while they insert into and form pores in anionic LUVs. CD studies indicated that these compounds adopt different conformations depending on the ratio of POPC to POPG lipids present in the liposome. There are detectable spectroscopic and thermodynamic differences between how each of these AMPs interacts with membranes, that is position and total charge density defines how these AMPs interact with specific membrane models and thus partially explain the resulting diversity of antibacterial activity of these compounds.

Nonadditivity in conformational entropy upon molecular rigidification reveals a universal mechanism affecting folding cooperativity

Previously, we employed a Maxwell counting distance constraint model (McDCM) to describe α-helix formation in polypeptides. Unlike classical helix-coil transition theories, the folding mechanism derives from nonadditivity in conformational entropy caused by rigidification of molecular structure as intramolecular cross-linking interactions form along the backbone. For example, when a hydrogen bond forms within a flexible region, both energy and conformational entropy decrease. However, no conformational entropy is lost when the region is already rigid because atomic motions are not constrained further. Unlike classical zipper models, the same mechanism also describes a coil-to-β-hairpin transition. Special topological features of the helix and hairpin structures allow the McDCM to be solved exactly. Taking full advantage of the fact that Maxwell constraint counting is a mean field approximation applied to the distribution of cross-linking interactions, we present an exact transfer matrix method that does not require any special topological feature. Upon application of the model to proteins, cooperativity within the folding transition is yet again appropriately described. Notwithstanding other contributing factors such as the hydrophobic effect, this simple model identifies a universal mechanism for cooperativity within polypeptide and protein-folding transitions, and it elucidates scaling laws describing hydrogen-bond patterns observed in secondary structure. In particular, the native state should have roughly twice as many constraints as there are degrees of freedom in the coil state to ensure high fidelity in two-state folding cooperativity, which is empirically observed.

Nonadditivity in the alpha-helix to coil transition

The Lifson-Roig Model (LRM) and all its variants describe the α-helix to coil transition in terms of additive component-free energies within a free energy decomposition scheme, and these contributions are interpreted through sequence-context dependent nucleation and propagation parameters. Although this phenomenological approach is able to adequately fit experimental data on helix content and heat capacity, the number of required parameters increases dramatically with additional sequence variation. Moreover, due to nonadditive competing microscopic effects that are difficult to disentangle within a LRM, large uncertainties within the parameters emerge. We offer an alternative view that removes the need for sequence-context parameterization by focusing on individual microsopic interactions within a free energy decomposition and explicitly account for nonadditivity in conformational entropy through network rigidity using a Distance Constraint Model (DCM). We apply a LRM and a DCM to previously published experimental heat capacity and helix content data for a series of heterogeneous polypeptides. Both models describe the experimental data well, and the parameters from both models are consistent with prior work. However, the number of DCM parameters is independent of sequence-variability, the parameter values exhibit better transferability, and the helix nucleation is predicted by the DCM explicitly through the nonadditive nature of conformational entropy. The importance of these results is that the DCM offers a system-independent approach for modeling stability within polypeptides and proteins, where the demonstrated accuracy for the α-helix to coil transition over a series of heterogeneous polypeptides described here is one case in point.

A second look at mini-protein stability: analysis of FSD-1 using circular dichroism, differential scanning calorimetry, and simulations

Mini-proteins that contain <50 amino acids often serve as model systems for studying protein folding because their small size makes long timescale simulations possible. However, not all mini-proteins are created equal. The stability and structure of FSD-1, a 28-residue mini-protein that adopted the betabetaalpha zinc-finger motif independent of zinc binding, was investigated using circular dichroism, differential scanning calorimetry, and replica-exchange molecular dynamics. The broad melting transition of FSD-1, similar to that of a helix-to-coil transition, was observed by using circular dichroism, differential scanning calorimetry, and replica-exchange molecular dynamics. The N-terminal beta-hairpin was found to be flexible. The FSD-1 apparent melting temperature of 41 degrees C may be a reflection of the melting of its alpha-helical segment instead of the entire protein. Thus, despite its attractiveness due to small size and purposefully designed helix, sheet, and turn structures, the status of FSD-1 as a model system for studying protein folding should be reconsidered.

Helix/coil nucleation: a local response to global demands

A complete description of protein structure and function must include a proper treatment of mechanisms that lead to cooperativity. The helix/coil transition serves as a simple example of a cooperative folding process, commonly described by a nucleation-propagation mechanism. The prevalent view is that coil structure must first form a short segment of helix in a localized region despite paying a free energy cost (nucleation). Afterward, helical structure propagates outward from the nucleation site. Both processes entail enthalpy-entropy compensation that derives from the loss in conformational entropy on helix formation with concomitant gain in favorable interactions. Nucleation-propagation models inherently assume that cooperativity arises from a sequential series of local events. An alternative distance constraint model asserts there is a direct link between available degrees of freedom and cooperativity through the nonadditivity in conformational entropy. That is, helix nucleation is a concerted manifestation of rigidity propagating through atomic structure. The link between network rigidity and nonadditivity of conformational entropy is shown in this study by solving the distance constraint model using a simple global constraint counting approximation. Cooperativity arises from competition between excess and deficiency in available degrees of freedom in the coil and helix states respectively.

Mechanism of formation of the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. II. Interplay of local backbone conformational dynamics and long-range hydrophobic interactions in hairpin formation

Two peptides, corresponding to the turn region of the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus, consisting of residues 51-56 [IG(51-56)] and 50-57 [IG(50-57)], respectively, were studied by circular dichroism and NMR spectroscopy at various temperatures and by differential scanning calorimetry. Our results show that the part of the sequence corresponding to the beta-turn in the native structure (DDATKT) of the B3 domain forms bent conformations similar to those observed in the native protein. The formation of a turn is observed for both peptides in a broad range of temperatures (T = 283-323 K), which confirms the conclusion drawn from our previous studies of longer sequences from the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G (16, 14, and 12 residues), that the DDATKT sequence forms a nucleation site for formation of the beta-hairpin structure of peptides corresponding to the C-terminal part of all the B domains of the immunoglobulin binding protein G. We also show and discuss the role of long-range hydrophobic interactions as well as local conformational properties of polypeptide chains in the mechanism of formation of the beta-hairpin structure.

Optimized molecular dynamics force fields applied to the helix-coil transition of polypeptides

Obtaining the correct balance of secondary structure propensities is a central priority in protein force-field development. Given that current force fields differ significantly in their alpha-helical propensities, a correction to match experimental results would be highly desirable. We have determined simple backbone energy corrections for two force fields to reproduce the fraction of helix measured in short peptides at 300 K. As validation, we show that the optimized force fields produce results in excellent agreement with nuclear magnetic resonance experiments for folded proteins and short peptides not used in the optimization. However, despite the agreement at ambient conditions, the dependence of the helix content on temperature is too weak, a problem shared with other force fields. A fit of the Lifson-Roig helix-coil theory shows that both the enthalpy and entropy of helix formation are too small: the helix extension parameter w agrees well with experiment, but its entropic and enthalpic components are both only about half the respective experimental estimates. Our structural and thermodynamic analyses point toward the physical origins of these shortcomings in current force fields, and suggest ways to address them in future force-field development.

Statistical estimation of statistical mechanical models: helix-coil theory and peptide helicity prediction

Analysis of biopolymer sequences and structures generally adopts one of two approaches: use of detailed biophysical theoretical models of the system with experimentally-determined parameters, or largely empirical statistical models obtained by extracting parameters from large datasets. In this work, we demonstrate a merger of these two approaches using Bayesian statistics. We adopt a common biophysical model for local protein folding and peptide configuration, the helix-coil model. The parameters of this model are estimated by statistical fitting to a large dataset, using prior distributions based on experimental data. L(1)-norm shrinkage priors are applied to induce sparsity among the estimated parameters, resulting in a significantly simplified model. Formal statistical procedures for evaluating support in the data for previously proposed model extensions are presented. We demonstrate the advantages of this approach including improved prediction accuracy and quantification of prediction uncertainty, and discuss opportunities for statistical design of experiments. Our approach yields a 39% improvement in mean-squared predictive error over the current best algorithm for this problem. In the process we also provide an efficient recursive algorithm for exact calculation of ensemble helicity including sidechain interactions, and derive an explicit relation between homo- and heteropolymer helix-coil theories and Markov chains and (non-standard) hidden Markov models respectively, which has not appeared in the literature previously.

Miranda cargo-binding domain forms an elongated coiled-coil homodimer in solution: implications for asymmetric cell division in Drosophila

Miranda is a multidomain adaptor protein involved in neuroblast asymmetric division in Drosophila melanogaster. The central domain of Miranda is necessary for cargo binding of the neural transcription factor Prospero, the Prospero-mRNA carrier Staufen, and the tumor suppressor Brat. Here, we report the first solution structure of Miranda central "cargo-binding" domain (residues 460-660) using small-angle X-ray scattering. Ab initio modeling of the scattering data yields an elongated "rod-like" molecule with a maximum linear dimension (D(max)) of approximately 22 nm. Moreover, circular dichroism and cross-linking experiments indicate that the cargo-binding domain is predominantly helical and forms a parallel coiled-coil homodimer in solution. Based on the results, we modeled the full-length Miranda protein as a double-headed, double-tailed homodimer with a long central coiled-coil region. We discuss the cargo-binding capacity of the central domain and propose a structure-based mechanism for cargo release and timely degradation of Miranda in developing neuroblasts.

Structure of the nucleocapsid-binding domain from the mumps virus polymerase; an example of protein folding induced by crystallization

The human pathogen mumps virus, like all paramyxoviruses, encodes a polymerase responsible for virally directed RNA synthesis. The template for the polymerase is the nucleocapsid, a filamentous protein-RNA complex harboring the viral genome. Interaction of the polymerase and the nucleocapsid is mediated by a small domain tethered to the end of the phosphoprotein (P), one of the polymerase subunits. We report the X-ray crystal structure of this region of mumps virus P (the nucleocapsid-binding domain, or NBD, amino acids 343-391). The mumps P NBD forms a compact bundle of three alpha-helices within the crystal, a fold apparently conserved across the Paramyxovirinae. In solution, however, the domain exists in the molten globule state. This is demonstrated through application of differential scanning calorimetry, circular dichroism spectroscopy, NMR spectroscopy, and dynamic light scattering. While the mumps P NBD is compact and has persistent secondary structure, it lacks a well-defined tertiary structure under normal solution conditions. It can, however, be induced to fold by addition of a stabilizing methylamine cosolute. The domain provides a rare example of a molten globule that can be crystallized. The structure that is stabilized in the crystal represents the fully folded state of the domain, which must be transiently realized during binding to the viral nucleocapsid. While the intermolecular forces that govern the polymerase-nucleocapsid interaction appear to be different in measles, mumps, and Sendai viruses, for each of these viruses, polymerase translocation involves the coupled binding and folding of protein domains. In all cases, we suggest that this will result in a weak-affinity protein complex with a short lifetime, which allows the polymerase to take rapid steps forward.

Conformational stability and DNA binding energetics of the rat thyroid transcription factor 1 homeodomain

The conformational stability of the rat thyroid transcription factor 1 homeodomain, TTF-1HD, has been investigated by means of circular dichroism (CD) and differential scanning calorimetry (DSC) measurements at pH 5.0 as a function of KCl concentration. Thermal unfolding of TTF-1HD is a reversible two-state transition. The protein is not stable against temperature, showing a denaturation temperature of 32 degrees C in the absence of salt and 50 degrees C at 75 mM KCl. The binding energetics of TTF-1HD to its target DNA sequence has been characterized by means of isothermal titration calorimetry (ITC) measurements, complemented with CD data. At 25 degrees C, pH 5.0 and 75 mM KCl, the binding constant amounts to 1.5 x 10(8)M(-1) and the binding enthalpy change amounts to -41 kJ mol(-1). The process is enthalpy driven, but also the entropy change is favorable to complex formation. To gain a molecular level understanding of the interactions determining the association of TTF-1HD to the target DNA sequence structural information would be requested, but it is not yet available. Therefore, structural models of two complexes, TTF-1HD with the target DNA sequence and TTF-1HD with a modified DNA sequence, have been constructed by using as a template the NMR structure of the complex between NK-2 HD and its target DNA, and by performing molecular dynamics simulations 3.5 ns long. Analysis of these models allows one to shed light on the origin of the DNA binding specificity characteristic of TTF-1HD.

An α-Helical Burst in the src SH3 Folding Pathway

Src SH3 is a small all-beta-sheet protein composed of a single domain. We studied the folding behavior of src SH3 at various conditions by circular dichroism (CD), fluorescence, and X-ray solution scattering methods. On the src SH3 folding pathway, an alpha-helix-rich intermediate appeared not only at subzero temperatures but also above 0 degrees C. The fraction of alpha-helix in the kinetically observed intermediate is ca. 26% based on the kinetic CD experiment. X-ray solution scattering revealed that the intermediate was compact but not fully packed. The analysis of CD implies that the amplitude of the burst phase is proportional to the helical fraction calculated according to the helix-coil transition theory. This strongly suggests that the initial folding core is formed by the collapse of much less stably existing alpha-helices.

Unfolding thermodynamics of Trp-cage, a 20 residue miniprotein, studied by differential scanning calorimetry and circular dichroism spectroscopy

Small proteins provide convenient models for computational studies of protein folding and stability, which are usually compared with experimental data. Until recently, the unfolding of Trp-cage was considered to be a two-state process. However, no direct experimental evidence for this has been presented, and in some cases, the contrary has been suggested. To elucidate a detailed unfolding mechanism, we studied the thermodynamics of unfolding of Trp-cage by differential scanning calorimetry (DSC) and circular dichroism (CD) spectroscopy. The observation that at low temperatures only approximately 90-95% of Trp-cage exists in the native conformation presented an analytical challenge. Nevertheless, it was found that the DSC and CD data can be fitted simultaneously to the same set of thermodynamic parameters. The major uncertainty in such a global fit is the heat capacity change upon unfolding, DeltaCp. This can be circumvented by obtaining DeltaCp directly from the difference between heat capacity functions of the native and unfolded states. Using such an analysis it is shown that Trp-cage unfolding can be represented by a two-state model with the following thermodynamic parameters: Tm = 43.9 +/- 0.8 degrees C, DeltaH(Tm) = 56 +/- 2 kJ/mol, DeltaCp = 0.3 +/- 0.1 kJ/(mol.K). Using these thermodynamic parameters it is estimated that Trp-cage is marginally stable at 25 degrees C, DeltaG(25 degrees C) = 3.2 +/- 0.2 kJ/mol, which is only 30% more than the thermal fluctuation energy at this temperature.

Effects of monovalent anions on a temperature-dependent heat capacity change for Escherichia coli SSB tetramer binding to single-stranded DNA

We have previously shown that the linkage of temperature-dependent protonation and DNA base unstacking equilibria contribute significantly to both the negative enthalpy change (DeltaH(obs)) and the negative heat capacity change (DeltaC(p,obs)) for Escherichia coli SSB homotetramer binding to single-stranded (ss) DNA. Using isothermal titration calorimetry we have now examined DeltaH(obs) over a much wider temperature range (5-60 degrees C) and as a function of monovalent salt concentration and type for SSB binding to (dT)(70) under solution conditions that favor the fully wrapped (SSB)(65) complex (monovalent salt concentration >or=0.20 M). Over this wider temperature range we observe a strongly temperature-dependent DeltaC(p,obs). The DeltaH(obs) decreases as temperature increases from 5 to 35 degrees C (DeltaC(p,obs) <0) but then increases at higher temperatures up to 60 degrees C (DeltaC(p,obs) >0). Both salt concentration and anion type have large effects on DeltaH(obs) and DeltaC(p,obs). These observations can be explained by a model in which SSB protein can undergo a temperature- and salt-dependent conformational transition (below 35 degrees C), the midpoint of which shifts to higher temperature (above 35 degrees C) for SSB bound to ssDNA. Anions bind weakly to free SSB, with the preference Br(-) > Cl(-) > F(-), and these anions are then released upon binding ssDNA, affecting both DeltaH(obs) and DeltaC(p,obs). We conclude that the experimentally measured values of DeltaC(p,obs) for SSB binding to ssDNA cannot be explained solely on the basis of changes in accessible surface area (ASA) upon complex formation but rather result from a series of temperature-dependent equilibria (ion binding, protonation, and protein conformational changes) that are coupled to the SSB-ssDNA binding equilibrium. This is also likely true for many other protein-nucleic acid interactions.

A coarse-grained model for the formation of α helix with a noninteger period on simple cubic lattices

Periodicity is an important parameter in the characterization of a helix in proteins. In this work, a coarse-grained model for a homopolypeptide in simple cubic lattices has been extended to build an alpha helix with a noninteger period. The lattice model is based on the bond fluctuation algorithm in which bond lengths and orientations are altered in a quasicontinuous way, and the simulation is performed via dynamic Monte Carlo simulation. Hydrogen bonds are assumed to be formed between a virtual-carbonyl group in a residue and a downstream virtual-imino group in another residue. Consequently, the period of the formed alpha helix is a noninteger. An improved spatial correlation function has been suggested to quantitatively describe the periodicity of the helical conformation, by which helical period and helical persistent length can be calculated by statistics. The resultant periods are very close to 3.6 residues; the persistent length based upon the improved definition can be larger or smaller than the chain length and reflect the inherent regularity of a chain including one or multiple helical blocks. The simulation outputs agree with the calculation of the Zimm-Bragg theory based upon the associated nucleation parameter and propagation parameter as well.

Elimination of the C-cap in ubiquitin - structure, dynamics and thermodynamic consequences

Single amino acid substitutions rarely produce substantial changes in protein structure. Here we show that substitution of the C-cap residue in the alpha-helix of ubiquitin with proline (34P variant) leads to dramatic structural changes. The resulting conformational perturbation extends over the last two turns of the alpha-helix and leads to enhanced flexibility for residues 27-37. Thermodynamic analysis of this ubiquitin variant using differential scanning calorimetry reveals that the thermal unfolding transition remains highly cooperative, exhibiting two-state behavior. Similarities with the wild type in the thermodynamic parameters (heat capacity change upon unfolding and m-value) of unfolding monitored by DSC and chemical denaturation suggests that the 34P variant has comparable buried surface area. The hydrophobic core of 34P variant is not packed as well as that of the wild type protein as manifested by a lower enthalpy of unfolding. The increased mobility of the polypeptide chain of this ubiquitin variant allows the transient opening of the hydrophobic core as evidenced by ANS binding. Taken together, these results suggest exceptional robustness of cooperativity in protein structures.