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.

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.

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.

Design of three-dimensional domain-swapped dimers and fibrous oligomers

Three-dimensional (3D) domain-swapped proteins are intermolecularly folded analogs of monomeric proteins; both are stabilized by the identical interactions, but the individual domains interact intramolecularly in monomeric proteins, whereas they form intermolecular interactions in 3D domain-swapped structures. The structures and conditions of formation of several domain-swapped dimers and trimers are known, but the formation of higher order 3D domain-swapped oligomers has been less thoroughly studied. Here we contrast the structural consequences of domain swapping from two designed three-helix bundles: one with an up-down-up topology, and the other with an up-down-down topology. The up-down-up topology gives rise to a domain-swapped dimer whose structure has been determined to 1.5 Å resolution by x-ray crystallography. In contrast, the domain-swapped protein with an up-down-down topology forms fibrils as shown by electron microscopy and dynamic light scattering. This demonstrates that design principles can predict the oligomeric state of 3D domain-swapped molecules, which should aid in the design of domain-swapped proteins and biomaterials.

Measuring the topology of the universe

Observations of microwave background fluctuations can yield information not only about the geometry of the universe but potentially about the topology of the universe. If the universe is negatively curved, then the characteristic scale for the topology of the universe is the curvature radius. Thus, if we are seeing the effects of the geometry of the universe, we can hope to soon see signatures of the topology of the universe. The cleanest signature of the topology of the universe is written on the microwave sky: There should be thousands of pairs of matched circles. These circles can be used to determine the precise topology and volume of the universe. Because we see hundreds of slices through the fundamental domain of the universe, we can use the microwave observations to reconstruct the initial conditions of the entire universe on the scale of a few megaparsecs.

Retroposon analysis of major cetacean lineages: The monophyly of toothed whales and the paraphyly of river dolphins

SINE (short interspersed element) insertion analysis elucidates contentious aspects in the phylogeny of toothed whales and dolphins (Odontoceti), especially river dolphins. Here, we characterize 25 informative SINEs inserted into unique genomic loci during evolution of odontocetes to construct a cladogram, and determine a total of 2.8 kb per taxon of the flanking sequences of these SINE loci to estimate divergence times among lineages. We demonstrate that: (i) Odontocetes are monophyletic; (ii) Ganges River dolphins, beaked whales, and ocean dolphins diverged (in this order) after sperm whales; (iii) three other river dolphin taxa, namely the Amazon, La Plata, and Yangtze river dolphins, form a monophyletic group with Yangtze River dolphins being the most basal; and (iv) the rapid radiation of extant cetacean lineages occurred some 28–33 million years B.P., in strong accord with the fossil record. The combination of SINE and flanking sequence analysis suggests a topology and set of divergence times for odontocete relationships, offering alternative explanations for several long-standing problems in cetacean evolution.

Deficiency of rds/peripherin causes photoreceptor death in mouse models of digenic and dominant retinitis pigmentosa

Retinitis pigmentosa (RP) is a group of inherited blinding diseases caused by mutations in multiple genes including RDS. RDS encodes rds/peripherin (rds), a 36-kDa glycoprotein in the rims of rod and cone outer-segment (OS) discs. Rom1 is related to rds with similar membrane topology and the identical distribution in OS. In contrast to RDS, no mutations in ROM1 alone have been associated with retinal disease. However, an unusual digenic form of RP has been described. Affected individuals in several families were doubly heterozygous for a mutation in RDS causing a leucine 185 to proline substitution in rds (L185P) and a null mutation in ROM1. Neither mutation alone caused clinical abnormalities. Here, we generated transgenic/knockout mice that duplicate the amino acid substitutions and predicted levels of rds and rom1 in patients with RDS-mediated digenic and dominant RP. Photoreceptor degeneration in the mouse model of digenic RP was faster than in the wild-type and monogenic controls by histological, electroretinographic, and biochemical analysis. We observed a positive correlation between the rate of photoreceptor loss and the extent of OS disorganization in mice of several genotypes. Photoreceptor degeneration in RDS-mediated RP appears to be caused by a simple deficiency of rds and rom1. The critical threshold for the combined abundance of rds and rom1 is ≈60% of wild type. Below this value, the extent of OS disorganization results in clinically significant photoreceptor degeneration.

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.

Cells that register logical relationships among proteins

Two-hybrid methods have augmented the classical genetic techniques biologists use to assign function to genes. Here, we describe construction of a two-bait interaction trap that uses yeast cells to register more complex protein relationships than those detected in existing two-hybrid systems. We show that such cells can identify bridge or connecting proteins and peptide aptamers that discriminate between closely related allelic variants. The protein relationships detected by these cells are analogous to classical genetic relationships, but lend themselves to systematic application to the products of entire genomes and combinatorial libraries. We show that, by performing logical operations on the phenotypic outputs of these complex cells and existing two-hybrid cells, we can make inferences about the topology and order of protein interactions. Finally, we show that cells that register such relationships can perform logical operations on protein inputs. Thus these cells will be useful for analysis of gene and allele function, and may also define a path for construction of biological computational devices.

Topological testing of the mechanism of homology search promoted by RecA protein

To initiate homologous recombination, sequence similarity between two DNA molecules must be searched for and homology recognized. How the search for and recognition of homology occurs remains unproven. We have examined the influences of DNA topology and the polarity of RecA–single-stranded (ss)DNA filaments on the formation of synaptic complexes promoted by RecA. Using two complementary methods and various ssDNA and duplex DNA molecules as substrates, we demonstrate that topological constraints on a small circular RecA–ssDNA filament prevent it from interwinding with its duplex DNA target at the homologous region. We were unable to detect homologous pairing between a circular RecA–ssDNA filament and its relaxed or supercoiled circular duplex DNA targets. However, the formation of synaptic complexes between an invading linear RecA–ssDNA filament and covalently closed circular duplex DNAs is promoted by supercoiling of the duplex DNA. The results imply that a triplex structure formed by non-Watson–Crick hydrogen bonding is unlikely to be an intermediate in homology searching promoted by RecA. Rather, a model in which RecA-mediated homology searching requires unwinding of the duplex DNA coupled with local strand exchange is the likely mechanism. Furthermore, we show that polarity of the invading RecA–ssDNA does not affect its ability to pair and interwind with its circular target duplex DNA.

Solution structure of the antiapoptotic protein bcl-2

The structures of two isoforms of Bcl-2 that differ by two amino acids have been determined by NMR spectroscopy. Because wild-type Bcl-2 behaved poorly in solution, the structures were determined by using Bcl-2/Bcl-xL chimeras in which part of the putative unstructured loop of Bcl-2 was replaced with a shortened loop from Bcl-xL. These chimeric proteins have a low pI compared with the wild-type protein and are soluble. The structures of the two Bcl-2 isoforms consist of 6 α-helices with a hydrophobic groove on the surface similar to that observed for the homologous protein, Bcl-xL. Comparison of the Bcl-2 structures to that of Bcl-xL shows that although the overall fold is the same, there are differences in the structural topology and electrostatic potential of the binding groove. Although the structures of the two isoforms of Bcl-2 are virtually identical, differences were observed in the ability of the proteins to bind to a 25-residue peptide from the proapoptotic Bad protein and a 16-residue peptide from the proapoptotic Bak protein. These results suggest that there are subtle differences in the hydrophobic binding groove in Bcl-2 that may translate into differences in antiapoptotic activity for the two isoforms.

A brief survey of symmetry in mathematics

This paper presents a brief survey of the idea of symmetry in mathematics, as exemplified by some particular developments in algebra, differential equations, topology, and number theory.

Rapid refolding of a proline-rich all-beta-sheet fibronectin type III module.

Fibronectin type III modules contain approximately 90 residues and are an extremely common building block of animal proteins. Despite containing a complex all-beta-sheet topology and eight prolines, the refolding of the 10th type III module of human fibronectin has been found to be very rapid, with native core packing, amide hydrogen bonding, and backbone conformation all recovered within 1 s at 5 degrees C. These observations indicate that this domain can overcome many structural characteristics often thought to slow the folding process.