Evidence for nonrandom hydrophobicity structures in protein chains.

The question of whether proteins originate from random sequences of amino acids is addressed. A statistical analysis is performed in terms of blocked and random walk values formed by binary hydrophobic assignments of the amino acids along the protein chains. Theoretical expectations of these variables from random distributions of hydrophobicities are compared with those obtained from functional proteins. The results, which are based upon proteins in the SWISS-PROT data base, convincingly show that the amino acid sequences in proteins differ from what is expected from random sequences in a statistically significant way. By performing Fourier transforms on the random walks, one obtains additional evidence for nonrandomness of the distributions. We have also analyzed results from a synthetic model containing only two amino acid types, hydrophobic and hydrophilic. With reasonable criteria on good folding properties in terms of thermodynamical and kinetic behavior, sequences that fold well are isolated. Performing the same statistical analysis on the sequences that fold well indicates similar deviations from randomness as for the functional proteins. The deviations from randomness can be interpreted as originating from anticorrelations in terms of an Ising spin model for the hydrophobicities. Our results, which differ from some previous investigations using other methods, might have impact on how permissive with respect to sequence specificity protein folding process is-only sequences with nonrandom hydrophobicity distributions fold well. Other distributions give rise to energy landscapes with poor folding properties and hence did not survive the evolution.

The primary structures of the Archaeon Halobacterium salinarium blue light receptor sensory rhodopsin II and its transducer, a methyl-accepting protein.

Recently, a large family of transducer proteins in the Archaeon Halobacterium salinarium was identified. On the basis of the comparison of the predicted structural domains of these transducers, three distinct subfamilies of transducers were proposed. Here we report isolation, complete gene sequences, and analysis of the encoded primary structures of transducer gene htrII, a member of family B, and its blue light receptor gene (sopII) of sensory rhodopsin II (SRII). The start codon ATG of the 714-bp sopII gene is one nucleotide beyond the termination codon TGA of the 2298-bp htrII gene. The deduced protein sequence of HtrII predicts a eubacterial chemotaxis transducer type with two hydrophobic membrane-spanning segments connecting sizable domains in the periplasm and cytoplasm. HtrII has a common feature with HtrI, the sensory rhodopsin I transducer; like HtrI, HtrII possesses a hydrophilic loop structure just after the second transmembrane segment. The C-terminal 299 residues (765 amino acid residues total) of HtrII show strong homology to the signaling and methylation domain of eubacterial transducer Tsr. The hydropathy plot of the primary structure of SRII indicates seven membrane-spanning alpha-helical segments, a characteristic feature of retinylidene proteins ("rhodopsins") from a widespread family of photoactive pigments. SRII shows high identity with SRI (42%), bacteriorhodopsin (BR) (32%), and halorhodopsin (24%). The crucial positions for retinal binding sites in these proteins are nearly identical, with the exception of Met-118 (numbering according to the mature BR sequence), which is replaced by Val in SRII. In BR, residues Asp-85 and Asp-96 are crucial in proton pumping. In SRII, the position corresponding to Asp-85 in BR is conserved, but the corresponding position of Asp-96 is replaced by an aromatic Tyr. Coexpression of the htrII and sopII genes restores SRII phototaxis to a mutant (Pho81) that contains a deletion in the htrI/sopI and insertion in htrII/sopII regions. This paper describes the first example that both HtrI and HtrII exist in the same halobacterial cell, confirming that different sensory rhodopsins SRI and SRII in the same organism have their own distinct transducers.

Three quaternary structures for a single protein.

The structure of a multisubunit protein (immunoglobulin light chain) was solved in three crystal forms, differing only in the solvent of crystallization. The three structures were obtained at high ionic strength and low pH, high ionic strength and high pH, and low ionic strength and neutral pH. The three resulting "snapshots" of possible structures show that their variable-domain interactions differ, reflecting their stabilities under specific solvent conditions. In the three crystal forms, the variable domains had different rotational and translational relationships, whereas no alteration of the constant domains was found. The critical residues involved in the observed effect of the solvent are tryptophans and histidines located between the two variable domains in the dimeric structure. Tryptophan residues are commonly found in interfaces between proteins and their subunits, and histidines have been implicated in pH-dependent conformation changes. The quaternary structure observed for a multisubunit protein or protein complex in a crystal may be influenced by the interactions of the constituents within the molecule or complex and/or by crystal packing interactions. The comparison of buried surface areas and hydrogen bonds between the domains forming the molecule and between the molecules forming the crystals suggest that, for this system, the interactions within the molecule are most likely the determining factors.

Crystal structures of the trimeric human immunodeficiency virus type 1 matrix protein: implications for membrane association and assembly.

The human immunodeficiency virus type 1 (HIV-1) matrix protein forms a structural shell associated with the inner viral membrane and performs other essential functions throughout the viral life cycle. The crystal structure of the HIV-1 matrix protein, determined at 2.3 angstrom resolution, reveals that individual matrix molecules are composed of five major helices capped by a three-stranded mixed beta-sheet. Unexpectedly, the protein assembles into a trimer in three different crystal lattices, burying 1880 angstrom2 of accessible surface area at the trimer interfaces. Trimerization appears to create a large, bipartite membrane binding surface in which exposed basic residues could cooperate with the N-terminal myristoyl groups to anchor the protein on the acidic inner membrane of the virus.

Visualization of intermediate and transition-state structures in protein-tyrosine phosphatase catalysis.

Engineering site-specific amino acid substitutions into the protein-tyrosine phosphatase (PTPase) PTP1 and the dual-specific vaccinia H1-related phosphatase (VHR), has kinetically isolated the two chemical steps of the reaction and provided a rare opportunity for examining transition states and directly observing the phosphoenzyme intermediate. Changing serine to alanine in the active-site sequence motif HCXXGXXRS shifted the rate-limiting step from intermediate formation to intermediate hydrolysis. Using phosphorus 31P NMR, the covalent thiol-phosphate intermediate was directly observed during catalytic turnover. The importance of the conserved aspartic acid (D92 in VHR and D181 in PTP1) in both chemical steps was established. Kinetic analysis of D92N and D181N mutants indicated that aspartic acid acts as a general acid by protonating the leaving-group phenolic oxygen. Structure-reactivity experiments with native and aspartate mutant enzymes established that proton transfer is concomitant with P-O cleavage, such that no charge develops on the phenolic oxygen. Steady- and presteady-state kinetics, as well as NMR analysis of the double mutant D92N/S131A (VHR), suggested that the conserved aspartic acid functions as a general base during intermediate hydrolysis. As a general base, aspartate would activate a water molecule to facilitate nucleophilic attack. The amino acids involved in transition-state stabilization for cysteinylphosphate hydrolysis were confirmed by the x-ray structure of the Yersinia PTPase complexed with vanadate, a transition-state mimic that binds covalently to the active-site cysteine. Consistent with the NMR, x-ray, biochemical, and kinetic data, a unifying mechanism for catalysis is proposed.

Core binding factor beta-smooth muscle myosin heavy chain chimeric protein involved in acute myeloid leukemia forms unusual nuclear rod-like structures in transformed NIH 3T3 cells.

Patients with the M4Eo subtype of acute myeloid leukemia almost invariably are found to have an inversion of chromosome 16 in their leukemic cells, which results in a gene fusion between the transcription factor called core binding factor beta (CBFbeta) on 16q and a smooth muscle myosin heavy chain (SMMHC) gene on 16p. Subcellular localizations of the wild-type CBFbeta and the CBFbeta-SMMHC fusion protein were determined by immunofluorescence of NIH 3T3 cells that overexpress wild-type or fusion protein. Normal CBFbeta showed an unexpected perinuclear pattern consistent with primary localization in the Golgi complex. The CBFbeta-SMMHC fusion protein had a very different pattern. Nuclear staining included rod-like crystalline structures as long as 11 microm. The heterodimeric partner of CBFbeta, CBFalpha, formed part of this complex. Cytoplasmic staining included stress fibers that colocalized with actin, probably as a consequence of the myosin heavy chain component of the fusion protein. Deletion of different regions of the CBFbeta portion of the fusion protein showed that binding to CBFalpha was not required for nuclear translocation. However, deletion of parts of the SMMHC domain of the fusion protein involved in myosin-mediated filament formation resulted in proteins that did not form rod-like structures. These observations confirm previous indirect evidence that the CBFbeta-SMMHC fusion protein is capable of forming macromolecular nuclear aggregates and suggests possible models for the mechanism of leukemic transformation.

Definition of amide protection factors for early kinetic intermediates in protein folding

Hydrogen–deuterium exchange experiments have been used previously to investigate the structures of well defined states of a given protein. These include the native state, the unfolded state, and any intermediates that can be stably populated at equilibrium. More recently, the hydrogen–deuterium exchange technique has been applied in kinetic labeling experiments to probe the structures of transiently formed intermediates on the kinetic folding pathway of a given protein. From these equilibrium and nonequilibrium studies, protection factors are usually obtained. These protection factors are defined as the ratio of the rate of exchange of a given backbone amide when it is in a fully solvent-exposed state (usually obtained from model peptides) to the rate of exchange of that amide in some state of the protein or in some intermediate on the folding pathway of the protein. This definition is straightforward for the case of equilibrium studies; however, it is less clear-cut for the case of transient kinetic intermediates. To clarify the concept for the case of burst-phase intermediates, we have introduced and mathematically defined two different types of protection factors: one is Pstruc, which is more related to the structure of the intermediate, and the other is Papp, which is more related to the stability of the intermediate. Kinetic hydrogen–deuterium exchange data from disulfide-intact ribonuclease A and from cytochrome c are discussed to explain the use and implications of these two definitions.

Region-dependent dynamics of cAMP response element-binding protein phosphorylation in the basal ganglia

The cAMP response element-binding protein (CREB) is an activity-dependent transcription factor that is involved in neural plasticity. The kinetics of CREB phosphorylation have been suggested to be important for gene activation, with sustained phosphorylation being associated with downstream gene expression. If so, the duration of CREB phosphorylation might serve as an indicator for time-sensitive plastic changes in neurons. To screen for regions potentially involved in dopamine-mediated plasticity in the basal ganglia, we used organotypic slice cultures to study the patterns of dopamine- and calcium-mediated CREB phosphorylation in the major subdivisions of the striatum. Different durations of CREB phosphorylation were evoked in the dorsal and ventral striatum by activation of dopamine D1-class receptors. The same D1 stimulus elicited (i) transient phosphorylation (≤15 min) in the matrix of the dorsal striatum; (ii) sustained phosphorylation (≤2 hr) in limbic-related structures including striosomes, the nucleus accumbens, the fundus striati, and the bed nucleus of the stria terminalis; and (iii) prolonged phosphorylation (up to 4 hr or more) in cellular islands in the olfactory tubercle. Elevation of Ca2+ influx by stimulation of L-type Ca2+ channels, NMDA, or KCl induced strong CREB phosphorylation in the dorsal striatum but not in the olfactory tubercle. These findings differentiate the response of CREB to dopamine and calcium signals in different striatal regions and suggest that dopamine-mediated CREB phosphorylation is persistent in limbic-related regions of the neonatal basal ganglia. The downstream effects activated by persistent CREB phosphorylation may include time-sensitive neuroplasticity modulated by dopamine.

Folding protein models with a simple hydrophobic energy function: The fundamental importance of monomer inside/outside segregation

The present study explores a “hydrophobic” energy function for folding simulations of the protein lattice model. The contribution of each monomer to conformational energy is the product of its “hydrophobicity” and the number of contacts it makes, i.e., E(h⃗, c⃗) = −Σi=1N cihi = −(h⃗.c⃗) is the negative scalar product between two vectors in N-dimensional cartesian space: h⃗ = (h1, … , hN), which represents monomer hydrophobicities and is sequence-dependent; and c⃗ = (c1, … , cN), which represents the number of contacts made by each monomer and is conformation-dependent. A simple theoretical analysis shows that restrictions are imposed concomitantly on both sequences and native structures if the stability criterion for protein-like behavior is to be satisfied. Given a conformation with vector c⃗, the best sequence is a vector h⃗ on the direction upon which the projection of c⃗ − c̄⃗ is maximal, where c̄⃗ is the diagonal vector with components equal to c̄, the average number of contacts per monomer in the unfolded state. Best native conformations are suggested to be not maximally compact, as assumed in many studies, but the ones with largest variance of contacts among its monomers, i.e., with monomers tending to occupy completely buried or completely exposed positions. This inside/outside segregation is reflected on an apolar/polar distribution on the corresponding sequence. Monte Carlo simulations in two dimensions corroborate this general scheme. Sequences targeted to conformations with large contact variances folded cooperatively with thermodynamics of a two-state transition. Sequences targeted to maximally compact conformations, which have lower contact variance, were either found to have degenerate ground state or to fold with much lower cooperativity.

Complete resolution of the solid-state NMR spectrum of a uniformly 15N-labeled membrane protein in phospholipid bilayers

Complete resolution of the amide resonances in a three-dimensional solid-state NMR correlation spectrum of a uniformly 15N-labeled membrane protein in oriented phospholipid bilayers is demonstrated. The three orientationally dependent frequencies, 1H chemical shift, 1H–15N dipolar coupling, and 15N chemical shift, associated with each amide resonance are responsible for resolution among resonances and provide sufficient angular restrictions for protein structure determination. Because the protein is completely immobilized by the phospholipids on the relevant NMR time scales (10 kHz), the linewidths will not degrade in the spectra of larger proteins. Therefore, these results demonstrate that solid-state NMR experiments can overcome the correlation time problem and extend the range of proteins that can have their structures determined by NMR spectroscopy to include uniformly 15N-labeled membrane proteins in phospholipid bilayers.

mRNA binding protein mrnp 41 localizes to both nucleus and cytoplasm

We have identified and molecularly characterized a human protein with a Mr of 40,880 Da. After UV irradiation of HeLa cells, this protein was cross-linked to poly(A)-containing mRNA and was therefore designated mrnp 41 (for mRNA binding protein of 41 kDa). Cell fractionation and immunoblotting showed mrnp 41 in both the cytoplasm and the nucleus and particularly in the nuclear envelope. Immunofluorescence microscopy localized mrnp 41 to distinct foci in the nucleoplasm, to the nuclear rim, and to meshwork-like structures throughout the cytoplasm. The cytoplasmic meshwork staining was disrupted by prior treatment of cells with the actin filament- or microtubule-disrupting drugs cytochalasin or nocodazole, respectively, suggesting association of mrnp 41 with the cytoskeleton. Double immunofluorescence with antibodies against mrnp 41 and the cytoplasmic poly(A) binding protein showed colocalization to the cytoplasmic meshwork. Immunogold electronmicroscopy confirmed mrnp 41’s cytoplasmic and nucleoplasmic localization and revealed a striking labeling of nuclear pore complexes. Together these data suggest that mrnp 41 may function in nuclear export of mRNPs and/or in cytoplasmic transport on, or attachment to, the cytoskeleton. Consistent with a role of mrnp 41 in nuclear export are previous reports that mutations in homologs of mrnp 41 in Schizosaccharomyces pombe, designated Rae1p, or in Saccharomyces cerevisiae, designated Gle2p, result in mRNA accumulation in the nucleus although it is presently not known whether these homologs are mRNA binding proteins as well.

Linking topography of its potential surface with the dynamics of folding of a protein model

The “3-color, 46-bead” model of a folding polypeptide is the vehicle for adapting to proteins a mode of analysis used heretofore for atomic clusters, to relate the topography of the potential surface to the dynamics that lead to formation of selected structures. The analysis is based on sequences of stationary points—successive minima, joined by saddles—that rise monotonically in energy from basin bottoms. Like structure-seeking clusters, the potential surface of the model studied here is staircase-like, rather than sawtooth-like, with highly collective motions required for passage from one minimum to the next. The surface has several deep basins whose minima correspond to very similar structures, but which are separated by high energy barriers.

Crystal structures of the copper and nickel complexes of RNase A: Metal-induced interprotein interactions and identification of a novel copper binding motif

We report the crystal structures of the copper and nickel complexes of RNase A. The overall topology of these two complexes is similar to that of other RNase A structures. However, there are significant differences in the mode of binding of copper and nickel. There are two copper ions per molecule of the protein, but there is only one nickel ion per molecule of the protein. Significant changes occur in the interprotein interactions as a result of differences in the coordinating groups at the common binding site around His-105. Consequently, the copper- and nickel-ion-bound dimers of RNase A act as nucleation sites for generating different crystal lattices for the two complexes. A second copper ion is present at an active site residue His-119 for which all the ligands are from one molecule of the protein. At this second site, His-119 adopts an inactive conformation (B) induced by the copper. We have identified a novel copper binding motif involving the α-amino group and the N-terminal residues.

A possible role for kinase-associated protein phosphatase in the Arabidopsis CLAVATA1 signaling pathway

Continuous growth and development in plants are accomplished by meristems, groups of undifferentiated cells that persist as stem cells and initiate organs. While the structures of the apical and floral meristems in dicotyledonous plants have been well described, little is known about the underlying molecular mechanisms controlling cell proliferation and differentiation in these structures. We have shown previously that the CLAVATA1 (CLV1) gene in Arabidopsis encodes a receptor kinase-like protein that controls the size of the apical and floral meristems. Here, we show that KAPP, a gene encoding a kinase-associated protein phosphatase, is expressed in apical and young floral meristems, along with CLV1. Overexpression of KAPP mimics the clv1 mutant phenotype. Furthermore, CLV1 has kinase activity: it phosphorylates both itself and KAPP. Finally, KAPP binds and dephosphorylates CLV1. We present a model where KAPP functions as a negative regulator of the CLAVATA1 signal transduction pathway.

A structural census of the current population of protein sequences

We examine the occurrence of the ≈300 known protein folds in different groups of organisms. To do this, we characterize a large fraction of the currently known protein sequences (≈140,000) in structural terms, by matching them to known structures via sequence comparison (or by secondary-structure class prediction for those without structural homologues). Overall, we find that an appreciable fraction of the known folds are present in each of the major groups of organisms (e.g., bacteria and eukaryotes share 156 of 275 folds), and most of the common folds are associated with many families of nonhomologous sequences (i.e., >10 sequence families for each common fold). However, different groups of organisms have characteristically distinct distributions of folds. So, for instance, some of the most common folds in vertebrates, such as globins or zinc fingers, are rare or absent in bacteria. Many of these differences in fold usage are biologically reasonable, such as the folds of metabolic enzymes being common in bacteria and those associated with extracellular transport and communication being common in animals. They also have important implications for database-based methods for fold recognition, suggesting that an unknown sequence from a plant is more likely to have a certain fold (e.g., a TIM barrel) than an unknown sequence from an animal.