Commitment to folded and aggregated states occurs late in interleukin-1 beta folding.

A point mutation, lysine 97 to isoleucine, in the all-beta cytokine interleukin-1 beta (IL-1 beta) exhibits an increased propensity to form inclusion bodies in vivo and aggregates in vitro. In an effort to better understand the aggregation reaction and determine when intervention may allow rescue of protein from aggregation during renaturation, we developed a novel application of mass spectrometry using isotopic labeling to determine the step(s) at which K97I commits to either the native or aggregated state. Interestingly, despite the early formation of a folding intermediate ensemble at an observed rate lambda(2) of 4.0 s(-1), K97I commits to folding at a significantly slower rate lambda(CF) of 0.021 s(-1). This rate of commitment to folding is in excellent agreement with the observed rate of K97I native state formation (lambda(1) = 0.018 s(-1)). K97I also commits slowly to aggregation at an observed rate lambda(CA) of 0.023 s(-1). Earlier folding species and aggregates present prior to these commitment steps are likely to be in a reversible equilibrium between monomeric folding intermediates and higher-order oligomers. Kinetic and equilibrium experimental measurements of folding and aggregation processes are consistent with a nucleation-dependent model of aggregation.

Aggregation events occur prior to stable intermediate formation during refolding of interleukin 1beta.

A point mutation, lysine 97 --> isoleucine (K97I), in a surface loop in the beta-sheet protein interleukin 1beta (IL-1beta), exhibits increased levels of inclusion body (IB) formation relative to the wild-type protein (WT) when expressed in Escherichia coli. Despite the common observation that less stable proteins are often found in IBs, K97I is more stable than WT. We examined the folding pathway of the mutant and wild-type proteins at pH 6.5 and 25 degrees C with manual-mixing and stopped-flow optical spectroscopy to determine whether changes in the properties of transiently populated species in vitro correlate with the observation of increased aggregation in vivo. The refolding reactions of the WT and K97I proteins are both described by three exponential processes. Two exponential processes characterize fast events (0.1-1.0 s) in folding while the third exponential process correlates with a slow (70 s) single pathway to and from the native state. The K97I replacement affects the earlier steps in the refolding pathway. Aggregation, absent in the WT refolding reaction, occurs in K97I above a critical protein concentration of 18 microM. This observation is consistent with an initial nucleation step mediating protein aggregation. Stopped-flow kinetic studies of the K97I aggregation process demonstrate that K97I aggregates most rapidly during the earliest refolding times, when unfolded protein conformers remain highly populated and the concentration of folding intermediates is low. Folding and aggregation studies together support a model in which the formation of stable folding intermediates afford protection against further K97I aggregation.

Conformation-invariant structures of the alpha1beta1 human hemoglobin dimer.

Analysis of the conformational differences between the oxy and deoxy forms of hemoglobin is complicated by shifting coordinate systems and correlated motions between different parts of the molecule. Methods independent of any frame of reference were used to study the differences in structure between the oxy and deoxy forms of the human hemoglobin alphabeta dimer. Differences between the deoxy and oxy dimer structures can be characterized as rearrangements of 15 substructures persisting between the two conformations. Such substructures are of two kinds, either rigid domains or tertiary substructures. Rigid domains are groups of residues for which all inter-residue distances are conformationally invariant. Residues belonging to a rigid domain do not have to be spatially contiguous nor must they have consecutive sequence numbers. The largest such substructure is a rigid core that spans both the alpha and beta monomers and includes 44% of the dimer. Other rigid domains exist within the heme pockets. An alternative but closely related view of the molecule is based on tertiary substructures. Unlike a rigid domain, a tertiary substructure must have consecutively numbered residues and the residue that ends one tertiary substructure begins the next. The decomposition of the dimer into tertiary substructures represents the dimer as a framework of connected stiff structural elements. Viewed as a set of tertiary substructures, the hemoglobin dimer has the same three principal functional elements: the dimer core and the alpha and beta heme pockets, with the heme pockets held to the dimer core by CD and FG corners. The tertiary substructures that comprise the dimer core include 51% of the molecule. When ligands bind at the hemes, the FG corners communicate structural changes in the hemes to the dimer cores, which may mediate heme-heme cooperativity.

A gel as an array of channels.

We consider the theory of charged point molecules ('probes') being pulled by an electric field through a two-dimensional net of channels that represents a piece of gel. Associated with the position in the net is a free energy of interaction between the probe and the net; this free energy fluctuates randomly with the position of the probe in the net. The free energy is intended to represent weak interactions between the probe and the gel, such as entropy associated with the restriction of the freedom of motion of the probe by the gel, or electrostatic interactions between the probe and charges fixed to the gel. The free energy can be thought of as a surface with the appearance of a rough, hilly landscape spread over the net; the roughness is measured by the standard deviation of the free-energy distribution. Two variations of the model are examined: (1) the net is assumed to have all channels open, or (2) only channels parallel to the electric field are open and all the cross-connecting channels are closed. Model (1) is more realistic but presents a two-dimensional mathematical problem which can only be solved by slow iteration methods, while model (2) is less realistic but presents a one-dimensional problem that can be reduced to simple quadratures and is easy to solve by numerical integration. In both models the mobility of the probe decreases as the roughness parameter is increased, but the effect is larger in the less realistic model (2) if the same free-energy surface is used in both. The mobility in model (2) is reduced both by high points in the rough surface ('bumps') and by low points ('traps'), while in model (1) only the traps are effective, since the probes can flow around the bumps through the cross channels. The mobility in model (2) can be made to agree with model (1) simply by cutting off the bumps of the surface. Thus the simple model (2) can be used in place of the more realistic model (1) that is more difficult to compute.

Rigid domains in proteins: an algorithmic approach to their identification.

A rigid domain, defined here as a tertiary structure common to two or more different protein conformations, can be identified numerically from atomic coordinates by finding sets of residues, one in each conformation, such that the distance between any two residues within the set belonging to one conformation is the same as the distance between the two structurally equivalent residues within the set belonging to any other conformation. The distance between two residues is taken to be the distance between their respective alpha carbon atoms. With the methods of this paper we have found in the deoxy and oxy conformations of the human hemoglobin alpha 1 beta 1 dimer a rigid domain closely related to that previously identified by Baldwin and Chothia (J. Mol. Biol. 129: 175-220, 1979). We provide two algorithms, both using the difference-distance matrix, with which to search for rigid domains directly from atomic coordinates. The first finds all rigid domains in a protein but has storage and processing demands that become prohibitively large with increasing protein size. The second, although not necessarily finding every rigid domain, is computationally tractable for proteins of any size. Because of its efficiency we are able to search protein conformations recursively for groups of non-intersecting domains. Different protein conformations, when aligned by superimposing their respective domain structures, can be examined for structural differences in regions complementing a rigid domain.

The degradation of T7 DNA in converging flow.

The flow-induced degradation of T7 DNA (Molecular Size = 40 Kbp) was studied in a flow device that generates converging flow rather than simple shear flow. We discovered that the sizes of the degradation products were very broadly distributed, covering the range from 10 Kbp to 36 Kbp. An explanation for the broadness of the distribution is given based on a computer simulation of the experiment. The significance of converging flow to the routine handling of large DNA is emphasized.

Understanding the anomalous electrophoresis of bent DNA molecules: a reptation model.

In polyacrylamide gel electrophoresis, the retardation of DNA molecules containing regions of intrinsic curvature can be explained by a novel reptation model that includes the elastic free energy of the DNA chain. Computer simulations based on this model give results that reproduce the dependence of anomalous mobility on gel concentration, which is quantified by new experimental data on the mobilities of circularly permuted isomers of kinetoplast DNA fragments. Fitting of the data required allowing for the elasticity of the gel.

Separations of open-circular DNA using pulsed-field electrophoresis.

The effect of high electric fields on the gel-electrophoretic mobility of open-circular DNA in agarose differs dramatically from that on linear molecules of the same molecular weight. At high fields, sufficiently large circular forms are prevented from migrating into the gel whereas linear molecules and smaller circular DNAs migrate normally. This effect is strongly field dependent, affecting circular molecules of decreasing size with increasing field strength. We have studied this effect with a series of plasmid DNAs ranging from 2.9 to 56 kilobase pairs using continuous and reversing-pulse electric fields. Application of reversing pulses abolishes the effect under certain conditions and supports the model for the gel electrophoresis of open-circular DNA where circular forms are trapped by engaging the free end of an agarose gel fiber.

Counter-ion condensation and system dimensionality.

Condensation of the counter-ions around a highly charged infinitely long cylindrical molecule, such as DNA, can be described in detail in terms of the solutions of the Poisson-Boltzmann (Gouy-Chapman) equation. By using the Alfrey-Berg-Morawetz (1951) solution of this equation one can show that a certain fraction of the counter-ions remain within finite distances of the poly-ion even when the volume of the system is expanded indefinitely; these ions can be appropriately called "condensed". The fraction of the macromolecule's charge represented by these ions is just 1-1/xi, where xi is the linear charge-density parameter of the macromolecule; this is also the value given by Manning's theory. The question arises: Is this property unique to the infinite cylinder? Using the same PB equation, we can consider the infinite charged plane and a large finite charged sphere for comparison. In the case of the plane all of the counter-ions are condensed in the above sense, not just a fraction, for any surface charge density of the plane. These ions form the classical Gouy double layer. On the other hand, none of the counter-ions of the charged sphere are condensed in the above sense, no matter how high the surface charge density. Thus the cylinder is a unique intermediate case in which a fraction of the counter-ions are condensed if the linear charge density is higher than the critical value of unity.