A yeast genetic system for selecting small molecule inhibitors of protein–protein interactions in nanodroplets

Cellular processes are mediated by complex networks of molecular interactions. Dissection of their role most commonly is achieved by using genetic mutations that alter, for example, protein–protein interactions. Small molecules that accomplish the same result would provide a powerful complement to the genetic approach, but it generally is believed that such molecules are rare. There are several natural products, however, that illustrate the feasibility of this approach. Split-pool synthesis now provides a simple mechanical means to prepare vast numbers of complex, even natural product-like, molecules individually attached to cell-sized polymer beads. Here, we describe a genetic system compatible with split-pool synthesis that allows the detection of cell-permeable, small molecule inhibitors of protein–protein interactions in 100- to 200-nl cell culture droplets, prepared by a recently described technique that arrays large numbers of such droplets. These “nanodroplets” contain defined media, cells, and one or more beads containing ≈100 pmol of a photoreleasable small molecule and a controlled number of cells. The engineered Saccharomyces cerevisiae cells used in this study express two interacting proteins after induction with galactose whose interaction results in cell death in the presence of 5-fluoroorotic acid (inducible reverse two-hybrid assay). Disruption of the interaction by a small molecule allows growth, and the small molecule can be introduced into the system hours before induction of the toxic interaction. We demonstrate that the interaction between the activin receptor R1 and the immunophilin protein FKBP12 can be disrupted by the small molecule FK506 at nanomolar concentrations in nanodroplets. This system should provide a general method for selecting cell-permeable ligands that can be used to study the relevance of protein–protein interactions in living cells or organisms.

Identification of a second aryl phosphate-binding site in protein-tyrosine phosphatase 1B: A paradigm for inhibitor design

The structure of the catalytically inactive mutant (C215S) of the human protein-tyrosine phosphatase 1B (PTP1B) has been solved to high resolution in two complexes. In the first, crystals were grown in the presence of bis-(para-phosphophenyl) methane (BPPM), a synthetic high-affinity low-molecular weight nonpeptidic substrate (Km = 16 μM), and the structure was refined to an R-factor of 18.2% at 1.9 Å resolution. In the second, crystals were grown in a saturating concentration of phosphotyrosine (pTyr), and the structure was refined to an R-factor of 18.1% at 1.85 Å. Difference Fourier maps showed that BPPM binds PTP1B in two mutually exclusive modes, one in which it occupies the canonical pTyr-binding site (the active site), and another in which a phosphophenyl moiety interacts with a set of residues not previously observed to bind aryl phosphates. The identification of a second pTyr molecule at the same site in the PTP1B/C215S–pTyr complex confirms that these residues constitute a low-affinity noncatalytic aryl phosphate-binding site. Identification of a second aryl phosphate binding site adjacent to the active site provides a paradigm for the design of tight-binding, highly specific PTP1B inhibitors that can span both the active site and the adjacent noncatalytic site. This design can be achieved by tethering together two small ligands that are individually targeted to the active site and the proximal noncatalytic site.

Computer-based analysis of the binding steps in protein complex formation

Computer models were used to examine whether and under what conditions the multimeric protein complex is inhibited by high concentrations of one of its components—an effect analogous to the prozone phenomenon in precipitin tests. A series of idealized simple “ball-and-stick” structures representing small oligomeric complexes of protein molecules formed by reversible binding reactions were analyzed to determine the binding steps leading to each structure. The equilibrium state of each system was then determined over a range of starting concentrations and Kds and the steady-state concentration of structurally complete oligomer calculated for each situation. A strong inhibitory effect at high concentrations was shown by any protein molecule forming a bridge between two or more separable parts of the complex. By contrast, proteins linked to the outside of the complex by a single bond showed no inhibition whatsoever at any concentration. Nonbridging, multivalent proteins in the body of the complex could show an inhibitory effect or not depending on the structure of the complex and the strength of its bonds. On the basis of this study, we suggest that the prozone phenomenon will occur widely in living cells and that it could be a crucial factor in the regulation of protein complex formation.

Functional separation of pre-rRNA processing steps revealed by truncation of the U3 small nucleolar ribonucleoprotein component, Mpp10

The U3 small nucleolar ribonucleoprotein (snoRNP) is required for three cleavage events that generate the mature 18S rRNA from the pre-rRNA. In Saccharomyces cerevisiae, depletion of Mpp10, a U3 snoRNP-specific protein, halts 18S rRNA production and impairs cleavage at the three U3 snoRNP-dependent sites: A0, A1, and A2. We have identified truncation mutations of Mpp10 that affect 18S rRNA synthesis and confer cold-sensitivity and slow growth. However, distinct from yeast cells depleted of Mpp10, the mutants carrying these truncated Mpp10 proteins accumulate a novel precursor, resulting from cleavage at only A0. The Mpp10 truncations do not alter association of Mpp10 with the U3 snoRNA, nor do they affect snoRNA or protein stability. Thus, the role in processing of the U3 snoRNP can be separated into cleavage at the A0 site, which occurs in the presence of truncated Mpp10, and cleavage at the A1/A2 sites, which occurs only with intact Mpp10. These results strongly argue for a role for Mpp10 in processing at the A1/A2 sites.

Mode matches and their locations in the hydrophobic free energy sequences of peptide ligands and their receptor eigenfunctions

Patterns in sequences of amino acid hydrophobic free energies predict secondary structures in proteins. In protein folding, matches in hydrophobic free energy statistical wavelengths appear to contribute to selective aggregation of secondary structures in “hydrophobic zippers.” In a similar setting, the use of Fourier analysis to characterize the dominant statistical wavelengths of peptide ligands’ and receptor proteins’ hydrophobic modes to predict such matches has been limited by the aliasing and end effects of short peptide lengths, as well as the broad-band, mode multiplicity of many of their frequency (power) spectra. In addition, the sequence locations of the matching modes are lost in this transformation. We make new use of three techniques to address these difficulties: (i) eigenfunction construction from the linear decomposition of the lagged covariance matrices of the ligands and receptors as hydrophobic free energy sequences; (ii) maximum entropy, complex poles power spectra, which select the dominant modes of the hydrophobic free energy sequences or their eigenfunctions; and (iii) discrete, best bases, trigonometric wavelet transformations, which confirm the dominant spectral frequencies of the eigenfunctions and locate them as (absolute valued) moduli in the peptide or receptor sequence. The leading eigenfunction of the covariance matrix of a transmembrane receptor sequence locates the same transmembrane segments seen in n-block-averaged hydropathy plots while leaving the remaining hydrophobic modes unsmoothed and available for further analyses as secondary eigenfunctions. In these receptor eigenfunctions, we find a set of statistical wavelength matches between peptide ligands and their G-protein and tyrosine kinase coupled receptors, ranging across examples from 13.10 amino acids in acid fibroblast growth factor to 2.18 residues in corticotropin releasing factor. We find that the wavelet-located receptor modes in the extracellular loops are compatible with studies of receptor chimeric exchanges and point mutations. A nonbinding corticotropin-releasing factor receptor mutant is shown to have lost the signatory mode common to the normal receptor and its ligand. Hydrophobic free energy eigenfunctions and their transformations offer new quantitative physical homologies in database searches for peptide-receptor matches.

Production, specificity, and functionality of monoclonal antibodies to specific peptide–major histocompatibility complex class II complexes formed by processing of exogenous protein

Several unanswered questions in T cell immunobiology relating to intracellular processing or in vivo antigen presentation could be approached if convenient, specific, and sensitive reagents were available for detecting the peptide–major histocompatibility complex (MHC) class I or class II ligands recognized by αβ T cell receptors. For this reason, we have developed a method using homogeneously loaded peptide–MHC class II complexes to generate and select specific mAb reactive with these structures using hen egg lysozyme (HEL) and I-Ak as a model system. mAbs specific for either HEL-(46–61)–Ak or HEL-(116–129)–Ak have been isolated. They cross-react with a small subset of I-Ak molecules loaded with self peptides but can nonetheless be used for flow cytometry, immunoprecipitation, Western blotting, and intracellular immunofluorescence to detect specific HEL peptide–MHC class II complexes formed by either peptide exposure or natural processing of native HEL. An example of the utility of these reagents is provided herein by using one of the anti-HEL-(46–61)–Ak specific mAbs to visualize intracellular compartments where I-Ak is loaded with HEL-derived peptides early after antigen administration. Other uses, especially for in vivo tracking of specific ligand-bearing antigen-presenting cells, are discussed.

A FYVE-finger-containing protein, Rabip4, is a Rab4 effector involved in early endosomal traffic

The small GTPase Rab4 is implicated in endocytosis in all cell types, but also plays a specific role in some regulated processes. To better understand the role of Rab4 in regulation of vesicular trafficking, we searched for an effector(s) that specifically recognizes its GTP-bound form. We cloned a ubiquitous 69-kDa protein, Rabip4, that behaves as a Rab4 effector in the yeast two-hybrid system and in the mammalian cell. Rabip4 contains two coiled-coil domains and a FYVE-finger domain. When expressed in CHO cells, Rabip4 is present in early endosomes, because it is colocated with endogenous Early Endosome Antigen 1, although it is absent from Rab11-positive recycling endosomes and Rab-7 positive late endosomes. The coexpression of Rabip4 with active Rab4, but not with inactive Rab4, leads to an enlargement of early endosomes. It strongly increases the degree of colocalization of markers of sorting (Rab5) and recycling (Rab11) endosomes with Rab4. Furthermore, the expression of Rabip4 leads to the intracellular retention of a recycling molecule, the glucose transporter Glut 1. We propose that Rabip4, an effector of Rab4, controls early endosomal traffic possibly by activating a backward transport step from recycling to sorting endosomes.

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.

Molecular cloning of Ian4: a BCR/ABL-induced gene that encodes an outer membrane mitochondrial protein with GTP-binding activity

Using the representation difference analysis technique, we have identified a novel gene, Ian4, which is preferentially expressed in hematopoietic precursor 32D cells transfected with wild-type versus mutant forms of the Bcr/Abl oncogene. Ian4 expression was undetectable in 32D cells transfected with v-src, oncogenic Ha-ras or v-Abl. Murine Ian4 maps to chromosome 6, 25 cM from the centromere. The Ian4 mRNA contains two open reading frames (ORFs) separated by 5 nt. The first ORF has the potential to encode for a polypeptide of 67 amino acids without apparent homology to known proteins. The second ORF encodes a protein of 301 amino acids with a GTP/ATP-binding site in the N-terminus and a hydrophobic domain in the extreme C-terminus. The IAN-4 protein resides in the mitochondrial outer membrane and the last 20 amino acids are necessary for this localization. The IAN-4 protein has GTP-binding activity and shares sequence homology with a novel family of putative GTP-binding proteins: the immuno-associated nucleotide (IAN) family.

Supernatant protein factor, which stimulates the conversion of squalene to lanosterol, is a cytosolic squalene transfer protein and enhances cholesterol biosynthesis

Squalene epoxidase, a membrane-associated enzyme that converts squalene to squalene 2,3-oxide, plays an important role in the maintenance of cholesterol homeostasis. In 1957, Bloch and colleagues identified a factor from rat liver cytosol termed “supernatant protein factor (SPF),” which promotes the squalene epoxidation catalyzed by rat liver microsomes with oxygen, NADPH, FAD, and phospholipid [Tchen, T. T. & Bloch, K. (1957) J. Biol. Chem. 226, 921–930]. Although purification of SPF by 11,000-fold was reported, no information is so far available on the primary structure or biological function of SPF. Here we report the cDNA cloning and expression of SPF from rat and human. The encoded protein of 403 amino acids belongs to a family of cytosolic lipid-binding/transfer proteins such as α-tocopherol transfer protein, cellular retinal binding protein, yeast phosphatidylinositol transfer protein (Sec14p), and squid retinal binding protein. Recombinant SPF produced in Escherichia coli enhances microsomal squalene epoxidase activity and promotes intermembrane transfer of squalene in vitro. SPF mRNA is expressed abundantly in the liver and small intestine, both of which are important sites of cholesterol biosynthesis. SPF is expressed significantly in isolated hepatocytes, but the expression level was markedly decreased after 48 h of in vitro culture. Moreover, SPF was not detectable in most of the cell lines tested, including HepG2 and McARH7777 hepatomas. Transfection of SPF cDNA in McARH7777 significantly stimulated de novo cholesterol biosynthesis. These data suggest that SPF is a cytosolic squalene transfer protein capable of regulating cholesterol biosynthesis.

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.

Synechocystis HSP17 is an amphitropic protein that stabilizes heat-stressed membranes and binds denatured proteins for subsequent chaperone-mediated refolding

The small heat shock proteins (sHSPs) are ubiquitous stress proteins proposed to act as molecular chaperones to prevent irreversible protein denaturation. We characterized the chaperone activity of Synechocystis HSP17 and found that it has not only protein-protective activity, but also a previously unrecognized ability to stabilize lipid membranes. Like other sHSPs, recombinant Synechocystis HSP17 formed stable complexes with denatured malate dehydrogenase and served as a reservoir for the unfolded substrate, transferring it to the DnaK/DnaJ/GrpE and GroEL/ES chaperone network for subsequent refolding. Large unilamellar vesicles made of synthetic and cyanobacterial lipids were found to modulate this refolding process. Investigation of HSP17-lipid interactions revealed a preference for the liquid crystalline phase and resulted in an elevated physical order in model lipid membranes. Direct evidence for the participation of HSP17 in the control of thylakoid membrane physical state in vivo was gained by examining an hsp17− deletion mutant compared with the isogenic wild-type hsp17+ revertant Synechocystis cells. We suggest that, together with GroEL, HSP17 behaves as an amphitropic protein and plays a dual role. Depending on its membrane or cytosolic location, it may function as a “membrane stabilizing factor” as well as a member of a multichaperone protein-folding network. Membrane association of sHSPs could antagonize the heat-induced hyperfluidization of specific membrane domains and thereby serve to preserve structural and functional integrity of biomembranes.

The Cellular Level of PR500, a Protein Complex Related to the 19S Regulatory Particle of the Proteasome, Is Regulated in Response to Stresses in Plants

In Arabidopsis seedlings and cauliflower florets, Rpn6 (a proteasome non-ATPase regulatory subunit) was found in two distinct protein complexes of ∼800 and 500 kDa, respectively. The large complex likely represents the proteasome 19S regulator particle (RP) because it displays the expected subunit composition and all characteristics. The small complex, designated PR500, shares at least three subunits with the “lid” subcomplex of 19S RP and is loosely associated with an hsp70 protein. In Arabidopsis COP9 signalosome mutants, PR500 was specifically absent or reduced to an extent that correlates with the severity of the mutations. Furthermore, PR500 was also diminished in response to potential protein-misfolding stresses caused by the heat shock and canavanine treatment. Immunofluorescence studies suggest that PR500 has a distinct localization pattern and is enriched in specific nuclear foci. We propose that PR500 may be evolved in higher plants to cope with the frequently encountered environmental stresses.

A Family of Small Coiled-Coil–forming Proteins Functioning at the Late Endosome in Yeast

The multispanning membrane protein Ste6, a member of the ABC-transporter family, is transported to the yeast vacuole for degradation. To identify functions involved in the intracellular trafficking of polytopic membrane proteins, we looked for functions that block Ste6 transport to the vacuole upon overproduction. In our screen, we identified several known vacuolar protein sorting (VPS) genes (SNF7/VPS32, VPS4, and VPS35) and a previously uncharacterized open reading frame, which we named MOS10 (more of Ste6). Sequence analysis showed that Mos10 is a member of a small family of coiled-coil–forming proteins, which includes Snf7 and Vps20. Deletion mutants of all three genes stabilize Ste6 and show a “class E vps phenotype.” Maturation of the vacuolar hydrolase carboxypeptidase Y was affected in the mutants and the endocytic tracer FM4-64 and Ste6 accumulated in a dot or ring-like structure next to the vacuole. Differential centrifugation experiments demonstrated that about half of the hydrophilic proteins Mos10 and Vps20 was membrane associated. The intracellular distribution was further analyzed for Mos10. On sucrose gradients, membrane-associated Mos10 cofractionated with the endosomal t-SNARE Pep12, pointing to an endosomal localization of Mos10. The growth phenotypes of the mutants suggest that the “Snf7-family” members are involved in a cargo-specific event.

Computational method to reduce the search space for directed protein evolution

We introduce a computational method to optimize the in vitro evolution of proteins. Simulating evolution with a simple model that statistically describes the fitness landscape, we find that beneficial mutations tend to occur at amino acid positions that are tolerant to substitutions, in the limit of small libraries and low mutation rates. We transform this observation into a design strategy by applying mean-field theory to a structure-based computational model to calculate each residue's structural tolerance. Thermostabilizing and activity-increasing mutations accumulated during the experimental directed evolution of subtilisin E and T4 lysozyme are strongly directed to sites identified by using this computational approach. This method can be used to predict positions where mutations are likely to lead to improvement of specific protein properties.