Effects of histidine protonation and rotameric states on virtual screening of M. tuberculosis RmlC

While it is well established that protonation and tautomeric states of ligands can significantly affect the results of virtual screening, such effects of ionizable residues of protein receptors are less well understood. In this study, we focus on histidine protonation and rotameric states and their impact on virtual screening of Mycobacterium tuberculosis enzyme RmlC. Depending on the net charge and the location of proton(s), a histidine can adopt three states: HIP (+1 charged, both δ- and ε-nitrogens protonated), HID (neutral, δ-nitrogen protonated), and HIE (neutral, ε-nitrogen protonated). Due to common ambiguities in X-ray crystal structures, a histidine may also be resolved as three additional states with its imidazole ring flipped. Here, we systematically investigate the predictive power of 36 receptor models with different protonation and rotameric states of two histidines in the RmlC active site by using results from a previous high-throughput screening. By measuring enrichment factors and area under the receiver operating characteristic curves, we show that virtual screening results vary depending on hydrogen bonding networks provided by the histidines, even in the cases where the ligand does not obviously interact with the side chain. Our results also suggest that, even with the help of widely used pKa prediction software, assigning histidine protonation and rotameric states for virtual screening can still be challenging and requires further examination and systematic characterization of the receptor-ligand complex.

[1]  Ramu Anandakrishnan,et al.  H++ 3.0: automating pK prediction and the preparation of biomolecular structures for atomistic molecular modeling and simulations , 2012, Nucleic Acids Res..

[2]  Cen Gao,et al.  Estimating binding affinities by docking/scoring methods using variable protonation states , 2011, Proteins.

[3]  Jan H. Jensen,et al.  Very fast prediction and rationalization of pKa values for protein–ligand complexes , 2008, Proteins.

[4]  Jonathan W. Essex,et al.  Ensemble Docking into Multiple Crystallographically Derived Protein Structures: An Evaluation Based on the Statistical Analysis of Enrichments , 2010, J. Chem. Inf. Model..

[5]  Junjun Mao,et al.  MCCE2: Improving protein pKa calculations with extensive side chain rotamer sampling , 2009, J. Comput. Chem..

[6]  D. Silverman,et al.  Proton transfer by histidine 67 in site-directed mutants of human carbonic anhydrase III. , 1995, Biochemistry.

[7]  Manfred Eigen,et al.  Proton Transfer, Acid-Base Catalysis, and Enzymatic Hydrolysis. Part I: ELEMENTARY PROCESSES†‡ , 1964 .

[8]  H. Carlson.  Protein flexibility and drug design: how to hit a moving target. , 2002, Current opinion in chemical biology.

[9]  S. Diamond,et al.  Identification of triazinoindol-benzimidazolones as nanomolar inhibitors of the Mycobacterium tuberculosis enzyme TDP-6-deoxy-d-xylo-4-hexopyranosid-4-ulose 3,5-epimerase (RmlC). , 2010, Bioorganic & medicinal chemistry.

[10]  Yvonne C. Martin,et al.  Let’s not forget tautomers , 2009, J. Comput. Aided Mol. Des..

[11]  B. X. Carlson,et al.  A single glycine residue at the entrance to the first membrane-spanning domain of the gamma-aminobutyric acid type A receptor beta(2) subunit affects allosteric sensitivity to GABA and anesthetics. , 2000, Molecular pharmacology.

[12]  Matthew P. Repasky,et al.  Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. , 2006, Journal of medicinal chemistry.

[13]  M. Hong,et al.  Protonation, tautomerization, and rotameric structure of histidine: a comprehensive study by magic-angle-spinning solid-state NMR. , 2011, Journal of the American Chemical Society.

[14]  C. Soares,et al.  Constant-pH molecular dynamics using stochastic titration , 2002 .

[15]  D. Case,et al.  Constant pH molecular dynamics in generalized Born implicit solvent , 2004, J. Comput. Chem..

[16]  Antti Poso,et al.  The Effect of Ligand-Based Tautomer and Protomer Prediction on Structure-Based Virtual Screening. , 2010 .

[17]  D. D. Perrin,et al.  pKa prediction for organic acids and bases , 1981 .

[18]  Jean-Paul Renaud,et al.  The role of the distal histidine in myoglobin and haemoglobin , 1988, Nature.

[19]  Lenwood S. Heath,et al.  H++: a server for estimating pKas and adding missing hydrogens to macromolecules , 2005, Nucleic Acids Res..

[20]  Hege S. Beard,et al.  Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. , 2004, Journal of medicinal chemistry.

[21]  Jan H. Jensen,et al.  PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. , 2011, Journal of chemical theory and computation.

[22]  Tom Fawcett,et al.  An introduction to ROC analysis , 2006, Pattern Recognit. Lett..

[23]  Jeremy R. Greenwood,et al.  Epik: a software program for pKa prediction and protonation state generation for drug-like molecules , 2007, J. Comput. Aided Mol. Des..

[24]  Rommie E. Amaro,et al.  An improved relaxed complex scheme for receptor flexibility in computer-aided drug design , 2008, J. Comput. Aided Mol. Des..

[25]  Andreas Bender,et al.  A Discussion of Measures of Enrichment in Virtual Screening: Comparing the Information Content of Descriptors with Increasing Levels of Sophistication , 2005, J. Chem. Inf. Model..

[26]  J. Errey,et al.  RmlC, a C3' and C5' carbohydrate epimerase, appears to operate via an intermediate with an unusual twist boat conformation. , 2007, Journal of molecular biology.

[27]  Anthony Nicholls,et al.  What do we know and when do we know it? , 2008, J. Comput. Aided Mol. Des..

[28]  E. Alexov,et al.  Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins. , 2002, Biophysical journal.

[29]  Alexander D. MacKerell,et al.  Consideration of Molecular Weight during Compound Selection in Virtual Target-Based Database Screening , 2003, J. Chem. Inf. Comput. Sci..

[30]  W. Cleland.  Low-barrier hydrogen bonds and enzymatic catalysis. , 2000, Archives of biochemistry and biophysics.

[31]  Patrick G. Blachly,et al.  Measuring the successes and deficiencies of constant pH molecular dynamics: A blind prediction study , 2011, Proteins.

[32]  Jan H. Jensen,et al.  Very fast empirical prediction and rationalization of protein pKa values , 2005, Proteins.

[33]  Woody Sherman,et al.  Analysis and comparison of 2D fingerprints: insights into database screening performance using eight fingerprint methods , 2010, J. Cheminformatics.

[34]  F. Cohen,et al.  Prion protein selectively binds copper(II) ions. , 1998, Biochemistry.

[35]  David G. Lloyd,et al.  Considerations in Compound Database Preparation-"Hidden" Impact on Virtual Screening Results , 2005, J. Chem. Inf. Model..

[36]  Thomas E. Exner,et al.  Influence of Protonation, Tautomeric, and Stereoisomeric States on Protein-Ligand Docking Results , 2009, J. Chem. Inf. Model..

[37]  Scott G. Franzblau,et al.  Drug Targeting Mycobacterium tuberculosis Cell Wall Synthesis: Genetics of dTDP-Rhamnose Synthetic Enzymes and Development of a Microtiter Plate-Based Screen for Inhibitors of Conversion of dTDP-Glucose to dTDP-Rhamnose , 2001, Antimicrobial Agents and Chemotherapy.

[38]  J. Mccammon,et al.  Accounting for Receptor Flexibility and Enhanced Sampling Methods in Computer‐Aided Drug Design , 2013, Chemical biology & drug design.

[39]  Matthew P. Repasky,et al.  Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. , 2004, Journal of medicinal chemistry.

[40]  J A McCammon,et al.  Accommodating protein flexibility in computational drug design. , 2000, Molecular pharmacology.

[41]  P. Paufler.  Crystal Structure Analysis for Chemists and Biologists. , 1995 .

[42]  Gregory A. Grothaus,et al.  A simple clustering algorithm can be accurate enough for use in calculations of pKs in macromolecules , 2006, Proteins.

[43]  D. Silverman,et al.  Role of histidine 64 in the catalytic mechanism of human carbonic anhydrase II studied with a site-specific mutant. , 1989, Biochemistry.

[44]  Jan H. Jensen,et al.  Improved Treatment of Ligands and Coupling Effects in Empirical Calculation and Rationalization of pKa Values. , 2011, Journal of chemical theory and computation.

[45]  Woody Sherman,et al.  Large-Scale Systematic Analysis of 2D Fingerprint Methods and Parameters to Improve Virtual Screening Enrichments , 2010, J. Chem. Inf. Model..

[46]  Xin Li,et al.  Assignment of polar states for protein amino acid residues using an interaction cluster decomposition algorithm and its application to high resolution protein structure modeling , 2006, Proteins.

[47]  E. Alexov,et al.  Incorporating protein conformational flexibility into the calculation of pH-dependent protein properties. , 1997, Biophysical journal.