Solvation shell dynamics studied by molecular dynamics simulation in relation to the translational and rotational dynamics of supercritical water and benzene

The solvation shell dynamics of supercritical water is analyzed by molecular dynamics simulation with emphasis on its relationship to the translational and rotational dynamics. The relaxation times of the solvation number (tau S), the velocity autocorrelation function (tau D), the angular momentum correlation function (tau J), and the second-order reorientational correlation function (tau 2R) are studied at a supercritical temperature of 400 degrees C over a wide density region of 0.01-1.5 g cm(-3). The relaxation times are decomposed into those conditioned by the solvation number n, and the effect of the short-ranged structure is examined in terms of its probability Pn of occurrence. In the low to medium-density range of 0.01-0.4 g cm(-3), the time scales of water dynamics are in the following sequence: tau D>tau S approximately or > tau J approximately or > tau 2R. This means that the rotation in supercritical water is of the "in-shell" type while the translational diffusion is not. The comparison to supercritical benzene is also performed and the effect of hydrogen bonding is examined. The water diffusion is not of the in-shell type up to the ambient density of 1.0 g cm(-3), which corresponds to the absence of the transition from the collision to the Brownian picture, whereas such transition is present in the case of benzene. The absence of the transition in water comes from the fast reorganization of the hydrogen bonds and the enhanced mobility of the solvation shell in supercritical conditions.     

Dynamics and conformations of PEO chains chemically bonded on silica: comparison between 1H and 2H NMR

1H NMR was used to study the motion of monomer units in a layer of poly(ethylene oxide) chains grafted on silica. First, the dependence of the relaxation times on the grafting ratios is discussed qualitatively from a phenomenological point of view. Next, the NMR line narrowing effect by high-speed rotation is observed in the same samples with different grafting ratios. The magic angle spinning technique permits determination of two correlation times for each grafting ratio: tau(c) characteristic of an environment with a fast motion and tau(l) characteristic of an environment with a slow motion. In addition, the dynamics of these grafted chains are investigated by deuterium NMR (2H NMR), which is sensitive to the anisotropy of molecular motion. The evolution has been studied for two extreme grafting ratios and each time as a function of temperature. The anisotropy is more marked at low temperatures and for a low grafting ratio. The results are consistent with the 1H NMR relaxation times measured as a function of temperature.     

NMR studies of the dynamics of a bifunctional rhodamine probe attached to troponin C

Fluorescence polarization measurements of bifunctional rhodamine (BR) probes provide a powerful approach to determine the in situ orientation of proteins within ordered complexes such as muscle fibers. For accurate interpretation of fluorescence measurements, it is important to understand the probe dynamics relative to the protein to which it is attached. We previously determined the structure of the N-domain of chicken skeletal troponin C, BR-labeled on the C helix, in complex with the switch region of troponin I, and demonstrated that the probe does not perturb the structure or dynamics of the protein. In this study, the motion of the fluorescence label relative to the protein has been characterized using NMR relaxation measurements of 13C-labeled methyl groups on the BR probe and 15N-labeled backbone amides of the protein. Probe dynamics were monitored using off-resonance 13C-R(1rho), 13C-R(1) and {1H}-13C NOE at magnetic field strengths of 500, 600, and 800 MHz. Relaxation data were interpreted in terms of the overall rotational correlation time of the protein and a two-time scale model for internal motion of the BR methyl groups, using a numerical optimization with Monte Carlo parameter error estimation. The analysis yields a 1.5 +/- 0.4 ps correlation time for rotation around the three-fold methyl symmetry axis, and a 0.8 +/- 0.4 ns rotational correlation time for reorientation of the 13C-14N bond with an associated S2s of 0.79 +/- 0.03. Order parameters of the backbone NH vectors in the helix to which the probe is attached average S2 approximately 0.85, implying that the amplitude of independent reorientation of the BR probe is small in magnitude, consistent with results from fluorescence polarization measurements in reconstituted muscle fibers.     

Thermally activated escape rate for the Brownian motion of a fixed axis rotator in an asymmetrical double-well potential for all values of the dissipation

The Kramers theory of the escape rate of a Brownian particle from a potential well as extended by Mel'nikov and Meshkov is used to evaluate the relaxation times and the dynamic susceptibility for the rotational Brownian motion of fixed axis rotators in an asymmetric double-well potential. An expression for the escape rate valid for all values of the dissipation including the very low damping (VLD), very high damping (VHD), and crossover regimes is derived. It is shown that this expression provides a good asymptotic estimate of the inverse of the smallest nonvanishing eigenvalue lambda(1) of the underlying Fokker-Planck operator calculated by using the matrix-continued fraction method. For low barriers, where the Mel'nikov and Meshkov approach is not applicable, analytic equations for the correlation time tau( parallel) of the longitudinal dipole correlation function in the VLD and VHD limits are derived and a simple extrapolating equation valid for all values of the damping is proposed.     

The chlorobenzene-argon ground state intermolecular potential energy surface

Using the coupled cluster singles and doubles including connected triple excitations model with the augmented correlation consistent polarized valence double zeta basis set extended with a set of 3s3p2d1f1g midbond functions, we evaluate the ground state intermolecular potential energy surface of the chlorobenzene-argon van der Waals complex. The minima of 420 cm(-1) are characterized by Ar atom position vectors of the length 3.583 A, forming an angle of 9.87 degrees with respect to the axis perpendicular to the chlorobenzene plane. These results are compared to those obtained for similar complexes and to the experimental data available. From the potential the three-dimensional vibrational eigenfunctions and eigenvalues are calculated and the results allow to correct and complete the experimental assignment.     

The quantitative measurement of rotational motion of the subfragment-1 region of myosin by saturation transfer epr spectroscopy.

According to current models of muscle contraction (Huxley, H. E., Science 164: 1356-1366 [1969]), motion of flexible myosin crossbridges is essential to the contractile cycle. Using a spin-label analog of iodoacetamide bound to the subfragment No. 1 (S1) region of myosin, we have obtained rotational correlation times (tau 2) for this region of the molecule with the ultimate goal of making quantitative measurements of the motion of the crossbridges under conditions comparable to those in living, contracting muscle. We used the newly developed technique of saturation transfer electron paramagnetic resonance spectroscopy (Hyde, J.S., and Thomas, D.D., Ann. N.Y. Acad Sci. 22:680-692 [1973]), which is uniquely sensitive to rotational motion in the range of 10(-7)-10(-3) sec. Our results indicate that the spin label is rigidly bound to S1 (tau 2 for isolated S1 is 2 X 10(-7) sec) and that the motion of the label reflects the motion of the S1 region of myosin. the value of tau 2 for the S1 segment of myosin is less than twice that for isolated S1, while the molecular weights differ by a factor of 4, indicating flexibility of myosin in agreement with the conclusions of Mendelson et al. (Biochemistry 12:2250-2255 [1973]). Adding F-actin increses tau 2 in either myosin or isolated S1 by a factor of mearly 103, indicating rigid immobilization of S1 by actin. Formation of myosin filaments (at an ionic strength of 0.15 or less) increses tau 2 by a factor of 10-30, depending on the ionic strength, indicating a decrease of the rotational mobility of S1 in these agregates. The remaining motion is at least a factor of 10 faster than would be expected for the filament itself, suggesting motion of the S1 region independent of the filament backbone but slower than in a single molecule. F-actin has a strong immobilizing effect on labeled S1 in myosin filaments (in 0.137 M KC1), but the immobilization is less complete than that observed when F-actin is added to labeled myosin monomers (in 0.5 M KC1). A spin-label analog of maleimide, attached to the SH-2 thiol groups of S1, is immobilized to a much lesser extent by F-actin than is the label on SH-1 groups. The maleimide label also was attached directly to F-actin and was sufficiently immobilized to suggest rigid binding to actin.     

How does internal motion influence the relaxation of the water protons in Ln(III)DOTA-like complexes?

Taking advantage of the Curie contribution to the relaxation of the protons in the Tb(III) complex, and the quadrupolar relaxation of the 17O and 2H nuclei on the Eu(III) complex, the effect of the internal motion of the water molecule bound to [Ln(DOTAM)(H2O)]3+ complexes was quantified. The determination of the quadrupolar coupling constant of the bound water oxygen chi(Omicron)(1 + eta(Omicron)2/3)1/2 = 5.2 +/- 0.5 MHz allows a new analysis of the 17O and 1H NMR data of the [Gd(DOTA)(H2O)]- complex with different rotational correlation times for the Gd(III)-O(water) and Gd(III)-H(water) vectors. The ratio of the rotational correlation times for the Ln(III)-H(water) vector and the overall rotational correlation time is calculated tau(RH)/tau(RO) = 0.65 +/- 0.2. This could have negative consequences on the water proton relaxivity, which we discuss in particular for macromolecular systems. It appears that the final effect is actually attenuated and should be around 10% for such large systems undergoing local motion of the chelating groups.     

Activation energy of catalysis-related domain motion in E. coli adenylate kinase

Adenylate kinase from E. coli (AKeco), folded into domains CORE, AMPbd, and LID, catalyzes the reaction AMP + ATP <--> 2ADP. Previous X-ray crystallography and optical solution methods showed that the domains AMPbd and LID, and the conserved P-loop, execute large-amplitude catalysis-related motions. We used (15)N NMR spin relaxation methods to find that the simplified model-free (MF) analysis does not, whereas our general Slowly Relaxing Local Structure analysis does, detect catalytic domain motion. SRLS set for the first time the correlation time for domain motion at tau(L)perpendicular = 8.2 ns, to be compared with tau(m) = 15.1 ns for global tumbling. These results were obtained at 303 K. Herein we conduct a temperature-dependent investigation of tau(L)perpendicular and tau(m) in the range of 288-310 K. We found that the activation energy for global tumbling is Ea = 16.9 +/- 0.5 kJ/mol, the hydrodynamic volume of hydrated AKeco is 65.6 +/- 2.1 nm3, its radius is 2.50 +/- 0.03 nm, and the number of hydration layers is 1.77. The average tau(L)perpendicular value decreases from 11 ns at 288 K to 4 ns at 310 K, with activation energies of 29.7 +/- 3.3, 32.1 +/- 4.3, and 30.4 +/- 4.3 kJ/mol for the domains AMPbd and LID, and the catalytic P-loop, respectively. These values are two-to-three times smaller than typical activation energies of enzymatic reactions. Hence kinase catalysis appears not to be controlled by domain motion in the ligand-free enzyme. However, the latter process clearly facilitates important mechanical aspects such as steric recognition and capturing of the AMP and ATP substrates, their proper positioning for phosphorylation, and the release of the ADP product.     

Rotational dynamics account for pH-dependent relaxivities of PAMAM dendrimeric, Gd-based potential MRI contrast agents

The EPTPA5) chelate, which ensures fast water exchange in GdIII complexes, has been coupled to three different generations (5, 7, and 9) of polyamidoamine (PAMAM) dendrimers through benzylthiourea linkages (H5EPTPA = ethylenepropylenetriamine-N,N,N',N',N'-pentaacetic acid). The proton relaxivities measured at pH 7.4 for the dendrimer complexes G5-(GdEPTPA)111, G7-(GdEPTPA)253 and G9-(GdEPTPA)1157 decrease with increasing temperature, indicating that, for the first time for dendrimers, slow water exchange does not limit relaxivity. At a given field and temperature, the relaxivity increases from G5 to G7, and then slightly decreases for G9 (r1 = 20.5, 28.3 and 27.9 mM(-1) s(-1), respectively, at 37 degrees C, 30 MHz). The relaxivities show a strong and reversible pH dependency for all three dendrimer complexes. This originates from the pH-dependent rotational dynamics of the dendrimer skeleton, which was evidenced by a combined variable-temperature and multiple-field 17O NMR and 1H relaxivity study performed at pH 6.0 and 9.9 on G5-(GdEPTPA)111. The longitudinal 17O and 1H relaxation rates of the dendrimeric complex are strongly pH-dependent, whereas they are not for the [Gd(EPTPA)(H2O)]2- monomer chelate. The longitudinal 17O and 1H relaxation rates have been analysed by the Lipari-Szabo spectral density functions and correlation times have been calculated for the global motion of the entire macromolecule (tau(gO)) and the local motion of the GdIII chelates on the surface (tau(lO)), correlated by means of an order parameter S2. The dendrimer complex G5-(GdEPTPA)111 has a considerably higher tau(gO) under acidic than under basic conditions (tau(298)gO = 4040 ps and 2950 ps, respectively), while local motions are less influenced by pH (tau(298)lO = 150 and 125 ps). The order parameter, characterizing the rigidity of the macromolecule, is also higher at pH 6.0 than at pH 9.9 (S2 = 0.43 vs 0.36, respectively). The pH dependence of the global correlation time can be related to the protonation of the tertiary amine groups in the PAMAM skeleton, which leads to an expanded and more rigid dendrimeric structure at lower pH. The increase of tau(gO) with decreasing pH is responsible for the pH dependent proton relaxivities. The water exchange rate on G5-(GdEPTPA)111(k(298)ex = 150 x 10(6) s(-1)) shows no significant pH dependency and is similar to the one measured for the monomer [Gd(EPTPA)(H2O)]2-. The proton relaxivity of G5-(GdEPTPA)111 is mainly limited by the important flexibility of the dendrimer structure, and to a small extent, by a faster than optimal water exchange rate.     

Solvent-dependent ultrafast internal conversion dynamics of n'-apo-beta-carotenoic-n'-acids (n = 8, 10, 12)

The ultrafast internal conversion dynamics of 12'-apo-beta-carotenoic-12'-acid (12'CA), 10'-apo-beta-carotenoic-10'-acid (10'CA) and 8'-apo-beta-carotenoic-8'-acid (8'CA) have been investigated by femtosecond pump-probe spectroscopy. The three apocarotenoic acids were excited to the S(2) state with different excess energies. Time constants tau(1) for the IC process S(1)/ICT --> S(0) were measured by probing the dynamics at 390 nm (S(0) --> S(2)), 575 nm (S(1)/ICT --> S(n)), 850, 860 and 890 nm (S(2) --> S(n) and S(1)/ICT --> S(0)). In nonpolar solvents, the observed reduction of the tau(1) values with increasing conjugation length of the acids is consistent with a reduction of the energy gap between the S(1)/ICT and S(0) states. The values are in good agreement with those of the corresponding apocarotenals studied previously in our groups. In polar solvents, a pronounced reduction of tau(1) values was observed for 12'CA, however the behavior was different from that observed for the respective aldehyde 12'-apo-beta-caroten-12'-al studied previously: First, the degree of tau(1) reduction in methanol was milder for 12'CA (218 --> 55 ps) than for 12'-apo-beta-caroten-12'-al (220 --> 8 ps). Secondly, for 12'CA the plateau of solvent independent tau(1) values extended further into the mid-polar range (up to 0.5 on the Deltaf scale) than previously observed for the 12'-aldehyde. For 10'CA the polarity effect on the tau(1) values was weaker ( approximately 71 ps in n-hexane and 34 ps in methanol) and for 8'CA it disappeared completely ( approximately 24 ps averaged over all solvents). The polarity-induced reduction of tau(1) is likely due to the stabilization of an intramolecular charge transfer state in polar solvents. This S(1)/ICT state is also responsible for the stimulated emission in the near IR, which has been observed in this specific class of carotenoids with a terminal carboxyl group for the first time. The occurrence of stimulated emission in the near IR region is also consistent with the steady-state fluorescence spectra which are reported along with the absorption spectra of these species. Possible reasons for the different behavior of the apocarotenoic acids compared to the respective aldehydes are discussed.     

Hydrophobic distal pocket affects NO-heme geminate recombination dynamics in dehaloperoxidase and H64V myoglobin

The recombination dynamics of NO with dehaloperoxidase (DHP) from Amphitrite ornata following photolysis were measured by femtosecond time-resolved absorption spectroscopy. Singular value decomposition (SVD) analysis reveals two important basis spectra. The first SVD basis spectrum reports on the population of photolyzed NO molecules and has the appearance of the equilibrium difference spectrum between the deoxy and NO forms of DHP. The first basis time course has two kinetic components with time constants of tau(11) approximately 9 ps and tau(12) approximately 50 ps that correspond to geminate recombination. The fast geminate process tau(11) arises from a contact pair with the heme iron in a bound state with S = 3/2 spin. The slow geminate process tau(12) corresponds to the recombination from a more remote docking site >3 A from the heme iron with the greater barrier corresponding to a S = 5/2 spin state. The second SVD basis spectrum represents a time-dependent Soret band shift indicative of heme photophysical processes and protein relaxation with time constants of tau(21) approximately 3 ps and tau(22) approximately 17 ps, respectively. A comparison between the more rapid rate constant of the slow geminate phase in DHP-NO and horse heart myoglobin (HHMbNO) or sperm whale myoglobin (SWMbNO) suggests that protein interactions with photolyzed NO are weaker in DHP than in the wild-type MbNOs, consistent with the hydrophobic distal pocket of DHP. The slower protein relaxation rate tau(22) in DHP-NO relative to HHMbNO implies less effective trapping in the docking site of the distal pocket and is consistent with a greater yield for the fast geminate process. The trends observed for DHP-NO also hold for the H64V mutant of SWMb (H64V MbNO), consistent with a more hydrophobic distal pocket for that protein as well. We examine the influence of solution viscosity on NO recombination by varying the glycerol content in the range from 0% to 90% (v/v). The dominant effect of increasing viscosity is the increase of the rate of the slow geminate process, tau(12), coupled with a population decrease of the slow geminate component. Both phenomena are similar to the effect of viscosity on wild-type Mb due to slowing of protein relaxation resulting from an increased solution viscosity and protein surface dehydration.     

Thermally activated escape rate for the Brownian motion of a fixed axis rotator in a double well potential for all values of the dissipation

The extension of the Kramers theory of the escape rate of a Brownian particle from a potential well to the entire range of damping proposed by Mel'nikov and Meshkov [J. Chem, Phys. 85, 1018 (1986)] is applied to the rotational Brownian motion of fixed axis rotators in a double well cosine potential. The procedure yields an expression for the Kramers escape rate valid for all values of the dissipation including the very low damping (VLD), very high damping (VHD), and crossover regimes. This equation provides a good asymptotic estimate of the correlation time tau per pendicular of the longitudinal dipole moment correlation function calculated by solving the underlying Langevin equation using the matrix-continued fraction method. Moreover, for low barriers, where the Mel'nikov and Meshkov approach is not applicable, analytic equations for tau in the VLD and VHD limits are derived and a simple extrapolating equation that is valid for all values of the damping is proposed.     

Site-Specific Hyperphosphorylation Inhibits,Rather than Promotes,Tau Fibrillization,Seeding Capacity,and Its Microtubule Binding

The consistent observation of phosphorylated tau in the pathology of Alzheimer's disease has contributed to the emergence of a model where hyperphosphorylation triggers both tau disassociation from microtubules and its subsequent aggregation. Herein, we applied a total chemical synthetic approach to site-specifically phosphorylate the microtubule binding repeat domain of tau (K18) at single (pS356) or multiple (pS356/pS262 and pS356/pS262/pS258) residues. We show that hyperphosphorylation of K18 inhibits 1) its aggregation in vitro, 2) its seeding activity in cells, 3) its binding to microtubules, and 4) its ability to promote microtubule polymerization. The inhibition increased with increasing the number of phosphorylated sites, with phosphorylation at S262 having the strongest effect. Our results argue against the hyperphosphorylation hypothesis and underscore the importance of revisiting the role of site-specific hyperphosphorylation in regulating tau functions in health and disease.     

Site‐Specific Hyperphosphorylation Inhibits,Rather than Promotes,Tau Fibrillization,Seeding Capacity,and Its Microtubule Binding

The consistent observation of phosphorylated tau in the pathology of Alzheimer's disease has contributed to the emergence of a model where hyperphosphorylation triggers both tau disassociation from microtubules and its subsequent aggregation. Herein, we applied a total chemical synthetic approach to site‐specifically phosphorylate the microtubule binding repeat domain of tau (K18) at single (pS356) or multiple (pS356/pS262 and pS356/pS262/pS258) residues. We show that hyperphosphorylation of K18 inhibits 1) its aggregation in vitro, 2) its seeding activity in cells, 3) its binding to microtubules, and 4) its ability to promote microtubule polymerization. The inhibition increased with increasing the number of phosphorylated sites, with phosphorylation at S262 having the strongest effect. Our results argue against the hyperphosphorylation hypothesis and underscore the importance of revisiting the role of site‐specific hyperphosphorylation in regulating tau functions in health and disease.     

Comparing ab initio density-functional and wave function theories: the impact of correlation on the electronic density and the role of the correlation potential

The framework of ab initio density-functional theory (DFT) has been introduced as a way to provide a seamless connection between the Kohn-Sham (KS) formulation of DFT and wave-function based ab initio approaches [R. J. Bartlett, I. Grabowski, S. Hirata, and S. Ivanov, J. Chem. Phys. 122, 034104 (2005)]. Recently, an analysis of the impact of dynamical correlation effects on the density of the neon atom was presented [K. Jankowski, K. Nowakowski, I. Grabowski, and J. Wasilewski, J. Chem. Phys. 130, 164102 (2009)], contrasting the behaviour for a variety of standard density functionals with that of ab initio approaches based on second-order M?ller-Plesset (MP2) and coupled cluster theories at the singles-doubles (CCSD) and singles-doubles perturbative triples [CCSD(T)] levels. In the present work, we consider ab initio density functionals based on second-order many-body perturbation theory and coupled cluster perturbation theory in a similar manner, for a range of small atomic and molecular systems. For comparison, we also consider results obtained from MP2, CCSD, and CCSD(T) calculations. In addition to this density based analysis, we determine the KS correlation potentials corresponding to these densities and compare them with those obtained for a range of ab initio density functionals via the optimized effective potential method. The correlation energies, densities, and potentials calculated using ab initio DFT display a similar systematic behaviour to those derived from electronic densities calculated using ab initio wave function theories. In contrast, typical explicit density functionals for the correlation energy, such as VWN5 and LYP, do not show behaviour consistent with this picture of dynamical correlation, although they may provide some degree of correction for already erroneous explicitly density-dependent exchange-only functionals. The results presented here using orbital dependent ab initio density functionals show that they provide a treatment of exchange and correlation contributions within the KS framework that is more consistent with traditional ab initio wave function based methods.     

Local hydrogen bonding dynamics and collective reorganization in water: ultrafast infrared spectroscopy of HOD/D(2)O

We present an investigation into hydrogen bonding dynamics and kinetics in water using femtosecond infrared spectroscopy of the OH stretching vibration of HOD in D(2)O. Infrared vibrational echo peak shift and polarization-selective pump-probe experiments were performed with mid-IR pulses short enough to capture all relevant dynamical processes. The experiments are self-consistently analyzed with a nonlinear response function expressed in terms of three dynamical parameters for the OH stretching vibration: the frequency correlation function, the lifetime, and the second Legendre polynomial dipole reorientation correlation function. It also accounts for vibrational-relaxation-induced excitation of intermolecular motion that appears as heating. The long time, picosecond behavior is consistent with previous work, but new dynamics are revealed on the sub-200 fs time scale. The frequency correlation function is characterized by a 50 fs decay and 180 fs beat associated with underdamped intermolecular vibrations of hydrogen bonding partners prior to 1.4 ps exponential relaxation. The reorientational correlation function observes a 50 fs librational decay prior to 3 ps diffusive reorientation. Both of these correlation functions compare favorably with the predictions from classical molecular dynamics simulations. The time-dependent behavior can be separated into short and long time scales by the 340 fs correlation time for OH frequency shifts. The fast time scales arise from dynamics that are mainly local: fluctuations in hydrogen bond distances and angles within relatively fixed intermolecular configurations. On time scales longer than the correlation time, dephasing and reorientations reflect collective reorganization of the liquid structure. Since the OH transition frequency and dipole are only weakly sensitive to these collective coordinates, this is a kinetic regime which gives an effective rate for exchange of intermolecular structures.     

Is hairpin formation in single-stranded polynucleotide diffusion-controlled?

An intriguing puzzle in biopolymer science is the observation that single-stranded DNA and RNA oligomers form hairpin structures on time scales of tens of microseconds, considerably slower than the estimated time for loop formation for a semiflexible polymer of similar length. To address the origin of the slow kinetics and to determine whether hairpin dynamics are diffusion-controlled, the effect of solvent viscosity (eta) on hairpin kinetics was investigated using laser temperature-jump techniques. The viscosity was varied by addition of glycerol, which significantly destabilizes hairpins. A previous study on the viscosity dependence of hairpin dynamics, in which all the changes in the measured rates were attributed to a change in solvent viscosity, reported an apparent scaling of relaxation times (tau(r)) on eta as tau(r) approximately eta(0.8). In this study, we demonstrate that if the effect of viscosity on the measured rates is not deconvoluted from the inevitable effect of change in stability, then separation of tau(r) into opening (tau(o)) and closing (tau(c)) times yields erroneous behavior, with different values (and opposite signs) of the apparent scaling exponents, tau(o) approximately eta(-0.4) and tau(c) approximately eta(1.5). Under isostability conditions, obtained by varying the temperature to compensate for the destabilizing effect of glycerol, both tau(o) and tau(c) scale as approximately eta(1.1+/-0.1). Thus, hairpin dynamics are strongly coupled to solvent viscosity, indicating that diffusion of the polynucleotide chain through the solvent is involved in the rate-determining step.     

NMR and monte carlo simulation studies of water adsorption dynamics

Nuclear magnetic resonance water proton spin-spin relaxation time T(2) was measured in wood samples with moisture contents ranging from 0.50 to 26.4%. The experimental results are discussed in terms of Monte Carlo simulations, which determine the correlation times tau for reorientation of the water molecule proton-proton vectors. We demonstrate that 1/T(2) and tau qualitatively follow the same behavior with surface hydration. The common application of the multisite exchange model to hydrated systems is discussed in light of the new results.     

Coupled cluster study of the energetic properties of S2x (x=0,+1,1)

Ab initio electronic structure calculations are reported for S2, and its ions S2+ and S2-. Geometric parameters are calculated using the singles and doubles coupled cluster method, including a perturbational correction for connected triple excitation, together with systematic sequences of correlation consistent basis sets extrapolated to the complete basis set (CBS) limit. Energetic and structural properties of S2 and the S2 cation and anion are reported. The heat of formation of S2 (3Sigmag-) at 0 K in the gas phase is predicted to be 29.8 kcal/mol from the average of CBS two extrapolation procedures, less than the experimental heat of formation of S2 of 30.66+/-0.07 kcal/mol. The 0 K adiabatic ionization potential and electron affinity are predicted to be 9.37 and 1.68 eV, respectively.     

Ab initio potential energy surface and intermolecular vibrations of the naphthalene-argon van der Waals complex

The intermolecular potential energy surface (PES) of the naphthalene-argon (NpAr) complex is constructed using an ab initio method. The molecule-argon interaction energy is computed at the level of the second-order M?ller-Plesset (MP2) theory combined with the augmented correlation consistent polarized valence double-ζ basis set. The analytical PES fitted to a large set of single energy values is further improved with the help of correction functions determined by calculations of the interaction energy at the coupled cluster level including single and double excitations supplemented by triple excitations performed for a limited set of intermolecular configurations. The PES determined is very flat near its four equivalent global minima of -493 cm(-1) located from both sides of the Np plane at a distance of 3.435 A? and shifted from the center of Np by ±0.43 A? along its long symmetry axis. The large-amplitude motion of Ar in the complex is investigated, and dynamical consequence of a strong intermode coupling is discovered in the excited vibrational states. The theoretical results obtained allow for the reassignment of the spectral bands observed in the electronic transition S(1) ← S(0) of the NpAr complex.