Orchestrating bacterial cell morphogenesis

MreB proteins are bacterial homologues of actin that directly determine cell shape and are involved in a range of other cellular processes in non-spherical bacteria. Like F-actin in eukaryotes, MreBs self-assemble into dynamic filamentous structures that are essential for cell viability. Recent studies have demonstrated that the MreB cytoskeletal scaffold governs shape determination by controlling functions related to the bacterial cell wall (probably by recruiting and directing peptidoglycan-synthesizing and modifying proteins). Here I consider general implications for bacterial morphogenesis, and the basis for differences in wall expansion and cylindrical cell shape, based on recent studies aimed to determine the role of MreBs in bacteria with different modes of growth.     

The dispanins: a novel gene family of ancient origin that contains 14 human members

The Interferon induced transmembrane proteins (IFITM) are a family of transmembrane proteins that is known to inhibit cell invasion of viruses such as HIV-1 and influenza. We show that the IFITM genes are a subfamily in a larger family of transmembrane (TM) proteins that we call Dispanins, which refers to a common 2TM structure. We mined the Dispanins in 36 eukaryotic species, covering all major eukaryotic groups, and investigated their evolutionary history using Bayesian and maximum likelihood approaches to infer a phylogenetic tree. We identified ten human genes that together with the known IFITM genes form the Dispanin family. We show that the Dispanins first emerged in eukaryotes in a common ancestor of choanoflagellates and metazoa, and that the family later expanded in vertebrates where it forms four subfamilies (A-D). Interestingly, we also find that the family is found in several different phyla of bacteria and propose that it was horizontally transferred to eukaryotes from bacteria in the common ancestor of choanoflagellates and metazoa. The bacterial and eukaryotic sequences have a considerably conserved protein structure. In conclusion, we introduce a novel family, the Dispanins, together with a nomenclature based on the evolutionary origin.     

NMR structure and dynamics of a designed water-soluble transmembrane domain of nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is an important therapeutic target for a wide range of pathophysiological conditions, for which rational drug designs often require receptor structures at atomic resolution. Recent proof-of-concept studies demonstrated a water-solubilization approach to structure determination of membrane proteins by NMR (Slovic et al., PNAS, 101: 1828-1833, 2004; Ma et al., PNAS, 105: 16537-42, 2008). We report here the computational design and experimental characterization of WSA, a water-soluble protein with ~83% sequence identity to the transmembrane (TM) domain of the nAChR α1 subunit. Although the design was based on a low-resolution structural template, the resulting high-resolution NMR structure agrees remarkably well with the recent crystal structure of the TM domains of the bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC), demonstrating the robustness and general applicability of the approach. NMR T(2) dispersion measurements showed that the TM2 domain of the designed protein was dynamic, undergoing conformational exchange on the NMR timescale. Photoaffinity labeling with isoflurane and propofol photolabels identified a common binding site in the immediate proximity of the anesthetic binding site found in the crystal structure of the anesthetic-GLIC complex. Our results illustrate the usefulness of high-resolution NMR analyses of water-solubilized channel proteins for the discovery of potential drug binding sites.     

An automatic method for predicting transmembrane protein structures using cryo-EM and evolutionary data

The transmembrane (TM) domains of many integral membrane proteins are composed of alpha-helix bundles. Structure determination at high resolution (<4 A) of TM domains is still exceedingly difficult experimentally. Hence, some TM-protein structures have only been solved at intermediate (5-10 A) or low (>10 A) resolutions using, for example, cryo-electron microscopy (cryo-EM). These structures reveal the packing arrangement of the TM domain, but cannot be used to determine the positions of individual amino acids. The observation that typically, the lipid-exposed faces of TM proteins are evolutionarily more variable and less charged than their core provides a simple rule for orienting their constituent helices. Based on this rule, we developed score functions and automated methods for orienting TM helices, for which locations and tilt angles have been determined using, e.g., cryo-EM data. The method was parameterized with the aim of retrieving the native structure of bacteriorhodopsin among near- and far-from-native templates. It was then tested on proteins that differ from bacteriorhodopsin in their sequences, architectures, and functions, such as the acetylcholine receptor and rhodopsin. The predicted structures were within 1.5-3.5 A from the native state in all cases. We conclude that the computational method can be used in conjunction with cryo-EM data to obtain approximate model structures of TM domains of proteins for which a sufficiently heterogeneous set of homologs is available. We also show that in those proteins in which relatively short loops connect neighboring helices, the scoring functions can discriminate between near- and far-from-native conformations even without the constraints imposed on helix locations and tilt angles that are derived from cryo-EM.     

Structural/functional homology between the bacterial and eukaryotic cytoskeletons

Structural proteins are now known to be as necessary for controlling cell division and cell shape in prokaryotes as they are in eukaryotes. Bacterial ParM and MreB not only have atomic structures that resemble eukaryotic actin and form similar filaments, but they are also equivalent in function: the assembly of ParM drives intracellular motility and MreB maintains the shape of the cell. FtsZ resembles tubulin in structure and in its dynamic assembly, and is similarly controlled by accessory proteins. Bacterial MinD and eukaryotic dynamin appear to have similar functions in membrane control. In dividing eukaryotic organelles of bacterial origin, bacterial and eukaryotic proteins work together.     

Strategies for dealing with conformational sampling in structural calculations of flexible or kinked transmembrane peptides.

Peptides corresponding to transmembrane (TM) segments from membrane proteins provide a potential route for the determination of membrane protein structure. We have determined that 2 functionally critical TM segments from the mammalian Na+/H+ exchanger display well converged structure in regions separated by break points. The flexibility of these break points results in conformational sampling in solution. A brief review of available NMR structures of helical membrane proteins demonstrates that there are a number of published structures showing similar properties. Such flexibility is likely indicative of kinks in the full-length protein. This minireview focuses on methods and protocols for NMR structure calculation and analysis of peptide structures under conditions of conformational sampling. The methods outlined allow the identification and analysis of structured peptides containing break points owing to conformational sampling and the differentiation between oligomerization and ensemble-averaged observation of multiple peptide conformations.     

A putative prokaryote voltage-gated Ca(2+) channel with only one 6TM motif per subunit

Until now, voltage-gated Ca(2+) channel proteins have been found only in eukaryotes. Here we report that a gene recently discovered in the eubacterium Bacillus halodurans codes for a protein closely related to eukaryotic Ca(2+) channels, but that has only one 6-transmembrane-segement (6TM) motif, instead of four, in its pore-forming subunit. This is supported by the comparison of consensus sequences, which, along with the patterns of residue conservation, indicates a similar structure in the membrane to voltage-gated K(+) channels. From this we hypothesize that Ca(2+) channels originally evolved in bacteria, and that the specific eubacteria protein highlighted here is an ideal candidate for structure determination efforts.     

Error tolerant NMR backbone resonance assignment and automated structure generation

Error tolerant backbone resonance assignment is the cornerstone of the NMR structure determination process. Although a variety of assignment approaches have been developed, none works sufficiently well on noisy fully automatically picked peaks to enable the subsequent automatic structure determination steps. We have designed an integer linear programming (ILP) based assignment system (IPASS) that has enabled fully automatic protein structure determination for four test proteins. IPASS employs probabilistic spin system typing based on chemical shifts and secondary structure predictions. Furthermore, IPASS extracts connectivity information from the inter-residue information and the (automatically picked) (15)N-edited NOESY peaks which are then used to fix reliable fragments. When applied to automatically picked peaks for real proteins, IPASS achieves an average precision and recall of 82% and 63%, respectively. In contrast, the next best method, MARS, achieves an average precision and recall of 77% and 36%, respectively. The assignments generated by IPASS are then fed into our protein structure calculation system, FALCON-NMR, to determine the 3D structures without human intervention. The final models have backbone RMSDs of 1.25?, 0.88?, 1.49?, and 0.67? to the reference native structures for proteins TM1112, CASKIN, VRAR, and HACS1, respectively. The web server is publicly available at http://monod.uwaterloo.ca/nmr/ipass.     

Aspartate dehydrogenase,a novel enzyme identified from structural and functional studies of TM1643

The open reading frame TM1643 of Thermotoga maritima belongs to a large family of proteins, with homologues in bacteria, archaea, and eukaryotes. TM1643 is found in an operon with two other genes that encode enzymes involved in the biosynthesis of NAD. In several bacteria, the gene in the position occupied by TM1643 encodes an aspartate oxidase (NadB), which synthesizes iminoaspartate as a substrate for NadA, the next enzyme in the pathway. The amino acid sequence of TM1643 does not share any recognizable homology with aspartate oxidase or with other proteins of known functions or structures. To help define the biological functions of TM1643, we determined its crystal structure at 2.6A resolution and performed a series of screens for enzymatic function. The structure reveals the presence of an N-terminal Rossmann fold domain with a bound NAD(+) cofactor and a C-terminal alpha+beta domain. The structural information suggests that TM1643 may be a dehydrogenase and the active site of the enzyme is located at the interface between the two domains. The enzymatic characterization of TM1643 revealed that it possesses NAD or NADP-dependent dehydrogenase activity toward l-aspartate but no aspartate oxidase activity. The product of the aspartate dehydrogenase activity is also iminoaspartate. Therefore, our studies demonstrate that two different enzymes, an oxidase and a dehydrogenase, may have evolved to catalyze the first step of NAD biosynthesis in prokaryotes. TM1643 establishes a new class of amino acid dehydrogenases.     

Modeling Transmembrane Domain Dimers/Trimers of Plexin Receptors: Implications for Mechanisms of Signal Transmission across the Membrane

Single-pass transmembrane (TM) receptors transmit signals across lipid bilayers by helix association or by configurational changes within preformed dimers. The structure determination for such TM regions is challenging and has mostly been accomplished by NMR spectroscopy. Recently, the computational prediction of TM dimer structures is becoming recognized for providing models, including alternate conformational states, which are important for receptor regulation. Here we pursued a strategy to predict helix oligomers that is based on packing considerations (using the PREDDIMER webserver) and is followed by a refinement of structures, utilizing microsecond all-atom molecular dynamics simulations. We applied this method to plexin TM receptors, a family of 9 human proteins, involved in the regulation of cell guidance and motility. The predicted models show that, overall, the preferences identified by PREDDIMER are preserved in the unrestrained simulations and that TM structures are likely to be diverse across the plexin family. Plexin-B1 and –B3 TM helices are regular and tend to associate, whereas plexin-A1, -A2, –A3, -A4, -C1 and –D1 contain sequence elements, such as poly-Glycine or aromatic residues that distort helix conformation and association. Plexin-B2 does not form stable dimers due to the presence of TM prolines. No experimental structural information on the TM region is available for these proteins, except for plexin-C1 dimeric and plexin-B1 – trimeric structures inferred from X-ray crystal structures of the intracellular regions. Plexin-B1 TM trimers utilize Ser and Thr sidechains for interhelical contacts. We also modeled the juxta-membrane (JM) region of plexin-C1 and plexin-B1 and show that it synergizes with the TM structures. The structure and dynamics of the JM region and TM-JM junction provide determinants for the distance and distribution of the intracellular domains, and for their binding partners relative to the membrane. The structures suggest experimental tests and will be useful for the interpretation of future studies.     

Compaction of a Prokaryotic Signal-Anchor Transmembrane Domain Begins within the Ribosome Tunnel and Is Stabilized by SRP during Targeting

Cotranslational targeting of membrane proteins is mediated by the universally conserved signal recognition particle (SRP). In eukaryotes, SRP attenuates translation during targeting; however, in prokaryotes, a simplified SRP is believed to carry out targeting during continuing translation. Here, we show a detailed stepwise analysis of the targeting of subunit c of the F₀ component of the bacterial ATP synthase (F₀c) to the inner membrane. We show that the first transmembrane (TM) signal-anchor domain of F₀c forms a compacted structure within the distal portion of the ribosome tunnel. This structure is formed just prior to the interaction with SRP. In the absence of SRP this structure is lost as the TM domain exits the tunnel; however in the presence of SRP it is stabilized. Our results suggest differences in early protein folding of substrates for prokaryotic SRP‐dependent membrane protein targeting pathways, from that of eukaryotic SRP targeting. These results imply that early TM domain recognition by targeting factors acts to ensure that the efficiency of membrane targeting is maintained.     

Histidine phosphorylation in biological systems

Histidine phosphorylation is important in prokaryotes and occurs to the extent of 6% of total phosphorylation in eukaryotes. Nevertheless phosphohistidine residues are not normally observed in proteins due to rapid hydrolysis of the phosphoryl group under acidic conditions. Many rapid processes employ phosphohistidines, including the bacterial phosphoenolpyruvate:sugar phosphotransferase system (PTS), the bacterial two-component systems and reactions catalyzed by enzymes such as nucleoside diphosphate kinase and succinyl-CoA synthetase. In the PTS, the NMR structure of the phosphohistidine moiety of the phosphohistidine-containing protein was determined but no X-ray structures of phosphohistidine forms of PTS proteins have been elucidated. There have been crystal structures of a few phosphohistidine-containing proteins determined: nucleoside diphosphate kinase, succinyl-CoA synthetase, a cofactor-dependent phosphoglycerate mutase and the protein PAE2307 from the hyperthermophilic archaeon Pyrobaculum aerophilum. A common theme for these stable phosphohistidines is the occurrence of ion-pair hydrogen bonds (salt bridges) involving the non-phosphorylated nitrogen atom of the histidine imidazole ring with an acidic amino acid side chain.     

In situ proteolysis for protein crystallization and structure determination

We tested the general applicability of in situ proteolysis to form protein crystals suitable for structure determination by adding a protease (chymotrypsin or trypsin) digestion step to crystallization trials of 55 bacterial and 14 human proteins that had proven recalcitrant to our best efforts at crystallization or structure determination. This is a work in progress; so far we determined structures of 9 bacterial proteins and the human aminoimidazole ribonucleotide synthetase (AIRS) domain.     

Progress in structure prediction of alpha-helical membrane proteins

Transmembrane (TM) proteins comprise 20-30% of the genome but, because of experimental difficulties, they represent less than 1% of the Protein Data Bank. The dearth of membrane protein structures makes computational prediction a potentially important means of obtaining novel structures. Recent advances in computational methods have been combined with experimental data to constrain the modeling of three-dimensional structures. Furthermore, threading and ab initio modeling approaches that were effective for soluble proteins have been applied to TM domains. Surprisingly, experimental structures, proteomic analyses and bioinformatics have revealed unexpected architectures that counter long-held views on TM protein structure and stability. Future computational and experimental studies aimed at understanding the thermodynamic and evolutionary bases of these architectural details will greatly enhance predictive capabilities.     

Comparing function and structure between entire proteomes

More than 30 organisms have been sequenced entirely. Here, we applied a variety of simple bioinformatics tools to analyze 29 proteomes for representatives from all three kingdoms: eukaryotes, prokaryotes, and archaebacteria. We confirmed that eukaryotes have relatively more long proteins than prokaryotes and archaes, and that the overall amino acid composition is similar among the three. We predicted that approximately 15%-30% of all proteins contained transmembrane helices. We could not find a correlation between the content of membrane proteins and the complexity of the organism. In particular, we did not find significantly higher percentages of helical membrane proteins in eukaryotes than in prokaryotes or archae. However, we found more proteins with seven transmembrane helices in eukaryotes and more with six and 12 transmembrane helices in prokaryotes. We found twice as many coiled-coil proteins in eukaryotes (10%) as in prokaryotes and archaes (4%-5%), and we predicted approximately 15%-25% of all proteins to be secreted by most eukaryotes and prokaryotes. Every tenth protein had no known homolog in current databases, and 30%-40% of the proteins fell into structural families with >100 members. A classification by cellular function verified that eukaryotes have a higher proportion of proteins for communication with the environment. Finally, we found at least one homolog of experimentally known structure for approximately 20%-45% of all proteins; the regions with structural homology covered 20%-30% of all residues. These numbers may or may not suggest that there are 1200-2600 folds in the universe of protein structures. All predictions are available at http://cubic.bioc.columbia.edu/genomes.     

DNA replication in the archaea.

The archaeal DNA replication machinery bears striking similarity to that of eukaryotes and is clearly distinct from the bacterial apparatus. In recent years, considerable advances have been made in understanding the biochemistry of the archaeal replication proteins. Furthermore, a number of structures have now been obtained for individual components and higher-order assemblies of archaeal replication factors, yielding important insights into the mechanisms of DNA replication in both archaea and eukaryotes.     

DNA Replication in the Archaea

The archaeal DNA replication machinery bears striking similarity to that of eukaryotes and is clearly distinct from the bacterial apparatus. In recent years, considerable advances have been made in understanding the biochemistry of the archaeal replication proteins. Furthermore, a number of structures have now been obtained for individual components and higher-order assemblies of archaeal replication factors, yielding important insights into the mechanisms of DNA replication in both archaea and eukaryotes.     

NMR structure of the transmembrane domain of the n-acetylcholine receptor β2 subunit

Nicotinic acetylcholine receptors (nAChRs) are involved in fast synaptic transmission in the central and peripheral nervous system. Among the many different types of subunits in nAChRs, the β2 subunit often combines with the α4 subunit to form α4β2 pentameric channels, the most abundant subtype of nAChRs in the brain. Besides computational predictions, there is limited experimental data available on the structure of the β2 subunit. Using high-resolution NMR spectroscopy, we solved the structure of the entire transmembrane domain (TM1234) of the β2 subunit. We found that TM1234 formed a four-helix bundle in the absence of the extracellular and intracellular domains. The structure exhibited many similarities to those previously determined for the Torpedo nAChR and the bacterial ion channel GLIC. We also assessed the influence of the fourth transmembrane helix (TM4) on the rest of the domain. Although secondary structures and tertiary arrangements were similar, the addition of TM4 caused dramatic changes in TM3 dynamics and subtle changes in TM1 and TM2. Taken together, this study suggests that the structures of the transmembrane domains of these proteins are largely shaped by determinants inherent in their sequence, but their dynamics may be sensitive to modulation by tertiary and quaternary contacts.     

Prioritizing targets for structural biology through the lens of proteomics: The archaeal protein TGAM_1934 from Thermococcus gammatolerans

ORFans are hypothetical proteins lacking any significant sequence similarity with other proteins. Here, we highlighted by quantitative proteomics the TGAM_1934 ORFan from the hyperradioresistant Thermococcus gammatolerans archaeon as one of the most abundant hypothetical proteins. This protein has been selected as a priority target for structure determination on the basis of its abundance in three cellular conditions. Its solution structure has been determined using multidimensional heteronuclear NMR spectroscopy. TGAM_1934 displays an original fold, although sharing some similarities with the 3D structure of the bacterial ortholog of frataxin, CyaY, a protein conserved in bacteria and eukaryotes and involved in iron–sulfur cluster biogenesis. These results highlight the potential of structural proteomics in prioritizing ORFan targets for structure determination based on quantitative proteomics data. The proteomic data and structure coordinates have been deposited to the ProteomeXchange with identifier PXD000402 ( http://proteomecentral.proteomexchange.org/dataset/PXD000402 ) and Protein Data Bank under the accession number 2mcf, respectively.     

Coverage of whole proteome by structural genomics observed through protein homology modeling database

We have been developing FAMSBASE, a protein homology-modeling database of whole ORFs predicted from genome sequences. The latest update of FAMSBASE (), which is based on the protein three-dimensional (3D) structures released by November 2003, contains modeled 3D structures for 368,724 open reading frames (ORFs) derived from genomes of 276 species, namely 17 archaebacterial, 130 eubacterial, 18 eukaryotic and 111 phage genomes. Those 276 genomes are predicted to have 734,193 ORFs in total and the current FAMSBASE contains protein 3D structure of approximately 50% of the ORF products. However, cases that a modeled 3D structure covers the whole part of an ORF product are rare. When portion of an ORF with 3D structure is compared in three kingdoms of life, in archaebacteria and eubacteria, approximately 60% of the ORFs have modeled 3D structures covering almost the entire amino acid sequences, however, the percentage falls to about 30% in eukaryotes. When annual differences in the number of ORFs with modeled 3D structure are calculated, the fraction of modeled 3D structures of soluble protein for archaebacteria is increased by 5%, and that for eubacteria by 7% in the last 3 years. Assuming that this rate would be maintained and that determination of 3D structures for predicted disordered regions is unattainable, whole soluble protein model structures of prokaryotes without the putative disordered regions will be in hand within 15 years. For eukaryotic proteins, they will be in hand within 25 years. The 3D structures we will have at those times are not the 3D structure of the entire proteins encoded in single ORFs, but the 3D structures of separate structural domains. Measuring or predicting spatial arrangements of structural domains in an ORF will then be a coming issue of structural genomics.