Systematic derivation of partition functions for ligand binding to two-dimensional lattices

The Ising problem consists in finding the analytical solution of the partition function of a lattice once the interaction geometry among its elements is specified. No general analytical solution is available for this problem, except for the one-dimensional case. Using site-specific thermodynamics, it is shown that the partition function for ligand binding to a two-dimensional lattice can be obtained from those of one-dimensional lattices with known solution. The complexity of the lattice is reduced recursively by application of a contact transformation that involves a relatively small number of steps. The transformation implemented in a computer code solves the partition function of the lattice by operating on the connectivity matrix of the graph associated with it. This provides a powerful new approach to the Ising problem, and enables a systematic analysis of two-dimensional lattices that model many biologically relevant phenomena. Application of this approach to finite two-dimensional lattices with positive cooperativity indicates that the binding capacity per site diverges as Na (N = number of sites in the lattice) and experiences a phase-transition-like discontinuity in the thermodynamic limit N → ∞. The zeroes of the partition function tend to distribute on a slightly distorted unit circle in complex plane and approach the positive real axis already for a 5×5 square lattice. When the lattice has negative cooperativity, its properties mimic those of a system composed of two classes of independent sites with the apparent population of low-affinity binding sites increasing with the size of the lattice, thereby accounting for a phenomenon encountered in many ligand-receptor interactions.

Identification of Two Type V Myosins in Fission Yeast, One of Which Functions in Polarized Cell Growth and Moves Rapidly in the Cell

We characterized the novel Schizosaccharomyces pombe genes myo4+ and myo5+, both of which encode myosin-V heavy chains. Disruption of myo4 caused a defect in cell growth and led to an abnormal accumulation of secretory vesicles throughout the cytoplasm. The mutant cells were rounder than normal, although the sites for cell polarization were still established. Elongation of the cell ends and completion of septation required more time than in wild-type cells, indicating that Myo4 functions in polarized growth both at the cell ends and during septation. Consistent with this conclusion, Myo4 was localized around the growing cell ends, the medial F-actin ring, and the septum as a cluster of dot structures. In living cells, the dots of green fluorescent protein-tagged Myo4 moved rapidly around these regions. The localization and movement of Myo4 were dependent on both F-actin cables and its motor activity but seemed to be independent of microtubules. Moreover, the motor activity of Myo4 was essential for its function. These results suggest that Myo4 is involved in polarized cell growth by moving with a secretory vesicle along the F-actin cables around the sites for polarization. In contrast, the phenotype of myo5 null cells was indistinguishable from that of wild-type cells. This and other data suggest that Myo5 has a role distinct from that of Myo4.

How optimization of potential functions affects protein folding.

The relationship between the optimization of the potential function and the foldability of theoretical protein models is studied based on investigations of a 27-mer cubic-lattice protein model and a more realistic lattice model for the protein crambin. In both the simple and the more complicated systems, optimization of the energy parameters achieves significant improvements in the statistical-mechanical characteristics of the systems and leads to foldable protein models in simulation experiments. The foldability of the protein models is characterized by their statistical-mechanical properties--e.g., by the density of states and by Monte Carlo folding simulations of the models. With optimized energy parameters, a high level of consistency exists among different interactions in the native structures of the protein models, as revealed by a correlation function between the optimized energy parameters and the native structure of the model proteins. The results of this work are relevant to the design of a general potential function for folding proteins by theoretical simulations.

Analysis of exocytosis mutants indicates close coupling between regulated secretion and transcription activation in Tetrahymena

Stimulation of regulated secretory cells promotes protein release via the fusion of cytoplasmic storage vesicles with the plasma membrane. In Tetrahymena thermophila, brief exposure to secretagogue results in synchronous fusion of the entire set of docked dense-core granules with the plasma membrane. We show that stimulation is followed by rapid new dense-core granule synthesis involving gene induction. Two genes encoding granule matrix proteins, GRL1 and GRL4, are shown to undergo induction following stimulation, resulting in ≈10-fold message accumulation within 1 h. The mechanism of induction involves transcriptional regulation, and the upstream region of GRL1 functions in vivo as an inducible promoter in a heterologous reporter construct using the gene encoding green fluorescent protein. Taking advantage of the characterized exocytosis (exo−) mutants available in this system, we asked whether the signals for regranulation were generated directly by the initial stimulation, or whether downstream events were required for transcription activation. Three mutants, with defects at three distinct stages in the regulated secretory pathway, failed to show induction of GRL1 and GRL4 after exposure to secretagogue. These results argue that regranulation depends upon signals generated by the final steps in exocytosis.

In vivo functions of the Saccharomyces cerevisiae Hsp90 chaperone

In the highly concentrated environment of the cell, polypeptide chains are prone to aggregation during synthesis (as nascent chains await the emergence of the remainder of their folding domain), translocation, assembly, and exposure to stresses that cause previously folded proteins to unfold. A large and diverse group of proteins, known as chaperones, transiently associate with such folding intermediates to prevent aggregation, but in many cases the specific functions of individual chaperones are still not clear. In vivo, Hsp90 (heat shock protein 90) plays a role in the maturation of components of signal transduction pathways but also exhibits chaperone activity with diverse proteins in vitro, suggesting a more general function. We used a unique temperature-sensitive mutant of Hsp90 in Saccharomyces cerevisiae, which rapidly and completely loses activity on shift to high temperatures, to examine the breadth of Hsp90 functions in vivo. The data suggest that Hsp90 is not required for the de novo folding of most proteins, but it is required for a specific subset of proteins that have greater difficulty reaching their native conformations. Under conditions of stress, Hsp90 does not generally protect proteins from thermal inactivation but does enhance the rate at which a heat-damaged protein is reactivated. Thus, although Hsp90 is one of the most abundant chaperones in the cell, its in vivo functions are highly restricted.

Identification of a novel p53 functional domain that is necessary for efficient growth suppression

Activation of the p53 tumor suppressor protein has been demonstrated to block cell growth by inducing either a transient cell cycle arrest or programmed cell death (apoptosis). Although evidence exists linking p53’s function as an activator of transcription to its ability to effect cell cycle arrest, the role of this activity in the induction of apoptosis remains unclear. To gain insight into the molecular mechanisms underlying p53-mediated antiproliferative pathways, a study was initiated to explore the functions of a putative p53 signaling domain. This region of the human p53 protein is localized between amino acids 61 and 94 (out of 393) and is noteworthy in that it contains five repeats of the sequence PXXP (where P represents proline and X any amino acid). This motif has been shown to play a role in signal transduction via its SH3 domain binding activity. A p53 cDNA deletion mutant (ΔproAE), which lacks this entire proline-rich domain (deleted for amino acids 62–91), was created and characterized for a variety of p53 functions. The entire domain has been shown to be completely dispensable for transcriptional activation. On the other hand, this deletion of the p53 proline-rich domain impairs p53’s ability to suppress tumor cell growth in culture. Amino acid substitution mutations at residues 22 and 23 of p53 (eliminates transcriptional activity) also impair p53-mediated inhibition of cell growth in culture. Unlike wild-type p53, the ΔproAE mutant cDNA can be stably expressed in tumor derived cell lines with few immediate detrimental effects. These cells express physiologic levels of p53 protein that are induced normally in response to DNA damage, indicating that removal of the proline-rich domain does not disrupt p53’s upstream regulation by DNA damage. These data indicate that, in addition to the transcriptional activation domain, the p53 proline-rich domain plays a critical role in the transmission of antiproliferative signals downstream of the p53 protein and may link p53 to a direct signal transduction pathway.

Advances in random matrix theory, zeta functions, and sphere packing

Over four hundred years ago, Sir Walter Raleigh asked his mathematical assistant to find formulas for the number of cannonballs in regularly stacked piles. These investigations aroused the curiosity of the astronomer Johannes Kepler and led to a problem that has gone centuries without a solution: why is the familiar cannonball stack the most efficient arrangement possible? Here we discuss the solution that Hales found in 1998. Almost every part of the 282-page proof relies on long computer verifications. Random matrix theory was developed by physicists to describe the spectra of complex nuclei. In particular, the statistical fluctuations of the eigenvalues (“the energy levels”) follow certain universal laws based on symmetry types. We describe these and then discuss the remarkable appearance of these laws for zeros of the Riemann zeta function (which is the generating function for prime numbers and is the last special function from the last century that is not understood today.) Explaining this phenomenon is a central problem. These topics are distinct, so we present them separately with their own introductory remarks.

Optical imaging of functional domains in the cortex of the awake and behaving monkey

As demonstrated by anatomical and physiological studies, the cerebral cortex consists of groups of cortical modules, each comprising populations of neurons with similar functional properties. This functional modularity exists in both sensory and association neocortices. However, the role of such cortical modules in perceptual and cognitive behavior is unknown. To aid in the examination of this issue we have applied the high spatial resolution optical imaging methodology to the study of awake, behaving animals. In this paper, we report the optical imaging of orientation domains and blob structures, approximately 100–200 μm in size, in visual cortex of the awake and behaving monkey. By overcoming the spatial limitations of other existing imaging methods, optical imaging will permit the study of a wide variety of cortical functions at the columnar level, including motor and cognitive functions traditionally studied with positron-emission tomography or functional MRI techniques.

PRIM: Proximity imaging of green fluorescent protein-tagged polypeptides

We report a serendipitous discovery that extends the impressive catalog of reporter functions performed by green fluorescent protein (GFP) or its derivatives. When two GFP molecules are brought into proximity, changes in the relative intensities of green fluorescence emitted upon excitation at 395 vs. 475 nm result. These spectral changes provide a sensitive ratiometric index of the extent of self-association that can be exploited to quantitatively image homo-oligomerization or clustering processes of GFP-tagged proteins in vivo. The method, which we term proximity imaging (PRIM), complements fluorescence resonance energy transfer between a blue fluorescent protein donor and a GFP acceptor, a powerful method for imaging proximity relationships between different proteins. However, unlike fluorescence resonance energy transfer (which is a spectral interaction), PRIM depends on direct contact between two GFP modules, which can lead to structural perturbations and concomitant spectral changes within a module. Moreover, the precise spatial arrangement of the GFP molecules within a given dimer determines the magnitude and direction of the spectral change. We have used PRIM to detect FK1012-induced dimerization of GFP fused to FK506-binding protein and clustering of glycosylphosphatidylinositol-anchored GFP at cell surfaces.

GIPC, a PDZ domain containing protein, interacts specifically with the C terminus of RGS-GAIP

We have identified a mammalian protein called GIPC (for GAIP interacting protein, C terminus), which has a central PDZ domain and a C-terminal acyl carrier protein (ACP) domain. The PDZ domain of GIPC specifically interacts with RGS-GAIP, a GTPase-activating protein (GAP) for Gαi subunits recently localized on clathrin-coated vesicles. Analysis of deletion mutants indicated that the PDZ domain of GIPC specifically interacts with the C terminus of GAIP (11 amino acids) in the yeast two-hybrid system and glutathione S-transferase (GST)-GIPC pull-down assays, but GIPC does not interact with other members of the RGS (regulators of G protein signaling) family tested. This finding is in keeping with the fact that the C terminus of GAIP is unique and possesses a modified C-terminal PDZ-binding motif (SEA). By immunoblotting of membrane fractions prepared from HeLa cells, we found that there are two pools of GIPC–a soluble or cytosolic pool (70%) and a membrane-associated pool (30%). By immunofluorescence, endogenous and GFP-tagged GIPC show both a diffuse and punctate cytoplasmic distribution in HeLa cells reflecting, respectively, the existence of soluble and membrane-associated pools. By immunoelectron microscopy the membrane pool of GIPC is associated with clusters of vesicles located near the plasma membrane. These data provide direct evidence that the C terminus of a RGS protein is involved in interactions specific for a given RGS protein and implicates GAIP in regulation of additional functions besides its GAP activity. The location of GIPC together with its binding to GAIP suggest that GAIP and GIPC may be components of a G protein-coupled signaling complex involved in the regulation of vesicular trafficking. The presence of an ACP domain suggests a putative function for GIPC in the acylation of vesicle-bound proteins.

Opening of a monomer-monomer interface of the trimeric bacteriophage T4-coded GP45 sliding clamp is required for clamp loading onto DNA

The replication system of bacteriophage T4 uses a trimeric ring-shaped processivity clamp (gp45) to tether the replication polymerase (gp43) to the template-primer DNA. This ring is placed onto the DNA by an ATPase-driven clamp-loading complex (gp44/62) where it then transfers, in closed form, to the polymerase. It generally has been assumed that one of the functions of the loading machinery is to open the clamp to place it around the DNA. However, the mechanism by which this occurs has not been fully defined. In this study we design and characterize a double-mutant gp45 protein that contains pairs of cysteine residues located at each monomer-monomer interface of the trimeric clamp. This mutant protein is functionally equivalent to wild-type gp45. However, when all three monomer-monomer interfaces are tethered by covalent crosslinks formed (reversibly or irreversibly) between the cysteine pairs these closed clamps can no longer be loaded onto the DNA nor onto the polymerase, effectively eliminating processive strand-displacement DNA synthesis. Analysis of the individual steps of the clamp-loading process shows that the ATPase-dependent interactions between the clamp and the clamp loader that precede DNA binding are hyperstimulated by the covalently crosslinked ring, suggesting that binding of the closed ring induces a futile, ATP-driven, ring-opening cycle. These findings and others permit further characterization and ordering of the steps involved in the T4 clamp-loading process.

The mouse pale ear (ep) mutation is the homologue of human Hermansky–Pudlak syndrome

The recessive mutation at the pale ear (ep) locus on mouse chromosome 19 was found to be the homologue of human Hermansky–Pudlak syndrome (HPS). A positional cloning strategy using yeast artificial chromosomes spanning the HPS locus was used to identify the HPS gene and its murine counterpart. These genes and their predicted proteins are highly conserved at the nucleotide and amino acid levels. Sequence analysis of the mutant ep gene revealed the insertion of an intracisternal A particle element in a protein-coding 3′ exon. Here we demonstrate that mice with the ep mutation exhibit abnormalities similar to human HPS patients in melanosomes and platelet-dense granules. These results establish an animal model of HPS and will facilitate biochemical and molecular analyses of the functions of this protein in the membranes of specialized intracellular organelles.

Interactions among enzymes of the Arabidopsis flavonoid biosynthetic pathway

Flavonoids are secondary metabolites derived from phenylalanine and acetate metabolism that perform a variety of essential functions in higher plants. Studies over the past 30 years have supported a model in which flavonoid metabolism is catalyzed by an enzyme complex localized to the endoplasmic reticulum [Hrazdina, G. & Wagner, G. J. (1985) Arch. Biochem. Biophys. 237, 88–100]. To test this model further we assayed for direct interactions between several key flavonoid biosynthetic enzymes in developing Arabidopsis seedlings. Two-hybrid assays indicated that chalcone synthase, chalcone isomerase (CHI), and dihydroflavonol 4-reductase interact in an orientation-dependent manner. Affinity chromatography and immunoprecipitation assays further demonstrated interactions between chalcone synthase, CHI, and flavonol 3-hydroxylase in lysates from Arabidopsis seedlings. These results support the hypothesis that the flavonoid enzymes assemble as a macromolecular complex with contacts between multiple proteins. Evidence was also found for posttranslational modification of CHI. The importance of understanding the subcellular organization of elaborate enzyme systems is discussed in the context of metabolic engineering.

β-Catenin mutations in cell lines established from human colorectal cancers

β-catenin has functions as both an adhesion and a signaling molecule. Disruption of these functions through mutations of the β-catenin gene (CTNNB1) may be important in the development of colorectal tumors. We examined the entire coding sequence of β-catenin by reverse transcriptase–PCR (RT-PCR) and direct sequencing of 23 human colorectal cancer cell lines from 21 patients. In two cell lines, there was apparent instability of the β-catenin mRNA. Five different mutations (26%) were found in the remaining 21cell lines (from 19 patients). A three-base deletion (codon 45) was identified in the cell line HCT 116, whereas cell lines SW 48, HCA 46, CACO 2, and Colo 201 each contained single-base missense mutations (codons 33, 183, 245, and 287, respectively). All 23 cell lines had full-length β-catenin protein that was detectable by Western blotting and that coprecipitated with E-cadherin. In three of the cell lines with CTNNB1 mutations, complexes of β-catenin with α-catenin and APC were detectable. In SW48 and HCA 46, however, we did not detect complexes of β-catenin protein with α-catenin and APC, respectively. These results show that selection of CTNNB1 mutations occurs in up to 26% of colorectal cancers from which cell lines are derived. In these cases, mutation selection is probably for altered β-catenin function, which may significantly alter intracellular signaling and intercellular adhesion and may serve as a complement to APC mutations in the early stages of tumorigenesis.

Probing the assembly of transcription initiation complexes through changes in σN protease sensitivity

The alternative bacterial σN RNA polymerase holoenzyme binds promoters as a transcriptionally inactive complex that is activated by enhancer-binding proteins. Little is known about how sigma factors respond to their ligands or how the responses lead to transcription. To examine the liganded state of σN, the assembly of end-labeled Klebsiella pneumoniae σN into holoenzyme, closed promoter complexes, and initiated transcription complexes was analyzed by enzymatic protein footprinting. V8 protease-sensitive sites in free σN were identified in the acidic region II and bordering or within the minimal DNA binding domain. Interaction with core RNA polymerase prevented cleavage at noncontiguous sites in region II and at some DNA binding domain sites, probably resulting from conformational changes. Formation of closed complexes resulted in further protections within the DNA binding domain, suggesting close contact to promoter DNA. Interestingly, residue E36 becomes sensitive to proteolysis in initiated transcription complexes, indicating a conformational change in holoenzyme during initiation. Residue E36 is located adjacent to an element involved in nucleating strand separation and in inhibiting polymerase activity in the absence of activation. The sensitivity of E36 may reflect one or both of these functions. Changing patterns of protease sensitivity strongly indicate that σN can adjust conformation upon interaction with ligands, a property likely important in the dynamics of the protein during transcription initiation.