Extension of the thermodynamics of small systems to open metastable states: An example

By using a simplified model of small open liquid-like clusters with surface effects, in the gas phase, it is shown how the statistical thermodynamics of small systems can be extended to include metastable supersaturated gaseous states not too far from the gas–liquid equilibrium transition point. To accomplish this, one has to distinguish between mathematical divergence and physical convergence of the open-system partition function.

Folding thermodynamics of a model three-helix-bundle protein

The calculated folding thermodynamics of a simple off-lattice three-helix-bundle protein model under equilibrium conditions shows the experimentally observed protein transitions: a collapse transition, a disordered-to-ordered globule transition, a globule to native-state transition, and the transition from the active native state to a frozen inactive state. The cooperativity and physical origin of the various transitions are explored with a single “optimization” parameter and characterized with the Lindemann criterion for liquid versus solid-state dynamics. Below the folding temperature, the model has a simple free energy surface with a single basin near the native state; the surface is similar to that calculated from a simulation of the same three-helix-bundle protein with an all-atom representation [Boczko, E. M. & Brooks III, C. L. (1995) Science 269, 393–396].

A meanfield approach to the thermodynamics of a protein-solvent system with application to the oligomerization of the tumor suppressor p53

The thermodynamic stability and oligomerization status of the tumor suppressor p53 tetramerization domain have been studied experimentally and theoretically. A series of hydrophilic mutations at Met-340 and Leu-344 of human p53 were designed to disrupt the hydrophobic dimer–dimer interface of the tetrameric oligomerization domain of p53 (residues 325–355). Meanfield calculations of the free energy of the solvated mutants as a function of interdimer distance were compared with experimental data on the thermal stability and oligomeric state (tetramer, dimer, or equilibrium mixture of both) of each mutant. The calculations predicted a decreasing stability and oligomeric state for the following amino acids at residue 340: Met (tetramer) > Ser Asp, His, Gln, > Glu, Lys (dimer), whereas the experimental results showed the following order: Met (tetramer) > Ser > Gln > His, Lys > Asp, Glu (dimers). For residue 344, the calculated trend was Leu (tetramer) > Ala > Arg, Gln, Lys (dimer), and the experimental trend was Leu (tetramer) > Ala, Arg, Gln, Lys (dimer). The discrepancy for the lysine side chain at residue 340 is attributed to the dual nature of lysine, both hydrophobic and charged. The incorrect prediction of stability of the mutant with Asp at residue 340 is attributed to the fact that within the meanfield approach, we use the wild-type backbone configuration for all mutants, but low melting temperatures suggest a softening of the α-helices at the dimer–dimer interface. Overall, this initial application of meanfield theory toward a protein-solvent system is encouraging for the application of the theoretical model to more complex systems.

Brownian flocculation of polymer colloids in the presence of a secondary minimum

Most analyses of Brownian flocculation apply to conditions where London–van der Waals attractive forces cause particles to be strongly bound in a deep interparticle potential well. In this paper, results are reported that show the interaction between primary- and secondary-minimum flocculation when the interparticle potential curve reflects both attractive and electrostatic repulsive forces. The process is highly time-dependent because of transfer of particles from secondary- to primary-minimum flocculation. Essential features of the analysis are corroborated by experiments with 0.80-μm polystyrene spheres suspended in aqueous solutions of NaCl over a range of ionic strengths. In all cases, experiments were restricted to the initial stage of coagulation, where singlets and doublets predominate.

Kinetics and thermodynamics of T cell receptor– autoantigen interactions in murine experimental autoimmune encephalomyelitis

In the current study, cellular and molecular approaches have been used to analyze the biophysical nature of T cell receptor (TCR)–peptide MHC (pMHC) interactions for two autoreactive TCRs. These two TCRs recognize the N-terminal epitope of myelin basic protein (MBP1–11) bound to the MHC class II protein, I-Au, and are associated with murine experimental autoimmune encephalomyelitis. Mice transgenic for the TCRs have been generated and characterized in other laboratories. These analyses indicate that the mice either develop encephalomyelitis spontaneously (172.10 TCR) or only if immunized with autoantigen in adjuvant (1934.4 TCR). Here, we show that the 172.10 TCR binds MBP1–11:I-Au with a 4–5-fold higher affinity than the 1934.4 TCR. Consistent with the higher affinity, 172.10 T hybridoma cells are significantly more responsive to autoantigen than 1934.4 cells. The interaction of the 172.10 TCR with cognate ligand is more entropically unfavorable than that of the 1934.4 TCR, indicating that the 172.10 TCR undergoes greater conformational rearrangements upon ligand binding. The studies therefore suggest a correlation between the strength and plasticity of a TCR–pMHC interaction and the frequency of spontaneous disease in the corresponding TCR transgenic mice. The comparative analysis of these two TCRs has implications for understanding autoreactive T cell recognition and activation.

Large electrostatic differences in the binding thermodynamics of a cationic peptide to oligomeric and polymeric DNA.

Results presented here demonstrate that the thermodynamics of oligocation binding to polymeric and oligomeric DNA are not equivalent because of long-range electrostatic effects. At physiological cation concentrations (0.1-0.3 M) the binding of an oligolysine octacation KWK6-NH2 (+8 charge) to single-stranded poly(dT) is much stronger per site and significantly more salt concentration dependent than the binding of the same ligand to an oligonucleotide, dT(pdT)10 (-10 charge). These large differences are consistent with Poisson-Boltzmann calculations for a model that characterizes the charge distributions with key preaveraged structural parameters. Therefore, both the experimental and the theoretical results presented here show that the polyelectrolyte character of a polymeric nucleic acid makes a large contribution to both the magnitude and the salt concentration dependence of its binding interactions with simple oligocationic ligands.

Polymer Models of Meiotic and Mitotic Chromosomes

Polymers tied together by constraints exhibit an internal pressure; this idea is used to analyze physical properties of the bottle-brush–like chromosomes of meiotic prophase that consist of polymer-like flexible chromatin loops, attached to a central axis. Using a minimal number of experimental parameters, semiquantitative predictions are made for the bending rigidity, radius, and axial tension of such brushes, and the repulsion acting between brushes whose bristles are forced to overlap. The retraction of lampbrush loops when the nascent transcripts are stripped away, the oval shape of diplotene bivalents between chiasmata, and the rigidity of pachytene chromosomes are all manifestations of chromatin pressure. This two-phase (chromatin plus buffer) picture that suffices for meiotic chromosomes has to be supplemented by a third constituent, a chromatin glue to understand mitotic chromosomes, and explain how condensation can drive the resolution of entanglements. This process resembles a thermal annealing in that a parameter (the affinity of the glue for chromatin and/or the affinity of the chromatin for buffer) has to be tuned to achieve optimal results. Mechanical measurements to characterize this protein–chromatin matrix are proposed. Finally, the propensity for even slightly chemically dissimilar polymers to phase separate (cluster like with like) can explain the apparent segregation of the chromatin into A+T- and G+C-rich regions revealed by chromosome banding.

Folding protein models with a simple hydrophobic energy function: The fundamental importance of monomer inside/outside segregation

The present study explores a “hydrophobic” energy function for folding simulations of the protein lattice model. The contribution of each monomer to conformational energy is the product of its “hydrophobicity” and the number of contacts it makes, i.e., E(h⃗, c⃗) = −Σi=1N cihi = −(h⃗.c⃗) is the negative scalar product between two vectors in N-dimensional cartesian space: h⃗ = (h1, … , hN), which represents monomer hydrophobicities and is sequence-dependent; and c⃗ = (c1, … , cN), which represents the number of contacts made by each monomer and is conformation-dependent. A simple theoretical analysis shows that restrictions are imposed concomitantly on both sequences and native structures if the stability criterion for protein-like behavior is to be satisfied. Given a conformation with vector c⃗, the best sequence is a vector h⃗ on the direction upon which the projection of c⃗ − c̄⃗ is maximal, where c̄⃗ is the diagonal vector with components equal to c̄, the average number of contacts per monomer in the unfolded state. Best native conformations are suggested to be not maximally compact, as assumed in many studies, but the ones with largest variance of contacts among its monomers, i.e., with monomers tending to occupy completely buried or completely exposed positions. This inside/outside segregation is reflected on an apolar/polar distribution on the corresponding sequence. Monte Carlo simulations in two dimensions corroborate this general scheme. Sequences targeted to conformations with large contact variances folded cooperatively with thermodynamics of a two-state transition. Sequences targeted to maximally compact conformations, which have lower contact variance, were either found to have degenerate ground state or to fold with much lower cooperativity.

Dynamics of fluorescence fluctuations in green fluorescent protein observed by fluorescence correlation spectroscopy

We have investigated the pH dependence of the dynamics of conformational fluctuations of green fluorescent protein mutants EGFP (F64L/S65T) and GFP-S65T in small ensembles of molecules in solution by using fluorescence correlation spectroscopy (FCS). FCS utilizes time-resolved measurements of fluctuations in the molecular fluorescence emission for determination of the intrinsic dynamics and thermodynamics of all processes that affect the fluorescence. Fluorescence excitation of a bulk solution of EGFP decreases to zero at low pH (pKa = 5.8) paralleled by a decrease of the absorption at 488 nm and an increase at 400 nm. Protonation of the hydroxyl group of Tyr-66, which is part of the chromophore, induces these changes. When FCS is used the fluctuations in the protonation state of the chromophore are time resolved. The autocorrelation function of fluorescence emission shows contributions from two chemical relaxation processes as well as diffusional concentration fluctuations. The time constant of the fast, pH-dependent chemical process decreases with pH from 300 μs at pH 7 to 45 μs at pH 5, while the time-average fraction of molecules in a nonfluorescent state increases to 80% in the same range. A second, pH-independent, process with a time constant of 340 μs and an associated fraction of 13% nonfluorescent molecules is observed between pH 8 and 11, possibly representing an internal proton transfer process and associated conformational rearrangements. The FCS data provide direct measures of the dynamics and the equilibrium properties of the protonation processes. Thus FCS is a convenient, intrinsically calibrated method for pH measurements in subfemtoliter volumes with nanomolar concentrations of EGFP.

Branched co-polymers of histidine and lysine are efficient carriers of plasmids

We previously determined that a linear co-polymer of histidine and lysine (HK) in combination with liposomes enhanced the transfection efficiency of cationic liposomes. In the current study, we designed a series of HK polymers with increased branching and/or histidine/lysine ratio to determine if either variable affects transfection efficiency. In the presence of liposomes, the branched polymer with the highest number of histidines, HHK4b, was the most effective at enhancing gene expression. Furthermore, when serum was added to the medium during transfection, the combination of HHK4b and liposomes as a gene-delivery vehicle increased luciferase expression 400-fold compared to liposomes alone. In contrast to linear HK polymers, the higher branched HHK polymers were effective carriers of plasmids in the absence of liposomes. Without liposomes, the HHK4b carrier enhanced luciferase expression 15-fold in comparison with the lesser branched HHK2b carrier and increased expression by 5-logs in comparison with the HHK or HK carrier. The interplay of several parameters including increased condensation of DNA, buffering of acidic endosomes and differential binding affinities of polymer with DNA have a role in the enhancement of transfection by the HK polymers. In addition to suggesting that branched HK polymers are promising gene-delivery vehicles, this study provides a framework for the development of more efficient peptide-bond-based polymers of histidine and lysine.

Binding of the cellulose-binding domain of exoglucanase Cex from Cellulomonas fimi to insoluble microcrystalline cellulose is entropically driven.

Isothermal titration microcalorimetry is combined with solution-depletion isotherm data to analyze the thermodynamics of binding of the cellulose-binding domain (CBD) from the beta-1,4-(exo)glucanase Cex of Cellulomonas fimi to insoluble bacterial microcrystalline cellulose. Analysis of isothermal titration microcalorimetry data against two putative binding models indicates that the bacterial microcrystalline cellulose surface presents two independent classes of binding sites, with the predominant high-affinity site being characterized by a Langmuir-type Ka of 6.3 (+/-1.4) x 10(7) M-1 and the low-affinity site by a Ka of 1.1 (+/-0.6) x 10(6) M-1. CBDCex binding to either site is exothermic, but is mainly driven by a large positive change in entropy. This differs from protein binding to soluble carbohydrates, which is usually driven by a relatively large exothermic standard enthalpy change for binding. Differential heat capacity changes are large and negative, indicating that sorbent and protein dehydration effects make a dominant contribution to the driving force for binding.

Thermodynamic parameters for the helix-coil transition of oligopeptides: molecular dynamics simulation with the peptide growth method.

The helix-coil transition equilibrium of polypeptides in aqueous solution was studied by molecular dynamics simulation. The peptide growth simulation method was introduced to generate dynamic models of polypeptide chains in a statistical (random) coil or an alpha-helical conformation. The key element of this method is to build up a polypeptide chain during the course of a molecular transformation simulation, successively adding whole amino acid residues to the chain in a predefined conformation state (e.g., alpha-helical or statistical coil). Thus, oligopeptides of the same length and composition, but having different conformations, can be incrementally grown from a common precursor, and their relative conformational free energies can be calculated as the difference between the free energies for growing the individual peptides. This affords a straightforward calculation of the Zimm-Bragg sigma and s parameters for helix initiation and helix growth. The calculated sigma and s parameters for the polyalanine alpha-helix are in reasonable agreement with the experimental measurements. The peptide growth simulation method is an effective way to study quantitatively the thermodynamics of local protein folding.