Effect of the environment on the protein dynamical transition: a neutron scattering study

We performed an elastic neutron scattering investigation of the molecular dynamics of lysozyme solvated in glycerol, at different water contents h (grams of water/grams of lysozyme). The marked non-Gaussian behavior of the elastic intensity was studied in a wide experimental momentum transfer range, as a function of the temperature. The internal dynamics is well described in terms of the double-well jump model. At low temperature, the protein total mean square displacements exhibit an almost linear harmonic trend irrespective of the hydration level, whereas at the temperature T(d) a clear changeover toward an anharmonic regime marks a protein dynamical transition. The decrease of T(d) from approximately 238 K to approximately 195 K as a function of h is reminiscent of that found in the glass transition temperature of aqueous solutions of glycerol, thus suggesting that the protein internal dynamics as a whole is slave to the environment properties. Both T(d) and the total mean square displacements indicate that the protein flexibility strongly rises between 0.1 and 0.2h. This hydration-dependent dynamical activation, which is similar to that of hydrated lysozyme powders, is related to the specific interplay of the protein with the surrounding water and glycerol molecules.     

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.     

Association of sugar-based carboranes with cationic liposomes: an electron spin resonance and light scattering study

The possibility of cationic (di-oleoyltrimethylammonium propane, DOTAP)/(l-α-dioleoylphosphatidyl-ethanolamine, DOPE) liposomes to act as carriers of boronated compounds such as 1,2-dicarba-closo-dodecaboran(12)-1-ylmethyl](β-d-galactopyranosyl)-(1→4)-β-d-glucopyranoside and 1,2-di-(β-d-gluco-pyranosyl-ox)methyl-1,2-dicarba-closo-dodeca-borane(12) has been investigated by Electron Spin Resonance (ESR) of n-doxyl stearic acids (n-DSA) and Quasi-Elastic Light Scattering (QELS). Both these carboranes have potential use in Boron Neutron Capture Therapy (BNCT), which is a targeted therapy for the treatment of radiation resistant tumors. They were shown to give aggregation both in plain water and in saline solution. Carborane aggregates were, however, disrupted when DOTAP/DOPE liposome solutions were used as dispersing agents. The computer analysis of the ESR spectra from carborane-loaded liposomes allowed to establish an increase of the order degree in the liposome bilayer with increasing carborane concentration, together with a decreased mobility. The same discontinuities of both correlation time and order parameter with respect to temperature variations were observed in carborane-containing and carborane-free liposomes. This suggested that a homogeneous dispersion of nitroxides and carboranes occurred in the liposome bilayer. The ESR line shape analysis proved that no dramatic changes were induced in the liposome environment by carborane insertion. QELS data showed that the overall liposome structure was preserved, with a slight decrease in the mean hydrodynamic radius and increase in polydispersity caused by the guest molecules.     

Association of sugar-based carboranes with cationic liposomes: an electron spin resonance and light scattering study

The possibility of cationic (di-oleoyltrimethylammonium propane, DOTAP)/(L-alpha-dioleoylphosphatidyl-ethanolamine, DOPE) liposomes to act as carriers of boronated compounds such as 1,2-dicarba-closo-dodecaboran(12)-1-ylmethyl](beta-D-galactopyranosyl)-(1-->4)-beta-D-glucopyranoside and 1,2-di-(beta-D-gluco-pyranosyl-ox)methyl-1,2-dicarba-closo-dodeca-borane(12) has been investigated by Electron Spin Resonance (ESR) of n-doxyl stearic acids (n-DSA) and Quasi-Elastic Light Scattering (QELS). Both these carboranes have potential use in Boron Neutron Capture Therapy (BNCT), which is a targeted therapy for the treatment of radiation resistant tumors. They were shown to give aggregation both in plain water and in saline solution. Carborane aggregates were, however, disrupted when DOTAP/DOPE liposome solutions were used as dispersing agents. The computer analysis of the ESR spectra from carborane-loaded liposomes allowed to establish an increase of the order degree in the liposome bilayer with increasing carborane concentration, together with a decreased mobility. The same discontinuities of both correlation time and order parameter with respect to temperature variations were observed in carborane-containing and carborane-free liposomes. This suggested that a homogeneous dispersion of nitroxides and carboranes occurred in the liposome bilayer. The ESR line shape analysis proved that no dramatic changes were induced in the liposome environment by carborane insertion. QELS data showed that the overall liposome structure was preserved, with a slight decrease in the mean hydrodynamic radius and increase in polydispersity caused by the guest molecules.     

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.     

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.     

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.     

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.     

Solution behavior of a sugar-based carborane for boron neutron capture therapy: a nuclear magnetic resonance investigation

The physico-chemical properties of beta-lactosyl-closo-orthocarborane in water solution were investigated by multinuclear NMR, (13)C NMR relaxation, and ab-initio calculations. This molecule represents a potentially selective boron carrier in Boron Neutron Cancer Therapy (BNCT) and exhibits amphiphilic characteristics. Its structural and dynamic features were studied comparing NMR data acquired in both aggregating and non-aggregating conditions. Aggregates are characterized by rapid exchange with the bulk and by high sensitivity to temperature conditions. An unusually stable intra-molecular CH...O hydrogen bond was found to persist in water solution both for the free molecules and after aggregate formation. At the same time, inter-molecular specific CH()O interactions do not seem effective in the aggregate formation process, which appears to take place only on non-specific hydrophobic basis.     

Deuterated iodoacetanilide: A covalent protein label for neutron scattering studies

A simple preparation of a deuterated alkylating reagent is described. Its reaction with the β⁹³-SH groups of hemoglobin was used to illustrate the preparation of a protein derivative which contains 5 covalently attached non-exchangeble deuterium atoms at a unique locus on each β chain.     

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.     

The dynamical Matryoshka model: 1. Incoherent neutron scattering functions for lipid dynamics in bilayers

Fluid lipid bilayers are the building blocks of biological membranes. Although there is a large amount of experimental data using incoherent quasi-elastic neutron scattering (QENS) techniques to study membranes, very little theoretical works have been developed to study the local dynamics of membranes. The main objective of this work is to build a theoretical framework to study and describe the local dynamics of lipids and derive analytical expressions of intermediate scattering functions (ISF) for QENS. As results, we developed the dynamical Matryoshka model which describes the local dynamics of lipid molecules in membrane layers as a nested hierarchical convolution of three motional processes: (i) individual motions described by the vibrational motions of H-atoms; (ii) internal motions including movements of the lipid backbone, head groups and tails, and (iii) molecule movements of the lipid molecule as a whole. The analytical expressions of the ISF associated with these movements are all derived. For use in analyzing the QENS experimental data, we also derived an analytical expression for the aggregate ISF of the Matryoshka model which involves an elastic term plus three inelastic terms of well-separated time scales and whose amplitudes and rates are functions of the lipid motions. And as an illustrative application, we used the aggregated ISF to analyze the experimental QENS data on a lipid sample of multilamellar bilayers of DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine). It is clear from this analysis that the dynamical Matryoshka model describes very well the experimental data and allow extracting the dynamical parameters of the studied system.     

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.     

Conformational dynamics of a multidomain protein by neutron scattering and computational analysis

The flexible conformations of a multidomain protein are responsible for its biological functions. Although MurD, a 47-kDa protein that consists of three domains, sequentially changes its domain conformation from an open form to a closed form through a semiclosed form in its enzymatic reaction, the domain dynamics in each conformation remains unclear. In this study, we verify the conformational dynamics of MurD in the corresponding three states (apo and ATP- and inhibitor-bound states) with a combination of small-angle x-ray and neutron scattering (SAXS and SANS), dynamic light scattering (DLS), neutron backscattering (NBS), neutron spin echo (NSE) spectroscopy, and molecular dynamics (MD) simulations. Applying principal component analysis of the MD trajectories, twisting and open-closed domain modes are identified as the major collective coordinates. The deviations of the experimental SAXS profiles from the theoretical calculations based on the known crystal structures become smaller in the ATP-bound state than in the apo state, and a further decrease is evident upon inhibitor binding. These results suggest that domain motions of the protein are suppressed step by step of each ligand binding. The DLS and NBS data yield collective and self-translational diffusion constants, respectively, and we used them to extract collective domain motions in nanometer and nanosecond scales from the NSE data. In the apo state, MurD shows both twisting and open-closed domain modes, whereas an ATP binding suppresses twisting domain motions, and a further reduction of open-closed mode is seen in the inhibitor-binding state. These observations are consistent with the structure modifications measured by the small-angle scattering as well as the MD simulations. Such changes in the domain dynamics associated with the sequential enzymatic reactions should be related to the affinity and reaction efficiency with a ligand that binds specifically to each reaction state.     

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%.     

The location of DNA in complexes of recA protein with double-stranded DNA. A neutron scattering study

Purified recA protein is found as rodlike homopolymers, and it forms filamentous complexes with double-stranded DNA that are stable in the presence of ATP gamma S, a nonhydrolyzable analogue of ATP. The structure of these filaments has been described in some detail by electron microscopy. Here we confirm the mass per length of 6.5 recA/100 A in solution by small-angle neutron scattering and extend the analysis to homopolymers of recA protein, finding a mass per length of about 7 recA/100 A and a radial mass distribution (cross-sectional radius of gyration) significantly different for the two filaments. The models proposed so far for the structure of the complex have placed the DNA in the center of the filament. Here we verify this assumption using small-angle neutron scattering to locate the DNA in the complexes, exploiting the contrast variation method in D2O/H2O mixtures. Model calculations show that the natural contrast difference between DNA and protein is not sufficient to locate the DNA (which accounts for only 4.7% of the mass in the complex). When deuterated DNA is used, the contrast difference is enhanced, and model calculations and experiment then converge, indicating that the DNA is indeed near the axis of the complex.     

Sequence-dependent dynamical instability of the human prion protein: a comparative simulation study

The present study aimed to explore the most probable regions of the human prion protein backbone for which the initial steps of conformational transitions as a result of intrinsic and extrinsic perturbing factors on the protein structure can be assigned. A total of 0.3-μs molecular dynamics simulations on several analog structures of the protein have been performed. To mimic the impact of the extrinsic and intrinsic destructive parameters on the dynamical characteristics of the protein, mild acidic conditions and R208H mutation have been simulated. The findings indicated that distribution of conformational flexibilities along the protein chain was almost independent of the induced perturbing factors, and was mostly centralized on certain distinct parts of the structure comprising residues 132–145 and 187–203. Analyses also revealed that the segment comprising residues 187–203 may be considered as a peptide sequence, possessing high potential to start the initial steps of conformational rearrangements due to the induced physicochemical alterations. Sequence alignment and molecular dynamics data also revealed that segment 178–203 prefers to accommodate in extended structures rather than α-helices. Region 178–203 may be considered as a peptide switch capable of initiating the conformational transitions due to the introduced modifications and perturbing parameters.     

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.