Molecular modeling and dynamics simulations of PNP from Streptococcus agalactiae

Docking with GemDock

GEMDOCK is a protein-ligand docking software that makes use of an elegant biologically inspired computational methodology based on the differential evolution algorithm. As any docking program, GEMDOCK has two major features to predict the binding of a small-molecule ligand to the binding site of a protein target: the search algorithm and the scoring function to evaluate the generated poses. The GEMDOCK scoring function uses a piecewise potential energy function integrated into the differential evolutionary algorithm. GEMDOCK has been applied to a wide range of protein systems with docking accuracy similar to other docking programs such as Molegro Virtual Docker, AutoDock4, and AutoDock Vina. In this chapter, we explain how to carry out protein-ligand docking simulations with GEMDOCK. We focus this tutorial on the protein target cyclin-dependent kinase 2.

SAnDReS: A Computational Tool for Docking

Since the early 1980s, we have witnessed considerable progress in the development and application of docking programs to assess protein-ligand interactions. Most of these applications had as a goal the identification of potential new binders to protein targets. Another remarkable progress is taking place in the determination of the structures of protein-ligand complexes, mostly using X-ray diffraction crystallography. Considering these developments, we have a favorable scenario for the creation of a computational tool that integrates into one workflow all steps involved in molecular docking simulations. We had these goals in mind when we developed the program SAnDReS. This program allows the integration of all computational features related to modern docking studies into one workflow. SAnDReS not only carries out docking simulations but also evaluates several docking protocols allowing the selection of the best approach for a given protein system. SAnDReS is a free and open-source (GNU General Public License) computational environment for running docking simulations. Here, we describe the combination of SAnDReS and AutoDock4 for protein-ligand docking simulations. AutoDock4 is a free program that has been applied to over a thousand receptor-ligand docking simulations. The dataset described in this chapter is available for downloading at https://github.com/azevedolab/sandres.

Van der Waals Potential in Protein Complexes

Van der Waals forces are determinants of the formation of protein-ligand complexes. Physical models based on the Lennard-Jones potential can estimate van der Waals interactions with considerable accuracy and with a computational complexity that allows its application to molecular docking simulations and virtual screening of large databases of small organic molecules. Several empirical scoring functions used to evaluate protein-ligand interactions approximate van der Waals interactions with the Lennard-Jones potential. In this chapter, we present the main concepts necessary to understand van der Waals interactions relevant to molecular recognition of a ligand by the binding pocket of a protein target. We describe the Lennard-Jones potential and its application to calculate potential energy for an ensemble of structures to highlight the main features related to the importance of this interaction for binding affinity.

Homology Modeling of Protein Targets with MODELLER

Homology modeling is a computational approach to generate three-dimensional structures of protein targets when experimental data about similar proteins are available. Although experimental methods such as X-ray crystallography and nuclear magnetic resonance spectroscopy successfully solved the structures of nearly 150,000 macromolecules, there is still a gap in our structural knowledge. We can fulfill this gap with computational methodologies. Our goal in this chapter is to explain how to perform homology modeling of protein targets for drug development. We choose as a homology modeling tool the program MODELLER. To illustrate its use, we describe how to model the structure of human cyclin-dependent kinase 3 using MODELLER. We explain the modeling procedure of CDK3 apoenzyme and the structure of this enzyme in complex with roscovitine.

Hydrogen Bonds in Protein-Ligand Complexes

Fast and reliable evaluation of the hydrogen bond potential energy has a significant impact in the drug design and development since it allows the assessment of large databases of organic molecules in virtual screening projects focused on a protein of interest. Semi-empirical force fields implemented in molecular docking programs make it possible the evaluation of protein-ligand binding affinity where the hydrogen bond potential is a common term used in the calculation. In this chapter, we describe the concepts behind the programs used to predict hydrogen bond potential energy employing semi-empirical force fields as the ones available in the programs AMBER, AutoDock4, TreeDock, and ReplicOpter. We described here the 12-10 potential and applied it to evaluate the binding affinity for an ensemble of crystallographic structures for which experimental data about binding affinity are available.

Molecular Dynamics Simulation of Cholera Toxin A-1 Polypeptide

A molecular dynamics (MD) simulation study of the enzymatic portion of cholera toxin; cholera toxin A-1 polypeptide (CTA1) was performed at 283, 310 and 323 K. From total energy analysis it was observed that this toxin is stable thermodynamically and these outcomes were likewise confirmed by root mean square deviations (RMSD) investigations. The Cα root mean square fluctuation (RMSF) examinations revealed that there are a number of residues inside CTA1, which can be used as target for designing and synthesizing inhibitory drugs, in order to inactivate cholera toxin inside the human body. The fluctuations in the radius of gyration and hydrogen bonding in CTA1 proved that protein unfolding and refolding were normal routine phenomena in its structure at all temperatures. Solvent accessible surface area study identified the hydrophilic nature of the CTA1, and due to this property it can be a potential biological weapon. The structural identification (STRIDE) algorithm for proteins was successfully used to determine the partially disordered secondary structure of CTA1. On account of this partially disordered secondary structure, it can easily deceive the proteolytic enzymes of the endoplasmic reticulum of host cells.

Virtual screening: an in silico tool for interlacing the chemical universe with the proteome

In silico screening both in the forward (traditional virtual screening) and reverse sense (inverse virtual screening (IVS)) are helpful techniques for interlacing the chemical universe of small molecules with the proteome. The former, which is using a protein structure and a large chemical database, is well-known by the scientific community. We have chosen here to provide an overview on the latter, focusing on validation and target prioritization strategies. By comparing it to complementary or alternative wet-lab approaches, we put IVS in the broader context of chemical genomics, target discovery and drug design. By giving examples from the literature and an own example on how to validate the approach, we provide guidance on the issues related to IVS.

Crystal structure and molecular dynamics studies of purine nucleoside phosphorylase from Mycobacterium tuberculosis associated with acyclovir

Consumption has been a scourge of mankind since ancient times. This illness has charged a high price to human lives. Many efforts have been made to defeat Mycobacterium tuberculosis (Mt). The M. tuberculosis purine nucleoside phosphorylase (MtPNP) is considered an interesting target to pursuit new potential inhibitors, inasmuch it belongs to the purine salvage pathway and its activity might be involved in the mycobacterial latency process. Here we present the MtPNP crystallographic structure associated with acyclovir and phosphate (MtPNP:ACY:PO(4)) at 2.10 Å resolution. Molecular dynamics simulations were carried out in order to dissect MtPNP:ACY:PO(4) structural features, and the influence of the ligand in the binding pocket stability. Our results revealed that the ligand leads to active site lost of stability, in agreement with experimental results, which demonstrate a considerable inhibitory activity against MtPNP (K(i) = 150 nM). Furthermore, we observed that some residues which are important in the proper ligand's anchor into the human homologous enzyme do not present the same importance to MtPNP. Therewithal, these findings contribute to the search of new specific inhibitors for MtPNP, since peculiarities between the mycobacterial and human enzyme binding sites have been identified, making a structural-based drug design feasible.

Exploring the molecular basis for selective binding of Mycobacterium tuberculosis Asp kinase toward its natural substrates and feedback inhibitors: a docking and molecular dynamics study

Tuberculosis (TB) is still a major public health problem, compounded by the human immunodeficiency virus (HIV)-TB co-infection and recent emergence of multidrug-resistant (MDR) and extensively drug resistant (XDR)-TB. In this context, aspartokinase of mycobacterium tuberculosis has drawn attention for designing novel anti-TB drugs. Asp kinase is an enzyme responsible for the synthesis of 4-phospho-L-aspartate from L-aspartate and involved in the branched biosynthetic pathway leading to the synthesis of amino acids lysine, threonine, methionine and isoleucine. An intermediate of lysine biosynthetic branch, mesodiaminopimelate is also a component of the peptidoglycan which is a component of bacterial cell wall. To interfere with the production of all these amino acids and cell wall, it is possible to inhibit Asp kinase activity. This can be achieved using Asp kinase inhibitors. In order to design novel Asp kinase inhibitors as effective anti-TB drugs, it is necessary to have an understanding of the binding sites of Asp kinase. As no crystal structure of the enzyme has yet been published, we built a homology model of Asp kinase using the crystallized Asp kinase from M. Jannaschii, as template structures (2HMF and 3C1M). After the molecular dynamics refinement, the optimized homology model was assessed as a reliable structure by PROCHECK, ERRAT, WHAT-IF, PROSA2003 and VERIFY-3D. The results of molecular docking studies with natural substrates, products and feedback inhibitors are in agreement with the published data and showed that ACT domain plays an important role in binding to ligands. Based on the docking conformations, pharmacophore model can be developed by probing the common features of ligands. By analyzing the results, ACT domain architecture, certain key residues that are responsible for binding to feedback inhibitors and natural substrates were identified. This would be very helpful in understanding the blockade mechanism of Asp kinase and providing insights into rational design of novel Asp kinase inhibitors for M.tuberculosis.

Crystal structure and molecular dynamics studies of human purine nucleoside phosphorylase complexed with 7-deazaguanine

In humans, purine nucleoside phosphorylase (HsPNP) is responsible for degradation of deoxyguanosine, and genetic deficiency of this enzyme leads to profound T-cell mediated immunosuppression. HsPNP is a target for inhibitor development aiming at T-cell immune response modulation. Here we report the crystal structure of HsPNP in complex with 7-deazaguanine (HsPNP:7DG) at 2.75 A. Molecular dynamics simulations were employed to assess the structural features of HsPNP in both free form and in complex with 7DG. Our results show that some regions, responsible for entrance and exit of substrate, present a conformational variability, which is dissected by dynamics simulation analysis. Enzymatic assays were also carried out and revealed that 7-deazaguanine presents a lower inhibitory activity against HsPNP (K(i)=200 microM). The present structure may be employed in both structure-based design of PNP inhibitors and in development of specific empirical scoring functions.

Molecular dynamics studies of a hexameric purine nucleoside phosphorylase

Purine nucleoside phosphorylase (PNP) (EC.2.4.2.1) is an enzyme that catalyzes the cleavage of N-ribosidic bonds of the purine ribonucleosides and 2-deoxyribonucleosides in the presence of inorganic orthophosphate as a second substrate. This enzyme is involved in purine-salvage pathway and has been proposed as a promising target for design and development of antimalarial and antibacterial drugs. Recent elucidation of the three-dimensional structure of PNP by X-ray protein crystallography left open the possibility of structure-based virtual screening initiatives in combination with molecular dynamics simulations focused on identification of potential new antimalarial drugs. Most of the previously published molecular dynamics simulations of PNP were carried out on human PNP, a trimeric PNP. The present article describes for the first time molecular dynamics simulations of hexameric PNP from Plasmodium falciparum (PfPNP). Two systems were simulated in the present work, PfPNP in ligand free form, and in complex with immucillin and sulfate. Based on the dynamical behavior of both systems the main results related to structural stability and protein-drug interactions are discussed.

Molecular modeling, dynamics and docking studies of purine nucleoside phosphorylase from Streptococcus pyogenes

Purine Nucleoside Phosphorylase (PNP) catalyzes the reversible phosphorolysis of N-glycosidic bonds of purine nucleosides and deoxynucleosides, except for adenosine, to generate ribose 1-phosphate and the purine base. PNP has been submitted to intensive structural studies. This work describes for the first time a structural model of PNP from Streptococcus pyogenes (SpPNP). We modeled the complexes of SpPNP with six different ligands in order to determine the structural basis for specificity of these ligands against SpPNP. Molecular dynamics (MD) simulations were performed in order to evaluate the overall stability of SpPNP model. The analysis of the MD simulation was assessed mainly by principal component analysis (PCA) to explore the trimeric structure behavior. Structural comparison, between SpPNP and human PNP, was able to identify the main features responsible for differences in ligand-binding affinities, such as mutation in the purine-binding site and in the second phosphate-binding site. The PCA analysis suggests a different behavior for each subunit in the trimer structure.

Molecular modeling and dynamics simulation of human cyclin-dependent kinase 3 complexed with inhibitors

The complex CDK3-cyclin is involved in the control of the progression of G0. While the mechanisms governing early and late G1 progression are well understood, very little is known about the G0-G1 transition. Human CDK3 is closely related to cyclin-dependent kinase 2 (CDK2). Since there is no crystallographic structure of human CDK3, this work describes for the first time a molecular model of human CDK3 complexed with several inhibitors. Comparison of the binary complexes with different inhibitors strongly indicates that those inhibitors should inhibit CDK3 as well as CDK2.

Structural studies of shikimate dehydrogenase from Bacillus anthracis complexed with cofactor NADP

Bacillus anthracis has been employed as an agent of bioterrorism, with high mortality, despite anti-microbial treatment, which strongly indicates the need of new drugs to treat anthrax. Shikimate pathway is a seven step biosynthetic route which generates chorismic acid from phosphoenol pyruvate and erythrose-4-phosphate. Chorismic acid is the major branch point in the synthesis of aromatic amino acids, ubiquinone, and secondary metabolites. The shikimate pathway is essential for many pathological organisms, whereas it is absent in mammals. Therefore, these enzymes are potential targets for the development of nontoxic antimicrobial agents and herbicides and have been submitted to intensive structural studies. The forth enzyme of this pathway is responsible for the conversion of dehydroshikimate to shikimate in the presence of NADP. In order to pave the way for structural and functional efforts toward development of new antimicrobials we describe the molecular modeling of shikimate dehydrogenase from Bacillus anthracis complexed with the cofactor NADP. This study was able to identify the main residues of the NADP binding site responsible for ligand affinities. This structural study can be used in the design of more specific drugs against infectious diseases.