Temperature dependence of looping rates in a short peptide

Knowledge of the influence of chain length and amino acid sequence on the structural and dynamic properties of small peptides in solution provides essential information on protein folding pathways. The combination of time-resolved optical spectroscopy and molecular dynamics (MD) simulation methods has become a powerful tool to investigate the kinetics of end-to-end collisions (looping rates) in short peptides, which are relevant in early protein folding events. We applied the combination of both techniques to study temperature-dependent (280-340 K) looping rates of the Dbo-AlaGlyGln-Trp-NH2 peptide, where Dbo represents a 2,3-diazabicyclo[2.2.2]oct-2-ene-labeled asparagine, which served as a fluorescent probe in the time-resolved spectroscopic experiments. The experimental looping rates increased from 4.8 x 10(7) s(-1) at 283 K to 2.0 x 10(8) s(-1) at 338 K in H2O. The corresponding Arrhenius plot provided as activation parameters Ea = 21.5 +/- 1.0 kJ mol(-1) and ln(A/s-1) = 26.8 +/- 0.2 in H2O. The results in D2O were consistent with a slight solvent viscosity effect, i.e., the looping rates were 10-20% slower. MD simulations were performed with the GROMOS96 force field in a water solvent model, which required first a parametrization of the synthetic amino acid Dbo. After corrections for solvent viscosity effects, the calculated looping rates varied from 1.5 x 10(8) s(-1) at 280 K to 8.2 x 10(8) s(-1) at 340 K in H2O, which was about four times larger than the experimental data. The calculated activation parameters were Ea = 24.7 +/- 1.5 kJ mol(-1) and ln(A/s(-1)) = 29.4 +/- 0.1 in H2O.     

Kinetics of chain collapse in dilute polymer solutions: molecular weight and solvent dependences

The molecular weight and solvent dependences of the characteristic time of chain collapse were studied for poly(methyl methacrylate) (PMMA) of the molecular weight Mw=6.4x10(6) and 1.14x10(7) in pure acetonitrile (AcN) and in the mixed solvent of AcN+water (10 vol %). The size of PMMA chains was measured as a function of the time after the quench by static light scattering and the chain collapse processes were expressed by the plot of the expansion factor alpha2 vs ln t. The chain collapse in the mixed solvent AcN+water (10 vol %) was found to occur much faster than that in pure AcN, though the measurement of the former collapse process required several hours. In order to make a comparison between the rates of chain collapses, the fast chain collapse process was superposed on the slow one by scaling the time of the fast process as gammat. The scale factor gamma was determined by comparing the chain collapse processes of nearly the same equilibrium expansion factor with each other. Accordingly, the superposition of the collapse for Mw=6.4x10(6) on that for Mw=1.14x10(7) yielded gammam=4.0+/-0.6 for the process in AcN+water and 5.5+/-0.6 in AcN. The superposition of the chain collapse process in AcN+water on that in AcN yielded gammas=9.5+/-1.4 for Mw=6.4x10(6) and 12.0+/-1.8 for Mw=1.14x10(7). This analysis suggests that gammam and gammas are constant independent of each other. Thus, by assuming the molecular weight dependence of gammam approximately Mz, the characteristic time tauexp of chain collapse was conjectured as tauexp approximately kappaMz, where kappa reflects the nature of solvent species. The ratio of kappa for PMMA in AcN to that in AcN+water is given by gammas. The exponent was estimated to be z=2.4+/-0.7 for AcN+water and 3.0+/-0.7 for AcN. These values are compatible with the theoretical prediction z=3 based on a phenomenological model, though the observed characteristic times are longer by several orders of magnitude than those of the theoretical prediction.     

Numerical Simulation of the Effect of Solvent Viscosity on the Motions of a β‐Peptide Heptamer

This report examines the effect of a decrease in solvent viscosity on the simulated folding behaviour of a β‐peptide heptamer in methanol. Simulations of the molecular dynamics of the heptamer H‐β³‐HVal‐β³‐HAla‐β³‐HLeu‐(S,S)‐β³‐HAla(αMe)‐β³‐HVal‐β³‐HAla‐β³‐HLeu‐OH in methanol, with an explicit representation of the methanol molecules, were performed for 80 ns at various solvent viscosities. The simulations indicate that at a solvent viscosity of one third of that of methanol, only the dynamic aspects of the folding process are altered, and that the rate of folding is increased. At a viscosity of one tenth of that of methanol, insufficient statistics are obtained within the 80 ns period. We suggest that 80 ns is an insufficient time to reach conformational equilibrium at very low viscosity because the dependence of the folding rate of a β‐peptide on solvent viscosity has two regimes; a result that was observed in another computational study for α‐peptides.     

Viscosity scaling for the glassy phase of protein folding

Although commendable progress has been made in the understanding of the physics of protein folding, a key unresolved issue is whether Kramers' diffusion model of chemical reactions is generally applicable to activated barrier crossing events during folding. To examine the solvent viscosity effect on the folding transition of native-like trapped intermediates, laser flash photolysis has been used to measure the microsecond folding kinetics of a natively folded state of CO-liganded ferrocytochrome c (M-state) in the 1-250 cP range of glycerol viscosity at pH 7.0, 20 degrees C. The single rate coefficient for the folding of the M-state to the native state of the protein (i.e., the M --> N folding process) decreases initially when the solvent viscosity is low (<10 cP), but saturates at higher viscosity, indicating that Kramers model is not general enough for scaling the viscosity dependence of post-transition folding involving glassy dynamics. Analysis based on the Grote-Hynes idea of time dependent friction in conjunction with defect diffusion dynamics can account for the observed non-Kramers scaling.     

Effect of finite size on cooperativity and rates of protein folding

We analyze the dependence of cooperativity of the thermal denaturation transition and folding rates of globular proteins on the number of amino acid residues, N, using lattice models with side chains, off-lattice Go models, and the available experimental data. A dimensionless measure of cooperativity, Omega(c) (0 < Omega(c) < infinity), scales as Omega(c) approximately N(zeta). The results of simulations and the analysis of experimental data further confirm the earlier prediction that zeta is universal with zeta = 1 + gamma, where exponent gamma characterizes the susceptibility of a self-avoiding walk. This finding suggests that the structural characteristics in the denaturated state are manifested in the folding cooperativity at the transition temperature. The folding rates k(F) for the Go models and a dataset of 69 proteins can be fit using k(F) = k(F)0 exp(-cN(beta)). Both beta = 1/2 and 2/3 provide a good fit of the data. We find that k(F) = k(F)0 exp(-cN(1/2)), with the average (over the dataset of proteins) k(F)0 approximately (0.2 micros)(-1) and c approximately 1.1, can be used to estimate folding rates to within an order of magnitude in most cases. The minimal models give identical N dependence with c approximately 1. The prefactor for off-lattice Go models is nearly 4 orders of magnitude larger than the experimental value.     

Shear thinning and shear dilatancy of liquid n-hexadecane via equilibrium and nonequilibrium molecular dynamics simulations: Temperature, pressure, and density effects

Equilibrium and nonequilibrium molecular dynamics (MD) simulations have been performed in both isochoric-isothermal (NVT) and isobaric-isothermal (NPT) ensemble systems. Under steady state shearing conditions, thermodynamic states and rheological properties of liquid n-hexadecane molecules have been studied. Between equilibrium and nonequilibrium states, it is important to understand how shear rates (gamma) affect the thermodynamic state variables of temperature, pressure, and density. At lower shear rates of gamma<1 x 10(11) s(-1), the relationships between the thermodynamic variables at nonequilibrium states closely approximate those at equilibrium states, namely, the liquid is very near its Newtonian fluid regime. Conversely, at extreme shear rates of gamma>1 x 10(11) s(-1), specific behavior of shear dilatancy is observed in the variations of nonequilibrium thermodynamic states. Significantly, by analyzing the effects of changes in temperature, pressure, and density on shear flow system, we report a variety of rheological properties including the shear thinning relationship between viscosity and shear rate, zero-shear-rate viscosity, rotational relaxation time, and critical shear rate. In addition, the flow activation energy and the pressure-viscosity coefficient determined through Arrhenius and Barus equations acceptably agree with the related experimental and MD simulation results.     

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.     

Rheology of viscoelastic mixed surfactant solutions: effect of scission on nonlinear flow and rheochaos

The linear and nonlinear rheology of viscoelastic mixed anionic-zwitterionic surfactant solutions has been systematically investigated. In the linear viscoelastic regime, these systems display nearly Maxwellian behavior with a unique relaxation time, tau0, and a characteristic elastic plateau modulus, G0. Linear rheological data were used to calculate the repitation and breaking times of the micelles, tau(rep) and tau(b), respectively. Surprisingly, the elastic modulus G0 significantly increases with salt concentration c(s), whereas tau(b) decreases by 1 order of magnitude. The strong effect of c(s) on the material parameters and microstructure of rodlike micelles allowed for the systematic investigation of the effect of these parameters on nonlinear flow. For samples with relatively long tau(b), the quasi-static flow diagram (stress vs shear rate) shows a stress peak followed by a metastable branch (a region of decreasing shear stress), whereas for samples with relatively short tau(b), this phenomenon is not observed. Transient flow responses corroborate quasi-static flow findings and further reveal the significance of microscopic dynamic parameters on flow behavior. Shear stress time series were recorded at constant shear rates, and above a critical shear rate, gamma(c2), stress fluctuations are observed. The amplitude of these stress fluctuations, Delta sigma, was found to scale as Delta sigma approximately equal to G0(tau(b)| gamma - gamma(c2)|)beta with beta approximately 0.5. This scaling is observed for micellar systems with tau(b) ranging from 0.12 to 0.01 s and G0 ranging from 1 x 10(3) to 7 x 10(3) dyn/cm2.     

UV Raman spatially resolved melting dynamics of isotopically labeled polyalanyl peptide: slow alpha-helix melting follows 3(10)-helices and pi-bulges premelting

We used UV resonance Raman (UVRR) to examine the spatial dependence of the T-jump secondary structure relaxation of an isotopically labeled 21-residue mainly Ala peptide, AdP. The AdP penultimate Ala residues were perdeuterated, leaving the central residues hydrogenated, to allow separate monitoring of melting of the middle versus the end peptide bonds. For 5 to 30 degrees C T-jumps, the central peptide bonds show a approximately 2-fold slower relaxation time (189 +/- 31 ns) than do the exterior peptide bonds (97 +/- 15 ns). In contrast, for a 20 to 40 degrees C T-jump, the central peptide bond relaxation appears to be faster (56 +/- 6 ns) than that of the penultimate peptide bonds (131 +/- 46 ns). We show that, if the data are modeled as a two-state transition, we find that only exterior peptide bonds show anti-Arrhenius folding behavior; the middle peptide bonds show both normal Arrhenius-like folding and unfolding. This anti-Arrhenius behavior results from the involvement of pi-bulges/helices and 3(10)-helix states in the melting. The unusual temperature dependence of the (un)folding rates of the interior and exterior peptide bonds is due to the different relative (un)folding rates of 3(10)-helices, alpha-helices, and pi-bulges/helices. Pure alpha-helix unfolding rates are approximately 12-fold slower (approximately 1 micros) than that of pi-bulges and 3(10)-helices. In addition, we also find that the alpha-helix is most stable at the AdP N-terminus where eight consecutive Ala occur, whereas the three hydrophilic Arg located in the middle and at the C-terminus destabilize the alpha-helix in these regions and induce defects such as pi-bulges and 3(10)-helices.     

Bonding and (hyper)polarizability in the sodium dimer

We report a conventional ab initio and density functional theory study of the polarizability (alpha(alphabeta)/e(2)a(0) (2)E(h) (-1)) and hyperpolarizability (gamma(alphabetagammadelta)/e(4)a(0) (4)E(h) (-3)) of the sodium dimer. A large [18s14p9d2f1g] basis set is thought to yield near-Hartree-Fock values for both properties: alpha=272.28, Deltaalpha=127.22 and gamma=2157.6 x 10(3) at R(e)=3.078 87 A. Electron correlation has a remarkable effect on the Cartesian components of gamma(alphabetagammadelta). Our best value for the mean is gamma=1460.1 x 10(3). The (hyper)polarizability shows very strong bond-length dependence. The effect is drastically different for the longitudinal and transverse components of the hyperpolarizability. The following first derivatives were extracted from high-level coupled cluster calculations: (dalpha/dR)(e)=54.1, (dDeltaalpha/dR)(e)=88.1e(2)a(0)E(h) (-1), and (dgamma/dR)(e)=210 x 10(3)e(4)a(0) (3)E(h) (-3). We associate the (hyper)polarizability to bonding effects between the two sodium atoms by introducing the differential property per atom Q(diff)/2 identical with (Q[Na(2)(X (1)Sigma(g) (+))]/2-Q[Na((2)S)]). The differential (hyper)polarizability per atom is predicted to be strongly negative for the dimer at R(e), as [alpha(Na(2))/2-alpha(Na)]=-33.8 and [gamma(Na(2))/2-gamma(Na)]=-226.3 x 10(3). The properties calculated with the widely used B3LYP and B3PW91 density functional methods differ significantly. The B3PW91 results are in reasonable agreement with the conventional ab initio values. Last, we observe that low-level ab initio and density functional theory methods underestimate the dipole polarizability anisotropy. Experimental data on this important property are highly desirable.     

Numerical simulation of the effect of solvent viscosity on the motions of a beta-peptide heptamer

This report examines the effect of a decrease in solvent viscosity on the simulated folding behaviour of a beta-peptide heptamer in methanol. Simulations of the molecular dynamics of the heptamer H-beta3-HVal-beta3-HAla-beta3-HLeu-(S,S)-beta3-HAla(alphaMe)-beta3-HVal-beta3-HAla-beta3-HLeu-OH in methanol, with an explicit representation of the methanol molecules, were performed for 80 ns at various solvent viscosities. The simulations indicate that at a solvent viscosity of one third of that of methanol, only the dynamic aspects of the folding process are altered, and that the rate of folding is increased. At a viscosity of one tenth of that of methanol, insufficient statistics are obtained within the 80 ns period. We suggest that 80 ns is an insufficient time to reach conformational equilibrium at very low viscosity because the dependence of the folding rate of a beta-peptide on solvent viscosity has two regimes; a result that was observed in another computational study for alpha-peptides.     

Efficient rate enhancement due to intramolecular general base (IGB) assistance in the hydrolysis of N-(o-hydroxyphenyl)phthalimide

The rates of the hydrolyses of N-(o-hydroxyphenyl)phthalimide (1) and N-(o-methoxyphenyl)phthalimide (2), studied at different pH, show that the hydrolysis of 1 involves intramolecular general base (IGB) assistance where the o-O- group of ionized 1 acts as IGB and H2O as the reactant. The rate enhancement due to the IGB-assisted reaction of H2O with ionized 1 is>8x10(4)-fold. Pseudo-first-order rate constant for the reaction of water with 2 is approximately 2x10(3)-fold smaller than the first-order rate constant (0.10 s-1) for pH-independent hydrolysis of 1 within the pH range of 9.60-10.10. Second-order rate constants (kOH) for hydroxide ion-assisted hydrolysis of ionized 1 and 2 are 3.0 and 29.1 M-1 s-1, respectively. The solvent deuterium kinetic isotope effect (dKIE) on the rate of alkaline hydrolysis of 1 and 2 reveals that the respective values of kOH/kOD are 0.84 and 0.78, where kOD represents the second-order rate constant for DO--assisted cleavage of these imides (1 and 2). The value of kwH2O/kdD2O is 2.04, with kwH2O and kdD2O representing pseudo-first-order rate constants for the reactions of ionized 1 with H2O and D2O, respectively.     

Ozonolysis of oleic acid adsorbed to polar and nonpolar aerosol particles

Single-particle kinetic studies of the reaction between oleic acid and O 3 have been conducted on two different types of core particles: polystyrene latex (PSL) and silica. Oleic acid was found to adsorb to both particle types in multilayer islands that resulted in an adsorbed layer of a total volume estimated to be less than one monolayer. The rate of the surface reaction between surface-adsorbed oleic acid and O 3 has been shown for the first time to be influenced by the composition of the aerosol substrate in a mixed organic/inorganic particle. A Langmuir-Hinshelwood mechanism was applied to the observed dependence of the pseudo-first-order rate constant with [O 3], and the resulting fit parameters for the ozone partition coefficient ( K O 3 ) and maximum first order rate constant ( k 1,max ) suggest that the reaction proceeded faster on the less polar PSL core at lower [O 3] due to the increased residence time of O 3 on the PSL surface, but the reaction was ultimately more efficient on the silica surface at high [O 3]. Values for the uptake coefficient, gamma oleic , for reaction of oleic acid on PSL spheres decrease from 2.5 x 10 (-5) to 1 x 10 (-5) with increasing [O 3] from 4 to 25 ppm and overlap at high [O 3] with the estimated values for gamma oleic on silica, which decrease from 1.6 x 10 (-5) to 1.3 x 10 (-5). The relationship between gamma oleic and the more common expression for gamma O 3 is discussed.     

Solvent viscosity dependence of the protein folding dynamics

Solvent viscosity has been frequently adopted as an adjustable parameter in various computational studies (e.g., protein folding simulations) with implicit solvent models. A common approach is to use low viscosities to expedite simulations. While using viscosities lower than that of aqueous is unphysical, such treatment is based on observations that the viscosity affects the kinetics (rates) in a well-defined manner as described by Kramers' theory. Here, we investigate the effect of viscosity on the detailed dynamics (mechanism) of protein folding. On the basis of a simple mathematical model, we first show that viscosity may indeed affect the dynamics in a complex way. By applying the model to the folding of a small protein, we demonstrate that the detailed dynamics is affected rather pronouncedly especially at unphysically low viscosities, cautioning against using such viscosities. In this regard, our model may also serve as a diagnostic tool for validating low-viscosity simulations. It is also suggested that the viscosity dependence can be further exploited to gain information about the protein folding mechanism.     

The adsorption of cetyltrimethylammonium bromide and propanol mixtures with regard to wettability of polytetrafluoroethylene II. Adsorption at polytetrafluoroethylene-aqueous solution interface and wettability

Measurements of contact angles (theta) of aqueous solutions of cetyltrimethylammonium bromide (CTAB) and propanol mixtures at constant CTAB concentration equal to 1x10(-5), 1x10(-4), 6x10(-4) and 1x10(-3) M on polytetrafluoroethylene (PTFE) were carried out. The obtained results indicate that the wettability of PTFE by aqueous solutions of these mixtures depends on their composition and concentration. They also indicate that, contrary to Zisman, there is no linear relationship between cos theta and the surface tension (gamma(LV)), but a linear relationship exists between the adhesional (gamma(LV)cos theta) and surface tension of aqueous solutions of CTAB and propanol mixtures. Curve gamma(LV)cos theta vs gamma(LV) has a slope equal -1 suggesting that adsorption of CTAB and propanol mixtures and the orientation of their molecules at aqueous solution-air and PTFE-aqueous solution interfaces is the same. Extrapolating this curve to the value of gamma(LV)cos theta corresponding to theta=0, the value of the critical tension of PTFE wetting equal 23.4 mN/m was determined. This value was higher than that obtained from contact angles of n-alkanes on PTFE surface (20.24 mN/m). The difference between the critical surface tension values of wetting probably resulted from the fact that at cos theta=1 the PTFE-aqueous solution of CTAB and propanol mixture interface tension was not equal to zero. This tension was determined on the basis of the measured contact angles and Young equation. It appeared that the values of PTFE-aqueous solution of the CTAB and propanol mixtures interface tension can be satisfactorily determined by modified Szyszkowski equation only for solutions in which probably CTAB and propanol molecules are present in monomeric form. However, it appeared that using the equation of Miller et al., in which the possibility of aggregation of propanol molecules in the interface layer is taken into account, it is possible to describe the PTFE-solution interfacial tension for all systems studied in the same way as by the Young equation. On the basis of linear dependence between the adhesional and surface tension it was established that the work of adhesion of aqueous solution of CTAB and propanol mixtures does not depend on its composition and concentration, and the average value of this work was equal to 46.85 mJ/m(2), which was similar to that obtained for adhesion of aqueous solutions of two cationic surfactants mixtures to PTFE surface.     

The nucleation rate of crystalline ice in amorphous solid water

The kinetics of crystalline ice nucleation and growth in nonporous, molecular beam deposited amorphous solid water (ASW) films are investigated at temperatures near 140 K. We implement an experimental methodology and corresponding model of crystallization kinetics to decouple growth from nucleation and quantify the temperature dependence and absolute rates of both processes. Nucleation rates are found to increase from approximately 3x10(13) m(-3) s(-1) at 134 K to approximately 2x10(17) m(-3) s(-1) at 142 K, corresponding to an Arrhenius activation energy of 168 kJ/mol. Over the same temperature range, the growth velocity increases from approximately 0.4 to approximately 4 A s(-1), also exhibiting Arrhenius behavior with an activation energy of 47 kJ/mol. These nucleation rates are up to ten orders of magnitude larger than in liquid water near 235 K, while growth velocities are approximately 10(9) times smaller. Crystalline ice nucleation kinetics determined in this study differ significantly from those reported previously for porous, background vapor deposited ASW, suggesting the nucleation mechanism is dependent upon film morphology.     

Slowdown of H/D exchange reaction rate and water dynamics in ionic liquids: deactivation of solitary water solvated by small anions in 1-butyl-3-methyl-imidazolium chloride

The H/D exchange reaction and the rotational dynamics of heavy water (D2O) are studied at 50 degrees C in the ionic liquid, 1-butyl-3-methylimidazolium chloride ([bmim][Cl]), in the [D2O] range of 3-55 M. The initial H/D exchange rates are observed as 1.0 x 10(-7), 4.5 x 10(-6), 1.0 x 10(-5), 4.1 x 10(-5), 1.1 x 10(-4), and 3.7 x 10(-4) s(-1), respectively, at [D2O] of 2.8, 7.1, 8.1, 11, 15, and 25 M. The rate is very slow and less than 10(-5) s(-1) at [D2O] below approximately 7 M. It steeply increases to the order of 10(-4)s(-1) for 7 M < [D2O] < 10 M, and linearly increases with [D2O] in the more water-rich region. The intercept of the linear region at [D2O] = approximately 9 M is interpreted by considering that each chloride anion deactivates 1.6 equiv water molecules due to the strong solvation. Correspondingly, the rotational correlation time of D2O at [D2O] < 7 M is 1 order of magnitude larger than that in water-rich conditions.     

Gibbs monolayers at the air/water interface: surface partitioning and lateral mobility of an electrochemically active surfactant

Monolayer films of a water-soluble surfactant, 4-octaneamido-2,2,6,6-tetramethyl-1-piperidinyloxy (C8-TEMPO) were investigated at the air/water interface. An electrochemical, horizontal touch method was developed to measure the equilibrium surface concentrations (gamma) of C8TEMPO. The dependence of gamma on the solution concentration followed a Langmuir isotherm and yielded the partition constant K = (2.3+/-0.2) x 10(4) M(-1). These results were verified by surface tension measurements and Brewster angle microscopy. Within experimental error, the same K values were obtained. The lateral diffusion constants vs surface concentration of this molecule were measured by 2D voltammetry. In these experiments, the component of the oxidation current due to C8TEMPO in the bulk of the solution was subtracted from the total measured current to obtain the component due to the lateral surface diffusion. In the ange of mean molecular areas from 120 to 400 A2/molecule, the lateral diffusion constant of C8TEMPO increased from 1.0 x 10(-6) to 1.0 x 10(-5) cm2/s. The latter value is about 2.5 times larger than the C8TEMPO diffusion constant in bulk water. Comparison of the lateral mobilities of C8TEMPO and two longer alkane chain, water-insoluble homologues, C14TEMPO and C18TEMPO, showed no statistically significant differences.     

Measurement of the dependence of interfacial charge-transfer rate constants on the reorganization energy of redox species at n-ZnO/H2O interfaces

The interfacial energetic and kinetics behavior of n-ZnO/H2O contacts have been determined for a series of compounds, cobalt trisbipyridine (Co(bpy)3(3+/2+)), ruthenium pentaamine pyridine (Ru(NH3)5 py(3+/2+)), cobalt bis-1,4,7-trithiacyclononane (Co(TTCN)2(3+/2+)), and osmium bis-dimethyl bipyridine bis-imidazole (Os(Me2bpy)2(Im)2(3+/2+)), which have similar formal reduction potentials yet which have reorganization energies that span approximately 1 eV. Differential capacitance vs potential and current density vs potential measurements were used to measure the interfacial electron-transfer rate constants for this series of one-electron outer-sphere redox couples. Each interface displayed a first-order dependence on the concentration of redox acceptor species and a first-order dependence on the concentration of electrons in the conduction band at the semiconductor surface, in accord with expectations for the ideal model of a semiconductor/liquid contact. Rate constants varied from 1 x 10(-19) to 6 x 10(-17) cm4 s(-1). The interfacial electron-transfer rate constant decreased as the reorganization energy, lambda, of the acceptor species increased, and a plot of the logarithm of the electron-transfer rate constant vs (lambda + deltaG(o)')(2)/4lambda k(B)T (where deltaG(o)' is the driving force for interfacial charge transfer) was linear with a slope of approximately -1. The rate constant at optimal exoergicity was found to be approximately 5 x 10(-17) cm4 s(-1) for this system. These results show that interfacial electron-transfer rate constants at semiconductor electrodes are in good agreement with the predictions of a Marcus-type model of interfacial electron-transfer reactions.