A statistical mechanical model for β-hairpin kinetics

Understanding the mechanism of protein secondary structure formation is an essential part of the protein-folding puzzle. Here, we describe a simple statistical mechanical model for the formation of a β-hairpin, the minimal structural element of the antiparallel β-pleated sheet. The model accurately describes the thermodynamic and kinetic behavior of a 16-residue, β-hairpin-forming peptide, successfully explaining its two-state behavior and apparent negative activation energy for folding. The model classifies structures according to their backbone conformation, defined by 15 pairs of dihedral angles, and is further simplified by considering only the 120 structures with contiguous stretches of native pairs of backbone dihedral angles. This single sequence approximation is tested by comparison with a more complete model that includes the 215 possible conformations and 15 × 215 possible kinetic transitions. Finally, we use the model to predict the equilibrium unfolding curves and kinetics for several variants of the β-hairpin peptide.

Stretching single-domain proteins: Phase diagram and kinetics of force-induced unfolding

Single-molecule force spectroscopy reveals unfolding of domains in titin on stretching. We provide a theoretical framework for these experiments by computing the phase diagrams for force-induced unfolding of single-domain proteins using lattice models. The results show that two-state folders (at zero force) unravel cooperatively, whereas stretching of non-two-state folders occurs through intermediates. The stretching rates of individual molecules show great variations reflecting the heterogeneity of force-induced unfolding pathways. The approach to the stretched state occurs in a stepwise “quantized” manner. Unfolding dynamics and forces required to stretch proteins depend sensitively on topology. The unfolding rates increase exponentially with force f till an optimum value, which is determined by the barrier to unfolding when f = 0. A mapping of these results to proteins shows qualitative agreement with force-induced unfolding of Ig-like domains in titin. We show that single-molecule force spectroscopy can be used to map the folding free energy landscape of proteins in the absence of denaturants.

Construction of a high-affinity receptor site for dihydropyridine agonists and antagonists by single amino acid substitutions in a non-l-type Ca2+ channel

The activity of l-type Ca2+ channels is increased by dihydropyridine (DHP) agonists and inhibited by DHP antagonists, which are widely used in the therapy of cardiovascular disease. These drugs bind to the pore-forming α1 subunits of l-type Ca2+ channels. To define the minimal requirements for DHP binding and action, we constructed a high-affinity DHP receptor site by substituting a total of nine amino acid residues from DHP-sensitive l-type α1 subunits into the S5 and S6 transmembrane segments of domain III and the S6 transmembrane segment of domain IV of the DHP-insensitive P/Q-type α1A subunit. The resulting chimeric α1A/DHPS subunit bound DHP antagonists with high affinity in radioligand binding assays and was inhibited by DHP antagonists with high affinity in voltage clamp experiments. Substitution of these nine amino acid residues yielded 86% of the binding energy of the l-type α1C subunit and 92% of the binding energy of the l-type α1S subunit for the high-affinity DHP antagonist PN200–110. The activity of chimeric Ca2+ channels containing α1A/DHPS was increased 3.5 ± 0.7-fold by the DHP agonist (−)Bay K8644. The effect of this agonist was stereoselective as in l-type Ca2+ channels since (+) Bay K8644 inhibited the activity of α1A/DHPS. The results show conclusively that DHP agonists and antagonists bind to a single receptor site at which they have opposite effects on Ca2+ channel activity. This site contains essential components from both domains III and IV, consistent with a domain interface model for binding and allosteric modulation of Ca2+ channel activity by DHPs.

Single-molecule protein folding: Diffusion fluorescence resonance energy transfer studies of the denaturation of chymotrypsin inhibitor 2

We report single-molecule folding studies of a small, single-domain protein, chymotrypsin inhibitor 2 (CI2). CI2 is an excellent model system for protein folding studies and has been extensively studied, both experimentally (at the ensemble level) and theoretically. Conformationally assisted ligation methodology was used to synthesize the proteins and site-specifically label them with donor and acceptor dyes. Folded and denatured subpopulations were observed by fluorescence resonance energy transfer (FRET) measurements on freely diffusing single protein molecules. Properties of these subpopulations were directly monitored as a function of guanidinium chloride concentration. It is shown that new information about different aspects of the protein folding reaction can be extracted from such subpopulation properties. Shifts in the mean transfer efficiencies are discussed, FRET efficiency distributions are translated into potentials, and denaturation curves are directly plotted from the areas of the FRET peaks. Changes in stability caused by mutation also are measured by comparing pseudo wild-type CI2 with a destabilized mutant (K17G). Current limitations and future possibilities and prospects for single-pair FRET protein folding investigations are discussed.

Structure and function of the Bacillus hybrid enzyme GluXyn-1: Native-like jellyroll fold preserved after insertion of autonomous globular domain

The 1,3–1,4-β-glucanase from Bacillus macerans (wtGLU) and the 1,4-β-xylanase from Bacillus subtilis (wtXYN) are both single-domain jellyroll proteins catalyzing similar enzymatic reactions. In the fusion protein GluXyn-1, the two proteins are joined by insertion of the entire XYN domain into a surface loop of cpMAC-57, a circularly permuted variant of wtGLU. GluXyn-1 was generated by protein engineering methods, produced in Escherichia coli and shown to fold spontaneously and have both enzymatic activities at wild-type level. The crystal structure of GluXyn-1 was determined at 2.1 Å resolution and refined to R = 17.7% and R(free) = 22.4%. It shows nearly ideal, native-like folding of both protein domains and a small, but significant hinge bending between the domains. The active sites are independent and accessible explaining the observed enzymatic activity. Because in GluXyn-1 the complete XYN domain is inserted into the compact folding unit of GLU, the wild-type-like activity and tertiary structure of the latter proves that the folding process of GLU does not depend on intramolecular interactions that are short-ranged in the sequence. Insertion fusions of the GluXyn-1 type may prove to be an easy route toward more stable bifunctional proteins in which the two parts are more closely associated than in linear end-to-end protein fusions.

Limited internal friction in the rate-limiting step of a two-state protein folding reaction

Small, single-domain proteins typically fold via a compact transition-state ensemble in a process well fitted by a simple, two-state model. To characterize the rate-limiting conformational changes that underlie two-state folding, we have investigated experimentally the effects of changing solvent viscosity on the refolding of the IgG binding domain of protein L. In conjunction with numerical simulations, our results indicate that the rate-limiting conformational changes of the folding of this domain are strongly coupled to solvent viscosity and lack any significant “internal friction” arising from intrachain collisions. When compared with the previously determined solvent viscosity dependencies of other, more restricted conformational changes, our results suggest that the rate-limiting folding transition involves conformational fluctuations that displace considerable amounts of solvent. Reconciling evidence that the folding transition state ensemble is comprised of highly collapsed species with these and similar, previously reported results should provide a significant constraint for theoretical models of the folding process.

Structural determinants of oxidative folding in proteins

A method for determining the kinetic fate of structured disulfide species (i.e., whether they are preferentially oxidized or reshuffle back to an unstructured disulfide species) is introduced. The method relies on the sensitivity of unstructured disulfide species to low concentrations of reducing agents. Because a structured des species that preferentially reshuffles generally first rearranges to an unstructured species, a small concentration of reduced DTT (e.g., 260 μM) suffices to distinguish on-pathway intermediates from dead-end species. We apply this method to the oxidative folding of bovine pancreatic ribonuclease A (RNase A) and show that des[40–95] and des[65–72] are productive intermediates, whereas des[26–84] and des[58–110] are metastable dead-end species that preferentially reshuffle. The key factor in determining the kinetic fate of these des species is the relative accessibility of both their thiol groups and disulfide bonds. Productive intermediates tend to be disulfide-secure, meaning that their structural fluctuations preferentially expose their thiol groups, while keeping their disulfide bonds buried. By contrast, dead-end species tend to be disulfide-insecure, in that their structural fluctuations expose their disulfide bonds in concert with their thiol groups. This distinction leads to four generic types of oxidative folding pathways. We combine these results with those of earlier studies to suggest a general three-stage model of oxidative folding of RNase A and other single-domain proteins with multiple disulfide bonds.

Cladistic association analysis of Y chromosome effects on alcohol dependence and related personality traits

Association between Y chromosome haplotype variation and alcohol dependence and related personality traits was investigated in a large sample of psychiatrically diagnosed Finnish males. Haplotypes were constructed for 359 individuals using alleles at eight loci (seven microsatellite loci and a nucleotide substitution in the DYZ3 alphoid satellite locus). A cladogram linking the 102 observed haplotype configurations was constructed by using parsimony with a single-step mutation model. Then, a series of contingency tables nested according to the cladogram hierarchy were used to test for association between Y haplotype and alcohol dependence. Finally, using only alcohol-dependent subjects, we tested for association between Y haplotype and personality variables postulated to define subtypes of alcoholism—antisocial personality disorder, novelty seeking, harm avoidance, and reward dependence. Significant association with alcohol dependence was observed at three Y haplotype clades, with significance levels of P = 0.002, P = 0.020, and P = 0.010. Within alcohol-dependent subjects, no relationship was revealed between Y haplotype and antisocial personality disorder, novelty seeking, harm avoidance, or reward dependence. These results demonstrate, by using a fully objective association design, that differences among Y chromosomes contribute to variation in vulnerability to alcohol dependence. However, they do not demonstrate an association between Y haplotype and the personality variables thought to underlie the subtypes of alcoholism.

Identification of a cytochrome b-type NAD(P)H oxidoreductase ubiquitously expressed in human cells

Cytochrome b-type NAD(P)H oxidoreductases are involved in many physiological processes, including iron uptake in yeast, the respiratory burst, and perhaps oxygen sensing in mammals. We have identified a cytosolic cytochrome b-type NAD(P)H oxidoreductase in mammals, a flavohemoprotein (b5+b5R) containing cytochrome b5 (b5) and b5 reductase (b5R) domains. A genetic approach, using blast searches against dbest for FAD-, NAD(P)H-binding sequences followed by reverse transcription–PCR, was used to clone the complete cDNA sequence of human b5+b5R from the hepatoma cell line Hep 3B. Compared with the classical single-domain b5 and b5R proteins localized on endoplasmic reticulum membrane, b5+b5R also has binding motifs for heme, FAD, and NAD(P)H prosthetic groups but no membrane anchor. The human b5+b5R transcript was expressed at similar levels in all tissues and cell lines that were tested. The two functional domains b5* and b5R* are linked by an approximately 100-aa-long hinge bearing no sequence homology to any known proteins. When human b5+b5R was expressed as c-myc adduct in COS-7 cells, confocal microscopy revealed a cytosolic localization at the perinuclear space. The recombinant b5+b5R protein can be reduced by NAD(P)H, generating spectrum typical of reduced cytochrome b with alpha, beta, and Soret peaks at 557, 527, and 425 nm, respectively. Human b5+b5R flavohemoprotein is a NAD(P)H oxidoreductase, demonstrated by superoxide production in the presence of air and excess NAD(P)H and by cytochrome c reduction in vitro. The properties of this protein make it a plausible candidate oxygen sensor.

Kinetic equilibrium of forces and molecular events in muscle contraction

In biomolecular systems, the mechanical transfer of free energy occurs with both high efficiency and high speed. It is shown here that such a transfer can be achieved only if the participating free-energy-storing elements exhibit opposing relationships between their content of free energy and the force they exert in the transfer direction. A kinetic equilibrium of forces (KEF) results, in which the transfer of free energy is mediated essentially by thermal molecular motion. On the basis of present evidence, KEF is used as a guiding principle in developing a mechanical model of the crossbridge cycle in muscle contraction. The model allows the basic features of molecular events to be visualized in terms of plausible structures. Real understanding of the process will require identification of the elements that perform the functions described here. Besides chemomechanical energy transduction, KEF may have a role in other biomolecular processes in which free energy is transferred mechanically over large distances.

Truncated hexa-octahedral magnetite crystals in ALH84001: Presumptive biosignatures

McKay et al. [(1996) Science 273, 924–930] suggested that carbonate globules in the meteorite ALH84001 contained the fossil remains of Martian microbes. We have characterized a subpopulation of magnetite (Fe3O4) crystals present in abundance within the Fe-rich rims of these carbonate globules. We find these Martian magnetites to be both chemically and physically identical to terrestrial, biogenically precipitated, intracellular magnetites produced by magnetotactic bacteria strain MV-1. Specifically, both magnetite populations are single-domain and chemically pure, and exhibit a unique crystal habit we describe as truncated hexa-octahedral. There are no known reports of inorganic processes to explain the observation of truncated hexa-octahedral magnetites in a terrestrial sample. In bacteria strain MV-1 their presence is therefore likely a product of Natural Selection. Unless there is an unknown and unexplained inorganic process on Mars that is conspicuously absent on the Earth and forms truncated hexa-octahedral magnetites, we suggest that these magnetite crystals in the Martian meteorite ALH84001 were likely produced by a biogenic process. As such, these crystals are interpreted as Martian magnetofossils and constitute evidence of the oldest life yet found.

Computation of Surface Electrical Potentials of Plant Cell Membranes : Correspondence to Published Zeta Potentials from Diverse Plant Sources

A Gouy-Chapman-Stern model has been developed for the computation of surface electrical potential (ψ0) of plant cell membranes in response to ionic solutes. The present model is a modification of an earlier version developed to compute the sorption of ions by wheat (Triticum aestivum L. cv Scout 66) root plasma membranes. A single set of model parameters generates values for ψ0 that correlate highly with published ζ potentials of protoplasts and plasma membrane vesicles from diverse plant sources. The model assumes ion binding to a negatively charged site (R− = 0.3074 μmol m−2) and to a neutral site (P0 = 2.4 μmol m−2) according to the reactions R− + IΖ ⇌ RIΖ−1 and P0 + IΖ ⇌ PIΖ, where IΖ represents an ion of charge Ζ. Binding constants for the negative site are 21,500 m−1 for H+, 20,000 m−1 for Al3+, 2,200 m−1 for La3+, 30 m−1 for Ca2+ and Mg2+, and 1 m−1 for Na+ and K+. Binding constants for the neutral site are 1/180 the value for binding to the negative site. Ion activities at the membrane surface, computed on the basis of ψ0, appear to determine many aspects of plant-mineral interactions, including mineral nutrition and the induction and alleviation of mineral toxicities, according to previous and ongoing studies. A computer program with instructions for the computation of ψ0, ion binding, ion concentrations, and ion activities at membrane surfaces may be requested from the authors.

Compressibility as a means to detect and characterize globular protein states.

We report compressibility data on single-domain, globular proteins which suggest a general relationship between protein conformational transitions and delta kzeroS, the change in the partial specific adiabatic compressibility which accompanies the transition. Specifically, we find transitions between native and compact intermediate states to be accompanied by small increases in kzeroS of +(1-4) x 10(-6) cm3.g-1.bar-1 (1 bar = 100 kPa). By contrast, transitions between native and partially unfolded states are accompanied by small decreases in kzeroS of -(3-7) x 10(-6) cm3.g-1.bar-1, while native-to-fully unfolded transitions result in large decreases in kzeroS of -(18-20) x 10(-6) cm3.g-1.bar-1. Thus, for the single-domain, globular proteins studied here, changes in kzeroS correlate with the type of transition being monitored, independent of the specific protein. Consequently, kzeroS measurements may provide a convenient approach for detecting the existence of and for defining the nature of protein transitions, while also characterizing the hydration properties of individual protein states.

Crystal structure of the cell cycle-regulatory protein suc1 reveals a beta-hinge conformational switch.

The Schizosaccharomyces pombe cell cycle-regulatory protein suc1, named as the suppressor of cdc2 temperature-sensitive mutations, is essential for cell cycle progression. To understand suc1 structure-function relationships and to help resolve conflicting interpretations of suc1 function based on genetic studies of suc1 and its functional homologs in both lower and higher eukaryotes, we have determined the crystal structure of the beta-interchanged suc1 dimer. Each domain consists of three alpha-helices and a four-stranded beta-sheet, completed by the interchange of terminal beta-strands between the two subunits. This beta-interchanged suc1 dimer, when compared with the beta-hairpin single-domain folds of suc1, reveals a beta-hinge motif formed by the conserved amino acid sequence HVPEPH. This beta-hinge mediates the subunit conformation and assembly of suc1: closing produces the intrasubunit beta-hairpin and single-domain fold, whereas opening leads to the intersubunit beta-strand interchange and interlocked dimer assembly reported here. This conformational switch markedly changes the surface accessibility of sequence-conserved residues available for recognition of cyclin-dependent kinase, suggesting a structural mechanism for beta-hinge-mediated regulation of suc1 biological function. Thus, suc1 belongs to the family of domain-swapping proteins, consisting of intertwined and dimeric protein structures in which the dual assembly modes regulate their function.

Mechanical and chemical unfolding of a single protein: A comparison

Is the mechanical unraveling of protein domains by atomic force microscopy (AFM) just a technological feat or a true measurement of their unfolding? By engineering a protein made of tandem repeats of identical Ig modules, we were able to get explicit AFM data on the unfolding rate of a single protein domain that can be accurately extrapolated to zero force. We compare this with chemical unfolding rates for untethered modules extrapolated to 0 M denaturant. The unfolding rates obtained by the two methods are the same. Furthermore, the transition state for unfolding appears at the same position on the folding pathway when assessed by either method. These results indicate that mechanical unfolding of a single protein by AFM does indeed reflect the same event that is observed in traditional unfolding experiments. The way is now open for the extensive use of AFM to measure folding reactions at the single-molecule level. Single-molecule AFM recordings have the added advantage that they define the reaction coordinate and expose rare unfolding events that cannot be observed in the absence of chemical denaturants.