Dynamics of apomyoglobin in the α-to-β transition and of partially unfolded aggregated protein

Direct Evidence for the Effect of Glycerol on Protein Hydration and Thermal Structural Transition

The mechanisms of protein stabilization by uncharged solutes, such as polyols and sugars, have been intensively studied with respect to the chemical thermodynamics of molecular crowding. In particular, many experimental and theoretical studies have been conducted to explain the mechanism of the protective action on protein structures by glycerol through the relationship between hydration and glycerol solvation on protein surfaces. We used wide-angle x-ray scattering (WAXS), small-angle neutron scattering, and theoretical scattering function simulation to quantitatively characterize the hydration and/or solvation shell of myoglobin in aqueous solutions of up to 75% v/v glycerol. At glycerol concentrations below ∼40% v/v, the preservation of the hydration shell was dominant, which was reasonably explained by the preferential exclusion of glycerol from the protein surface (preferential hydration). In contrast, at concentrations above 50% v/v, the partial penetration or replacement of glycerol into or with hydration-shell water (neutral solvation by glycerol) was gradually promoted. WAXS results quantitatively demonstrated the neutral solvation, in which the replacement of hydrated water by glycerol was proportional to the volume fraction of glycerol in the solvent multiplied by an exchange rate (β ≤ 1). These phenomena were confirmed by small-angle neutron scattering measurements. The observed WAXS data covered the entire hierarchical structure of myoglobin, ranging from tertiary to secondary structures. We separately analyzed the effect of glycerol on the thermal stability of myoglobin at each hierarchical structural level. The thermal transition midpoint temperature at each hierarchical structural level was raised depending on the glycerol concentration, with enhanced transition cooperativeness between different hierarchical structural levels. The onset temperature of the helix-to-cross β-sheet transition (the initial process of amyloid formation) was evidently elevated. However, oligomerization connected to fibril formation was suppressed, even at a low glycerol concentration.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

Complex kinetics and residual structure in the thermal unfolding of yeast triosephosphate isomerase

BACKGROUND

Saccharomyces cerevisiae triosephosphate isomerase (yTIM) is a dimeric protein that shows noncoincident unfolding and refolding transitions (hysteresis) in temperature scans, a phenomenon indicative of the slow forward and backward reactions of the native-unfolded process. Thermal unfolding scans suggest that no stable intermediates appear in the unfolding of yTIM. However, reported evidence points to the presence of residual structure in the denatured monomer at high temperature.

RESULTS

Thermally denatured yTIM showed a clear trend towards the formation of aggregation-prone, β-strand-like residual structure when pH decreased from 8.0 to 6.0, even though thermal unfolding profiles retained a simple monophasic appearance regardless of pH. However, kinetic studies performed over a relatively wide temperature range revealed a complex unfolding mechanism comprising up to three observable phases, with largely different time constants, each accompanied by changes in secondary structure. Besides, a simple sequential mechanism is unlikely to explain the observed variation of amplitudes and rate constants with temperature. This kinetic complexity is, however, not linked to the appearance of residual structure. Furthermore, the rate constant for the main unfolding phase shows small, rather unvarying values in the pH region where denatured yTIM gradually acquires a β-strand-like conformation. It appears, therefore, that the residual structure has no influence on the kinetic stability of the native protein. However, the presence of residual structure is clearly associated with increased irreversibility.

CONCLUSIONS

The slow temperature-induced unfolding of yeast TIM shows three kinetic phases. Rather than a simple sequential pathway, a complex mechanism involving off-pathway intermediates or even parallel pathways may be operating. β-strand-type residual structure, which appears below pH 8.0, is likely to be associated with increased irreversible aggregation of the unfolded protein. However, this denatured form apparently accelerates the refolding process.

Correlation between supercoiling and conformational motions of the bacterial flagellar filament

The bacterial flagellar filament is a very large macromolecular assembly of a single protein, flagellin. Various supercoiled states of the filament exist, which are formed by two structurally different conformations of flagellin in different ratios. We investigated the correlation between supercoiling of the protofilaments and molecular dynamics in the flagellar filament using quasielastic and elastic incoherent neutron scattering on the picosecond and nanosecond timescales. Thermal fluctuations in the straight L- and R-type filaments were measured and compared to the resting state of the wild-type filament. Amplitudes of motion on the picosecond timescale were found to be similar in the different conformational states. Mean-square displacements and protein resilience on the 0.1 ns timescale demonstrate that the L-type state is more flexible and less resilient than the R-type, whereas the wild-type state lies in between. Our results provide strong support that supercoiling of the protofilaments in the flagellar filament is determined by the strength of molecular forces in and between the flagellin subunits.

Unfolding simulations of holomyoglobin from four mammals: identification of intermediates and β-sheet formation from partially unfolded states

Myoglobin (Mb) is a centrally important, widely studied mammalian protein. While much work has investigated multi-step unfolding of apoMb using acid or denaturant, holomyoglobin unfolding is poorly understood despite its biological relevance. We present here the first systematic unfolding simulations of holoMb and the first comparative study of unfolding of protein orthologs from different species (sperm whale, pig, horse, and harbor seal). We also provide new interpretations of experimental mean molecular ellipticities of myoglobin intermediates, notably correcting for random coil and number of helices in intermediates. The simulated holoproteins at 310 K displayed structures and dynamics in agreement with crystal structures (R g ~1.48-1.51 nm, helicity ~75%). At 400 K, heme was not lost, but some helix loss was observed in pig and horse, suggesting that these helices are less stable in terrestrial species. At 500 K, heme was lost within 1.0-3.7 ns. All four proteins displayed exponentially decaying helix structure within 20 ns. The C- and F-helices were lost quickly in all cases. Heme delayed helix loss, and sperm whale myoglobin exhibited highest retention of heme and D/E helices. Persistence of conformation (RMSD), secondary structure, and ellipticity between 2-11 ns was interpreted as intermediates of holoMb unfolding in all four species. The intermediates resemble those of apoMb notably in A and H helices, but differ substantially in the D-, E- and F-helices, which interact with heme. The identified mechanisms cast light on the role of metal/cofactor in poorly understood holoMb unfolding. We also observed β-sheet formation of several myoglobins at 500 K as seen experimentally, occurring after disruption of helices to a partially unfolded, globally disordered state; heme reduced this tendency and sperm-whale did not display any sheet propensity during the simulations.

Misfolding and amyloid aggregation of apomyoglobin

Apomyoglobin is an excellent example of a monomeric all α-helical globular protein whose folding pathway has been extensively studied and well characterized. Structural perturbation induced by denaturants or high temperature as well as amino acid substitution have been described to induce misfolding and, in some cases, aggregation. In this article, we review the molecular mechanism of the aggregation process through which a misfolded form of a mutated apomyoglobin aggregates at physiological pH and room temperature forming an amyloid fibril. The results are compared with data showing that either amyloid or aggregate formation occurs under particular denaturing conditions or upon cleavage of the residues corresponding to the C-terminal helix of apomyoglobin. The results are discussed in terms of the sequence regions that are more important than others in determining the amyloid aggregation process.

Change of dynamics of raft-model membrane induced by amyloid-β protein binding

While the steady-state existence in the size and shape of liquid-ordered microdomains in cell membranes, the so-called "lipid rafts", still remain the subject of debate, glycosphingolipid-cholesterol rich regions in plasma membranes have been considered to have a function as platforms for signaling and sorting. In addition, recent spectroscopic studies show that the interaction between monosialoganglioside and amyloid beta (Aβ protein promotes the transition of Aβ from the native structure to the cross-beta fold in amyloid aggregates. However, there is few evidence on the dynamics of "lipid rafts" membranes. As the neutron spin-echo (NSE) technique is well known to detect directly slow dynamics of membrane systems in situ, by the combination of NSE and small-angle X-ray scattering we have studied the effect of the interaction between raft-model membrane and amyloid Aβ proteins on the structure and dynamics of a large uni-lamellar vesicle (LUV) consisting of monosialoganglioside-cholesterol-phospholipid ternary mixtures as a model of lipid-raft membrane. We have found that the interaction between the Aβ proteins and the model membrane at the liquid crystal phase significantly suppresses a bending-diffusion motion with a minor effect on the LUV structure. The present results would suggest a possibility of non-receptor-mediated disorder in signaling through a modulation of a membrane dynamics induced by the association of amyloidogenic peptides on a plasma membrane.

Neutron scattering: a tool to detect in vivo thermal stress effects at the molecular dynamics level in micro-organisms

In vivo molecular dynamics in Halobacterium salinarum cells under stress conditions was measured by neutron scattering experiments coupled with microbiological characterization. Molecular dynamics alterations were detected with respect to unstressed cells, reflecting a softening of protein structures consistent with denaturation. The experiments indicated that the neutron scattering method provides a promising tool to study molecular dynamics modifications in the proteome of living cells induced by factors altering protein folds.

Thermal fluctuations of haemoglobin from different species: adaptation to temperature via conformational dynamics

Thermodynamic stability, configurational motions and internal forces of haemoglobin (Hb) of three endotherms (platypus, Ornithorhynchus anatinus; domestic chicken, Gallus gallus domesticus and human, Homo sapiens) and an ectotherm (salt water crocodile, Crocodylus porosus) were investigated using circular dichroism, incoherent elastic neutron scattering and coarse-grained Brownian dynamics simulations. The experimental results from Hb solutions revealed a direct correlation between protein resilience, melting temperature and average body temperature of the different species on the 0.1 ns time scale. Molecular forces appeared to be adapted to permit conformational fluctuations with a root mean square displacement close to 1.2 Å at the corresponding average body temperature of the endotherms. Strong forces within crocodile Hb maintain the amplitudes of motion within a narrow limit over the entire temperature range in which the animal lives. In fully hydrated powder samples of human and chicken, Hb mean square displacements and effective force constants on the 1 ns time scale showed no differences over the whole temperature range from 10 to 300 K, in contrast to the solution case. A complementary result of the study, therefore, is that one hydration layer is not sufficient to activate all conformational fluctuations of Hb in the pico- to nanosecond time scale which might be relevant for biological function. Coarse-grained Brownian dynamics simulations permitted to explore residue-specific effects. They indicated that temperature sensing of human and chicken Hb occurs mainly at residues lining internal cavities in the β-subunits.

Dynamics-stability relationships in apo- and holomyoglobin: a combined neutron scattering and molecular dynamics simulations study

The removal of the heme group from myoglobin (Mb) results in a destabilization of the protein structure. The dynamic basis of the destabilization was followed by comparative measurements on holo- (holo-Mb) and apomyoglobin (apo-Mb). Mean-squared displacements (MSD) and protein resilience on the picosecond-to-nanosecond timescale were measured by elastic incoherent neutron scattering. Differences in thermodynamic parameters, MSD, and resilience were observed for both proteins. The resilience of holo-Mb was significantly lower than that of apo-Mb, indicating entropic stabilization by a higher degree of conformational sampling in the heme-bound folded protein. Molecular dynamics simulations provided site-specific information. Averaged over the whole structure, the molecular dynamics simulations yielded similar MSD and resilience values for the two proteins. The mobility of residues around the heme group in holo-Mb showed a smaller MSD and higher resilience compared to the same residue group in apo-Mb. It is of interest that in holo-Mb, higher MSD values are observed for the residues outside the heme pocket, indicating an entropic contribution to protein stabilization by heme binding, which is in agreement with experimental results.

Neutron Scattering Reveals Enhanced Protein Dynamics in Concanavalin A Amyloid Fibrils

Protein aggregation is one of the most challenging topics in life sciences, and it is implicated in several human pathologies. The nature and the role of toxic species is highly debated, with amyloid fibrils being among the most relevant species for their peculiar structural and functional properties. Protein dynamics and in particular the ability to fluctuate through a large number of conformational substates are closely related to protein function. This Letter focuses on amyloid fibril dynamics, and, to our knowledge, it is the first neutron scattering study on a protein (Concanavalin A) isolated in its fibril state. Our results reveal enhanced atomic fluctuations in amyloid fibrils and indicate that the protein is "softer" in the fibril state with respect to the native and amorphous aggregate states. We discuss this finding in terms of a structural interpretation and suggest that the paradigm ordered structure ↔ lower flexibility can be questioned when considering the local fast side-chain protein dynamics.

Study of Protein Dynamics vs. Amyloid Formation

Protein fibril formation has been often associated to manifestation of serious and devastating amyloyd diseases, including Alzheimer and BSA. The main mechanism for the formation of amyloid fibrils is the accumulation of protein aggregates in body´s organs.In this paper, we try to compare the dynamical behaviour of two amyloidogenic proteins, the Insulin and the Myoglobin. Insulin has been chosen for its pharmacological extensive use in diabete´s therapy, while Myoglobin is used as control, since its dynamics is now largely known. The investigation has been performed through incoherent elastic neutron scattering over a wide temperature range. Our results suggest an enhanced stiffness of Insulin with respect to Myoglobin.

Dynamics in Biological Systems as seen by QENS

Quasielastic incoherent neutron scattering is a well suited and established experimental method to study protein and water dynamics in the picosecond to nanosecond time- and Ångstrom length-scale. Using deuterium labelling either protein or water motions can be selected and brought into focus. Protein and cell water dynamics were separately studied in red blood cells. A consistent picture of cytoplasmic water and protein dynamics in whole cells is emerging from recent experimental results.