Molecular insight into pseudolysin inhibition using the MM-PBSA and LIE methods.

Pseudolysin, the extracellullar elastase of Pseudomonas aeruginosa (EC: 3.4.24.26) plays an important role in the pathogenesis of P. aeruginosa infections. In the present study, molecular dynamics simulations and theoretical affinity predictions were used to gain molecular insight into pseudolysin inhibition. Four low molecular weight inhibitors were docked at their putative binding sites and molecular dynamics (MD) simulations were performed for 5.0 ns, and the free energy of binding was calculated by the linear interaction energy method. The number and the contact surface area of stabilizing hydrophobic, aromatic, and hydrogen bonding interactions appears to reflect the affinity differences between the inhibitors. The proteinaceous inhibitor, Streptomyces metalloproteinase inhibitor (SMPI) was docked in three different binding positions and MD simulations were performed for 3.0 ns. The MD trajectories were used for molecular mechanics-Poisson-Boltzmann surface area analysis of the three binding positions. Computational alanine scanning of the average pseudolysin-SMPI complexes after MD revealed residues at the pseudolysin-SMPI interface giving the main contribution to the free energy of binding. The calculations indicated that SMPI interacts with pseudolysin via the rigid active site loop, but that also contact sites outside this loop contribute significantly to the free energy of association.

[1]  M. Sanner,et al.  Reduced surface: an efficient way to compute molecular surfaces. , 1996, Biopolymers.

[2]  P A Kollman,et al.  What determines the van der Waals coefficient β in the LIE (linear interaction energy) method to estimate binding free energies using molecular dynamics simulations? , 1999, Proteins.

[3]  Brian W. Matthews,et al.  Structural basis of the action of thermolysin and related zinc peptidases , 1988 .

[4]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[5]  M. Kainosho,et al.  Mutational analysis of the reactive site loop of Streptomyces metalloproteinase inhibitor, SMPI. , 1999, Journal of biochemistry.

[6]  John Marelius,et al.  Calculation of Ligand Binding Free Energies from Molecular Dynamics Simulations , 1998 .

[7]  D. Johnson,et al.  Inactivation of human bronchial mucosal proteinase inhibitor by Pseudomonas aeruginosa elastase. , 1982, The American review of respiratory disease.

[8]  A. Szabó,et al.  Langevin modes of macromolecules , 1986 .

[9]  T. Clackson,et al.  Structural and functional analysis of the 1:1 growth hormone:receptor complex reveals the molecular basis for receptor affinity. , 1998, Journal of molecular biology.

[10]  K. Hiraga,et al.  Resynthesis of reactive site peptide bond and temporary inhibition of Streptomyces metalloproteinase inhibitor. , 1997, Journal of biochemistry.

[11]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[12]  P A Kollman,et al.  Calculation and prediction of binding free energies for the matrix metalloproteinases. , 2000, Journal of medicinal chemistry.

[13]  J. Aqvist,et al.  A new method for predicting binding affinity in computer-aided drug design. , 1994, Protein engineering.

[14]  B. Honig,et al.  Classical electrostatics in biology and chemistry. , 1995, Science.

[15]  J M Thornton,et al.  LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. , 1995, Protein engineering.

[16]  P. Kollman,et al.  Application of RESP charges to calculate conformational energies, hydrogen bond energies, and free energies of solvation , 1993 .

[17]  H M Holden,et al.  The binding of L-valyl-L-tryptophan to crystalline thermolysin illustrates the mode of interaction of a product of peptide hydrolysis. , 1989, The Journal of biological chemistry.

[18]  E. Kessler,et al.  Inhibition by phosphoramidon of Pseudomonas aeruginosa elastase injected intracorneally in rabbit eyes. , 1984, Current eye research.

[19]  Xiaojie Xu,et al.  Predictions of Binding of a Diverse Set of Ligands to Gelatinase-A by a Combination of Molecular Dynamics and Continuum Solvent Models , 2002 .

[20]  J. Powers,et al.  Active site directed irreversible inhibition of thermolysin. , 1978, Biochemistry.

[21]  M. Kainosho,et al.  Elucidation of the mode of interaction of thermolysin with a proteinaceous metalloproteinase inhibitor, SMPI, based on a model complex structure and a structural dynamics analysis. , 1998, Journal of molecular biology.

[22]  C. Hodge,et al.  Fitting an inhibitor into the active site of thermolysin: A molecular dynamics case study , 1996, Proteins.

[23]  H M Holden,et al.  Structures of two thermolysin-inhibitor complexes that differ by a single hydrogen bond. , 1987, Science.

[24]  M. Tashiro,et al.  Identification of reactive site of a proteinaceous metalloproteinase inhibitor from Streptomyces nigrescens TK-23. , 1997, Journal of biochemistry.

[25]  R. Stein,et al.  Mechanistic studies on the inhibition of thermolysin by a peptide hydroxamic acid , 1992 .

[26]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[27]  Xiaojie Xu,et al.  Binding Affinities for a Series of Selective Inhibitors of Gelatinase-A Using Molecular Dynamics with a Linear Interaction Energy Approach , 2001 .

[28]  D. Case,et al.  Langevin modes of macromolecules: Applications to crambin and DNA hexamers , 1990, Biopolymers.

[29]  B. Nordén,et al.  Computational modelling of inhibitor binding to human thrombin. , 2001, European journal of pharmaceutical sciences : official journal of the European Federation for Pharmaceutical Sciences.

[30]  K. Botzenhart,et al.  [Extracellular toxins of Pseudomonas aeruginosa. II. Effect of two proteases on human immunoglobulins IgG, IgA and secretory IgA (author's transl)]. , 1981, Zentralblatt fur Bakteriologie. 1. Abt. Originale. A: Medizinische Mikrobiologie, Infektionskrankheiten und Parasitologie.

[31]  E. Kessler,et al.  In Vitro Inhibition of Pseudomonas aeruginosa Elastase by Metal-Chelating Peptide Derivatives , 1982, Infection and immunity.

[32]  Y. Komori,et al.  Hemorrhagic activity and muscle damaging effect of Pseudomonas aeruginosa metalloproteinase (elastase). , 2001, Toxicon.

[33]  K. Morihara,et al.  Pseudolysin and other pathogen endopeptidases of thermolysin family. , 1995, Methods in enzymology.

[34]  Emil Alexov,et al.  Rapid grid‐based construction of the molecular surface and the use of induced surface charge to calculate reaction field energies: Applications to the molecular systems and geometric objects , 2002, J. Comput. Chem..

[35]  D. J. Price,et al.  Assessing scoring functions for protein-ligand interactions. , 2004, Journal of medicinal chemistry.

[36]  E. J. Miller,et al.  Specific cleavage of human type III and IV collagens by Pseudomonas aeruginosa elastase , 1986, Infection and immunity.

[37]  K. Flaherty,et al.  Three-dimensional structure of the elastase of Pseudomonas aeruginosa at 1.5-A resolution. , 1993, The Journal of biological chemistry.

[38]  P A Kollman,et al.  Molecular dynamics and continuum solvent studies of the stability of polyG-polyC and polyA-polyT DNA duplexes in solution. , 1998, Journal of biomolecular structure & dynamics.

[39]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[40]  R. Gray,et al.  Inhibition of Pseudomonas aeruginosa elastase and Pseudomonas keratitis using a thiol-based peptide , 1990, Antimicrobial Agents and Chemotherapy.

[41]  B. Honig,et al.  Calculation of the total electrostatic energy of a macromolecular system: Solvation energies, binding energies, and conformational analysis , 1988, Proteins.

[42]  J. Åqvist,et al.  Calculation of absolute binding free energies for charged ligands and effects of long‐range electrostatic interactions , 1996 .

[43]  J. Alcorn,et al.  Pseudomonas aeruginosa elastase degrades surfactant proteins A and D. , 2003, American journal of respiratory cell and molecular biology.

[44]  S. Murao,et al.  Purification and properties of a proteinaceous metallo-proteinase inhibitor from Streptomyces nigrescens TK-23. , 1979, Biochimica et biophysica acta.

[45]  F. Jørgensen,et al.  Structural Differences of Matrix Metalloproteinases. Homology Modeling and Energy Minimization of Enzyme-Substrate Complexes , 2000, Journal of biomolecular structure & dynamics.

[46]  A. Schmidtchen,et al.  Elastase-producing Pseudomonas aeruginosa degrade plasma proteins and extracellular products of human skin and fibroblasts, and inhibit fibroblast growth. , 2003, Microbial pathogenesis.

[47]  I. A. Holder,et al.  Experimental studies of the pathogenesis of infections owing to Pseudomonas aeruginosa: elastase, an IgG protease. , 1984, Canadian journal of microbiology.

[48]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[49]  Ingebrigt Sylte,et al.  The Protein-protein Interactions Between SMPI and Thermolysin Studied by Molecular Dynamics and MM/PBSA Calculations , 2005, Journal of biomolecular structure & dynamics.

[50]  E. Kessler,et al.  Specific inhibitors of Pseudomonas aeruginosa elastase as potential drugs for the treatment of Pseudomonas keratitis. , 1987, Antibiotics and chemotherapy.

[51]  A. Berger,et al.  On the size of the active site in proteases. I. Papain. , 1967, Biochemical and biophysical research communications.

[52]  J. Åqvist,et al.  Ligand binding affinities from MD simulations. , 2002, Accounts of chemical research.

[53]  J. Powers,et al.  Pseudomonas aeruginosa elastase. Development of a new substrate, inhibitors, and an affinity ligand. , 1980, Journal of Biological Chemistry.

[54]  D E Tronrud,et al.  Crystallographic structural analysis of phosphoramidates as inhibitors and transition-state analogs of thermolysin. , 1989, European journal of biochemistry.

[55]  Johan Åqvist,et al.  Ligand binding affinity prediction by linear interaction energy methods , 1998, J. Comput. Aided Mol. Des..

[56]  D. Baker,et al.  A simple physical model for binding energy hot spots in protein–protein complexes , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[57]  H. Tsuzuki,et al.  Protease and elastase of Pseudomonas aeruginosa: inactivation of human plasma alpha 1-proteinase inhibitor , 1979, Infection and immunity.

[58]  H. Berendsen,et al.  ALGORITHMS FOR MACROMOLECULAR DYNAMICS AND CONSTRAINT DYNAMICS , 1977 .

[59]  Ruben Abagyan,et al.  ICM—A new method for protein modeling and design: Applications to docking and structure prediction from the distorted native conformation , 1994, J. Comput. Chem..

[60]  J. A. Hobden Pseudomonas aeruginosa proteases and corneal virulence. , 2002, DNA and cell biology.

[61]  K. Sharp,et al.  Accurate Calculation of Hydration Free Energies Using Macroscopic Solvent Models , 1994 .

[62]  Jaime Prilusky,et al.  Automated analysis of interatomic contacts in proteins , 1999, Bioinform..

[63]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[64]  K Nadassy,et al.  Analysis of zinc binding sites in protein crystal structures , 1998, Protein science : a publication of the Protein Society.

[65]  David E. Kim,et al.  Computational Alanine Scanning of Protein-Protein Interfaces , 2004, Science's STKE.

[66]  M. Swindells,et al.  NMR structure of the Streptomyces metalloproteinase inhibitor, SMPI, isolated from Streptomyces nigrescens TK-23: another example of an ancestral beta gamma-crystallin precursor structure. , 1998, Journal of molecular biology.

[67]  W. L. Jorgensen Revised TIPS for simulations of liquid water and aqueous solutions , 1982 .