Use of SANS and biophysical techniques to reveal subtle conformational differences between native apo-calmodulin and its unfolded states

Apo-calmodulin, a small, mainly α, soluble protein is a calcium-dependent protein activator. It is made of two N- and C-terminal domains having a sequence homology of 70%, an identical folding but different stabilities, and is thus an interesting system for unfolding studies. The use of small angle neutron scattering (SANS) and other biophysical techniques has permitted to reveal conformational difference between native and thermal denatured states of apo-calmodulin. The results show that secondary and tertiary structures of apo-calmodulin evolve in a synchronous way, indicating the absence in the unfolding pathway of molten-globule state sufficiently stable to affect transition curves. From SANS experiments, at 85 °C, apo-calmodulin adopts a polymer chain conformation with some residual local structures. After cooling down, apo-calmodulin recovers a compact state, with a secondary structure close to the native one but with a higher radius of gyration and a different tyrosine environment. In fact on a timescale of few minutes, heat denaturation of apo-calmodulin is partially reversible, but on a time scale of hours (for SANS experiments), the long exposure to heat may lead to a non-reversibility due to some chemical perturbation of the protein. In fact, from Mass Spectrometry measurements, we got evidence of dehydration and deamidation of heated apo-calmodulin.     

Dynamics and dimension of an amyloidogenic disordered state of human β2-microglobulin

Human β₂-microglobulin (β₂m) aggregation is implicated in dialysis-related amyloidosis. Previously, it has been shown that β₂m adopts an ensemble of partially unfolded states at low pH. Here we provide detailed structural and dynamical insights into the acid unfolded and yet compact state of β₂m at pH 2.5 using a host of fluorescence spectroscopic tools. These tools allowed us to investigate protein conformational dynamics at low micromolar protein concentrations in an amyloid-forming condition. Our equilibrium fluorescence data in combination with circular dichroism data provide support in favor of progressive structural dissolution of β₂m with lowering pH. The acid unfolded intermediate at pH 2.5 has high 8-anilinonaphthalene, 1-sulfonic acid (ANS)-binding affinity and is devoid of significant secondary structural elements. Using fluorescence lifetime measurements, we have been able to monitor the conformational transition during the pH transition from the native to the compact disordered state. Additionally, using time-resolved fluorescence anisotropy measurements, we have been able to distinguish this compact disordered state from the canonical denatured state of the protein by identifying unique dynamic signatures pertaining to the segmental chain mobility. Taken together, our results demonstrate that β₂m at pH 2.5 adopts a compact noncanonical unfolded state resembling a collapsed premolten globule state. Additionally, our stopped-flow fluorescence kinetics results provide mechanistic insights into the formation of a compact disordered state from the native form.     

Thermodynamics and Kinetics of Single-Chain Monellin Folding with Structural Insights into Specific Collapse in the Denatured State Ensemble

Proteins, which behave as random coils in high denaturant concentrations undergo collapse transition similar to polymers on denaturant dilution. We study collapse in the denatured ensemble of single-chain monellin (MNEI) using a coarse-grained protein model and molecular dynamics simulations. The model is validated by quantitatively comparing the computed guanidinium chloride and pH-dependent thermodynamic properties of MNEI folding with the experiments. The computed properties such as the fraction of the protein in the folded state and radius of gyration (R_g) as function of [GuHCl] are in good agreement with the experiments. The folded state of MNEI is destabilized with an increase in pH due to the deprotonation of the residues Glu24 and Cys42. On decreasing [GuHCl], the protein in the unfolded ensemble showed specific compaction. The R_g of the protein decreased steadily with [GuHCl] dilution due to increase in the number of native contacts in all the secondary structural elements present in the protein. MNEI folding kinetics is complex with multiple folding pathways and transiently stable intermediates are populated in these pathways. In strong stabilizing conditions, the protein in the unfolded ensemble showed transition to a more compact unfolded state where R_g decreased by ≈ 17% due to the formation of specific native contacts in the protein. The intermediate populated in the dominant MNEI folding pathway satisfies the structural features of the dry molten globule inferred from experiments.     

Translational hydration water dynamics drives the protein glass transition

Experimental and computer simulation studies have revealed the presence of a glass-like transition in the internal dynamics of hydrated proteins at approximately 200 K involving an increase of the amplitude of anharmonic dynamics. This increase in flexibility has been correlated with the onset of protein activity. Here, we determine the driving force behind the protein transition by performing molecular dynamics simulations of myoglobin surrounded by a shell of water. A dual heat bath method is used with which, in any given simulation, the protein and solvent are held at different temperatures, and sets of simulations are performed varying the temperature of the components. The results show that the protein transition is driven by a dynamical transition in the hydration water that induces increased fluctuations primarily in side chains in the external regions of the protein. The water transition involves activation of translational diffusion and occurs even in simulations where the protein atoms are held fixed.     

Structural characteristics of thermostable immunogenic outer membrane protein from Salmonella enterica serovar Typhi

In this work, we explored the acid-induced unfolding pathway of non-porin outer membrane protein (OMP), an immunogenic protein from Salmonella Typhi, by monitoring the conformational changes over a pH range of 1.0–7.0 by circular dichroism, intrinsic fluorescence, ANS binding, acrylamide quenching, and dynamic light scattering. The spectroscopic measurements showed that OMP in its native state at pH 7.0 exists in more stable and compact conformation. In contrast, at pH 2.0, OMP retains substantial amount of secondary structure, disrupted side chain interactions, increased hydrodynamic radii, and nearly four-fold increase in ANS fluorescence with respect to the native state, indicating that MG state exists at pH 2.0. Quenching of tryptophan fluorescence by acrylamide further confirmed the accumulation of a partially unfolded state between native and unfolded state. The effect of pH on the conformation and thermostability of OMP points towards its heat resistance at neutral pH (T _m?~?69 °C at pH 7.0, monitored by change in MRE_{222 nm}). Acid unfolded state was also characterized by the lack of a cooperative thermal transition. All these results suggested that acid-induced unfolded state of OMP at pH 2.0 represented the molten globule state. The chemical denaturation studies with GuHCl and urea as denaturants showed dissimilar results. The chemical unfolding experiments showed that in both far-UV CD and fluorescence measurements, GuHCl is more efficient than urea. GuHCl is characterized by low C _m (~1 M), while urea is characterized by high C _m (~3 M). The fully unfolded states were reached at 2 M GuHCl and 4 M urea concentration, respectively. This study adds to several key considerations of importance in the development of therapeutic agents against typhoid fever for clinical purposes.     

Controlling the protein dynamical transition with sugar-based bioprotectant matrices: a neutron scattering study

Through elastic neutron scattering we measured the mean-square displacements of the hydrogen atoms of lysozyme embedded in a glucose-water glassy matrix as a function of the temperature and at various water contents. The elastic intensity of all the samples has been interpreted in terms of the double-well model in the whole temperature range. The dry sample shows an onset of anharmonicity at approximately 100 K, which can be attributed to the activation of methyl group reorientations. Such a protein intrinsic dynamics is decoupled from the external environment on the whole investigated temperature range. In the hydrated samples an additional and larger anharmonic contribution is provided by the protein dynamical transition, which appears at a higher temperature Td. As hydration increases the coupling between the protein internal dynamics and the surrounding matrix relaxations becomes more effective. The behavior of Td that, as a function of the water content, diminishes by approximately 60 K, supports the picture of the protein dynamics as driven by solvent relaxations. A possible connection between the protein dynamical response versus T and the thermal stability in glucose-water bioprotectant matrices is proposed.     

Exploration of the correlation between solvation dynamics and internal dynamics of a protein

Protein function is intimately related to the dynamics of the protein as well as to the dynamics of the solvent shell around the protein. Although it has been argued extensively that protein dynamics is slaved to solvent dynamics, experimental support for this hypothesis is scanty. In this study, measurements of fluorescence anisotropy decay kinetics have been used to determine the motional dynamics of the fluorophore acrylodan linked to several locations in a small protein barstar in its various structural forms, including the native and unfolded states as well as the acid and protofibril forms. Fluorescence upconversion and streak camera measurements have been used to determine the solvation dynamics around the fluorophore. Both the motional dynamics and solvent dynamics were found to be dependent upon the location of the probe as well as on the structural form of the protein. While the (internal) motional dynamics of the fluorophore occur in the 0.1-3 ns time domain, the observed mean solvent relaxation times are in the range of 20-300 ps. A strong positive correlation between these two dynamical modes was found in spite of the significant difference in their time scales. This observed correlation is a strong indicator of the coupling between solvent dynamics and the dynamics in the protein.     

Dependence of Protein Folding Stability and Dynamics on the Density and Composition of Macromolecular Crowders

We investigate the effect of macromolecular crowding on protein folding, using purely repulsive crowding particles and a self-organizing polymer model of protein folding. We find that the variation in folding stability with crowder size for typical α-, β-, and α/β-proteins is well described by an adaptation of the scaled particle theory. The native state, the transition state, and the unfolded protein are treated as effective hard spheres, with the folded and transition state radii independent of the size and concentration of the crowders. Remarkably, we find that, as the effective unfolded state radius is very weakly dependent on the crowder concentration, it can also be approximated by a single size. The same model predicts the effect of crowding on the folding barrier and therefore refolding rates with no adjustable parameters. A simple extension of the scaled-particle theory model, assuming additivity, can also describe the behavior of mixtures of crowding particles.     

Structure of stable protein folding intermediates by equilibrium phi-analysis: the apoflavodoxin thermal intermediate

Protein intermediates in equilibrium with native states may play important roles in protein dynamics but, in cases, can initiate harmful aggregation events. Investigating equilibrium protein intermediates is thus important for understanding protein behaviour (useful or pernicious) but it is hampered by difficulties in gathering structural information. We show here that the phi-analysis techniques developed to investigate transition states of protein folding can be extended to determine low-resolution three-dimensional structures of protein equilibrium intermediates. The analysis proposed is based solely on equilibrium data and is illustrated by determination of the structure of the apoflavodoxin thermal unfolding intermediate. In this conformation, a large part of the protein remains close to natively folded, but a 40 residue region is clearly unfolded. This structure is fully consistent with the NMR data gathered on an apoflavodoxin mutant designed specifically to stabilise the intermediate. The structure shows that the folded region of the intermediate is much larger than the proton slow-exchange core at 25 degrees C. It also reveals that the unfolded region is made of elements whose packing surface is more polar than average. In addition, it constitutes a useful guide to rationally stabilise the native state relative to the intermediate state, a far from trivial task.     

Recombination of S-peptide with S-protein during folding of ribonuclease S. II. Kinetic characterization of a stable folding intermediate shown by S-protein at pH 1.7.

At pH 1.7 S-peptide dissociates from S-protein but S-protein remains partly folded below 30 °C. A folded form of S-protein, labeled I₃, is detected and measured by its ability to combine rapidly with S-peptide at pH 6.8 and then to form native ribonuclease S. The second-order combination reaction (k = 0.7 × 10⁶m^?1s^?1 at 20 °C) can be monitored either by tyrosine absorbance or fluorescence emission; the subsequent first-order folding reaction (half-time, 68 ms; 20 °C) is monitored by 2′CMP ² binding. Combination with S-peptide and folding to form native RNAase S is considerably slower for both classes of unfolded S-protein (see preceding paper).I₃ shows a thermal folding transition at pH 1.7: it is completely unfolded above 32 °C and reaches a limiting low-temperature value of 65% below 10 °C. The 35% S-protein remaining at 10 °C is unfolded as judged by its refolding behavior in forming native RNAase S at pH 6.8. The folding transition of S-protein at pH 1.7 is a broad, multi-state transition. This is shown both by the large fraction of unfolded S-protein remaining at low temperatures and by the large differences between the folding transition curves monitored by I₃ and by tyrosine absorbance.The fact that S-protein remains partly folded after dissociation of S-peptide at pH 1.7 but not at pH 6.8 may be explained by two earlier observations. (1) Native RNAase A is stable in the temperature range of the S-protein folding transition at pH 1.7, and (2) the binding constant of S-protein for S-peptide falls steadily as the pH is lowered, by more than four orders of magnitude between pH 8.3 and pH 2.7, at 0 °C. The following explanation is suggested for why folding intermediates are observed easily in the transition of S-protein but not of RNAase A. The S-protein transition is shifted to lower temperatures, where folding intermediates should be more stable: consequently, intermediates in the folding of RNAase A which do not involve the S-peptide moiety and which are populated to almost detectable levels can be observed at the lower temperatures of the S-protein transition.     

Deconvoluting Protein (Un)folding Structural Ensembles Using X-Ray Scattering,Nuclear Magnetic Resonance Spectroscopy and Molecular Dynamics Simulation

The folding and unfolding of protein domains is an apparently cooperative process, but transient intermediates have been detected in some cases. Such (un)folding intermediates are challenging to investigate structurally as they are typically not long-lived and their role in the (un)folding reaction has often been questioned. One of the most well studied (un)folding pathways is that of Drosophila melanogaster Engrailed homeodomain (EnHD): this 61-residue protein forms a three helix bundle in the native state and folds via a helical intermediate. Here we used molecular dynamics simulations to derive sample conformations of EnHD in the native, intermediate, and unfolded states and selected the relevant structural clusters by comparing to small/wide angle X-ray scattering data at four different temperatures. The results are corroborated using residual dipolar couplings determined by NMR spectroscopy. Our results agree well with the previously proposed (un)folding pathway. However, they also suggest that the fully unfolded state is present at a low fraction throughout the investigated temperature interval, and that the (un)folding intermediate is highly populated at the thermal midpoint in line with the view that this intermediate can be regarded to be the denatured state under physiological conditions. Further, the combination of ensemble structural techniques with MD allows for determination of structures and populations of multiple interconverting structures in solution.     

Cation-induced toroidal condensation of DNA studies with Co3+(NH3)6

The unfolding and refolding of Staphylococcus aureus penicillinase have been followed by urea-gradient electrophoresis. Unfolding of the native state proceeds by an all-or-none transition to fully unfolded protein, with no detectable accumulation of partially unfolded states. In contrast, refolding is complex and proceeds by very rapid, reversible formation of a partially folded state, H, which had been detected and characterized previously, as it is the most stable conformation at intermediate denaturant concentrations. At very low urea concentrations, a more compact conformational state was observed as a transient intermediate in refolding. There was little kinetic heterogeneity of the unfolded protein, as is normally observed with proteins containing proline residues.     

Picosecond Internal Dynamics of Lysozyme as Affected by Thermal Unfolding in Nonaqueous Environment

A neutron-scattering investigation of the internal picosecond dynamics of lysozyme solvated in glycerol as a function of temperature in the range 200–410 K has been undertaken. The inelastic contribution to the measured intensity is characterized by the presence of a bump generally known as “boson peak”, clearly distinguishable at low temperature. When the temperature is increased the quasielastic component of the spectrum becomes more and more intrusive and progressively overwhelms the vibrational bump. This happens especially for T > 345 K when the protein goes through an unfolding process, which leads to the complete denaturation. The quasielastic term is the superposition of two components whose intensities and linewidths have been studied as a function of temperature. The slower component describes motions with characteristic times of ~4 ps corresponding to reorientations of polypeptide side chains. Both the intensity and linewidth of this kind of relaxations show two distinct regimes with a crossover in the temperature range where the melting process occurs, thus suggesting the presence of a dynamical transition correlated to the protein unfolding. Conversely the faster component might be ascribed to the local dynamics of hydrogen atoms caged by the nearest neighbors with characteristic time of ~0.3 ps.     

Characterizing intermolecular interactions that initiate native-like protein aggregation

Folded proteins can access aggregation-prone states without the need for transitions that cross the energy barriers for unfolding. In this study we characterized the initial steps of aggregation from a native-like state of the acylphosphatase from Sulfolobus solfataricus (Sso AcP). Using computer simulations restrained by experimental hydrogen/deuterium (H/D) exchange data, we provide direct evidence that under aggregation-promoting conditions Sso AcP populates a conformational ensemble in which native-like structure is retained throughout the sequence in the absence of local unfolding (N¹), although the protein exhibits an increase in hydrodynamic radius and dynamics. This transition leads an edge strand to experience an increased affinity for a specific unfolded segment of the protein. Direct measurements by means of H/D exchange rates, isothermal titration calorimetry, and intermolecular relaxation enhancements show that after formation of N¹, an intermolecular interaction with an antiparallel arrangement is established between the edge strand and the unfolded segment of the protein. However, under conditions that favor the fully native state of Sso AcP, such an interaction is not established. Thus, these results reveal a novel (to our knowledge) self-assembly mechanism for a folded protein that is based on the increased flexibility of highly aggregation-prone segments in the absence of local unfolding.     

The protein-solvent glass transition

The protein dynamical transition and its connection with the liquid-glass transition (GT) of hydration water and aqueous solvents are reviewed. The protein solvation shell exhibits a regular glass transition, characterized by steps in the specific heat and the thermal expansion coefficient at the calorimetric glass temperature T_G ≈ 170 K. It implies that the time scale of the structural α-relaxation has reached the experimental time window of 1–100 s. The protein dynamical transition, identified from elastic neutron scattering experiments by enhanced amplitudes of molecular motions exceeding the vibrational level [1], probes the α-process on a shorter time scale. The corresponding liquid-glass transition occurs at higher temperatures, typically 240 K. The GT is generally associated with diverging viscosities, the freezing of long-range translational diffusion in the supercooled liquid. Due to mutual hydrogen bonding, both, protein- and solvent relaxational degrees of freedom slow down in paralled near the GT. However, the freezing of protein motions, where surface-coupled rotational and librational degrees of freedom are arrested, is better characterized as a rubber-glass transition. In contrast, internal protein modes such as the rotation of side chains are not affected. Moreover, ligand binding experiments with myoglobin in various glass-forming solvents show, that only ligand entry and exit rates depend on the local viscosity near the protein surface, but protein-internal ligand migration is not coupled to the solvent. The GT leads to structural arrest on a macroscopic scale due to the microscopic cage effect on the scale of the intermolecular distance. Mode coupling theory provides a theoretical framework to understand the microcopic nature of the GT even in complex systems. The role of the α- and β-process in the dynamics of protein hydration water is evaluated. The protein-solvent GT is triggered by hydrogen bond fluctuations, which give rise to fast β-processes. High-frequency neutron scattering spectra indicate increasing hydrogen bond braking above T_G.     

Ligand Entry into Fatty Acid Binding Protein via Local Unfolding Instead of Gap Widening

Fatty acid binding proteins play an important role in the transportation of fatty acids. Despite intensive studies, how fatty acids enter the protein cavity for binding is still controversial. Here, a gap-closed variant of human intestinal fatty acid binding protein was generated by mutagenesis, in which the gap is locked by a disulfide bridge. According to its structure determined here by NMR, this variant has no obvious openings as the ligand entrance and the gap cannot be widened by internal dynamics. Nevertheless, it still takes up fatty acids and other ligands. NMR relaxation dispersion, chemical exchange saturation transfer, and hydrogen-deuterium exchange experiments show that the variant exists in a major native state, two minor native-like states, and two locally unfolded states in aqueous solution. Local unfolding of either βB–βD or helix 2 can generate an opening large enough for ligands to enter the protein cavity, but only the fast local unfolding of helix 2 is relevant to the ligand entry process.     

(Un)Folding Mechanisms of the FBP28 WW Domain in Explicit Solvent Revealed by Multiple Rare Event Simulation Methods

We report a numerical study of the (un)folding routes of the truncated FBP28 WW domain at ambient conditions using a combination of four advanced rare event molecular simulation techniques. We explore the free energy landscape of the native state, the unfolded state, and possible intermediates, with replica exchange molecular dynamics. Subsequent application of bias-exchange metadynamics yields three tentative unfolding pathways at room temperature. Using these paths to initiate a transition path sampling simulation reveals the existence of two major folding routes, differing in the formation order of the two main hairpins, and in hydrophobic side-chain interactions. Having established that the hairpin strand separation distances can act as reasonable reaction coordinates, we employ metadynamics to compute the unfolding barriers and find that the barrier with the lowest free energy corresponds with the most likely pathway found by transition path sampling. The unfolding barrier at 300 K is ∼17 k_BT ≈ 42 kJ/mol, in agreement with the experimental unfolding rate constant. This work shows that combining several powerful simulation techniques provides a more complete understanding of the kinetic mechanism of protein folding.     

Design and Structure of an Equilibrium Protein Folding Intermediate: A Hint into Dynamical Regions of Proteins

Partly unfolded protein conformations close to the native state may play important roles in protein function and in protein misfolding. Structural analyses of such conformations which are essential for their fully physicochemical understanding are complicated by their characteristic low populations at equilibrium. We stabilize here with a single mutation the equilibrium intermediate of apoflavodoxin thermal unfolding and determine its solution structure by NMR. It consists of a large native region identical with that observed in the X-ray structure of the wild-type protein plus an unfolded region. Small-angle X-ray scattering analysis indicates that the calculated ensemble of structures is consistent with the actual degree of expansion of the intermediate. The unfolded region encompasses discontinuous sequence segments that cluster in the 3D structure of the native protein forming the FMN cofactor binding loops and the binding site of a variety of partner proteins. Analysis of the apoflavodoxin inner interfaces reveals that those becoming destabilized in the intermediate are more polar than other inner interfaces of the protein. Natively folded proteins contain hydrophobic cores formed by the packing of hydrophobic surfaces, while natively unfolded proteins are rich in polar residues. The structure of the apoflavodoxin thermal intermediate suggests that the regions of natively folded proteins that are easily responsive to thermal activation may contain cores of intermediate hydrophobicity.