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.     

Effect of conformational states on protein dynamical transition

In order to examine the properties specific to the folded protein, the effect of the conformational states on protein dynamical transition was studied by incoherent elastic neutron scattering for both wild type and a deletion mutant of staphylococcal nuclease. The deletion mutant of SNase which lacks C-terminal 13 residues takes a compact denatured structure, and can be regarded as a model of intrinsic unstructured protein. Incoherent elastic neutron scattering experiments were carried out at various temperature between 10 K and 300 K on IN10 and IN13 installed at ILL. Temperature dependence of mean-square displacements was obtained by the q-dependence of elastic scattering intensity. The measurements were performed on dried and hydrated powder samples. No significant differences were observed between wild type and the mutant for the hydrated samples, while significant differences were observed for the dried samples. A dynamical transition at ∼ 140 K observed for both dried and hydrated samples. The slopes of the temperature dependence of MSD before transition and after transition are different between wild type and the mutant, indicating the folding induces hardening. The hydration water activates a further transition at ∼ 240 K. The behavior of the temperature dependence of MSD is indistinguishable for wild type and the mutant, indicating that hydration water dynamics dominate the dynamical properties.     

Conditioning action of the environment on the protein dynamics studied through elastic neutron scattering

The dynamics of lysozyme in the picosecond timescale has been studied when it is in dry and hydrated powder form and when it is embedded in glycerol, glycerol–water, glucose and glucose–water matrices. The investigation has been undertaken through elastic neutron scattering technique on the backscattering spectrometer IN13. The dynamics of dry powder and embedded-in-glucose lysozyme can be considered purely vibrational up to 100 K, where the onset of an anharmonic contribution takes place. This contribution can be attributed to the activation of methyl group reorientations and is described with an Arrhenius trend. An additional source of anharmonic dynamics appears at higher temperatures for lysozyme in hydrated powders and embedded in glycerol, glycerol–water and glucose–water matrices. This second process, also represented with an Arrhenius trend, corresponds to the so-called protein dynamical transition. Both the temperature where such a transition takes place and the magnitude of the protein mean square displacements depend on the environment. The dynamical response of the protein to temperature is put in relationship with its thermal stability.     

Small-angle neutron scattering study of the structure of protein/detergent complexes

Small-angle neutron scattering (SANS) was used to study the structure of protein/sodium dodecylsulfate complexes. Two water soluble proteins, bovine serum albumin (BSA) and ovalbumin (OVA), were used. The protein concentration was kept constant at 1 wt %, and protein/detergent wt ratio varied between 1/1, 1/1.5, 1/2 and 1/3. Absolute intensities of SANS distributions were analyzed by a fractal model. Analyses of large Q portions of SANS distributions established that sodium dodecylsulfate (SDS) molecules bound to a protein/SDS complex form micelle-like clusters. On the other hand, analyses of small Q portions of SANS distributions clearly showed that the arrangement of micelle-like clusters resembles a fractal packing of spheres. We showed that a protein/SDS complex can be characterized by four parameters extracted from the scattering experiment, namely, the average micelle size and its aggregation number, the fractal dimension characterizing the conformation of the micellar chains, the correlation length giving the extent of the unfolded polypeptide chains, and the numbers of micelle-like clusters in the complex.     

Interaction of the full-length Bax protein with biomimetic mitochondrial liposomes: a small-angle neutron scattering and fluorescence study

In response to apoptotic stimuli, the pro-apoptotic protein Bax inserts in the outer mitochondrial membrane, resulting in the formation of pores and the release of several mitochondrial components, and sealing the cell's fate. To study the binding of Bax to membranes, we used an in vitro system consisting of 50nm diameter liposomes prepared with a lipid composition mimicking that of mitochondrial membranes in which recombinant purified full-length Bax was inserted via activation with purified tBid. We detected the association of the protein with the membrane using fluorescence fluctuation methods, and found that it could well be described by an equilibrium between soluble and membrane-bound Bax and that at a high protein-to-liposome ratio the binding seemed to saturate at about 15 Bax proteins per 50nm diameter liposome. We then obtained structural data for samples in this saturated binding regime using small-angle neutron scattering under different contrast matching conditions. Utilizing a simple model to fit the neutron data, we observed that a significant amount of the protein mass protrudes above the membrane, in contrast to the conjecture that all of the membrane-associated Bax states are umbrella-like. Upon protein binding, we also observed a thinning of the lipid bilayer accompanied by an increase in liposome radius, an effect reminiscent of the action of antimicrobial peptides on membranes.     

Hydration of DNA: a neutron scattering study of oriented NaDNA

A preliminary neutron-scattering study of wet-spun oriented NaDNA samples with various contents of H₂O or D₂O has been performed to obtain information about the hydration of NaDNA. The preferred orientation was in the direction of the crystallographic c axis. Perpendicular to this direction the samples could be regarded as poly-crystalline. It has been found that the preferred crystallite orientation was also reflected in the water structure. This fact, taken together with the small energy widths of the measured near-elastic peaks indicates that water molecules in the sample are relatively strongly bound to the DNA helix either directly or via other water molecules. Angular distribution measurements of elastically scattered neutrons were performed with the momentum transfer vector directed along and perpendicular to the DNA helix axis. Different Debye-Waller factors were obtained, indicating different water dynamics in the two directions. Comparison between the H₂O and D₂O measurements suggests that the difference is caused by motions of rotational nature. The shape of the frequency distribution function for the H₂O molecules in the NaDNA sample was somewhat similar to that which has been measured in ice.     

Neutron frequency windows and the protein dynamical transition

Proteins undergo an apparent dynamical transition on temperature variation that has been correlated with the onset of function. The transition in the mean-square displacement, , that is observed using a spectrometer or computer simulation, depends on the relationship between the timescales of the relaxation processes activated and the timescale accessible to the instrument or simulation. Models are described of two extreme situations---an "equilibrium" model, in which the long-time dynamics changes with temperature and all motions are resolved by the instrument used; and a "frequency window" model, in which there is no change in the long-time dynamics but as the temperature increases, the relaxation frequencies move into the instrumental range. Here we demonstrate that the latter, frequency-window model can describe the temperature and timescale dependences of both the intermediate neutron scattering function and derived from molecular dynamics simulations of a small protein in a cryosolution. The frequency-window model also describes the energy-resolution and temperature-dependences of obtained from experimental neutron scattering on glutamate dehydrogenase in the same solvent. Although equilibrium effects should also contribute to dynamical transitions in proteins, the present results suggests that frequency-window effects can play a role in the simulations and experiments examined. Finally, misquotations of previous findings are discussed in the context of solvent activation of protein dynamics and the possible relationship of this to activity.     

Picosecond dynamical changes on denaturation of yeast phosphoglycerate kinase revealed by quasielastic neutron scattering

Quasielastic neutron scattering experiments performed on yeast phosphoglycerate kinase in the native form and denatured in 1.5 M guanidinium chloride reveal a change in the fast (picosecond time scale) diffusive internal dynamics of the protein. The momentum and energy transfer dependences of the scattering for both states are fitted by an analytical model in which, on the experimentally accessible picosecond time scale and angstrom length scale, the dynamics of a fraction of the nonexchangeable hydrogens in the protein is described as a superposition of vibrations with uniform diffusion in a sphere, the rest of the hydrogens undergoing only vibrational motion. The fraction diffusing changes, from ≈60% in the native protein to ≈82% in the denatured protein. The radius of the sphere also changes slightly, from ≈1.8 Å in the native protein to ≈2.2 Å in the denatured protein. Possible implications of these results for the general protein folding problem are discussed. Proteins 28:380–387, 1997 © 1997 Wiley-Liss, Inc.     

Morphology of fast-tumbling bicelles: a small angle neutron scattering and NMR study

Bilayered micelles, or bicelles, which consist of a mixture of long- and short-chain phospholipids, are a popular model membrane system. Depending on composition, concentration, and temperature, bicelle mixtures may adopt an isotropic phase or form an aligned phase in magnetic fields. Well-resolved (1)H NMR spectra are observed in the isotropic or so-called fast-tumbling bicelle phase, over the range of temperatures investigated (10-40 degrees C), for molar ratios of long-chain lipid to short-chain lipid between 0.20 and 1.0. Small angle neutron scattering data of this phase are consistent with the model in which bicelles were proposed to be disk-shaped. The experimentally determined dimensions are roughly consistent with the predictions of R.R. Vold and R.S. Prosser (J. Magn. Reson. B 113 (1996)). Differential paramagnetic shifts of head group resonances of dimyristoylphosphatidylcholine (DMPC) and dihexanoylphosphatidylcholine (DHPC), induced by the addition of Eu(3+), are also consistent with the bicelle model in which DHPC is believed to be primarily sequestered to bicelle rims. Selective irradiation of the DHPC aliphatic methyl resonances results in no detectable magnetization transfer to the corresponding DMPC methyl resonances (and vice versa) in bicelles, which also suggests that DHPC and DMPC are largely sequestered in the bicelle. Finally, (1)H spectra of the antibacterial peptide indolicidin (ILPWKWPWWPWRR-NH(2)) are compared, in a DPC micellar phase and the above fast-tumbling bicellar phases for a variety of compositions. The spectra exhibit adequate resolution and improved dispersion of amide and aromatic resonances in certain bicelle mixtures.     

Molecular dynamics decomposition of temperature-dependent elastic neutron scattering by a protein solution

Molecular dynamics simulations are performed of bovine pancreatic trypsin inhibitor in a cryosolution over a range of temperatures from 80 to 300 K and the origins identified of elastic dynamic neutron scattering from the solution. The elastic scattering and mean-square displacement calculated from the molecular dynamics trajectories are in reasonable agreement with experiments on a larger protein in the same solvent. The solvent and protein contributions to the scattering from the simulation model are determined. At lower temperatures (< approximately 200 K) or on shorter timescales ( approximately 10 ps) the scattering contributions are proportional to the isotopic nuclear scattering cross-sections of each component. However, for T > 200 K marked deviations from these cross-sections are seen due to differences in the dynamics of the components of the solution. Rapid activation of solvent diffusion leads to the variation with temperature of the total elastic intensity being determined largely by that of the solvent. At higher temperatures (>240 K) and longer times ( approximately 100 ps) the protein makes the only significant contribution to the scattering, the solvent scattering having moved out of the accessible time-space window. Decomposition of the protein mean-square displacement shows that the observed dynamical transition in the solution at 200-220 K involves activation of both internal motions and external whole-molecule rotational and translational diffusion. The proportion that the external dynamics contributes to the protein mean-square displacement increases to approximately 30 and 60% at 300 K on the 10- and 100-ps timescales, respectively.     

Complexes of RecA protein in solution. A study by small angle neutron scattering

RecA complexes on DNA and self-polymers were analysed by small-angle neutron scattering in solution. By Guinier analysis at small angles and by model analysis of a subsidiary peak at wider angles, we find that the filaments fall into two groups: the DNA complex in the presence of ATP gamma S, an open helix with pitch 95 A, a cross-sectional radius of gyration of 33 A and a mass per length of about six RecA units per turn, which corresponds to the state of active enzyme; and the compact form (bound to single-stranded DNA in the absence of ATP, or binding ATP gamma S in the absence of DNA, or just the protein on its own), a helical structure with pitch 70 A, cross-sectional radius of gyration 40 A and mass per length about five RecA units per turn, which corresponds to the conditions of inactive enzyme. The results are discussed in the perspective of unifying previous conflicting structural results obtained by electron microscopy.     

Cancer as a dynamical phase transition

Background

We review and extend the work of Rosen and Casti who discuss category theory with regards to systems biology and manufacturing systems, respectively.

Results

We describe anticipatory systems, or long-range feed-forward chemical reaction chains, and compare them to open-loop manufacturing processes. We then close the loop by discussing metabolism-repair systems and describe the rationality of the self-referential equation f = f (f). This relationship is derived from some boundary conditions that, in molecular systems biology, can be stated as the cardinality of the following molecular sets must be about equal: metabolome, genome, proteome. We show that this conjecture is not likely correct so the problem of self-referential mappings for describing the boundary between living and nonliving systems remains an open question. We calculate a lower and upper bound for the number of edges in the molecular interaction network (the interactome) for two cellular organisms and for two manufacturomes for CMOS integrated circuit manufacturing.

Conclusions

We show that the relevant mapping relations may not be Abelian, and that these problems cannot yet be resolved because the interactomes and manufacturomes are incomplete.     

Effects of soman inhibition and of structural differences on cholinesterase molecular dynamics: a neutron scattering study

Incoherent elastic neutron scattering experiments on members of the cholinesterase family were carried out to investigate how molecular dynamics is affected by covalent inhibitor binding and by differences in primary and quaternary structure. Tetrameric native and soman-inhibited human butyrylcholinesterase (HuBChE) as well as native dimeric Drosophila melanogaster acetylcholinesterase (DmAChE) hydrated protein powders were examined. Atomic mean-square displacements (MSDs) were found to be identical for native HuBChE and for DmAChE in the whole temperature range examined, leading to the conclusion that differences in activity and substrate specificity are not reflected by a global modification of subnanosecond molecular dynamics. MSDs of native and soman-inhibited HuBChE were identical below the thermal denaturation temperature of the native enzyme, indicating a common mean free-energy surface. Denaturation of the native enzyme is reflected by a relative increase of MSDs consistent with entropic stabilization of the unfolded state. The results suggest that the stabilization of HuBChE phosphorylated by soman is due to an increase in free energy of the unfolded state due to a decrease in entropy.     

Small-angle neutron scattering by a strongly denatured protein: analysis using random polymer theory.

Small-angle neutron scattering profiles are presented from phosphoglycerate kinase, in the native form and strongly denatured in 4 M guanidinium chloride (GdnHCl) solution. The data are interpreted using a model in which the excess scattering density associated with the protein is represented as a finite freely jointed chain of spheres. The similarity of the model-derived scattering function to experiment increases asymptotically with the number of spheres. The improvement of the fit obtained with more than approximately 200 spheres (i.e., two residues per sphere) is insignificant. The effects of finite size of the scattering units and of scattering length variation along the polypeptide chain are examined. Improved agreement with experiment is obtained when these effects are taken into account. A method for rapid calculation of the scattering profile of a full, all-atom configuration is examined. It is found that a representation of the chain containing two scattering units per residue, placed at the backbone and side-chain scattering length centroids, reproduces the full, all-atom profile to within 2%.     

Protein flexibility from the dynamical transition: a force constant analysis

A standard analysis of the scattered neutron incoherent elastic intensity measured with very good energy resolution yields elastic scans, i.e., mean-square displacements of atomic motions (in a pico to nanosecond time scale) in a sample as a function of temperature. This provides a quick way for characterizing the dynamical behavior of biological macromolecules, such behavior being correlated with biological function and activity. Elastic scans of proteins exhibit a dynamical transition at approximately 200 K, marking a cross-over in molecular fluctuations between harmonic and nonharmonic dynamical regimes. This paper presents an approach allowing analysis of the elastic scan in terms of force constants and related parameters, such as the free energy barrier DeltaG at the transition. We find that the increased protein flexibility beyond the dynamical transition is associated with DeltaG approximately equals RT and effective force constants of the order of 0.1-3 N/m. The analysis provides a set of parameters for characterizing molecular resilience and exploring relations among dynamics, function, and activity in proteins.     

Picosecond dynamics of T and R forms of aspartate transcarbamylase: a neutron scattering study

E. coli aspartate transcarbamylase (ATCase) is a 310 kDa allosteric enzyme which catalyses the first committed step in pyrimidine biosynthesis. The binding of its substrates, carbamylphosphate and aspartate, induces significant conformational changes. This enzyme shows homotropic cooperative interactions between the catalytic sites for the binding of aspartate. This property is explained by a quaternary structure transition from T state (aspartate low affinity) to R state (aspartate high affinity) accompanied by a 5% increase of radius of gyration of ATCase. The same quaternary structure change is observed upon binding of the bisubstrate analogue PALA (N-(phosphonacetyl)-L-aspartate. Owing to the large incoherent neutron scattering cross-section of the hydrogen atom and the abundance of this element in proteins, inelastic neutron scattering gives a global view of protein dynamics as sensed via the individual motions of its hydrogen atoms. We present neutron scattering results of the local dynamics (few angstroms), at short time (few tens of picoseconds), of ATCase in T and R forms. Compared to the T form, we observe an increased mobility of the protein in the R form that we associate to an increase of accessible surface area to the solvent. Beyond this specific result, this highlights the key role of the accessible surface area (ASA) in dynamic contribution to inelastic neutron data in the picosecond time scale. In particular, we want to stress out (i) that a difference at the picosecond time scale does not allow to conclude to a difference in the dynamics at a longer time scale and to address whether the T state is looser than the R state (ii) how challenging is, any comparison in terms of general dynamics (tense or relaxed) between dynamic values deduced from experimental neutron data on proteins with different sequences and therefore ASA. This caveat holds particularly when comparing dynamics of a mesophile with the corresponding extremophile.     

Temperature dependence of protein dynamics: computer simulation analysis of neutron scattering properties

The temperature dependence of the internal dynamics of an isolated protein, bovine pancreatic trypsin inhibitor, is examined using normal mode analysis and molecular dynamics (MD) simulation. It is found that the protein exhibits marked anharmonic dynamics at temperatures of approximately 100-120 K, as evidenced by departure of the MD-derived average mean square displacement from that of the harmonic model. This activation of anharmonic dynamics is at lower temperatures than previously detected in proteins and is found in the absence of solvent molecules. The simulation data are also used to investigate neutron scattering properties. The effects are determined of instrumental energy resolution and of approximations commonly used to extract mean square displacement data from elastic scattering experiments. Both the presence of a distribution of mean square displacements in the protein and the use of the Gaussian approximation to the dynamic structure factor lead to quantified underestimation of the mean square displacement obtained.     

Dynamics of well-folded and natively disordered proteins in solution: a time-of-flight neutron scattering study

Casein proteins belong to the class of natively disordered proteins. The existence of disordered biologically active proteins questions the assumption that a well-folded structure is required for function. A hypothesis generally put forward is that the unstructured nature of these proteins results from the functional need of a higher flexibility. This interplay between structure and dynamics was investigated in a series of time-of-flight neutron scattering experiments, performed on casein proteins, as well as on three well-folded proteins with distinct secondary structures, namely, myoglobin (alpha), lysozyme (alpha/beta) and concanavalin A (beta). To illustrate the subtraction of the solvent contribution from the scattering spectra, we used the dynamic susceptibility spectra emphasizing the high frequency part of the spectrum, where the solvent dominates. The quality of the procedure is checked by comparing the corrected spectra to those of the dry and hydrated protein with negligible solvent contamination. Results of spectra analysis reveal differences in motional amplitudes of well-folded proteins, where beta-sheet structures appear to be more rigid than a cluster of alpha-helices. The disordered caseins display the largest conformational displacements. Moreover their global diffusion rates deviate from the expected dependence, suggesting further large-scale conformational motions.