Search

In addition to amide hydrogen bonds and the hydrophobic effect, interactions involving π-bonded sp 2 atoms of amides, aromatics, and other groups occur in protein self-assembly processes including folding, oligomerization, and condensate formation. These interactions also occur in aqueous solutions of amide and aromatic compounds, where they can be quantified. Previous analysis of thermodynamic coefficients quantifying net-favorable interactions of amide compounds with other amides and aromatics revealed that interactions of amide sp 2 O with amide sp 2 N unified atoms (presumably C═O···H-N hydrogen bonds) and amide/aromatic sp 2 C (lone pair π, n-π*) are particularly favorable. Sp 3 C-sp 3 C (hydrophobic), sp 3 C-sp 2 C (hydrophobic, CH-π), sp 2 C-sp 2 C (hydrophobic, π-π), and sp 3 C-sp 2 N interactions are favorable, sp 2 C-sp 2 N interactions are neutral, while sp 2 O-sp 2 O and sp 2 N-sp 2 N self-interactions and sp 2 O-sp 3 C interactions are unfavorable. Here, from determinations of favorable effects of 14 amides on naphthalene solubility at 10, 25, and 45 °C, we dissect amide-aromatic interaction free energies into enthalpic and entropic contributions and find these vary systematically with amide composition. Analysis of these results yields enthalpic and entropic contributions to intrinsic strengths of interactions of amide sp 2 O, sp 2 N, sp 2 C, and sp 3 C unified atoms with aromatic sp 2 C atoms. For each interaction, enthalpic and entropic contributions have the same sign and are much larger in magnitude than the interaction free energy itself. The amide sp 2 O-aromatic sp 2 C interaction is enthalpy-driven and entropically unfavorable, consistent with direct chemical interaction (e.g., lone pair-π), while amide sp 3 C- and sp 2 C-aromatic sp 2 C interactions are entropy-driven and enthalpically unfavorable, consistent with hydrophobic effects. These findings are relevant for interactions involving π-bonded sp 2 atoms in protein processes.

In addition to amide hydrogen bonds and the hydrophobic effect, interactions involving π-bonded sp 2 atoms of amides, aromatics and other groups occur in protein self-assembly processes including folding, oligomerization and condensate formation. These interactions also occur in aqueous solutions of amide and aromatic compounds, where they can be quantified. Previous analysis of thermodynamic coefficients quantifying net-favorable interactions of amide compounds with other amides and aromatics revealed that interactions of amide sp 2 O with amide sp 2 N unified atoms (presumably C=O···H-N hydrogen bonds) and amide/aromatic sp 2 C (lone pair-π, n-π * ) are particularly favorable. Sp 3 C-sp 3 C (hydrophobic), sp 3 C-sp 2 C (hydrophobic, CH-π), sp 2 C-sp 2 C (hydrophobic, π-π) and sp 3 C-sp 2 N interactions are favorable, sp 2 C-sp 2 N interactions are neutral, while sp 2 O-sp 2 O and sp 2 N-sp 2 N self-interactions and sp 2 O-sp 3 C interactions are unfavorable. Here, from determinations of favorable effects of fourteen amides on naphthalene solubility at 10, 25 and 45 °C, we dissect amide-aromatic interaction free energies into enthalpic and entropic contributions and find these vary systematically with amide composition. Analysis of these results yields enthalpic and entropic contributions to intrinsic strengths of interactions of amide sp 2 O, sp 2 N, sp 2 C and sp 3 C unified atoms with aromatic sp 2 C atoms. For each interaction, enthalpic and entropic contributions have the same sign and are much larger in magnitude than the interaction free energy itself. The amide sp 2 O-aromatic sp 2 C interaction is enthalpy-driven and entropically unfavorable, consistent with direct chemical interaction (e.g. lone pair-π) while amide sp 3 C- and sp 2 C-aromatic sp 2 C interactions are entropy-driven and enthalpically unfavorable, consistent with hydrophobic effects. These findings are relevant for interactions involving π-bonded sp 2 atoms in protein processes.

E. coli single-stranded-DNA binding protein (EcSSB) displays nearest-neighbor (NN) and non-nearest-neighbor (NNN)) cooperativity in binding ssDNA during genome maintenance. NNN cooperativity requires the intrinsically-disordered linkers (IDL) of the C-terminal tails. Potassium glutamate (KGlu), the primary E. coli salt, promotes NNN-cooperativity, while KCl inhibits it. We find that KGlu promotes compaction of a single polymeric SSB-coated ssDNA beyond what occurs in KCl, indicating a link of compaction to NNN-cooperativity. EcSSB also undergoes liquid-liquid phase separation (LLPS), inhibited by ssDNA binding. We find that LLPS, like NNN-cooperativity, is promoted by increasing [KGlu] in the physiological range, while increasing [KCl] and/or deletion of the IDL eliminate LLPS, indicating similar interactions in both processes. From quantitative determinations of interactions of KGlu and KCl with protein model compounds, we deduce that the opposing effects of KGlu and KCl on SSB LLPS and cooperativity arise from their opposite interactions with amide groups. KGlu interacts unfavorably with the backbone (especially Gly) and side chain amide groups of the IDL, promoting amide-amide interactions in LLPS and NNN-cooperativity. By contrast, KCl interacts favorably with these amide groups and therefore inhibits LLPS and NNN-cooperativity. These results highlight the importance of salt interactions in regulating the propensity of proteins to undergo LLPS., Competing Interests: Declaration of interests The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

Transcription initiation is highly regulated by promoter sequence, transcription factors, and ligands. All known transcription inhibitors, an important class of antibiotics, act in initiation. To understand regulation and inhibition, the biophysical mechanisms of formation and stabilization of the "open" promoter complex (OC), of synthesis of a short RNA-DNA hybrid upon nucleotide addition, and of escape of RNA polymerase (RNAP) from the promoter must be understood. We previously found that RNAP forms three different OC with λP R promoter DNA. The 37 °C RNAP-λP R OC (RP O ) is very stable. At lower temperatures, RP O is less stable and in equilibrium with an intermediate OC (I 3 ). Here, we report step-by-step rapid quench-flow kinetic data for initiation and growth of the RNA-DNA hybrid at 25 and 37 °C that yield rate constants for each step of productive nucleotide addition. Analyzed together, with previously published data at 19 °C, our results reveal that I 3 and not RP O is the productive initiation complex at all temperatures. From the strong variations of rate constants and activation energies and entropies for individual steps of hybrid extension, we deduce that contacts of RNAP with the bubble strands are disrupted stepwise as the hybrid grows and translocates. Stepwise disruption of RNAP-strand contacts is accompanied by stepwise bubble collapse, base stacking, and duplex formation, as the hybrid extends to a 9-mer prior to disruption of upstream DNA-RNAP contacts and escape of RNAP from the promoter., Competing Interests: The authors declare no competing interest.

Folding and other protein self-assembly processes are driven by favorable interactions between O, N, and C unified atoms of the polypeptide backbone and side chains. These processes are perturbed by solutes that interact with these atoms differently than water does. Amide NH···O=C hydrogen bonding and various π-system interactions have been better characterized structurally or by simulations than experimentally in water, and unfavorable interactions are relatively uncharacterized. To address this situation, we previously quantified interactions of alkyl ureas with amide and aromatic compounds, relative to interactions with water. Analysis yielded strengths of interaction of each alkylurea with unit areas of different hybridization states of unified O, N, and C atoms of amide and aromatic compounds. Here, by osmometry, we quantify interactions of 10 pairs of amides selected to complete this dataset. An analysis yields intrinsic strengths of six favorable and four unfavorable atom-atom interactions, expressed per unit area of each atom and relative to interactions with water. The most favorable interactions are sp 2 O-sp 2 C (lone pair-π, presumably n -π*), sp 2 C-sp 2 C (π-π and/or hydrophobic), sp 2 O-sp 2 N (hydrogen bonding) and sp 3 C-sp 2 C (CH-π and/or hydrophobic). Interactions of sp 3 C with itself (hydrophobic) and with sp 2 N are modestly favorable, while sp 2 N interactions with sp 2 N and with amide/aromatic sp 2 C are modestly unfavorable. Amide sp 2 O-sp 2 O interactions and sp 2 O-sp 3 C interactions are more unfavorable, indicating the preference of amide sp 2 O to interact with water. These intrinsic interaction strengths are used to predict interactions of amides with proteins and chemical effects of amides (including urea, N -ethylpyrrolidone [NEP], and polyvinylpyrrolidone [PVP]) on protein stability., Competing Interests: The authors declare no competing interest.

FRET (fluorescence resonance energy transfer) between far-upstream (-100) and downstream (+14) cyanine dyes (Cy3, Cy5) showed extensive bending and wrapping of λP R promoter DNA on Escherichia coli RNA polymerase (RNAP) in closed and open complexes (CC and OC, respectively). Here we determine the kinetics and mechanism of DNA bending and wrapping by FRET and of formation of RNAP contacts with -100 and +14 DNA by single-dye protein-induced fluorescence enhancement (PIFE). FRET and PIFE kinetics exhibit two phases: rapidly reversible steps forming a CC ensemble ({CC}) of four intermediates [initial (RP C ), early (I 1E ), mid (I 1M ), and late (I 1L )], followed by conversion of {CC} to OC via I 1L . FRET and PIFE are first observed for I 1E , not RP c . FRET and PIFE together reveal large-scale bending and wrapping of upstream and downstream DNA as RP C advances to I 1E , decreasing the Cy3-Cy5 distance to ∼75 Å and making RNAP-DNA contacts at -100 and +14. We propose that far-upstream DNA wraps on the upper β'-clamp while downstream DNA contacts the top of the β-pincer in I 1E . Converting I 1E to I 1M (∼1 s time scale) reduces FRET efficiency with little change in -100 or +14 PIFE, interpreted as clamp opening that moves far-upstream DNA (on β') away from downstream DNA (on β) to increase the Cy3-Cy5 distance by ∼14 Å. FRET increases greatly in converting I 1M to I 1L , indicating bending of downstream duplex DNA into the clamp and clamp closing to reduce the Cy3-Cy5 distance by ∼21 Å. In the subsequent rate-determining DNA-opening step, in which the clamp may also open, I 1L is converted to the initial unstable OC (I 2 ). Implications for facilitation of CC-to-OC isomerization by upstream DNA and upstream binding, DNA-bending transcription activators are discussed.

To determine the step-by-step kinetics and mechanism of transcription initiation and escape by E. coli RNA polymerase from the λP R promoter, we quantify the accumulation and decay of transient short RNA intermediates on the pathway to promoter escape and full-length (FL) RNA synthesis over a wide range of NTP concentrations by rapid-quench mixing and phosphorimager analysis of gel separations. Experiments are performed at 19 °C, where almost all short RNAs detected are intermediates in FL-RNA synthesis by productive complexes or end-products in nonproductive (stalled) initiation complexes and not from abortive initiation. Analysis of productive-initiation kinetic data yields composite second-order rate constants for all steps of NTP binding and hybrid extension up to the escape point (11-mer). The largest of these rate constants is for incorporation of UTP into the dinucleotide pppApU in a step which does not involve DNA opening or translocation. Subsequent steps, each of which begins with reversible translocation and DNA opening, are slower with rate constants that vary more than 10-fold, interpreted as effects of translocation stress on the translocation equilibrium constant. Rate constants for synthesis of 4- and 5-mer, 7-mer to 9-mer, and 11-mer are particularly small, indicating that RNAP-promoter interactions are disrupted in these steps. These reductions in rate constants are consistent with the previously determined ∼9 kcal cost of escape from λP R . Structural modeling and previous results indicate that the three groups of small rate constants correspond to sequential disruption of in-cleft, -10, and -35 interactions. Parallels to escape by T7 RNAP are discussed.

Alkylureas display hydrocarbon and amide groups, the primary functional groups of proteins. To obtain the thermodynamic information that is needed to analyze interactions of amides and proteins with nucleobases and nucleic acids, we quantify preferential interactions of alkylureas with nucleobases differing in the amount and composition of water-accessible surface area (ASA) by solubility assays. Using an established additive ASA-based analysis, we interpret these thermodynamic results to determine interactions of each alkylurea with five types of nucleobase unified atoms (carbonyl sp 2 O, amino sp 3 N, ring sp 2 N, methyl sp 3 C, and ring sp 2 C). All alkylureas interact favorably with nucleobase sp 2 C and sp 3 C atoms; these interactions become more favorable with an increasing level of alkylation of urea. Interactions with nucleobase sp 2 O are most favorable for urea, less favorable for methylurea and ethylurea, and unfavorable for dialkylated ureas. Contributions to overall alkylurea-nucleobase interactions from interactions with each nucleobase atom type are proportional to the ASA of that atom type with proportionality constant (interaction strength) α, as observed previously for urea. Trends in α-values for interactions of alkylureas with nucleobase atom types parallel those for corresponding amide compound atom types, offset because nucleobase α-values are more favorable. Comparisons between ethylated and methylated ureas show interactions of amide compound sp 3 C with nucleobase sp 2 C, sp 3 C, sp 2 N, and sp 3 N atoms are favorable while amide sp 3 C-nucleobase sp 2 O interactions are unfavorable. Strongly favorable interactions of urea with nucleobase sp 2 O but weakly favorable interactions with nucleobase sp 3 N indicate that amide sp 2 N-nucleobase sp 2 O and nucleobase sp 3 N-amide sp 2 O hydrogen bonding (NH···O═C) interactions are favorable while amide sp 2 N-nucleobase sp 3 N interactions are unfavorable. These favorable amide-nucleobase hydrogen bonding interactions are prevalent in specific protein-nucleotide complexes.

Significant, otherwise-unavailable information about mechanisms and transition states (TS) of protein folding and binding is obtained from solute effects on rate constants. Here we characterize TS for lac repressor(R)-lac operator(O) binding by analyzing effects of RO-stabilizing and RO-destabilizing solutes on association (ka) and dissociation (kd) rate constants. RO-destabilizing solutes (urea, KCl) reduce ka comparably (urea) or more than (KCl) they increase kd, demonstrating that they destabilize TS relative to reactants and RO, and that TS exhibits most of the Coulombic interactions between R and O. Strikingly, three solutes which stabilize RO by favoring burial/dehydration of amide oxygens and anionic phosphate oxygens all reduce kd without affecting ka significantly. The lack of stabilization of TS by these solutes indicates that O phosphates remain hydrated in TS and that TS preferentially buries aromatic carbons and amide nitrogens while leaving amide oxygens exposed. In our proposed mechanism, DNA-binding-domains (DBD) of R insert in major grooves of O pre-TS, forming most Coulombic interactions of RO and burying aromatic carbons. Nucleation of hinge helices creates TS, burying sidechain amide nitrogens. Post-TS, hinge helices assemble and the DBD-hinge helix-O-DNA module docks on core repressor, partially dehydrating phosphate oxygens and tightening all interfaces to form RO., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Quantitative information about amide interactions in water is needed to understand their contributions to protein folding and amide effects on aqueous processes and to compare with computer simulations. Here we quantify interactions of urea, alkylated ureas, and other amides by osmometry and amide-aromatic hydrocarbon interactions by solubility. Analysis of these data yields strengths of interaction of ureas and naphthalene with amide sp 2 O, amide sp 2 N, aliphatic sp 3 C, and amide and aromatic sp 2 C unified atoms in water. Interactions of amide sp 2 O with urea and naphthalene are favorable, while amide sp 2 O-alkylurea interactions are unfavorable, becoming more unfavorable with increasing alkylation. Hence, amide sp 2 O-amide sp 2 N interactions (proposed n-σ* hydrogen bond) and amide sp 2 O-aromatic sp 2 C (proposed n-π*) interactions are favorable in water, while amide sp 2 O-sp 3 C interactions are unfavorable. Interactions of all ureas with sp 3 C and amide sp 2 N are favorable and increase in strength with increasing alkylation, indicating favorable sp 3 C-amide sp 2 N and sp 3 C-sp 3 C interactions. Naphthalene results show that aromatic sp 2 C-amide sp 2 N interactions in water are unfavorable while sp 2 C-sp 3 C interactions are favorable. These results allow interactions of amide and hydrocarbon moieties and effects of urea and alkylureas on aqueous processes to be predicted or interpreted in terms of structural information. We predict strengths of favorable urea-benzene and N-methylacetamide interactions from experimental information to compare with simulations and indicate how amounts of hydrocarbon and amide surfaces buried in protein folding and other biopolymer processes and transition states can be determined from analysis of urea and diethylurea effects on equilibrium and rate constants.

To investigate roles of the discriminator and open complex (OC) lifetime in transcription initiation by Escherichia coli RNA polymerase (RNAP; α 2 ββ'ωσ 70 ), we compare productive and abortive initiation rates, short RNA distributions, and OC lifetime for the λP R and T7A1 promoters and variants with exchanged discriminators, all with the same transcribed region. The discriminator determines the OC lifetime of these promoters. Permanganate reactivity of thymines reveals that strand backbones in open regions of long-lived λP R -discriminator OCs are much more tightly held than for shorter-lived T7A1-discriminator OCs. Initiation from these OCs exhibits two kinetic phases and at least two subpopulations of ternary complexes. Long RNA synthesis (constrained to be single round) occurs only in the initial phase (<10 s), at similar rates for all promoters. Less than half of OCs synthesize a full-length RNA; the majority stall after synthesizing a short RNA. Most abortive cycling occurs in the slower phase (>10 s), when stalled complexes release their short RNA and make another without escaping. In both kinetic phases, significant amounts of 8-nt and 10-nt transcripts are produced by longer-lived, λP R -discriminator OCs, whereas no RNA longer than 7 nt is produced by shorter-lived T7A1-discriminator OCs. These observations and the lack of abortive RNA in initiation from short-lived ribosomal promoter OCs are well described by a quantitative model in which ∼1.0 kcal/mol of scrunching free energy is generated per translocation step of RNA synthesis to overcome OC stability and drive escape. The different length-distributions of abortive RNAs released from OCs with different lifetimes likely play regulatory roles.

Potassium glutamate (KGlu) is the primary Escherichia coli cytoplasmic salt. After sudden osmotic upshift, cytoplasmic KGlu concentration increases, initially because of water efflux and subsequently by K + transport and Glu - synthesis, allowing water uptake and resumption of growth at high osmolality. In vitro, KGlu ranks with Hofmeister salts KF and K 2 SO 4 in driving protein folding and assembly. Replacement of KCl by KGlu stabilizes protein-nucleic acid complexes. To interpret and predict KGlu effects on protein processes, preferential interactions of KGlu with 15 model compounds displaying six protein functional groups-sp 3 (aliphatic) C; sp 2 (aromatic, amide, carboxylate) C; amide and anionic (carboxylate) O; and amide and cationic N-were determined by osmometry or solubility assays. Analysis of these data yields interaction potentials (α-values) quantifying non-Coulombic chemical interactions of KGlu with unit area of these six groups. Interactions of KGlu with the 15 model compounds predicted from these six α-values agree well with experimental data. KGlu interactions with all carbon groups and with anionic (carboxylate) and amide oxygen are unfavorable, while KGlu interactions with cationic and amide nitrogen are favorable. These α-values, together with surface area information, provide quantitative predictions of why KGlu is an effective E. coli cytoplasmic osmolyte (because of the dominant effect of unfavorable interactions of KGlu with anionic and amide oxygens and hydrocarbon groups on the water-accessible surface of cytoplasmic biopolymers) and why KGlu is a strong stabilizer of folded proteins (because of the dominant effect of unfavorable interactions of KGlu with hydrocarbon groups and amide oxygens exposed in unfolding)., (Copyright © 2016 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

In vitro, replacing KCl with potassium glutamate (KGlu), the Escherichia coli cytoplasmic salt and osmolyte, stabilizes folded proteins and protein-nucleic acid complexes. To understand the chemical basis for these effects and rank Glu- in the Hofmeister anion series for protein unfolding, we quantify and interpret the strong stabilizing effect of KGlu on the ribosomal protein domain NTL9, relative to the effects of other stabilizers (KCl, KF, and K2SO4) and destabilizers (GuHCl and GuHSCN). GuHSCN titrations at 20 ° C, performed as a function of the concentration of KGlu or another salt and monitored by NTL9 fluorescence, are analyzed to obtain R-values quantifying the Hofmeister salt concentration (m3) dependence of the unfolding equilibrium constant K(obs) [r-value = −d ln K(obs)/dm3 = (1/RT) dΔG(obs) ° /dm3 = m-value/RT]. r-Values for both stabilizing K+ salts and destabilizing GuH+ salts are compared with predictions from model compound data. For two-salt mixtures, we find that contributions of stabilizing and destabilizing salts to observed r-values are additive and independent. At 20 ° C, we determine a KGlu r-value of 3.22 m(−1) and K2SO4, KF, KCl, GuHCl, and GuHSCN r-values of 5.38, 1.05, 0.64, −1.38, and −3.00 m(−1), respectively. The KGlu r-value represents a 25-fold (1.9 kcal) stabilization per molal KGlu added. KGlu is much more stabilizing than KF, and the stabilizing effect of KGlu is larger in magnitude than the destabilizing effect of GuHSCN. Interpretation of the data reveals good agreement between predicted and observed relative r-values and indicates the presence of significant residual structure in GuHSCN-unfolded NTL9 at 20 ° C.

Initial recognition of promoter DNA by RNA polymerase (RNAP) is proposed to trigger a series of conformational changes beginning with bending and wrapping of the 40-50 bp of DNA immediately upstream of the -35 region. Kinetic studies demonstrated that the presence of upstream DNA facilitates bending and entry of the downstream duplex (to +20) into the active site cleft to form an advanced closed complex (CC), prior to melting of ∼13 bp (-11 to +2), including the transcription start site (+1). Atomic force microscopy and footprinting revealed that the stable open complex (OC) is also highly wrapped (-60 to +20). To test the proposed bent-wrapped model of duplex DNA in an advanced RNAP-λP(R) CC and compare wrapping in the CC and OC, we use fluorescence resonance energy transfer (FRET) between cyanine dyes at far-upstream (-100) and downstream (+14) positions of promoter DNA. Similarly large intrinsic FRET efficiencies are observed for the CC (0.30 ± 0.07) and the OC (0.32 ± 0.11) for both probe orientations. Fluorescence enhancements at +14 are observed in the single-dye-labeled CC and OC. These results demonstrate that upstream DNA is extensively wrapped and the start site region is bent into the cleft in the advanced CC, reducing the distance between positions -100 and +14 on promoter DNA from >300 to <100 Å. The proximity of upstream DNA to the downstream cleft in the advanced CC is consistent with the proposed mechanism for facilitation of OC formation by upstream DNA.

Osmosensing transporters mediate osmolyte accumulation to forestall cellular dehydration as the extracellular osmolality increases. ProP is a bacterial osmolyte-H(+) symporter, a major facilitator superfamily member, and a paradigm for osmosensing. ProP activity is a sigmoid function of the osmolality. It is determined by the osmolality, not the magnitude or direction of the osmotic shift, in cells and salt-loaded proteoliposomes. The activation threshold varies directly with the proportion of anionic phospholipid in cells and proteoliposomes. The osmosensory mechanism was probed by varying the salt composition and concentration outside and inside proteoliposomes. Data analysis was based on the hypothesis that the fraction of maximal transporter activity at a particular luminal salt concentration reflects the proportion of ProP molecules in an active conformation. ProP attained the same activity at the same osmolality when diverse, membrane-impermeant salts were added to the external medium. Contributions of Coulombic and/or Hofmeister salt effects to ProP activation were examined by varying the luminal salt cation (K(+) and Na(+)) and anion (chloride, phosphate, and sulfate) composition and then systematically increasing the luminal salt concentration by increasing the external osmolality. ProP activity increased with the sixth power of the univalent cation concentration, independent of the type of anion. This indicates that salt activation of ProP is a Coulombic, cation effect resulting from salt cation accumulation and not site-specific cation binding. Possible origins of this Coulombic effect include folding or assembly of anionic cytoplasmic ProP domains, an increase in local membrane surface charge density, and/or the juxtaposition of anionic protein and membrane surfaces during activation.

Small and large PEGs greatly increase chemical potentials of globular proteins (μ2), thereby favoring precipitation, crystallization, and protein-protein interactions that reduce water-accessible protein surface and/or protein-PEG excluded volume. To determine individual contributions of PEG-protein chemical and excluded volume interactions to μ2 as functions of PEG molality m3 , we analyze published chemical potential increments μ23  = dμ2/dm3 quantifying unfavorable interactions of PEG (PEG200-PEG6000) with BSA and lysozyme. For both proteins, μ23 increases approximately linearly with the number of PEG residues (N3). A 1 molal increase in concentration of PEG -CH2 OCH2 - groups, for any chain-length PEG, increases μBSA by ∼2.7 kcal/mol and μlysozyme by ∼1.0 kcal/mol. These values are similar to predicted chemical interactions of PEG -CH2 OCH2 - groups with these protein components (BSA ∼3.3 kcal/mol, lysozyme ∼0.7 kcal/mol), dominated by unfavorable interactions with amide and carboxylate oxygens and counterions. While these chemical effects should be dominant for small PEGs, larger PEGS are expected to exhibit unfavorable excluded volume interactions and reduced chemical interactions because of shielding of PEG residues in PEG flexible coils. We deduce that these excluded volume and chemical shielding contributions largely compensate, explaining why the dependence of μ23 on N3 is similar for both small and large PEGs., (© 2015 Wiley Periodicals, Inc.)

In transcription initiation by Escherichia coli RNA polymerase (RNAP), initial binding to promoter DNA triggers large conformational changes, bending downstream duplex DNA into the RNAP cleft and opening 13bp to form a short-lived open intermediate (I2). Subsequent conformational changes increase lifetimes of λPR and T7A1 open complexes (OCs) by >10(5)-fold and >10(2)-fold, respectively. OC lifetime is a target for regulation. To characterize late conformational changes, we determine effects on OC dissociation kinetics of deletions in RNAP mobile elements σ(70) region 1.1 (σ1.1), β' jaw and β' sequence insertion 3 (SI3). In very stable OC formed by the wild type WT RNAP with λPR (RPO) and by Δσ1.1 RNAP with λPR or T7A1, we conclude that downstream duplex DNA is bound to the jaw in an assembly with SI3, and bases -4 to +2 of the nontemplate strand discriminator region are stably bound in a positively charged track in the cleft. We deduce that polyanionic σ1.1 destabilizes OC by competing for binding sites in the cleft and on the jaw with the polyanionic discriminator strand and downstream duplex, respectively. Examples of σ1.1-destabilized OC are the final T7A1 OC and the λPR I3 intermediate OC. Deleting σ1.1 and either β' jaw or SI3 equalizes OC lifetimes for λPR and T7A1. DNA closing rates are similar for both promoters and all RNAP variants. We conclude that late conformational changes that stabilize OC, like early ones that bend the duplex into the cleft, are primary targets of regulation, while the intrinsic DNA opening/closing step is not., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

Transcription initiation is a highly regulated step of gene expression. Here, we discuss the series of large conformational changes set in motion by initial specific binding of bacterial RNA polymerase (RNAP) to promoter DNA and their relevance for regulation. Bending and wrapping of the upstream duplex facilitates bending of the downstream duplex into the active site cleft, nucleating opening of 13 bp in the cleft. The rate-determining opening step, driven by binding free energy, forms an unstable open complex, probably with the template strand in the active site. At some promoters, this initial open complex is greatly stabilized by rearrangements of the discriminator region between the -10 element and +1 base of the nontemplate strand and of mobile in-cleft and downstream elements of RNAP. The rate of open complex formation is regulated by effects on the rapidly-reversible steps preceding DNA opening, while open complex lifetime is regulated by effects on the stabilization of the initial open complex. Intrinsic DNA opening-closing appears less regulated. This noncovalent mechanism and its regulation exhibit many analogies to mechanisms of enzyme catalysis.

Small solutes are useful probes of large conformational changes in RNA polymerase-promoter interactions and other biopolymer processes. In general, a large effect of a solute on an equilibrium constant (or rate constant) indicates a large change in water-accessible biopolymer surface area in the corresponding step (or transition state), resulting from conformational changes, interface formation, or both. Here, we describe nitrocellulose filter binding assays from series used to determine the urea dependence of open complex formation and dissociation with Escherichia coli RNA polymerase and phage λPR promoter DNA. Then, we describe the subsequent data analysis and interpretation of these solute effects.

Protein folding has been extensively studied, but many questions remain regarding the mechanism. Characterizing early unstable intermediates and the high-free-energy transition state (TS) will help answer some of these. Here, we use effects of denaturants (urea, guanidinium chloride) and temperature on folding and unfolding rate constants and the overall equilibrium constant as probes of surface area changes in protein folding. We interpret denaturant kinetic m-values and activation heat capacity changes for 13 proteins to determine amounts of hydrocarbon and amide surface buried in folding to and from TS, and for complete folding. Predicted accessible surface area changes for complete folding agree in most cases with structurally determined values. We find that TS is advanced (50-90% of overall surface burial) and that the surface buried is disproportionately amide, demonstrating extensive formation of secondary structure in early intermediates. Models of possible pre-TS intermediates with all elements of the native secondary structure, created for several of these proteins, bury less amide and hydrocarbon surface than predicted for TS. Therefore, we propose that TS generally has both the native secondary structure and sufficient organization of other regions of the backbone to nucleate subsequent (post-TS) formation of tertiary interactions. The approach developed here provides proof of concept for the use of denaturants and other solutes as probes of amount and composition of the surface buried in coupled folding and other large conformational changes in TS and intermediates in protein processes.

To quantify interactions of the osmolyte l-proline with protein functional groups and predict their effects on protein processes, we use vapor pressure osmometry to determine chemical potential derivatives dμ2/dm3 = μ23, quantifying the preferential interactions of proline (component 3) with 21 solutes (component 2) selected to display different combinations of aliphatic or aromatic C, amide, carboxylate, phosphate or hydroxyl O, and amide or cationic N surface. Solubility data yield μ23 values for four less-soluble solutes. Values of μ23 are dissected using an ASA-based analysis to test the hypothesis of additivity and obtain α-values (proline interaction potentials) for these eight surface types and three inorganic ions. Values of μ23 predicted from these α-values agree with the experiment, demonstrating additivity. Molecular interpretation of α-values using the solute partitioning model yields partition coefficients (Kp) quantifying the local accumulation or exclusion of proline in the hydration water of each functional group. Interactions of proline with native protein surfaces and effects of proline on protein unfolding are predicted from α-values and ASA information and compared with experimental data, with results for glycine betaine and urea, and with predictions from transfer free energy analysis. We conclude that proline stabilizes proteins because of its unfavorable interactions with (exclusion from) amide oxygens and aliphatic hydrocarbon surfaces exposed in unfolding and that proline is an effective in vivo osmolyte because of the osmolality increase resulting from its unfavorable interactions with anionic (carboxylate and phosphate) and amide oxygens and aliphatic hydrocarbon groups on the surface of cytoplasmic proteins and nucleic acids.

Urea destabilizes helical and folded conformations of nucleic acids and proteins, as well as protein-nucleic acid complexes. To understand these effects, extend previous characterizations of interactions of urea with protein functional groups, and thereby develop urea as a probe of conformational changes in protein and nucleic acid processes, we obtain chemical potential derivatives (μ23 = dμ2/dm3) quantifying interactions of urea (component 3) with nucleic acid bases, base analogues, nucleosides, and nucleotide monophosphates (component 2) using osmometry and hexanol-water distribution assays. Dissection of these μ23 values yields interaction potentials quantifying interactions of urea with unit surface areas of nucleic acid functional groups (heterocyclic aromatic ring, ring methyl, carbonyl and phosphate O, amino N, sugar (C and O); urea interacts favorably with all these groups, relative to interactions with water. Interactions of urea with heterocyclic aromatic rings and attached methyl groups (as on thymine) are particularly favorable, as previously observed for urea-homocyclic aromatic ring interactions. Urea m-values determined for double helix formation by DNA dodecamers near 25 °C are in the range of 0.72-0.85 kcal mol(-1)m(-1) and exhibit little systematic dependence on nucleobase composition (17-42% GC). Interpretation of these results using the urea interaction potentials indicates that extensive (60-90%) stacking of nucleobases in the separated strands in the transition region is required to explain the m-value. Results for RNA and DNA dodecamers obtained at higher temperatures, and literature data, are consistent with this conclusion. This demonstrates the utility of urea as a quantitative probe of changes in surface area (ΔASA) in nucleic acid processes.

Differences in kinetics of transcription initiation by RNA polymerase (RNAP) at different promoters tailor the pattern of gene expression to cellular needs. After initial binding, large conformational changes occur in promoter DNA and RNAP to form initiation-capable complexes. To understand the mechanism and regulation of transcription initiation, the nature and sequence of these conformational changes must be determined. Escherichia coli RNAP uses binding free energy to unwind and separate 13 base pairs of λP(R) promoter DNA to form the unstable open intermediate I(2), which rapidly converts to much more stable open complexes (I(3), RP(o)). Conversion of I(2) to RP(o) involves folding/assembly of several mobile RNAP domains on downstream duplex DNA. Here, we investigate effects of a 42-residue deletion in the mobile β' jaw (ΔJAW) and truncation of promoter DNA beyond +12 (DT+12) on the steps of initiation. We find that in stable ΔJAW open complexes the downstream boundary of hydroxyl radical protection shortens by 5-10 base pairs, as compared to wild-type (WT) complexes. Dissociation kinetics of open complexes formed with ΔJAW RNAP and/or DT+12 DNA resemble those deduced for the structurally uncharacterized intermediate I(3). Overall rate constants (k(a)) for promoter binding and DNA opening by ΔJAW RNAP are much smaller than for WT RNAP. Values of k(a) for WT RNAP with DT+12 and full-length λP(R) are similar, though contributions of binding and isomerization steps differ. Hence, the jaw plays major roles both early and late in RP(o) formation, while downstream DNA functions primarily as the assembly platform after DNA opening.

We investigate how the coulombic Gibbs free energy and salt ion association per phosphate charge of DNA oligomers vary with oligomer size ( i.e. number of charged residues ∣ Z D ∣) at 0.15 M univalent salt by non-linear Poisson Boltzmann (NLPB) analysis of all-atom DNA models. Calculations of these quantities ([Formula: see text], [Formula: see text]) are performed for short and long double-stranded (ds) and single-stranded (ss) DNA oligomers, ranging from 4 to 118 phosphates (ds) and from 2 to 59 phosphates (ss). Behaviors of [Formula: see text] and [Formula: see text] as functions of ∣ Z D ∣ provide a measure of the range of the coulombic end effect and determine the size of an oligomer at which an interior region with the properties (per charge) of the infinite-length polyelectrolyte first appears. This size (10-11 phosphates at each end for ds DNA and 6-9 for ss DNA at 0.15 M salt) is in close agreement with values obtained previously by Monte Carlo and NLPB calculations for cylindrical models of polyions, and by analysis of binding of oligocations to DNA oligomers. Differences in [Formula: see text] and in [Formula: see text] between ss and ds DNA are used to predict effects of oligomeric size and salt concentration on duplex stability in the vicinity of 0.15 M salt. Results of all-atom calculations are compared with results of less structurally detailed models and with experimental data.

To explain the large, opposite effects of urea and glycine betaine (GB) on stability of folded proteins and protein complexes, we quantify and interpret preferential interactions of urea with 45 model compounds displaying protein functional groups and compare with a previous analysis of GB. This information is needed to use urea as a probe of coupled folding in protein processes and to tune molecular dynamics force fields. Preferential interactions between urea and model compounds relative to their interactions with water are determined by osmometry or solubility and dissected using a unique coarse-grained analysis to obtain interaction potentials quantifying the interaction of urea with each significant type of protein surface (aliphatic, aromatic hydrocarbon (C); polar and charged N and O). Microscopic local-bulk partition coefficients K(p) for the accumulation or exclusion of urea in the water of hydration of these surfaces relative to bulk water are obtained. K(p) values reveal that urea accumulates moderately at amide O and weakly at aliphatic C, whereas GB is excluded from both. These results provide both thermodynamic and molecular explanations for the opposite effects of urea and glycine betaine on protein stability, as well as deductions about strengths of amide NH--amide O and amide NH--amide N hydrogen bonds relative to hydrogen bonds to water. Interestingly, urea, like GB, is moderately accumulated at aromatic C surface. Urea m-values for protein folding and other protein processes are quantitatively interpreted and predicted using these urea interaction potentials or K(p) values.

Initiation of RNA synthesis from DNA templates by RNA polymerase (RNAP) is a multi-step process, in which initial recognition of promoter DNA by RNAP triggers a series of conformational changes in both RNAP and promoter DNA. The bacterial RNAP functions as a molecular isomerization machine, using binding free energy to remodel the initial recognition complex, placing downstream duplex DNA in the active site cleft and then separating the nontemplate and template strands in the region surrounding the start site of RNA synthesis. In this initial unstable "open" complex the template strand appears correctly positioned in the active site. Subsequently, the nontemplate strand is repositioned and a clamp is assembled on duplex DNA downstream of the open region to form the highly stable open complex, RP(o). The transcription initiation factor, σ(70), plays critical roles in promoter recognition and RP(o) formation as well as in early steps of RNA synthesis., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Small solutes affect protein and nucleic acid processes because of favorable or unfavorable chemical interactions of the solute with the biopolymer surface exposed or buried in the process. Large solutes also exclude volume and affect processes where biopolymer molecularity and/or shape changes. Here, we develop an analysis to separate and interpret or predict excluded volume and chemical effects of a flexible coil polymer on a process. We report a study of the concentration-dependent effects of the full series from monomeric to polymeric PEG on intramolecular hairpin and intermolecular duplex formation by 12-nucleotide DNA strands. We find that chemical effects of PEG on these processes increase in proportion to the product of the amount of DNA surface exposed on melting and the amount of PEG surface that is accessible to this DNA, and these effects are completely described by two interaction terms that quantify the interactions between this DNA surface and PEG end and interior groups. We find that excluded volume effects, once separated from these chemical effects, are quantitatively described by the analytical theory of Hermans, which predicts the excluded volume between a flexible polymer and a rigid molecule. From this analysis, we show that at constant concentration of PEG monomer, increasing PEG size increases the excluded volume effect but decreases the chemical interaction effect, because in a large PEG coil a smaller fraction of the monomers are accessible to the DNA. Volume exclusion by PEG has a much larger effect on intermolecular duplex formation than on intramolecular hairpin formation.

Previous isothermal titration calorimetry (ITC) and Förster resonance energy transfer studies demonstrated that Escherichia coli HU(αβ) binds nonspecifically to duplex DNA in three different binding modes: a tighter-binding 34-bp mode that interacts with DNA in large (>34 bp) gaps between bound proteins, reversibly bending it by 140(o) and thereby increasing its flexibility, and two weaker, modestly cooperative small site-size modes (10 bp and 6 bp) that are useful for filling gaps between bound proteins shorter than 34 bp. Here we use ITC to determine the thermodynamics of these binding modes as a function of salt concentration, and we deduce that DNA in the 34-bp mode is bent around-but not wrapped on-the body of HU, in contrast to specific binding of integration host factor. Analyses of binding isotherms (8-bp, 15-bp, and 34-bp DNA) and initial binding heats (34-bp, 38-bp, and 160-bp DNA) reveal that all three modes have similar log-log salt concentration derivatives of the binding constants (Sk(i)) even though their binding site sizes differ greatly; the most probable values of Sk(i) on 34-bp DNA or larger DNA are -7.5±0.5. From the similarity of Sk(i) values, we conclude that the binding interfaces of all three modes involve the same region of the arms and saddle of HU. All modes are entropy-driven, as expected for nonspecific binding driven by the polyelectrolyte effect. The bent DNA 34-bp mode is most endothermic, presumably because of the cost of HU-induced DNA bending, while the 6-bp mode is modestly exothermic at all salt concentrations examined. Structural models consistent with the observed Sk(i) values are proposed., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Though opening of the start site (+1) region of promoter DNA is required for transcription by RNA polymerase (RNAP), surprisingly little is known about how and when this occurs in the mechanism. Early events at the lambdaP(R) promoter load this region of duplex DNA into the active site cleft of Escherichia coli RNAP, forming the closed, permanganate-unreactive intermediate I(1). Conversion to the subsequent intermediate I(2) overcomes a large enthalpic barrier. Is I(2) open? Here we create a burst of I(2) by rapidly destabilizing open complexes (RP(o)) with 1.1 M NaCl. Fast footprinting reveals that thymines at positions from -11 to +2 in I(2) are permanganate-reactive, demonstrating that RNAP opens the entire initiation bubble in the cleft in a single step. Rates of decay of all observed thymine reactivities are the same as the I(2) to I(1) conversion rate determined by filter binding. In I(2), permanganate reactivity of the +1 thymine on the template (t) strand is the same as the RP(o) control, whereas nontemplate (nt) thymines are significantly less reactive than in RP(o). We propose that: (i) the +1(t) thymine is in the active site in I(2); (ii) conversion of I(2) to RP(o) repositions the nt strand in the cleft; and (iii) movements of the nt strand are coupled to the assembly and DNA binding of the downstream clamp and jaw that occurs after DNA opening and stabilizes RP(o). We hypothesize that unstable open intermediates at the lambdaP(R) promoter resemble the unstable, transcriptionally competent open complexes formed at ribosomal promoters.

Transcription by all RNA polymerases (RNAPs) requires a series of large-scale conformational changes to form the transcriptionally competent open complex RP(o). At the lambdaP(R) promoter, Escherichia coli sigma(70) RNAP first forms a wrapped, closed 100 bp complex I(1). The subsequent step opens the entire DNA bubble, creating the relatively unstable (open) complex I(2). Additional conformational changes convert I(2) to the stable RP(o). Here we probe these events by dissecting the effects of Na(+) salts of Glu(-), F(-), and Cl(-) on each step in this critical process. Rapid mixing and nitrocellulose filter binding reveal that the binding constant for I(1) at 25 degrees C is approximately 30-fold larger in Glu(-) than in Cl(-) at the same Na(+) concentration, with the same log-log salt concentration dependence for both anions. In contrast, both the rate constant and equilibrium constant for DNA opening (I(1) to I(2)) are only weakly dependent on salt concentration, and the opening rate constant is insensitive to replacement of Cl(-) with Glu(-). These very small effects of salt concentration on a process (DNA opening) that is strongly dependent on salt concentration in solution may indicate that the backbones of both DNA strands interact with polymerase throughout the process and/or that compensation is present between ion uptake and release. Replacement of Cl(-) with Glu(-) or F(-) at 25 degrees C greatly increases the lifetime of RP(o) and greatly reduces its salt concentration dependence. By analogy to Hofmeister salt effects on protein folding, we propose that the excluded anions Glu(-) and F(-) drive the folding and assembly of the RNAP clamp/jaw domains in the conversion of I(2) to RP(o), while Cl(-) does not. Because the Hofmeister effect of Glu(-) or F(-) largely compensates for the destabilizing Coulombic effect of any salt on the binding of this assembly to downstream promoter DNA, RP(o) remains long-lived even at 0.5 M Na(+) in Glu(-) or F(-) salts. The observation that Esigma(70) RP(o) complexes are exceedingly long-lived at moderate to high Glu(-) concentrations argues that Esigma(70) RNAP does not dissociate from strong promoters in vivo when the cytoplasmic glutamate concentration increases during osmotic stress.

The majority ( approximately 70%) of surface buried in protein folding is hydrocarbon, whereas in DNA helix formation, the majority ( approximately 65%) of surface buried is relatively polar nitrogen and oxygen. Our previous quantification of salt exclusion from hydrocarbon (C) accessible surface area (ASA) and accumulation at amide nitrogen (N) and oxygen (O) ASA leads to a prediction of very different Hofmeister effects on processes that bury mostly polar (N, O) surface compared to the range of effects commonly observed for processes that bury mainly nonpolar (C) surface, e.g., micelle formation and protein folding. Here we quantify the effects of salts on folding of the monomeric DNA binding domain (DBD) of lac repressor (lac DBD) and on formation of an oligomeric DNA duplex. In accord with this prediction, no salt investigated has a stabilizing Hofmeister effect on DNA helix formation. Our ASA-based analyses of model compound data and estimates of the surface area buried in protein folding and DNA helix formation allow us to predict Hofmeister effects on these processes. We observe semiquantitative to quantitative agreement between these predictions and the experimental values, obtained from a novel separation of coulombic and Hofmeister effects. Possible explanations of deviations, including salt-dependent unfolded ensembles and interactions with other types of surface, are discussed.

To improve our understanding of the effects of small solutes on protein stability, we conducted atomistic simulations to quantitatively characterize the interactions between two broadly used small solutes, urea and glycine betaine (GB), and a triglycine peptide, which is a good model for a protein backbone. Multiple solute concentrations were analyzed, and each solute-peptide-water ternary system was studied with approximately 200-300 ns of molecular dynamics simulations with the CHARMM force field. The comparison between calculated preferential interaction coefficients (Gamma(23)) and experimentally measured values suggests that semiquantitative agreement with experiments can be obtained if care is exercised to balance interactions among the solute, protein, and water. On the other hand, qualitatively incorrect (i.e., wrong sign in Gamma(23)) results can be obtained if a solute model is constructed by directly taking parameters for chemically similar groups from an existing force field. Such sensitivity suggests that small solute thermodynamic data can be valuable in the development of accurate force field models of biomolecules. Further decomposition of Gamma(23) into group contributions leads to additional insights regarding the effects of small solutes on protein stability. For example, use of the CHARMM force field predicts that urea preferentially interacts with not only amide groups in the peptide backbone but also aliphatic groups, suggesting a role for these interactions in urea-induced protein denaturation; quantitatively, however, it is likely that the CHARMM force field overestimates the interaction between urea and aliphatic groups. The results with GB support a simple thermodynamic model that assumes additivity of preferential interaction between GB and various biomolecular surfaces.

Interactions of histidine with nucleic acid phosphates and histidine pK(a) shifts make important contributions to many protein-nucleic acid binding processes. To characterize these phenomena in simplified systems, we quantified binding of a histidine-containing model peptide HWKK ((+)NH(3)-His-Trp-Lys-Lys-NH(2)) and its lysine analogue KWKK ((+)NH(3)-Lys-Trp-Lys-Lys-NH(2)) to a single-stranded RNA model, polyuridylate (polyU), by changes in tryptophan fluorescence as a function of salt concentration and pH. For both HWKK and KWKK, equilibrium binding constants, K(obs), and magnitudes of log-log salt derivatives, SK(obs) identical with (partial differential logK(obs)/partial differential log[Na(+)]), decreased with increasing pH in the manner expected for a titration curve model in which deprotonation of the histidine and alpha-amino groups weakens binding and reduces its salt-dependence. Fully protonated HWKK and KWKK exhibit the same K(obs) and SK(obs) within uncertainty, and these SK(obs) values are consistent with limiting-law polyelectrolyte theory for +4 cationic oligopeptides binding to single-stranded nucleic acids. The pH-dependence of HWKK binding to polyU provides no evidence for pK(a) shifts nor any requirement for histidine protonation, in stark contrast to the thermodynamics of coupled protonation often seen for these cationic residues in the context of native protein structure where histidine protonation satisfies specific interactions (e.g., salt-bridge formation) within highly complementary binding interfaces. The absence of pK(a) shifts in our studies indicates that additional Coulombic interactions across the nonspecific-binding interface between RNA and protonated histidine or the alpha-amino group are not sufficient to promote proton uptake for these oligopeptides. We present our findings in the context of hydration models for specific vs nonspecific nucleic acid binding.

Noncovalent self-assembly of biopolymers is driven by molecular interactions between functional groups on complementary biopolymer surfaces, replacing interactions with water. Since individually these interactions are comparable in strength to interactions with water, they have been difficult to quantify. Solutes (osmolytes, denaturants) exert often large effects on these self-assembly interactions, determined in sign and magnitude by how well the solute competes with water to interact with the relevant biopolymer surfaces. Here, an osmometric method and a water-accessible surface area (ASA) analysis are developed to quantify and interpret the interactions of the remarkable osmolyte glycine betaine (GB) with molecular surfaces in water. We find that GB, lacking hydrogen bond donors, is unable to compete with water to interact with anionic and amide oxygens; this explains its effectiveness as an osmolyte in the Escherichia coli cytoplasm. GB competes effectively with water to interact with amide and cationic nitrogens (hydrogen bonding) and especially with aromatic hydrocarbon (cation-pi). The large stabilizing effect of GB on lac repressor-lac operator binding is predicted quantitatively from ASA information and shown to result largely from dehydration of anionic DNA phosphate oxygens in the protein-DNA interface. The incorporation of these results into theoretical and computational analyses will likely improve the ability to accurately model intra- and interprotein interactions. Additionally, these results pave the way for development of solutes as kinetic/mechanistic and thermodynamic probes of conformational changes and formation/disruption of molecular interfaces that occur in the steps of biomolecular self-assembly processes.

Bacterial RNA polymerase and a "sigma" transcription factor form an initiation-competent "open" complex at a promoter by disruption of about 14 base pairs. Strand separation is likely initiated at the highly conserved -11 A-T base pair. Amino acids in conserved region 2.3 of the main Escherichia coli sigma factor, sigma(70), are involved in this process, but their roles are unclear. To monitor the fates of particular bases upon addition of RNA polymerase, promoters bearing single substitutions of the fluorescent A-analog 2-aminopurine (2-AP) at -11 and two other positions in promoter DNA were examined. Evidence was obtained for an open intermediate on the pathway to open complex formation, in which these 2-APs are no longer stacked onto their neighboring bases. The tyrosine at residue 430 in region 2.3 of sigma(70) was shown to be involved in quenching the fluorescence of a 2-AP substituted at -11, presumably through a stacking interaction. These data refine the structural model for open complex formation and reveal a novel interaction involved in DNA melting by RNA polymerase.

Salts and uncharged solutes in aqueous solution exert effects on a wide range of processes in which large amounts ofbiopolymer surface are buried or exposed (folding/unfolding, complexation/dissociation, or precipitation/dissolution). A simple two-state solute partitioning model (SPM, where the solute is partitioned between the bulk and surface water) allows the interpretation and prediction of the thermodynamic effects of various uncharged solutes (e.g., urea, glycine betaine) on protein and nucleic acid processes in terms of structural information. The correlation of solute effects with various coarse-grained types of biopolymer surface exposed or buried in a process provides a novel probe for investigation of large-scale conformational changes. Solutes that are fully excluded from one or more types of biopolymer surface are useful to quantify changes in water of hydration of these surfaces in biopolymer processes. Additionally, application of the SPM to the analysis of non-Coulombic salt effects on various model processes provides an estimate for the hydration layer at surfaces and shows that ion effects are additive and independent of the nature of the counterion.

Little is known about the thermodynamic forces that drive the folding pathways of higher-order RNA structure. In this study, we employ calorimetric [isothermal titration calorimetry (ITC) and differential scanning calorimetry (DSC)] and spectroscopic (NMR and UV) methods to characterize the thermodynamics of the GAAA tetraloop-receptor interaction, utilizing a previously described bivalent construct. ITC studies indicate that the bivalent interaction is enthalpy driven and highly stable, with a binding constant (K(obs)) of 5.5x10(6) M(-1) and enthalpy (DeltaH(obs)(o)) of -33.8 kcal/mol at 45 degrees C in 20 mM KCl and 2 mM MgCl(2). Thus, we derive the DeltaH(obs)(o) for a single tetraloop-receptor interaction to be -16.9 kcal/mol at these conditions. UV absorbance data indicate that an increase in base stacking quality contributes to the enthalpy of complex formation. These highly favorable thermodynamics are consistent with the known critical role for the tetraloop-receptor motif in the folding of large RNAs. Additionally, a significant heat capacity change (DeltaC(p,obs)(o)) of -0.24 kcal mol(-1) K(-1) was determined by ITC. DSC and UV-monitored thermal denaturation experiments indicate that the bivalent tetraloop-receptor construct follows a minimally five-state unfolding pathway and suggest the observed DeltaC(p,obs)(o) for the interaction results from a temperature-dependent unbound receptor RNA structure.

Recently, surface spectroscopies and simulations have begun to characterize the non-uniform distributions of salt ions near macroscopic and molecular surfaces. The thermodynamic consequences of these non-uniform distributions determine the often-large ion-specific effects of Hofmeister salts on a very wide range of processes in water. For uncharged surfaces, where these nonuniform ion distributions are confined to the first few layers of water at the surface, a two-state approximation to the distributions of water and ions, called the salt ion partitioning model (SPM) has both molecular and thermodynamic signiicance. Here, we summarize SPM results quantifying the local accumulation of H + , exclusion of HO - , and range of partitioning behavior of Hofmeister anions and cations near macroscopic and molecular interfaces. These results provide a database to interpret or predict Hofmeister salt effects on aqueous processes in terms of structural information regarding amount and composition of the surface exposed or buried in these processes.

Escherichia coli HU(alphabeta), a major nucleoid-associated protein, organizes chromosomal DNA and facilitates numerous DNA transactions. Using isothermal titration calorimetry, fluorescence resonance energy transfer and a series of DNA lengths (8 bp, 15 bp, 34 bp, 38 bp and 160 bp) we established that HU(alphabeta) interacts with duplex DNA using three different nonspecific binding modes. Both the HU to DNA molar ratio ([HU]/[DNA]) and DNA length dictate the dominant HU binding mode. On sufficiently long DNA (> or =34 bp), at low [HU]/[DNA], HU populates a noncooperative 34 bp binding mode with a binding constant of 2.1+/-0.4x10(6) M(-1), and a binding enthalpy of +7.7+/-0.6 kcal/mol at 15 degrees C and 0.15 M Na(+). With increasing [HU]/[DNA], HU bound in the noncooperative 34 bp mode progressively converts to two cooperative (omega approximately 20) modes with site sizes of 10 bp and 6 bp. These latter modes exhibit smaller binding constants (1.1+/-0.2x10(5) M(-1) for the 10 bp mode, 3.5+/-1.4x10(4) M(-1) for the 6 bp mode) and binding enthalpies (4.2+/-0.3 kcal/mol for the 10 bp mode, -1.6+/-0.3 kcal/mol for the 6 bp mode). As DNA length increases to 34 bp or more at low [HU]/[DNA], the small modes are replaced by the 34 bp binding mode. Fluorescence resonance energy transfer data demonstrate that the 34 bp mode bends DNA by 143+/-6 degrees whereas the 6 bp and 10 bp modes do not. The model proposed in this study provides a novel quantitative and comprehensive framework for reconciling previous structural and solution studies of HU, including single molecule (force extension measurement), fluorescence, and electrophoretic gel mobility-shift assays. In particular, it explains how HU condenses or extends DNA depending on the relative concentrations of HU and DNA.

Quantitative interpretation and prediction of Hofmeister ion effects on protein processes, including folding and crystallization, have been elusive goals of a century of research. Here, a quantitative thermodynamic analysis, developed to treat noncoulombic interactions of solutes with biopolymer surface and recently extended to analyze the effects of Hofmeister salts on the surface tension of water, is applied to literature solubility data for small hydrocarbons and model peptides. This analysis allows us to obtain a minimum estimate of the hydration b1 (H2O A(-2)), of hydrocarbon surface and partition coefficients Kp, characterizing the distribution of salts and salt ions between this hydration water and bulk water. Assuming that Na+ and SO4(2-) ions of Na2SO4 (the salt giving the largest reduction in hydrocarbon solubility as well as the largest increase in surface tension) are fully excluded from the hydration water at hydrocarbon surface, we obtain the same b1 as for air-water surface (approximately 0.18 H2O A(-2)). Rank orders of cation and anion partition coefficients for nonpolar surface follow the Hofmeister series for protein processes, but are strongly offset for cations in the direction of exclusion (preferential hydration). By applying a coarse-grained decomposition of water accessible surface area (ASA) into nonpolar, polar amide, and other polar surface and the same hydration b1 to interpret peptide solubility increments, we determine salt partition coefficients for amide surface. These partition coefficients are separated into single-ion contributions based on the observation that both Cl- and Na+ (also K+) occupy neutral positions in the middle of the anion and cation Hofmeister series for protein folding. Independent of this assignment, we find that all cations investigated are strongly accumulated at amide surface while most anions are excluded. Cation and anion effects are independent and additive, allowing successful prediction of Hofmeister salt effects on micelle formation and other processes from structural information (ASA).

To characterize driving forces and driven processes in formation of a large-interface, wrapped protein-DNA complex analogous to the nucleosome, we have investigated the thermodynamics of binding the 34-base pair (bp) H' DNA sequence to the Escherichia coli DNA-remodeling protein integration host factor (IHF). Isothermal titration calorimetry and fluorescence resonance energy transfer are applied to determine effects of salt concentration [KCl, KF, K glutamate (KGlu)] and of the excluded solute glycine betaine (GB) on the binding thermodynamics at 20 degrees C. Both the binding constant K(obs) and enthalpy Delta H degrees (obs) depend strongly on [salt] and anion identity. Formation of the wrapped complex is enthalpy driven, especially at low [salt] (e.g., Delta H(o)(obs)=-20.2 kcal x mol(-1) in 0.04 M KCl). Delta H degrees (obs) increases linearly with [salt] with a slope (d Delta H degrees (obs)/d[salt]), which is much larger in KCl (38+/-3 kcal x mol(-1) M(-1)) than in KF or KGlu (11+/-2 kcal x mol(-1) M(-1)). At 0.33 M [salt], K(obs) is approximately 30-fold larger in KGlu or KF than in KCl, and the [salt] derivative SK(obs)=dlnK(obs)/dln[salt] is almost twice as large in magnitude in KCl (-8.8+/-0.7) as in KF or KGlu (-4.7+/-0.6). A novel analysis of the large effects of anion identity on K(obs), SK(obs) and on Delta H degrees (obs) dissects coulombic, Hofmeister, and osmotic contributions to these quantities. This analysis attributes anion-specific differences in K(obs), SK(obs), and Delta H degrees (obs) to (i) displacement of a large number of water molecules of hydration [estimated to be 1.0(+/-0.2)x10(3)] from the 5340 A(2) of IHF and H' DNA surface buried in complex formation, and (ii) significant local exclusion of F(-) and Glu(-) from this hydration water, relative to the situation with Cl(-), which we propose is randomly distributed. To quantify net water release from anionic surface (22% of the surface buried in complexation, mostly from DNA phosphates), we determined the stabilizing effect of GB on K(obs): dlnK(obs)/d[GB]=2.7+/-0.4 at constant KCl activity, indicating the net release of ca. 150 H(2)O molecules from anionic surface.

The formation of the transcriptionally competent open complex (RP(o)) by Escherichia coli RNA polymerase at the lambda P(R) promoter involves at least three steps and two kinetically significant intermediates (I(1) and I(2)). Understanding the sequence of conformational changes (rearrangements in the jaws of RNA polymerase, DNA opening) that occur in the conversion of I(1) to RP(o) requires: (1) dissecting the rate constant k(d) for the dissociation of RP(o) into contributions from individual steps and (2) isolating and characterizing I(2). To deconvolute k(d), we develop experiments involving rapid upshifts to elevated concentrations of RP(o)-destabilizing solutes ("perturbants": urea and KCl) to create a burst in the population of I(2). At high concentrations of either perturbant, k(d) approaches the same [perturbant]-independent value, interpreted as the elementary rate constant k(-2) for I(2)-->I(1). The large effects of [urea] and [salt] on K(3) (the equilibrium constant for I(2) is in equilibrium with RP(o)) indicate that a large-scale folding transition in polymerase occurs and a new interface with the DNA forms late in the mechanism. We deduce that I(2) at the lambda P(R) promoter is always unstable relative to RP(o), even at 0 degrees C, explaining previous difficulties in detecting it by using temperature downshifts. The division of the large positive enthalpy change between the late steps of the mechanism suggests that an additional unstable intermediate (I(3)) may exist between I(2) and RP(o).

Two fundamentally different thermodynamic approaches are in use to interpret or predict the effects of urea on biopolymer processes: one is a synthesis of transfer free energies obtained from measurements of the effects of urea on the solubilities of small, model compounds; the other is an analysis of preferential interactions of urea with a range of folded and unfolded biopolymer surfaces. Here, we compare the predictions of these two approaches for the contribution of urea-amide (peptide) interactions to destabilization of folded proteins by urea. For these comparisons, we develop independent thermodynamic analyses of osmometric and solubility data characterizing interactions of a model compound with urea (or any other solute) and apply them to all five model compounds (glycine, alanine, diglycine, glycylalanine, and triglycine) where both isopiestic distillation (ID) and solubility data in aqueous urea solutions are available. We use model-independent expressions to calculate mu ex 23, the derivative of the "excess" chemical potential of solute "2" (either a model compound or a biopolymer) with respect to the molality of solute "3" (urea). Analyses of ID data for these systems reveal significant dependences of mu ex 23 on both m2 and m3, which must be taken into account in making comparisons with values of mu ex 23 obtained from solubility studies or from analyses of urea-biopolymer preferential interactions. Values of mu ex 23 calculated from model compound ID data at low m2 and m3 are directly proportional to the amount of polar amide (N, O) surface area, and not to any other type of surface. The proportionality constant in this limit, mu ex 23 /(RT x ASA) = (1.0 +/- 0.1) x 10(-3) A(-2), is very similar to that previously obtained by analysis of urea-biopolymer preferential interactions ((1.4 +/- 0.3) x 10(-3) A(-2)). This level of agreement for amide surface in the low concentration limit, as well as the absence of any significant preferential interaction of urea with Gly and Ala, reinforces the conclusion that the primary preferential interaction of urea with protein surface is a favorable interaction (resulting in local accumulation of urea) at polar amide surface, located mostly on the peptide backbone. However, mu ex 23 for interactions of urea with these model amides is found from both ID and solubility data to be urea concentration-dependent, in contrast to the urea concentration independence of the analogous quantity for protein unfolding.

We apply a recently developed surface-bulk partitioning model to interpret the effects of individual Hofmeister cations and anions on the surface tension of water. The most surface-excluded salt (Na2SO4) provides a minimum estimate for the number of water molecules per unit area of the surface region of 0.2 H2O A-2. This corresponds to a lower bound thickness of the surface region of approximately 6 A, which we assume is a property of this region and not of the salt investigated. At salt concentrations < or = 1 m, single-ion partition coefficients Kp,i, defined relative to Kp,Na+ = Kp,SO42- = 0, are found to be independent of bulk salt concentration and additive for different salt ions. Semiquantitative agreement with surface-sensitive spectroscopy data and molecular dynamics simulations is attained. In most cases, the rank orders of Kp,i for both anions and cations follow the conventional Hofmeister series, qualitative rankings of ions based on their effects on protein processes (folding, precipitation, assembly). Most anions that favor processes that expose protein surface to water (e.g., SCN-), and hence must interact favorably with (i.e., accumulate at) protein surface, are also accumulated at the air-water interface (Kp >1, e.g., Kp,SCN- =1.6). Most anions that favor processes that remove protein surface from water (e.g., F-), and hence are excluded from protein surface, are also excluded from the air-water interface (Kp,F- = 0.5). The guanidinium cation, a strong protein denaturant and therefore accumulated at the protein surface exposed in unfolding, is somewhat excluded from the air-water surface (Kp,GuH+ = 0.7), but is much less excluded than alkali metal cations (e.g., Kp,Na+ identical with 0, Kp,K+ = 0.1). Hence, cation Kp values for the air-water surface appear shifted (toward exclusion) as compared with values inferred for interactions of these cations with protein surface.

The architecture of cellular RNA polymerases (RNAPs) dictates that transcription can begin only after promoter DNA bends into a deep channel and the start site nucleotide (+1) binds in the active site located on the channel floor. Formation of this transcriptionally competent "open" complex (RP(o)) by Escherichia coli RNAP at the lambdaP(R) promoter is greatly accelerated by DNA upstream of base pair -47 (with respect to +1). Here we report real-time hydroxyl radical (*OH) and potassium permanganate (KMnO4) footprints obtained under conditions selected for optimal characterization of the first kinetically significant intermediate (I(1)) in RP(o) formation. .OH footprints reveal that the DNA backbone from -71 to -81 is engulfed by RNAP in I(1) but not in RP(o); downstream protection extends to approximately +20 in both complexes. KMnO4 footprinting detects solvent-accessible thymine bases in RP(o), but not in I(1). We conclude that upstream DNA wraps more extensively on RNAP in I(1) than in RP(o) and that downstream DNA (-11 to +20) occupies the active-site channel in I(1) but is not yet melted. Mapping of the footprinting data onto available x-ray structures provides a detailed model of a kinetic intermediate in bacterial transcription initiation and suggests how transient contacts with upstream DNA in I(1) might rearrange the channel to favor entry of downstream duplex DNA.

In vitro changes in polymer volume fraction (macromolecular crowding) and changes in solute or salt concentration typically have large effects on protein and nucleic acid processes (e.g., folding, binding, assembly, precipitation, crystallization). However, the large changes in these concentration variables, which occur in vivo as part of cellular responses to osmotic stress, appear to have much less dramatic effects on cellular biopolymer processes. Methods of changing intracellular concentrations by varying the extracellular osmolality or the concentration of a permeable solute or by titrating cells with an impermeable solute (plasmolysis) under conditions where an active response is suppressed are reviewed. The first in vivo biophysical studies of protein folding and protein diffusion performed as a function of these variables are also discussed.

In this issue of Cell, Hsu et al. (2006) report on the binding activity of a variant of the bacterial transcriptional specificity factor sigma (sigma) to promoter DNA. This study demonstrates that the sigma variant induces a large distortion in the transcriptional start site in the absence of core RNA polymerase, raising intriguing new questions about the roles of sigma and core RNA polymerase in transcription initiation.

Recently, surface-sensitive spectroscopy data and molecular dynamics simulations have generated intense interest in the distribution of electrolyte ions between bulk water and the air/water interface. A partitioning model for cations and anions developed for biopolymer surface is extended here to interpret the effects of selected acids, bases, and salts on the surface tension of water. Data for electrolytes were analyzed by using a lower-bound value for the number of water molecules in the surface region [0.2 H(2)O A(-2) (approximately two layers of water)], obtained by assuming that both Na(+) and SO(4)(2-) (i.e., Na(2)SO(4)) are fully excluded from this region. Surface-bulk partition coefficients of atmospherically relevant anions and the proton are determined. Notably, we find that H(+) is most strongly surface-accumulated, I(-) is modestly accumulated, NO(3)(-) is evenly distributed, and OH(-) is weakly excluded.