Functional analysis of predicted coiled-coil regions in the Escherichia coli K-12 O-antigen polysaccharide chain length determinant Wzz

Wzz is a membrane protein that determines the chain length distribution of the O-antigen lipopolysaccharide by an unknown mechanism. Wzz proteins consist of two transmembrane helices separated by a large periplasmic loop. The periplasmic loop of Escherichia coli K-12 Wzz (244 amino acids from K65 to A308) was purified and found to be a monomer with an extended conformation, as determined by gel filtration chromatography and analytical ultracentrifugation. Circular dichroism showed that the loop has a 60% helical content. The Wzz periplasmic loop also contains three regions with predicted coiled coils. To probe the function of the predicted coiled coils, we constructed amino acid replacement mutants of the E. coli K-12 Wzz protein, which were designed so that the coiled coils could be separate without compromising the helicity of the individual molecules. Mutations in one of the regions, spanning amino acids 108 to 130 (region I), were associated with a partial defect in O-antigen chain length distribution, while mutants with mutations in the region spanning amino acids 209 to 223 (region III) did not have an apparent functional defect. In contrast, mutations in the region spanning amino acids 153 to 173 (region II) eliminated the Wzz function. This phenotype was associated with protein instability, most likely due to conformational changes caused by the amino acid replacements, which was confirmed by limited trypsin proteolysis. Additional mutagenesis based on a three-dimensional model of region I demonstrated that the amino acids implicated in function are all located at the same face of a predicted α-helix, suggesting that a coiled coil actually does not exist in this region. Together, our results suggest that the regions predicted to be coiled coils are important for Wzz function because they maintain the native conformation of the protein, although the existence of coiled coils could not be demonstrated experimentally.     

The hinge region of type VII collagen is intrinsically disordered

Type VII collagen (Col7) is important for skin integrity. As a major component of the anchoring fibrils, Col7 is essential for linking different skin layers together. The central collagenous domain of Col7 contains several interruptions of the collagen triple helix. The longest interruption is 39 amino acids long and referred to as the hinge region. The hinge region is highly conserved between species. This region was predicted to adopt a coiled coil structure and to serve as the trimerization domain of Col7.To gain insight into the potential function of the hinge region we investigated a heterologous expressed peptide by CD and NMR spectroscopy. CD spectroscopy implies that the hinge region is intrinsically disordered. Resonance assignment was performed and allowed secondary structure analysis based on the chemical shift values. Seven amino acids in the N-terminal moiety show residual α-helical conformation. Subsequent investigation of temperature dependency of amide chemical shifts indicated participation in hydrogen bonding of amino acid residues in the C-terminal moiety of the hinge region. Therefore, the hinge region does not form a coiled coil structure under the employed experimental conditions. The intrinsic disorder of the hinge region might be desired for flexibility to serve as a “hinge” or the hinge region is an important interaction site as typically observed for intrinsically disordered proteins.     

In silico investigation of the ATP7B gene: insights from functional prediction of non-synonymous substitution to protein structure

ATP7B is a copper-transporting ATPase that plays a key role in the regulation of copper homeostasis. Mutations in the ATP7B gene are causative for Wilson’s disease, and recent reports have suggested that genetic variants are associated with susceptibility to Alzheimer’s disease. Unfortunately, it is difficult to profile experimentally novel genetic variants in the ATP7B gene, because the human protein X-ray structure is not yet entirely understood. In order to investigate ATP7B non-synonymous substitutions, we used an in silico amino acid sequence-based approach. Specifically, we analyzed 337 ATP7B non-synonymous substitutions, which included Wilson’s disease-causing mutations (DVs) and non Wilson’s disease-causing variants (NDVs), with an algorithm that estimated a combined probability (cP_del) of an amino acidic change to be deleterious for the protein function. This approach appeared to reliably indentify the probability of DVs and NDVs to be deleterious and to profile still unknown gene variants. Specifically, after analyzing ATP7B protein domains with the cP_del method, we found results in line with the predicted–modeled domains and some new suggestions. In conclusion, a functional survey of amino acid changes in the ATP7B protein is provided herein, and we suggest that this bioinformatic method can furnish information about novel ATP7B mutations. Furthermore, the same approach can be applied to other uncharacterized proteins.     

Polymorphism in a histone H1 subtype with a short N-terminal domain in three legume species (Fabaceae, Fabaeae)

A number of alleles of an orthologous gene His6 encoding histone H1 subtype f (H1-6 in pea) accumulated in chromatin of old tissues were sequenced in three legume species: seven alleles in Pisum sativum, four in Vicia unijuga and eight in Lathyrus gmelinii. In the total of 19 alleles sequenced in the three species, 29 non-synonymous substitutions and six indels were found in the coding region; most of amino acid substitutions (26 of 29) and all indels occurred in the C-terminal hydrophilic domain of the encoded protein. All species were polymorphic for some non-synonymous substitutions, V. unijuga was also polymorphic for one and P. sativum for two indels. Three near-isogenic lines of P. sativum bearing different alleles showed differences in many quantitative traits; that in the growth dynamic could be tentatively attributed to the allelic substitution of subtype H1-6. The frequencies of four electromorphs in a sampled locality of V. unijuga were found to be close to those observed 25?years ago, although their rapid change in the past was supposed in the previous study.     

Functional Coding Variants in SLC6A15, a Possible Risk Gene for Major Depression

SLC6A15 is a neuron-specific neutral amino acid transporter that belongs to the solute carrier 6 gene family. This gene family is responsible for presynaptic re-uptake of the majority of neurotransmitters. Convergent data from human studies, animal models and pharmacological investigations suggest a possible role of SLC6A15 in major depressive disorder. In this work, we explored potential functional variants in this gene that could influence the activity of the amino acid transporter and thus downstream neuronal function and possibly the risk for stress-related psychiatric disorders. DNA from 400 depressed patients and 400 controls was screened for genetic variants using a pooled targeted re-sequencing approach. Results were verified by individual re-genotyping and validated non-synonymous coding variants were tested in an independent sample (N = 1934). Nine variants altering the amino acid sequence were then assessed for their functional effects by measuring SLC6A15 transporter activity in a cellular uptake assay. In total, we identified 405 genetic variants, including twelve non-synonymous variants. While none of the non-synonymous coding variants showed significant differences in case-control associations, two rare non-synonymous variants were associated with a significantly increased maximal ³H proline uptake as compared to the wildtype sequence. Our data suggest that genetic variants in the SLC6A15 locus change the activity of the amino acid transporter and might thus influence its neuronal function and the risk for stress-related psychiatric disorders. As statistically significant association for rare variants might only be achieved in extremely large samples (N >70,000) functional exploration may shed light on putatively disease-relevant variants.     

Large-scale analysis of non-synonymous coding region single nucleotide polymorphisms

MOTIVATION: Single nucleotide polymorphisms (SNPs) are the most common form of genetic variant in humans. SNPs causing amino acid substitutions are of particular interest as candidates for loci affecting susceptibility to complex diseases, such as diabetes and hypertension. To efficiently screen SNPs for disease association, it is important to distinguish neutral variants from deleterious ones. RESULTS: We describe the use of Pfam protein motif models and the HMMER program to predict whether amino acid changes in conserved domains are likely to affect protein function. We find that the magnitude of the change in the HMMER E-value caused by an amino acid substitution is a good predictor of whether it is deleterious. We provide internet-accessible display tools for a genomewide collection of SNPs, including 7391 distinct non-synonymous coding region SNPs in 2683 genes. AVAILABILITY: http://lpgws.nci.nih.gov/cgi-bin/GeneViewer.cgi     

Identification of Site-Specific Adaptations Conferring Increased Neural Cell Tropism during Human Enterovirus 71 Infection

Enterovirus 71 (EV71) is one of the most virulent enteroviruses, but the specific molecular features that enhance its ability to disseminate in humans remain unknown. We analyzed the genomic features of EV71 in an immunocompromised host with disseminated disease according to the different sites of infection. Comparison of five full-length genomes sequenced directly from respiratory, gastrointestinal, nervous system, and blood specimens revealed three nucleotide changes that occurred within a five-day period: a non-conservative amino acid change in VP1 located within the BC loop (L97R), a region considered as an immunogenic site and possibly important in poliovirus host adaptation; a conservative amino acid substitution in protein 2B (A38V); and a silent mutation in protein 3D (L175). Infectious clones were constructed using both BrCr (lineage A) and the clinical strain (lineage C) backgrounds containing either one or both non-synonymous mutations. In vitro cell tropism and competition assays revealed that the VP1₉₇ Leu to Arg substitution within the BC loop conferred a replicative advantage in SH-SY5Y cells of neuroblastoma origin. Interestingly, this mutation was frequently associated in vitro with a second non-conservative mutation (E167G or E167A) in the VP1 EF loop in neuroblastoma cells. Comparative models of these EV71 VP1 variants were built to determine how the substitutions might affect VP1 structure and/or interactions with host cells and suggest that, while no significant structural changes were observed, the substitutions may alter interactions with host cell receptors. Taken together, our results show that the VP1 BC loop region of EV71 plays a critical role in cell tropism independent of EV71 lineage and, thus, may have contributed to dissemination and neurotropism in the immunocompromised patient.     

Angiotensinogen T174M and M235T Variants and Hypertension in the Hani and Yi Minority Groups of China

A case–control study of 538 individuals investigated whether the angiotensinogen gene (AGT) might be implicated in the pathogenesis of essential hypertension in the Hani and Yi populations of China. Genotypes for two diallelic DNA polymorphisms observed at amino acid residues 174 (T174M) and 235 (M235T) within the coding sequence were determined. M235T and T174M genotyping with PCR-RFLP was performed in 267 normotensive subjects and 271 hypertensive subjects. No significant difference was found between normotensives and hypertensives in genotype distribution and allele frequency for either M235T or T174M in the Hani or the Yi populations (P > 0.05). Relative to carriers of the 235T/235T and 174T/174T combination, the others had a significantly elevated risk of hypertension (OR = 1.62, 95% CI 1.02–2.59; P = 0.043) in the Hani population. The AGT M235T and T174M variants in combination may play a role in the genetic predisposition to develop essential hypertension in the Hani minority of China.     

Nucleotide variation of sFRP5 gene is not associated with obesity in children and adolescents

Because sFRP5 was shown to be an important extracellular modulator of the Wnt pathway, regulating adipogenesis, we wanted to investigate the role of sFRP5 variants in human, monogenic obesity by performing mutation analysis. We screened the complete sFRP5 coding region in 622 obese children and adolescents and 503 lean control individuals by high-resolution melting curve analysis and direct sequencing. We found a total of 15 sequence variants in sFRP5, 10 of which resulted in a non-synonymous amino acid change. Five of these variants were, to our knowledge, not previously reported. For one of the variants (c.-3G>A), we identified a trend towards association between the variant frequency and the obese phenotype. We argue that, when looking at conservation and location inside known protein domains, several of the identified variants (D103N, A113V, K212N and H317L), may affect sFRP5 protein function. In addition, we found c.-3G>A, residing in the Kozak sequence, with a lower frequency in cases compared to controls. However, functional studies investigating the effect of sFRP5 variants on protein function are necessary to determine the true role of sFRP5 genetic variation in human, monogenic obesity.     

Genetic diversity of the toll-like receptor 2 (TLR2) in hare (Lepus capensis) populations from Tunisia

Toll-like receptors (TLRs) are a major group of proteins that recognize molecular components of infectious agents, known as pathogen associated molecular patterns (PAMPs). The structure of these genes is similar and characterized by the presence of an ectodomain, a signal transmembrane segment and a highly conserved cytoplasmic domain. The latter domain is homologous to the human interleukin-1 receptor (IL1R) and human IL-18 receptor (IL-18R) and designated TIR domain. The latter domain of the TLR genes was suggested to be very conservative and its evolution is driven by purifying selection. Variability and evolution of the TIR sequences of TLR2 gene were studied in three hare populations from Tunisia with different ecological characteristics (NT–North Tunisia with Mediterranean, CT–Central Tunisia with semi-arid, and ST–South Tunisia with arid climate). Sequencing of a 372 bp fragment of TIR2 revealed 25 alleles among 110 hares. Twenty variable nucleotide positions were detected, of which 7 were non-synonymous. The highest variability was observed in CT, with 16 polymorphic positions. In ST, only 4 polymorphic nucleotide positions were detected with all diversity values lower than those recorded for the other two populations. By using several approaches, no positive selection was detected. However, evidence of purifying selection was found at two positions. The logistic models of the most common TIR2 protein variant that we run to examine whether its occurrence was affected by climatic variation independent of the geographic sample location suggested only a longitudinal effect. Finally, the mapping of the non-synonymous mutations to the inferred tertiary protein structure showed that they were all localized in the different loop regions. Among all non-synonymous substitutions, three were suggested to be deleterious as evidenced by PROVEAN analysis. The observed patterns of variability characterized by low genetic diversity in ST might suggest that the TIR region was more affected, than other markers, by genetic drift or/and that these patterns were shaped by different selective pressures under different ecological conditions. Notably, this low diversity was not detected by other (putatively neutral) microsatellite markers analysed in the course of other studies. But low diversity was also found for two MHC class II adaptive immune genes. As expected from functionally important regions, the evolution of the TIR2 domain is mainly driven by purifying selection. However, the occurrence of deleterious non-synonymous substitutions might highlight the flexible evolution of the TIR genes and/or their interactions with other proteins.     

Association of the OCA2 Polymorphism His615Arg with Melanin Content in East Asian Populations: Further Evidence of Convergent Evolution of Skin Pigmentation

The last decade has witnessed important advances in our understanding of the genetics of pigmentation in European populations, but very little is known about the genes involved in skin pigmentation variation in East Asian populations. Here, we present the results of a study evaluating the association of 10 Single Nucleotide Polymorphisms (SNPs) located within 5 pigmentation candidate genes (OCA2, DCT, ADAM17, ADAMTS20, and TYRP1) with skin pigmentation measured quantitatively in a sample of individuals of East Asian ancestry living in Canada. We show that the non-synonymous polymorphism rs1800414 (His615Arg) located within the OCA2 gene is significantly associated with skin pigmentation in this sample. We replicated this result in an independent sample of Chinese individuals of Han ancestry. This polymorphism is characterized by a derived allele that is present at a high frequency in East Asian populations, but is absent in other population groups. In both samples, individuals with the derived G allele, which codes for the amino acid arginine, show lower melanin levels than those with the ancestral A allele, which codes for the amino acid histidine. An analysis of this non-synonymous polymorphism using several programs to predict potential functional effects provides additional support for the role of this SNP in skin pigmentation variation in East Asian populations. Our results are consistent with previous research indicating that evolution to lightly-pigmented skin occurred, at least in part, independently in Europe and East Asia.     

Selective pressures at a codon-level predict deleterious mutations in human disease genes

Deleterious mutations affecting biological function of proteins are constantly being rejected by purifying selection from the gene pool. The non-synonymous/synonymous substitution rate ratio (omega) is a measure of selective pressure on amino acid replacement mutations for protein-coding genes. Different methods have been developed in order to predict non-synonymous changes affecting gene function. However, none has considered the estimation of selective constraints acting on protein residues. Here, we have used codon-based maximum likelihood models in order to estimate the selective pressures on the individual amino acid residues of a well-known model protein: p53. We demonstrate that the number of residues under strong purifying selection in p53 is much higher than those that are strictly conserved during the evolution of the species. In agreement with theoretical expectations, residues that have been noted to be of structural relevance, or in direct association with DNA, were among those showing the highest signals of purifying selection. Conversely, those changing according to a neutral, or nearly neutral mode of evolution, were observed to be irrelevant for protein function. Finally, using more than 40 human disease genes, we demonstrate that residues evolving under strong selective pressures (omega<0.1) are significantly associated (p<0.01) with human disease. We hypothesize that non-synonymous change on amino acids showing omega<0.1 will most likely affect protein function. The application of this evolutionary prediction at a genomic scale will provide an a priori hypothesis of the phenotypic effect of non-synonymous coding single nucleotide polymorphisms (SNPs) in the human genome.     

Disorder predictors also predict backbone dynamics for a family of disordered proteins

Several algorithms have been developed that use amino acid sequences to predict whether or not a protein or a region of a protein is disordered. These algorithms make accurate predictions for disordered regions that are 30 amino acids or longer, but it is unclear whether the predictions can be directly related to the backbone dynamics of individual amino acid residues. The nuclear Overhauser effect between the amide nitrogen and hydrogen (NHNOE) provides an unambiguous measure of backbone dynamics at single residue resolution and is an excellent tool for characterizing the dynamic behavior of disordered proteins. In this report, we show that the NHNOE values for several members of a family of disordered proteins are highly correlated with the output from three popular algorithms used to predict disordered regions from amino acid sequence. This is the first test between an experimental measure of residue specific backbone dynamics and disorder predictions. The results suggest that some disorder predictors can accurately estimate the backbone dynamics of individual amino acids in a long disordered region.     

Defining the Erythrocyte Binding Domains of Plasmodium vivax Tryptophan Rich Antigen 33.5

Tryptophan-rich antigens play important role in host-parasite interaction. One of the Plasmodium vivax tryptophan-rich antigens called PvTRAg33.5 had earlier been shown to be predominantly of alpha helical in nature with multidomain structure, induced immune responses in humans, binds to host erythrocytes, and its sequence is highly conserved in the parasite population. In the present study, we divided this protein into three different parts i.e. N-terminal (amino acid position 24–106), middle (amino acid position 107–192), and C-terminal region (amino acid position 185–275) and determined the erythrocyte binding activity of these fragments. This binding activity was retained by the middle and C-terminal fragments covering 107 to 275 amino acid region of the PvTRAg33.5 protein. Eight non-overlapping peptides covering this 107 to 275 amino acid region were then synthesized and tested for their erythrocyte binding activity to further define the binding domains. Only two peptides, peptide P4 (at 171–191 amino acid position) and peptide P8 (at 255–275 amino acid position), were found to contain the erythrocyte binding activity. Competition assay revealed that each peptide recognizes its own erythrocyte receptor. These two peptides were found to be located on two parallel helices at one end of the protein in the modelled structure and could be exposed on its surface to form a suitable site for protein-protein interaction. Natural antibodies present in the sera of the P. vivax exposed individuals or the polyclonal rabbit antibodies against this protein were able to inhibit the erythrocyte binding activity of PvTRAg33.5, its fragments, and these two synthetic peptides P4 and P8. Further studies on receptor-ligand interaction might lead to the development of the therapeutic reagent.     

Relationship between structure and amino acid sequence of strongly twisted and coiled β-hairpins in globular proteins

β-Hairpins are widespread in proteins, and it is possible to find them both within β-sheets and separately. In this work, a comparative analysis of amino acid sequences of β-strands within strongly twisted β-hairpins from different structural protein subclasses has been conducted. Strongly twisted and coiled β-hairpin generates in the space a right double helix out of β-strands that are connected by a loop region (connections). The frequencies of amino acid residues on the internal (concave) and external (convex) surfaces of strongly twisted β-hairpins have been determined (220 β-hairpins from nonhomologous proteins were studied). The concave surface of these β-hairpins is mainly generated by hydrophobic residues, while the convex surface by hydrophilic residues; accordingly, the alternation of hydrophobic internal and hydrophilic external residues is observed in their amino acid sequences. Amino acid residues of glycine and alanine (especially in places of the largest twisting of the strands) were anomalously frequently found in internal positions of strongly twisted and coiled β-hairpins. It was established that internal positions never contain the proline residues, while external positions in the twisting region contain them in a relatively large amount. It was demonstrated that at least one amino acid residue in α_L- or ε-conformation is required for generation of relatively short (up to 7 amino acid residues) connection. As a rule, these positions are occupied by glycines. Thus, not only the alternation of hydrophobic and hydrophilic amino acid residues, but also the presence of one or two glycine residues in the connection region and the excess of glycines and alanines in the places of the largest strand twisting on the concave surface, as well as the presence of prolines on the convex surface, are required to generate a strongly twisted and coiled β-hairpin.     

Rise and fall of the delta globin gene

The complete nucleotide sequence of the gene phoE, which codes for the phosphate limitation inducible outer membrane pore protein of Escherichia coli K12 was established. The results show that PhoE protein is synthesized in a precursor form with a 21 amino acid residue amino-terminal extension. This peptide has the general characteristics of a signal sequence. The promoter region of phoE has no homlogy with the consensus sequence of E. coli promoter regions, but homologous sequences with the promoter region of phoA, the structural gene for alkaline phosphatase, were observed. The deduced amino acid sequence showed that the mature PhoE protein is composed of 330 amino acid residues with a calculated molecular weight of 36,782. A number of 81 charged amino acids was found scattered throughout the protein while no large stretches of hydrophobic amino acids were observed. Hydrophobicity and hydration profiles of PhoE protein showed five pronounced hydrophilic maxima which are all located in the region from the amino terminus to residue 212.When the deduced amino acid sequence of PhoE protein was compared with the established sequence of the OmpF pore protein, a number of 210 identical residues was found. Some aspects of the structure-function relationship of PhoE protein are discussed in view of the hydrophobicity and hydration profiles, and the homology between PhoE protein and OmpF protein.     

Belief and Behavior Aspects of the EAT-26: The Case of Schoolgirls in Belize

This study investigates components of eating attitudes in a sample of Belizean schoolgirls and argues for separate analysis of eating beliefs and eating behaviors using the EAT-26 in populations undergoing rapid cultural change. The EAT-26 was utilized in a novel manner, preserving the ethnographic and empirical distinction between belief and behavior components of eating attitudes. Participants included a sample of secondary schoolgirls (n = 80) undergoing acculturative stress. Participants reported more disordered eating beliefs than behaviors. Respondents having higher belief scores than behavior scores were more likely to prefer thinner body build and to be concerned about boys’ assessments of their bodies. Girls with higher behavior scores were less likely to report eating when hungry and stopping when full. In conclusion, discriminant validity was found between attitudinal and behavioral aspects of the EAT-26 as evidenced by face validity and patterns in predicting body image preference and desired weight change. Such a distinction has implications for assessing risk for disordered eating among populations undergoing acculturative stress. Among such populations, while behavioral symptoms might be absent or present in subclinical levels, disordered beliefs associated with psychological distress or potential precursors to eating-disordered behavior might be detected and should be investigated further.     

Characterization of kinectin, a kinesin-binding protein: primary sequence and N-terminal topogenic signal analysis.

Kinectin is a kinesin-binding protein (Toyoshima et al., 1992) that is required for kinesin-based motility (Kumar et al., 1995). A kinectin cDNA clone containing a 4.7-kilobase insert was isolated from an embryonic chick brain cDNA library by immunoscreening with a panel of monoclonal antibodies. The cDNA contained an open reading frame of 1364 amino acids encoding a protein of 156 kDa. A bacterially expressed product of the full length cDNA bound purified kinesin. Transient expression in CV-1 cells gave an endoplasmic reticulum distribution that depended upon the N-terminal domain. Analysis of the predicted amino acid sequence indicated a highly hydrophobic near N-terminal stretch of 28 amino acids and a large portion (326-1248) of predicted alpha helical coiled coils. The 30-kDa fragment containing the N-terminal hydrophobic region was produced by cell-free in vitro translation and found to assemble with canine pancreas rough microsomes. Cleavage of the N terminus was not observed confirming its role as a potential transmembrane domain. Thus, the kinectin cDNA encodes a cytoplasmic-oriented integral membrane protein that binds kinesin and is likely to be a coiled-coil dimer.     

Role of the buried glutamate in the alpha-helical coiled coil domain of the macrophage scavenger receptor

The macrophage scavenger receptor exhibits a pH-dependent conformational change around the carboxy-terminal half of the alpha-helical coiled coil domain, which has a representative amino acid sequence of a (defgabc)n heptad. We previously demonstrated that a peptide corresponding to this region formed a random coil structure at pH 7 and an alpha-helical coiled coil structure at pH 5 [Suzuki, K., Doi, T., Imanishi, T., Kodama, T., and Tanaka, T. (1997) Biochemistry 36, 15140-15146]. To determine the amino acid responsible for the conformational change, we prepared several peptides in which the acidic amino acids were replaced with neutral amino acids. Analyses of their structures by circular dichroism and sedimentation equilibrium gave the result that the presence of Glu242 at the d position was sufficient to induce the pH-dependent conformational change of the alpha-helical coiled coil domain. Furthermore, we substituted a Glu residue for the Ile residue at the d or a position of a de novo designed peptide (IEKKIEA)4, which forms a highly stable triple-stranded coiled coil. These peptides exhibited a pH-dependent conformational change similar to that of the scavenger receptor. Therefore, we conclude that a buried Glu residue in the hydrophobic core of a triple-stranded coiled coil has the potential to induce the pH-dependent conformational change. This finding makes it possible to elucidate the functions of natural proteins and to create a de novo protein designed to undergo a pH-dependent conformational change.     

The double helix coiled coil structure of murein lipoprotein from Escherichia coli.

A rod-like structure is proposed for the murein lipoprotein of Escherichia coli, built of two parallel unbroken α-helices arranged in a coiled coil of the same type as in the muscle protein tropomyosin. The amino acid sequence has the required regular pattern of hydrophobic amino acids at intervals of three and four residues and the secondary structure predicted from the sequence is 80% helical. A space-filling model confirms that the coiled coil model is stereochemically reasonable, and energy calculations for a series of coils with different radii suggest that the best structure is one with the helix axes 8.25 Å apart. Energyrefined atomic co-ordinates have been calculated which show that the hydrophobic side-chains form a series of close-packed unstrained contacts between the two helices along the entire length of the sequence. On the basis of this study the hexagonal membrane pore model and the segmented helix model proposed by others seem unlikely. The coiled coil has a strongly hydrophilic outer surface, suggesting that the protein has a watery environment within the E. coli cell envelope and is not strictly a membrane protein. Probably only the fatty acid portion of the lipoprotein penetrates into the lipid region of the outer membrane, so that the protein may act as a tie or a spacer between the lipid and the murein wall.