Phosphate closes the solution structure of the 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) from Mycobacterium tuberculosis

Discovery of new inhibitors of Mycobacterium tuberculosis EPSP synthase - A computational study

Tuberculosis (TB) is a highly infectious disease caused by the pathogen Mycobacterium tuberculosis (Mtb). EPSP Synthase (MtEPSPS), the enzyme responsible for the sixth step of the shikimate pathway, is a potential target for the development of new drugs for the treatment of TB, as it is essential in mycobacteria but absent in humans. In this work, we performed virtual screening using sets of molecules from two databases and three crystallographic structures of MtEPSPS. The initial hits obtained from molecular docking were filtered based on predicted binding affinity and interactions with binding site residues. Subsequently, molecular dynamics simulations were carried out to analyze the stability of protein-ligand complexes. We have found that MtEPSPS forms stable interactions with several candidates, including already approved pharmaceutical drugs such as Conivaptan and Ribavirin monophosphate. In particular, Conivaptan had the highest estimated binding affinity with the open conformation of the enzyme. The complex formed between MtEPSPS and Ribavirin monophosphate was also energetically stable as shown by RMSD, Rg and FEL analyses, and the ligand was stabilized by hydrogen bonds with important residues of the binding site. The findings reported in this work could serve as the basis of promising scaffolds for the discovery, design, and development of new anti-TB drugs.

Using modern approaches to sedimentation velocity to detect conformational changes in proteins

It has been known for decades that proteins undergo conformational changes in response to binding ligands. Such changes are usually accompanied by a loss of entropy by the protein, and thus conformational changes are integral to the thermodynamics of ligand association. Methods to detect these alterations are numerous; here, we focus on the sedimentation velocity (SV) mode of AUC, which has several advantages, including ease of use and rigorous data-selection criteria. In SV, it is assumed that conformational changes manifest primarily as differences in the sedimentation coefficient (the s-value). Two methods of determining s-value differences were assessed. The first method used the widely adopted c(s) distribution to gather statistics on the s-value differences to determine whether the observed changes were reliable. In the second method, a decades-old technique called "difference SV" was revived and updated to address its viability in this era of modern instrumentation. Both methods worked well to determine the extent of conformational changes to three model systems. Both simulations and experiments were used to explore the strengths and limitations of the methods. Finally, software incorporating these methodologies was produced.

Mycobacterium tuberculosis Shikimate Pathway Enzymes as Targets for the Rational Design of Anti-Tuberculosis Drugs

Roughly a third of the world’s population is estimated to have latent Mycobacterium tuberculosis infection, being at risk of developing active tuberculosis (TB) during their lifetime. Given the inefficacy of prophylactic measures and the increase of drug-resistant M. tuberculosis strains, there is a clear and urgent need for the development of new and more efficient chemotherapeutic agents, with selective toxicity, to be implemented on patient treatment. The component enzymes of the shikimate pathway, which is essential in mycobacteria and absent in humans, stand as attractive and potential targets for the development of new drugs to treat TB. This review gives an update on published work on the enzymes of the shikimate pathway and some insight on what can be potentially explored towards selective drug development.

How Docking Programs Work

Protein-ligand docking simulations are of central interest for computer-aided drug design. Docking is also of pivotal importance to understand the structural basis for protein-ligand binding affinity. In the last decades, we have seen an explosion in the number of three-dimensional structures of protein-ligand complexes available at the Protein Data Bank. These structures gave further support for the development and validation of in silico approaches to address the binding of small molecules to proteins. As a result, we have now dozens of open source programs and web servers to carry out molecular docking simulations. The development of the docking programs and the success of such simulations called the attention of a broad spectrum of researchers not necessarily familiar with computer simulations. In this scenario, it is essential for those involved in experimental studies of protein-ligand interactions and biophysical techniques to have a glimpse of the basics of the protein-ligand docking simulations. Applications of protein-ligand docking simulations to drug development and discovery were able to identify hits, inhibitors, and even drugs. In the present chapter, we cover the fundamental ideas behind protein-ligand docking programs for non-specialists, which may benefit from such knowledge when studying molecular recognition mechanism.

Electrostatic Energy in Protein-Ligand Complexes

Computational analysis of protein-ligand interactions is of pivotal importance for drug design. Assessment of ligand binding energy allows us to have a glimpse of the potential of a small organic molecule as a ligand to the binding site of a protein target. Considering scoring functions available in docking programs such as AutoDock4, AutoDock Vina, and Molegro Virtual Docker, we could say that they all rely on equations that sum each type of protein-ligand interactions to model the binding affinity. Most of the scoring functions consider electrostatic interactions involving the protein and the ligand. In this chapter, we present the main physics concepts necessary to understand electrostatics interactions relevant to molecular recognition of a ligand by the binding pocket of a protein target. Moreover, we analyze the electrostatic potential energy for an ensemble of structures to highlight the main features related to the importance of this interaction for binding affinity.

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.

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.

Structural basis of natural promoter recognition by a unique nuclear receptor, HNF4alpha. Diabetes gene product

HNF4alpha (hepatocyte nuclear factor 4alpha) plays an essential role in the development and function of vertebrate organs, including hepatocytes and pancreatic beta-cells by regulating expression of multiple genes involved in organ development, nutrient transport, and diverse metabolic pathways. As such, HNF4alpha is a culprit gene product for a monogenic and dominantly inherited form of diabetes, known as maturity onset diabetes of the young (MODY). As a unique member of the nuclear receptor superfamily, HNF4alpha recognizes target genes containing two hexanucleotide direct repeat DNA-response elements separated by one base pair (DR1) by exclusively forming a cooperative homodimer. We describe here the 2.0 angstroms crystal structure of human HNF4alpha DNA binding domain in complex with a high affinity promoter element of another MODY gene, HNF1alpha, which reveals the molecular basis of unique target gene selection/recognition, DNA binding cooperativity, and dysfunction caused by diabetes-causing mutations. The predicted effects of MODY mutations have been tested by a set of biochemical and functional studies, which show that, in contrast to other MODY gene products, the subtle disruption of HNF4alpha molecular function can cause significant effects in afflicted MODY patients.

Structural studies of shikimate 5-dehydrogenase from Mycobacterium tuberculosis

Tuberculosis (TB) remains the leading cause of mortality due to a single bacterial pathogen, Mycobacterium tuberculosis. The reemergence of TB as a potential public health threat, the high susceptibility of human immunodeficiency virus-infected persons to the disease, the proliferation of multi-drug-resistant strains (MDR-TB) and, more recently, of extensively drug resistant isolates (XDR-TB) have created a need for the development of new antimycobacterial agents. Amongst the several proteins and/or enzymes to be studied as potential targets to develop novel drugs against M. tuberculosis, the enzymes of the shikimate pathway are attractive targets because they are essential in algae, higher plants, bacteria, and fungi, but absent from mammals. The mycobacterial shikimate pathway leads to the biosynthesis of chorismate, which is a precursor of aromatic amino acids, naphthoquinones, menaquinones, and mycobactins. Here we report the structural studies by homology modeling and circular dichroism spectroscopy of the shikimate dehydrogenase from M. tuberculosis (MtSDH), which catalyses the fourth step of the shikimate pathway. Our structural models show that the MtSDH has similar structure to other shikimate dehydrogenase structures previously reported either in presence or absence of NADP, despite the low amino acid sequence identity. The circular dichroism spectra corroborate the secondary structure content observed in the MtSDH models developed. The enzyme was stable up to 50 degrees C presenting a cooperative unfolding profile with the midpoint of the unfolding temperature value of approximately 63-64 degrees C, as observed in the unfolding experiment followed by circular dichroism. Our MtSDH structural models and circular dichroism data showed small conformational changes induced by NADP binding. We hope that the data presented here will assist the rational design of antitubercular agents.

Molecular modeling and dynamics studies of Shikimate Kinase from Bacillus anthracis

Bacillus anthracis has been used as weapon in bioterrorist activities, with high mortality, despite anti-microbial treatment, which strongly indicates a need of new drugs to treat anthrax. Shikimate Pathway is a seven-step biosynthetic route which generates chorismic acid. 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 non-toxic anti-microbial agents and herbicides and have been submitted to intensive structural studies. Shikimate Kinase is the fifth enzyme of shikimate pathway and catalyzes the specific phosphorylation of the 3-hydroxyl group of shikimate using ATP as a co-substrate, resulting in shikimate-3-phosphate and ADP. The present work describes for the first time a structural model for the Shikimate Kinase from B. anthracis using molecular modeling approach and molecular dynamics simulations. This study was able to identify the main residues of the ATP-binding and the shikimate pockets responsible for ligand affinities. Analysis of the molecular dynamics simulations indicates the structural features responsible for the stability of the structure. This study may help in the identification of new inhibitors for this enzyme.