A computationally directed screen identifying interacting coiled coils from Saccharomyces cerevisiae

Computational methods can frequently identify protein-interaction motifs in otherwise uncharacterized open reading frames. However, the identification of candidate ligands for these motifs (e.g., so that partnering can be determined experimentally in a directed manner) is often beyond the scope of current computational capabilities. One exception is provided by the coiled-coil interaction motif, which consists of two or more α helices that wrap around each other: the ligands for coiled-coil sequences are generally other coiled-coil sequences, thereby greatly simplifying the motif/ligand recognition problem. Here, we describe a two-step approach to identifying protein–protein interactions mediated by two-stranded coiled coils that occur in Saccharomyces cerevisiae. Coiled coils from the yeast genome are first predicted computationally, by using the multicoil program, and associations between coiled coils are then determined experimentally by using the yeast two-hybrid assay. We report 213 unique interactions between 162 putative coiled-coil sequences. We evaluate the resulting interactions, focusing on associations identified between components of the spindle pole body (the yeast centrosome).

High-sensitivity array analysis of gene expression for the early detection of disseminated breast tumor cells in peripheral blood

Early detection is an effective means of reducing cancer mortality. Here, we describe a highly sensitive high-throughput screen that can identify panels of markers for the early detection of solid tumor cells disseminated in peripheral blood. The method is a two-step combination of differential display and high-sensitivity cDNA arrays. In a primary screen, differential display identified 170 candidate marker genes differentially expressed between breast tumor cells and normal breast epithelial cells. In a secondary screen, high-sensitivity arrays assessed expression levels of these genes in 48 blood samples, 22 from healthy volunteers and 26 from breast cancer patients. Cluster analysis identified a group of 12 genes that were elevated in the blood of cancer patients. Permutation analysis of individual genes defined five core genes (P ≤ 0.05, permax test). As a group, the 12 genes generally distinguished accurately between healthy volunteers and patients with breast cancer. Mean expression levels of the 12 genes were elevated in 77% (10 of 13) untreated invasive cancer patients, whereas cluster analysis correctly classified volunteers and patients (P = 0.0022, Fisher's exact test). Quantitative real-time PCR confirmed array results and indicated that the sensitivity of the assay (1:2 × 108 transcripts) was sufficient to detect disseminated solid tumor cells in blood. Expression-based blood assays developed with the screening approach described here have the potential to detect and classify solid tumor cells originating from virtually any primary site in the body.

Dynamic properties of the Ras switch I region and its importance for binding to effectors

We have investigated the dynamic properties of the switch I region of the GTP-binding protein Ras by using mutants of Thr-35, an invariant residue necessary for the switch function. Here we show that these mutants, previously used as partial loss-of-function mutations in cell-based assays, have a reduced affinity to Ras effector proteins without Thr-35 being involved in any interaction. The structure of Ras(T35S)⋅GppNHp was determined by x-ray crystallography. Whereas the overall structure is very similar to wildtype, residues from switch I are completely invisible, indicating that the effector loop region is highly mobile. 31P-NMR data had indicated an equilibrium between two rapidly interconverting conformations, one of which (state 2) corresponds to the structure found in the complex with the effectors. 31P-NMR spectra of Ras mutants (T35S) and (T35A) in the GppNHp form show that the equilibrium is shifted such that they occur predominantly in the nonbinding conformation (state 1). On addition of Ras effectors, Ras(T35S) but not Ras(T35A) shift to positions corresponding to the binding conformation. The structural data were correlated with kinetic experiments that show two-step binding reaction of wild-type and (T35S)Ras with effectors requires the existence of a rate-limiting isomerization step, which is not observed with T35A. The results indicate that minor changes in the switch region, such as removing the side chain methyl group of Thr-35, drastically affect dynamic behavior and, in turn, interaction with effectors. The dynamics of the switch I region appear to be responsible for the conservation of this threonine residue in GTP-binding proteins.

Cyanogenesis in Cassava1 : The Role of Hydroxynitrile Lyase in Root Cyanide Production

In the cyanogenic crop cassava (Manihot esculenta, Crantz), the final step in cyanide production is the conversion of acetone cyanohydrin, the deglycosylation product of linamarin, to cyanide plus acetone. This process occurs spontaneously at pH greater than 5.0 or enzymatically and is catalyzed by hydroxynitrile lyase (HNL). Recently, it has been demonstrated that acetone cyanohydrin is present in poorly processed cassava root food products. Since it has generally been assumed that HNL is present in all cassava tissues, we reinvestigated the enzymatic properties and tissue-specific distribution of HNL in cassava. We report the development of a rapid two-step purification protocol for cassava HNL, which yields an enzyme that is catalytically more efficient than previously reported (Hughes, J., Carvalho, F., and Hughes, M. [1994] Arch Biochem Biophys 311: 496–502). Analyses of the distribution of HNL activity and protein indicate that the accumulation of acetone cyanohydrin in roots is due to the absence of HNL, not to inhibition of the enzyme. Furthermore, the absence of HNL in roots and stems is associated with very low steady-state HNL transcript levels. It is proposed that the lack of HNL in cassava roots accounts for the high acetone cyanohydrin levels in poorly processed cassava food products.

Rapid activation of Na+/H+ exchange by aldosterone in renal epithelial cells requires Ca2+ and stimulation of a plasma membrane proton conductance.

There is increasing evidence for an additional acute, nongenomic action of the mineralocorticoid hormone aldosterone on renal epithelial cells, leading to a two-step model of mineralocorticoid action on electrolyte excretion. We investigated the acute effect of aldosterone on intracellular free Ca2+ and on intracellular pH in an aldosterone-sensitive Madin-Darby canine kidney cell clone. Within seconds of application of aldosterone, but not of the glucocorticoid hydrocortisone, there was a 3-fold sustained increase of intracellular Ca2+ at a half-maximal concentration of 10(-10) mol/liter. Omission of extracellular Ca2+ prevented this hormone response. In the presence of extracellular Ca2+ aldosterone led to intracellular alkalinization. The Na+/H+ exchange inhibitor ethyl-isopropanol-amiloride (EIPA) prevented the aldosterone-induced alkalinization but not the aldosterone-induced increase of intracellular Ca2+. Omission of extracellular Ca2+ also prevented aldosterone-induced alkalinization. Instead, aldosterone led to a Zn(2+)-dependent intracellular acidification in the presence of EIPA, indicative of an increase of plasma membrane proton conductance. Under control conditions, Zn2+ prevented the aldosterone-induced alkalinization completely. We conclude that aldosterone stimulated net-entry of Ca2+ from the extracellular compartment and a plasma membrane H+ conductance as prerequisites for the stimulation of plasma membrane Na+/H+ exchange which in turn modulates K+ channel acitivity. It is probable that the aldosterone-sensitive H+ conductance maintains Na+/H+ exchange activity by providing an acidic environment in the vicinity of the exchanger. Thus, genomic action of aldosterone determines cellular transport equipment, whereas the nongenomic action regulates transporter activity that requires responses within seconds or minutes, which explains the rapid effects on electrolyte excretion.

Molecular cloning of violaxanthin de-epoxidase from romaine lettuce and expression in Escherichia coli.

Plants need to avoid or dissipate excess light energy to protect photosystem II (PSII) from photoinhibitory damage. Higher plants have a conserved system that dissipates excess energy as heat in the light-harvesting complexes of PSII that depends on the transthylakoid delta pH and violaxanthin de-epoxidase (VDE) activity. To our knowledge, we report the first cloning of a cDNA encoding VDE and expression of functional enzyme in Escherichia coli. VDE is nuclear encoded and has a transit peptide with characteristic features of other lumen-localized proteins. The cDNA encodes a putative polypeptide of 473 aa with a calculated molecular mass of 54,447 Da. Cleavage of the transit peptide results in a mature putative polypeptide of 348 aa with a calculated molecular mass of 39,929 Da, close to the apparent mass of the purified enzyme (43 kDa). The protein has three interesting domains including (i) a cysteine-rich region, (ii) a lipocalin signature, and (iii) a highly charged region. The E. coli expressed enzyme de-epoxidizes violaxanthin sequentially to antheraxanthin and zeaxanthin, and is inhibited by dithiothreitol, similar to VDE purified from chloroplasts. This confirms that the cDNA encodes an authentic VDE of a higher plant and is unequivocal evidence that the same enzyme catalyzes the two-step mono de-epoxidation reaction. The cloning of VDE opens new opportunities for examining the function and evolution of the xanthophyll cycle, and possibly enhancing light-stress tolerance of plants.

A loop-loop "kissing" complex is the essential part of the dimer linkage of genomic HIV-1 RNA.

RNA-RNA interactions govern a number of biological processes. Several RNAs, including natural sense and antisense RNAs, interact by means of a two-step mechanism: recognition is mediated by a loop-loop complex, which is then stabilized by formation of an extended intermolecular duplex. It was proposed that the same mechanism holds for dimerization of the genomic RNA of human immunodeficiency virus type 1 (HIV-1), an event thought to control crucial steps of HIV-1 replication. However, whereas interaction between the partially self-complementary loop of the dimerization initiation site (DIS) of each monomer is well established, formation of the extended duplex remained speculative. Here we first show that in vitro dimerization of HIV-1 RNA is a specific process, not resulting from simple annealing of denatured molecules. Next we used mutants of the DIS to test the formation of the extended duplex. Four pairs of transcomplementary mutants were designed in such a way that all pairs can form the loop-loop "kissing" complex, but only two of them can potentially form the extended duplex. All pairs of mutants form heterodimers whose thermal stability, dissociation constant, and dynamics were analyzed. Taken together, our results indicate that, in contrast with the interactions between natural sense and antisense RNAs, no extended duplex is formed during dimerization of HIV-1 RNA. We also showed that 55-mer sense RNAs containing the DIS are able to interfere with the preformed HIV-1 RNA dimer.

Isolation of HLA-DR1.(staphylococcal enterotoxin A)2 trimers in solution.

Mutational studies indicate that the superantigen staphylococcal enterotoxin A (SEA) has two separate binding sites for major histocompatibility complex (MHC) class II molecules. Direct evidence is provided here for the formation of SEA-MHC class II trimers in solution. Isoelectric focusing separated SEA-HLA-DR1 complexes into both dimers and HLA-DR1.SEA2 trimers. The molar ratio of components was determined by dual isotope labeling. The SEA mutant SEA-F47S, L48S, Y92A, which is deficient in MHC class II alpha-chain binding, formed only dimers with HLA-DR1, whereas a second SEA mutant, SEA-H225A, which lacks high-affinity MHC class II beta-chain binding was incapable of forming any complexes. Thus SEA binding to its MHC receptor is a two-step process involving initial beta-chain binding followed by cooperative binding of a second SEA molecule to the class II alpha chain.

An amino acid sequence motif sufficient for subnuclear localization of an arginine/serine-rich splicing factor.

We have identified an amino acid sequence in the Drosophila Transformer (Tra) protein that is capable of directing a heterologous protein to nuclear speckles, regions of the nucleus previously shown to contain high concentrations of spliceosomal small nuclear RNAs and splicing factors. This sequence contains a nucleoplasmin-like bipartite nuclear localization signal (NLS) and a repeating arginine/serine (RS) dipeptide sequence adjacent to a short stretch of basic amino acids. Sequence comparisons from a number of other splicing factors that colocalize to nuclear speckles reveal the presence of one or more copies of this motif. We propose a two-step subnuclear localization mechanism for splicing factors. The first step is transport across the nuclear envelope via the nucleoplasmin-like NLS, while the second step is association with components in the speckled domain via the RS dipeptide sequence.

Microsatellite variability and genetic distances.

We analyze the within- and between-population dynamics of the distribution of the number of repeats at multiple microsatellite DNA loci subject to stepwise mutation. Analytical expressions for moments up to the fourth order within a locus and the variance of between-locus variance at mutation-drift equilibrium have been obtained. These statistics may be used to test the appropriateness of the one-step mutation model and to detect between-locus variation in the mutation rate. Published data are compatible with the one-step mutation model, although they do not reject the two-step model. Using both multinomial sampling and diffusion approximations for the analysis of the genetic distance introduced by Goldstein et al. [Goldstein, D. B., Linares, A. R., Cavalli-Sforza, L. L. & Feldman, M. W. (1995) Proc. Natl. Acad. Sci. USA 92, 6723-6727], we show that this distance follows a chi 2 distribution with degrees of freedom equal to the number of loci when there is no variation in mutation rates among the loci. In the presence of such variation, the variance of the distance is obtained. We conclude that the number of microsatellite loci required for the construction of phylogenetic trees with reliable branch lengths may be several hundred. Also, mutations that change repeat scores by several units, even though extremely rare, may dramatically influence estimates of population parameters.

Human immunodeficiency virus reverse transcriptase substrate-induced conformational changes and the mechanism of inhibition by nonnucleoside inhibitors.

A combination of transient kinetic and equilibrium titration methods has been used to show that both primer/template and nucleotide binding to human immunodeficiency virus type 1 (HIV-1) reverse transcriptase are two-step processes. In both cases, after initial formation of relatively weakly bound states, isomerization reactions lead to tightly bound states. In the case of deoxynucleotide binding to the reverse transcriptase-primer/template complex, the second step in the interaction is rate-limiting in the overall reaction during processive polymerization. Discrimination against incorrect nucleotides occurs both in the initial weak binding and in the second step but is purely kinetic in the second step (as opposed to thermodynamic in the first step). Nonnucleoside inhibitors have a relatively small effect on nucleotide-binding steps (overall affinity is reduced by a factor of ca. 10), while the affinity of the primer/template duplex is increased by at least a factor of 10. The major effect of nonnucleoside inhibitors is on the chemical step (nucleotide transfer).

Low concentrations of diacylglycerol promote the binding of apolipophorin III to a phospholipid bilayer: a surface plasmon resonance spectroscopy study.

The binding of the exchangeable apolipoprotein apolipophorin III (apoLp-III) to an egg phosphatidylcholine bilayer as a function of the concentration of diacylglycerol (DG) in the bilayer was studied by surface plasmon resonance spectroscopy. At a DG concentration of 2 mol % in the bilayer, the binding of apoLp-III reached saturation. Under saturating conditions, apoLp-III forms a closely packed monolayer approximately 55 A thick, in which each molecule of protein occupies approximately 500 A2 at the membrane surface. These dimensions are consistent with the molecular size of the apoLp-III molecule determined by x-ray crystallography, if apoLp-III binds to the bilayer with the long axis of the apoLp-III normal to the membrane surface. In the absence of protein, the overall structure of the lipid bilayer was not significantly changed up to 2.5 mol% DG. However, at 4 and 6 mol % DG, the presence of nonbilayer structures was observed. The addition of apoLp-III to a membrane containing 6 mol % DG promoted the formation of large lipid-protein complexes. These data support a two-step sequential binding mechanism for binding of apoLp-III to a lipid surface. The first step is a recognition process, consisting of the adsorption of apoLp-III to a nascent hydrophobic defect in the phospholipid bilayer caused by the presence of DG. This recognition process might depend on the presence of a hydrophobic sensor located at one of the ends of the long axis of the apoLp-III molecule but would be consolidated through H-bond and electrostatic interactions. Once primary binding is achieved, subsequent enlargement of the hydrophobic defect in the lipid surface would trigger the unfolding of the apolipoprotein and binding via the amphipathic alpha-helices. This two-step sequential binding mechanism could be a general mechanism for all exchangeable apolipoproteins. A possible physiological role of the ability of apoLp-III to bind to lipid structures in two orientations is also proposed.

The First Step of Gibberellin Biosynthesis in Pumpkin Is Catalyzed by at Least Two Copalyl Diphosphate Synthases Encoded by Differentially Regulated Genes

The first step in gibberellin biosynthesis is catalyzed by copalyl diphosphate synthase (CPS) and ent-kaurene synthase. We have cloned from pumpkin (Cucurbita maxima L.) two cDNAs, CmCPS1 and CmCPS2, that each encode a CPS. Both recombinant fusion CmCPS proteins were active in vitro. CPS are translocated into plastids and processed by cleavage of transit peptides. For CmCPS1 and CmCPS2, the putative transit peptides cannot exceed the first 99 and 107 amino acids, respectively, because longer N-terminal deletions abolished activity. Levels of both CmCPS transcripts were strictly regulated in an organ-specific and developmental manner. Both transcripts were almost undetectable in leaves and were abundant in petioles. CmCPS1 transcript levels were high in young cotyledons and low in roots. In contrast, CmCPS2 transcripts were undetectable in cotyledons but present at significant levels in roots. In hypocotyls, apices, and petioles, CmCPS1 transcript levels decreased with age much more rapidly than those of CmCPS2. We speculate that CmCPS1 expression is correlated with the early stages of organ development, whereas CmCPS2 expression is correlated with subsequent growth. In contrast, C. maxima ent-kaurene synthase transcripts were detected in every organ at almost constant levels. Thus, ent-kaurene biosynthesis may be regulated through control of CPS expression.

Analysis of a peptide hormone–receptor interaction in the yeast two-hybrid system

Interaction between a peptide hormone and extracellular domains of its receptor is a crucial step for initiation of hormone action. We have developed a modification of the yeast two-hybrid system to study this interaction and have used it to characterize the interaction of insulin-like growth factor 1 (IGF-1) with its receptor by using GAL4 transcriptional regulation with a β-galactosidase assay as readout. In this system, IGF-1 and proIGF-1 bound to the cysteine-rich domain, extracellular domain, or entire IGF-1 proreceptor. This interaction was specific. Thus, proinsulin showed no significant interaction with the IGF-1 receptor, while a chimeric proinsulin containing the C-peptide of IGF-1 had an intermediate interaction, consistent with its affinity for the IGF-1 receptor. Over 2000 IGF-1 mutants were generated by PCR and screened for interaction with the color assay. About 40% showed a strong interaction, 20% showed an intermediate interaction, and 40% give little or no signal. Of 50 mutants that were sequenced, several (Leu-5 → His, Glu-9 → Val, Arg-37 → Gly, and Met-59 → Leu) appeared to enhance receptor association, others resulted in weaker receptor interaction (Tyr-31 → Phe and Ile-43 → Phe), and two gave no detectable signal (Leu-14 → Arg and Glu-46 → Ala). Using PCR-based mutagenesis with proinsulin, we also identified a gain of function mutant (proinsulin Leu-17 → Pro) that allowed for a strong IGF-1–receptor interaction. These data demonstrate that the specificity of the interaction between a hormone and its receptor can be characterized with high efficiency in the two-hybrid system and that novel hormone analogues may be found by this method.

Limited internal friction in the rate-limiting step of a two-state protein folding reaction

Small, single-domain proteins typically fold via a compact transition-state ensemble in a process well fitted by a simple, two-state model. To characterize the rate-limiting conformational changes that underlie two-state folding, we have investigated experimentally the effects of changing solvent viscosity on the refolding of the IgG binding domain of protein L. In conjunction with numerical simulations, our results indicate that the rate-limiting conformational changes of the folding of this domain are strongly coupled to solvent viscosity and lack any significant “internal friction” arising from intrachain collisions. When compared with the previously determined solvent viscosity dependencies of other, more restricted conformational changes, our results suggest that the rate-limiting folding transition involves conformational fluctuations that displace considerable amounts of solvent. Reconciling evidence that the folding transition state ensemble is comprised of highly collapsed species with these and similar, previously reported results should provide a significant constraint for theoretical models of the folding process.