ABSTRACT 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 (<~200 K) or on shorter timescales (~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 (~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 ~30 and 60% at 300 K on the 10- and 100-ps timescales, respectively.
INTRODUCTION
Neutron scattering is a major technique for the examination of picosecond dynamics in biological systems. However, due to the complexity of the motions present and the large number of disparate scattering nuclei, direct interpretation of experimental results with analytical models can be fraught with difficulty. Molecular dynamics (MD) simulation can be used to overcome this problem. The dynamic structure factor can be calculated from the simulation trajectories, compared with experiment, and used to decompose the contributions to the measured scattering (Smith, 1991, 2000; Kneller and Smith, 1994; Souaille et al., 1996; Morelon et al., 1998; Tarek and Tobias, 2000; Hayward and Smith, 2002).
In recent years neutron scattering and other techniques have been used extensively in research stimulated by the finding that the temperature-dependence of the average mean-square displacement of protein internal motion exhibits a change in gradient at temperatures ranging from ~170-240 K (Knapp et al., 1982; Parak and Knapp, 1984; Doster et al., 1989; Smith, 1991, 2000; Smith et al., 1990; Rasmussen et al., 1992; Demmel et al., 1997; Fitter et al., 1997; Cordone et al., 1998; Ostermann et al., 2000). Below the transition temperature vibrational motions dominate. Above the transition temperature, the motions may involve continuous and/or jump diffusion and thus allow sampling of different conformational substates existing in different energy wells (Elber and Karplus, 1987; Karplus and Petsko, 1990; Frauenfelder et al., 1991; Kneller and Smith, 1994). It has also been observed that solvent, when present, has an important controlling effect on the dynamical transition (Ferrand et al., 1993; Demmel et al., 1997; Fitter, 1999; Fitter et al., 1997, 1998a,b; Cordone et al., 1998; Vitkup et al., 2000; Reat et al., 2000). Strong timescale dependence has also been found both theoretically (Hayward and Smith, 2002), and in elastic neutron scattering experiments of an enzyme, glutamate dehydrogenase (GDH) in CD^sub 3^OD/D^sub 2^O, 70:30, v/v cryosolvent (Daniel et al., 1998, 1999). Some studies have found a relationship between the onset of anharmonic motion and activity of proteins (Parak et al., 1980; Rasmussen et al., 1992; Ferrand et al., 1993; Ding et al., 1994; Ostermann et al., 2000).
Much of the dynamical transition research has involved the use of elastic neutron scattering to derive the average mean-square displacement of the system. In a previous article, MD simulations of an isolated protein, bovine pancreatic trypsin inhibitor (BPTI), were performed as a function of temperature and the simulation data were used to evaluate approximations used in deriving mean-square displacements from experiment (Hayward and Smith, 2002). Here we extend the MD analysis to BPTI in solution. The solvent chosen is CD^sub 3^OD/D^sub 2^O, 70:30, v/v, the same as was used in the experiments reported in Daniel et al. (1998, 1999) and Reat et al. (2000). The primary aim of this work is to decompose the contributions to the temperature-dependent elastic scattering.
An important question concerns under what conditions the isotopic coherent and incoherent cross-sections of the sample can be used as a guide to determine the proportions of scattering from each component. Isotopic cross-sections, which are commonly used to estimate the proportions of scattering from the components of a sample, are determined by nuclear properties that are independent of the nuclear dynamics. Formally, they give the cross-section for energy-integrated scattering. They also give an approximate estimate of the relative scattering intensities for components of a system in which the dynamics of the components are similar. However, in cases where the dynamics are significantly different, deviation from this approximation is expected.
Following the previous study of elastic neutron scattering (Hayward and Smith, 2002), here mean-square displacements suitable for comparison with experiment are derived by taking the atomic trajectories, calculating the dynamic structure factor, and subjecting this to the same data treatment process as was used experimentally. A comparison with experiment is made, and the components of the elastic scattering are examined. The scattering is decomposed into solvent and protein contributions, and the protein is further decomposed into external (whole-molecule rotation and translation) diffusion and internal dynamics. The simulation analysis permits the identification of dynamical signals responsible for elastic neutron scattering intensities from a small protein in solution as the timescale and temperature probed vary.
METHODS
Molecular dynamics simulations
The model system consists of one BPTI molecule in CD^sub 3^OD/D^sub 2^O, 70:30, v/v in an orthorhombic box of size 41.0 x 37.5 x 46.5 [Angstrom]^sup 3^ (on average at 300 K) with periodic boundary conditions. BPTI has 892 atoms and four internal water molecules. There are 658 CD^sub 3^OD and 661 D^sub 2^O molecules in the solvent, providing at least three solvent shells around the protein. This is adequate for simulating a protein in a bulk solvent environment.
As in typical neutron scattering experiments, the exchangeable BPTI hydrogen atoms were replaced by deuterium, leaving 324 hydrogen atoms. This system models a solution of a protein in CD^sub 3^OD/D^sub 2^O, 70:30 v/v cryosolvent, as has been used in several neutron experiments (Daniel et al., 1998, 1999; Reat et al., 2000). However, the protein used in the above experiments is GDH. The reason for choosing BPTI for the simulations and not GDH is that due to computational limitations, only a small protein could be simulated for a time long enough to cover the experimental timescale and at the required temperatures. Therefore, the assumption is made that the difference between the scattering from GDH and BPTI in the same solvent can be neglected. That this is reasonable is suggested from comparisons of the scattering of different proteins in the same CD^sub 3^OD/D^sub 2^O, 70:30 v/v solution, which was the same to within experimental error (Daniel et al., 1998, 1999; Reat et al., 2000). Moreover, the focus of the present work is on decomposing contributions to elastic scattering, a quantity that depends largely on the average mean-square displacement and is relatively insensitive to dynamical details.
The system was simulated using CHARMM (Brooks et al., 1983) version 27 with all-atom parameter set 22 for the protein and methanol (Mackerell et al., 1998). All water molecules were represented by the TIP3P potential (Jorgensen et al., 1983). The simulations were performed in the NPT ensemble. A timestep of 0.001 ps was used with SHAKE (Ryckaert et al., 1977) applied to constrain bonds containing hydrogen or deuterium atoms. Nonbonded and electrostatic interactions were truncated using a shifting function (Steinbach and Brooks, 1994) at 13.0 [Angstrom]. The pressure of the system was set to be 1 atm.
Simulations were performed at 18 different temperatures: 80 K, 100 K, 120 K, 140 K, and 160 K, then in steps of 10-280 K, and finally 300 K. The starting structure for the first temperature simulated, 80 K, was the energy-minimized BPTI crystal structure (Parkin et al., 1996), Protein Data Bank reference (IBPI) (Berman et al., 2000), equilibrated for 500 ps in a fully equilibrated box of solvent. The starting structures for the rest of the simulations at increasing temperatures were the final structures from the preceding temperature. The systems at each temperature were equilibrated for 150 ps and then data collected every 0.1 ps for 520 ps; i.e., a total simulation time of 670 ps per temperature. The total simulation time was 12.060 ns requiring 597 CPU h running in parallel on 64 processors on an IBM SP2 computer; i.e., 38,208 CPU h. After the final simulation at 300 K the RMS deviation of the backbone heavy atoms from the energy-minimized crystal structure was 1.3 [Angstrom] and that calculated using only the secondary structural elements was 1.0 [Angstrom]. Therefore the protein structure remained stable.
Neutron scattering properties
We used the nMOLDYN program (Kneller et al., 1995) to calculate neutron scattering properties from the atomic trajectories.
Dynamic structure factor
Timescale dependence
We also examine ther the timescale dependence of the scattering. This is of particular interest because it was demonstrated, using elastic neutron scattering from instruments of different energy resolutions, that dynamical transition behavior of two enzyme/cryosolvent solutions is strongly time-dependent (Daniel et al., 1999). The time-dependent mean-square displacement, [left angle bracket]u^sup 2^^sup [is proportional to]^ (t)[right angle bracket] = [left angle bracket](R ^sub [is proportional to]^(t) - R ^sub [is proportional to]^ (0))^sup 2^[right angle bracket], were calculated from the molecular dynamics configurations as follows,
External and internal dynamics
The protein trajectory was decomposed into exernal (whole-melecule diffusion) and internal components. The internal motions were extracted by superimposing every frame from the atomic trajectory with and RMS coordinate fit onto the first frame. The result is a new atomic trajectory of internal motions. The translational mean-square displacement was calculated by creating a new trajectory of the center-of-mass of the protein. The center-of-mass was weighted by the scattering cross-sections. The mean-square displacement for the internal protein motion was subtracted from that of the total protein dynamics. The rotational components of the protein motion were separated by subtraction of the translational mean-square displacement form the external mean-square displacement.
Results
S^sub int^(q ,0) and [left angle bracket]u^sup 2^[right angle bracket] are shown in Fig. 1 from experiments on GDH in CD^sub 3^OD/D^sub 2^O, 70:30 v/v and from the model system of BPTI in the same cryosolvent. As described in Methods, for the purposes of comparison the simulation-derived dynamic structure factor has been calculated at the experimental protein concentration and processed in exactly the save way as the experimental quantity. Experiment and simulation are in close agreement up to temperatures of 200 K. The small peak in the experimental [left angle bracket]u^sup 2^[right angle bracket] at 190 K is of unknown origin-it may be due to contamination by diffraction at this temperature. Above 220 K deviation from linearity, i.e., dynamical transition behavior, is seen.
This occurs somewhat lower in temperature and more gradually in the simulation than in the experiment. The differences between experiment and simulation at these temperatures are small but significant and my be due to errors in the simulation (e.g., force field, sampling) or in the experiment (e.g., contamination from small-angle or multiple scattering) as well as differences expected between GDH and BPTI (for example, in the whole-molecule diffusional properties). At 300 K the simulation-derived and experimental [left angle bracket]u^sup 2^[right angle bracket] agree closely.
We now undertake a decomposition of the dynamical contributions to the scattering cross-sections for the protein, solvent, and total system are shown for the simulation model in Table 1. Based on these values, 60% of the total scattering originates from the solvent.
The relationship between the dynamic structure factor and the isotopic cross-section is characterized by the zeroth-order sum rule
Thus the energy-integrated scattering is given by the scattering cross-sections. However, the proportions of scattering in any given subspace of (q,w) also depend on the microscopic dynamics of the components. To estimate this effect, the relative contributions of solvent and protein to the elastic scattering intensity were calculated. This was performed using the Gaussian approximation in the limit q [arrow right] 0 with Eq. 8. Equation 8 requires the input of [left angle bracket]u^sub a^^sup 2^[right angle bracket], calculated here using Eq. 11 with t= 10 ps or 200 ps. Thus, the scattering contributions on two different timescales are evaluated. The 10-ps timescale approximates an instrumental energy resolution of 100 [mu]eV (FWHM) whereas the 200-ps timescale approximates 5 [mu]eV.
The results for [left angle bracket]u^sup 2^[right angle bracket] calculated at 10 ps are shown in Fig. 2 A. At each temperature, the sum of the protein and solvent contributions gives the intensity for the whole system. At 80 K the proportions of the scattering are given by the ration of the isotopic cross-sections in Table 1. This remains true up to 180 K and reflects the fact that at low temperatures the solvent and protein both vibrate approximately harmonically with roughly similar dynamics. Above 180 K one ovserves the interesting effect that the variation of the solvent intensity with temperature is much stronger than that of the protein, the consequence of which is that that of the solvent largely determines the change with temperature of the scattering from the solution.
The long-time results (200 ps) are shown in Fig. 2 B. Once again, the shape of the profile for the whole system is largely determined by that of the solvent. The relative contributions for T < 180 K are similar to those of Fig. 2 A and are again given by the proportions of the isotopic cross-sections. At higher temperatures diffusive motion in the solvent is
quickly activated, such that at T > 270 K the solvent scattering is outside the (q, [omega]) window examined, and the solvent intensity is zero. Diffusive motion of the protein is also activated but is more confined, remaining within the accessible (q, [omega]) window at high temperatures. Thus, at T > 270 K only the protein motion is seen.
To decompose the various contributions to the protein [left angle bracket]u^sup 2^[right angle bracket], the internal and external (whole-molecule diffusive) dynamics were determined, the latter being further subdivided into rotational and translational components. The protein [left angle bracket]u^sup 2^[right angle bracket] were calculated directly from the atomic trajectories, on timescales of 10 and 100 ps.
The results at 10 and 100 ps are shown in Fig. 3. On both timescales, at T < 200 K the internal motions make the only significant contribution to the protein (u^sup 2^). Deviation in linearity is seen in the internal dynamics in Fig. 3 B at a relatively low temperature (150 K), a phenomenon that was also seen in simulations of isolated BPTI (Hayward and Smith, 2002). There is a further change in slope of the total protein (u^sup 2^) at -200-220 K on both timescales which involves a change in the internal (u^sup 2^) together with activation of the external motions. At 10 ps and above -200 K the global diffusive and internal [left angle bracket]u^sup 2^[right angle bracket] increase with temperature at similar rates, and the external motion contributes 30% of (u^sup 2^) for T > 280 K. However, at 100 ps and above 200 K the external (u^sup 2^) increases at a much faster rate than the internal (u^sup 2^). Both (u^sup 2^) are equal in value at 260 K. At 300 K the external motions contribute 60% to the overall displacement. Rotational and translational motions contribute about equally on both timescales. Self-diffusion coefficients were estimated from the slope of the mean-square displacement curves over 40-120 ps according to the Einstein relation D = lim^sub t[arrow right][infinity]^(1/6t) [left angle bracket]u^sup 2^(t)[right angle bracket]. The 300 K rotational and translational diffusion coefficients are 3-4 X 10 7 cm^sup 2^ s^sup -1^ and are in the range of experimental values for small proteins in solution.
CONCLUSIONS
This work demonstrates the usefulness of MD simulation in understanding dynamical contributions to neutron scattering from condensed-phase biological systems. When comparing experimental with simulation-derived elastic scattering intensities it must be realized that both contain errors. For example, the simulation will contain force field and sampling errors whereas the experiment may contain inaccuracies originating from counting statistics, multiple scattering, and small-angle scattering. Moreover, the fact that in the present case the simulation and experiment were performed on different proteins may also affect the scattering, particularly as concerns the rotational and translation diffusive component. In a previous article it was shown that both the use of the Gaussian approximation and the existence of a distribution of mean-square displacements lead to quantified errors in the mean-square displacement derived from an experiment on a protein (Hayward and Smith, 2002). However, in the present comparison in Fig. 1, as the same method was used to derive the mean-square displacements from the simulation as from experiment these two sources of error will largely cancel. Nevertheless, given the various remaining potential sources of experimental and simulation error, the temperature-dependence of the experimental and simulation-derived mean-square displacement and integrated elastic intensity in Fig. 1 are in reasonable agreement.
The MD simulations here are decomposed so as to identify the dynamical components of the scattering from the simulation model. At lower temperatures (<200 K) or on shorter timescales (10 ps) the scattering contributions are given by the isotopic scattering cross-sections. This is probably due to the fact that diffusive motion is absent and the whole system vibrates approximately harmonically. In this case the mean-square displacements of the protein and solvent, and thus the elastic scattering, will be approximately the same. However, for T > 200 K, marked deviations are seen due to differences in the dynamics of the components of the solution. Above 200 K, three classes of dynamics are activated: diffusive solvent motion, external protein diffusion (rotational and translational), and internal protein anharmonic motions. The solvent diffusion rapidly increases with temperature, such that the reduction with temperature of the total elastic intensity is largely dominated by the reduction in solvent self-coherent scattering. The solvent diffusion is such that at higher temperatures (>240 K) and longer times (100 ps) the protein makes the only significant contribution to the scattering. This result provides some theoretical support for the recent intriguing report that protein dynamics can be measured using neutron scattering even in H2O solution if an instrument of appropriate energy resolution is used (Tehei et al., 2001).
Decomposition of the protein mean-square displacement shows that internal motion dominates below the dynamical transition at 200 K, whereas above it, both anharmonic internal dynamics and external rotational whole-molecule diffusion are activated at 180-210 K. The finding that external motion contributes a significant proportion of the picosecond scattering of a protein solution even with deuterated solvent is in harmony with the results of Perez et al. (1999) in which a careful decomposition of the external and internal contributions to scattering from lysozyme in aqueous solution was performed and was checked with classical measurements of diffusion constants. In the present work it is found that the proportion that the external dynamics contributes to the mean-square displacement depends on the timescale observed, with a correspondingly higher external contribution on longer timescales. This contribution will clearly also be modified by the size of the protein treated, with BPTI, which is relatively small, diffusing relatively quickly. However, a useful experimental strategy for the future would be the development of techniques to reduce the contribution of protein external diffusion (e.g., by tethering on a solid support) such that the "time-space" window in which internal protein motions dominate observed scattering can be fully utilized to observe internal, possibly functional, protein fluctuations.
We acknowledge the High-Performance Computing Center in Karlsruhe for time on their IBM RS/6000 SP-256. We thank Valerie Real and Erika Balog for assistance with the experimental data and Torsten Becker and Alex Tournier for useful discussions.
[Sidebar]
Biophysical Journal Volume 85 August 2003 679-685
[Reference]
REFERENCES
Berman, H. M., J. Westbrook, Z. Feng, G. Gilliland, T. N. Bhat, H. Weissig, I. N. Shindyalov, and P. E. Bourne. 2000. The protein data bank. Nucl. Acids Res. 28:235-242.
Brooks, B. R., R. E. Bruccoleri, B. D. Olafson, D. J. States, S. Swaminathan, and M. Karplus. 1983. CHARMM. J. Comp. Chem. 4:187-217.
Cordone, L., P. Galajda, E. Vilrano, A. Gassmann, A. Ostermann, and F. Parak. 1998. A reduction of protein specific motions in co-ligated myoglobin embedded in a lrehalose glass. Eur. Biophys. J. 27:173-176.
Daniel, R. M., J. C. Smith, M. Ferrand, S. Hery, R. Dunn, and J. L. Finney. 1998. Enzyme activity below the dynamical transition at 220K. Biophys. J. 75:2504-2507.
Daniel, R. M., J. L. Finney, V. Reat, R. Dunn, M. Ferrand, and J. C. Smith. 1999. Enzyme dynamics and activity: time-scale dependence of dynamical transitions in glutamate dehydrogenase solution. Biophys. J. 77:2184-2190.
Demmel, F., W. Doster, W. Petry, and A. Schulte. 1997. Vibrational frequency shifts as a probe of hydrogen bonds: thermal expansion and glass transition of myoglobin in mixed solvents. Eur. Biophys. J. 26:327335.
Ding, X., B. Rasmussen, G. A. Petesko, and D. Ringe. 1994. Direct structural observation of an acyl-enzyme intermediate in the hydrolysis of an ester substrate by elastase. Biochemistry. 33:9285-9293.
Doster, W., S. Cusack, and W. Petry. 1989. Dynamical transition of myoglobin revealed by inelastic neutron scattering. Nature. 337:754756.
Elber, R., and M. Karplus. 1987. Multiple conformation states of proteins: a molecular dynamics analysis of myoglobin. Science. 235:318-321.
Ferrand, M., A. J. Dianoux, W. Petry, and G. Zaccai. 1993. Thermal motions and function of bacteriorhodopsin in purple membranes: effects of temperature and hydration studied by neutron scattering. Proc. Natl. Acad. Sci. USA. 90:9668-9672.
Fitter, J., R. E. Lechner, and N. A. Dencher. 1997. Picosecond molecular motions in bacteriorhodopsin from neutron scattering. Biophys. J. 73:2126-2137.
Fitter, J., O. P. Ernst, T. Hauss, R. E. Lechner, K. P. Hoffmann, and N. A. Dencher. 1998a. Molecular motions and hydration of purple membranes and disk membranes studied by neutron scattering. Eur. Biophys. J. 27:638-645.
Fitter, J., S. A. W. Verclas, R. E. Lechner, H. Seelert, and N. A. Dencher. 1998b. Function and picosecond dynamics of bacteriorhodopsin in purple membrane at different lipidation and hydration. FEBS Lett. 433:321-325.
Fitter, J. 1999. The temperature dependence of internal molecular motions in hydrated and dry [alpha]-amylase: the role of hydration water in the dynamical transition of proteins. Biophys. J. 76:1034-1042.
Frauenfelder, H., S. G. Sligar, and P. G. Wolynes. 1991. The energy landscapes and motions of proteins. Science. 254:1598-1603.
Hayward, J. A., and J. C. Smith. 2002. Temperature dependence of protein dynamics: computer simulation analysis of neutron scattering properties. Biophys. J. 82:1216-1225.
Jorgensen, W. L., J. Chandrasekhar, J. D. Madura, R. W. Impey, and M. L. Klein. 1983. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 79:926-935.
Karplus, M., and G. A. Petsko. 1990. Molecular dynamics simulations in biology. Nature. 347:631-639.
Knapp, E. W., S. F. Fischer, and F. Parak. 1982. Protein dynamics from Mossbauer. J. Phys. Chem. 86:5042-5047.
Kneller, G. R., and J. C. Smith. 1994. Liquid-like side-chain dynamics in myoglobin. J. Mol. Biol. 242:181-185.
Kneller, G. R., V. Keiner, M. Kneller, and M. Schiller. 1995. nMOLDYN: a program package for a neutron scattering oriented analysis of molecular dynamics simulations. Comp. Phys. Comm. 91:191-214.
Lovesey, S. W. 1984. Theory of Neutron Scattering from Condensed Matter. Clarendon Press, Oxford.
Mackerell, A. D., D. Bashford, M. Bellot, J. R. Dunbrack, R. L. Evenseck, M. J. Field, S. Fischer, J. Gao, H. Guo, S. Ha, D. Joseph, L. Kuchnir, K. Kuczera, F. T. K. Lau, C. Mattos, S. Michnick, T. Ngo, D. T. Nguyen, B. Prodhom, I. W. E. Reiher, B. Roux, M. Schlenkrich, J. C. Smith, R. Stote, J. Straub, M. Watanabe, J. Wiorkiewicz-Kuczera, J. Yin, and M. Karplus. 1998. All-atom empirical potential for molecular modeling and dynamics studies of proteins. J. Phys. Chem. B. 102:3586-3616.
Morelon, N.-D., G. Kneller, M. Ferrand, A. Grand, J. C. Smith, and M. Bee. 1998. Dynamics of an alkane chain included in an organic matrix. J. Chem. Phys. 109:2883-2904.
Ostermann, A., R. Waschipky, F. G. Parak, and G. U. Nienhaus. 2000. Ligand binding and conformational motions in myoglobin. Nature. 404:205-208.
Parak, F., and E. W. Knapp. 1984. A consistent picture of protein dynamics. Proc. Natl. Acad. Sci. USA. 81:7088-7092.
Parak, F., E. N. Frolov, A. A. Kononenko, R. L. Mossbauer, V. I. Goldanskii, and A. B. Rubin. 1980. Evidence for a correlation between the photoinduced electron transfer and dynamic properties of the chromatophore membranes from Rhodospirillum rubrum. FEBS Lett. 117:368-372.
Parkin, S., B. Rupp, and H. Hope. 1996. Structure of bovine pancreatic trypsin inhibitor at 125K definition of carboxyl-terminal residues Gly-57 and Ala-58. Crystallogr. D. Biol. Crystallogr. 52:18-29.
Perez, J., J.-M. Zanotti, and D. Durand. 1999. Evolution of the internal dynamics of two globular proteins from dry powder to solution. Biophys. J. 77:454-469.
Rahman, A., K. S. Singwi, and A. Sjolander. 1962. Theory of slow neutron scattering by liquids. I. Phys. Rev. 126:986-996.
Rasmussen, B. F., A. M. Stock, D. Ringe, and G. A. Petsko. 1992. Crystalline ribonuclease A loses function below the dynamical transition at 220K. Nature. 357:423-424.
Reat, V., R. Dunn, M. Ferrand, J. L. Finney, R. M. Daniel, and J. C. Smith. 2000. Solvent dependence of dynamic transitions in protein solutions. Proc. Natl. Acad. Sci. USA. 97:9961-9966.
Ryckaert, J. P., G. Ciccotti, and H. J. C. Berendsen. 1977. Numerical integration of the Cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes (SHAKE). J. Comp. Phys. 23:327-341.
Smith, J. C., K. Kuczera, and M. Karplus. 1990. Temperature-dependence of myoglobin dynamics: neutron spectra calculated from molecular dynamics simulations of myoglobin. Proc. Natl. Acad. Sci. USA. 87:1601-1605.
Smith, J. C. 1991. Protein dynamics: comparison of simulations with inelastic neutron scattering experiments. Q. Rev. Biophys. 24:227-291.
Smith, J. C. 2000. Inelastic and quasielastic neutron scattering: complementarity with biomolecular simulation. In Structure and Dynamics of Biomolecules: Neutron Synchroton Radiation for Condensed Matter Studies. Oxford University Press, Oxford, pp. 161-180.
Souaille, M., F. Guillaume, and J. C. Smith. 1996. Dynamics of an alkane in a urea host. J. Chem. Phys. 105:1529-1546.
Steinbach, P. J., and B. R. Brooks. 1994. New spherical-cutoff methods for long-range forces in macromolecular simulation. J. Comp. Chem. 15:667-683.
Tarek, M., and D. J. Tobias. 2000. Enzyme dynamics and activity: time-scale dependence of dynamical transitions in glutamate dehydrogenase solution. Biophys. J. 79:3244-3257.
Tehei, M., D. Madern, C. Pfister, and G. Zaccai. 2001. Fast dynamics of malate dehydrogenase and bovine serum albumin measured by neutron scattering under various solvent conditions influencing protein stability. Proc. Natl. Acad. Sci. USA. 98:14356-14361.
Vitkup, D., D. Ringe, G. A. Petsko, and M. Karplus. 2000. Solvent mobility and the protein glass transition. Nat. Struct. Biol. 7:34-38.
[Author Affiliation]
Jennifer A. Hayward,* John L. Finney,[dagger] Roy M. Daniel,[double dagger] and Jeremy C. Smith*
*Interdisciplinary Centre for Computational Science (JWR), University of Heidelberg, Heidelberg, Germany; [dagger]Department of Physics and Astronomy, University College London, London, England; and [dagger]Department of Biological Sciences, University of Waikato, Hamilton, New Zealand
[Author Affiliation]
Submitted October 22, 2002, and accepted for publication February 14, 2003.
Address reprint requests to Jeremy C. Smith, Tel.: 49-622-154-8857; Fax: 49-622-154-8868; E-mail: biocomputing@iwr.uni-heidelberg.de.
(C) 2003 by the Biophysical Society
0006-3495/03/08/679/07 $2.00
