Modeling and simulations of a bacterial outer membrane protein: OprF from Pseudomonas aeruginosa

OprF is a major outer membrane protein from Pseudomonas aeruginosa, a homolog of OmpA from Escherichia coli. The N-terminal domains of both proteins have been demonstrated to form low conductance channels in lipid bilayers. Homology models, consisting of an eight-stranded beta-barrel, of the N-terminal domain OprF have been constructed based on the crystal structure of the corresponding domain from E. coli OmpA. OprF homology models have been evaluated via a set (6 x 10 ns) of simulations of the beta-barrel embedded within a solvated dimyristoyl-phosphatidylcholine (DMPC) bilayer. The conformational stability of the models is similar to that of the crystal structure of OmpA in comparable simulations. There is a degree of water penetration into the pore-like center of the OprF barrel. The presence of an acidic/basic (E8/K121) side-chain interaction within the OprF barrel may form a "gate" able to close/open a central pore. Lipid-protein interactions within the simulations were analyzed and revealed that aromatic side-chains (Trp, Tyr) of OprF interact with lipid headgroups. Overall, the behavior of the OprF model in simulations supports the suggestion that this molecule is comparable to OmpA. The simulations help to explain the mechanism of formation of low conductance pores within the outer membrane.     

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.     

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.     

Pseudomonas aeruginosa porin OprF exists in two different conformations

The major nonspecific porin of Pseudomonas aeruginosa, OprF, produces a large channel yet allows only a slow diffusion of various solutes. Here we provide an explanation of this apparent paradox. We first show, by introduction of tobacco etch virus protease cleavage site in the middle of OprF protein, that most of OprF population folds as a two-domain protein with an N-terminal beta-barrel domain and a C-terminal periplasmic domain rich in alpha-helices. However, sedimentation of unilamellar proteoliposomes through an iso-osmotic gradient showed that only about 5% of the OprF population produced open channels. Gel filtration showed that the open channel conformers tended to occur in oligomeric associations. Because the open channel conformer is likely to fold as a single domain protein with a large beta-barrel, we reasoned that residues near the C terminus may be exposed on cell surface in this conformer. Introduction of a cysteine residue at position 312 produced a functional mutant protein. By using bulky biotinylation reagents on intact cells, we showed that this cysteine residue was not exposed on cell surface in most of the OprF population. However, the minority OprF population that was biotinylated in such experiments was enriched for the conformer with pore-forming activity and had a 10-fold higher pore-forming specific activity than the bulk OprF population. Finally trypsin treatment, which preferentially cleaves the C-terminal domain of the two-domain conformer, did not affect the pore-forming activity of OprF nor did it digest the minority conformer whose residue 312 is exposed on cell surface.     

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.     

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.     

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.     

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.     

The N-terminal beta-barrel structure of lipid body lipoxygenase mediates its binding to liposomes and lipid bodies.

Phospholipase A2 and a particular isoform of lipoxygenase are synthesized and transferred to lipid bodies during the stage of triacylglycerol mobilization in germinating cucumber seedlings. Lipid body lipoxygenase (LBLOX) is post-translationally transported to lipid bodies without proteolytic modification. Fractionation of homogenates from cucumber cotyledons or transgenic tobacco leaves expressing LBLOX showed that a small but significant amount was detectable in the microsomal fraction. A beta-barrel-forming N-terminal domain in the structure of LBLOX, as deduced from sequence data, was shown to be crucial for selective intracellular transport from the cytosol to lipid bodies. Although a specific signal sequence for targeting protein domains to the lipid bodies could not be established, it was evident that the beta-barrel represents a membrane-binding domain that is functionally comparable with the C2 domains of mammalian phospholipases. The intact beta-barrel of LBLOX was demonstrated to be sufficient to target in vitro a fusion protein of LBLOX beta-barrel with glutathione S-transferase (GST) to lipid bodies. In addition, binding experiments on liposomes using lipoxygenase isoforms, LBLOX deletions and the GST-fusion protein confirmed the role of the beta-barrel as the membrane-targeting domain. In this respect, the cucumber LBLOX differs from cytosolic isoforms in cucumber and from the soybean LOX-1. When the beta-barrel of LBLOX was destroyed by insertion of an additional peptide sequence, its ability to target proteins to membranes was abolished.     

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.     

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.     

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.     

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.     

Limited proteolysis and biophysical characterization of the lipovitellin homology region in apolipoprotein B

Apolipoprotein B (apoB) is the essential nonexchangeable protein in chylomicrons and very low-density lipoprotein-derived lipoprotein particles, including low-density lipoprotein (LDL). ApoB has been a key target for cardiovascular research because of its essential role in the assembly, secretion, delivery, and receptor binding of LDL. The three-dimensional structure of apoB has not been determined. However, the N-terminal region of apoB is homologous to the lipid storage protein lipovitellin, which allows the modeling of this region based on the X-ray structure of lipovitellin. The model of the N-terminal 17% of apoB (B17) suggests that, like lipovitellin, B17 consists of an N-terminal beta-barrel domain, a helical domain, and a beta-sheet domain (C-sheet). Here we test the validity of this model by limited proteolysis of B17 and the characterization of individual domains expressed in Escherichia coli and insect cell systems that are consistent with the model and proteolysis data. Circular dichroism studies of the individual domains indicate that they are folded and their secondary structures are in agreement with the model. We find that the helical domain and C-sheet of apoB interact with each other in vitro, suggesting a strong interaction between these two domains, even without a covalent peptide bond linkage. Our data suggest that the three lipovitellin-like domains exist in B17. Furthermore, the domains fold independently with secondary structures and stabilities like those of intact B17.     

A minimal transmembrane beta-barrel platform protein studied by nuclear magnetic resonance

In this study, we were concerned with the structural role of the surface-exposed extracellular loops of the N-terminal transmembrane (TM) domain of OmpA. A variant of the TM domain of outer membrane protein A (OmpA) with all four such loops shortened, which we call the beta-barrel platform (BBP), was successfully refolded. This indicates that the removed parts of the surface-exposed loops indeed do not contain amino acid sequences critical for this membrane protein's refolding in vitro. BBP has the potential to be used as a template beta-barrel membrane protein structure for the development of novel functions, although our results also highlight the potential difficulties that can arise when functionality is being engineered into the loop regions of membrane proteins. We have used solution nuclear magnetic resonance spectroscopy to determine the global fold of BBP+EF, BBP with a metal ion-binding EF-hand inserted in one of the shortened loops. BBP and BBP+EF in dihexanoylphosphatidylcholine micelles are eight-stranded antiparallel beta-barrels, and BBP represents the smallest beta-structured integral membrane protein known to date.     

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).     

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.     

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.     

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.     

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.