Molecular analysis of the membrane insertion domain of proteinase 3, the Wegener's autoantigen, in RBL cells: implication for its pathogenic activity.

PR3, also called myeloblastin, is a neutrophil serine protease that promotes myeloid cell proliferation by cleaving the cyclin-dependent kinase inhibitor p21(cip1/waf1). In addition, it is the target of ANCA in GPA, a necrotizing vasculitis. Anti-PR3 ANCA binding to membrane-expressed PR3 triggers neutrophil activation, potentiating vascular inflammation. This study performed in RBL cells identifies the structural motifs of PR3 membrane anchorage and examines its impact on PR3 proinflammatory and proliferative functions. With the use of MD simulations and mutagenesis, we demonstrate that the mutations of four hydrophobic (F180, F181, L228, F229) or four basic (R193, R194, K195, R227) amino acids abrogated PR3 membrane anchorage. The hydrophobic patch-deficient PR3 mutant (PR34H4A) was still able to cleave the synthetic substrate Boc-Ala-Pro-Val in cell lysates. However, in contrast to WT PR3, PR34H4A was not expressed at the plasma membrane after degranulation and failed to cleave extracellular fibronectin, was not externalized after apoptosis and did not impair macrophage phagocytosis of apoptotic cells, did not promote myeloid cell proliferation and failed to cleave p21/waf1. PR3 membrane insertion appears to be pivotal for its proinflammatory activities, such as extracellular proteolysis and impairment of apoptotic cell clearance, but also for myeloid cell proliferation. Targeting membrane-associated PR3 might constitute a novel, anti-inflammatory therapeutic strategy in inflammatory disease especially in vasculitis, but this approach has to be validated in mature neutrophils.

Molecular analysis of antimicrobial agent translocation through the membrane porin BpsOmp38 from an ultraresistant Burkholderia pseudomallei strain.

Burkholderia pseudomallei (Bps) is a Gram-negative bacterium that causes melioidosis, an infectious disease of animals and humans common in northern and north-eastern parts of Thailand. Successful treatment of melioidosis is difficult due to intrinsic resistance of Bps to various antibacterial agents. It has been suggested that the antimicrobial resistance of this organism may result from poor permeability of the active compounds through porin channels located in the outer membrane (OM) of the bacterium. In previous work, a 38-kDa protein, named "BpsOmp38", was isolated from the OM of Bps. A topology prediction and liposome-swelling assay suggested that BpsOmp38 comprises a ß-barrel structure and acts as a general diffusion porin. The present study employed black lipid membrane (BLM) reconstitution to demonstrate the single-channel conductance of the trimeric BpsOmp38 to be 2.7±0.3 nS in 1M KCl. High-time resolution BLM measurements displayed ion current blockages of seven antimicrobial agents in a concentration-dependent manner with the translocation on-rate (kon) following the order: norfloxacin≫ertapenem>ceftazidime>cefepime>imipenem>meropenem>penicillin G. The dwell time of a selected antimicrobial agent (ertapenem) decayed exponentially with increasing temperature. The energy barrier for the ertapenem binding to the affinity site inside the BpsOmp38 channel was estimated from the Arrhenius plot to be 12 kT and for the ertapenem release to be 13 kT at +100 mV. The BLM data obtained from this study provide the first insight into antimicrobial agent translocation through the BpsOmp38 channel.

Molecular simulations reveal the mechanism and the determinants for ampicillin translocation through OmpF.

We use a multiscale approach, combining molecular dynamics simulations with metadynamics, to simulate the translocation of ampicillin through OmpF from Escherichia coli (E. coli). In-depth analysis has allowed us to reveal the complete picture of the translocation process in terms of both energetics and physicochemical properties. We have demonstrated the existence of a unique affinity site at the constriction region, accessible from both sides and defined by specific pore-antibiotic interactions. By providing optimal binding, the constriction region works like an enzyme toward the permeation of ampicillin. We find reduction in entropy to be compensated by enthalpic contributions from a favorable network of interactions (hydrogen bonds and hydrophobic contacts) which is also mediated by two slow water molecules bridging the antibiotic-pore interactions. Finally, as ampicillin assumes a preferential value for a torsional angle when at the constriction region, we investigated the consequence of the conformational preorganization of ampicillin toward its translocation. As a whole, our analysis opens the way to chemical modifications of antibiotics to allow improving uptake through porins contributing to combat bacterial resistance.

Toward screening for antibiotics with enhanced permeation properties through bacterial porins.

Gram-negative bacteria are protected by an outer membrane barrier, and to reach their periplasmic target, penicillins have to diffuse through outer membrane porins such as OmpF. Here we propose a structure-dynamics-based strategy for improving such antibiotic uptake. Using a variety of experiments (high-resolution single channel recording, Minimum Inhibitory Concentration (MIC), liposome swelling assay) and accelerated molecular simulations, we decipher the subtle balance of interactions governing ampicillin diffusion through the porin OmpF. This suggests mutagenesis of a hot spot residue of OmpF for which additional simulations reveal drastic changes in the molecular and energetic pathway of ampicillin's diffusion. Inverting the problem, we predict and describe how benzylpenicillin diffuses with a lower effective energy barrier by interacting differently with OmpF. The thorough comparison between the theoretical predictions and the three independent experiments, which were set up to measure the kinetics of transport and biological activity, gives insights on how to combine such different investigation techniques with the aim of providing complementary validation. Our study illustrates the importance of microscopic interactions at the constriction region of the biological channel to control the antibiotic flux through it. We conclude by providing a complete inventory of the channel and antibiotic hot spots and discuss the implications in terms of antibacterial screening and design.

Implication of porins in beta-lactam resistance of Providencia stuartii.

An integrative approach combining biophysical and microbiological methods was used to characterize the antibiotic translocation through the outer membrane of Providencia stuartii. Two novel members of the General Bacterial Porin family of Enterobacteriaceae, named OmpPst1 and OmpPst2, were identified in P. stuartii. In the presence of ertapenem (ERT), cefepime (FEP), and cefoxitin (FOX) in growth media, several resistant derivatives of P. stuartii ATCC 29914 showed OmpPst1-deficiency. These porin-deficient strains showed significant decrease of susceptibility to ß-lactam antibiotics. OmpPst1 and OmpPst2 were purified to homogeneity and reconstituted into planar lipid bilayers to study their biophysical characteristics and their interactions with ß-lactam molecules. Determination of ß-lactam translocation through OmpPst1 and OmpPst2 indicated that the strength of interaction decreased in the order of ertapenem ≫ cefepime > cefoxitin. Moreover, the translocation of these antibiotics through OmpPst1 was more efficient than through OmpPst2. Heterologous expression of OmpPst1 in the porin-deficient E. coli strain BL21(DE3)omp8 was associated with a higher antibiotic susceptibility of the E. coli cells to ß-lactams compared with expression of OmpPst2. All our data enlighten the involvement of porins in the resistance of P. stuartii to ß-lactam antibiotics.

Investigating reaction pathways in rare events simulations of antibiotics diffusion through protein channels.

In Gram-negative bacteria, outer-membrane protein channels, such as OmpF of Escherichia coli, constitute the entry point of various classes of antibiotics. While antibacterial research and development is declining, bacterial resistance to antibiotics is rising and there is an emergency call for a new way to develop potent antibacterial agents and to bring them to the market faster and at reduced cost. An emerging strategy is to follow a bottom-up approach based on microscopically founded computational based screening, however such strategy needs better-tuned methods. Here we propose to use molecular dynamics (MD) simulations combined with the metadynamics algorithm, to study antibiotic translocation through OmpF at a molecular scale. This recently designed algorithm overcomes the time scale problem of classical MD by accelerating some reaction coordinates. It is expected that the initial assumption of the reaction coordinates is a key determinant for the efficiency and accuracy of the simulations. Previous studies using different computational schemes for a similar process only used one reaction coordinate, which is the directionality. Here we go further and see how it is possible to include more informative reaction coordinates, accounting explicitly for: (i) the antibiotic flexibility and (ii) interactions with the channel. As model systems, we select two compounds covering the main classes of antibiotics, ampicillin and moxifloxacine. We decipher the molecular mechanism of translocation of each antibiotic and highlight the important parameters that should be taken into account for improving further simulations. This will benefit the screening and design for antibiotics with better permeation properties.

Structures of human proteinase 3 and neutrophil elastase--so similar yet so different.

Proteinase 3 and neutrophil elastase are serine proteinases of the polymorphonuclear neutrophils, which are considered to have both similar localization and ligand specificity because of their high sequence similarity. However, recent studies indicate that they might have different and yet complementary physiologic roles. Specifically, proteinase 3 has intracellular specific protein substrates resulting in its involvement in the regulation of intracellular functions such as proliferation or apoptosis. It behaves as a peripheral membrane protein and its membrane expression is a risk factor in chronic inflammatory diseases. Moreover, in contrast to human neutrophil elastase, proteinase 3 is the preferred target antigen in Wegener's granulomatosis, a particular type of vasculitis. We review the structural basis for the different ligand specificities and membrane binding mechanisms of both enzymes, as well as the putative anti-neutrophil cytoplasm autoantibody epitopes on human neutrophil elastase 3. We also address the differences existing between murine and human enzymes, and their consequences with respect to the development of animal models for the study of human proteinase 3-related pathologies. By integrating the functional and the structural data, we assemble many pieces of a complicated puzzle to provide a new perspective on the structure-function relationship of human proteinase 3 and its interaction with membrane, partner proteins or cleavable substrates. Hence, precise and meticulous structural studies are essential tools for the rational design of specific proteinase 3 substrates or competitive ligands that modulate its activities.

Molecular basis of enrofloxacin translocation through OmpF, an outer membrane channel of Escherichia coli--when binding does not imply translocation.

The molecular pathway of enrofloxacin, a fluoroquinolone antibiotic, through the outer membrane channel OmpF of Escherichia coli is investigated. High-resolution ion current fluctuation analysis reveals a strong affinity for enrofloxacin to OmpF, the highest value ever recorded for an antibiotic-channel interaction. A single point mutation in the constriction zone of OmpF, replacing aspartic acid at the 113 position with asparagine (D113N), lowers the affinity to a level comparable to other antibiotics. All-atom molecular dynamics simulations allow rationalizing the translocation pathways: wild-type OmpF has two symmetric binding sites for enrofloxacin located at each channel entry separated by a large energy barrier in the center, which inhibits antibiotic translocation. In this particular case, our simulations suggest that the ion current blockages are caused by molecules occupying either one of these peripheral binding sites. Removal of the negative charge on position 113 removes the central barrier and shifts the two peripheral binding sites to a unique central site, which facilitates translocation. Fluorescence steady-state measurements agree with the different location of binding sites for wild-type OmpF and the mutant. Our results demonstrate how a single-point mutation of the porin, and the resulting intrachannel shift of the affinity site, may substantially modify translocation.

Bridging timescales and length scales: from macroscopic flux to the molecular mechanism of antibiotic diffusion through porins.

Our aim in this study was to provide an atomic description of ampicillin translocation through OmpF, the major outer membrane channel in Escherichia coli and main entry point for beta-lactam antibiotics. By applying metadynamics simulations, we also obtained the energy barriers along the diffusion pathway. We then studied the effect of mutations that affect the charge and size at the channel constriction zone, and found that in comparison to the wild-type, much lower energy barriers are required for translocation. The expected higher translocation rates were confirmed on the macroscopic scale by liposome-swelling assays. A microscopic view on the millisecond timescale was obtained by analysis of temperature-dependent ion current fluctuations in the presence of ampicillin and provide the enthalpic part of the energy barrier. By studying antibiotic translocation over various timescales and length scales, we were able to discern its molecular mechanism and rate-limiting interactions, and draw biologically relevant conclusions that may help in the design of drugs with enhanced permeation rates.

Structural and dynamical properties of the porins OmpF and OmpC: insights from molecular simulations.

In this paper we investigate the structural and dynamical properties of the two major porins (OmpF and OmpC) in Escherichia coli, using molecular dynamics (MD) simulations. In particular we characterized the atomic fluctuations, correlated motions, temperature dependence, solvent-accessible cross-sectional area and water dynamics in the key regions of the two channels. Our in-depth analysis allows us to highlight the importance of both the key conserved and substituted residues between OmpF and OmpC. The latter is characterized by a narrower and longer constriction region with respect to OmpF. OmpC also showed a higher stability upon increasing temperature. We then present the results of transport properties by using accelerated MD simulations to probe the diffusion of norfloxacin (a fluoroquinolone antibiotic) through the two porins OmpF/OmpC. Our study constitutes a step forward towards understanding the structure-function relationship of the two porins' channels. This will benefit the research of antibacterials with improved permeation properties and nanopores that aim to use these porins as sensing systems.

Challenges in pKa predictions for proteins: the case of Asp213 in human proteinase 3.

Knowledge of the protonation states of the ionizable residues in an enzyme is a prerequisite to an accurate description of its structure and mechanism. In practice, the use of the inappropriate protonation state for an amino acid in a molecular modeling computation (e.g., molecular dynamics simulation) is likely to lead to unrealistic results. Although methods using solvers of the linearized Poisson-Boltzmann equation have proven to yield accurate pK(a) predictions, they bear a number of limitations. They are quite demanding in terms of computational power and are sensitive to representation of the charges and their position (force field and protein conformation). Moreover they depend on the choice of a dielectric constant for the protein interior. In this manuscript, we describe the difficulties met when trying to predict the protonation state of a buried amino acid, located in a protein for which very little biochemical data is available. Such a case is highly representative of the challenges faced in theoretical biology studies. Proteinase 3 (PR3) is an enzyme involved in proteolytic events associated with inflammation. It is a potential target in the development of new anti-inflammatory therapeutic strategies. We report the results of pK(a) predictions of the aspartic acid 213 of PR3 with a FDPB solver. We probed the influence of the choice of the dielectric constant for the protein interior epsilon(p) and the benefits of conformational sampling by molecular dynamics (MD) on the pK(a) prediction of this carboxylate group. Using only the FDPB calculations, we could not conclude on the protonation state of Asp213. MD simulations confronted to knowledge of the ligand-binding and reaction mechanism led us to decide on a protonated form of this aspartic acid. We also demonstrate that the use of the wrong protonation state leads to an unreliable structural model for PR3. pK(a) prediction with a fast empirical method yielded a pK(a) of 8.4 for Asp213, which is in agreement with our choice of protonation state based on MD simulations.

Computational prediction of the binding site of proteinase 3 to the plasma membrane.

Proteinase 3 (PR3) is a neutrophil-derived serine proteinase localized within cytoplasmic granules which can be released upon activation. PR3 is exposed at the neutrophil plasma membrane where it can mediate proinflammatory effects. Moreover, PR3 membrane expression is of special relevance in patients with Wegener's granulomatosis, a systemic vasculitis presenting anticytoplasmic neutrophil autoantibodies (ANCA) against PR3, which can bind to PR3 expressed at the surface of neutrophils and amplify their activation state. Therefore, it is of special relevance to unravel the molecular mechanisms governing its association with the membrane to be able to modulate it. To this end, we performed molecular dynamics (MD) simulations of PR3 with the implicit membrane model IMM1-GC to identify its interfacial binding site (IBS). Both the energies and structures resulting from the MD suggest that PR3 associates strongly with anionic membranes. We observe a unique IBS consisting of five basic (R177, R186A, R186B, K187, R222) and six hydrophobic (F165, F166, F224, L223, F184, W218) amino acids. The basic residues provide the driving force to orient PR3 at the membrane surface, so that the hydrophobic residues can anchor into the hydrocarbon region. Energy decomposition and in silico mutations show that only a few residues account for the membrane association. Similar calculations with HNE suggest a different membrane-binding mechanism. Our results agree with previous experimental observations and this work predicts, for the first time, the structural determinants of the binding of PR3 to membranes.

Differences in the substrate binding sites of murine and human proteinase 3 and neutrophil elastase.

Understanding the differences between murine (m) and human (h) proteinase 3 (PR3) and neutrophil elastase (NE) is crucial for the interpretation of in vivo studies of inflammatory processes. We built structural models of mPR3 and mNE and analyzed their surface properties. We performed molecular dynamics (MD) simulations on several enzyme-peptide complexes to investigate their interaction patterns. The analysis of trajectories confirms that murine and human complexes have different interaction patterns with peptidic substrates. We provide a map of the binding sites of the murine proteases and suggest sequence motifs that we predict to be specific for mPR3 or mNE.

Influence of charge distribution at the active site surface on the substrate specificity of human neutrophil protease 3 and elastase. A kinetic and molecular modeling analysis.

The biological functions of human neutrophil protease 3 (Pr3) differ from those of neutrophil elastase despite their close structural and functional resemblance. Although both proteases are strongly cationic, their sequences differ mainly in the distribution of charged residues. We have used these differences in electrostatic surface potential in the vicinity of their active site to produce fluorescence resonance energy transfer (FRET) peptide substrates for investigating individual Pr3 subsites. The specificities of subsites S5 to S3' were investigated both kinetically and by molecular dynamic simulations. Subsites S2, S1', and S2' were the main definers of Pr3 specificity. Combinations of results for each subsite were used to deduce a consensus sequence that was complementary to the extended Pr3 active site and was not recognized by elastase. Similar sequences were identified in natural protein substrates such as NFkappaB and p21 that are specifically cleaved by Pr3. FRET peptides derived from these natural sequences were specifically hydrolyzed by Pr3 with specificity constants k(cat)/K(m) in the 10(6) m(-1) s(-1) range. The consensus Pr3 sequence may also be used to predict cleavage sites within putative protein targets like the proform of interleukin-18, or to develop specific Pr3 peptide-derived inhibitors, because none is available for further studies on the physiopathological function of this protease.

A novel locust (Schistocerca gregaria) serine protease inhibitor with a high affinity for neutrophil elastase.

We have purified to homogeneity two forms of a new serine protease inhibitor specific for elastase/chymotrypsin from the ovary gland of the desert locust Schistocerca gregaria. This protein, greglin, has 83 amino acid residues and bears putative phosphorylation sites. Amino acid sequence alignments revealed no homology with pacifastin insect inhibitors and only a distant relationship with Kazal-type inhibitors. This was confirmed by computer-based structural studies. The most closely related homologue is a putative gene product from Ciona intestinalis with which it shares 38% sequence homology. Greglin is a fast-acting and tight binding inhibitor of human neutrophil elastase (k(ass)=1.2x10(7) M(-1) x s(-1), K(i)=3.6 nM) and subtilisin. It also binds neutrophil cathepsin G, pancreatic elastase and chymotrypsin with a lower affinity (26 nM< or =K(i)< or =153 nM), but does not inhibit neutrophil protease 3 or pancreatic trypsin. The capacity of greglin to inhibit neutrophil elastase was not significantly affected by exposure to acetonitrile, high temperature (90 degrees C), low or high pH (2.5-11.0), N-chlorosuccinimide-mediated oxidation or the proteolytic enzymes trypsin, papain and pseudolysin from Pseudomonas aeruginosa. Greglin efficiently inhibits the neutrophil elastase activity of sputum supernatants from cystic fibrosis patients. Its biological function in the locust ovary gland is currently unknown, but its physicochemical properties suggest that it can be used as a template to design a new generation of highly resistant elastase inhibitors for treating inflammatory diseases.

Odorant binding and conformational dynamics in the odorant-binding protein.

In mammals, the olfactory epithelium secretes odorant-binding proteins (OBPs), which are lipocalins found freely dissolved in the mucus layer protecting the olfactory neurons. OBPs may act as passive transporters of predominantly hydrophobic odorant molecules across the aqueous mucus layer, or they may play a more active role in which the olfactory neuronal receptor recognizes the OBP-ligand complex. To better understand the molecular events accompanying the initial steps in the olfaction process, we have performed molecular dynamics studies of rat and pig OBPs with the odorant molecule thymol. These calculations provide an atomic level description of conformational changes and pathway intermediates that remain difficult to study directly. A series of eight independent molecular dynamics trajectories of rat OBP permitted the observation of a consensus pathway for ligand unbinding and the calculation of the potential of mean force (PMF) along this path. Titration microcalorimetry confirmed the specific binding of thymol to this protein with a strong hydrophobic component. In both rat and pig OBPs we observed lipocalin strand pair opening in the presence of ligand, consistent with potential roles of these proteins in olfactive receptor recognition.

Inspection of the binding sites of proteinase3 for the design of a highly specific substrate.

Proteinase3 (PR3) and human neutrophil elastase (HNE) are homologous proteases from the polymorphonuclear neutrophils and have been thought for a long time to have close enzymatic specificity. We have used molecular dynamics simulations to investigate and compare the interactions between different peptides and the two enzymes. The important role played especially by the C-terminal part of the peptides is confirmed. We provide a map of the subsites of PR3 and a description of the interaction scheme for six ligands. The main difference between HNE and PR3 concerns S2, S1', S2', and S3'. The recognition subsites in PR3 are interconnected; in particular, Lys99 participates to a hydrophobic (S4) and a polar (S2) pocket. On the basis of the simulations, we suggest that VADVKDR is a highly specific sequence for PR3; enzymatic assays confirm that it is cleaved by PR3 with a high specificity constant (k(cat)/K(m) = 3,400,000 M(-1) s(-1)) and not by HNE.

Cleavage of p21/WAF1/CIP1 by proteinase 3 modulates differentiation of a monocytic cell line. Molecular analysis of the cleavage site.

Proteinase 3 (PR3), also called myeloblastin, is involved in the control of myeloid cell growth, but the underlying molecular mechanisms have not been elucidated. In U937/PR3, stably transfected with PRCRSV/PR3 to overexpress PR3, PMA-induced p21 expression was significantly decreased as compared with control U937, and this phenomenon was reversed in the presence of the serine proteinase inhibitor, pefabloc. Conversely, when PR3 was inactivated by small interfering RNA, p21 protein was increased, and PMA-induced monocytic differentiation was potentiated. Mass spectrometry analysis identified Ala45 as the primary cleavage site on p21, and the recombinant mutated p21A45R, generated by site-directed mutagenesis and expressed in Escherichia coli, was resistant to in vitro PR3 cleavage. The U937 cells were then stably transfected with either PRCRSV/p21 or PRCRSV/p21A45R, to ectopically express wild type p21 or PR3-resistant p21, respectively. In U937/p21A45R treated with PS-341, a selective proteasome inhibitor, a significant decrease in the S phase and a blockade in the G0-G1 phase of cell cycle were observed when compared with U937/p21 or control U937. This suggested that both PR3 and the proteasome are efficiently involved in the proteolytic regulation of p21 expression in myeloid cells. Moreover, PMA-induced p21 expression was more pronounced in U937/p21A45R compared with U937/p21 and was concomitant with the morphological features of early differentiation. Our data demonstrated that p21 is one specific target of PR3 and that PR3-mediated p21 cleavage prevents monocytic differentiation.