Crystal structure of Omp32, the anion-selective porin from Comamonas acidovorans, in complex with a periplasmic peptide at 2.1 A resolution

Motifs in outer membrane protein sequences: Applications for discrimination

Discriminating outer membrane proteins (OMPs) from other folding types of globular and membrane proteins is an important problem for predicting their secondary and tertiary structures and detecting outer membrane proteins from genomic sequences as well. In this work, we have systematically analyzed the distribution of amino acid residues in the sequences of globular and outer membrane proteins with several motifs, such as A*B, A**B, etc. We observed that the motifs E*L, A*K and L*E occur frequently in globular proteins while S*S, N*S and R*D predominantly occur in OMPs. We have devised a statistical method based on frequently occurring motifs in globular and OMPs and obtained an accuracy of 96% and 82% for correctly identifying OMPs and excluding globular proteins, respectively. Further, we noticed that the motifs of transmembrane helical (TMH) proteins are different from that of OMPs. While I*A, I*L and L*I prefer in TMH proteins S*S, N*S and N*N predominantly occur in OMPs. The information about the occurrence of A*B motifs in TMH and OMPs could discriminate them with an accuracy of 80% for excluding OMPs and 100% for identifying OMPs. The influence of protein size and structural class for discrimination is discussed.

Application of residue distribution along the sequence for discriminating outer membrane proteins

Discriminating outer membrane proteins from other folding types of globular and membrane proteins is an important problem both for detecting outer membrane proteins from genomic sequences and for the successful prediction of their secondary and tertiary structures. In this work, we have systematically analyzed the distribution of amino acid residues in the sequences of globular and outer membrane proteins. We observed that the occurrence of two neighboring aliphatic and polar residues is significantly higher in outer membrane proteins than in globular proteins. From the information about the dipeptide composition we have devised a statistical method for discriminating outer membrane proteins from other globular and membrane proteins. Our approach correctly picked up the outer membrane proteins with an accuracy of 95% for the training set of 337 proteins. On the other hand, our method has correctly excluded the globular proteins at an accuracy of 79% in a non-redundant dataset of 674 proteins. Furthermore, the present method is able to correctly exclude alpha-helical membrane proteins up to an accuracy of 87%. These accuracy levels are comparable to other methods in the literature. The influence of protein size and structural class for discrimination is discussed.

Ionic transport phenomena in nanofluidics: experimental and theoretical study of the exclusion-enrichment effect on a chip

In nanometer-sized apertures with charged surfaces, the extension of the electrical double layer results in the electrostatic exclusion of co-ions and enrichment in counterions, which affects the permselectivity of such structures. A modeling of this phenomenon is proposed and is compared with quantitative measurements of the ionic permeability change of a Pyrex nanoslit at low ionic strength. The comparison of experimental results with theoretical predictions justifies that electrostatic forces are the governing forces in nanofluidics.

Expression and refolding of Omp38 from Burkholderia pseudomallei and Burkholderia thailandensis, and its function as a diffusion porin

In the present paper, we describe cloning and expression of two outer membrane proteins, BpsOmp38 (from Burkholderia pseudomallei) and BthOmp38 (from Burkholderia thailandensis) lacking signal peptide sequences, using the pET23d(+) expression vector and Escherichia coli host strain Origami(DE3). The 38 kDa proteins, expressed as insoluble inclusion bodies, were purified, solubilized in 8 M urea, and then subjected to refolding experiments. As seen on SDS/PAGE, the 38 kDa band completely migrated to approximately 110 kDa when the purified monomeric proteins were refolded in a buffer system containing 10% (w/v) Zwittergent 3-14, together with a subsequent heating to 95 degrees C for 5 min. CD spectroscopy revealed that the 110 kDa proteins contained a predominant beta-sheet structure, which corresponded completely to the structure of the Omp38 proteins isolated from B. pseudomallei and B. thailandensis. Immunoblot analysis using anti-BpsOmp38 polyclonal antibodies and peptide mass analysis by MALDI-TOF (matrix-assisted laser-desorption ionization-time-of-flight) MS confirmed that the expressed proteins were BpsOmp38 and BthOmp38. The anti-BpsOmp38 antibodies considerably exhibited the inhibitory effects on the permeation of small sugars through the Omp38-reconstituted liposomes. A linear relation between relative permeability rates and M(r) of neutral sugars and charged antibiotics suggested strongly that the in vitro re-assembled Omp38 functioned fully as a diffusion porin.

Evaluation of methods for predicting the topology of beta-barrel outer membrane proteins and a consensus prediction method

Background

Prediction of the transmembrane strands and topology of β-barrel outer membrane proteins is of interest in current bioinformatics research. Several methods have been applied so far for this task, utilizing different algorithmic techniques and a number of freely available predictors exist. The methods can be grossly divided to those based on Hidden Markov Models (HMMs), on Neural Networks (NNs) and on Support Vector Machines (SVMs). In this work, we compare the different available methods for topology prediction of β-barrel outer membrane proteins. We evaluate their performance on a non-redundant dataset of 20 β-barrel outer membrane proteins of gram-negative bacteria, with structures known at atomic resolution. Also, we describe, for the first time, an effective way to combine the individual predictors, at will, to a single consensus prediction method.

Results

We assess the statistical significance of the performance of each prediction scheme and conclude that Hidden Markov Model based methods, HMM-B2TMR, ProfTMB and PRED-TMBB, are currently the best predictors, according to either the per-residue accuracy, the segments overlap measure (SOV) or the total number of proteins with correctly predicted topologies in the test set. Furthermore, we show that the available predictors perform better when only transmembrane β-barrel domains are used for prediction, rather than the precursor full-length sequences, even though the HMM-based predictors are not influenced significantly. The consensus prediction method performs significantly better than each individual available predictor, since it increases the accuracy up to 4% regarding SOV and up to 15% in correctly predicted topologies.

Conclusions

The consensus prediction method described in this work, optimizes the predicted topology with a dynamic programming algorithm and is implemented in a web-based application freely available to non-commercial users at .

Folding and assembly of beta-barrel membrane proteins

Beta-barrel membrane proteins occur in the outer membranes of Gram-negative bacteria, mitochondria and chloroplasts. The membrane-spanning sequences of beta-barrel membrane proteins are less hydrophobic than those of alpha-helical membrane proteins, which is probably the main reason why completely different folding and membrane assembly pathways have evolved for these two classes of membrane proteins. Some beta-barrel membrane proteins can be spontaneously refolded into lipid bilayer model membranes in vitro. They may also have this ability in vivo although lipid and protein chaperones likely assist with their assembly in appropriate target membranes. This review summarizes recent work on the thermodynamic stability and the mechanism of membrane insertion of beta-barrel membrane proteins in lipid model and biological membranes. How lipid compositions affect folding and assembly of beta-barrel membrane proteins is also reviewed. The stability of these proteins in membranes is not as large as previously thought (<10 kcal/mol) and is modulated by elastic forces of the lipid bilayer. Detailed kinetic studies indicate that beta-barrel membrane proteins fold in distinct steps with several intermediates that can be characterized in vitro. Formation of the barrel is synchronized with membrane insertion and all beta-hairpins insert simultaneously in a concerted pathway.

Synthesis and characterisation of a ligand that forms a stable tetrahedral intermediate in the active site of the Aureobacterium species (–)γ-lactamase

The crystal structure of a (-) gamma-lactamase from an Aureobacterium species showed a molecule bound covalently to the active site serine residue. This enzyme complex represented the first structure of a stably bound tetrahedral intermediate for an alpha/beta hydrolase fold enzyme. The structural elucidation of tetrahedral intermediates is important for the understanding of enzymatic mechanism, substrate recognition and enzyme inhibition. In this paper, we report the synthesis and subsequent characterisation of (3aR,7aS)-3a,4,7,7a-tetrahydrobenzo-[1,3]-dioxol-2-one (BD1), the molecule modelled into the Aureobacterium (-) gamma-lactamase active site. This molecule has been confirmed to be an inhibitor and to be displaced from the enzyme by the racemic gamma-lactam substrate.

Prediction of transmembrane regions of beta-barrel proteins using ANN- and SVM-based methods

This article describes a method developed for predicting transmembrane beta-barrel regions in membrane proteins using machine learning techniques: artificial neural network (ANN) and support vector machine (SVM). The ANN used in this study is a feed-forward neural network with a standard back-propagation training algorithm. The accuracy of the ANN-based method improved significantly, from 70.4% to 80.5%, when evolutionary information was added to a single sequence as a multiple sequence alignment obtained from PSI-BLAST. We have also developed an SVM-based method using a primary sequence as input and achieved an accuracy of 77.4%. The SVM model was modified by adding 36 physicochemical parameters to the amino acid sequence information. Finally, ANN- and SVM-based methods were combined to utilize the full potential of both techniques. The accuracy and Matthews correlation coefficient (MCC) value of SVM, ANN, and combined method are 78.5%, 80.5%, and 81.8%, and 0.55, 0.63, and 0.64, respectively. These methods were trained and tested on a nonredundant data set of 16 proteins, and performance was evaluated using "leave one out cross-validation" (LOOCV). Based on this study, we have developed a Web server, TBBPred, for predicting transmembrane beta-barrel regions in proteins (available at http://www.imtech.res.in/raghava/tbbpred).

A simple statistical method for discriminating outer membrane proteins with better accuracy

MOTIVATION

Discriminating outer membrane proteins from other folding types of globular and membrane proteins is an important task both for identifying outer membrane proteins from genomic sequences and for the successful prediction of their secondary and tertiary structures.

RESULTS

We have systematically analyzed the amino acid composition of globular proteins from different structural classes and outer membrane proteins. We found that the residues, Glu, His, Ile, Cys, Gln, Asn and Ser, show a significant difference between globular and outer membrane proteins. Based on this information, we have devised a statistical method for discriminating outer membrane proteins from other globular and membrane proteins. Our approach correctly picked up the outer membrane proteins with an accuracy of 89% for the training set of 337 proteins. On the other hand, our method has correctly excluded the globular proteins at an accuracy of 79% in a non-redundant dataset of 674 proteins. Furthermore, the present method is able to correctly exclude alpha-helical membrane proteins up to an accuracy of 80%. These accuracy levels are comparable to other methods in the literature, and this is a simple method, which could be used for dissecting outer membrane proteins from genomic sequences. The influence of protein size, structural class and specific residues for discrimination is discussed.

Separation methods in the analysis of protein membrane complexes

The separation of membrane protein complexes can be divided into two categories. One category, which is operated on a relatively large scale, aims to purify the membrane protein complex from membrane fractions while retaining its native form, mainly to characterize its nature. The other category aims to analyze the constituents of the membrane protein complex, usually on a small scale. Both of these face the difficulty of isolating the membrane protein complex without interference originating from the hydrophobic nature of membrane proteins or from the close association with membrane lipids. To overcome this difficulty, many methods have been employed. Crystallized membrane protein complexes are the most successful example of the former category. In these purification methods, special efforts are made in the steps prior to the column chromatography to enrich the target membrane protein complexes. Although there are specific aspects for each complex, the most popular method for isolating these membrane protein complexes is anion-exchange column chromatography, especially using weak anion-exchange columns. Another remarkable trend is metal affinity column chromatography, which purifies the membrane protein complex as an intact complex in one step. Such protein complexes contain subunit proteins which are genetically engineered so as to include multiple-histidine tags at carboxyl- or amino-termini. The key to these successes for multi-subunit complex isolation is the idea of keeping the expression at its physiological level, rather than overexpression. On the other hand, affinity purification using the Fv fragment, in which a Strep tag is genetically introduced, is ideal because this method does not introduce any change to the target protein. These purification methods supported by affinity interaction can be applied to minor membrane protein complexes in the membrane system. Isoelectric focusing (IEF) and blue native (BN) electrophoresis have also been employed to prepare membrane protein complexes. Generally, a combination of two or more chromatographic and/or electrophoretic methods is conducted to separate membrane protein complexes. IEF or BN electrophoresis followed by 2nd dimension electrophoresis serve as useful tools for analytical demand. However, some problems still exist in the 2D electrophoresis using IEF. To resolve such problems, many attempts have been made, e.g. introduction of new chaotropes, surfactants, reductants or supporting matrices. This review will focus in particular on two topics: the preparative methods that achieved purification of membrane protein complexes in the native (intact) form, and the analytical methods oriented to resolve the membrane proteins. The characteristics of these purification and analytical methods will be discussed along with plausible future developments taking into account the nature of membrane protein complexes.

BOMP: a program to predict integral beta-barrel outer membrane proteins encoded within genomes of Gram-negative bacteria

This work describes the development of a program that predicts whether or not a polypeptide sequence from a Gram-negative bacterium is an integral beta-barrel outer membrane protein. The program, called the beta-barrel Outer Membrane protein Predictor (BOMP), is based on two separate components to recognize integral beta-barrel proteins. The first component is a C-terminal pattern typical of many integral beta-barrel proteins. The second component calculates an integral beta-barrel score of the sequence based on the extent to which the sequence contains stretches of amino acids typical of transmembrane beta-strands. The precision of the predictions was found to be 80% with a recall of 88% when tested on the proteins with SwissProt annotated subcellular localization in Escherichia coli K 12 (788 sequences) and Salmonella typhimurium (366 sequences). When tested on the predicted proteome of E.coli, BOMP found 103 of a total of 4346 polypeptide sequences to be possible integral beta-barrel proteins. Of these, 36 were found by BLAST to lack similarity (E-value score < 1e-10) to proteins with annotated subcellular localization in SwissProt. BOMP predicted the content of integral beta-barrels per predicted proteome of 10 different bacteria to range from 1.8 to 3%. BOMP is available at http://www.bioinfo.no/tools/bomp.

The crystal structure of a (-) gamma-lactamase from an Aureobacterium species reveals a tetrahedral intermediate in the active site

The structure of the recombinant (-) gamma-lactamase from an Aureobacterium species has been solved at 1.73A resolution in the cubic space group F23 with unit cell parameters a=b=c=240.6A. The trimeric enzyme has an alpha/beta hydrolase fold and closely resembles the cofactor free haloperoxidases. The structure has been solved in complex with a covalently bound ligand originating from the host cell and also in the unligated form. The associated density in the former structure has been interpreted as the two-ring ligand (3aR,7aS)-3a,4,7,7a-tetrahydro-benzo [1,3] dioxol-2-one which forms a tetrahedral complex with OG of the catalytic Ser98. Soaks of these crystals with the industrial substrate gamma-lactam or its structural analogue, norcamphor, result in the displacement of the ligand from the enzyme active site, thereby allowing determination of the unligated structure. The presence of the ligand in the active site protects the enzyme from serine hydrolase inhibitors. Cyclic ethylene carbonate, the first ring of the ligand, was shown to be a substrate of the enzyme.

Structure of the OmpA-like domain of RmpM from Neisseria meningitidis

RmpM is a putative peptidoglycan binding protein from Neisseria meningitidis that has been shown to interact with integral outer membrane proteins such as porins and TonB-dependent transporters. Here we report the 1.9 A crystal structure of the C-terminal domain of RmpM. The 150-residue domain adopts a betaalphabetaalphabetabeta fold, as first identified in Bacillus subtilis chorismate mutase. The C-terminal RmpM domain is homologous to the periplasmic, C-terminal domain of Escherichia coli OmpA; these domains are thought to be responsible for non-covalent interactions with peptidoglycan. From the structure of the OmpA-like domain of RmpM, we suggest a putative peptidoglycan binding site and identify residues that may be essential for binding. Both the crystal structure and solution experiments indicate that RmpM may exist as a dimer. This would promote more efficient peptidoglycan binding, by allowing RmpM to interact simultaneously with two glycan chains through its C-terminal, OmpA-like binding domain, while its (structurally uncharacterized) N-terminal domain could stabilize oligomers of porins and TonB-dependent transporters in the outer membrane.

Porines bactériennes et sensibilité aux antibiotiques

In Gram negative bacteria, hydrophilic antibiotics such as beta-lactams and fluoroquinolons used the bacterial porin channel during their entry. The balance of the porin expression level and the molecular parameters which govern the molecule diffusion through the pore are important physiological points. Acquired in vivo beta-lactam resistance is often associated with porin loss, and recently clinical resistant strains synthetizing mutated porin have been described. These data highlight both the importance of the channel characteristics and the amino acid residues involved in the drug diffusion process. In addition, several mechanisms, including various repressors or activators as well as molecules inhibiting the pore synthesis or activity, argue for the complexity and plasticity of the bacterial control of porin function. All these aspects play a key role in both membrane permeability and efficiency of the antibiotic resistance process.

Functional reconstitution, gene isolation and topology modelling of porins from Burkholderia pseudomallei and Burkholderia thailandensis

The sequences for Omp38 from Burkholderia pseudomallei and Burkholderia thailandensis have been deposited in the DDBJ, EMBL, GenBank(R) and GSDB Nucleotide Sequence Databases under the accession numbers AY312416 and AY312417 respectively. The intracellular pathogen Burkholderia pseudomallei is the causative agent of tropical melioidosis, and Burkholderia thailandensis is a closely-related Gram-negative bacterium that does not cause serious disease. Like other bacteria, the major outer membrane (OM) porins of Burkholderia strains, Bps Omp38 and Bth Omp38 may have roles in antibiotic resistance and immunity. We purified both proteins and found them to be immunologically related, SDS-resistant, heat-sensitive trimers with M (r) of approx. 110000. In functional liposome-swelling assays, both proteins showed similar permeabilities for small sugar molecules, compatible with a pore diameter of between 1.2 and 1.6 nm. Secondary structure analysis by FTIR (Fourier-transform infrared) spectroscopy revealed almost identical spectra with predominantly beta-sheet structures, typical of bacterial porins. MALDI-TOF (matrix-assisted laser-desorption ionization-time of flight) MS and ESI/MS (electrospray ionization MS) analysis of each protein showed extensive sequence similarities to the OpcP1 porin from Burkholderia cepacia (later found to be 76.5% identical). Based on information from the incomplete B. pseudomallei genome-sequencing project, the genes encoding Omp38 were identified and amplified by PCR from B. pseudomallei and B. thailandensis genomic DNA. The nucleotide sequences are 99.7% identical, and the predicted processed proteins are 100% identical. Topology prediction and molecular modelling suggest that this newly-isolated and cloned porin is a 16-stranded beta-barrel and the external loops of the protein could be important determinants of the immune response to infection.

Neural network-based prediction of transmembrane β-strand segments in outer membrane proteins

Prediction of transmembrane beta-strands in outer membrane proteins (OMP) is one of the important problems in computational chemistry and biology. In this work, we propose a method based on neural networks for identifying the membrane-spanning beta-strands. We introduce the concept of "residue probability" for assigning residues in transmembrane beta-strand segments. The performance of our method is evaluated with single-residue accuracy, correlation, specificity, and sensitivity. Our predicted segments show a good agreement with experimental observations with an accuracy level of 73% solely from amino acid sequence information. Further, the predictive power of N- and C-terminal residues in each segments, number of segments in each protein, and the influence of cutoff probability for identifying membrane-spanning beta-strands will be discussed. We have developed a Web server for predicting the transmembrane beta-strands from the amino acid sequence, and the prediction results are available at http://psfs.cbrc.jp/tmbeta-net/.

On the Role of Structural Information in Remote Homology Detection and Sequence Alignment: New Methods Using Hybrid Sequence Profiles

Structural alignments often reveal relationships between proteins that cannot be detected using sequence alignment alone. However, profile search methods based entirely on structural alignments alone have not been found to be effective in finding remote homologs. Here, we explore the role of structural information in remote homolog detection and sequence alignment. To this end, we develop a series of hybrid multidimensional alignment profiles that combine sequence, secondary and tertiary structure information into hybrid profiles. Sequence-based profiles are profiles whose position-specific scoring matrix is derived from sequence alignment alone; structure-based profiles are those derived from multiple structure alignments. We compare pure sequence-based profiles to pure structure-based profiles, as well as to hybrid profiles that use combined sequence-and-structure-based profiles, where sequence-based profiles are used in loop/motif regions and structural information is used in core structural regions. All of the hybrid methods offer significant improvement over simple profile-to-profile alignment. We demonstrate that both sequence-based and structure-based profiles contribute to remote homology detection and alignment accuracy, and that each contains some unique information. We discuss the implications of these results for further improvements in amino acid sequence and structural analysis.

Molecular basis of bacterial outer membrane permeability revisited

Gram-negative bacteria characteristically are surrounded by an additional membrane layer, the outer membrane. Although outer membrane components often play important roles in the interaction of symbiotic or pathogenic bacteria with their host organisms, the major role of this membrane must usually be to serve as a permeability barrier to prevent the entry of noxious compounds and at the same time to allow the influx of nutrient molecules. This review summarizes the development in the field since our previous review (H. Nikaido and M. Vaara, Microbiol. Rev. 49:1-32, 1985) was published. With the discovery of protein channels, structural knowledge enables us to understand in molecular detail how porins, specific channels, TonB-linked receptors, and other proteins function. We are now beginning to see how the export of large proteins occurs across the outer membrane. With our knowledge of the lipopolysaccharide-phospholipid asymmetric bilayer of the outer membrane, we are finally beginning to understand how this bilayer can retard the entry of lipophilic compounds, owing to our increasing knowledge about the chemistry of lipopolysaccharide from diverse organisms and the way in which lipopolysaccharide structure is modified by environmental conditions.

Variation of amino acid properties in all-beta globular and outer membrane protein structures

Discriminating outer membrane (OM) proteins from globular proteins is an important task. The structural analysis of beta-strands dominating globular (all-beta) proteins and OM proteins provides useful insight to distinguish between them. In this work, we analyze the characteristic features of the 20 amino acid residues in all-beta and OM proteins. We set up numerical indices for several properties of amino acid residues, such as, conformational parameters, surrounding hydrophobicity, accessible surface area and reduction in accessibility, and inter-residue contacts. We found that all the aromatic residues prefer to be in beta-strands of both globular and OM proteins. The surrounding hydrophobicity of aromatic and non-polar amino acid residues in globular proteins is significantly higher than that of OM proteins. The residues Trp, Arg, Phe and Gln show a remarkable difference of reduction in accessibility between all-beta globular (betaG) and OM proteins. The positively charged residues, Lys and Arg in the membrane part of OM proteins have more number of contacts than globular proteins. Further, the behavior of the 20 amino acid residues in beta-strand segments of globular and OM proteins have been discussed. The parameters developed in this work can be used for identifying transmembrane beta-strands in OM proteins and for discriminating betaG proteins from OM proteins.

Multistep mechanism of chloride translocation in a strongly anion-selective porin channel

The strongly anion-selective porin channel Omp32 from the bacterium Delftia acidovorans differs from other unspecific porins by its pronounced selectivity for anions and its particularly small channel cross-section. Multinanosecond molecular dynamics simulations of chloride ion movement in this pore protein suggest that translocated anions interact intimately with the charges of a "basic ladder", whose dynamics lead the anions in a stepwise manner through the constriction zone of the channel. The ladder-steps comprise the central clustered arginine groups and flanking basic residues at its exoplasmic and periplasmic sides. The computed free energy profile of ion movement in and around the constriction zone shows a corresponding succession of free energy minima and barriers. A number of polar atoms from other amino acids contribute to the coordination of Cl(-) at certain sites and to its temporary immobilization in the channel. A special binding site occurs at the transition of the constriction zone to the periplasmic funnel, binding the chloride ion over significant lengths of time. The results from our MD study offer a possible explanation for the nonlinear conductance properties and unusual salt-dependent characteristics of Omp32 observed earlier in experimental measurements.

Crystal structure of Neisserial surface protein A (NspA), a conserved outer membrane protein with vaccine potential

The neisserial surface protein A (NspA) from Neisseria meningitidis is a promising vaccine candidate because it is highly conserved among meningococcal strains and induces bactericidal antibodies. NspA is a homolog of the Opa proteins, which mediate adhesion to host cells. Here, we present the crystal structure of NspA, determined to 2.55-A resolution. NspA forms an eight-stranded antiparallel beta-barrel. The four loops at the extracellular side of the NspA molecule form a long cleft, which contains mainly hydrophobic residues and harbors a detergent molecule, suggesting that the protein might function in the binding of hydrophobic ligands, such as lipids. In addition, the structure provides a starting point for structure-based vaccine design.

Role of surface-exposed loops of Haemophilus influenzae protein P2 in the mitogen-activated protein kinase cascade

The outer membrane of gram-negative bacteria contains several proteins, and some of these proteins, the porins, have numerous biological functions in the interaction with the host; porins are involved in the activation of signal transduction pathways and, in particular, in the activation of the Raf/MEK1-MEK2/mitogen-activated protein kinase (MAPK) cascade. The P2 porin is the most abundant outer membrane protein of Haemophilus influenzae type b. A three-dimensional structural model for P2 was constructed based on the crystal structures of Klebsiella pneumoniae OmpK36 and Escherichia coli PhoE and OmpF. The protein was readily assembled into the beta-barrel fold characteristic of porins, despite the low sequence identity with the template proteins. The model provides information on the structural features of P2 and insights relevant for prediction of domains corresponding to surface-exposed loops, which could be involved in the activation of signal transduction pathways. To identify the role of surface-exposed loops, a set of synthetic peptides were synthesized according to the proposed model and were assayed for MEK1-MEK2/MAPK pathway activation. Our results show that synthetic peptides corresponding to surface loops of protein P2 are able to activate the MEK1-MEK2/MAPK pathways like the entire protein, while peptides modeled on internal beta strands are unable to induce significant phosphorylation of the MEK1-MEK2/MAPK pathways. In particular, the peptides corresponding to loops L5 (Lys206 to Gly219), L6B (Ser239 to Lys253), and L7 (Thr280 to Lys287) activate, as the whole protein, essentially JNK and p38.

The core of the tetrameric mycobacterial porin MspA is an extremely stable beta-sheet domain

MspA is the major porin of Mycobacterium smegmatis mediating the exchange of hydrophilic solutes across the cell wall and is the prototype of a new family of tetrameric porins with a single central pore of 10 nm in length. Infrared and circular dichroism spectroscopy revealed that MspA consists mainly of antiparallel beta-strands organized in a coherent domain. Heating to 92 and 112 degrees C was required to dissociate the MspA tetramer and to unfold the beta-sheet domain in the monomer, respectively. The stability of the MspA tetramer exceeded the remarkable stability of the porins of Gram-negative bacteria for every condition tested and was not reduced in the presence of 2% SDS and at any pH from 2 to 14. These results indicated that the interactions between the MspA subunits are different from those in the porins of Gram-negative bacteria and are discussed in the light of a channel-forming beta-barrel as a core structure of MspA. Surprisingly, the channel activity of MspA in 2% SDS and 7.6 m urea at 50 degrees C was reduced 13- and 30-fold, respectively, although the MspA tetramer and the beta-sheet domain were stable under those conditions. Channel closure by conformational changes of extracellular loops under those conditions is discussed to explain these observations. This study presents the first experimental evidence that outer membrane proteins not only from Gram-negative bacteria but also from mycobacteria are beta-sheet proteins and demonstrates that MspA constitutes the most stable transmembrane channel protein known so far. Thus, MspA may be of special interest for biotechnological applications.

Understanding pH-dependent selectivity of alamethicin K18 channels by computer simulation

Alamethicin K18 is a covalently linked alamethicin dimer in which the glutamine residue at position 18 in each helix has been replaced by a lysine residue. As described in previous work, channels formed by this peptide show pH-dependent selectivity. The maximum anion selectivity of the putative octameric conducting state is obtained at pH 7 or lower. Inasmuch as no change in selectivity is seen between pH 7 and pH 3, and because protons are expected to be in equilibrium with the open state of the channel during a selectivity measurement, the channel is believed to be fully charged (i.e., all eight lysines protonated) at pH 7. In an effort to understand how such a highly charged channel structure is stable in membranes and why it is not more selective for anions, we have performed a number of computer simulations of the system. Molecular dynamics simulations of 10 ns each of the octameric bundle in a lipid bilayer environment are presented, with either zero, four, or eight lysines charged in the absence of salt, and with eight lysines charged in the presence of 0.5 M and 1 M KCl. When no salt is present and all lysines are charged, on average 1.9 Cl(-) ions are inside the channel and the channel significantly deforms. With 0.5 M KCl present, 2.9 Cl(-) ions are inside the channel. With 1 M KCl present, four Cl(-) ions are present and the channel maintains a regular structure. Poisson-Boltzmann calculations on models of the octameric channel also predict an average of 2-4 Cl(-) ions near the lysine residues as a function of ionic strength. These counterions lower the apparent charge of the channel, which may underlie the decrease in selectivity observed experimentally with increasing salt concentrations. We suggest that to increase the selectivity of Alm K18 channels, positive charges could be engineered in a narrower part of the channel.

The expression of outer membrane proteins for crystallization

The production of sufficient amounts of chemically and conformationally homogenous protein is a major requirement for successful crystallization and structure determination. With membrane proteins, this constitutes a particular problem because the membrane volume is limited and the organisms are usually very sensitive to changes in membrane properties brought about by massive protein insertion. Moreover, the extraction of membrane proteins from the membrane with detergents is generally a harsh treatment, which gives rise to conformational aberrations. A number of successful procedures for functional expression followed by purification are reviewed here together with nonfunctional expression into inclusion bodies and subsequent (re)folding to produce functional proteins. Most of the data are for prokaryotic outer membrane proteins, but the outer membrane proteins of eukaryotic organelles are also considered as they do show similar features.

Structural characterization of the fusobacterial non-specific porin FomA suggests a 14-stranded topology, unlike the classical porins

Native and recombinant FomA proteins were extracted by detergent from the cell envelopes of Fusobacterium nucleatum and Escherichia coli, and purified to near homogeneity by chromatography. Circular dichroism analysis revealed that the FomA protein consists predominantly of beta-sheets, in line with the previously proposed 16-stranded beta-barrel topology model. Results obtained by trypsin treatment of intact cells and cell envelopes of F. nucleatum, and from limited proteolysis of purified FomA protein, indicated that the N-terminal part of the FomA protein is not an integral part of the beta-barrel, but forms a periplasmic domain. Based on these results a new topology model is proposed for the FomA protein, where the C-terminal part forms a 14-stranded beta-barrel separate from the periplasmic N-terminal domain.

The structure of bacterial outer membrane proteins

Integral membrane proteins come in two types, alpha-helical and beta-barrel proteins. In both types, all hydrogen bonding donors and acceptors of the polypeptide backbone are completely compensated and buried while nonpolar side chains point to the membrane. The alpha-helical type is more abundant and occurs in cytoplasmic (or inner) membranes, whereas the beta-barrels are known from outer membranes of bacteria. The beta-barrel construction is described by the number of strands and the shear number, which is a measure for the inclination angle of the beta-strands against the barrel axis. The common right-handed beta-twist requires shear numbers slightly larger than the number of strands. Membrane protein beta-barrels contain between 8 and 22 beta-strands and have a simple topology that is probably enforced by the folding process. The smallest barrels form inverse micelles and work as enzymes or they bind to other macromolecules. The medium-range barrels form more or less specific pores for nutrient uptake, whereas the largest barrels occur in active Fe(2+) transporters. The beta-barrels are suitable objects for channel engineering, because the structures are simple and because many of these proteins can be produced into inclusion bodies and recovered therefrom in the exact native conformation.

Phylogenetic Evidence for Frequent Positive Selection and Recombination in the Meningococcal Surface Antigen PorB

Previous estimates of rates of synonymous (d(S)) and nonsynonymous (d(N)) substitution among Neisseria meningitidis gene sequences suggested that the surface loops of the variable outer membrane protein PorB were under only weak selection pressure from the host immune response. These findings were consistent with studies indicating that PorB variants were not always protective in immunological and microbiological assays and questioned the suitability of this protein as a vaccine component. PorB, which is expressed at high levels on the surface of the meningococcus, has been implicated in mechanisms of pathogenesis and has also been used as a typing target in epidemiological investigations. In this work, using more precise estimates of selection pressures and recombination rates, we have shown that some residues in the surface loops of PorB are under very strong positive selection, as great as that observed in human immunodeficiency virus-1 surface glycoproteins, whereas amino acids within the loops and the membrane-spanning regions of the protein are under purifying selection, presumably because of structural constraints. Congruence tests showed that recombination occurred at a rate that was not sufficient to erase all phylogenetic similarity and did not greatly bias selection analysis. Homology models of PorB structure indicated that many strongly selected sites encoded residues that were predicted to be exposed to host immune responses, implying that this protein is under strong immune selection and requires further examination as a potential vaccine candidate. These data show that phylogenetic inference can be used to complement immunological and biochemical data in the choice of vaccine candidates.

Identification and analysis of amino acid mutations in porin IB that mediate intermediate-level resistance to penicillin and tetracycline in Neisseria gonorrhoeae

PenB is the third resistance determinant in the stepwise acquisition of multiple resistance genes in chromosomally mediated resistant Neisseria gonorrhoeae (CMRNG). Alterations in por(IB), one of two alleles at the por locus that encodes the outer membrane protein porin IB (PIB), were recently reported to be responsible for the increased resistance to penicillin and tetracycline conferred by penB, but the specific mutations conferring antibiotic resistance were not identified experimentally. To determine which amino acids in PIB confer increased resistance, we transformed a recipient strain with chimeras of the por(IB) genes from strains FA1090 and FA140 (penB2). These studies revealed that two amino acid changes, G120D and A121D, were both necessary and sufficient to confer increased resistance to penicillin and tetracycline. Site-saturation and site-directed mutagenesis of Gly-120 and Ala-121 revealed that both a single mutation, G120K, and the double mutations G120R A121H and G120P A121P also conferred antibiotic resistance to the recipient strain. The identical mutations in PIA increased penicillin and tetracycline resistance either moderately or not at all. Analysis of por(IB) genes present in the GenBank database from 51 clinical isolates demonstrated that lysine and aspartate mutations at positions 120 and/or 121 also occur in nature. These studies demonstrate that charged amino acids at positions 120 and 121 in PIB are highly preferential for conferring resistance to penicillin and tetracycline in N. gonorrhoeae.

Electrostatic properties of the anion selective porin Omp32 from Delftia acidovorans and of the arginine cluster of bacterial porins

The functional properties of the anion-selective porin Omp32 from the bacterium Delftia acidovorans, formerly Comamonas acidovorans, are determined by the particularly narrow channel constriction and the electrostatic field inside and outside the pore. A cluster of arginines (Arg 38, Arg 75, and Arg 133) determines the electrostatic field close to the constriction zone. Stacked amino acids carrying charges are prone to drastic pK(a) shifts. However, optimized calculations of the titration behavior of charged groups, based on the finite-difference Poisson-Boltzmann technique, suggest that all the arginines are charged at physiological pH. Protonation of the clustered arginines is stabilized by one buried glutamate residue (Glu 58), which is strongly interacting with Arg 75 and Arg 38. This functional arrangement of three charged amino acid residues is of general significance because it is found in the constriction zones of all known 16-stranded porins from the alpha-, beta-, and gamma-proteobacteria.

Membrane protein complexes

A wide range of membrane protein structures have been published during the past two years. These have included proteins from both eucaryotic and heterologously overexpressed sources. Whereas some of these proteins were crystallised using conventional techniques, others employed the new methods of lipidic cubic phase crystallisation and antibody fragment co-crystallisation. These and other new approaches to the expression and crystallisation of membrane proteins are accelerating structural studies of membrane protein complexes.

Toward genomic identification of beta-barrel membrane proteins: composition and architecture of known structures

The amino acid composition and architecture of all beta-barrel membrane proteins of known three-dimensional structure have been examined to generate information that will be useful in identifying beta-barrels in genome databases. The database consists of 15 nonredundant structures, including several novel, recent structures. Known structures include monomeric, dimeric, and trimeric beta-barrels with between 8 and 22 membrane-spanning beta-strands each. For this analysis the membrane-interacting surfaces of the beta-barrels were identified with an experimentally derived, whole-residue hydrophobicity scale, and then the barrels were aligned normal to the bilayer and the position of the bilayer midplane was determined for each protein from the hydrophobicity profile. The abundance of each amino acid, relative to the genomic abundance, was calculated for the barrel exterior and interior. The architecture and diversity of known beta-barrels was also examined. For example, the distribution of rise-per-residue values perpendicular to the bilayer plane was found to be 2.7 +/- 0.25 A per residue, or about 10 +/- 1 residues across the membrane. Also, as noted by other authors, nearly every known membrane-spanning beta-barrel strand was found to have a short loop of seven residues or less connecting it to at least one adjacent strand. Using this information we have begun to generate rapid screening algorithms for the identification of beta-barrel membrane proteins in genomic databases. Application of one algorithm to the genomes of Escherichia coli and Pseudomonas aeruginosa confirms its ability to identify beta-barrels, and reveals dozens of unidentified open reading frames that potentially code for beta-barrel outer membrane proteins.

Continuum electrostatic analysis of irregular ionization and proton allocation in proteins

Irregular (nonsigmoidal) ionization behavior of titratable groups in proteins is analyzed theoretically, using a computational algorithm designed to count explicitly for tautomers of titratable groups and different locations of polar hydrogens. On the basis of calculations for different model systems (acid-acid, base-base, acid-base pairs, and cluster of three strongly interacting groups), it is demonstrated that the pK values, extracted from nonsigmoidal titration curves by fitting to a sum of Henderson-Hasselbalch equations, do not describe the ionization equilibrium correctly. The conditions for observation of irregular titration curves are derived analytically for the case of arbitrary couple of interacting ionizable groups. A possible relation between irregularly shaped titration curves and tautomerization is also illustrated. The protonation-deprotonation equilibrium of Asp76 in ribonuclease T1 is shown to be coupled to dipole reorientation of a water molecule bound at the protein-solvent interface. This finding provides a new interpretation of the experimentally observed chemical shift of this residue.

The role of the TolC family in protein transport and multidrug efflux. From stereochemical certainty to mechanistic hypothesis

Gram-negative bacteria are enveloped by a system of two membranes, and they use specialized multicomponent, energy-driven pumps to transport molecules directly across this double-layered partition from the cell interior to the extra-cellular environment. One component of these pumps is embedded in the outer-membrane, and the paradigm for its structure and function is the TolC protein from Escherichia coli. A common component of a wide variety of efflux pumps, TolC and its homologues are involved in the export of chemically diverse molecules ranging from large protein toxins, such as alpha-hemolysin, to small toxic compounds, such as antibiotics. TolC family members thus play important roles in conferring pathogenic bacteria with both virulence and multidrug resistance. These pumps assemble reversibly in a transient process that brings together TolC or its homologue, an inner-membrane-associated periplasmic component, an integral inner-membrane translocase and the substrate itself. TolC can associate in this fashion with a variety of different partners to participate in the transport of diverse substrates. We review here the structure and function of TolC and the other components of the efflux/transport pump.

Prediction of the transmembrane regions of beta-barrel membrane proteins with a neural network-based predictor

A method based on neural networks is trained and tested on a nonredundant set of beta-barrel membrane proteins known at atomic resolution with a jackknife procedure. The method predicts the topography of transmembrane beta strands with residue accuracy as high as 78% when evolutionary information is used as input to the network. Of the transmembrane beta-strands included in the training set, 93% are correctly assigned. The predictor includes an algorithm of model optimization, based on dynamic programming, that correctly models eight out of the 11 proteins present in the training/testing set. In addition, protein topology is assigned on the basis of the location of the longest loops in the models. We propose this as a general method to fill the gap of the prediction of beta-barrel membrane proteins.