Molecular docking between the RNA polymerase of the Moniliophthora perniciosa mitochondrial plasmid and Rifampicin produces a highly stable complex
© Andrade et al.; licensee BioMed Central Ltd. 2013
Received: 22 August 2012
Accepted: 21 February 2013
Published: 26 February 2013
Moniliophthora perniciosa (Stahel) Aime & Phillips-Mora is the causal agent of witches’ broom disease (WBD) in cacao (Theobroma cacao). When the mitochondrial genome of this fungus had been completely sequenced, an integrated linear-type plasmid that encodes viral-like RNA polymerases was found. The structure of this polymerase was previously constructed using a homology modeling approach.
Using a virtual screening process, accessing the Kegg, PubChem and ZINC databases, we selected the eight most probable macrocyclic polymerase inhibitors to test against M. perniciosa RNA polymerase (RPO). AutoDock Vina was used to perform docking calculations for each molecule. This software returned affinity energy values for several ligand conformations. Subsequently, we used PyMOL 1.4 and Ligand Scout 3.1 to check the stereochemistry of chiral carbons, substructure, superstructure, number of rotatable bonds, number of rings, number of donor groups, and hydrogen bond receptors.
On the basis of this evidence we selected Rifampicin, a bacterial RNA polymerase inhibitor, and then AMBER 12 was used to simulate the behavior of the RPO-Rifampicin complex after a set of 5000 ps and up to 300 K in water. This calculation returned a graph of potential energy against simulation time and showed that the ligand remained inside the active site after the simulation was complete, with an average energy of -15 x 102 Kcal/Mol.
The results indicate that Rifampicin could be a good inhibitor for testing in vitro and in vivo against M. perniciosa.
KeywordsMoniliophthora perniciosa RNA polymerase Rifampicin Docking MM/PBSA
Moniliophthora perniciosa (Stahel) Aime & Phillips-Mora is the causal agent of witches’ broom disease (WBD) in cacao (Theobroma cacao). When the mitochondrial genome (http://www.ncbi.nlm.nih.gov/genomes/GenomesGroup.cgi?taxid=4751opt=organelle/) of this fungus had been completely sequenced, an integrated linear-type plasmid that encodes viral-like DNA and RNA polymerases was found .
RNA polymerase (RNAP) is responsible for transcribing DNA and is the direct or indirect target of most regulators of transcription . The enzyme from the M. perniciosa mitochondrial plasmid (RPO) is a 766 amino acid DNA-dependent RNA polymerase belonging to the single chain family of polymerases, which occur in viruses and cellular organelles . Its active site is distributed between two domains, Palm (Asp457 and Asp695) and Fingers (Tyr537 and Lys529), involved in transcription [3–5]. The mechanism of transcription by this enzyme shares several similarities with other multichain RNA polymerases [3, 6], so it could have inhibitors in common with other polymerases.
Rifampicin is a macrocyclic molecule the ansamycin family. It contains a methyl-piperazinyliminomethyl side chain at position 3, a cyclopentyl-piperazinyliminomethyl side chain at position 3, and a cyclic spiro-piperidyl side chain at positions 3 and 4 . This drug has been used since 1968 to combat Mycobacterium tuberculosis but is considered a broad-spectrum antibiotic [8, 9]. It has a high capacity to bind and inhibit DNA-dependent RNA polymerases (RNAP) from bacteria through its specific interaction is with the polymerase β subunit . The essential catalytic core of RNAP is evolutionarily conserved among all cellular organisms .
Virtual structure-based screening has become prominent in drug discovery, using protein targets [11, 12]. Several free ligand databases are widely available today. Searching for molecules that can complex with target proteins can be done either by keywords (e.g. Kegg and PubChem databases) or by using a structure-activity relationship, available in the Zinc  and PubChem databases.
One of the most important techniques for receptor-based drug design is molecular docking . Using crystallographic or modeled protein structures, molecular docking is often employed to screen compound libraries and to predict the conformation of a protein-ligand complex and calculate its affinity energy . In general, docking programs such as AutoDock Vina  generate multiple protein-ligand conformations by sampling the ligand’s probable conformations in the binding pocket of the target protein, using flexible ligand-rigid receptor docking [11, 14]. Scoring functions are used for docking calculations by these programs in an attempt to approximate the standard chemical potentials of the system . AutoDock Vina uses a force-field-based scoring function approach to estimate binding affinities by calculating the non-bonded interactions based on traditional force fields, identify the correct binding position of a ligand, and rank ligands by their predicted binding affinities [11, 14]. On the other hand, the problems of molecular docking as a screening tool have been widely discussed: the scoring functions are in general inaccurate and neglect solvent-related terms, and protein flexibility is ignored . Coupled Molecular Docking and Molecular Dynamics is a good way to solve this problem because it can treat both proteins and ligands in a flexible manner, allowing the binding site to be relaxed around the ligand [11, 12, 15].
Molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) combines molecular mechanics energy and implicit solvation models and is more rigorous than most empirical or knowledge-based scoring functions [11, 12]. It allows for rigorous free-energy decomposition into contributions originating from different groups of atoms or types of interaction [11, 16]. In the MM/PB-SA method the free energy is calculated using snapshots of solute molecules obtained from explicit-solvent MD simulation .
The aim of this study was to search a series of molecules likely to form complexes with the RPO from M. perniciosa, available at the Kegg, PubChem and Zinc databases, and to select a potential inhibitor using a coupled Molecular Docking and Molecular Dynamics - MM / PBSA - approach.
We searched initially by keywords and by nucleoside molecules in the Kegg (http://www.genome.jp/kegg/) and PubChem (http://pubchem.ncbi.nlm.nih.gov/) databases. Only molecules described as inhibitors of RNA polymerases were selected from the outputs. All 2D structures were copied in Similes format for comparison with other Zinc database (http://zinc.docking.org/) molecules; this increases the number of molecules that can be used for docking, starting from a known ligand. The 3D structures of these molecules were downloaded in mol2 and pdb formats for use in Virtual Screening, which was carried out by Molecular Docking and Molecular Dynamics. Following a protocol described by Irwin and Schoichet (2005) , molecules obtained from the Zinc database were selected for comparison, given 95-99% similarity, with the structures found in the Kegg and PubChem databases. In addition, selected structures were downloaded in mol2 and pdb formats for subsequent Docking Studies and Molecular Dynamics.
The structures downloaded from the Kegg, Pubchem and Zinc databases were first checked in PyMOL 1.4  to evaluate the presence or absence of hydrogen, the stereochemistry of chiral carbons, substructure, superstructure, number of rotatable bonds, number of rings and number of hydrogen acceptors and donors.
The ligand and receptor (RPO) molecules were prepared in AutoDockTools 1.5.6 . All polar hydrogens were added to the receptor and Kollman United Atomic Charges were computed. For all ligands we added polar hydrogens and computed the Gasteiger charges. The grid definition, adjusted to the RPO active site, was set up manually by following the recommendations of the program manual [14, 18]. The structures of the ligand and receptor were then saved in pqbqt format to be used for docking calculations. AutoDock Vina was used to perform Docking Scoring for each ligand-receptor complex . Before running each Docking calculation, a configuration file was generated with information about grid size and coordinates and indicating the ligand and receptor files. The reports for each calculation were analyzed to obtain affinity energy (Kcal/mol) values for each ligand conformation in its respective complex. In addition, we used PyMOL 1.4 to verify the number of hydrogen bonds and non-covalent interactions between each ligand conformation and the catalytic residues of RPO that are involved in the recognition and polymerization mechanisms. In order to optimize the choice of an ideal complex we selected just one ligand that fit best in the RPO active site, considering all stereochemical aspects previously evaluated and the free-energy results.
Molecular dynamics of complex
In this study we used the MM/PBSA protocol to calculate affinity and stability in the ligand-receptor RPO complex interaction, using the package Amber 12 . Initially, we used the Antechamber program to make the ff99 force field recognize the types of atom in both ligand and receptor and avoid errors during the calculations.
tLEaP was used to neutralize charges (ff99 force field) and the RPO–ligand complex was immersed in a rectangular box of TIP3P water molecules. Following the protocol we used Sander to carry out a Molecular Dynamics (MD) equilibrium, restricted to a region of the protein that contains the active site (amino acids 457-695), according to the following parameters: 1000 cycles of steepest descent and 1000 cycles of conjugate gradient minimization, heating MD for 200 picoseconds (ps), density equilibrium for 200 ps, followed by Equilibrium Dynamics for 600 ps at constant pressure and 300 K. After the system equilibrated we followed this with the MM/PBSA protocol [11, 19, 20]. We then simulated a total of 4000 ps production steps of Molecular Dynamics, divided into four sets of 1000 ps (prod1, prod2, prod3 and prod4), saving the coordinates every 10 ps. Furthermore, we used the mm_pbsa.pl script to extract snapshots (without the water) from our production runs and obtain its trajectories. In addition, we checked the stability of the complex by plotting Potential Energy × Time (ps) and RMSD × Time (ps) graphs from all simulation trajectories. As a final step, we used the ambpdb command to generate a pdb file of the complex after the last stage of the Molecular Dynamics, and this structure was analyzed in PyMOL 1.4 to verify whether the ligand remained in the active site after the process was complete. In addition, Ligand Scout 3.1  was used to generate 2D and 3D interaction maps of RPO-Rifampicin, presenting Hydrogen Bond Acceptors and Donors (HBA and HBD) and all hydrophobic interactions in the active site.
Results and discussion
Structures and binding energies of RPO complexes from AutoDock Vina
Potential RPO inhibitors selected from the Kegg, PubChem and Zinc databases used in docking studies
The use of RNA polymerases as molecular targets for virtual screening is not restricted to prokaryotes. An RNA-dependent RNA polymerase (RdRp) is an attractive target for anti-HCV agents . However, we found no articles that dealt specifically with the use of inhibitors of fungal polymerases or polymerases encoded by mitochondrial genes. On the other hand, many authors report that all cellular RNA polymerases are relatively conserved in amino acid sequence and catalytic mechanism [2, 7, 9, 23], We can therefore understand why we generally find the same class (macrocyclic) of RNA polymerase inhibitors acting on different groups of organisms. In addition, Rifampicin can probably act in vitro and in vivo, inhibiting mitochondrial transcription by RPO and thus blocking the mitochondrial metabolism of M. perniciosa.
Molecular Dynamics MM/PBSA of RPO-Rifampicin complex
The pdb of the complex, generated after molecular dynamics, showed that Rifampicin remains within the active site of RPO after 5000 ps of simulation. Furthermore, we can infer that this simulation time was sufficient to show that Rifampicin could block M. perniciosa RPO activity.
In this article we describe the selection of one potential inhibitor of the M. perniciosa mitochondrial plasmid RNA polymerase from among eight molecules found in public databases, using a virtual screening approach. Rifampicin is a bacterial RNAPs inhibitor; nevertheless it forms a very stable complex with RPO, perhaps because this type of enzyme is highly conserved among organisms. Rifampicin forms complexes with exactly those amino acids in the active site pocket that are involved in transcription by RPO. This integration remained stable throughout the 5000 ps of Molecular Dynamics.
In a further step, we could analyze different mechanisms of Biomolecular Simulation to describe the mechanism by which Rifampicin inhibits the RPO, and whether Rifampicin will acquire different conformations within this process that can effectively interact within the active site of this enzyme.
We cannot discard use of the other ligands described in Table 1 as potential inhibitors of RPO or of others described as macrocyclic. Rifampicin was selected because it forms the best ligand-complex interaction. In future work we will perform biochemical tests in vitro and in vivo to verify whether our selected inhibitor can act effectively against M. perniciosa by making RPO-dependent transcription unfeasible.
We thank the Graduate program in Biotechnology (PPGBiotec / UEFS-Fiocruz) for logistical support during the execution of this work; The Bioinformatics Laboratory (LAPEM–UEFS) for the infrastructure available to perform calculations; and The State University of Feira de Santana.
- Formighieri EF, Tiburcio RA, Armas ED, Medrano FJ, Shimo H, Carels N, Góes-Neto A, Cotomacci C, Carazzolle MF, Sardinha-Pinto N, Thomazella DP, Rincones J, Digiampietri L, Carraro DM, Azeredo-Espin AM, Reis SF, Deckmann AC, Gramacho K, Gonçalves MS, Moura Neto JP, Barbosa LV, Meinhardt LW, Cascardo JC, Pereira GA: The mitochondrial genome of the phytopathogenic basidiomycete Moniliophthora perniciosa is 109 kb in size and contains a stable integrated plasmid. Mycol Res. 2008, 112 (10): 1136-1152. 10.1016/j.mycres.2008.04.014.View ArticlePubMedGoogle Scholar
- Tuske S, Sarafianos SG, Wang X, Hudson B, Sineva E, Mukhopadhyay J, Birktoft JJ, Leroy O, Ismail S, Clark AD, Dharia C, Napoli A, Laptenko O, Lee J, Borukhov S, Ebright RH, Arnold E: Inhibition of bacterial RNA polymerase by streptolydigin: stabilization of a straight-bridge-helix active-center conformation. Cell. 2005, 122 (4): 541-552. 10.1016/j.cell.2005.07.017.PubMed CentralView ArticlePubMedGoogle Scholar
- Andrade BS, Taranto AG, Góes-Neto A, Duarte AA: Comparative modeling of DNA and RNA polymerases from Moniliophthora perniciosa mitochondrial plasmid. Theor Biol Med Model. 2009, 6: 22-10.1186/1742-4682-6-22.PubMed CentralView ArticlePubMedGoogle Scholar
- Cheetham GMT, Jeruzalmi D, Steitz T: Structural basis for initiation of transcription from an RNA polymerase-promoter complex. Nature. 1999, 399 (6731): 80-83. 10.1038/19999.View ArticlePubMedGoogle Scholar
- Bonner G, Patra D, Lafer EM, Sousa R: Mutations in T7 RNA polymerase that support the proposal for a common polymerase active site structure. EMBO J. 1992, 11 (10): 3767-3775.PubMed CentralPubMedGoogle Scholar
- Sousa R, Chung YJ, Rose JP, Wang BC: Structure of bacteriophage T7 RNA polymerase at 3.3 A resolution. Nature. 1993, 364 (6438): 593-599. 10.1038/364593a0.View ArticlePubMedGoogle Scholar
- Ho MX, Hudson BP, Das K, Arnold E, Ebright RH: Structures of RNA polymerase-antibiotic complexes. Curr Opin Struct Biol. 2009, 19 (6): 715-723. 10.1016/j.sbi.2009.10.010.PubMed CentralView ArticlePubMedGoogle Scholar
- Mick V, Domínguez MA, Tubau F, Liñares J, Pujol M, Martín R: Molecular characterization of resistance to Rifampicin in an emerging hospital-associated Methicillin-resistant Staphylococcus aureus clone ST228, Spain. BMC Microbiol. 2010, 10: 68-10.1186/1471-2180-10-68.PubMed CentralView ArticlePubMedGoogle Scholar
- Campbell EA, Korzheva N, Mustaev A, Murakami K, Nair S, Goldfarb A, Darst SA: Structural mechanism for rifampicin inhibition of bacterial RNA polymerase. Cell. 2001, 104 (6): 901-912. 10.1016/S0092-8674(01)00286-0.View ArticlePubMedGoogle Scholar
- Aboshkiwa M, Rowland G, Coleman G: Nucleotide sequence of the Staphylococcus aureus RNA polymerase rpoB gene and comparison of its predicted amino acid sequence with those of other bacteria. Biochem. 1995, 1262: 73-78.Google Scholar
- Hou T, Wang J, Li Y, Wang W: Assessing the performance of the molecular mechanics/Poisson Boltzmann surface area and molecular mechanics/generalized Born surface area methods II: The accuracy of ranking poses generated from docking. J Comput Chem. 2011, 32 (5): 866-877. 10.1002/jcc.21666.PubMed CentralView ArticlePubMedGoogle Scholar
- Okimoto N, Futatsugi N, Fuji H, Suenaga A, Morimoto G, Yanai R, Ohno Y, Narumi T, Taiji M: High-performance drug discovery: computational screening by combining docking and molecular dynamics simulations. PLoS Comput Biol. 2009, 5 (10): e1000528-10.1371/journal.pcbi.1000528.PubMed CentralView ArticlePubMedGoogle Scholar
- Irwin JJ, Shoichet BK: ZINC: A free database of commercially available compounds for virtual screening. J Chem Inf Model. 2005, 45 (1): 177-182. 10.1021/ci049714+.PubMed CentralView ArticlePubMedGoogle Scholar
- Trott O, Olson AJ: AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. 2010, 31 (2): 455-461.PubMed CentralPubMedGoogle Scholar
- Aqvist J, Luzhkov VB, Brandsdal BO: Ligand binding affinities from MD simulations. Acc Chem Res. 2002, 35 (6): 358-365. 10.1021/ar010014p.View ArticlePubMedGoogle Scholar
- Hou T, Zhang W, Case DA, Wang W: Characterization of domain-peptide interaction interface: a case study on the amphiphysin-1 SH3 domain. J Mol Biol. 2008, 376 (4): 1201-1214. 10.1016/j.jmb.2007.12.054.View ArticlePubMedGoogle Scholar
- DeLano WL: The PyMOL Molecular Graphics System. 2002, San Carlos, CA: DeLano ScientificGoogle Scholar
- Sanner MF: Python: A Programming Language for Software Integration and Development. J Mol Graphics Mod. 1999, 17: 57-61.Google Scholar
- Case DA, Darden TA, Cheatham TE, Simmerling CL, Wang J, Duke RE, Luo R, Walker RC, Zhang W, Merz KM, Roberts B, Wang B, Hayik S, Roitberg A, Seabra G, Kolossvai I, Wong KF, Paesani F, Vanicek J, Liu J, Wu X, Brozell SR, Steinbrecher T, Gohlke H, Cai Q, Ye X, Wang J, Hsieh M-J, Cui G, Roe DR, Mathews DH, Seetin MG, Sagui C, Babin V, Luchko T, Gusarov S, Kovalenko A, Kollman PA: AMBER 12. 2012, San Francisco: University of CaliforniaGoogle Scholar
- Fogolari F, Brigo A, Molinari H: Protocol for MM/PBSA molecular dynamics simulations of proteins. Biophys J. 2003, 85 (1): 159-166. 10.1016/S0006-3495(03)74462-2.PubMed CentralView ArticlePubMedGoogle Scholar
- Wolber G, Langer T: LigandScout: 3-D Pharmacophores Derived from Protein-Bound Ligands and Their Use as Virtual Screening Filters. J Chem Inf Model. 2005, 45 (1): 160-169. 10.1021/ci049885e.View ArticlePubMedGoogle Scholar
- Zaheer-ul-Haq , Halim SA, Uddin R, Madura JD: Benchmarking docking and scoring protocol for the identification of potential acetylcholinesterase inhibitors. J Mol Graph Model. 2010, 28 (8)): 870-882.View ArticlePubMedGoogle Scholar
- Holtje HD, Sippl W, Rognan D, Folkers G: Molecular modeling: basic principles and applications. 2008, Weinheim, Cambridge: Wiley-VCH, 3Google Scholar
- Campbell EA, Pavlova O, Zenkin N, Leon F, Irschik H, Jansen R, Severinov K, Darst SA: Structural, functional, and genetic analysis of sorangicin inhibition of bacterial RNA polymerase. EMBO J. 2005, 24 (4): 674-682. 10.1038/sj.emboj.7600499.PubMed CentralView ArticlePubMedGoogle Scholar
- Korzheva N, Mustaev A, Kozlov M, Malhotra A, Nikiforov V, Goldfarb A, Darst SA: A structural model of transcription elongation. Science. 2000, 289: 619-625. 10.1126/science.289.5479.619.View ArticlePubMedGoogle Scholar
- Lee JC, Tseng CK, Chen KJ, Huang KJ, Lin CK, Lin YT: A cell-based reporter assay for inhibitor screening of hepatitis C virus RNA-dependent RNA polymerase. Anal Biochem. 2010, 403 (1–2): 52-62.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.