Characterization of the outer membrane protein OprF of Pseudomonas aeruginosa in a lipopolysaccharide membrane by computer simulation

The N-terminal domain of outer membrane protein OprF of Pseudomonas aeruginosa forms a membrane spanning eight-stranded antiparallel beta-barrel domain that folds into a membrane channel with low conductance. The structure of this protein has been modeled after the crystal structure of the homologous protein OmpA of Escherichia coli. A number of molecular dynamics simulations have been carried out for the homology modeled structure of OprF in an explicit molecular model for the rough lipopolysaccharide (LPS) outer membrane of P. aeruginosa. The structural stability of the outer membrane model as a result of the strong electrostatic interactions compared with simple lipid bilayers is restricting both the conformational flexibility and the lateral diffusion of the porin in the membrane. Constricting side-chain interactions within the pore are similar to those found in reported simulations of the protein in a solvated lipid bilayer membrane. Because of the strong interactions between the loop regions of OprF and functional groups in the saccharide core of the LPS, the entrance to the channel from the extracellular space is widened compared with the lipid bilayer simulations in which the loops are extruding in the solvent. The specific electrostatic signature of the LPS membrane, which results in a net intrinsic dipole across the membrane, is found to be altered by the presence of OprF, resulting in a small electrically positive patch at the position of the channel.     

OmpA: a pore or not a pore? Simulation and modeling studies

The bacterial outer membrane protein OmpA is composed of an N-terminal 171-residue beta-barrel domain (OmpA(171)) that spans the bilayer and a periplasmic, C-terminal domain of unknown structure. OmpA has been suggested to primarily serve a structural role, as no continuous pore through the center of the barrel can be discerned in the crystal structure of OmpA(171). However, several groups have recorded ionic conductances for bilayer-reconstituted OmpA(171). To resolve this apparent paradox we have used molecular dynamics (MD) simulations on OmpA(171) to explore the conformational dynamics of the protein, in particular the possibility of transient formation of a central pore. A total of 19 ns of MD simulations of OmpA(171) have been run, and the results were analyzed in terms of 1) comparative behavior of OmpA(171) in different bilayer and bilayer-mimetic environments, 2) solvation states of OmpA(171), and 3) pore characteristics in different MD simulations. Significant mobility was observed for residues and water molecules within the beta-barrel. A simulation in which putative gate region side chains of the barrel interior were held in a non-native conformation led to an open pore, with a predicted conductance similar to experimental measurements. The OmpA(171) pore has been shown to be somewhat more dynamic than suggested by the crystal structure. A gating mechanism is proposed to explain its documented channel properties, involving a flickering isomerization of Arg138, forming alternate salt bridges with Glu52 (closed state) and Glu128 (open state).     

A multidomain outer membrane protein from Pasteurella multocida: modelling and simulation studies of PmOmpA

PmOmpA is a two-domain outer membrane protein from Pasteurella multocida. The N-terminal domain of PmOmpA is a homologue of the transmembrane beta-barrel domain of OmpA from Escherichia coli, whilst the C-terminal domain of PmOmpA is a homologue of the extra-membrane Neisseria meningitidis RmpM C-terminal domain. This enables a model of a complete two domain PmOmpA to be constructed and its conformational dynamics explored via MD simulations of the protein embedded within two different phospholipid bilayers (DMPC and DMPE). The conformational stability of the transmembrane beta-barrel is similar to that of a homology model of OprF from Pseudomonas aeruginosa in bilayer simulations. There is a degree of water penetration into the interior of the beta-barrel, suggestive of a possible transmembrane pore. Although the PmOmpA model is stable over 20 ns simulations, retaining its secondary structure and fold integrity throughout, substantial flexibility is observed in a short linker region between the N- and the C-terminal domains. At low ionic strength, the C-terminal domain moves to interact electrostatically with the lipid bilayer headgroups. This study demonstrates that computational approaches may be applied to more complex, multi-domain outer membrane proteins, rather than just to transmembrane beta-barrels, opening the possibility of in silico proteomics approaches to such proteins.     

The amino terminus of Pseudomonas aeruginosa outer membrane protein OprF forms channels in lipid bilayer membranes: correlation with a three-dimensional model

Pseudomonas aeruginosa OprF forms 0.36-nS channels and, rarely, 2- to 5-nS channels in lipid bilayer membranes. We show that a protein comprising only the N-terminal 162-amino-acid domain of OprF formed the smaller, but not the larger, channels in lipid bilayers. Circular dichroism spectroscopy indicated that this protein folds into a beta-sheet-rich structure, and three-dimensional comparative modeling revealed that it shares significant structural similarity with the amino terminus of the orthologous protein Escherichia coli OmpA, which has been shown to form a beta-barrel. OprF and OmpA share only 15% identity in this domain, yet these results support the utility of modeling such widely divergent beta-barrel domains in three dimensions in order to reveal similarities not readily apparent through primary sequence comparisons. The model is used to further hypothesize why porin activity differs for the N-terminal domains of OprF and OmpA.     

Refolded outer membrane protein A of Escherichia coli forms ion channels with two conductance states in planar lipid bilayers

Outer membrane protein A (OmpA), a major structural protein of the outer membrane of Escherichia coli, consists of an N-terminal 8-stranded beta-barrel transmembrane domain and a C-terminal periplasmic domain. OmpA has served as an excellent model for studying the mechanism of insertion, folding, and assembly of constitutive integral membrane proteins in vivo and in vitro. The function of OmpA is currently not well understood. Particularly, the question whether or not OmpA forms an ion channel and/or nonspecific pore for uncharged larger solutes, as some other porins do, has been controversial. We have incorporated detergent-purified OmpA into planar lipid bilayers and studied its permeability to ions by single channel conductance measurements. In 1 M KCl, OmpA formed small (50-80 pS) and large (260-320 pS) channels. These two conductance states were interconvertible, presumably corresponding to two different conformations of OmpA in the membrane. The smaller channels are associated with the N-terminal transmembrane domain, whereas both domains are required to form the larger channels. The two channel activities provide a new functional assay for the refolding in vitro of the two respective domains of OmpA. Wild-type and five single tryptophan mutants of urea-denatured OmpA are shown to refold into functional channels in lipid bilayers.     

Escherichia coli outer membrane protein A adheres to human brain microvascular endothelial cells

Escherichia coli K1 is the most common gram-negative bacterium causing neonatal meningitis. The outer membrane protein A (OmpA) assembles a beta-barrel structure having four surface-exposed loops in E. coli outer membrane. OmpA of meningitis-causing E. coli K1 is shown to contribute to invasion of the human brain microvascular endothelial cells (HBMEC), the main cellular component of the blood-brain barrier (BBB). However, the direct evidence of OmpA protein interacting with HBMEC is not clear. In this study, we showed that OmpA protein, solubilized from the outer membrane of E. coli, adhered to HBMEC surface. To verify OmpA interaction with the HBMEC, we purified N-terminal membrane-anchoring beta-barrel domain of OmpA and all surface-exposed loops deleted OmpA proteins, and showed that the surface-exposed loops of OmpA were responsible for adherence to HBMEC. These findings indicate that the OmpA is the adhesion molecule with HBMEC and the surface-exposed loops of OmpA are the determinant of this interaction.     

Molecular dynamics simulations of a bacterial autotransporter: NalP from Neisseria meningitidis

NalP is an autotransporter secretory protein found in the outer membrane of Neisseria meningitidis. The crystal structure of the NalP translocator domain revealed a transmembrane beta-barrel containing a central alpha-helix. The role of this alpha-helix, and of the conformational dynamics of the beta-barrel pore have been studied via atomistic molecular dynamics simulations. Three simulations, each of 10 ns duration, of NalP embedded within a solvated DMPC bilayer were performed. The helix was removed from the barrel interior in one simulation. The conformational stability of the protein is similar to that of other outer membrane proteins, e.g., OmpA, in comparable simulations. The transmembrane beta-barrel is stable even in the absence of the alpha-helix. Removal of the helix results in an influx of water into the pore region, suggesting the helix acts as a 'plug'. Water molecules entering the resultant pore form hydrogen bonds with the barrel lining that compensate for the loss of helix-barrel hydrogen bonds. The dimensions of the pore fluctuate over the course of the simulation revealing it to be flexible, but only wide enough to allow transport of the passenger domain in an unfolded or extended conformation. The simulations help us to understand the role of the central helix in plugging the pore and in maintaining the width of the barrel, and show that the NalP monomer is sufficient for the transport of the passenger domain in an unfolded or extended conformation.     

The function of OmpA in Escherichia coli

Outer membrane protein A (OmpA) is a major protein in the Escherichia coli outer membrane. In this study, the function of OmpA in E. coli stress survival was examined. An E. coli K1 ompA-deletion mutant was significantly more sensitive than that of its parent strain to sodium dodecyl sulfate (SDS), cholate, acidic environment, high osmolarity, and pooled human serum. A number of amino acid changes at the extracellular loops of OmpA did not affect the viability of E. coli, while short peptide insertions in the periplasmic turns of the OmpA beta-barrel decreased E. coli resistance to environmental stresses. Moreover, ompA mutants were found to survive much better within brain microvascular endothelial cells than the wild-type strain, supporting that OmpA is a major target in mammalian host cell defense. These results indicated that OmpA plays a vital structural role in E. coli, and suggested that a perfect beta-barrel structure of OmpA is important for outer membrane stability. Based on these results and the published OmpA structural analyses, I propose that OmpA is composed of three functional domains including a hydrophilic extracellular mass, a beta-barrel transmembrane structure, and a peptidoglycan binding domain.     

Folding kinetics of the outer membrane proteins OmpA and FomA into phospholipid bilayers

The folding mechanism of outer membrane proteins (OMPs) of Gram-negative bacteria into lipid bilayers has been studied using OmpA of E. coli and FomA of F. nucleatum as examples. Both, OmpA and FomA are soluble in unfolded form in urea and insert and fold into phospholipid bilayers upon strong dilution of the denaturant urea. OmpA is a structural protein and forms a small ion channel, composed of an 8-stranded transmembrane beta-barrel domain. FomA is a voltage-dependent porin, predicted to form a 14 stranded beta-barrel. Both OMPs fold into a range of model membranes of very different phospholipid compositions. Three membrane-bound folding intermediates of OmpA were discovered in folding studies with dioleoylphosphatidylcholine bilayers that demonstrated a highly synchronized mechanism of secondary and tertiary structure formation of beta-barrel membrane proteins. A study on FomA folding into lipid bilayers indicated the presence of parallel folding pathways for OMPs with larger transmembrane beta-barrels.     

Ion channel formation by N-terminal domain: a common feature of OprFs of Pseudomonas and OmpA of Escherichia coli

The proteolytic fragments of OprFs of Pseudomonas aeruginosa and Pseudomonas fluorescens were identified, respectively, as the first 175 and 177 amino acids from the N-terminal domain. They induced ion channels after reincorporation into planar lipid bilayers (85 and 75 pS, respectively, in 1 M NaCl). A similar conductance value (72 pS) was found for the eight beta-strand OmpA N-terminal domain (OmpA171) of Escherichia coli. We conclude that the N-terminal domain of OprFs is sufficient to induce ion channels and the comparison with OmpA171, provides strong evidence of the existence of an eight-stranded beta-barrel in the N-terminal domain of OprFs.     

OmpA: Gating and dynamics via molecular dynamics simulations

Outer membrane proteins (OMPs) of Gram-negative bacteria have a variety of functions including passive transport, active transport, catalysis, pathogenesis and signal transduction. Whilst the structures of ∼ 25 OMPs are currently known, there is relatively little known about their dynamics in different environments. The outer membrane protein, OmpA from Escherichia coli has been studied extensively in different environments both experimentally and computationally, and thus provides an ideal test case for the study of the dynamics and environmental interactions of outer membrane proteins. We review molecular dynamics simulations of OmpA and its homologues in a variety of different environments and discuss possible mechanisms of pore gating. The transmembrane domain of E. coli OmpA shows subtle differences in dynamics and interactions between a detergent micelle and a lipid bilayer environment. Simulations of the crystallographic unit cell reveal a micelle-like network of detergent molecules interacting with the protein monomers. Simulation and modelling studies emphasise the role of an electrostatic-switch mechanism in the pore-gating mechanism. Simulation studies have been extended to comparative models of OmpA homologues from Pseudomonas aeruginosa (OprF) and Pasteurella multocida (PmOmpA), the latter model including the periplasmic C-terminal domain.     

The major outer membrane protein of Fusobacterium nucleatum (FomA) folds and inserts into lipid bilayers via parallel folding pathways

Membrane protein insertion and folding was studied for the major outer membrane protein of Fusobacterium nucleatum (FomA), which is a voltage-dependent general diffusion porin. The transmembrane domain of FomA forms a beta-barrel that is predicted to consist of 14 beta-strands. Here, unfolded FomA is shown to insert and fold spontaneously and quantitatively into phospholipid bilayers upon dilution of the denaturant urea, which was shown previously only for outer membrane protein A (OmpA) of Escherichia coli. Folding of FomA is demonstrated by circular dichroism and fluorescence spectroscopy, by SDS-polyacrylamide gel electrophoresis, and by single-channel recordings. Refolded FomA had a single-channel conductance of 1.1 nS at 1 M KCl, in agreement with the conductance of FomA isolated from membranes in native form. In contrast to OmpA, which forms a smaller eight-stranded beta-barrel domain, folding kinetics of the larger FomA were slower and provided evidence for parallel folding pathways of FomA into lipid bilayers. Two pathways were observed independent of membrane thickness with two different lipid bilayers, which were either composed of dicapryl phosphatidylcholine or dioleoyl phosphatidylcholine. This is the first observation of parallel membrane insertion and folding pathways of a beta-barrel membrane protein from an unfolded state in urea into lipid bilayers. The kinetics of both folding pathways depended on the chain length of the lipid and on temperature with estimated activation energies of 19 kJ/mol (dicapryl phosphatidylcholine) and 70 kJ/mol (dioleoyl phosphatidylcholine) for the faster pathways.     

Secondary and tertiary structure formation of the beta-barrel membrane protein OmpA is synchronized and depends on membrane thickness

The mechanism of membrane insertion and folding of a beta-barrel membrane protein has been studied using the outer membrane protein A (OmpA) as an example. OmpA forms an eight-stranded beta-barrel that functions as a structural protein and perhaps as an ion channel in the outer membrane of Escherichia coli. OmpA folds spontaneously from a urea-denatured state into lipid bilayers of small unilamellar vesicles. We have used fluorescence spectroscopy, circular dichroism spectroscopy, and gel electrophoresis to investigate basic mechanistic principles of structure formation in OmpA. Folding kinetics followed a second-order rate law and is strongly depended on the hydrophobic thickness of the lipid bilayer. When OmpA was refolded into model membranes of dilaurylphosphatidylcholine, fluorescence kinetics were characterized by a rate constant that was about fivefold higher than the rate constants of formation of secondary and tertiary structure, which were determined by circular dichroism spectroscopy and gel electrophoresis, respectively. The formation of beta-sheet secondary structure and closure of the beta-barrel of OmpA were correlated with the same rate constant and coupled to the insertion of the protein into the lipid bilayer. OmpA, and presumably other beta-barrel membrane proteins therefore do not follow a mechanism according to the two-stage model that has been proposed for the folding of alpha-helical bundle membrane proteins. These different folding mechanisms are likely a consequence of the very different intramolecular hydrogen bonding and hydrophobicity patterns in these two classes of membrane proteins.     

OmpA is the principal nonspecific slow porin of Acinetobacter baumannii

Acinetobacter species show high levels of intrinsic resistance to many antibiotics. The major protein species in the outer membrane of Acinetobacter baumannii does not belong to the high-permeability trimeric porin family, which includes Escherichia coli OmpF/OmpC, and instead is a close homolog of E. coli OmpA and Pseudomonas aeruginosa OprF. We characterized the pore-forming function of this OmpA homolog, OmpA(Ab), by a reconstitution assay. OmpA(Ab) produced very low pore-forming activity, about 70-fold lower than that of OmpF and an activity similar to that of E. coli OmpA and P. aeruginosa OprF. The pore size of the OmpA(Ab) channel was similar to that of OprF, i.e., about 2 nm in diameter. The low permeability of OmpA(Ab) is not due to the inactivation of this protein during purification, because the permeability of the whole A. baumannii outer membrane was also very low. Furthermore, the outer membrane permeability to cephalothin and cephaloridine, measured in intact cells, was about 100-fold lower than that of E. coli K-12. The permeability of cephalothin and cephaloridine in A. baumannii was decreased 2- to 3-fold when the ompA(Ab) gene was deleted. These results show that OmpA(Ab) is the major nonspecific channel in A. baumannii. The low permeability of this porin, together with the presence of constitutive β-lactamases and multidrug efflux pumps, such as AdeABC and AdeIJK, appears to be essential for the high levels of intrinsic resistance to a number of antibiotics.     

Fate of predator and prey proteins during growth of Bdellovibrio bacteriovorus on Escherichia coli and Pseudomonas syringae prey

A two-dimensional electrophoretic analysis of protein distribution followed by identification of selected proteins by mass spectrometry was performed on fresh bdellovibrio cultures containing attack phase cells of the predatory bacterium Bdellovibrio bacteriovorus strain 109J-1 and the remains of an Escherichia coli or a Pseudomonas syringae pv. tomato prey. Cleavage of the peptidoglycan-associated outer membrane proteins (OMPs) OmpA in E. coli and OprF in P. syringae occurred in both prey. The tryptic peptides obtained from the cleavage products of OmpA and OprF were all located within the 19-kDa pronase-resistant N-terminal parts of the corresponding proteins. The predator cell fraction was separated from the prey ghosts in fresh bdellovibrio cultures by centrifugation on a Percoll-sucrose cushion. Proteins from each fraction were separated by two-dimensional electrophoresis and identified by mass spectrometric analysis. As no prey OMP could be detected in the predator cell fraction, it was concluded that prey OMPs are not transferred to the predator, as had been suggested previously. However, a protein from the predator was found bound to ghost cell envelopes. This protein may correspond to a protein earlier suggested to be associated with the prey outer or cytoplasmic membranes. Along with recently described polypeptides from B. bacteriovorus strains 100 and 114, it forms a new family of putative outer membrane proteins.     

In vivo membrane assembly of split variants of the E.coli outer membrane protein OmpA.

The two-domain, 325 residue outer membrane protein OmpA of Escherichia coli is a well-established model for the study of membrane assembly. The N-terminal domain, consisting of approximately 170 amino acid residues, is embedded in the membrane, presumably in the form of a beta-barrel consisting of eight antiparallel transmembrane beta-strands. A set of 16 gene variants carrying deletions in the membrane-embedded domain of OmpA was constructed. When pairs of these mutant genes were co-expressed in E.coli, it was found that a functional OmpA protein could be assembled efficiently from two complementary protein fragments. Assembly was found when the polypeptide chain was split at the second or third periplasmic turn. All four protein termini were located in the periplasmic space. Interestingly, duplication of transmembrane strands five and six led to a variant with an unusual topology: the N-terminus of one fragment and the C-terminus of the other fragment were exposed at the cell surface. This is the first demonstration of correct membrane assembly of split beta-structured membrane proteins. These findings are important for a better understanding of their folding/assembly pathway and may have implications for the development of artificial outer membrane proteins and for the cell surface display of heterologous peptides or proteins.     

Correct folding of the beta-barrel of the human membrane protein VDAC requires a lipid bilayer

Spontaneous membrane insertion and folding of beta-barrel membrane proteins from an unfolded state into lipid bilayers has been shown previously only for few outer membrane proteins of Gram-negative bacteria. Here we investigated membrane insertion and folding of a human membrane protein, the isoform 1 of the voltage-dependent anion-selective channel (hVDAC1) of mitochondrial outer membranes. Two classes of transmembrane proteins with either alpha-helical or beta-barrel membrane domains are known from the solved high-resolution structures. VDAC forms a transmembrane beta-barrel with an additional N-terminal alpha-helix. We demonstrate that similar to bacterial OmpA, urea-unfolded hVDAC1 spontaneously inserts and folds into lipid bilayers upon denaturant dilution in the absence of folding assistants or energy sources like ATP. Recordings of the voltage-dependence of the single channel conductance confirmed folding of hVDAC1 to its active form. hVDAC1 developed first beta-sheet secondary structure in aqueous solution, while the alpha-helical structure was formed in the presence of lipid or detergent. In stark contrast to bacterial beta-barrel membrane proteins, hVDAC1 formed different structures in detergent micelles and phospholipid bilayers, with higher content of beta-sheet and lower content of alpha-helix when inserted and folded into lipid bilayers. Experiments with mixtures of lipid and detergent indicated that the content of beta-sheet secondary structure in hVDAC1 decreased at increased detergent content. Unlike bacterial beta-barrel membrane proteins, hVDAC1 was not stable even in mild detergents such as LDAO or dodecylmaltoside. Spontaneous folding of outer membrane proteins into lipid bilayers indicates that in cells, the main purpose of membrane-inserted or associated assembly factors may be to select and target beta-barrel membrane proteins towards the outer membrane instead of actively assembling them under consumption of energy as described for the translocons of cytoplasmic membranes.     

Pseudomonas aeruginosa porin OprF: properties of the channel

Using ion channel reconstitution in planar lipid bilayers, we examined the channel-forming activity of subfractions of Pseudomonas aeruginosa OprF, which was shown to exist in two different conformations: a minority single domain conformer and a majority two-domain conformer (Sugawara, E., Nestorovich, E. M., Bezrukov, S. M., and Nikaido, H. (2006) J. Biol. Chem. 281, 16220-16229). With the fraction depleted for the single domain conformer, we were unable to detect formation of any channels with well defined conductance levels. With the unfractionated OprF, we saw only rare channel formation. However, with the single domain-enriched fraction of OprF, we observed regular insertion of channels with highly reproducible conductances. Single OprF channels demonstrate rich kinetic behavior exhibiting spontaneous transitions between several subconformations that differ in ionic conductance and radius measured in polymer exclusion experiments. Although we showed that the effective radius of the most conductive conformation exceeds that of the general outer membrane porin of Escherichia coli, OmpF, we also found that a single OprF channel mainly exists in weakly conductive subconformations and switches to the fully open state for a short time only. Therefore, the low permeability of OprF reported earlier may be due to two factors: mainly to the paucity of the single domain conformer in the OprF population and secondly to the predominance of weakly conductive subconformations within the single domain conformer.     

Secondary structure of the outer membrane proteins OmpA of Escherichia coli and OprF of Pseudomonas aeruginosa.

When purified without the use of ionic detergents, both OmpA and OprF proteins contained nearly 20% alpha-helical structures, which disappeared completely upon the addition of sodium dodecyl sulfate. This result suggests that the proteins fold in a similar manner, with an N-terminal, membrane-spanning beta-barrel domain and a C-terminal, globular, periplasmic domain.     

A multidomain outer membrane protein from Pasteurella multocida: Modelling and simulation studies of PmOmpA

PmOmpA is a two-domain outer membrane protein from Pasteurella multocida. The N-terminal domain of PmOmpA is a homologue of the transmembrane β-barrel domain of OmpA from Escherichia coli, whilst the C-terminal domain of PmOmpA is a homologue of the extra-membrane Neisseria meningitidis RmpM C-terminal domain. This enables a model of a complete two domain PmOmpA to be constructed and its conformational dynamics explored via MD simulations of the protein embedded within two different phospholipid bilayers (DMPC and DMPE). The conformational stability of the transmembrane β-barrel is similar to that of a homology model of OprF from Pseudomonas aeruginosa in bilayer simulations. There is a degree of water penetration into the interior of the β-barrel, suggestive of a possible transmembrane pore. Although the PmOmpA model is stable over 20 ns simulations, retaining its secondary structure and fold integrity throughout, substantial flexibility is observed in a short linker region between the N- and the C-terminal domains. At low ionic strength, the C-terminal domain moves to interact electrostatically with the lipid bilayer headgroups. This study demonstrates that computational approaches may be applied to more complex, multi-domain outer membrane proteins, rather than just to transmembrane β-barrels, opening the possibility of in silico proteomics approaches to such proteins.