Slip complexity in earthquake fault models.

We summarize studies of earthquake fault models that give rise to slip complexities like those in natural earthquakes. For models of smooth faults between elastically deformable continua, it is critical that the friction laws involve a characteristic distance for slip weakening or evolution of surface state. That results in a finite nucleation size, or coherent slip patch size, h*. Models of smooth faults, using numerical cell size properly small compared to h*, show periodic response or complex and apparently chaotic histories of large events but have not been found to show small event complexity like the self-similar (power law) Gutenberg-Richter frequency-size statistics. This conclusion is supported in the present paper by fully inertial elastodynamic modeling of earthquake sequences. In contrast, some models of locally heterogeneous faults with quasi-independent fault segments, represented approximately by simulations with cell size larger than h* so that the model becomes "inherently discrete," do show small event complexity of the Gutenberg-Richter type. Models based on classical friction laws without a weakening length scale or for which the numerical procedure imposes an abrupt strength drop at the onset of slip have h* = 0 and hence always fall into the inherently discrete class. We suggest that the small-event complexity that some such models show will not survive regularization of the constitutive description, by inclusion of an appropriate length scale leading to a finite h*, and a corresponding reduction of numerical grid size.

Exploring the folding free energy surface of a three-helix bundle protein

The multidimensional free energy surface for a small fast folding helical protein is explored based on first-principle calculations. The model represents the 46-residue segment from fragment B of staphylococcal protein A. The relationship between collapse and tertiary structure formation, and the order of collapse and secondary structure formation, are investigated. We find that the initial collapse process gives rise to a transition state with about 30% of the native tertiary structure and 50–70% of the native helix content. We also observe two distinct distributions of native helix in this collapsed state (Rg ≈ 12 Å), one with about 20% of the native helical hydrogen bonds, the other with near 70%. The former corresponds to a local minimum. The barrier from this metastable state to the native state is about 2 kBT. In the latter case, folding is essentially a downhill process involving topological assembly. In addition, the order of formation of secondary structure among the three helices is examined. We observe cooperative formation of the secondary structure in helix I and helix II. Secondary structure in helix III starts to form following the formation of certain secondary structure in both helix I and helix II. Comparisons of our results with those from theory and experiment are made.

Partial molar volume, surface area, and hydration changes for equilibrium unfolding and formation of aggregation transition state: High-pressure and cosolute studies on recombinant human IFN-γ

The equilibrium dissociation of recombinant human IFN-γ was monitored as a function of pressure and sucrose concentration. The partial molar volume change for dissociation was −209 ± 13 ml/mol of dimer. The specific molar surface area change for dissociation was 12.7 ± 1.6 nm2/molecule of dimer. The first-order aggregation rate of recombinant human IFN-γ in 0.45 M guanidine hydrochloride was studied as a function of sucrose concentration and pressure. Aggregation proceeded through a transition-state species, N*. Sucrose reduced aggregation rate by shifting the equilibrium between native state (N) and N* toward the more compact N. Pressure increased aggregation rate through increased solvation of the protein, which exposes more surface area, thus shifting the equilibrium away from N toward N*. The changes in partial molar volume and specific molar surface area between the N* and N were −41 ± 9 ml/mol of dimer and 3.5 ± 0.2 nm2/molecule, respectively. Thus, the structural change required for the formation of the transition state for aggregation is small relative to the difference between N and the dissociated state. Changes in waters of hydration were estimated from both specific molar surface area and partial molar volume data. From partial molar volume data, estimates were 25 and 128 mol H2O/mol dimer for formation of the aggregation transition state and for dissociation, respectively. From surface area data, estimates were 27 and 98 mol H2O/mol dimer. Osmotic stress theory yielded values ≈4-fold larger for both transitions.

DNA-mediated immunization in a transgenic mouse model of the hepatitis B surface antigen chronic carrier state.

Transgenic mice expressing the sequences coding for the envelope proteins of the hepatitis B virus (HBV) in the liver have been used as a model of the HBV chronic carrier state. We evaluated the possibility of inducing a specific immune response to the viral envelope antigens and thus potentially controlling chronic HBV infection. Using HBV-specific DNA-mediated immunization in this transgenic model, we show that the immune response induced after a single intramuscular injection of DNA resulted in the complete clearance of circulating hepatitis B surface antigen and in the long-term control of transgene expression in hepatocytes. This response does not involve a detectable cytopathic effect in the liver. Adoptive transfer of fractionated primed spleen cells from DNA-immunized mice shows that T cells are responsible for the down-regulation of HBV mRNA in the liver of transgenic mice. To our knowledge, this is the first demonstration of a potential immunotherapeutic application of DNA-mediated immunization against an infectious disease and raises the possibility of designing more effective ways of treating HBV chronic carriers.

Membrane-bound state of the colicin E1 channel domain as an extended two-dimensional helical array

Atomic level structures have been determined for the soluble forms of several colicins and toxins, but the structural changes that occur after membrane binding have not been well characterized. Changes occurring in the transition from the soluble to membrane-bound state of the C-terminal 190-residue channel polypeptide of colicin E1 (P190) bound to anionic membranes are described. In the membrane-bound state, the α-helical content increases from 60–64% to 80–90%, with a concomitant increase in the average length of the helical segments from 12 to 16 or 17 residues, close to the length required to span the membrane bilayer in the open channel state. The average distance between helical segments is increased and interhelix interactions are weakened, as shown by a major loss of tertiary structure interactions, decreased efficiency of fluorescence resonance energy transfer from an energy donor on helix V of P190 to an acceptor on helix IX, and decreased resonance energy transfer at higher temperatures, not observed in soluble P190, implying freedom of motion of helical segments. Weaker interactions are also shown by a calorimetric thermal transition of low cooperativity, and the extended nature of the helical array is shown by a 3- to 4-fold increase in the average area subtended per molecule to 4,200 Å2 on the membrane surface. The latter, with analysis of the heat capacity changes, implies the absence of a developed hydrophobic core in the membrane-bound P190. The membrane interfacial layer thus serves to promote formation of a highly helical extended two-dimensional flexible net. The properties of the membrane-bound state of the colicin channel domain (i.e., hydrophobic anchor, lengthened and loosely coupled α-helices, and close association with the membrane interfacial layer) are plausible structural features for the state that is a prerequisite for voltage gating, formation of transmembrane helices, and channel opening.

Unique molecular surface features of in vivo tolerized T cells

Differential expression of surface markers can frequently be used to distinguish functional subsets of T cells, yet a surface phenotype unique to T cells induced into an anergic state has not been described. Here, we report that CD4 T cells rendered anergic in vivo by superantigen can be identified by loss of the 6C10 T cell marker. Inoculation of Vβ8.1 T cell antigen receptor (TCR) transgenic mice with a Vβ8.1-reactive minor lymphocyte-stimulating superantigen (Mls-1a) induces tolerance to Mls-1a by clonal anergy. CD4 lymph node T cells from Mls-1a inoculated transgenic mice enriched for the 6C10− phenotype neither proliferate nor produce interleukin-2 upon TCR engagement, whereas 6C10+ CD4 T cells retain responsiveness. Analysis of T cell memory markers demonstrate that 6C10− T cells remain 3G11hi but express heterogeneous levels of CD45RB, CD62L, CD44, and the CD69 early activation marker, suggesting that T cells at various degrees of activation can be functionally anergic. These studies demonstrate that anergic T cells can be purified based on 6C10 expression permitting examination of issues concerning biochemical and biological features specific to T cell anergy.

A peptide that inhibits hydroxyapatite growth is in an extended conformation on the crystal surface

Proteins play an important role in the biological mechanisms controlling hard tissue development, but the details of molecular recognition at inorganic crystal interfaces remain poorly characterized. We have applied a recently developed homonuclear dipolar recoupling solid-state NMR technique, dipolar recoupling with a windowless sequence (DRAWS), to directly probe the conformation of an acidic peptide adsorbed to hydroxyapatite (HAP) crystals. The phosphorylated hexapeptide, DpSpSEEK (N6, where pS denotes phosphorylated serine), was derived from the N terminus of the salivary protein statherin. Constant-composition kinetic characterization demonstrated that, like the native statherin, this peptide inhibits the growth of HAP seed crystals when preadsorbed to the crystal surface. The DRAWS technique was used to measure the internuclear distance between two 13C labels at the carbonyl positions of the adjacent phosphoserine residues. Dipolar dephasing measured at short mixing times yielded a mean separation distance of 3.2 ± 0.1 Å. Data obtained by using longer mixing times suggest a broad distribution of conformations about this average distance. Using a more complex model with discrete α-helical and extended conformations did not yield a better fit to the data and was not consistent with chemical shift analysis. These results suggest that the peptide is predominantly in an extended conformation rather than an α-helical state on the HAP surface. Solid-state NMR approaches can thus be used to determine directly the conformation of biologically relevant peptides on HAP surfaces. A better understanding of peptide and protein conformation on biomineral surfaces may provide design principles useful for the modification of orthopedic and dental implants with coatings and biological growth factors that are designed to enhance biocompatibility with surrounding tissue.

Exploring the origins of topological frustration: Design of a minimally frustrated model of fragment B of protein A

Topological frustration in an energetically unfrustrated off-lattice model of the helical protein fragment B of protein A from Staphylococcus aureus was investigated. This Gō-type model exhibited thermodynamic and kinetic signatures of a well-designed two-state folder with concurrent collapse and folding transitions and single exponential kinetics at the transition temperature. Topological frustration is determined in the absence of energetic frustration by the distribution of Fersht φ values. Topologically unfrustrated systems present a unimodal distribution sharply peaked at intermediate φ, whereas highly frustrated systems display a bimodal distribution peaked at low and high φ values. The distribution of φ values in protein A was determined both thermodynamically and kinetically. Both methods yielded a unimodal distribution centered at φ = 0.3 with tails extending to low and high φ values, indicating the presence of a small amount of topological frustration. The contacts with high φ values were located in the turn regions between helices I and II and II and III, intimating that these hairpins are in large part required in the transition state. Our results are in good agreement with all-atom simulations of protein A, as well as lattice simulations of a three- letter code 27-mer (which can be compared with a 60-residue helical protein). The relatively broad unimodal distribution of φ values obtained from the all-atom simulations and that from the minimalist model for the same native fold suggest that the structure of the transition state ensemble is determined mostly by the protein topology and not energetic frustration.

Preferential exclusion of sucrose from recombinant interleukin-1 receptor antagonist: Role in restricted conformational mobility and compaction of native state

Understanding the mechanism for sucrose-induced protein stabilization is important in many diverse fields, ranging from biochemistry and environmental physiology to pharmaceutical science. Timasheff and Lee [Lee, J. C. & Timasheff, S. N. (1981) J. Biol. Chem. 256, 7193–7201] have established that thermodynamic stabilization of proteins by sucrose is due to preferential exclusion of the sugar from the protein’s surface, which increases protein chemical potential. The current study measures the preferential exclusion of 1 M sucrose from a protein drug, recombinant interleukin 1 receptor antagonist (rhIL-1ra). It is proposed that the degree of preferential exclusion and increase in chemical potential are directly proportional to the protein surface area and that, hence, the system will favor the protein state with the smallest surface area. This mechanism explains the observed sucrose-induced restriction of rhIL-1ra conformational fluctuations, which were studied by hydrogen–deuterium exchange and cysteine reactivity measurements. Furthermore, infrared spectroscopy of rhlL-1ra suggested that a more ordered native conformation is induced by sucrose. Electron paramagnetic resonance spectroscopy demonstrated that in the presence of sucrose, spin-labeled cysteine 116 becomes more buried in the protein’s interior and that the hydrodynamic diameter of the protein is reduced. The preferential exclusion of sucrose from the protein and the resulting shift in the equilibrium between protein states toward the most compact conformation account for sucrose-induced effects on rhIL-1ra.

Equilibrium distributions of topological states in circular DNA: Interplay of supercoiling and knotting

Two variables define the topological state of closed double-stranded DNA: the knot type, K, and ΔLk, the linking number difference from relaxed DNA. The equilibrium distribution of probabilities of these states, P(ΔLk, K), is related to two conditional distributions: P(ΔLk|K), the distribution of ΔLk for a particular K, and P(K|ΔLk) and also to two simple distributions: P(ΔLk), the distribution of ΔLk irrespective of K, and P(K). We explored the relationships between these distributions. P(ΔLk, K), P(ΔLk), and P(K|ΔLk) were calculated from the simulated distributions of P(ΔLk|K) and of P(K). The calculated distributions agreed with previous experimental and theoretical results and greatly advanced on them. Our major focus was on P(K|ΔLk), the distribution of knot types for a particular value of ΔLk, which had not been evaluated previously. We found that unknotted circular DNA is not the most probable state beyond small values of ΔLk. Highly chiral knotted DNA has a lower free energy because it has less torsional deformation. Surprisingly, even at |ΔLk| > 12, only one or two knot types dominate the P(K|ΔLk) distribution despite the huge number of knots of comparable complexity. A large fraction of the knots found belong to the small family of torus knots. The relationship between supercoiling and knotting in vivo is discussed.

A transient expansion of the native state precedes aggregation of recombinant human interferon-γ

Aggregation of proteins, even under conditions favoring the native state, is a ubiquitous problem in biotechnology and biomedical engineering. Providing a mechanistic basis for the pathways that lead to aggregation should allow development of rational approaches for its prevention. We have chosen recombinant human interferon-γ (rhIFN-γ) as a model protein for a mechanistic study of aggregation. In the presence of 0.9 M guanidinium hydrochloride, rhIFN-γ aggregates with first order kinetics, a process that is inhibited by addition of sucrose. We describe a pathway that accounts for both the observed first-order aggregation of rhIFN-γ and the effect of sucrose. In this pathway, aggregation proceeds through a transient expansion of the native state. Sucrose shifts the equilibrium within the ensemble of rhIFN-γ native conformations to favor the most compact native species over more expanded ones, thus stabilizing rhIFN-γ against aggregation. This phenomenon is attributed to the preferential exclusion of sucrose from the protein surface. In addition, kinetic analysis combined with solution thermodynamics shows that only a small (9%) expansion surface area is needed to form the transient native state that precedes aggregation. The approaches used here link thermodynamics and aggregation kinetics to provide a powerful tool for understanding both the pathway of protein aggregation and the rational use of excipients to inhibit the process.

The active oligomeric state of the minimalistic influenza virus M2 ion channel is a tetramer

The influenza A virus M2 integral membrane protein is an ion channel that permits protons to enter virus particles during uncoating of virions in endosomes and also modulates the pH of the trans-Golgi network in virus-infected cells. The M2 protein is a homo-oligomer of 97 residues, and analysis by chemical cross-linking and SDS/PAGE indicates M2 forms a tetramer. However, a higher order molecular form is sometimes observed and, thus, it is necessary to determine the active form of the molecule. This was done by studying the currents of oocytes that expressed mixtures of the wild-type M2 protein (epitope tagged) and the mutant protein M2-V27S, which is resistant to the inhibitor amantadine. The composition of mixed oligomers of the two proteins expressed at the plasma membrane of individual oocytes was quantified after antibody capture of the cell surface expressed molecules and it was found that the subunits mixed freely. When the ratio of wild-type to mutant protein subunits was 0.85:0.15, the amantadine sensitivity was reduced to 50% and for a ratio of 0.71:0.29 to 20%. These results are consistent with the amantadine-resistant mutant being dominant and the oligomeric state being a tetramer.

Cell surface expression of mouse macrosialin and human CD68 and their role as macrophage receptors for oxidized low density lipoprotein

We have previously identified a 94- to 97-kDa oxidized low density lipoprotein (LDL)-binding protein in mouse macrophages as macrosialin (MS), a member of the lamp family. Earlier immunostaining studies have shown that MS and its human homolog, CD68, are predominantly intracellular proteins. However, using sensitive techniques such as flow cytometry (FACS) and cell-surface-specific biotinylation, we now show that there is significant surface expression of these proteins. FACS analysis of intact cells using mAb FA/11 showed small but definite surface expression of MS in resident mouse peritoneal macrophages but this was greatly enhanced with thioglycollate elicitation. Biotinylation of intact cells and detergent-solubilized cell preparations followed by immunoprecipitation revealed 10–15% of the total MS content of elicited macrophages on the plasma membrane. Similar results were obtained with untreated RAW 264.7 cells. FACS analysis of intact THP-1 monocytic cells showed minimal surface expression of CD68 on unactivated cells (4% of total cell content). Stimulation with phorbol 12-myristate 13-acetate increased both surface and total CD68 expression considerably. Furthermore, the specific binding at 4°C and uptake at 37°C of 125I-labeled oxidized LDL by activated THP-1 cells was inhibited by 30–50% by CD68 mAbs KP-1 and EBM-11. Thus, although the surface expression of MS/CD68 at steady-state represents only a small percentage of their total cellular content, these proteins can play a significant role in oxidized LDL uptake by activated macrophages in vitro and could contribute to foam cell formation in atherosclerotic lesions.

The photoisomerization of retinal in bacteriorhodopsin: Experimental evidence for a three-state model

The primary events in the all-trans to 13-cis photoisomerization of retinal in bacteriorhodopsin have been investigated with femtosecond time-resolved absorbance spectroscopy. Spectra measured over a broad range extending from 7000 to 22,400 cm−1 reveal features whose dynamics are inconsistent with a model proposed earlier to account for the highly efficient photoisomerization process. Emerging from this work is a new three-state model. Photoexcitation of retinal with visible light accesses a shallow well on the excited state potential energy surface. This well is bounded by a small barrier, arising from an avoided crossing that separates the Franck–Condon region from the nearby reactive region of the photoisomerization coordinate. At ambient temperatures, the reactive region is accessed with a time constant of ≈500 fs, whereupon the retinal rapidly twists and encounters a second avoided crossing region. The protein mediates the passage into the second avoided crossing region and thereby exerts control over the quantum yield for forming 13-cis retinal. The driving force for photoisomerization resides in the retinal, not in the surrounding protein. This view contrasts with an earlier model where photoexcitation was thought to access directly a reactive region of the excited-state potential and thereby drive the retinal to a twisted conformation within 100–200 fs.

Effect of the protein import machinery at the mitochondrial surface on precursor stability

Many biological processes require proteins to undergo conformational changes at the surface of membranes. For example, some precursor proteins unfold at the surface of mitochondria and chloroplasts before translocation into the organelles, and toxins such as colicin A unfold to the molten globule state at bacterial surfaces before inserting into the cell membrane. It is commonly thought that the membrane surfaces and the associated protein machinery destabilize the substrate proteins and that this effect is required for membrane insertion or translocation. One of the best characterized translocation processes is protein import into mitochondria. By measuring the contributions of individual interactions within a model protein to its stability at the mitochondrial surface and in free solution, we show here that the mitochondrial surface neither induces the molten globule state in this protein nor preferentially destabilizes any type of interaction (e.g., hydrogen bonds, nonpolar, etc.) within the protein. Because it is not possible to measure absolute protein stability at the surface of mitochondria, we determined the stability of a tightly associated protein–protein complex at the mitochondrial import site as a model of the stability of a protein. We found the binding constants of the protein–protein complex at the mitochondrial surface and in free solution to be identical. Our results demonstrate that the mitochondrial surface does not destabilize importing precursor proteins in its vicinity.