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.

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.

Neuregulin-3 (NRG3): A novel neural tissue-enriched protein that binds and activates ErbB4

We describe the identification of Neuregulin-3 (NRG3), a novel protein that is structurally related to the neuregulins (NRG1). The NRG1/neuregulins are a diverse family of proteins that arise by alternative splicing from a single gene. These proteins play an important role in controlling the growth and differentiation of glial, epithelial, and muscle cells. The biological effects of NRG1 are mediated by receptor tyrosine kinases ErbB2, ErbB3, and ErbB4. However, genetic studies have suggested that the activity of ErbB4 may also be regulated in the central nervous system by a ligand distinct from NRG1. NRG3 is predicted to contain an extracellular domain with an epidermal growth factor (EGF) motif, a transmembrane domain, and a large cytoplasmic domain. We show that the EGF-like domain of NRG3 binds to the extracellular domain of ErbB4 in vitro. Moreover, NRG3 binds to ErbB4 expressed on cells and stimulates tyrosine phosphorylation of this receptor. The expression of NRG3 is highly restricted to the developing and adult nervous system. These data suggest that NRG3 is a novel, neural-enriched ligand for ErbB4.

Activation of the c-Jun N-terminal kinase pathway by a novel protein kinase related to human germinal center kinase

The c-Jun N-terminal kinase (JNK), or stress-activated protein kinase plays a crucial role in cellular responses stimulated by environmental stress and proinflammatory cytokines. However, the mechanisms that lead to the activation of the JNK pathway have not been elucidated. We have isolated a cDNA encoding a novel protein kinase that has significant sequence similarities to human germinal center kinase (GCK) and human hematopoietic progenitor kinase 1. The novel GCK-like kinase (GLK) has a nucleotide sequence that encodes an ORF of 885 amino acids with 11 kinase subdomains. Endogenous GLK could be activated by UV radiation and proinflammatory cytokine tumor necrosis factor α. When transiently expressed in 293 cells, GLK specifically activated the JNK, but not the p42/44MAPK/extracellular signal-regulated kinase or p38 kinase signaling pathways. Interestingly, deletion of amino acids 353–835 in the putative C-terminal regulatory region, or mutation of Lys-35 in the putative ATP-binding domain, markedly reduced the ability of GLK to activate JNK. This result indicates that both kinase activity and the C-terminal region of GLK are required for maximal activation of JNK. Furthermore, GLK-induced JNK activation could be inhibited by a dominant-negative mutant of mitogen-activated protein kinase kinase kinase 1 (MEKK1) or mitogen-activated protein kinase kinase 4/SAPK/ERK kinase 1 (SEK1), suggesting that GLK may function upstream of MEKK1 in the JNK signaling pathway.

Functional characterization of the C—C chemokine-like molecules encoded by molluscum contagiosum virus types 1 and 2

Many viruses have evolved mechanisms for evading the host immune system by synthesizing proteins that interfere with the normal immune response. The poxviruses are among the most accomplished at deceiving their hosts’ immune systems. The nucleotide sequence of the genome of the human cutaneous poxvirus, molluscum contagiosum virus (MCV) type 1, was recently reported to contain a region that resembles a human chemokine. We have cloned and expressed the chemokine-like genes from MCV type 1 and the closely related MCV type 2 to determine a potential role for these proteins in the viral life cycle. In monocyte chemotaxis assays, the viral proteins have no chemotactic activity but both viral proteins block the chemotactic response to the human chemokine, macrophage inflammatory protein (MIP)-1α. Like MIP-1α, both viral proteins also inhibit the growth of human hematopoietic progenitor cells, but the viral proteins are more potent in this activity than the human chemokine. These viral chemokines antagonize the chemotactic activity of human chemokines and have an inhibitory effect on human hematopoietic progenitor cells. We hypothesize that the inhibition of chemotaxis is an immune evasion function of these proteins during molluscum contagiosum virus infection. The significance of hematopoietic progenitor cell inhibition in viral pathogenesis is uncertain.

Structural mimicry of a native protein by a minimized binding domain

The affinity between molecules depends both on the nature and presentation of the contacts. Here, we observe coupling of functional and structural elements when a protein binding domain is evolved to a smaller functional mimic. Previously, a 38-residue form of the 59-residue B-domain of protein A, termed Z38, was selected by phage display. Z38 contains 13 mutations and binds IgG only 10-fold weaker than the native B-domain. We present the solution structure of Z38 and show that it adopts a tertiary structure remarkably similar to that observed for the first two helices of B-domain in the B-domain/Fc complex [Deisenhofer, J. (1981) Biochemistry 20, 2361–2370], although it is significantly less stable. Based on this structure, we have improved on Z38 by designing a 34-residue disulfide-bonded variant (Z34C) that has dramatically enhanced stability and binds IgG with 9-fold higher affinity. The improved stability of Z34C led to NMR spectra with much greater chemical shift dispersion, resulting in a more precisely determined structure. Z34C, like Z38, has a structure virtually identical to the equivalent region from native protein A domains. The well-defined hydrophobic core of Z34C reveals key structural features that have evolved in this small, functional domain. Thus, the stabilized two-helix peptide, about half the size and having one-third of the remaining residues altered, accurately mimics both the structure and function of the native domain.

Probing the role of packing specificity in protein design

By using a protein-design algorithm that quantitatively considers side-chain packing, the effect of specific steric constraints on protein design was assessed in the core of the streptococcal protein G β1 domain. The strength of packing constraints used in the design was varied, resulting in core sequences that reflected differing amounts of packing specificity. The structural flexibility and stability of several of the designed proteins were experimentally determined and showed a trend from well-ordered to highly mobile structures as the degree of packing specificity in the design decreased. This trend both demonstrates that the inclusion of specific packing interactions is necessary for the design of native-like proteins and defines a useful range of packing specificity for the design algorithm. In addition, an analysis of the modeled protein structures suggested that penalizing for exposed hydrophobic surface area can improve design performance.

The macrophage/endothelial cell mannose receptor cDNA encodes a protein that binds oligosaccharides terminating with SO4-4-GalNAcβ1,4GlcNAcβ or Man at independent sites

Lutropin (LH) and other glycoproteins bearing oligosaccharides with the terminal sequence SO4-4-GalNAcβ1,4GlcNAcβ1,4Man- (S4GGnM) are rapidly removed from the circulation by an S4GGnM-specific receptor (S4GGnM-R) expressed at the surface of hepatic endothelial cells. The S4GGnM-R isolated from rat liver is closely related to the macrophage mannose-specific receptor (Man-R) isolated from rat lung both antigenically and structurally. The S4GGnM-R and Man-R isolated from these tissues nonetheless differ in their ability to bind ligands bearing terminal GalNAc-4-SO4 or Man. In this paper, we have explored the structural relationship between the Man-R and the S4GGnM-R by examining the properties of the recombinant Man-R in the form of a transmembrane protein and a soluble chimeric fusion protein in which the transmembrane and cytosolic domains have been replaced by the Fc region of human IgG1. Like the S4GGnM-R isolated from liver, the chimeric fusion protein is able to bind ligands terminating with GalNAc-4-SO4 and Man at independent sites. When expressed in CHO cells the recombinant Man-R is able to mediate the uptake of ligands bearing either terminal GalNAc-4-SO4 or terminal Man. We propose that the Man-R be renamed the Man/S4GGnM receptor on the basis of its multiple and independent specificities.

Two sterol regulatory element-like sequences mediate up-regulation of caveolin gene transcription in response to low density lipoprotein free cholesterol

Caveolae form the terminus for a major pathway of intracellular free cholesterol (FC) transport. Caveolin mRNA levels in confluent human skin fibroblasts were up-regulated following increased uptake of low density lipoprotein (LDL) FC. The increase induced by FC was not associated with detectable change in mRNA stability, indicating that caveolin mRNA levels were mediated at the level of gene transcription. A total of 924 bp of 5′ flanking region of the caveolin gene were cloned and sequenced. The promoter sequence included three G+C-rich potential sterol regulatory elements (SREs), a CAAT sequence and a Sp1 consensus sequence. Deletional mutagenesis of individual SRE-like sequences indicated that of these two (at −646 and −395 bp) were essential for the increased transcription rates mediated by LDL-FC, whereas the third was inconsequential. Gel shift analysis of protein binding from nuclear extracts to these caveolin promoter DNA sequences, together with DNase I footprinting, confirmed nucleoprotein binding to the SRE-like elements as part of the transcriptional response to LDL-FC. A supershift obtained with antibody to SRE-binding protein 1 (SPEBP-1) indicated that this protein binds at −395 bp. There was no reaction at −395 bp with anti-Sp1 antibody nor with either antibody at −646 bp. The cysteine protease inhibitor N-acetyl-leu-leu-norleucinal (ALLN), which inhibits SREBP catabolism, superinhibited caveolin mRNA levels regardless of LDL-FC. This finding suggests that SREBP inhibits caveolin gene transcription in contrast to its stimulating effect on other promoters. The findings of this study are consistent with the postulated role for caveolin as a regulator of cellular FC homeostasis in quiescent peripheral cells, and the coordinate regulation by SREBP of FC influx and efflux.

CLARP, a death effector domain-containing protein interacts with caspase-8 and regulates apoptosis

We have identified and characterized CLARP, a caspase-like apoptosis-regulatory protein. Sequence analysis revealed that human CLARP contains two amino-terminal death effector domains fused to a carboxyl-terminal caspase-like domain. The structure and amino acid sequence of CLARP resemble those of caspase-8, caspase-10, and DCP2, a Drosophila melanogaster protein identified in this study. Unlike caspase-8, caspase-10, and DCP2, however, two important residues predicted to be involved in catalysis were lost in the caspase-like domain of CLARP. Analysis with fluorogenic substrates for caspase activity confirmed that CLARP is catalytically inactive. CLARP was found to interact with caspase-8 but not with FADD/MORT-1, an upstream death effector domain-containing protein of the Fas and tumor necrosis factor receptor 1 signaling pathway. Expression of CLARP induced apoptosis, which was blocked by the viral caspase inhibitor p35, dominant negative mutant caspase-8, and the synthetic caspase inhibitor benzyloxycarbonyl-Val-Ala-Asp-(OMe)-fluoromethylketone (zVAD-fmk). Moreover, CLARP augmented the killing ability of caspase-8 and FADD/MORT-1 in mammalian cells. The human clarp gene maps to 2q33. Thus, CLARP represents a regulator of the upstream caspase-8, which may play a role in apoptosis during tissue development and homeostasis.

Genetic definition of a protein-splicing domain: Functional mini-inteins support structure predictions and a model for intein evolution

Inteins are protein-splicing elements, most of which contain conserved sequence blocks that define a family of homing endonucleases. Like group I introns that encode such endonucleases, inteins are mobile genetic elements. Recent crystallography and computer modeling studies suggest that inteins consist of two structural domains that correspond to the endonuclease and the protein-splicing elements. To determine whether the bipartite structure of inteins is mirrored by the functional independence of the protein-splicing domain, the entire endonuclease component was deleted from the Mycobacterium tuberculosis recA intein. Guided by computer modeling studies, and taking advantage of genetic systems designed to monitor intein function, the 440-aa Mtu recA intein was reduced to a functional mini-intein of 137 aa. The accuracy of splicing of several mini-inteins was verified. This work not only substantiates structure predictions for intein function but also supports the hypothesis that, like group I introns, mobile inteins arose by an endonuclease gene invading a sequence encoding a small, functional splicing element.

Sterols regulate processing of carbohydrate chains of wild-type SREBP cleavage-activating protein (SCAP), but not sterol-resistant mutants Y298C or D443N

SREBP cleavage activating protein (SCAP), a membrane-bound glycoprotein, regulates the proteolytic activation of sterol regulatory element binding proteins (SREBPs), which are membrane-bound transcription factors that control lipid synthesis in animal cells. SCAP-stimulated proteolysis releases active fragments of SREBPs from membranes of the endoplasmic reticulum and allows them to enter the nucleus where they activate transcription. Sterols such as 25-hydroxycholesterol inactivate SCAP, suppressing SREBP proteolysis and turning off cholesterol synthesis. We here report the isolation of Chinese hamster ovary cells with a point mutation in SCAP (Y298C) that renders the protein resistant to inhibition by 25-hydroxycholesterol. Like the previously described D443N mutation, the Y298C mutation occurs within the putative sterol-sensing domain, which is part of the polytopic membrane attachment region of SCAP. Cells that express SCAP(Y298C) continued to process SREBPs in the presence of 25-hydroxycholesterol and hence they resisted killing by this sterol. In wild-type Chinese hamster ovary cells the N-linked carbohydrate chains of SCAP were mostly in the endoglycosidase H-sensitive form when cells were grown in medium containing 25-hydroxycholesterol. In contrast, when cells were grown in sterol-depleted medium, these chains were converted to an endoglycosidase H-resistant form. 25-Hydroxycholesterol had virtually no effect in cells expressing SCAP(D443N) or SCAP(Y298C). The relation between this regulated carbohydrate processing to the SCAP-regulated proteolysis of SREBP remains to be explored.

Structure of the thrombin complex with triabin, a lipocalin-like exosite-binding inhibitor derived from a triatomine bug

Triabin, a 142-residue protein from the saliva of the blood-sucking triatomine bug Triatoma pallidipennis, is a potent and selective thrombin inhibitor. Its stoichiometric complex with bovine α-thrombin was crystallized, and its crystal structure was solved by Patterson search methods and refined at 2.6-Å resolution to an R value of 0.184. The analysis revealed that triabin is a compact one-domain molecule essentially consisting of an eight-stranded β-barrel. The eight strands A to H are arranged in the order A-C-B-D-E-F-G-H, with the first four strands exhibiting a hitherto unobserved up-up-down-down topology. Except for the B-C inversion, the triabin fold exhibits the regular up-and-down topology of lipocalins. In contrast to the typical ligand-binding lipocalins, however, the triabin barrel encloses a hydrophobic core intersected by a unique salt-bridge cluster. Triabin interacts with thrombin exclusively via its fibrinogen-recognition exosite. Surprisingly, most of the interface interactions are hydrophobic. A prominent exception represents thrombin’s Arg-77A side chain, which extends into a hydrophobic triabin pocket forming partially buried salt bridges with Glu-128 and Asp-135 of the inhibitor. The fully accessible active site of thrombin in this complex is in agreement with its retained hydrolytic activity toward small chromogenic substrates. Impairment of thrombin’s fibrinogen converting activity or of its thrombomodulin-mediated protein C activation capacity upon triabin binding is explained by usage of overlapping interaction sites of fibrinogen, thrombomodulin, and triabin on thrombin. These data demonstrate that triabin inhibits thrombin via a novel and unique mechanism that might be of interest in the context of potential therapeutic applications.

The helical domain of a G protein α subunit is a regulator of its effector

The α subunit (Gα) of heterotrimeric G proteins is a major determinant of signaling selectivity. The Gα structure essentially comprises a GTPase “Ras-like” domain (RasD) and a unique α-helical domain (HD). We used the vertebrate phototransduction model to test for potential functions of HD and found that the HD of the retinal transducin Gα (Gαt) and the closely related gustducin (Gαg), but not Gαi1, Gαs, or Gαq synergistically enhance guanosine 5′-γ[-thio]triphosphate bound Gαt (GαtGTPγS) activation of bovine rod cGMP phosphodiesterase (PDE). In addition, both HDt and HDg, but not HDi1, HDs, or HDq attenuate the trypsin-activated PDE. GαtGDP and HDt attenuation of trypsin-activated PDE saturate with similar affinities and to an identical 38% of initial activity. These data suggest that interaction of intact Gαt with the PDE catalytic core may be caused by the HD moiety, and they indicate an independent site(s) for the HD moiety of Gαt within the PDE catalytic core in addition to the sites for the inhibitory Pγ subunits. The HD moiety of GαtGDP is an attenuator of the activated catalytic core, whereas in the presence of activated GαtGTPγS the independently expressed HDt is a potent synergist. Rhodopsin catalysis of Gαt activation enhances the PDE activation produced by subsaturating levels of Gαt, suggesting a HD-moiety synergism from a transient conformation of Gαt. These results establish HD-selective regulations of vertebrate retinal PDE, and they provide evidence demonstrating that the HD is a modulatory domain. We suggest that the HD works in concert with the RasD, enhancing the efficiency of G protein signaling.