Transition-state structure as a unifying basis in protein-folding mechanisms: Contact order, chain topology, stability, and the extended nucleus mechanism

I attempt to reconcile apparently conflicting factors and mechanisms that have been proposed to determine the rate constant for two-state folding of small proteins, on the basis of general features of the structures of transition states. Φ-Value analysis implies a transition state for folding that resembles an expanded and distorted native structure, which is built around an extended nucleus. The nucleus is composed predominantly of elements of partly or well-formed native secondary structure that are stabilized by local and long-range tertiary interactions. These long-range interactions give rise to connecting loops, frequently containing the native loops that are poorly structured. I derive an equation that relates differences in the contact order of a protein to changes in the length of linking loops, which, in turn, is directly related to the unfavorable free energy of the loops in the transition state. Kinetic data on loop extension mutants of CI2 and α-spectrin SH3 domain fit the equation qualitatively. The rate of folding depends primarily on the interactions that directly stabilize the nucleus, especially those in native-like secondary structure and those resulting from the entropy loss from the connecting loops, which vary with contact order. This partitioning of energy accounts for the success of some algorithms that predict folding rates, because they use these principles either explicitly or implicitly. The extended nucleus model thus unifies the observations of rate depending on both stability and topology.

How does a β-hairpin fold/unfold? Competition between topology and heterogeneity in a solvable model

We study the competition between topological effects and sequence inhomogeneities in determining the thermodynamics and the un/folding kinetics of a β-hairpin. Our work utilizes a new exactly solvable model that allows for arbitrary configurations of native contacts. In general, the competition between heterogeneity and topology results in a crossover of the dominant transition state. Interestingly, near this crossover, the single reaction coordinate picture can be seriously misleading. Our results also suggest that inferring the folding pathway from unfolding simulations is not always justified.

The Integral Membrane Protein Snl1p Is Genetically Linked to Yeast Nuclear Pore Complex Function

Integral membrane proteins are predicted to play key roles in the biogenesis and function of nuclear pore complexes (NPCs). Revealing how the transport apparatus is assembled will be critical for understanding the mechanism of nucleocytoplasmic transport. We observed that expression of the carboxyl-terminal 200 amino acids of the nucleoporin Nup116p had no effect on wild-type yeast cells, but it rendered the nup116 null strain inviable at all temperatures and coincidentally resulted in the formation of nuclear membrane herniations at 23°C. To identify factors related to NPC function, a genetic screen for high-copy suppressors of this lethal nup116-C phenotype was conducted. One gene (designated SNL1 for suppressor of nup116-C lethal) was identified whose expression was necessary and sufficient for rescuing growth. Snl1p has a predicted molecular mass of 18.3 kDa, a putative transmembrane domain, and limited sequence similarity to Pom152p, the only previously identified yeast NPC-associated integral membrane protein. By both indirect immunofluorescence microscopy and subcellular fractionation studies, Snl1p was localized to both the nuclear envelope and the endoplasmic reticulum. Membrane extraction and topology assays suggested that Snl1p was an integral membrane protein, with its carboxyl-terminal region exposed to the cytosol. With regard to genetic specificity, the nup116-C lethality was also suppressed by high-copy GLE2 and NIC96. Moreover, high-copy SNL1 suppressed the temperature sensitivity of gle2–1 and nic96-G3 mutant cells. The nic96-G3 allele was identified in a synthetic lethal genetic screen with a null allele of the closely related nucleoporin nup100. Gle2p physically associated with Nup116p in vitro, and the interaction required the N-terminal region of Nup116p. Therefore, genetic links between the role of Snl1p and at least three NPC-associated proteins were established. We suggest that Snl1p plays a stabilizing role in NPC structure and function.

Trapping of branched DNA in microfabricated structures.

We have observed electrostatic trapping of tribranched DNA molecules undergoing electrophoresis in a microfabricated pseudo-two-dimensional array of posts. Trapping occurs in a unique transport regimen in which the electrophoretic mobility is extremely sensitive to polymer topology. The arrest of branched polymers is explained by considering their center-of-mass motion; in certain conformations, owing to the constraints imposed by the obstacles a molecule cannot advance without the center of mass first moving a short distance backwards. The depth of the resulting local potential well can be much greater than the thermal energy so that escape of an immobilized molecule can be extremely slow. We summarize the expected behavior of the mobility as a function of field strength and topology and point out that the microfabricated arrays are highly suitable for detecting an extremely small number of branched molecules in a very large population of linear molecules.

Additional evidence for an eight-transmembrane-domain topology for Caenorhabditis elegans and human presenilins

Presenilins have been implicated in the genesis of Alzheimer’s disease and in facilitating LIN-12/Notch activity during development. All presenilins have multiple hydrophobic regions that could theoretically span a membrane, and a description of the membrane topology is a crucial step toward deducing the mechanism of presenilin function. Previously, we proposed an eight-transmembrane-domain model for presenilin, based on studies of the Caenorhabditis elegans SEL-12 presenilin. Here, we describe experiments that support the view that two of the hydrophobic regions of SEL-12 function as the seventh and eighth transmembrane domains. Furthermore, we have shown that human presenilin 1 behaves like SEL-12 presenilin when analyzed by our methods. Our results provide additional experimental support for the eight-transmembrane-domain model of presenilin topology.

Structure and function in rhodopsin: Topology of the C-terminal polypeptide chain in relation to the cytoplasmic loops*

Cysteine mutagenesis and site-directed spin labeling in the C-terminal region of rhodopsin have been used to probe the local structure and proximity of that region to the cytoplasmic loops. Each of the native amino acids in the sequence T335–T340 was replaced with Cys, one at a time. The sulfhydryl groups of all mutants reacted rapidly with the sulfhydryl reagent 4,4′-dithiodipyridine, which indicated a high degree of solvent accessibility. Furthermore, to probe the proximity relationships, a series of double Cys mutants was constructed. One Cys in all sets was at position 338 and the other was at a position in the sequence S240–V250 in the EF interhelical loop, at position 65 in the AB interhelical loop, or at position 140 in the CD interhelical loop. In the dark state, no significant disulfide formation was observed between C338 and C65 or C140 under the conditions used, whereas a relatively rapid disulfide formation was observed between C338 and C242 or C245. Spin labels in the double Cys mutants showed the strongest magnetic interactions between the nitroxides attached to C338 and C245 or C246. Light activation of the double mutant T242C/S338C resulted in slower disulfide formation, whereas interactions between nitroxides at C338 and C245 or C246 decreased. These results suggest the proximity of the C-terminal residue C338 to residues located on the outer face of a cytoplasmic helical extension of the F helix with an apparent increase of distance upon photoactivation.

Functional Dissection and Hierarchy of Tubulin-folding Cofactor Homologues in Fission Yeast

We describe the isolation of fission yeast homologues of tubulin-folding cofactors B (Alp11) and E (Alp21), which are essential for cell viability and the maintenance of microtubules. Alp11B contains the glycine-rich motif (the CLIP-170 domain) involved in microtubular functions, whereas, unlike mammalian cofactor E, Alp21E does not. Both mammalian and yeast cofactor E, however, do contain leucine-rich repeats. Immunoprecipitation analysis shows that Alp11B interacts with both α-tubulin and Alp21E, but not with the cofactor D homologue Alp1, whereas Alp21E also interacts with Alp1D. The cellular amount of α-tubulin is decreased in both alp1 and alp11 mutants. Overproduction of Alp11B results in cell lethality and the disappearance of microtubules, which is rescued by co-overproduction of α-tubulin. Both full-length Alp11B and the C-terminal third containing the CLIP-170 domain localize in the cytoplasm, and this domain is required for efficient binding to α-tubulin. Deletion of alp11 is suppressed by multicopy plasmids containing either alp21+ or alp1+, whereas alp21 deletion is rescued by overexpression of alp1+ but not alp11+. Finally, the alp1 mutant is not complemented by either alp11+ or alp21+. The results suggest that cofactors operate in a linear pathway (Alp11B-Alp21E-Alp1D), each with distinct roles.

Functional Hierarchy of Simultaneously Expressed Adhesion Receptors: Integrin α2β1 but Not CD44 Mediates MV3 Melanoma Cell Migration and Matrix Reorganization within Three-dimensional Hyaluronan-containing Collagen MatricesV⃞

Haptokinetic cell migration across surfaces is mediated by adhesion receptors including β1 integrins and CD44 providing adhesion to extracellular matrix (ECM) ligands such as collagen and hyaluronan (HA), respectively. Little is known, however, about how such different receptor systems synergize for cell migration through three-dimensionally (3-D) interconnected ECM ligands. In highly motile human MV3 melanoma cells, both β1 integrins and CD44 are abundantly expressed, support migration across collagen and HA, respectively, and are deposited upon migration, whereas only β1 integrins but not CD44 redistribute to focal adhesions. In 3-D collagen lattices in the presence or absence of HA and cross-linking chondroitin sulfate, MV3 cell migration and associated functions such as polarization and matrix reorganization were blocked by anti-β1 and anti-α2 integrin mAbs, whereas mAbs blocking CD44, α3, α5, α6, or αv integrins showed no effect. With use of highly sensitive time-lapse videomicroscopy and computer-assisted cell tracking techniques, promigratory functions of CD44 were excluded. 1) Addition of HA did not increase the migratory cell population or its migration velocity, 2) blocking of the HA-binding Hermes-1 epitope did not affect migration, and 3) impaired migration after blocking or activation of β1 integrins was not restored via CD44. Because α2β1-mediated migration was neither synergized nor replaced by CD44–HA interactions, we conclude that the biophysical properties of 3-D multicomponent ECM impose more restricted molecular functions of adhesion receptors, thereby differing from haptokinetic migration across surfaces.

Temperature, template topology, and factor requirements of archaeal transcription

Although Archaea are prokaryotic and resemble Bacteria morphologically, their transcription apparatus is remarkably similar to those of eukaryotic cell nuclei. Because some Archaea exist in environments with temperatures of around 100°C, they are likely to have evolved unique strategies for transcriptional control. Here, we investigate the effects of temperature and DNA template topology in a thermophilic archaeal transcription system. Significantly, and in marked contrast with characterized eucaryal systems, archaeal DNA template topology has negligible effect on transcription levels at physiological temperatures using highly purified polymerase and recombinant transcription factors. Furthermore, archaeal transcription does not require hydrolysis of the β-γ phosphoanhydride bond of ATP. However, at lower temperatures, negatively supercoiled templates are transcribed more highly than those that are positively supercoiled. Notably, the block to transcription on positively supercoiled templates at lowered temperatures is at the level of polymerase binding and promoter opening. These data imply that Archaea do not possess a functional homologue of transcription factor TFIIH, and that for the promoters studied, transcription is mediated by TATA box-binding protein, transcription factor TFB, and RNA polymerase alone. Furthermore, they suggest that the reduction of plasmid linking number by hyperthermophilic Archaea in vivo in response to cold shock is a mechanism to maintain gene expression under these adverse circumstances.

In Vivo Binding and Hierarchy of Assembly of the Yeast RNA Polymerase I Transcription Factors

Transcription by RNA polymerase I in Saccharomyces cerevisiae requires a series of transcription factors that have been genetically and biochemically identified. In particular, the core factor (CF) and the upstream activation factor (UAF) have been shown in vitro to bind the core element and the upstream promoter element, respectively. We have analyzed in vivo the DNAse I footprinting of the 35S promoter in wild-type and mutant strains lacking one specific transcription factor at the time. In this way we were able to unambiguously attribute the protections by the CF and the UAF to their respective putative binding sites. In addition, we have found that in vivo a binding hierarchy exists, the UAF being necessary for CF binding. Because the CF footprinting is lost in mutants lacking a functional RNA polymerase I, we also conclude that the final step of preinitiation-complex assembly affects binding of the CF, stabilizing its contact with DNA. Thus, in vivo, the CF is recruited to the core element by the UAF and stabilized on DNA by the presence of a functional RNA polymerase I.

Dynamical hierarchy in transition states: Why and how does a system climb over the mountain?

How a reacting system climbs through a transition state during the course of a reaction has been an intriguing subject for decades. Here we present and quantify a technique to identify and characterize local invariances about the transition state of an N-particle Hamiltonian system, using Lie canonical perturbation theory combined with microcanonical molecular dynamics simulation. We show that at least three distinct energy regimes of dynamical behavior occur in the region of the transition state, distinguished by the extent of their local dynamical invariance and regularity. Isomerization of a six-atom Lennard–Jones cluster illustrates this: up to energies high enough to make the system manifestly chaotic, approximate invariants of motion associated with a reaction coordinate in phase space imply a many-body dividing hypersurface in phase space that is free of recrossings even in a sea of chaos. The method makes it possible to visualize the stable and unstable invariant manifolds leading to and from the transition state, i.e., the reaction path in phase space, and how this regularity turns to chaos with increasing total energy of the system. This, in turn, illuminates a new type of phase space bottleneck in the region of a transition state that emerges as the total energy and mode coupling increase, which keeps a reacting system increasingly trapped in that region.

Multidrug resistance proteins QacA and QacB from Staphylococcus aureus: membrane topology and identification of residues involved in substrate specificity.

The closely related multidrug efflux pumps QacA and QacB, from the bacterial pathogen Staphylococcus aureus, both confer resistance to various toxic organic cations but differ in that QacB mediates lower levels of resistance to divalent cations. Cloning and nucleotide sequencing of the qacB gene revealed that qacB differs from qacA by only seven nucleotide substitutions. Random hydroxylamine mutagenesis of qacB was undertaken, selecting for variants that conferred increased resistance to divalent cations. Both QacA and the QacB mutants capable of conferring resistance to divalent cations contain an acidic residue at either amino acid 322 or 323, whereas QacB contains uncharged residues in these positions. Site-directed mutagenesis of qacA confirmed the importance of an acidic residue within this region of QacA in conferring resistance to divalent cations. Membrane topological analysis using alkaline phosphatase and beta-galactosidase fusions indicated that the QacA protein contains 14 transmembrane segments. Thus, QacA represents the first membrane transport protein shown to contain 14 transmembrane segments, and confirms that the major facilitator superfamily contains a family of proteins with 14 transmembrane segments.

Active site topology and reaction mechanism of GTP cyclohydrolase I.

GTP cyclohydrolase I of Escherichia coli is a torus-shaped homodecamer with D5 symmetry and catalyzes a complex ring expansion reaction conducive to the formation of dihydroneopterin triphosphate from GTP. The x-ray structure of a complex of the enzyme with the substrate analog, dGTP, bound at the active site was determined at a resolution of 3 A. In the decamer, 10 equivalent active sites are present, each of which contains a 10-A deep pocket formed by surface areas of 3 adjacent subunits. The substrate forms a complex hydrogen bond network with the protein. Active site residues were modified by site-directed mutagenesis, and enzyme activities of the mutant proteins were measured. On this basis, a mechanism of the enzyme-catalyzed reaction is proposed. Cleavage of the imidazole ring is initiated by protonation of N7 by His-179 followed by the attack of water at C8 of the purine system. Cystine Cys-110 Cys-181 may be involved in this reaction step. Opening of the imidazole ring may be in concert with cleavage of the furanose ring to generate a Schiff's base from the glycoside. The gamma-phosphate of GTP may be involved in the subsequent Amadori rearrangement of the carbohydrate side chain by activating the hydroxyl group of Ser-135.

Considerations on the folding topology and evolutionary origin of cadherin domains.

Cell-cell adhesion in zonula adherens and desmosomal junctions is mediated by cadherins, and recent crystal structures of the first domain from murine N-cadherin provide a plausible molecular basis for this adhesive action. A structure-based sequence analysis of this adhesive domain indicates that its fold is common to all extracellular cadherin domains. The cadherin folding topology is also shown to be similar to immunoglobulin-like domains and to other Greek-key beta-sandwich structures, as diverse as domains from plant cytochromes, bacterial cellulases, and eukaryotic transcription factors. Sequence similarities between cadherins and these other molecules are very low, however, and intron patterns are also different. On balance, independent origins for a favorable folding topology seem more likely than evolutionary divergence from an ancestor common to cadherins and immunoglobulins.