The Length and Viscosity Dependence of End-to-End Collision Rates in Single-Stranded DNA

Intramolecular dynamics play an essential role in the folding and function of biomolecules and, increasingly, in the operation of many biomimetic technologies. Thus motivated we have employed both experiment and simulation to characterize the end-to-end collision dynamics of unstructured, single-stranded DNAs ranging from 6 to 26 bases. We find that, because of the size and flexibility of the optical reporters employed experimentally, end-to-end collision dynamics exhibit little length dependence at length scales <11 bases. For longer constructs, however, the end-to-end collision rate exhibits a power-law relationship to polymer length with an exponent of −3.49 ± 0.13. This represents a significantly stronger length dependence than observed experimentally for unstructured polypeptides or predicted by polymer scaling arguments. Simulations indicate, however, that the larger exponent stems from electrostatic effects that become important over the rather short length scale of these highly charged polymers. Finally, we have found that the end-to-end collision rate also depends linearly on solvent viscosity, with an experimentally significant, nonzero intercept (the extrapolated rate at zero viscosity) that is independent of chain length—an observation that sheds new light on the origins of the “internal friction” observed in the dynamics of many polymer systems.     

Molecular-Dynamics Simulation of the Surface Diffusion of n-Alkanes on Pt(111)

Abstract

In the present study, the equilibrium adsorption and the dynamics of surface diffusion in a model of ethane and n-butane on a Pt(111) surface were simulated with molecular dynamics. At low temperatures, we found that both admolecules adsorb in a specific binding site. Through analysis of the trajectories, several features of the dynamics were resolved. At low temperature, we observed that diffusion occurs through a nearest-neighbor hopping mechanism involving both lateral rotation and axial translation. At high temperatures, the admolecule makes multiple-site hops and nonlocalized long flights. The temperature dependence of the diffusion coefficients was analyzed and was found to exhibit good Arrhenius behavior. The apparent diffusion coefficients follow trends seen in related experimental studies. In the case of ethane, a comparison between the diffusion barrier measured in the molecular-dynamics simulations and the theoretical barrier predicted by transition-state theory indicates that the simulated barrier is larger than the theoretical value. This finding is consistent with conclusions in recent studies of metal-atom diffusion on metal surfaces, where it was found that systems with low corrugation exhibit a non-unique relationship between the dynamical diffusion barrier and the potential-energy-surface topology.     

Calorimetric and spectroscopic studies of lipid thermotropic phase behavior in liver inner mitochondrial membranes from a mammalian hibernator

Arrhenius plots of various enzyme and transport systems associated with the liver mitochondrial inner membranes of ground squirrels exhibit changes in slope at temperatures of 20-25 degrees C in nonhibernating but not in hibernating animals. It has been proposed that the Arrhenius breaks observed in nonhibernating animals are the result of a gel to liquid-crystalline phase transition of the mitochondrial membrane lipids, which also occurs at 20-25 degrees C, and that the absence of such breaks in hibernating animals is due to a major depression of this lipid phase transition to temperatures below 4 degrees C. In order to test this hypothesis, we have examined the thermotropic phase behavior of liver inner mitochondrial membranes from hibernating and nonhibernating Richardson's ground squirrels, Spermophilus richardsonii, by differential scanning calorimetry and by 19F nuclear magnetic resonance and fluorescence polarization spectroscopy. Each of these techniques indicates that no lipid phase transition occurs in the membranes of either hibernating or nonhibernating ground squirrels within the physiological temperature range of this animal (4-37 degrees C). Moreover, differential scanning calorimetric measurements indicate that only a small depression of the lipid gel to liquid-crystalline phase transition, which is centered at about -5 degrees C in nonhibernating animals and at about -9 degrees C in hibernators, occurs. We thus conclude that the Arrhenius plot breaks observed in some membrane-associated enzymatic and transport activities of nonhibernating animals are not the result of a lipid phase transition and that a major shift in the gel to liquid-crystalline lipid phase transition temperature is not responsible for seasonal changes in the thermal behavior of these inner mitochondrial membrane proteins.     

Temperature dependence of rotational dynamics of protein and lipid in sarcoplasmic reticulum membranes

We have investigated the relationship between function and molecular dynamics of both the lipid and the Ca-ATPase protein in sarcoplasmic reticulum (SR), using temperature as a means of altering both activity and rotational dynamics. Conventional and saturation-transfer electron paramagnetic resonance (EPR) was used to probe rotational motions of spin-labels attached either to fatty acid hydrocarbon chains or to the Ca-ATPase sulfhydryl groups in SR. EPR studies were also performed on aqueous dispersions of extracted SR lipids, in order to study intrinsic lipid properties independent of the protein. While an Arrhenius plot of the Ca-ATPase activity exhibits a clear change in slope at 20 degrees C, Arrhenius plots of lipid hydrocarbon chain mobility are linear, indicating that an abrupt thermotropic change in the lipid hydrocarbon phase is not responsible for the Arrhenius break in enzymatic activity. The presence of protein was found to decrease the average hydrocarbon chain mobility, but linear Arrhenius plots were observed both in the intact SR and in extracted lipids. Lipid EPR spectra were analyzed by procedures that prevent the production of artifactual breaks in the Arrhenius plots. Similarly, using sample preparations and spectral analysis methods that minimize the temperature-dependent contribution of local probe mobility to the spectra of spin-labeled Ca-ATPase, we find that Arrhenius plots of overall protein rotational mobility also exhibit no change in slope. The activation energy for protein mobility is the same as that of ATPase activity above 20 degrees C; we discuss the possibility that overall protein mobility may be essential to the rate-limiting step above 20 degrees C.     

Langevin dynamics of peptides: the frictional dependence of isomerization rates of N-acetylalanyl-N'-methylamide.

The rate constant for the transition between the equatorial and axial conformations of N-acetylalanyl-N'-methylamide has been determined from Langevin dynamics (LD) simulations with no explicit solvent. The isomerization rate is maximum at collision frequency gamma = 2 ps-1, shows diffusive character for gamma greater than or equal to 10 ps-1, but does not approach zero even at gamma = 0.01 ps-1. This behavior differs from that found for a one-dimensional bistable potential and indicates that both collisional energy transfer with solvent and vibrational energy transfer between internal modes are important in the dynamics of barrier crossing for this system. It is suggested that conformational searches of peptides be carried out using LD with a collision frequency that maximizes the isomerization rate (i.e., gamma approximately 2 ps-1). This method is expected to be more efficient than either molecular dynamics in vacuo (which corresponds to LD with gamma = 0) or molecular dynamics in solvent (where dynamics is largely diffusive).     

Tritium exchange of spectrin versus temperature

Tritium-exchange experiments were performed to assess the dynamics of the spectrin molecule. From the kinetics presented, spectrin appears as a fast-exchanging protein compared to globular proteins. This is easily explained by its known elongated structure. Spectrin tritium exchange is also shown to be highly temperature sensitive, even in buffer conditions that greatly reduce the temperature dependence of hydroxyl catalyst concentration. At high ionic strength, the exchange rates of specific hydrogens deduced from a power-law simulation are shown to exhibit an Arrhenius behavior, with enthalpies ranging from 70 to 140 kjmol^?1. These results are discussed in relation to the known high α-helix content of spectrin. Finally, a molecular model of α-helix opening is proposed that provides satisfactory agreement with the totality of the spectrin behavior in tritium-exchange experiments.     

Temperature-dependencies of various catalytic activities of membrane-bound Na+/K+-ATPase from ox brain, ox kidney and shark rectal gland and of C12E8-solubilized shark Na+/K+-ATPase

The temperature dependence of ouabain-sensitive ATPase and phosphatase activities of membrane fragments containing the Na+/K+-ATPase were investigated in tissue from ox kidney, ox brain and from shark rectal glands. The shark enzyme was also tested in solubilized form. Arrhenius plots of the Na+/K+-ATPase activity seem to be linear up to about 20 degrees C, and non-linear above this temperature. The Arrhenius plots of mammalian enzyme (ox brain and kidney) were steeper, especially at temperatures below 20-30 degrees C, than that of shark enzyme. The Na+-ATPase activity showed a weaker temperature-dependence than the Na+/K+-ATPase activity. The phosphatase reactions measured, K+-stimulated, Na+/K+-stimulated and Na+/K+/ATP-stimulated, also showed a weaker temperature-dependence than the overall Na+/K+-ATPase activity. Among the phosphatase reactions, the largest change in slope of the Arrhenius plot was observed with the Na+/K+/ATP)-stimulated phosphatase reaction. The Arrhenius plots of the partial reactions were all non-linear. Solubilization of shark enzyme in C12E8 did not change the curvature of Arrhenius plots of the Na+/K+-ATPase activity or the K+-phosphatase activity. Since solubilization involves a disruption of the membrane and an 80% delipidation, the observed curvature of the Arrhenius plot can not be attributed to a property of the membrane as such.     

Temperature and Ca(2+)-dependence of the sarcoplasmic reticulum Ca(2+)-ATPase in haddock,salmon, rainbow trout and zebra cichlid

Temperature dependence of Ca(2+)-ATPase from the sarcoplasmic reticulum (SR) in rabbit muscle has been widely studied, and it is generally accepted that a break point in Arrhenius plot exist at approximately 20 degrees C. Whether the break point arises as a result of temperature dependent changes in the enzyme or its membrane lipid environment is still a matter of discussion. In this study we compared the temperature dependence and Ca(2+)-dependence of SR Ca(2+)-ATPase in haddock (Melanogrammus aeglefinus), salmon (Salmo salar), rainbow trout (Oncorhynchus mykiss) and zebra cichlid (Cichlasoma nigrofasciatum). The Arrhenius plot of zebra cichlid showed a break point at 20 degrees C, and the haddock Arrhenius plot was non-linear with pronounced changes in slope in the temperature area, 6-14 degrees C. In Arrhenius plot from both salmon and rainbow trout a plateau exists with an almost constant SR Ca(2+)-ATPase activity. The temperature range of the plateau was 14-21 and 18-25 degrees C in salmon and rainbow trout, respectively. Ca(2+)-dependence in the four different fish species investigated was very similar with half maximal activation (K(0.5)) between 0.2 and 0.6 micro M and half maximal inhibition (I(0.5)) between 60 and 250 micro M. Results indicated that interaction between SR Ca(2+)-ATPase and its lipid environment may play an important role for the different Arrhenius plot of the different types of fish species investigated.     

Structural diversity and extreme stability of unimolecular Oxytricha nova telomeric G-quadruplex

Oxytricha nova telomeric DNA contains guanine-rich short-tandem repeat sequences (GGGGTTTT) n and terminates as a single strand at the 3'-end. This single-stranded overhang forms a novel DNA structure, namely, G-quadruplex, comprising four quartets. In this study, we investigated the structures and dynamics of unimolecular Oxytricha nova ( O. nova) telomeric G-quadruplexes by performing single molecule fluorescence resonance energy transfer (FRET) spectroscopy and bulk circular dichroism (CD) measurements. We observed that unimolecular O. nova G-quadruplexes exhibit structural polymorphism according to monovalent cations. In the presence of Na (+), only antiparallel conformation is detected, which was demonstrated in previous studies; however, in the presence of K (+), they fold into two different conformations, a parallel conformation and an antiparallel one different from that induced by Na (+). Furthermore, these G-quadruplexes show extremely high stability in their dynamics when compared with human G-quadruplexes. While human telomeric G-quadruplexes that possess three quartets display fast dynamic behavior (<100 s) at low K (+) concentrations or high temperatures, O. nova G-quadruplexes maintain their conformational state for a long time (>1000 s), even at the lowest K (+) concentration and the highest temperature investigated. This high stability is primarily due to an extra quartet that results in additional cation coordination. In addition to cation coordination, we propose that other factors such as base stacking and the size of the thymine loop may contribute to the stability of O. nova G-quadruplexes; this is based on the fact that the O. nova G-quadruplexes were observed to be more stable than the human ones in the presence of Li (+), which is known to greatly destabilize G-quadruplexes because of imprecise coordination. This extreme stability of four-quartet G-quadruplexes enables telomere protection even in the absence of protective proteins or in the case of abrupt environmental changes, although only a single G-quadruplex structure can be derived from the short single-stranded overhang.     

Temperature dependence of the conductivity of individual potential-dependent K+-channels in mollusk neurons]

Using the patch-clamp method temperature dependences of the chord conductance of single potential--dependent slow and fast K+ channels in mollusk neurons were studied. Under control conditions (20 degrees C, 0 mV, [K+]o = 1.5 mM and [K+]i = 100 mM) the conductances of the fast and slow K+ channels were equal to 20-25 pS and 30-40 pS, respectively. Besides, the temperature dependences of the currents through the K+ channels of lesser conductance (5-20 pS) were studied. Some of these channels may be regarded as subtypes of the fast and slow K+ channels named above. It was found that for the channels of all types single channel currents arise with temperature. However, in the range of 10-20 degrees C an anomalous conductance decrease at temperature elevation was observed. For all channels except for the fast one at temperatures above 20 degrees C activation energy (delta Ea) calculated from the Arrhenius plots of the currents was about 4 kcal/mol. At the temperatures below 10 degrees C delta Ea was equal to about 12 kcal/mol. In this temperature range delta Ea had a pronounced potential dependency. Temperature dependences of the fast K+ channel conductance were opposite to those of the slow K+ channel to some extent.     

Enzyme activity below the dynamical transition at 220 K.

Enzyme activity requires the activation of anharmonic motions, such as jumps between potential energy wells. However, in general, the forms and time scales of the functionally important anharmonic dynamics coupled to motion along the reaction coordinate remain to be determined. In particular, the question arises whether the temperature-dependent dynamical transition from harmonic to anharmonic motion in proteins, which has been observed experimentally and using molecular dynamics simulation, involves the activation of motions required for enzyme function. Here we present parallel measurements of the activity and dynamics of a cryosolution of glutamate dehydrogenase as a function of temperature. The dynamical atomic fluctuations faster than approximately 100 ps were determined using neutron scattering. The results show that the enzyme remains active below the dynamical transition observed at approximately 220 K, i.e., at temperatures where no anharmonic motion is detected. Furthermore, the activity shows no significant deviation from Arrhenius behavior down to 190 K. The results indicate that the observed transition in the enzyme's dynamics is decoupled from the rate-limiting step along the reaction coordinate.     

Dielectric and calorimetric studies of hydrated purple membrane

Purple membranes (PM) from halobacteria were hydrated to approximately 0.4 and approximately 0.2 g H(2)O/g of PM and studied by dielectric spectroscopy and differential scanning calorimetry between 120 and 300 K. The dielectric process, attributed to a local (beta) relaxation of the confined supercooled water, shows an Arrhenius temperature behavior at low temperatures. In the case of the most hydrated PM a small deviation from the Arrhenius behavior occurs at 190-200 K together with a pronounced endothermic process and an increased activation energy. The observed crossover is accompanied by a reduction of the interlayer spacing due to the partial loss of the intermembrane water. All these effects at approximately 200 K are consistent with a scenario where the local relaxation process merges with a nonobservable alpha-relaxation of the interlayer water, giving rise to a more liquid-like behavior of the interfacial water. For the less hydrated sample the effects are less pronounced and shift to a slightly higher temperature.     

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.     

ON THE TEMPERATURE CHARACTERISTICS FOR FREQUENCY OF BREATHING MOVEMENTS IN INBRED STRAINS OF MICE AND IN THEIR HYBRID OFFSPRING. I

Young mice of a selected line of the dilute brown strain of mice exhibit over the range 15–25°C. (body temperature) a relation of frequency of breathing movements to temperature such that when fitted by the Arrhenius equation the data give a value for the constant µ of 24,000± calories or, less frequently, 28,000±. Young mice of an inbred albino strain show over the range 15–20°C. a value of µ = 34,000±, or, less frequently, 14,000±, with a critical temperature at about 20°C. and a value of µ = 14,000± above 20°C. The F₁ hybrids of these two strains, and the backcross generations to either parent strain, exhibit only those four values of the temperature characteristic observed in the parent strains and none other. One may therefore speak of the inheritance of the value of the constant µ, but the inheritance shows in this instance no Mendelian behavior. Furthermore there appears to be inherited the occurrence (or absence) of a critical temperature at 20°C. These experiments indicate the "biological reality" of the temperature characteristics.     

Complexity of aromatic ring-flip motions in proteins: Y97 ring dynamics in cytochrome c observed by cross-relaxation suppressed exchange NMR spectroscopy

Dynamics of large-amplitude conformational motions in proteins are complex and less understood, although these processes are intimately associated with structure, folding, stability, and function of proteins. Here, we use a large set of spectra obtained by cross-relaxation suppressed exchange NMR spectroscopy (EXSY) to study the 180° flipping motion of the Y97 ring of horse ferricytochrome c as a function of near-physiological temperature in the 288–308 K range. With rising temperature, the ring-flip rate constant makes a continuous transition from Arrhenius to anti-Arrhenius behavior through a narrow Arrhenius-like zone. This behavior is seen not only for the native state of the protein, but also for native-like states generated by adding subdenaturing amounts of guanidine deuterochloride (GdnDCl). Moderately destabilizing concentrations of the denaturant (1.5 M GdnDCl) completely removes the Arrhenius-like feature from the temperature window employed. The Arrhenius to anti-Arrhenius transition can be explained by the heat capacity model where temperature strengthens ground state interactions, perhaps hydrophobic in nature. The effect of the denaturant may appear to arise from direct protein-denaturant interactions that are structure-stabilizing under subdenaturing conditions. The temperature distribution of rate constants under different stability conditions also suggests that the prefactor in Arrhenius-like relations is temperature dependent. Although the use of the transition state theory (TST) offers several challenges associated with data interpretation, the present results and a consideration of others published earlier provide evidence for complexity of ring-flip dynamics in proteins.     

Environmentally coupled hydrogen tunneling. Linking catalysis to dynamics.

Many biological C-H activation reactions exhibit nonclassical kinetic isotope effects (KIEs). These nonclassical KIEs are too large (kH/kD > 7) and/or exhibit unusual temperature dependence such that the Arrhenius prefactor KIEs (AH/AD) fall outside of the semiclassical range near unity. The focus of this minireview is to discuss such KIEs within the context of the environmentally coupled hydrogen tunneling model. Full tunneling models of hydrogen transfer assume that protein or solvent fluctuations generate a reactive configuration along the classical, heavy-atom coordinate, from which the hydrogen transfers via nuclear tunneling. Environmentally coupled tunneling also invokes an environmental vibration (gating) that modulates the tunneling barrier, leading to a temperature-dependent KIE. These properties directly link enzyme fluctuations to the reaction coordinate for hydrogen transfer, making a quantum view of hydrogen transfer necessarily a dynamic view of catalysis. The environmentally coupled hydrogen tunneling model leads to a range of magnitudes of KIEs, which reflect the tunneling barrier, and a range of AH/AD values, which reflect the extent to which gating modulates hydrogen transfer. Gating is the primary determinant of the temperature dependence of the KIE within this model, providing insight into the importance of this motion in modulating the reaction coordinate. The potential use of variable temperature KIEs as a direct probe of coupling between environmental dynamics and the reaction coordinate is described. The evolution from application of a tunneling correction to a full tunneling model in enzymatic H transfer reactions is discussed in the context of a thermophilic alcohol dehydrogenase and soybean lipoxygenase-1.     

Conductance behavior of symmetrical tetraalkylammonium halides in aqueous sucrose solutions

The conductance behavior of some tetraalkylammonium halides (R₄NX) in saturated, aqueous solutions of sucrose has been investigated, and data on the conductance of these salts in water saturated by sucrose at 50° are reported at several temperatures within the range 25 to 70°. In these homogeneous, ternary systems, plots of —log K versus 1/T show a break at the saturation temperature, where two straight lines intersect one another. Divergence of the pair of straight lines has been found to decrease with increasing chain-length of the R₄N⁺ ions, in contrast to the structural behavior of common, alkali-metal ions. The results are interpreted in terms of the hydrophobic nature of the tetraalkylammonium halides, as well as the salting-in behavior of these salts towards sucrose molecules.     

H and C NMR studies on the temperature-dependent water and protein dynamics in hydrated elastin,myoglobin and collagen

²H NMR spin-lattice relaxation and line-shape analyses are performed to study the temperature-dependent dynamics of water in the hydration shells of myoglobin, elastin, and collagen. The results show that the dynamical behaviors of the hydration waters are similar for these proteins when using comparable hydration levels of h = 0.25–0.43. Since water dynamics is characterized by strongly nonexponential correlation functions, we use a Cole–Cole spectral density for spin-lattice relaxation analysis, leading to correlation times, which are in nice agreement with results for the main dielectric relaxation process observed for various proteins in the literature. The temperature dependence can roughly be described by an Arrhenius law, with the possibility of a weak crossover in the vicinity of 220 K. Near ambient temperatures, the results substantially depend on the exact shape of the spectral density so that deviations from an Arrhenius behavior cannot be excluded in the high-temperature regime. However, for the studied proteins, the data give no evidence for the existence of a sharp fragile-to-strong transition reported for lysozyme at about 220 K. Line-shape analysis reveals that the mechanism for the rotational motion of hydration waters changes in the vicinity of 220 K. For myoglobin, we observe an isotropic motion at high temperatures and an anisotropic large-amplitude motion at low temperatures. Both mechanisms coexist in the vicinity of 220 K. ¹³C CP MAS spectra show that hydration results in enhanced elastin dynamics at ambient temperatures, where the enhancement varies among different amino acids. Upon cooling, the enhanced mobility decreases. Comparison of ²H and ¹³C NMR data reveals that the observed protein dynamics is slower than the water dynamics.     

Temperature dependence of dynamics of hydrated myoglobin. Comparison of force field calculations with neutron scattering data

Molecular dynamics is used to probe the atomic motions of the carboxy-myoglobin protein as a function of temperature. Simulations of 150 picoseconds in length are carried out on the protein at 20, 60, 100, 180, 220, 240, 260, 280, 300, 320 and 340 K. The simulations attempt to mimic neutron scattering experiments very closely by including a partial hydration shell around the protein. Theoretical elastic, quasielastic and inelastic neutron scattering data are derived from the trajectories and directly compared with experiment. Compared to experiment, the simulation-derived elastic scattering curves show a decrease in intensity as a function of the scattering wavevector, q2. The inelastic and quasielastic spectra show that the inelastic peak is shifted to lower frequency than the experimental value, while quasielastic behavior is in good agreement with experiment. This suggests that the theoretical model is too flexible in the harmonic limit (low temperature), but accurately reproduces high-temperature behavior. Time correlation functions of the intermediate scattering function are determined. At low temperature there is one fast decay process, and at high temperatures there is an additional slow relaxation process that is due to quasielastic scattering. The average atomic fluctuations show that the protein behaves harmonically at low temperatures. At approximately 210 K, a glass-like transition in atomic fluctuations is seen. Above the transition temperature, the atomic fluctuations exhibit both harmonic and anharmonic behavior. Comparison of protein mobility behavior with experiment indicate the fluctuations derived from simulations are larger in the harmonic region. However, the anharmonic region agrees very well with experiment. The anharmonicity is large at all temperatures, with a gradual monotonic increase from 0.5 at 20 K to greater than 0.7 at 340 K without a noticeable change at the glass transition temperature. Heavy-atom dihedral transitions are monitored as a function of temperature. Trends in the type of dihedral transitions that occur with temperature are clearly visible. Dihedral transitions involving backbone atoms occur only above the glass transition temperature. The overall protein behavior results suggest that at low temperatures there is purely vibrational motion with one fast decay process, and above the glass transition temperature there is more anharmonic motion with a fast and a slower relaxation process occurring simultaneously.(ABSTRACT TRUNCATED AT 400 WORDS)