Mutational analysis of the mitochondrial Rieske iron-sulfur protein of Saccharomyces cerevisiae. III. Import, protease processing, and assembly into the cytochrome bc1 complex of iron-sulfur protein lacking the iron-sulfur cluster.

We have used site-directed mutagenesis of the Saccharomyces cerevisiae Rieske iron-sulfur protein gene (RIP 1) to convert cysteines 159, 164, 178, and 180 to serines, and to convert histidines 161 and 181 to arginines. These 4 cysteines and 2 histidines are conserved in all Rieske proteins sequenced to date, and 4 of these 6 residues are thought to ligate the iron-sulfur cluster to the apoprotein. We have also converted histidine 184 to arginine. This histidine is conserved only in respiring organisms. The site-directed mutations of the six fully conserved putative iron-sulfur cluster ligands result in an inactive iron-sulfur protein, lacking iron-sulfur cluster, and failure of the yeast to grow on nonfermentable carbon sources. In contrast, when histidine 184 is replaced by arginine, the iron-sulfur cluster is assembled properly and the yeast grow on nonfermentable carbon sources. The site-directed mutations of the 6 fully conserved residues do not prevent post-translational import of iron-sulfur protein precursor into mitochondria, nor do the mutations prevent processing of iron-sulfur protein precursor to mature size protein by mitochondrial proteases. Optical spectra of mitochondria from the six mutants indicate that cytochrome b is normal, in contrast to the deranged spectrum of cytochrome b which results when the iron-sulfur protein gene is deleted. In addition, mature size iron-sulfur apoprotein is associated with cytochrome bc1 complex purified from a site-directed mutant in which iron-sulfur cluster is not inserted. These results indicate that mature size iron-sulfur apoprotein, lacking iron-sulfur cluster, is inserted into the cytochrome bc1 complex, where it interacts with and preserves the optical properties of cytochrome b. Insertion of the iron-sulfur cluster is not an obligatory prerequisite to processing of the protein to its final size. Either the processing protease cannot distinguish between iron-sulfur protein with or without the iron-sulfur cluster, or insertion of the iron-sulfur cluster occurs after the protein is processed to its mature size, possibly after it is assembled in the cytochrome bc1 complex.     

The tryptophan cluster: a hypothetical structure of the DNA-binding domain of the myb protooncogene product

In the DNA-binding domain of the c-myb protooncogene product (c-Myb) which consists of three repeats of 51-52 amino acids, there are 3 perfectly conserved tryptophans in each repeat. Site-directed mutagenesis of these tryptophans showed that any single or multiple mutations of tryptophan to hydrophilic residues or alanine abolished or greatly reduced the sequence-specific DNA-binding activity, but mutations to hydrophobic amino acids retained considerable activity. Raman spectroscopic study showed that these tryptophans were buried in the protein core. These 3 tryptophans are proposed to form a cluster in the hydrophobic core in each repeat. This hypothetical structure is referred to as the "tryptophan cluster," and it may represent a characteristic property of a group of DNA-binding proteins including the myb- and ets-related proteins.     

Structural characterization of the second TSP1-module of human thrombospondin

The TSP1-module has been first identified as the type 1 repeat of thrombospondin-1. Members of this extracellular module-family have since been shown to be present in several hundred metazoan proteins as well as in proteins of some protists. Despite the widespread occurrence and biological importance of this module-type, relatively little is known about their three-dimensional structure. To define the structural features of this important module-family, we have expressed the second TSP1-domain of human thrombospondin 1 in Escherichia coli. Amino acid sequencing of proteolytic fragments of the recombinant protein have shown that its disulfide bonds connect the six conserved cysteines in a 1-5, 2-6, 3-4 pattern. Circular dichroism studies on the recombinant protein indicate that the disulfide-bonded TSP1-module consists primarily of distorted beta-strands.     

Fluorescence probing of the function-specific cysteines of rat microsomal NADPH-cytochrome P-450 reductase

Titration of NADPH-cytochrome P-450 reductase with a fluorigenic maleimide suggests that approximately four cysteines are initially accessible and in close proximity to four tryptophans. Perturbation of the cysteines and/or tryptophans results in concomitant decreases in enzymic activity. These cysteines were correlated with functional components by binding studies and subsequent tryptic peptide mapping on the acid mobile phase-reverse phase HPLC. Adenine nucleotides and cytochrome c block labelling of the more hydrophilic peptides, while detergents facilitate labelling of the more hydrophobic peptides. The more hydrophobic peptides contain the microsomal binding site of cytochrome P-450. Removal of the prosthetic flavins exposes more cysteines in the more hydrophilic and hydrophobic regions of the peptide map, associating the former with FAD and the latter with FMN binding sites.     

O-fucosylation stabilizes the TSR3 motif in thrombospondin-1 by interacting with nearby amino acids and protecting a disulfide bond

Thrombospondin type-1 repeats (TSRs) are small protein motifs containing six conserved cysteines forming three disulfide bonds that can be modified with an O-linked fucose. Protein O-fucosyltransferase 2 (POFUT2) catalyzes the addition of O-fucose to TSRs containing the appropriate consensus sequence, and the O-fucose modification can be elongated to a Glucose-Fucose disaccharide with the addition of glucose by β3-glucosyltransferase (B3GLCT). Elimination of Pofut2 in mice results in embryonic lethality in mice, highlighting the biological significance of O-fucose modification on TSRs. Knockout of POFUT2 in HEK293T cells has been shown to cause complete or partial loss of secretion of many proteins containing O-fucosylated TSRs. In addition, POFUT2 is localized to the endoplasmic reticulum (ER) and only modifies folded TSRs, stabilizing their structures. These observations suggest that POFUT2 is involved in an ER quality control mechanism for TSR folding and that B3GLCT also participates in quality control by providing additional stabilization to TSRs. However, the mechanisms by which addition of these sugars result in stabilization are poorly understood. Here, we conducted molecular dynamics (MD) simulations and provide crystallographic and NMR evidence that the Glucose-Fucose disaccharide interacts with specific amino acids in the TSR3 domain in thrombospondin-1 that are within proximity to the O-fucosylation modification site resulting in protection of a nearby disulfide bond. We also show that mutation of these amino acids reduces the stabilizing effect of the sugars in vitro. These data provide mechanistic details regarding the importance of O-fucosylation and how it participates in quality control mechanisms inside the ER.     

Physical characterization of the procollagen module of human thrombospondin 1 expressed in insect cells

Thrombospondin 1 (TSP1) is a homotrimeric glycoprotein composed of 150-kDa subunits connected by disulfide bridges. The procollagen module of thrombospondin 1 has been implicated in antiangiogenic activity. Procollagen modules are found in a number of extracellular proteins and are identifiable by 10 cysteines with characteristic spacing. We expressed and studied the procollagen module (C) of human TSP1, both by itself and in the context of the adjoining oligomerization sequence (o) and N-terminal module (N). The coding sequences were introduced into baculoviruses along with an N-terminal signal sequence and C-terminal polyhistidine tag. Proteins were purified from conditioned medium of infected insect cells by nickel-chelate chromatography. NoC is a disulfide bonded trimer and cleaves readily at a site of preferential proteolysis to yield monomeric N and trimeric oC. These are known properties of full-length TSP1. Mass spectroscopy indicated that C is N-glycosylated, and all 10 cysteine residues of C are in disulfides. By equilibrium ultracentrifugation, C is a monomer in physiological salt solution. Circular dichroism, intrinsic fluorescence, and differential scanning calorimetry experiments suggest that the stability of C is determined by the disulfides. The two tryptophans of C are in a polar, exposed environment as assessed by iodide fluorescence quenching and solvent perturbation. The oC far UV circular dichroism spectrum could be modeled as the sum of C and a coiled-coil oligomerization domain. The results indicate that the recombinant C folds autonomously into its native structure, and trimerization of the modules in TSP1 does not perturb their structures.     

Thrombospondin 1 and thrombospondin 2 are expressed as both homo- and heterotrimers.

There exist two distinct thrombospondin molecules (designated TSP1 and TSP2) which are encoded by separate genes. TSP1 is a trimeric cell surface and extracellular matrix molecule. Sequence comparison reveals that the 2 cysteines involved in interchain disulfide linkage and trimer assembly in TSP1 are conserved in TSP2 (Laherty, C. D., O'Rourke, K., Wolf, F. W., Katz, R., Seldin, M. F., and Dixit, V. M. (1992) J. Biol. Chem. 267, 3274-3281). Swiss 3T3 fibroblasts express both TSP1 and TSP2, and, therefore, an important question is whether TSP in such cells is expressed as homotrimers or as heterotrimers. We find that Swiss 3T3 cells and epithelial cells transfected with TSP expression vectors express both homo- and heterotrimeric forms of TSP. In addition, homotrimeric TSP2 has a lower affinity for heparin than homotrimeric TSP1. Thus, the heparin affinity of TSP can be modulated by the expression of TSP as homo- or heterotrimers.     

Structural organization of essential iron-sulfur clusters in the evolutionarily highly conserved ATP-binding cassette protein ABCE1

The ABC protein ABCE1, formerly named RNase L inhibitor RLI1, is one of the most conserved proteins in evolution and is expressed in all organisms except eubacteria. Because of its fundamental role in translation initiation and/or ribosome biosynthesis, ABCE1 is essential for life. Its molecular mechanism has, however, not been elucidated. In addition to two ABC ATPase domains, ABCE1 contains a unique N-terminal region with eight conserved cysteines, predicted to coordinate iron-sulfur clusters. Here we present detailed information on the type and on the structural organization of the Fe-S clusters in ABCE1. Based on biophysical, biochemical, and yeast genetic analyses, ABCE1 harbors two essential diamagnetic [4Fe-4S](2+) clusters with different electronic environments, one ferredoxin-like (CPX(n)CX(2)CX(2)C; Cys at positions 4-7) and one unique ABCE1-type cluster (CXPX(2)CX(3)CX(n)CP; Cys at positions 1, 2, 3, and 8). Strikingly, only seven of the eight conserved cysteines coordinating the Fe-S clusters are essential for cell viability. Mutagenesis of the cysteine at position 6 yielded a functional ABCE1 with the ferredoxin-like Fe-S cluster in a paramagnetic [3Fe-4S](+) state. Notably, a lethal mutation of the cysteine at position 4 can be rescued by ligand swapping with an adjacent, extra cysteine conserved among all eukaryotes.     

Structural and functional characterization of hBD-1(Ser35), a peptide deduced from a DEFB1 polymorphism

beta-Defensins are mammalian antimicrobial peptides that share a unique disulfide-bonding motif of six conserved cysteines. An intragenic polymorphism of the DEFB1 gene that changes a highly conserved Cys to Ser in the peptide coding region has recently been described. The deduced peptide cannot form three disulfide bonds, as one of the cysteines is unpaired. We have determined the cysteine connectivities of a corresponding synthetic hBD-1(Ser35) peptide, investigated the structure by circular dichroism spectroscopy, and assayed the in vitro antimicrobial activity. Despite a different arrangement of the disulfides, hBD-1(Ser35) proved as active as hBD-1 against the microorganisms tested. This activity likely depends on the ability of hBD-1(Ser35) to adopt an amphipathic conformation in hydrophobic environment, similar to the wild type peptide, as suggested by CD spectroscopy.     

Collagen-binding domains of gelatinase A and thrombospondin-derived peptides impede endocytic clearance of active gelatinase A and promote HT1080 fibrosarcoma cell invasion

Gelatinase A (matrix metalloproteinase-2, MMP-2) binds to several proteins through its collagen-binding domains (CBDs). Surface plasmon resonance analysis revealed a strong interaction between CBD123 and thrombospondin-1 (TSP-1), with a K(D) value of 2x10(-9) M. CBD123, as well as individual domains, behave as competitive inhibitors of the TSP-1-directed endocytic clearance of active MMP-2, but not of its latent form, by HT1080 fibrosarcoma cells. Enhanced level of active MMP-2 in conditioned medium was associated to increased matrigel invasion. Similarly, GGWSHWSPWSS and GGWSHW peptides, as tryptophan-rich peptides within properdin-repeat motifs (TSRs) of TSP-1, promoted MMP-2 accumulation and cell invasiveness. Our data document the importance of TSP-1 in promoting MMP-2-mediated cancer cell invasion through interaction between CBDs of the enzyme and TSRs motifs of TSP-1.     

Identification by mutagenesis of conserved arginine and tryptophan residues in rat liver carnitine palmitoyltransferase I important for catalytic activity

Carnitine palmitoyltransferase I catalyzes the conversion of long-chain acyl-CoA to acylcarnitines in the presence of l-carnitine. To determine the role of the conserved arginine and tryptophan residues on catalytic activity in the liver isoform of carnitine palmitoyltransferase I (L-CPTI), we separately mutated five conserved arginines and two tryptophans to alanine. Substitution of arginine residues 388, 451, and 606 with alanine resulted in loss of 88, 82, and 93% of L-CPTI activity, respectively. Mutants R601A and R655A showed less than 2% of the wild type L-CPTI activity. A change of tryptophan 391 and 452 to alanine resulted in 50 and 93% loss in carnitine palmitoyltransferase activity, respectively. The mutations caused decreases in catalytic efficiency of 80-98%. The residual activity in the mutant L-CPTIs was sensitive to malonyl-CoA inhibition. Mutants R388A, R451A, R606A, W391A, and W452A had no effect on the K(m) values for carnitine or palmitoyl-CoA. However, these mutations decreased the V(max) values for both substrates by 10-40-fold, suggesting that the main effect of the mutations was to decrease the stability of the enzyme-substrate complex. We suggest that conserved arginine and tryptophan residues in L-CPTI contribute to the stabilization of the enzyme-substrate complex by charge neutralization and hydrophobic interactions. The predicted secondary structure of the 100-amino acid residue region of L-CPTI, containing arginines 388 and 451 and tryptophans 391 and 452, consists of four alpha-helices similar to the known three-dimensional structure of the acyl-CoA-binding protein. We predict that this 100-amino acid residue region constitutes the putative palmitoyl-CoA-binding site in L-CPTI.     

Forced Dissociation of the Strand Dimer Interface between C-Cadherin Ectodomains

The force-induced dissociation of the strand dimer interface in C-cadherin has been studied using steered molecular dynamics simulations. The dissociation occurred, without domain unraveling, after the extraction of the conserved trypthophans (Trp2) from their respective hydrophobic pockets. The simulations revealed two stable positions for the Trp2 side chain inside the pocket. The most internal stable position involved a hydrogen bond between the ring Ne of Trp2 and the backbone carbonyl of Glu90. In the second stable position, the aromatic ring is located at the pocket entrance. After extracting the two tryptophans from their pockets, the complex exists in an intermediate bound state that involves a close packing of the tryptophans with residues Asp1 and Asp27 from both domains. Dissociation occurred after this residue association was broken. Simulations carried out with a complex formed between W2A mutants showed that the mutant complex dissociates more easily than the wild type complex does. These results correlate closely with the role of the conserved tryptophans suggested previously by site directed mutagenesis.     

Localization of disulfide bonds in the cystine knot domain of human von Willebrand factor

von Willebrand factor (VWF) is a multimeric glycoprotein that is required for normal hemostasis. After translocation into the endoplasmic reticulum, proVWF subunits dimerize through disulfide bonds between their C-terminal cystine knot-like (CK) domains. CK domains are characterized by six conserved cysteines. Disulfide bonds between cysteines 2 and 5 and between cysteines 3 and 6 define a ring that is penetrated by a disulfide bond between cysteines 1 and 4. Dimerization often is mediated by additional cysteines that differ among CK domain subfamilies. When expressed in a baculovirus system, recombinant VWF CK domains (residues 1957-2050) were secreted as dimers that were converted to monomers by selective reduction and alkylation of three unconserved cysteine residues: Cys(2008), Cys(2010), and Cys(2048). By partial reduction and alkylation, chemical and proteolytic digestion, mass spectrometry, and amino acid sequencing, the remaining intrachain disulfide bonds were characterized: Cys(1961)-Cys(2011) (), Cys(1987)-Cys(2041) (), Cys(1991)-Cys(2043) (), and Cys(1976)-Cys(2025). The mutation C2008A or C2010A prevented dimerization, whereas the mutation C2048A did not. Symmetry considerations and molecular modeling based on the structure of transforming growth factor-beta suggest that one or three of residues Cys(2008), Cys(2010), and Cys(2048) in each subunit mediate the covalent dimerization of proVWF.     

In vitro mutagenesis of Caenorhabditis elegans cuticle collagens identifies a potential subtilisin-like protease cleavage site and demonstrates that carboxyl domain disulfide bonding is required for normal function but not assembly.

The importance of conserved amino acids in the amino and carboxyl non-Gly-X-Y domains of Caenorhabditis elegans cuticle collagens was examined by analyzing site-directed mutations of the sqt-1 and rol-6 collagen genes in transgenic animals. Altered collagen genes on transgenic arrays were shown to produce appropriate phenotypes by injecting in vivo cloned mutant alleles. Equivalent alterations in sqt-1 and rol-6 generally produced the same phenotypes, indicating that conserved amino acids in these two collagens have similar functions. Serine substitutions for either of two conserved carboxyl domain cysteines produced LRol phenotypes. Substitution for both cysteines in sqt-1 also resulted in an LRol phenotype, demonstrating that disulfide bonding is important for normal function but not required for assembly. Arg-1 or Arg-4 to Cys mutations in homology block A (HBA; consensus, 1-RXRRQ-5; in the amino non-Gly-X-Y domain) caused RRol phenotypes, while the same alteration at Arg-3 had no effect, indicating that Arg-3 is functionally different from Arg-1 and Arg-4. Substitutions of Arg-4 with Ser, Leu, or Glu also produced the RRol phenotype, while Lys substitutions for Arg-1 or Arg-4 did not generate any abnormal phenotypes. His substitutions for Arg-1 or Arg-4 caused somewhat less severe RRol phenotypes. Therefore, strong positively charged residues, Arg or Lys, are required at positions 1 and 4 for normal function. The conserved pattern of arginines in HBA matches the cleavage sites of the subtilisin-like endoproteinases. HBA may be a cleavage site for a subtilisin-like protease, and cleavage may be important for cuticle collagen processing.     

The crystal structure of annexin Gh1 from Gossypium hirsutum reveals an unusual S3 cluster.

The three-dimensional crystal structure of recombinant annexin Gh1 from Gossypium hirsutum (cotton fibre) has been determined and refined to the final R-factor of 0.219 at the resolution of 2.1 A. This plant annexin consists of the typical 'annexin fold' and is similar to the previously solved bell pepper annexin Anx24(Ca32), but significant differences are seen when compared to the structure of nonplant annexins. A comparison with the structure of the mammalian annexin AnxA5 indicates that canonical calcium binding is geometrically possible within the membrane loops in domains I and II of Anx(Gh1) in their present conformation. All plant annexins possess a conserved tryptophan residue in the AB loop of the first domain; this residue was found to adopt both a loop-in and a loop-out conformation in the bell pepper annexin Anx24(Ca32). In Anx(Gh1), the conserved tryptophan residue is in a surface-exposed position, half way between both conformations observed in Anx24(Ca32). The present structure reveals an unusual sulfur cluster formed by two cysteines and a methionine in domains II and III, respectively. While both cysteines adopt the reduced thiolate forms and are separated by a distance of about 5.5 A, the sulfur atom of the methionine residue is in their close vicinity and apparently interacts with both cysteine sulfur atoms. While the cysteine residues are conserved in at least five plant annexins and in several mammalian members of the annexin family of proteins, the methionine residue is conserved only in three plant proteins. Several of these annexins carrying the conserved residues have been implicated in oxidative stress response. We therefore hypothesize that the cysteine motif found in the present structure, or possibly even the entire sulfur cluster, forms the molecular basis for annexin function in oxidative stress response.     

Isolation and amino-acid sequence analysis of human sperm protamines P1 and P2. Occurrence of two forms of protamine P2

The two protamines of human sperm cell nuclei, P1 and P2, were isolated in pure form after extraction with 6M guanidine/5% mercaptoethanol and alkylation with vinyl pyridine by reversed-phase high-performance liquid chromatography. The amino-acid sequence of protamine P1 was determined by analysing the intact protein and the fragments obtained by cyanogen bromide cleavage. Out of the 50 amino-acid residues 24 are arginines and 6 are cysteines. The sequence of protamine P2 was determined by analysing the intact protein and the fragments resulting from cleavage with endoproteinase Lys-C and thermolysin. Protamine P2 was found to occur in two forms which only differ in their N-terminal regions. The form P2' is three amino-acid residues longer at the N-terminus than the form P2'. Out of the 57 amino-acid residues in the longer form 27 are arginines and 5 are cysteines. Human protamine P1 is highly homologous with the protamines isolated from bull, boar, ram and mouse sperm cells, but human protamine P2 shows a novel type of structure, although also here the dominant amino acids are arginine and cysteine.     

Structural and Functional Similarities of Calcium Homeostasis Modulator 1 (CALHM1) Ion Channel with Connexins,Pannexins, and Innexins

CALHM1 (calcium homeostasis modulator 1) forms a plasma membrane ion channel that mediates neuronal excitability in response to changes in extracellular Ca²⁺ concentration. Six human CALHM homologs exist with no homology to other proteins, although CALHM1 is conserved across >20 species. Here we demonstrate that CALHM1 shares functional and quaternary and secondary structural similarities with connexins and evolutionarily distinct innexins and their vertebrate pannexin homologs. A CALHM1 channel is a hexamer, comprised of six monomers, each of which possesses four transmembrane domains, cytoplasmic amino and carboxyl termini, an amino-terminal helix, and conserved extracellular cysteines. The estimated pore diameter of the CALHM1 channel is ∼14 Å, enabling permeation of large charged molecules. Thus, CALHMs, connexins, and pannexins and innexins are structurally related protein families with shared and distinct functional properties.     

The amino acid sequence of boar H1t, a testis-specific H1 histone variant

H1t is a testis-specific H1 variant found in spermatocytes and spermatids of mammals. The complete amino acid sequence has been determined for H1t isolated from boar testes. The protein is composed of 211 amino acids with the composition 3 aspartic acids, 6 asparagines, 13 threonines, 21 serines, 7 glutamic acids, 7 glutamines, 15 prolines, 14 glycines, 37 alanines, 11 valines, 1 methionine, 3 isoleucines, 14 leucines, 1 tyrosine, 1 phenylalanine, 42 lysines, 15 arginines and a calculated molecular weight of 22,059 disregarding the post-translational acetylation of the amino-terminal alanine. The protein shows typical H1 domains with basic amino- and carboxyl-terminal regions separated by a conserved, presumably globular, core sequence. The core region is very homologous to the highly conserved core sequence found in somatic mammalian H1 histones but does differ from this sequence in 15 places. Accordingly, it may be appropriate to consider H1t as new variant category. The carboxyl-terminal half of H1t is distinguished from the standard somatic family by being somewhat shorter and by the presence of 10 arginine residues. In contrast to many H1 proteins, the carboxyl-terminal region of H1t does not show an obvious pattern of peptide repeats.     

Both aromatic and cationic residues contribute to the membrane-lytic and bactericidal activity of eosinophil cationic protein

Eosinophil cationic protein (ECP) and eosinophil derived neurotoxin (EDN) are proteins of the ribonuclease A (RNase A) superfamily that have developed biological properties related to the function of eosinophils. ECP is a potent cytotoxic molecule, and although the mechanism is still unknown this cytotoxic activity has been associated with its highly cationic character. Using liposome vesicles as a model, we have demonstrated that ECP tends to disrupt preferentially acidic membranes. On the basis of structure analysis, ECP variants modified at basic and hydrophobic residues have been constructed. Changes in the leakage of liposome vesicles by these ECP variants have indicated the role of both aromatic and basic specific amino acids in cellular membrane disruption. This is the case with the two tryptophans at positions 10 and 35, but not phenylalanine 76, and the two arginines 101 and 104. The bactericidal activity of both native ECP and point-mutated variants, tested against Escherichia coli and Staphylococcus aureus, suggests that basic amino acids play, in addition to the effect on the disruption of the cellular membrane, other roles such as specific binding on the surface of the bacteria cell.     

Computer modeling of the membrane interaction of FYVE domains

FYVE domains are membrane targeting domains that are found in proteins involved in endosomal trafficking and signal transduction pathways. Most FYVE domains bind specifically to phosphatidylinositol 3-phosphate (PI(3)P), a lipid that resides mainly in endosomal membranes. Though the specific interactions between FYVE domains and the headgroup of PI(3)P have been well characterized, principally through structural studies, the available experimental structures suggest several different models for FYVE/membrane association. Thus, the manner in which FYVE domains adsorb to the membrane surface remains to be elucidated. Towards this end, recent experiments have shown that FYVE domains bind PI(3)P in the context of phospholipid bilayers and that hydrophobic residues on a conserved loop are able to penetrate the membrane interface in a PI(3)P-dependent manner.Here, the finite difference Poisson-Boltzmann (FDPB) method has been used to calculate the energetic interactions of FYVE domains with phospholipid membranes. Based on the computational analysis, it is found that (1) recruitment to membranes is facilitated by non-specific electrostatic interactions that occur between basic residues on the domains and acidic phospholipids in the membrane, (2) the energetic analysis can quantitatively differentiate among the modes of membrane association proposed by the experimentally determined structures, (3) FDPB calculations predict energetically feasible models for the membrane-associated states of FYVE domains, (4) these models are consistent with the observation that conserved hydrophobic residues insert into the membrane interface, and (5) the calculations provide a molecular model for the hydrophobic partitioning: binding of PI(3)P significantly neutralizes positive potential in the region of the hydrophobic residues, which acts as an "electrostatic switch" by reducing the energetic barrier for membrane penetration. Finally, the computational results are extended to FYVE domains of unknown structure through the construction of high quality homology models for human FYVE sequences.