Expression, purification, and characterization of cytokinin signaling intermediates: Arabidopsis histidine phosphotransfer protein 1 (AHP1) and AHP2

Key message

We have expressed, purified, and biophysically characterized recombinant AHP1 and AHP2. Also, using computational homology models for AHP1, ARR7, and AHP1–ARR7 complex, we identified three-dimensional positioning of key amino acids.

Abstract

Cytokinin signaling involves activation of Arabidopsis Response Regulators (ARRs) by Arabidopsis Histidine Phosphotransfer Proteins (AHPs) by phosphorylation. Type-A ARRs are key regulators of several developmental pathways, but the mechanism underlying this phosphorylation and activation is not known in plants. In this study, we report the successful expression and purification of recombinant AHP1 and AHP2. Biophysical characterization shows that these two recombinant proteins were purified to homogeneity and possess well-defined secondary structures. Brief attempts to purify recombinant ARR7 posed problems during size-exclusion chromatography. Nevertheless, we generated computational homology models for AHP1, ARR7, and AHP1–ARR7 complex using crystal structures of homologous proteins from other organisms. The homology models helped to identify the three-dimensional positioning of the key conserved residues of AHP1 and ARR7 involved in phosphorylation. The similarity in positioning of these residues to other homologous proteins suggests that AHPs and type-A ARRs could be structurally conserved across kingdoms. Thus, our homology models can serve as valuable tools to gain structural insights into the phosphorylation and activation of cytokinin response regulators in plants.     

UMARS: Un-MAppable Reads Solution

Background

Few high-resolution structures of integral membranes proteins are available, as crystallization of such proteins needs yet to overcome too many technical limitations. Nevertheless, prediction of their transmembrane (TM) structure by bioinformatics tools provides interesting insights on the topology of these proteins.

Methods

We describe here how to extract new information from the analysis of hydrophobicity variations or hydrophobic pulses (HPulses) in the sequence of integral membrane proteins using the Hydrophobic Pulse Predictor, a new tool we developed for this purpose. To analyze the primary sequence of 70 integral membrane proteins we defined two levels of analysis: G1-HPulses for sliding windows of n = 2 to 6 and G2-HPulses for sliding windows of n = 12 to 16.

Results

The G2-HPulse analysis of 541 transmembrane helices allowed the definition of the new concept of transmembrane unit (TMU) that groups together transmembrane helices and segments with potential adjacent structures. In addition, the G1-HPulse analysis identified helix irregularities that corresponded to kinks, partial helices or unannotated structural events. These irregularities could represent key dynamic elements that are alternatively activated depending on the channel status as illustrated by the crystal structures of the lactose permease in different conformations.

Conclusions

Our results open a new way in the understanding of transmembrane secondary structures: hydrophobicity through hydrophobic pulses strongly impacts on such embedded structures and is not confined to define the transmembrane status of amino acids.     

Novel and unexpected bacterial diversity in an arsenic-rich ecosystem revealed by culture-dependent approaches

Background

BLAST is a commonly-used software package for comparing a query sequence to a database of known sequences; in this study, we focus on protein sequences. Position-specific-iterated BLAST (PSI-BLAST) iteratively searches a protein sequence database, using the matches in round i to construct a position-specific score matrix (PSSM) for searching the database in round i?+?1. Biegert and S?ding developed Context-sensitive BLAST (CS-BLAST), which combines information from searching the sequence database with information derived from a library of short protein profiles to achieve better homology detection than PSI-BLAST, which builds its PSSMs from scratch.

Results

We describe a new method, called domain enhanced lookup time accelerated BLAST (DELTA-BLAST), which searches a database of pre-constructed PSSMs before searching a protein-sequence database, to yield better homology detection. For its PSSMs, DELTA-BLAST employs a subset of NCBI??s Conserved Domain Database (CDD). On a test set derived from ASTRAL, with one round of searching, DELTA-BLAST achieves a ROC₅₀₀₀ of 0.270 vs. 0.116 for CS-BLAST. The performance advantage diminishes in iterated searches, but DELTA-BLAST continues to achieve better ROC scores than CS-BLAST.

Conclusions

DELTA-BLAST is a useful program for the detection of remote protein homologs. It is available under the ??Protein BLAST?? link at http://blast.ncbi.nlm.nih.gov.

Reviewers

This article was reviewed by Arcady Mushegian, Nick V. Grishin, and Frank Eisenhaber.     

Hydration studies on the archaeal protein Sso7d using NMR measurements and MD simulations

Background

Cytoplasmic class XI myosins are the fastest processive motors known. This class functions in high-velocity cytoplasmic streaming in various plant cells from algae to angiosperms. The velocities at which they process are ten times faster than its closest class V homologues.

Results

To provide sequence determinants and structural rationale for the molecular mechanism of this fast pace myosin, we have compared the sequences from myosin class V and XI through Evolutionary Trace (ET) analysis. The current study identifies class-specific residues of myosin XI spread over the actin binding site, ATP binding site and light chain binding neck region. Sequences for ET analysis were accumulated from six plant genomes, using literature based text search and sequence searches, followed by triple validation viz. CDD search, string-based searches and phylogenetic clustering. We have identified nine myosin XI genes in sorghum and seven in grape by sequence searches. Both the plants possess one gene product each belonging to myosin type VIII as well. During this process, we have re-defined the gene boundaries for three sorghum myosin XI genes using fgenesh program.

Conclusion

Molecular modelling and subsequent analysis of putative interactions involving these class-specific residues suggest a structural basis for the molecular mechanism behind high velocity of plant myosin XI. We propose a model of a more flexible switch I region that contributes to faster ADP release leading to high velocity movement of the algal myosin XI.     

Functional evolution of two subtly different (similar) folds

Background

The function of proteins is a direct consequence of their three-dimensional structure. The structural classification of proteins describes the ways of folding patterns all proteins could adopt. Although, the protein folds were described in many ways the functional properties of individual folds were not studied.

Results

We have analyzed two β-barrel folds generally adopted by small proteins to be looking similar but have different topology. On the basis of the topology they could be divided into two different folds named SH3-fold and OB-fold. There was no sequence homology between any of the proteins considered. The sequence diversity and loop variability was found to be important for various binding functions.

Conclusions

The function of Oligonucleotide/oligosaccharide-binding (OB) fold proteins was restricted to either DNA/RNA binding or sugar binding whereas the Src homology 3 (SH3) domain like proteins bind to a variety of ligands through loop modulations. A question was raised whether the evolution of these two folds was through DNA shuffling.     

Three dimensional structure directs T-cell epitope dominance associated with allergy

Background

CD4+ T-cell epitope immunodominance is not adequately explained by peptide selectivity in class II major histocompatibility proteins, but it has been correlated with adjacent segments of conformational flexibility in several antigens.

Methods

The published T-cell responses to two venom allergens and two aeroallergens were used to construct profiles of epitope dominance, which were correlated with the distribution of conformational flexibility, as measured by crystallographic B factors, solvent-accessible surface, COREX residue stability, and sequence entropy.

Results

Epitopes associated with allergy tended to be excluded from and lie adjacent to flexible segments of the allergen.

Conclusion

During the initiation of allergy, the N- and/or C-terminal ends of proteolytic processing intermediates were preferentially loaded into antigen presenting proteins for the priming of CD4+ T cells.     

Relationship Between Amino Acid Properties and Protein Stability: Buried Mutations

In order to understand the mechanism of protein stability and to develop a simple method for predicting mutation-induced stability changes, we analyzed the relationship between stability changes caused by buried mutations and changes in 48 amino acid properties. As expected from the importance of hydrophobicity, properties reflecting hydrophobicity are strongly correlated with the stability of proteins. We found that subgroup classification based on secondary structure increased correlations significantly, and mutations within -strand segments correlated better than did those in -helical segments, which may result from stronger hydrophobicity of the -strands. Multiple regression analyses incorporating combinations of three properties from among all possible combinations of the 48 properties increased the correlation coefficient to 0.88 and by an average of 13% for all data sets. Analyzing the stability of tryptophan synthase mutants with Glu49 replaced by all other residues except Arg revealed that combining buriedness, solvent-accessible surface area for denatured protein, and unfolding Gibbs free energy change increased the correlation to 0.95. Consideration of sequence and structural information (neighboring residues in sequence and in space) did not significantly strengthen the correlations in buried mutations, suggesting that nonspecific interactions dominate in the interior of proteins.     

One-step concentration of malarial parasite-infected red blood cells and removal of contaminating white blood cells

Background

The aim of this study was to develop site-specific antibodies as a tool to capture Plasmodium falciparum -dihydrofolate reductase (Pf-DHFR) from blood samples from P. falciparum infected individuals in order to detect, in a sandwich ELISA, structural alterations due to point mutations in the gene coding for Pf-DHFR. Furthermore, we wanted to study the potential use of homology models in general and of Pf-DHFR in particular in predicting antigenic malarial surface epitopes.

Methods

A homology model of Pf-DHFR domain was employed to define an epitope for the development of site-specific antibodies against Pf-DHFR. The homology model suggested an exposed loop encompassing amino acid residues 64–100. A synthetic peptide of 37-mers whose sequence corresponded to the sequence of amino acid residues 64–100 of Pf-DHFR was synthesized and used to immunize mice for antibodies. Additionally, polyclonal antibodies recognizing a recombinant DHFR enzyme were produced in rabbits.

Results and conclusions

Serum from mice immunized with the 37-mer showed strong reactivity against both the immunizing peptide, recombinant DHFR and a preparation of crude antigen from P. falciparum infected red blood cells. Five monoclonal antibodies were obtained, one of which showed reactivity towards crude antigen prepared from P. falciparum infected red cells. Western blot analysis revealed that both the polyclonal and monoclonal antibodies recognized Pf-DHFR. Our study provides insight into the potential use of homology models in general and of Pf-DHFR in particular in predicting antigenic malarial surface epitopes.     

Simulating non-small cell lung cancer with a multiscale agent-based model

Background

Abel and Trevors have delineated three aspects of sequence complexity, Random Sequence Complexity (RSC), Ordered Sequence Complexity (OSC) and Functional Sequence Complexity (FSC) observed in biosequences such as proteins. In this paper, we provide a method to measure functional sequence complexity.

Methods and Results

We have extended Shannon uncertainty by incorporating the data variable with a functionality variable. The resulting measured unit, which we call Functional bit (Fit), is calculated from the sequence data jointly with the defined functionality variable. To demonstrate the relevance to functional bioinformatics, a method to measure functional sequence complexity was developed and applied to 35 protein families. Considerations were made in determining how the measure can be used to correlate functionality when relating to the whole molecule and sub-molecule. In the experiment, we show that when the proposed measure is applied to the aligned protein sequences of ubiquitin, 6 of the 7 highest value sites correlate with the binding domain.

Conclusion

For future extensions, measures of functional bioinformatics may provide a means to evaluate potential evolving pathways from effects such as mutations, as well as analyzing the internal structural and functional relationships within the 3-D structure of proteins.     

Human T Lymphotropic Virus Type 1 protein Tax reduces histone levels

Background

Human immunodeficiency virus type 1 (HIV-1) infects cells by means of ligand-receptor interactions. This lentivirus uses the CD4 receptor in conjunction with a chemokine coreceptor, either CXCR4 or CCR5, to enter a target cell. HIV-1 is characterized by high sequence variability. Nonetheless, within this extensive variability, certain features must be conserved to define functions and phenotypes. The determination of coreceptor usage of HIV-1, from its protein envelope sequence, falls into a well-studied machine learning problem known as classification. The support vector machine (SVM), with string kernels, has proven to be very efficient for dealing with a wide class of classification problems ranging from text categorization to protein homology detection. In this paper, we investigate how the SVM can predict HIV-1 coreceptor usage when it is equipped with an appropriate string kernel.

Results

Three string kernels were compared. Accuracies of 96.35% (CCR5) 94.80% (CXCR4) and 95.15% (CCR5 and CXCR4) were achieved with the SVM equipped with the distant segments kernel on a test set of 1425 examples with a classifier built on a training set of 1425 examples. Our datasets are built with Los Alamos National Laboratory HIV Databases sequences. A web server is available at http://genome.ulaval.ca/hiv-dskernel.

Conclusion

We examined string kernels that have been used successfully for protein homology detection and propose a new one that we call the distant segments kernel. We also show how to extract the most relevant features for HIV-1 coreceptor usage. The SVM with the distant segments kernel is currently the best method described.     

Expression cartography of human tissues using self organizing maps

Background

Although homology-based methods are among the most widely used methods for predicting the structure and function of proteins, the question as to whether interface sequence conservation can be effectively exploited in predicting protein-protein interfaces has been a subject of debate.

Results

We studied more than 300,000 pair-wise alignments of protein sequences from structurally characterized protein complexes, including both obligate and transient complexes. We identified sequence similarity criteria required for accurate homology-based inference of interface residues in a query protein sequence. Based on these analyses, we developed HomPPI, a class of sequence homology-based methods for predicting protein-protein interface residues. We present two variants of HomPPI: (i) NPS-HomPPI (Non partner-specific HomPPI), which can be used to predict interface residues of a query protein in the absence of knowledge of the interaction partner; and (ii) PS-HomPPI (Partner-specific HomPPI), which can be used to predict the interface residues of a query protein with a specific target protein. Our experiments on a benchmark dataset of obligate homodimeric complexes show that NPS-HomPPI can reliably predict protein-protein interface residues in a given protein, with an average correlation coefficient (CC) of 0.76, sensitivity of 0.83, and specificity of 0.78, when sequence homologs of the query protein can be reliably identified. NPS-HomPPI also reliably predicts the interface residues of intrinsically disordered proteins. Our experiments suggest that NPS-HomPPI is competitive with several state-of-the-art interface prediction servers including those that exploit the structure of the query proteins. The partner-specific classifier, PS-HomPPI can, on a large dataset of transient complexes, predict the interface residues of a query protein with a specific target, with a CC of 0.65, sensitivity of 0.69, and specificity of 0.70, when homologs of both the query and the target can be reliably identified. The HomPPI web server is available at http://homppi.cs.iastate.edu/.

Conclusions

Sequence homology-based methods offer a class of computationally efficient and reliable approaches for predicting the protein-protein interface residues that participate in either obligate or transient interactions. For query proteins involved in transient interactions, the reliability of interface residue prediction can be improved by exploiting knowledge of putative interaction partners.     

A kingdom-specific protein domain HMM library for improved annotation of fungal genomes

Background

Mimivirus isolated from A. polyphaga is the largest virus discovered so far. It is unique among all the viruses in having genes related to translation, DNA repair and replication which bear close homology to eukaryotic genes. Nevertheless, only a small fraction of the proteins (33%) encoded in this genome has been assigned a function. Furthermore, a large fraction of the unassigned protein sequences bear no sequence similarity to proteins from other genomes. These sequences are referred to as ORFans. Because of their lack of sequence similarity to other proteins, they can not be assigned putative functions using standard sequence comparison methods. As part of our genome-wide computational efforts aimed at characterizing Mimivirus ORFans, we have applied fold-recognition methods to predict the structure of these ORFans and further functions were derived based on conservation of functionally important residues in sequence-template alignments.

Results

Using fold recognition, we have identified highly confident computational 3D structural assignments for 21 Mimivirus ORFans. In addition, highly confident functional predictions for 6 of these ORFans were derived by analyzing the conservation of functional motifs between the predicted structures and proteins of known function. This analysis allowed us to classify these 6 previously unannotated ORFans into their specific protein families: carboxylesterase/thioesterase, metal-dependent deacetylase, P-loop kinases, 3-methyladenine DNA glycosylase, BTB domain and eukaryotic translation initiation factor eIF4E.

Conclusion

Using stringent fold recognition criteria we have assigned three-dimensional structures for 21 of the ORFans encoded in the Mimivirus genome. Further, based on the 3D models and an analysis of the conservation of functionally important residues and motifs, we were able to derive functional attributes for 6 of the ORFans. Our computational identification of important functional sites in these ORFans can be the basis for a subsequent experimental verification of our predictions. Further computational and experimental studies are required to elucidate the 3D structures and functions of the remaining Mimivirus ORFans.     

Species-specific protein sequence and fold optimizations

Background

An organism's ability to adapt to its particular environmental niche is of fundamental importance to its survival and proliferation. In the largest study of its kind, we sought to identify and exploit the amino-acid signatures that make species-specific protein adaptation possible across 100 complete genomes.

Results

Environmental niche was determined to be a significant factor in variability from correspondence analysis using the amino acid composition of over 360,000 predicted open reading frames (ORFs) from 17 archae, 76 bacteria and 7 eukaryote complete genomes. Additionally, we found clusters of phylogenetically unrelated archae and bacteria that share similar environments by amino acid composition clustering. Composition analyses of conservative, domain-based homology modeling suggested an enrichment of small hydrophobic residues Ala, Gly, Val and charged residues Asp, Glu, His and Arg across all genomes. However, larger aromatic residues Phe, Trp and Tyr are reduced in folds, and these results were not affected by low complexity biases. We derived two simple log-odds scoring functions from ORFs (C_G) and folds (C_F) for each of the complete genomes. C_F achieved an average cross-validation success rate of 85 ± 8% whereas the C_G detected 73 ± 9% species-specific sequences when competing against all other non-redundant C_G. Continuously updated results are available at http://genome.mshri.on.ca.

Conclusion

Our analysis of amino acid compositions from the complete genomes provides stronger evidence for species-specific and environmental residue preferences in genomic sequences as well as in folds. Scoring functions derived from this work will be useful in future protein engineering experiments and possibly in identifying horizontal transfer events.     

Localization of the Carnation Italian ringspot virus replication protein p36 to the mitochondrial outer membrane is mediated by an internal targeting signal and the TOM complex

Background

Carnation Italian ringspot virus (CIRV) is a positive-strand RNA virus that causes massive structural alterations of mitochondria in infected host cells, the most conspicuous being the formation of numerous internal vesicles/spherules that are derived from the mitochondrial outer membrane and serve as the sites for viral RNA replication. While the membrane-bound components of the CIRV replication complex, including a 36-kD RNA-binding protein (p36), are known to be essential for these changes in mitochondrial morphology and are relatively well characterized in terms of their roles in nascent viral RNA synthesis, how these proteins are specifically targeted and inserted into mitochondria is poorly defined.

Results

Here we report on the molecular signal responsible for sorting p36 to the mitochondrial outer membrane. Using a combination of gain-of-function assays with portions of p36 fused to reporter proteins and domain-swapping assays with p36 and another closely-related viral RNA-binding protein, p33, that sorts specifically to the peroxisomal boundary membrane, we show that the mitochondrial targeting information in p36 resides within its two transmembrane domains (TMDs) and intervening hydrophilic loop sequence. Comprehensive mutational analysis of these regions in p36 revealed that the primary targeting determinants are the moderate hydrophobicity of both TMDs and the positively-charged face of an amphipathic helix within the intervening loop sequence. We show also using bimolecular fluorescence complementation (BiFC) that p36 interacts with certain components of the translocase complex in the mitochondrial outer membrane (TOM), but not with the sorting and assembly machinery (SAM).

Conclusion

Our results provide insight to how viruses, such as CIRV, exploit specific host-cell protein sorting pathways to facilitate their replication. The characterization of the targeting and insertion of p36 into the mitochondrial outer membrane also sheds light on the mechanisms involved in sorting of host-cell membrane proteins to mitochondria, a process that has been largely unexplored in plants.     

A homology model of restriction endonuclease SfiI in complex with DNA

Background

Restriction enzymes (REases) are commercial reagents commonly used in recombinant DNA technologies. They are attractive models for studying protein-DNA interactions and valuable targets for protein engineering. They are, however, extremely divergent: the amino acid sequence of a typical REase usually shows no detectable similarities to any other proteins, with rare exceptions of other REases that recognize identical or very similar sequences. From structural analyses and bioinformatics studies it has been learned that some REases belong to at least four unrelated and structurally distinct superfamilies of nucleases, PD-DxK, PLD, HNH, and GIY-YIG. Hence, they are extremely hard targets for structure prediction and homology-based inference of sequence-function relationships and the great majority of REases remain structurally and evolutionarily unclassified.

Results

SfiI is a REase which recognizes the interrupted palindromic sequence 5'GGCCNNNN^NGGCC3' and generates 3 nt long 3' overhangs upon cleavage. SfiI is an archetypal Type IIF enzyme, which functions as a tetramer and cleaves two copies of the recognition site in a concerted manner. Its sequence shows no similarity to other proteins and nothing is known about the localization of its active site or residues important for oligomerization. Using the threading approach for protein fold-recognition, we identified a remote relationship between SfiI and BglI, a dimeric Type IIP restriction enzyme from the PD-DxK superfamily of nucleases, which recognizes the 5'GCCNNNN^NGGC3' sequence and whose structure in complex with the substrate DNA is available. We constructed a homology model of SfiI in complex with its target sequence and used it to predict residues important for dimerization, tetramerization, DNA binding and catalysis.

Conclusions

The bioinformatics analysis suggest that SfiI, a Type IIF enzyme, is more closely related to BglI, an "orthodox" Type IIP restriction enzyme, than to any other REase, including other Type IIF REases with known structures, such as NgoMIV. NgoMIV and BglI belong to two different, very remotely related branches of the PD-DxK superfamily: the α-class (EcoRI-like), and the β-class (EcoRV-like), respectively. Thus, our analysis provides evidence that the ability to tetramerize and cut the two DNA sequences in a concerted manner was developed independently at least two times in the evolution of the PD-DxK superfamily of REases. The model of SfiI will also serve as a convenient platform for further experimental analyses.     

From genome to function: the Arabidopsis aquaporins

Background

In the post-genomic era newly sequenced genomes can be used to deduce organismal functions from our knowledge of other systems. Here we apply this approach to analyzing the aquaporin gene family in Arabidopsis thaliana. The aquaporins are intrinsic membrane proteins that have been characterized as facilitators of water flux. Originally termed major intrinsic proteins (MIPs), they are now also known as water channels, glycerol facilitators and aqua-glyceroporins, yet recent data suggest that they facilitate the movement of other low-molecular-weight metabolites as well.

Results

The Arabidopsis genome contains 38 sequences with homology to aquaporin in four subfamilies, termed PIP, TIP, NIP and SIP. We have analyzed aquaporin family structure and expression using the A. thaliana genome sequence, and introduce a new NMR approach for the purpose of analyzing water movement in plant roots in vivo.

Conclusions

Our preliminary data indicate a strongly transcellular component for the flux of water in roots.     

Genome-wide acceleration of protein evolution in flies (Diptera)

Background

The rate of molecular evolution varies widely between proteins, both within and among lineages. To what extent is this variation influenced by genome-wide, lineage-specific effects? To answer this question, we assess the rate variation between insect lineages for a large number of orthologous genes.

Results

When compared to the beetle Tribolium castaneum, we find that the stem lineage of flies and mosquitoes (Diptera) has experienced on average a 3-fold increase in the rate of evolution. Pairwise gene comparisons between Drosophila and Tribolium show a high correlation between evolutionary rates of orthologous proteins.

Conclusion

Gene specific divergence rates remain roughly constant over long evolutionary times, modulated by genome-wide, lineage-specific effects. Among the insects analysed so far, it appears that the Tribolium genes show the lowest rates of divergence. This has the practical consequence that homology searches for human genes yield significantly better matches in Tribolium than in Drosophila. We therefore suggest that Tribolium is better suited for comparisons between phyla than the widely employed dipterans.     

Optimized submerged batch fermentation strategy for systems scale studies of metabolic switching in Streptomyces coelicolor A3(2)

Background

Given the complex mechanisms underlying biochemical processes systems biology researchers tend to build ever increasing computational models. However, dealing with complex systems entails a variety of problems, e.g. difficult intuitive understanding, variety of time scales or non-identifiable parameters. Therefore, methods are needed that, at least semi-automatically, help to elucidate how the complexity of a model can be reduced such that important behavior is maintained and the predictive capacity of the model is increased. The results should be easily accessible and interpretable. In the best case such methods may also provide insight into fundamental biochemical mechanisms.

Results

We have developed a strategy based on the Computational Singular Perturbation (CSP) method which can be used to perform a "biochemically-driven" model reduction of even large and complex kinetic ODE systems. We provide an implementation of the original CSP algorithm in COPASI (a COmplex PAthway SImulator) and applied the strategy to two example models of different degree of complexity - a simple one-enzyme system and a full-scale model of yeast glycolysis.

Conclusion

The results show the usefulness of the method for model simplification purposes as well as for analyzing fundamental biochemical mechanisms. COPASI is freely available at http://www.copasi.org.     

An Experimentally Based Computer Search Identifies Unstructured Membrane-binding Sites in Proteins: APPLICATION TO CLASS I MYOSINS, PAKS, AND CARMIL

Programs exist for searching protein sequences for potential membrane-penetrating segments (hydrophobic regions) and for lipid-binding sites with highly defined tertiary structures, such as PH, FERM, C2, ENTH, and other domains. However, a rapidly growing number of membrane-associated proteins (including cytoskeletal proteins, kinases, GTP-binding proteins, and their effectors) bind lipids through less structured regions. Here, we describe the development and testing of a simple computer search program that identifies unstructured potential membrane-binding sites. Initially, we found that both basic and hydrophobic amino acids, irrespective of sequence, contribute to the binding to acidic phospholipid vesicles of synthetic peptides that correspond to the putative membrane-binding domains of Acanthamoeba class I myosins. Based on these results, we modified a hydrophobicity scale giving Arg- and Lys-positive, rather than negative, values. Using this basic and hydrophobic scale with a standard search algorithm, we successfully identified previously determined unstructured membrane-binding sites in all 16 proteins tested. Importantly, basic and hydrophobic searches identified previously unknown potential membrane-binding sites in class I myosins, PAKs and CARMIL (capping protein, Arp2/3, myosin I linker; a membrane-associated cytoskeletal scaffold protein), and synthetic peptides and protein domains containing these newly identified sites bound to acidic phospholipids in vitro.