Two-state Analysis Research Articles

Conformational Analysis of a Tetrasaccharide Based on NMR Spectroscopy and Molecular Dynamics Simulations

The conformational preference of the human milk oligosaccharide lacto-N-neotetraose, beta-d-Galp-(1 --> 4)-beta-d-GlcpNAc-(1 --> 3)-beta-d-Galp-(1 --> 4)-d-Glcp, has been analyzed using (1)H,(1)H T-ROESY and (1)H,(13)C trans-glycosidic J coupling experiments in isotropic solution and (1)H,(13)C residual dipolar couplings (RDCs) obtained in lyotropic liquid crystalline media. Molecular dynamics simulations of the tetrasaccharide with explicit water as the solvent revealed that two conformational states are significantly populated at the psi glycosidic torsion angle, defined by C(anomeric)-O-C-H, of the (1 --> 3)-linkage. Calculation of order parameters, related to the molecular shape, were based on the inertia tensor and fitting of experimental RDCs to different conformational states showed that psi(+) > 0 degrees is the major and psi(-) < 0 degrees is the minor conformation in solution, in complete agreement with a two-state analysis based on the T-ROESY data. Attention was also given to the effect of salt (200 mM NaCl) in the anisotropic medium, which was a ternary mixture of n-octyl-penta(ethylene glycol), n-octanol, and D(2)O.

Two-state analysis of single-molecule Raman spectra of crystal violet

Surface-enhanced Raman spectra of individual crystal violet molecules adsorbed on silver colloids were measured and analyzed. It was found that, because of the small number of hot spots on the surface of the colloids, the single-molecule limit is attained already at a molecule to colloid ratio of 10 3. Series of time-dependent spectra measured from individual molecules showed significant spectral jumps, involving essentially only the low-frequency part of the spectrum. Each spectrum could therefore be analyzed as belonging to one of two groups, depending on the ratio between its two parts. Spectral trajectories were converted into simplified two-state (high/low) trajectories. Probability distributions of the dwell times in the two states were constructed from the trajectories, and statistical analysis showed that the dynamic process underlying the spectral jumps is likely to be Markovian. We suggest that the jumps originate in fluctuations of molecule-surface charge transfer interactions, although our studies did not provide evidence for lateral motion of the molecules, as found previously with rhodamine 6G.

An Obligatory Intermediate in the Folding Pathway of Cytochromec552 from Hydrogenobacterthermophilus

The folding mechanism of many proteins involves the population of partially organized structures en route to the native state. Identification and characterization of these intermediates is particularly difficult, as they are often only transiently populated and may play different mechanistic roles, being either on-pathway productive species or off-pathway kinetic traps. Following different spectroscopic probes, and employing state-of-the-art kinetic analysis, we present evidence that the folding mechanism of the thermostable cytochrome c552 from Hydrogenobacter thermophilus does involve the presence of an elusive, yet compact, on-pathway intermediate. Characterization of the folding mechanism of this cytochrome c is particularly interesting for the purpose of comparative folding studies, because H. thermophilus cytochrome c552 shares high sequence identity and structural homology with its homologue from the mesophilic bacterium Pseudomonas aeruginosa cytochrome c551, which refolds through a broad energy barrier without the accumulation of intermediates. Analysis of the folding kinetics and correlation with the three-dimensional structure add new evidence for the validity of a consensus folding mechanism in the cytochrome c family.

Driving the Localized-to-Delocalized Transition in Unsymmetrical Dinuclear Ruthenium Mixed-Valence Complexes

The degree of delocalization in the symmetrical complexes [{Ru(bpy)2}2(μ-dpb′)]5+ and [{Ru(bpy)2}2(μ-dpb)]5+ (dpb = 2,3-bis(2-pyridyl)-1,4-benzoquinoxaline; dpb′ = dipyrido(2,3-a;3′,2′-c)benzophenazine; bpy = 2,2′-bipyridine) is diminished by the substitution of the terminal bpy ligands at one end of the complex. The results of a classical analysis for the diastereoisomeric forms of the series of complexes [{Ru(bpy)2}(μ-BL){Ru(pp)2}]5+ (pp = bpy, Me2bpy (4,4′-dimethyl-2,2′-bipyridine), Me4bpy (4,4′,5,5′-tetramethyl-2,2′-bipyridine)) indicate that a greater degree of ground-state delocalization exists in the complexes incorporating the bridging ligand dpb′ compared with the dpb analogue. A two-state analysis in which ΔE 0 (the redox asymmetry) is varied at constant H ab (the electronic coupling) and λ (the reorganizational energy) does not adequately describe the properties of the systems due to the importance of a third electronic state corresponding to the bridging ligand.

Frequency Analysis of Yield for Delineating Yield Response Zones

The yield in any given field or management zone is a product of interaction between many soil properties and production inputs. Therefore, multi-year yield maps may give better insight into determining potential management zones. This research was conducted to develop a methodology to delineate yield response zones by using two-state frequency analysis conducted on yield maps for 3 years on two commercial corn fields near Wiggins, Colorado. A zone was identified by the number of years that yield was equal and greater than the average yield in a given year. Classes producing statistically similar yield were combined resulting in three potential yield zones. Results indicated that the variability of yield over time and space could successfully be assessed at the same time without the drawbacks of averaging data from different years. Frequency analysis of multi-year yield data could be an effective way to establish yield response zones. Seventeen percent of the field #1 consistently produced lower yield than the mean while 43 of the field produced yield over the mean. Corresponding values for field #2 were 6% and 42%.The remainder of the fields produced fluctuating yields between years. These spatially and temporally sound yield response maps could be used to identify the yield-limiting factors in zones where yield is either low or fluctuating. Yield response maps could also be helpful to delineate potential management zones with the help of resource zones such as electrical conductivity and soil maps, along with the directed soil sampling results.

How delocalized is N,N,N',N'-tetraphenylphenylenediamine radical cation? An experimental and theoretical study on the electronic and molecular structure.

The electronic and molecular structure of N,N,N',N'-tetraphenylphenylenediamine radical cation 1(+) is in focus of this study. Resonance Raman experiments showed that at least eight vibrational modes are strongly coupled to the optical charge resonance band which is seen in the NIR. With the help of a DFT-based vibrational analysis, these eight modes were assigned to symmetric vibrations. The contribution of these symmetric modes to the total vibrational reorganization energy is dominant. These findings are in agreement with the conclusions from a simple two-state two-mode Marcus-Hush analysis which yields a tiny electron-transfer barrier. The excellent agreement of the X-ray crystal structure analysis and the DFT computed molecular structure of 1(+) on one hand as well as the solvent and solid-state IR spectra and the DFT-calculated IR active vibrations on the other hand prove 1(+) adopts a symmetrical delocalized Robin-Day class III structure both in the solid state and in solution.

Effect of drastic sequence alteration and D-amino acid incorporation on the membrane binding behavior of lytic peptides.

The amphipathic alpha-helix is a common motif found in many cell lytic peptides including antimicrobial peptides. We have recently shown that significantly altering the amphipathic structure of a lytic peptide by reshuffling its sequence and/or replacing a few l-amino acids with their D-enantiomers did not significantly affect the antimicrobial activity of the peptides nor their ability to bind and permeate negatively charged (PE/PG) membranes. However, a pronounced effect was observed regarding their hemolytic activity and their ability to bind and permeate zwitterionic (PC/Cho) membranes. To shed light on these findings, here we used surface plasmon resonance (SPR) with mono- and bilayer membranes. We found that the L-amino acid (aa) peptides bound 10-25-fold stronger to PC/Cho bilayers compared with monolayers, whereas the diastereomers bound similarly to both membranes. A two-state reaction model analysis of the data indicated that this difference is due to the insertion of the L-aa peptides into the PC/Cho bilayers, whereas the diastereomers are surface-localized. In contrast, only an approximately 2-fold difference was found with negatively charged membranes. Changes in the amphipathicity markedly affected only the insertion of the L-aa peptides into PC/Cho bilayers. Furthermore, whereas the all-L-aa peptides bound similarly to the PC/Cho and PE/PG membranes, the diastereomers bound approximately 100-fold better to PE/PG compared with PC/Cho membranes, and selectivity was determined only in the first binding step. The effect of the peptides on the lipid order determined by using ATR-FTIR studies supported these findings. Besides shedding light on the mode of action of these peptides, the present study demonstrates SPR as a powerful tool to differentiate between non-cell-selective compared with bacteria-selective peptides, based on differences in their membrane binding behavior.

The side chain of aspartic acid 69 dictates the folding mechanism of Bacillus subtilis HPr.

Many small, single-domain proteins show equilibrium and kinetic folding mechanisms that appear to be adequately described as two state. The two-state model makes several predictions that can be tested experimentally. First, the conformational stability determined at or extrapolated to a set of reference conditions should be independent of the measurement method (thermal or solvent denaturation or hydrogen exchange). Second, model-independent measures of the cardinal thermodynamic parameters (T(m), DeltaH) as determined from direct calorimetric means should be identical to those determined from the two-state analysis of thermal unfolding data. Third, the ratio of the kinetic folding and unfolding rate constants should be equal to K(eq) determined from an equilibrium measurement under the same conditions. Here, we show that the wild-type HPr protein from Bacillus subtilis does not meet all of these criteria under our standard conditions. However, if we replace the side chain of Asp69, or add moderate concentrations of salt, we find excellent two-state behavior in both equilibrium and kinetic folding. Thus, for this protein and possibly others, very subtle changes in the primary structure or in the solution conditions can dramatically alter the relative stabilities of the native intermediate, and unfolded ensembles can cause an observable change in the nature of the folding mechanism.

Folding of horse cytochrome c in the reduced state

Equilibrium and kinetic folding studies of horse cytochrome c in the reduced state have been carried out under strictly anaerobic conditions at neutral pH, 10°C, in the entire range of aqueous solubility of guanidinium hydrochloride (GdnHCl). Equilibrium unfolding transitions observed by Soret heme absorbance, excitation energy transfer from the lone tryptophan residue to the ferrous heme, and far-UV circular dichroism (CD) are all biphasic and superimposable, implying no accumulation of structural intermediates. The thermodynamic parameters obtained by two-state analysis of these transitions yielded ΔG(H2O) = 18.8(±1.45) kcal mol−1, and Cm = 5.1(±0.15) M GdnHCl, indicating unusual stability of reduced cytochrome c. These results have been used in conjunction with the redox potential of native cytochrome c and the known stability of oxidized cytochrome c to estimate a value of −164 mV as the redox potential of the unfolded protein. Stopped-flow kinetics of folding and unfolding have been recorded by Soret heme absorbance, and tryptophan fluorescence as observables. The refolding kinetics are monophasic in the transition region, but become biphasic as moderate to strongly native-like conditions are approached. There also is a burst folding reaction unobservable in the stopped-flow time window. Analyses of the two observable rates and their amplitudes indicate that the faster of the two rates corresponds to apparent two-state folding (U ↔ N) of 80–90% of unfolded molecules with a time constant in the range 190–550 μs estimated by linear extrapolation and model calculations. The remaining 10–20% of the population folds to an off-pathway intermediate, I, which is required to unfold first to the initial unfolded state, U, in order to refold correctly to the native state, N (I ↔ U ↔ N). The slower of the two observable rates, which has a positive slope in the linear functional dependence on the denaturant concentration indicating that an unfolding process under native-like conditions indeed exists, originates from the unfolding of I to U, which rate-limits the overall folding of these 10–20% of molecules. Both fast and slow rates are independent of protein concentration and pH of the refolding milieu, suggesting that the off-pathway intermediate is not a protein aggregate or trapped by heme misligation. The nature or type of unfolded-state heme ligation does not interfere with refolding. Equilibrium pH titration of the unfolded state yielded coupled ionization of the two non-native histidine ligands, H26 and H33, with a pKa value of 5.85. A substantial fraction of the unfolded population persists as the six-coordinate form even at low pH, suggesting ligation of the two methionine residues, M65 and M80. These results have been used along with the known ligand-binding properties of unfolded cytochrome c to propose a model for heme ligation dynamics. In contrast to refolding kinetics, the unfolding kinetics of reduced cytochrome c recorded by observation of Soret absorbance and tryptophan fluorescence are all slow, simple, and single-exponential. In the presence of 6.8 M GdnHCl, the unfolding time constant is ∼300(±125) ms. There is no burst unfolding reaction. Simulations of the observed folding-unfolding kinetics by numerical solutions of the rate equations corresponding to the three-state I ↔ U ↔ N scheme have yielded the microscopic rate constants.

Stabilization of apoflavodoxin by replacing hydrogen-bonded charged Asp or Glu residues by the neutral isosteric Asn or Gln.

Knowledge of protein stability principles provides a means to increase protein stability in a rational way. Here we explore the feasibility of stabilizing proteins by replacing solvent-exposed hydrogen-bonded charged Asp or Glu residues by the neutral isosteric Asn or GLN: The rationale behind this is a previous observation that, in some cases, neutral hydrogen bonds may be more stable that charged ones. We identified, in the apoflavodoxin from Anabaena PCC 7119, three surface-exposed aspartate or glutamate residues involved in hydrogen bonding with a single partner and we mutated them to asparagine or glutamine, respectively. The effect of the mutations on apoflavodoxin stability was measured by both urea and temperature denaturation. We observed that the three mutant proteins are more stable than wild-type (on average 0.43 kcal/mol from urea denaturation and 2.8 degrees C from a two-state analysis of fluorescence thermal unfolding data). At high ionic strength, where potential electrostatic repulsions in the acidic apoflavodoxin should be masked, the three mutants are similarly more stable (on average 0.46 kcal/mol). To rule out further that the stabilization observed is due to removal of electrostatic repulsions in apoflavodoxin upon mutation, we analysed three control mutants and showed that, when the charged residue mutated to a neutral one is not hydrogen bonded, there is no general stabilizing effect. Replacing hydrogen-bonded charged Asp or Glu residues by Asn or Gln, respectively, could be a straightforward strategy to increase protein stability.

Enthalpy distribution for the alpha-helix/random coil transition in a model peptide: a study of two-state behavior.

We use heat capacity data of Taylor et al. to calculate the enthalpy distribution of a model peptide using the moments/maximum-entropy method. The peptide was designed with small covalent loops at both ends of the molecule to nucleate alpha-helix thus giving a system that would be expected to show a helix-coil transition that is very close to being two state. If we subtract a background contribution from the heat capacity data, then the enthalpy distribution we obtain shows two distinct peaks representing helix and coil. The difference in the peak enthalpy values agrees closely with the DeltaH obtained from the two-state analysis. On the other hand, if we use the complete heat capacity without subtracting background we then obtain an enthalpy distribution that has only a single peak at all temperatures. We show that this result can be consistent with the existence of two states, helix and coil, but only if the range of variation of the enthalpy of each species is so large as to make the notion of a species fairly meaningless.

A change in the apparent m value reveals a populated intermediate under equilibrium conditions in Escherichia coli ribonuclease HI.

Experimental studies of protein stability often rely on the determination of an "m value", which describes the denaturant dependence of the free energy change between two states (DeltaG = DeltaG(H2O) - m[denaturant]). Changes in the m value accompanying site specific mutations are usually attributed to structural alterations in the native or unfolded ensemble. Here, we provide an example of significant reduction in the m value resulting from a subtle deviation in two-state behavior not detected by traditional methods. The protein that is studied is a variant of Escherchia coli RNase H in which three residues predicted to be involved in a partially buried salt bridge network were mutated to alanine (R46A, D102A, and D148A). Equilibrium denaturant profiles monitored by both fluorescence and circular dichroism appeared to be cooperative, and a two-state analysis yielded a DeltaG(UN) of approximately -3 kcal/mol with an m value of 1.4 kcal mol(-1) M(-1) (vs 2.3 for RNase H). Analysis of kinetic refolding experiments suggests that the system is actually three-state at equilibrium with an appreciable concentration of an intermediate state under low denaturant concentrations. The stability of the native state determined from a fit of these kinetic data is -6.7 kcal/mol, suggesting that the stability determined by traditional two-state equilibrium analysis is a gross underestimate. The only hint to this loss of two-state behavior was a decrease in the apparent m value, and the presence of the equilibrium intermediate was only identified by a kinetic analysis. Our work serves as a cautionary note; the possibility of a three-state system should be closely addressed before interpreting a change in the m value as a change in the native or unfolded state.

Melting studies of short DNA hairpins containing the universal base 5-nitroindole

Effects of the universal base 5-nitroindole on the thermodynamic stability of DNA hairpins having a 6 bp stem and four base loops were investigated by optical absorbance and differential scanning calorimetry techniques. Melting studies were conducted in buffer containing 115 mM Na(+). Five different modified versions of DNA hairpins containing a 5-nitroindole base or bases substituted at different positions in the stem and loop regions were examined. Thermo-dynamic parameters of the melting transitions estimated from a two-state analysis of optical melting curves and measured directly by calorimetry revealed that the presence of 5-nitroindole bases in the duplex stem or loop regions of short DNA hairpins significantly affects both their enthalpic and entropic melting components in a compensating manner, while the transition free energy varies linearly with the transition temperature. The calorimetrically determined enthalpy and entropy values of the modified hairpins were considerably smaller (43-53%) than the two-state optical parameters, suggesting that solvent effects may be significant in the melting processes of these hairpins. Results of circular dichroism measurements also revealed slight differences between the modified hairpins and the control in both the duplex and melted states, suggesting subtle structural differences between the control and DNA hairpins containing a 5-nitroindole base or bases.

Melting studies of short DNA hairpins: influence of loop sequence and adjoining base pair identity on hairpin thermodynamic stability.

Spectroscopic and calorimetric melting studies of 28 DNA hairpins were performed. These hairpins form by intramolecular folding of 16 base self-complementary DNA oligomer sequences. Sequence design dictated that the hairpin structures have a six base pair duplex linked by a four base loop and that the first five base pairs in the stem are the same in every molecule. Only loop sequence and identity of the duplex base pair closing the loop vary for the set of hairpins. For these DNA samples, melting studies were carried out to investigate effects of the variables on hairpin stability. Stability of the 28 oligomers was ascertained from their temperature-induced melting transitions in buffered 115 mM Na(+) solvent, monitored by ultraviolet absorbance and differential scanning calorimetry (DSC). Experiments revealed the melting temperatures of these molecules range from 32.4 to 60.5 degrees C and are concentration independent over strand concentrations of 0.5 to 260 microM; thus, as expected for hairpins, the melting transitions are apparently unimolecular. Model independent thermodynamic transition parameters, DeltaH(cal), DeltaS(cal), and DeltaG(cal), were determined from DSC measurements. Model dependent transition parameters, DeltaH(vH), DeltaS(vH), and DeltaG(vH) were estimated from a van't Hoff (two-state) analysis of optical melting transitions. Results of these studies reveal a significant sequence dependence to DNA hairpin stability. Thermodynamic parameters evaluated by either procedure reveal the transition enthalpy, DeltaH(cal) (DeltaH(vH)) can differ by as much as 20 kcal/mol depending on sequence. Similarly, values of the transition entropy DeltaS(cal) (DeltaS(vH)) can differ by as much as 60 cal/Kmol (eu) for different molecules. Differences in free energies DeltaG(cal) (DeltaG(vH)) are as large as 4 kcal/mol for hairpins with different sequences. Comparisons between the model independent calorimetric values and the thermodynamic parameters evaluated assuming a two-state model reveal that 10 of the 28 hairpins display non-two-state melting behavior. The database of sequence-dependent melting free energies obtained for the hairpins was employed to extract a set of n-n (nearest-neighbor) sequence dependent loop parameters that were able to reproduce the input data within error (with only two exceptions). Surprisingly, this suggests that the thermodynamic stability of the DNA hairpins can in large part be reasonably represented in terms of sums of appropriate nearest-neighbor loop sequence parameters.

Energetic and spatial parameters for gating of the bacterial large conductance mechanosensitive channel, MscL.

MscL is multimeric protein that forms a large conductance mechanosensitive channel in the inner membrane of Escherichia coli. Since MscL is gated by tension transmitted through the lipid bilayer, we have been able to measure its gating parameters as a function of absolute tension. Using purified MscL reconstituted in liposomes, we recorded single channel currents and varied the pressure gradient (P) to vary the tension (T). The tension was calculated from P and the radius of curvature was obtained using video microscopy of the patch. The probability of being open (Po) has a steep sigmoidal dependence on T, with a midpoint (T1/2) of 11.8 dyn/cm. The maximal slope sensitivity of Po/Pc was 0.63 dyn/cm per e-fold. Assuming a Boltzmann distribution, the energy difference between the closed and fully open states in the unstressed membrane was DeltaE = 18.6 kBT. If the mechanosensitivity arises from tension acting on a change of in-plane area (DeltaA), the free energy, TDeltaA, would correspond to DeltaA = 6.5 nm2. MscL is not a binary channel, but has four conducting states and a closed state. Most transition rates are independent of tension, but the rate-limiting step to opening is the transition between the closed state and the lowest conductance substate. This transition thus involves the greatest DeltaA. When summed over all transitions, the in-plane area change from closed to fully open was 6 nm2, agreeing with the value obtained in the two-state analysis. Assuming a cylindrical channel, the dimensions of the (fully open) pore were comparable to DeltaA. Thus, the tension dependence of channel gating is primarily one of increasing the external channel area to accommodate the pore of the smallest conducting state. The higher conducting states appear to involve conformational changes internal to the channel that don't involve changes in area.

Studies of DNA dumbbells. VI. Analysis of optical melting curves of dumbbells with a sixteen-base pair duplex stem and end-loops of variable size and sequence.

Optical melting curves of 22 DNA dumbbells with the 16-base pair duplex sequence 5'-G-C-A-T-C-A-T-C-G-A-T-G-A-T-G-C-3' linked on both ends by single-strand loops of A, or C, sequences (iota = 2, 3, 4, 6, 8, 10, 14). T sequences (iota = 2, 3, 4, 6, 8, 10), and G iota sequences (iota = 2, 4) were measured in phosphate buffered solvents containing 30, 70, and 120 mM Na+. For dumbbells with loops comprised of at least three nucleotides, stability is inversely proportional to end-loop size. Dumbbells with loops comprised of only two nucleotide bases generally have lower stabilities than dumbbells with three base nucleotide loops. Experimental melting curves were analyzed in terms of the numerically exact (multistate) statistical thermodynamic model of DNA dumbbell melting previously described (T. M. paner, M. Amaratunga & A. S. Benight (1992), Biopolymers 32, 881). Theoretically calculated melting curves were fitted to experimental curves by simultaneously adjusting model parameters representing statistical weights of intramolecular hairpin loop and single-strand circle states. The systematically determined empirical parameters provided evaluations of the energetics of hairpin loop formation as a function of loop size, sequence, and salt environment. Values of the free energies of hairpin loop formation delta Gloop(n > iota) and single-strand circles, delta Gcir(N) as a function of end-loop size, tau = 2-14, circle size, N = 32 + 2 iota, and loop sequence were obtained. These quantities were found to depend on end-loop size but not loop sequence. Their empirically determined values also varied with solvent ionic strength. Analytical expressions for the partition function Q(T) of the dumbbells were evaluated using the empirically evaluated best-fit loop parameters. From Q(T), the melting transition enthalpy delta H, entropy delta S, and free energy delta G, were evaluated for the dumbbells as a function of end-loop size, sequence, and [Na+]. Since the multistate analysis is based on the numerically exact model, and considers a statistically significant number of theoretically possible partially melted states, it does not require prior assumptions regarding the nature of the melting transition, i.e., whether or not it occurs in a two-state manner. For comparison with the multistate analysis, thermodynamic transition parameters were also evaluated directly from experimental melting curves assuming a two-state transition and using the graphical van't Hoff analysis. Comparisons between results of the multistate and two-state analyses suggested dumbbells with loops comprised of six or fewer residues melted in a two-state manner, while the melting processes for dumbbells with larger end-loops deviate from two-state behavior. Dependence of thermodynamic transition parameters on [Na+] as a function of loop size suggests single-strand end-loops have different counterion binding properties than the melted circle. Results are compared with those obtained in an earlier study of dumbbells with the slightly different stem sequence 5'-G-C-A-T-A-G-A-T-G-A-G-A-A-T-G-C-3' linked on the ends by T iota loops (iota = 2, 3, 4, 6, 8, 10, 14).

Studies of DNA dumbbells. VI. Analysis of optical melting curves of dumbbells with a sixteen-base pair duplex stem and end-loops of variable size and sequence

Optical melting curves of 22 DNA dumbbells with the 16-base pair duplex sequence 5'-G-C-A-T-C-A-T-C-G-A-T-G-A-T-G-C-3' linked on both ends by single-strand loops of A, or C, sequences (iota = 2, 3, 4, 6, 8, 10, 14). T sequences (iota = 2, 3, 4, 6, 8, 10), and G iota sequences (iota = 2, 4) were measured in phosphate buffered solvents containing 30, 70, and 120 mM Na+. For dumbbells with loops comprised of at least three nucleotides, stability is inversely proportional to end-loop size. Dumbbells with loops comprised of only two nucleotide bases generally have lower stabilities than dumbbells with three base nucleotide loops. Experimental melting curves were analyzed in terms of the numerically exact (multistate) statistical thermodynamic model of DNA dumbbell melting previously described (T. M. paner, M. Amaratunga & A. S. Benight (1992), Biopolymers 32, 881). Theoretically calculated melting curves were fitted to experimental curves by simultaneously adjusting model parameters representing statistical weights of intramolecular hairpin loop and single-strand circle states. The systematically determined empirical parameters provided evaluations of the energetics of hairpin loop formation as a function of loop size, sequence, and salt environment. Values of the free energies of hairpin loop formation delta Gloop(n > iota) and single-strand circles, delta Gcir(N) as a function of end-loop size, tau = 2-14, circle size, N = 32 + 2 iota, and loop sequence were obtained. These quantities were found to depend on end-loop size but not loop sequence. Their empirically determined values also varied with solvent ionic strength. Analytical expressions for the partition function Q(T) of the dumbbells were evaluated using the empirically evaluated best-fit loop parameters. From Q(T), the melting transition enthalpy delta H, entropy delta S, and free energy delta G, were evaluated for the dumbbells as a function of end-loop size, sequence, and [Na+]. Since the multistate analysis is based on the numerically exact model, and considers a statistically significant number of theoretically possible partially melted states, it does not require prior assumptions regarding the nature of the melting transition, i.e., whether or not it occurs in a two-state manner. For comparison with the multistate analysis, thermodynamic transition parameters were also evaluated directly from experimental melting curves assuming a two-state transition and using the graphical van't Hoff analysis. Comparisons between results of the multistate and two-state analyses suggested dumbbells with loops comprised of six or fewer residues melted in a two-state manner, while the melting processes for dumbbells with larger end-loops deviate from two-state behavior. Dependence of thermodynamic transition parameters on [Na+] as a function of loop size suggests single-strand end-loops have different counterion binding properties than the melted circle. Results are compared with those obtained in an earlier study of dumbbells with the slightly different stem sequence 5'-G-C-A-T-A-G-A-T-G-A-G-A-A-T-G-C-3' linked on the ends by T iota loops (iota = 2, 3, 4, 6, 8, 10, 14).

Doorway state enhanced intramolecular vibrational energy redistribution in the asymmetric =CH2 hydride stretch of methyl vinyl ether

The microwave-infrared double-resonance capabilities of an electric-resonance optothermal spectrometer have been used to assign the high resolution (5 MHz) infrared spectrum of the asymmetric =CH2 stretch of the cis conformer of methyl vinyl ether near 3130 cm−1. This vibrational state is anharmonically coupled to a near-resonant bath state by a 0.69 cm−1 matrix element resulting in two vibrational bands separated by about 1.44 cm−1. The two mixed states resulting from this interaction are further coupled to other near-resonant bath states with an average matrix element of about 0.01 cm−1. The coupled state density increases weakly with the total angular momentum, J, however, the intramolecular vibrational energy distribution (IVR) rate is approximately independent of the total angular momentum quantum number. Therefore, the rotationally mediated coupling mechanisms are weaker than the anharmonic terms in the redistribution process. A two-state analysis of the strong coupling, which includes a phenomenological IVR rate constant, suggests that the IVR rate in the two mixed states is dominated by the contribution from the coupled dark state. From the deconvolution of the IVR rates to remove the contribution from the dark state, the IVR lifetime of the asymmetric =CH2 stretch is determined to be 660 ps.

Na+ dependence of single-channel current and channel density generate saturation of Na+ uptake in A6 cells.

In high-resistance, salt-absorbing epithelia the apical amiloride-sensitive Na+ channel is the key site for regulation of salt and water balance. The saturation of macroscopic Na+ transport through these channels was investigated using A6 epithelial monolayers. The relation between transepithelial Na+ transport (INa) and apical Na+ concentration ([Na+]ap) under short-circuit conditions was studied. Michaelis-Menten analysis of the saturable short-circuit current (Isc) yielded an apparent Michaelis-Menten constant (KmI) of 5 mmol/l and a maximal current (Imax) of 8 microA/cm2. The microscopic parameters underlying INa, namely the single-channel current (i) and the open channel density (No), were investigated by the analysis of current fluctuations induced by the electroneutral amiloride analogue CDPC (6-chloro-3, 5-diaminopyrazine-2-carboxamide). A two-state model analysis yielded the absolute values of i (0.18 +/- 0.01 pA) and No (65.38 +/- 9.57 million channels/cm2 of epithelium) at [Na+]ap = 110 mmol/l containing 50 mumol/l CDPC. Our data indicate that in A6 cells both i and No depend on [Na+]ap. Between 3 and approximately 20 mmol/l the density of conducting pores, No, decreases sharply and behaves again as an almost [Na+]ap-independent parameter at higher [Na+]ap. The single-channel current clearly saturates with an apparent Michaelis-Menten constant, Kmi, of approximately 17 mmol/l. Thus, the [Na+]ap dependence of No as well as the limited transport capacity of the amiloride-sensitive Na+ channel are both responsible for the saturation of INa.

ALASKA: A Mathematica package for two-state kinetic analysis of protein folding reactions

A Mathematica package (ALASKA) has been developed to simplify the measurement of protein folding kinetics by analysis of 1H NMR lineshape analysis. This package reads NMR data in ASCII format and can simulate an aromatic 1 NMR spectrum with or without lineshape broadening from chemical exchange. We describe the analysis of a urea denaturation series of a fast-folding protein, the G46A/G48A variant of monomeric λ repressor.