Molecular modeling of the reductase domain to elucidate the reaction mechanism of reduction of peptidyl thioester into its corresponding alcohol in non-ribosomal peptide synthetases

BackgroundNonribosomal peptide synthetases (NRPSs) are multienzymatic, multidomain megasynthases involved in the biosynthesis of pharmaceutically important nonribosomal peptides. The peptaibol synthetase from Trichoderma virens (TPS) is an important member of the NRPS family that exhibits antifungal properties. The majority of the NRPSs terminate peptide synthesis with the thioesterase (TE) domain, which either hydrolyzes the thioester linkage, releasing the free peptic acid, or catalyzes the intramolecular macrocyclization to produce a macrolactone product. TPS is an important NRPS that does not encompass a TE domain, but rather a reductase domain (R domain) to release the mature peptide product reductively with the aid of a NADPH cofactor. However, the catalytic mechanism of the reductase domain has not yet been elucidated.ResultsWe present here a three-dimensional (3D) model of the reductase domain based on the crystal structure of vestitone reductase (VR). VR belongs to the short-chain dehydrogenase/reductase (SDR) superfamily and is responsible for the nicotinamide dinucleotide phosphate (NADPH)-dependent reduction of the substrate into its corresponding secondary alcohol product. The binding sites of the probable linear substrates, alamethicin, trichotoxin, antiamoebin I, chrysopermin C and gramicidin, were identified within the modeled R domain using multiple docking approaches. The docking results of the ligand in the active site of the R domain showed that reductase side chains have a high affinity towards ligand binding, while the thioester oxygen of each substrate forms a hydrogen bond with the OH group of Tyr176 and the thiol group of the substrate is closer to the Glu220. The modeling and docking studies revealed the reaction mechanism of reduction of thioester into a primary alcohol.ConclusionPeptaibol biosynthesis incorporates a single R domain, which appears to catalyze the four-electron reduction reaction of a peptidyl carrier protein (PCP)-bound peptide to its corresponding primary alcohol. Analysis of R domains present in the non-redundant (nr) database of the NCBI showed that the R domain always resides in the last NRPS module and is involved in either a two or four-electron reduction reaction.

[1]  D. Schomburg,et al.  Prediction of protein three-dimensional structures in insertion and deletion regions: a procedure for searching data bases of representative protein fragments using geometric scoring criteria. , 1995, Journal of molecular biology.

[2]  F. Stevens.  Efficient recognition of protein fold at low sequence identity by conservative application of Psi‐BLAST: validation , 2005, Journal of molecular recognition : JMR.

[3]  D. Ehmann,et al.  Lysine biosynthesis in Saccharomyces cerevisiae: mechanism of alpha-aminoadipate reductase (Lys2) involves posttranslational phosphopantetheinylation by Lys5. , 1999, Biochemistry.

[4]  Ian W. Davis,et al.  Structure validation by Cα geometry: ϕ,ψ and Cβ deviation , 2003, Proteins.

[5]  J. Madura,et al.  Structure modeling, ligand binding, and binding affinity calculation (LR‐MM‐PBSA) of human heparanase for inhibition and drug design , 2006, Proteins.

[6]  H Edelsbrunner,et al.  Analytical shape computation of macromolecules: II. Inaccessible cavities in proteins , 1998, Proteins.

[7]  M. Sansom.  Alamethicin and related peptaibols — model ion channels , 2004, European Biophysics Journal.

[8]  Michael J E Sternberg,et al.  Exploring the extremes of sequence/structure space with ensemble fold recognition in the program Phyre , 2008, Proteins.

[9]  Tudor I. Oprea,et al.  Property distribution of drug-related chemical databases* , 2000, J. Comput. Aided Mol. Des..

[10]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[11]  L. Tabernero,et al.  Substrate-induced closure of the flap domain in the ternary complex structures provides insights into the mechanism of catalysis by 3-hydroxy-3-methylglutaryl-CoA reductase. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[12]  Jeffry D Madura,et al.  Dopamine transporter comparative molecular modeling and binding site prediction using the LeuTAa leucine transporter as a template , 2008, Proteins.

[13]  W Patrick Walters,et al.  A detailed comparison of current docking and scoring methods on systems of pharmaceutical relevance , 2004, Proteins.

[14]  Mohamed A. Marahiel,et al.  Modular Peptide Synthetases Involved in Nonribosomal Peptide Synthesis. , 1997, Chemical reviews.

[15]  F. Kopp,et al.  Peptide macrocyclization: the reductase of the nostocyclopeptide synthetase triggers the self-assembly of a macrocyclic imine. , 2006, Journal of the American Chemical Society.

[16]  H. Wenschuh,et al.  The NMR solution structure of the ion channel peptaibol chrysospermin C bound to dodecylphosphocholine micelles. , 2000, European journal of biochemistry.

[17]  Peer Bork,et al.  SMART: a web-based tool for the study of genetically mobile domains , 2000, Nucleic Acids Res..

[18]  M. Levitt.  Accurate modeling of protein conformation by automatic segment matching. , 1992, Journal of molecular biology.

[19]  P. Leadlay,et al.  The biosynthetic gene cluster for the polyketide immunosuppressant rapamycin. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[20]  D. Eisenberg,et al.  Assessment of protein models with three-dimensional profiles , 1992, Nature.

[21]  Jie Liang,et al.  CASTp: Computed Atlas of Surface Topography of proteins , 2003, Nucleic Acids Res..

[22]  Thomas A. Halgren.  Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 , 1996, J. Comput. Chem..

[23]  M. Marahiel,et al.  Multimodular biocatalysts for natural product assembly , 2001, Naturwissenschaften.

[24]  J M Blaney,et al.  A geometric approach to macromolecule-ligand interactions. , 1982, Journal of molecular biology.

[25]  B. Shen,et al.  Biosynthesis of hybrid peptide-polyketide natural products. , 2001, Current opinion in drug discovery & development.

[26]  Liam J. McGuffin,et al.  The PSIPRED protein structure prediction server , 2000, Bioinform..

[27]  Jie Liang,et al.  CASTp: computed atlas of surface topography of proteins with structural and topographical mapping of functionally annotated residues , 2006, Nucleic Acids Res..

[28]  M. Marahiel,et al.  The tyrocidine biosynthesis operon of Bacillus brevis: complete nucleotide sequence and biochemical characterization of functional internal adenylation domains , 1997, Journal of bacteriology.

[29]  M. Marahiel,et al.  Synthesis of linear gramicidin requires the cooperation of two independent reductases. , 2005, Biochemistry.

[30]  Andrei L Lomize,et al.  Bmc Structural Biology , 2022 .

[31]  J. Thompson,et al.  The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. , 1997, Nucleic acids research.

[32]  W. Gerwick,et al.  Lyngbyatoxin biosynthesis: sequence of biosynthetic gene cluster and identification of a novel aromatic prenyltransferase. , 2004, Journal of the American Chemical Society.

[33]  E. Maser,et al.  Carbonyl Reductases and Pluripotent Hydroxysteroid Dehydrogenases of the Short-chain Dehydrogenase/reductase Superfamily , 2007, Drug metabolism reviews.

[34]  C. Walsh,et al.  The lyngbyatoxin biosynthetic assembly line: chain release by four-electron reduction of a dipeptidyl thioester to the corresponding alcohol. , 2007, Journal of the American Chemical Society.

[35]  T. Pollard,et al.  Annual review of biophysics and biomolecular structure , 1992 .

[36]  R. Durbin,et al.  Pfam: A comprehensive database of protein domain families based on seed alignments , 1997, Proteins.

[37]  S Subramaniam,et al.  Analytical shape computation of macromolecules: I. molecular area and volume through alpha shape , 1998, Proteins.

[38]  R. Müller,et al.  In vitro reconstitution of the myxochelin biosynthetic machinery of Stigmatella aurantiaca Sg a15: Biochemical characterization of a reductive release mechanism from nonribosomal peptide synthetases , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[39]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

[40]  D. Hopwood,et al.  Genetic Contributions to Understanding Polyketide Synthases. , 1997, Chemical reviews.

[41]  Akbar Nayeem,et al.  A comparative study of available software for high‐accuracy homology modeling: From sequence alignments to structural models , 2006, Protein science : a publication of the Protein Society.

[42]  R. Campbell,et al.  The first structure of UDP-glucose dehydrogenase reveals the catalytic residues necessary for the two-fold oxidation. , 2000, Biochemistry.

[43]  D. Marsh.  Peptide models for membrane channels. , 1996, The Biochemical journal.

[44]  P E Bourne,et al.  Protein structure alignment by incremental combinatorial extension (CE) of the optimal path. , 1998, Protein engineering.

[45]  H. Blöcker,et al.  Novel features in a combined polyketide synthase/non-ribosomal peptide synthetase: the myxalamid biosynthetic gene cluster of the myxobacterium Stigmatella aurantiaca Sga15. , 2001, Chemistry & biology.

[46]  Noriaki Hirayama,et al.  Ph4Dock: pharmacophore-based protein-ligand docking. , 2004, Journal of medicinal chemistry.

[47]  Gitanjali Yadav,et al.  NRPS-PKS: a knowledge-based resource for analysis of NRPS/PKS megasynthases , 2004, Nucleic Acids Res..

[48]  R. Dixon,et al.  Crystal structure of vestitone reductase from alfalfa (Medicago sativa L.). , 2007, Journal of molecular biology.

[49]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .

[50]  Stephen H. Bryant,et al.  CD-Search: protein domain annotations on the fly , 2004, Nucleic Acids Res..

[51]  J. Chugh,et al.  Model for a helical bundle channel based on the high-resolution crystal structure of trichotoxin_A50E. , 2002, Biochemistry.

[52]  C. Walsh,et al.  The parallel and convergent universes of polyketide synthases and nonribosomal peptide synthetases. , 1999, Chemistry & biology.

[53]  B. Wallace,et al.  Solution NMR studies of antiamoebin, a membrane channel-forming polypeptide. , 2003, Biophysical journal.

[54]  D. Ebbole,et al.  Identification of Peptaibols from Trichoderma virens and Cloning of a Peptaibol Synthetase* , 2002, The Journal of Biological Chemistry.

[55]  C. Walsh,et al.  Harnessing the biosynthetic code: combinations, permutations, and mutations. , 1998, Science.

[56]  R. Coppel,et al.  Identification of a peptide synthetase involved in the biosynthesis of glycopeptidolipids of Mycobacterium smegmatis , 1999, Molecular microbiology.

[57]  R. Müller,et al.  Myxochelin biosynthesis: direct evidence for two- and four-electron reduction of a carrier protein-bound thioester. , 2008, Journal of the American Chemical Society.

[58]  J R Jacobsen,et al.  Tolerance and specificity of polyketide synthases. , 1999, Annual review of biochemistry.

[59]  P. Leadlay,et al.  Organisation of the biosynthetic gene cluster for rapamycin in Streptomyces hygroscopicus: analysis of genes flanking the polyketide synthase. , 1996, Gene.

[60]  Natasja Brooijmans,et al.  Molecular recognition and docking algorithms. , 2003, Annual review of biophysics and biomolecular structure.

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

[62]  A. Pospiech,et al.  Two multifunctional peptide synthetases and an O-methyltransferase are involved in the biosynthesis of the DNA-binding antibiotic and antitumour agent saframycin Mx1 from Myxococcus xanthus. , 1996, Microbiology.

[63]  Xiaofeng Wu,et al.  Structure and Biosynthesis of the BT Peptide Antibiotic from Brevibacillus texasporus , 2005, Applied and Environmental Microbiology.

[64]  D. Eisenberg,et al.  A method to identify protein sequences that fold into a known three-dimensional structure. , 1991, Science.