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

Multidrug efflux in Gram-negative bacteria: structural modifications in active compounds leading to efflux pump avoidance

Gram-negative bacteria cause the majority of critically drug-resistant infections, necessitating the rapid development of new drugs with Gram-negative activity. However, drug design is hampered by the low permeability of the Gram-negative cell envelope and the function of drug efflux pumps, which extrude foreign molecules from the cell. A better understanding of the molecular determinants of compound recognition by efflux pumps is, therefore, essential. Here, we quantitatively analysed the activity of 73,737 compounds, recorded in the publicly accessible database CO-ADD, across three strains of E. coli - the wild-type, the efflux-deficient tolC variant, and the hyper-permeable lpxC variant, to elucidate the molecular principles of evading efflux pumps. We computationally investigated molecular features within this dataset that promote, or reduce, the propensity of being recognised by the TolC-dependent efflux systems in E. coli. Our results show that, alongside a range of physicochemical features, the presence or absence of specific chemical groups in the compounds substantially increases the probability of avoiding efflux. A comparison of our findings with inward permeability data further underscores the primary role of efflux in determining drug bioactivity in Gram-negative bacteria.

Design and directed evolution of noncanonical β-stereoselective metalloglycosidases

Metallohydrolases are ubiquitous in nearly all subclasses of hydrolases, utilizing metal elements to activate a water molecule and facilitate its subsequent dissociation of diverse chemical bonds. However, such a catalytic role of metal ions is rarely found with glycosidases that hydrolyze the glycosidic bonds in sugars. Herein, we design metalloglycosidases by constructing a hydrolytically active Zn-binding site within a barrel-shaped outer membrane protein OmpF. Structure- and mechanism-based redesign and directed evolution have led to the emergence of Zn-dependent glycosidases with catalytic proficiency of 2.8 × 10⁹ and high β-stereoselectivity. Biochemical characterizations suggest that the Zn-binding site constitutes a key catalytic motif along with at least one adjacent acidic residue. This work demonstrates that unprecedented metalloenzymes can be tailor-made, expanding the scope of inorganic reactivities in proteinaceous environments, resetting the structural and functional diversity of metalloenzymes, and providing the potential molecular basis of unidentified metallohydrolases and novel whole-cell biocatalysts.

Data-Driven Derivation of Molecular Substructures That Enhance Drug Activity in Gram-Negative Bacteria

The complex cell envelope of Gram-negative bacteria creates a formidable barrier to antibiotic influx. Reduced drug uptake impedes drug development and contributes to a wide range of drug-resistant bacterial infections, including those caused by extremely resistant species prioritized by the World Health Organization. To develop new and efficient treatments, a better understanding of the molecular features governing Gram-negative permeability is essential. Here, we present a data-driven approach, using matched molecular pair analysis and machine learning on minimal inhibitory concentration data from Gram-positive and Gram-negative bacteria to uncover chemical features that influence Gram-negative bioactivity. We find recurring chemical moieties, of a wider range than previously known, that consistently improve activity and suggest that this insight can be used to optimize compounds for increased Gram-negative uptake. Our findings may help to expand the chemical space of broad-spectrum antibiotics and aid the search for new antibiotic compound classes.

The impact of length variations in the L2 loop on the structure and thermal stability of non-specific porins: The case of OmpCs from the Yersinia pseudotuberculosis complex

Porins are integral proteins of the outer membranes of gram-negative bacteria. In membranes, they exist as homotrimers and the L2 loops contribute to their stability. Comparison of OmpC porins of the Yersinia pseudotuberculosis complex with other enterobacterial porins demonstrated L2 loop length diversity, which is caused by varying numbers of dipeptide/tripeptide repeats. The OmpC porins are highly homologous to each other, and they can be subdivided into five isoforms based on their L2 loop structure. Optical spectroscopy and SDS-PAGE experiments revealed that particularities of the L2 loops affected the structure and thermal stability of the porins. Thermal denaturation studies showed that porins with shorter loops, compared to porins with longer loops, had more stable tertiary and less stable secondary and quaternary structures. According to our comparative modeling results, the L2 loops differ in their structure by adopting different spatial positions and forming different polar bonds with a neighbor monomer. The replacement of asparagine with arginine at the C-terminus of the L2 loop shifts the loop upwards and causes the loss of contacts with the arginine clusters within the pores. The increase in the length of these loops ensures that they shift down toward the pore and restore their contacts with arginines on the channel wall, as is the case in classical nonspecific porins. Despite the fact that the surface charge density varies considerably among the OmpC porins, the L2 loops form a typical negatively charged region in the center of the trimer.

Excavating the surface-associated and secretory proteome of Mycobacterium leprae for identifying vaccines and diagnostic markers relevant immunodominant epitopes

For centuries, Mycobacterium leprae, etiological agent of leprosy, has been afflicting mankind regardless of extensive use of live-attenuated vaccines and antibiotics. Surface-associated and secretory proteins (SASPs) are attractive targets against bacteria. We have integrated biological knowledge with computational approaches and present a proteome-wide identification of SASPs. We also performed computational assignment of immunodominant epitopes as coordinates of prospective antigenic candidates in most important class of SASPs, the outer membrane proteins (OMPs). Exploiting the known protein sequence and structural characteristics shared by the SASPs from bacteria, 17 lipoproteins, 11 secretory and 19 novel OMPs (including 4 essential proteins) were identified in M. leprae As OMPs represent the most exposed antigens on the cell surface, their immunoinformatics analysis showed that the identified 19 OMPs harbor T-cell MHC class I epitopes and class II epitopes against HLA-DR alleles (54), while 15 OMPs present potential T-cell class II epitopes against HLA-DQ alleles (6) and 7 OMPs possess T-cell class II epitopes against HLA-DP alleles (5) of humans. Additionally, 11 M. leprae OMPs were found to have B-cell epitopes and these may be considered as prime candidates for the development of new immunotherapeutics against M. leprae.

Tuning the affinity of anion binding sites in porin channels with negatively charged residues: molecular details for OprP

The cell envelope of the Gram negative opportunistic pathogen Pseudomonas aeruginosa is poorly permeable to many classes of hydrophilic molecules including antibiotics due to the presence of the narrow and selective porins. Here we focused on one of the narrow-channel porins, that is, OprP, which is responsible for the high-affinity uptake of phosphate ions. Its two central binding sites for phosphate contain a number of positively charged amino acids together with a single negatively charged residue (D94). The presence of this negatively charged residue in a binding site for negatively charged phosphate ions is highly surprising due to the potentially reduced binding affinity. The goal of this study was to better understand the role of D94 in phosphate binding, selectivity, and transport using a combination of mutagenesis, electrophysiology, and free-energy calculations. The presence of a negatively charged residue in the binding site is critical for this specific porin OprP as emphasized by the evolutionary conservation of such negatively charged residue in the binding site of several anion-selective porins. Mutations of D94 in OprP to any positively charged or neutral residue increased the binding affinity of phosphate for OprP. Detailed analysis indicated that this anionic residue in the phosphate binding site of OprP, despite its negative charge, maintained energetically favorable phosphate binding sites in the central region of the channel and at the same time decreased residence time thus preventing excessively strong binding of phosphate that would oppose phosphate flux through the channel. Intriguingly mutations of D94 to positively charged residues, lysine and arginine, resulted in very different binding affinities and free energy profiles, indicating the importance of side chain conformations of these positively charged residues in phosphate binding to OprP.

Proteome-scale identification of outer membrane proteins in Mycobacterium avium subspecies paratuberculosis using a structure based combined hierarchical approach

The overall strategy used for the proteome-wide comprehensive computational investigation to identify outer membrane proteins fromMycobacterium aviumsubsp.paratuberculosis.

Ion permeation in the NanC porin from Escherichia coli: free energy calculations along pathways identified by coarse-grain simulations

Using the X-ray structure of a recently discovered bacterial protein, the N-acetylneuraminic acid-inducible channel (NanC), we investigate computationally K(+) and Cl(-) ions' permeation. We identify ion permeation pathways that are likely to be populated using coarse-grain Monte Carlo simulations. Next, we use these pathways as reaction coordinates for umbrella sampling-based free energy simulations. We find distinct tubelike pathways connecting specific binding sites for K(+) and, more pronounced, for Cl(-) ions. Both ions permeate the porin preserving almost all of their first hydration shell. The calculated free energy barriers are G(#) ≈ 4 kJ/mol and G(#) ≈ 8 kJ/mol for Cl(-) and K(+), respectively. Within the approximations associated with these values, discussed in detail in this work, we suggest that the porin is slightly selective for Cl(-) versus K(+). Our suggestion is consistent with the experimentally observed weak Cl(-) over K(+) selectivity. A rationale for the latter is suggested by a comparison with previous calculations on strongly anion selective porins.

Structure and function of the PorB porin from disseminating Neisseria gonorrhoeae

The outer membrane of Gram-negative bacteria contains a large number of channel-forming proteins, porins, for the uptake of small nutrient molecules. Neisseria gonorrhoeae PorBIA (PorB of serotype A) are associated with disseminating diseases and mediate a rapid bacterial invasion into host cells in a phosphate-sensitive manner. To gain insights into this structure-function relationship we analysed PorBIA by X-ray crystallography in the presence of phosphate and ATP. The structure of PorBIA in the complex solved at a resolution of 3.3 Å (1 Å=0.1 nm) displays a surplus of positive charges inside the channel. ATP ligand-binding in the channel is co-ordinated by the positively charged residues of the channel interior. These residues ligate the aromatic, sugar and pyrophosphate moieties of the ligand. Two phosphate ions were observed in the structure, one of which clamped by two arginine residues (Arg92 and Arg124) localized at the extraplasmic channel exit. A short β-bulge in β2-strand together with the long L3 loop narrow the barrel diameter significantly and further support substrate specificity through hydrogen bond interactions. Interestingly the structure also comprised a small peptide as a remnant of a periplasmic protein which physically links porin molecules to the peptidoglycan network. To test the importance of Arg92 on bacterial invasion the residue was mutated. In vivo assays of bacteria carrying a R92S mutation confirmed the importance of this residue for host-cell invasion. Furthermore systematic sequence and structure comparisons of PorBIA from Neisseriaceae indicated Arg92 to be unique in disseminating N. gonorrhoeae thereby possibly distinguishing invasion-promoting porins from other neisserial porins.

Continuum electrostatic calculations of the pKa of ionizable residues in an ion channel: dynamic vs. static input structure

We have computed the pK(a)'s of the ionizable residues of a protein ion channel, the Staphylococcus aureus toxin alpha-hemolysin, by using two types of input structures, namely the crystal structure of the heptameric alpha-hemolysin and a set of over four hundred snapshots from a 4.38 ns Molecular Dynamics simulation of the protein inserted in a phospholipid planar bilayer. The comparison of the dynamic picture provided by the Molecular Simulation with the static one based on the X-ray crystal structure of the protein embedded in a lipid membrane allows analyzing the influence of the fluctuations in the protein structure on its ionization properties. We find that the use of the dynamic structure provides interesting information about the sensitivity of the computed pK(a) of a given residue to small changes in the local structure. The calculated pK(a) are consistent with previous indirect estimations obtained from single-channel conductance and selectivity measurements.

The simulation approach to bacterial outer membrane proteins (Review)

The outer membrane of Gram-negative bacteria serves as a protective barrier against the external environment but is rendered selectively permeable to nutrients and waste by proteins called porins. Other outer membrane proteins (OMPs) provide the membrane with a variety of other functions including active transport, catalysis, pathogenesis and signal transduction. A relatively small number of crystal or NMR structures of these proteins are known, and it is therefore essential that the maximum possible information be extracted. In this respect, computational techniques enable extrapolation from time- and space-averaged static structures to dynamic, physiological events. Electrostatics approaches have been used to investigate the structures of porins. The stochastic simulation of ion trajectories through these channels has been possible with Brownian dynamics, which treats the membrane and solvent approximately, enabling the prediction of conduction properties. Molecular dynamics has also been applied, enabling fully atomistic descriptions of 'virtual outer membranes'. This has provided atomic resolution descriptions of solute permeation through porins. It has also yielded insights into the dynamics of gating in active transporters and ion channels, as well as providing clues to catalytic mechanisms in outer membrane enzymes. Additionally, simulations are beginning to reveal the common features of interactions between membrane proteins and lipids, with biological implications for OMP folding, stability and mechanism. Future prospects include the simulation of longer, larger and more complex outer membrane systems, with more accurate descriptions of inter-atomic forces.

Identification of outer membrane proteins of Mycobacterium tuberculosis

The cell wall of mycobacteria includes an unusual outer membrane of extremely low permeability. While Escherichia coli uses more than 60 proteins to functionalize its outer membrane, only two mycobacterial outer membrane proteins (OMPs) are known. The porin MspA of Mycobacterium smegmatis provided the proof of principle that integral mycobacterial OMPs share the beta-barrel structure, the absence of hydrophobic alpha-helices and the presence of a signal peptide with OMPs of gram-negative bacteria. These properties were exploited in a multi-step bioinformatic approach to predict OMPs of M. tuberculosis. A secondary structure analysis was performed for 587 proteins of M. tuberculosis predicted to be exported. Scores were calculated for the beta-strand content and the amphiphilicity of the beta-strands. Reference OMPs of gram-negative bacteria defined threshold values for these parameters that were met by 144 proteins of unknown function of M. tuberculosis. Two of them were verified as OMPs by a novel two-step experimental approach. Rv1698 and Rv1973 were detected only in the total membrane fraction of M. bovis BCG in Western blot experiments, while proteinase K digestion of whole cells showed the surface accessibility of these proteins. These findings established that Rv1698 and Rv1973 are indeed localized in the outer membrane and tripled the number of known OMPs of M. tuberculosis. Significantly, these results provide evidence for the usefulness of the bioinformatic approach to predict mycobacterial OMPs and indicate that M. tuberculosis likely has many OMPs with beta-barrel structure. Our findings pave the way to identify the set of proteins which functionalize the outer membrane of M. tuberculosis.

An arginine ladder in OprP mediates phosphate-specific transfer across the outer membrane

The outer membrane protein OprP mediates the transport of essential phosphate anions into the pathogenic bacterium Pseudomonas aeruginosa. Here we report the crystallographic structure of trimeric OprP at 1.9-A resolution, revealing an unprecedented 9-residue arginine 'ladder' that spans from the extracellular surface down through a constriction zone where phosphate is coordinated. Lysine residues coat the inner periplasmic surface, creating an 'electropositive sink' that pulls the phosphates through the eyelet and into the cell.

High resolution crystal structures and molecular dynamics studies reveal substrate binding in the porin Omp32

The porin Omp32 is the major outer membrane protein of the bacterium Delftia acidovorans. The crystal structures of the strongly anion-selective porin alone and in complex with the substrate malate were solved at 1.5 and 1.45 A resolution, respectively, and revealed a malate-binding motif adjacent to the channel constriction zone. Binding is mediated by interaction with a cluster of two arginine residues and two threonines. This binding site is specific for Omp32 and reflects the physiological adaptation of the organism to organic acids. Structural studies are combined with a 7-ns unbiased molecular dynamics simulation of the trimeric channel in a model membrane. Molecular dynamics trajectories show how malate ions are efficiently captured from the surrounding bulk solution by the electrostatic potential of the channel, translocated to the binding site region, and immobilized in the constriction zone. In accordance with these results, conductance measurements with Omp32 inserted in planar lipid membranes revealed binding of malate. The anion-selective channel Omp32 is the first reported example of a porin with a 16-stranded beta-barrel and proven substrate specificity. This finding suggests a new view on the correlation of porin structure with substrate binding in specific channels.

Topology of the porin MspA in the outer membrane of Mycobacterium smegmatis

MspA is the major porin of Mycobacterium smegmatis mediating the exchange of hydrophilic solutes across the outer membrane (OM). It is the prototype of a new family of octameric porins with a single central channel of 9.6 nm in length and consists of two hydrophobic beta-barrels of 3.7 nm in length and a more hydrophilic, globular rim domain. The length of the hydrophobic domain of MspA does not match the thicknesses of mycobacterial OMs of 5-12 nm as derived from electron micrographs. Further, the membrane topology of MspA is unknown as it is for any other mycobacterial OM protein. We used MspA as a molecular ruler to define the boundaries of the OM of M. smegmatis by surface labeling of single cysteine mutants. Seventeen mutants covered the surface of the rim domain and were biotinylated with a membrane-impermeable reagent. The label efficiencies in vitro were remarkably similar to the predicted accessibilities of the cysteines. By contrast, six of these mutants were protected from biotinylation in M. smegmatis cells. Tryptophan 21 defines a horizontal plane that dissects the surface-exposed versus the membrane-protected residues of MspA. The 8 phenylalanines at position 99 form a ring at the periplasmic end of the hydrophobic beta-barrel domain. These results indicated that (i) the membrane boundaries of MspA are defined by aromatic girdles as in porins of Gram-negative bacteria and (ii) loops and a 3.4-nm long part of the hydrophilic rim domain are embedded into the OM of M. smegmatis. This is the first report suggesting that elements other than hydrophobic alpha-helices or beta-sheets are integrated into a lipid membrane.

Chemical modification of the bacterial porin OmpF: gain of selectivity by volume reduction

OmpF is an essentially nonselective porin isolated from the outer membrane of Escherichia coli. Here we report on the manipulation of the ion selectivity of OmpF by chemical modification with MTS reagents (MTSET, MTSEA, and MTSES) and the (rather bulky) tripeptide glutathione, all cysteine specific. When recorded in a gradient of 0.1//1 M CaCl2 or 0.1//1 M NaCl, pH 7.4 solutions, measured reversal potentials of the most cation-selective modified mutants were (virtually) identical to the Nernst potential of Ca2+ or Na+. Compared to this full cation selectivity, the anion-selective modified mutants performed somewhat less but nevertheless showed high anion selectivity. We conclude that a low permanent charge in combination with a narrow pore can render the same selectivity as a highly charged but wider pore. These results favor the view that both the electrostatic potential arising form the fixed charge in the pore and the space available at the selectivity filter contribute to the charge selection (i.e., cation versus anion selectivity) of a biological ion channel.

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.

Improved 3D continuum calculations of ion flux through membrane channels

A continuum model, based on the Poisson-Nernst-Planck (PNP) theory, is applied to simulate steady-state ion flux through protein channels. The PNP equations are modified to explicitly account (1) for the desolvation of mobile ions in the membrane pore and (2) for effects related to ion sizes. The proposed algorithm for a three-dimensional self-consistent solution of PNP equations, in which final results are refined by a focusing technique, is shown to be suitable for arbitrary channel geometry and arbitrary protein charge distribution. The role of the pore shape and protein charge distribution in formation of basic electrodiffusion properties, such as channel conductivity and selectivity, as well as concentration distributions of mobile ions in the pore region, are illustrated by simulations on model channels. The influence of the ionic strength in the bulk solution and of the externally applied electric field on channel properties are also discussed.

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.