Ionizable Residues Research Articles

Mechanism of algal aggregation by Bacillus sp. strain RP1137.

Alga-derived biofuels are one of the best alternatives for economically replacing liquid fossil fuels with a fungible renewable energy source. Production of fuel from algae is technically feasible but not yet economically viable. Harvest of dilute algal biomass from the surrounding water remains one of the largest barriers to economic production of algal biofuel. We identified Bacillus sp. strain RP1137 in a previous study and showed that this strain can rapidly aggregate several biofuel-producing algae in a pH- and divalent-cation-dependent manner. In this study, we further characterized the mechanism of algal aggregation by RP1137. We show that aggregation of both algae and bacteria is optimal in the exponential phase of growth and that the density of ionizable residues on the RP1137 cell surface changes with growth stage. Aggregation likely occurs via charge neutralization with calcium ions at the cell surface of both algae and bacteria. We show that charge neutralization occurs at least in part through binding of calcium to negatively charged teichoic acid residues. The addition of calcium also renders both algae and bacteria more able to bind to hydrophobic beads, suggesting that aggregation may occur through hydrophobic interactions. Knowledge of the aggregation mechanism may enable engineering of RP1137 to obtain more efficient algal harvesting.

Study of pK Values and Effective Dielectric Constants of Ionizable Residues in Pentapeptides and in Staphylococcal Nuclease (SNase) Using a Mean-Field Approach

The determination of pK values of amino acid residues as a function of temperature and ionic concentration is crucial to understanding the dynamics of various biological processes such as adsorption of peptides and their interactions with active sites of enzymes. In this study we developed a mean-field model to calculate the position-dependent dielectric constants of ionizable groups and the mean electrostatic potential on the surface. Such potential, which takes into account the contributions exerted by neighboring groups and ions in solution, is responsible for the fine-tuning of the pK value of each residue. The proposed model was applied to the amino acids Asp, Glu, Lys, His, Tyr, and Cys, and since the results were consistent with experimentally obtained values, the model was extended and applied to computation of pK values of Gly and Ala pentapeptides and of ionizable residues of the enzyme staphylococcal nuclease (SNase). In this latter case, we used an approach similar to a first-neighbors approximation, and the results turned out to be in good agreement with previously reported data when considering only the interactions of charged groups located at distances of maximally 20 Å. These considerations and the little computational cost involved turn the suggested approach into a promising tool for the modeling of force fields in computational simulations.

Energy cost for the proper ionization of active site residues in 6-phosphogluconate dehydrogenase from T. brucei

The catalytic mechanism of 6-phosphogluconate dehydrogenase requires the inversion of a Lys/Glu couple from its natural ionization state. The pKa of these residues in free and substrate bound enzymes has been determined measuring by ITC the proton release/uptake induced by substrate binding at different pH values. Wt 6-phosphogluconate dehydrogenase from Trypanosoma brucei and two active site enzyme mutants, K185H and E192Q were investigated. Substrate binding was accompanied by proton release and was dependent on the ionization of a group with pKa 7.07 which was absent in the E192Q mutant. Kinetic data highlighted two pKa, 7.17 and 9.64, in the enzyme–substrate complex, the latter being absent in the E192Q mutant, suggesting that the substrate binding shifts Glu192 pKa from 7.07 to 9.64. A comparison of wt and E192Q mutant appears to show that the substrate binding shifts Lys185 pKa from 9.9 to 7.17. By comparing differences in proton release and the binding enthalpy of wt and mutant enzymes, the enthalpic cost of the change in the protonation state of Lys185 and Glu192 was estimated at ≈6.1kcal/mol. The change in protonation state of Lys185 and Glu192 has little effect on Gibbs free energy, 240–325cal/mol. However proton balance evidences the dissociation of other group(s) that can be collectively described by a single pKa shift from 9.1 to 7.54. This further change in ionization state of the enzyme causes an increase of free energy with a total cost of 1.2–2.3kcal/mol to set the enzyme into a catalytically competent form.

Constant pH molecular dynamics of proteins in explicit solvent with proton tautomerism

pH is a ubiquitous regulator of biological activity, including protein-folding, protein-protein interactions, and enzymatic activity. Existing constant pH molecular dynamics (CPHMD) models that were developed to address questions related to the pH-dependent properties of proteins are largely based on implicit solvent models. However, implicit solvent models are known to underestimate the desolvation energy of buried charged residues, increasing the error associated with predictions that involve internal ionizable residue that are important in processes like hydrogen transport and electron transfer. Furthermore, discrete water and ions cannot be modeled in implicit solvent, which are important in systems like membrane proteins and ion channels. We report on an explicit solvent constant pH molecular dynamics framework based on multi-site λ-dynamics (CPHMD(MSλD)). In the CPHMD(MSλD) framework, we performed seamless alchemical transitions between protonation and tautomeric states using multi-site λ-dynamics, and designed novel biasing potentials to ensure that the physical end-states are predominantly sampled. We show that explicit solvent CPHMD(MSλD) simulations model realistic pH-dependent properties of proteins such as the Hen-Egg White Lysozyme (HEWL), binding domain of 2-oxoglutarate dehydrogenase (BBL) and N-terminal domain of ribosomal protein L9 (NTL9), and the pKa predictions are in excellent agreement with experimental values, with a RMSE ranging from 0.72 to 0.84 pKa units. With the recent development of the explicit solvent CPHMD(MSλD) framework for nucleic acids, accurate modeling of pH-dependent properties of both major class of biomolecules-proteins and nucleic acids is now possible.

Steady-state kinetics with nitric oxide reductase (NOR): New considerations on substrate inhibition profile and catalytic mechanism

Nitric oxide reductase (NOR) from denitrifying bacteria is an integral membrane protein that catalyses the two electron reduction of NO to N2O, as part of the denitrification process, being responsible for an exclusive reaction, the NN bond formation, the key step of this metabolic pathway. Additionally, this class of enzymes also presents residual oxidoreductase activity, reducing O2 to H2O in a four electron/proton reaction. In this work we report, for the first time, steady-state kinetics with the Pseudomonas nautica NOR, either in the presence of its physiological electron donor (cyt. c552) or immobilised on a graphite electrode surface, in the presence of its known substrates, namely NO or O2. The obtained results show that the enzyme has high affinity for its natural substrate, NO, and different kinetic profiles according to the electron donor used. The kinetic data, as shown by the pH dependence, is modelled by ionisable amino acid residues nearby the di-nuclear catalytic site. The catalytic mechanism is revised and a mononitrosyl-non-heme Fe complex (FeBII-NO) species is favoured as the first catalytic intermediate involved on the NO reduction.

PH Dependence of Amylin Fibrillization

In type 2 diabetics, the hormone amylin misfolds into amyloid plaques implicated in the destruction of the pancreatic β-cells that make insulin and amylin. The aggregative misfolding of amylin is pH-dependent, and exposure of the hormone to acidic and basic environments could be physiologically important. Amylin has two ionizable residues between pH 3 and 9: the α-amino group and His18. Our approach to measuring the pKa values for these sites has been to look at the pH dependence of fibrillization in amylin variants that have only one of the two groups. The α-amino group at the unstructured N-terminus of amylin has a pKa near 8.0, similar to the value in random coil models. By contrast, His18, which is involved in the intermolecular β-sheet structure of the fibrils, has a pKa that is lowered to 5.0 in the fibrils compared to the random coil value of 6.5. The lowered pKa of His18 is due to the hydrophobic environment of the residue, and electrostatic repulsion between positively charged His18 residues on neighboring amylin molecules in the fibril. His18 acts as an electrostatic switch inhibiting fibrillization in its charged state. The presence of a charged side chain at position 18 also affects fibril morphology and lowers amylin cytotoxicity toward a MIN6 mouse model of pancreatic β-cells. In addition to the two expected pKa values, we detected an apparent pKa of ~4.0 for the amylin-derived peptide NAc-SNNFGAILSS-NH2, which has no titratable groups. This pKa is due to the pH-induced ionization of the dye thioflavin T. By using alternative methods to follow fibrillization such as the dye Nile Red or turbidimetry, we were able to distinguish between the titration of the dye and groups on the peptide. Large differences in reaction kinetics were observed between the different methods at acidic pH, because of charges on the ThT dye, which hinder fibril formation much like the charges on the protein.

Computational Evaluations of Charge Coupling and Hydrogen Bonding in the Active Site of a Family 7 Cellobiohydrolase

Solution pH and the pKa values of ionizable residues are critical factors known to influence enzyme catalysis, structural stability, and dynamical fluctuations. Presented here is an exhaustive computational study utilizing long time constant pH molecular dynamics, pH replica exchange simulations, and kinetic modeling to evaluate pH-dependent conformations, charge dynamics, residue pKa values, and the catalytic activity-pH profile for cellobiohydrolase Cel7B from Melanocarpus albomyces . The predicted pKa values support the role of Glu212 as the catalytic nucleophile and Glu217 as the acid-base residue. The presence of a charge-correlated active site and an extensive hydrogen bonding network is found to be critical in enabling favorable residue orientations for catalysis and shuttling excess protons around the active site. Clusters of amino acids are identified that act in concert to effectively modulate the optimal pH for catalysis while elevating the overall catalytic rate with respect to a noncoupled system. The work presented here demonstrates the complex and critical role of coupled ionizable residues to the proper functioning of cellobiohydrolase Cel7B, functionally related glycosyl hydrolases, and enzymes in general. The simulations also support the use of the CpHMD for the accurate prediction of residue pKa values and to evaluate the impact of pH on protein structure and charge dynamics.

Low-Temperature Trapping of Photointermediates of the Rhodopsin E181Q Mutant.

Three active-site components in rhodopsin play a key role in the stability and function of the protein: 1) the counter-ion residues which stabilize the protonated Schiff base, 2) water molecules, and 3) the hydrogen-bonding network. The ionizable residue Glu-181, which is involved in an extended hydrogen-bonding network with Ser-186, Tyr-268, Tyr-192, and key water molecules within the active site of rhodopsin, has been shown to be involved in a complex counter-ion switch mechanism with Glu-113 during the photobleaching sequence of the protein. Herein, we examine the photobleaching sequence of the E181Q rhodopsin mutant by using cryogenic UV-visible spectroscopy to further elucidate the role of Glu-181 during photoactivation of the protein. We find that lower temperatures are required to trap the early photostationary states of the E181Q mutant compared to native rhodopsin. Additionally, a Blue Shifted Intermediate (BSI, λmax = 498 nm, 100 K) is observed after the formation of E181Q Bathorhodopsin (Batho, λmax = 556 nm, 10 K) but prior to formation of E181Q Lumirhodopsin (Lumi, λmax = 506 nm, 220 K). A potential energy diagram of the observed photointermediates suggests the E181Q Batho intermediate has an enthalpy value 7.99 KJ/mol higher than E181Q BSI, whereas in rhodopsin, the BSI is 10.02 KJ/mol higher in enthalpy than Batho. Thus, the Batho to BSI transition is enthalpically driven in E181Q and entropically driven in native rhodopsin. We conclude that the substitution of Glu-181 with Gln-181 results in a significant perturbation of the hydrogen-bonding network within the active site of rhodopsin. In addition, the removal of a key electrostatic interaction between the chromophore and the protein destabilizes the protein in both the dark state and Batho intermediate conformations while having a stabilizing effect on the BSI conformation. The observed destabilization upon this substitution further supports that Glu-181 is negatively charged in the early intermediates of the photobleaching sequence of rhodopsin.

Redox specificity of 2-hydroxyacid-coupled NAD(+)/NADH dehydrogenases: a study exploiting "reactive" arginine as a reporter of protein electrostatics.

With “reactive” arginine as a kinetic reporter, 2-hydroxyacid dehydrogenases are assessed in basis of their specialization as NAD+-reducing or NADH-oxidizing enzymes. Specifically, M4 and H4 lactate dehydrogenases (LDHs) and cytoplasmic and mitochondrial malate dehydrogenases (MDHs) are compared to assess if their coenzyme specificity may involve electrostatics of cationic or neutral nicotinamide structure as the basis. The enzymes from diverse eukaryote and prokaryote sources thus are assessed in “reactivity” of functionally-critical arginine as a function of salt concentration and pH. Electrostatic calculations were performed on “reactive” arginines and found good correspondence with experiment. The reductive and oxidative LDHs and MDHs are assessed in their count over ionizable residues and in placement details of the residues in their structures as proteins. The variants found to be high or low in ΔpKa of “reactive” arginine are found to be also strong or weak cations that preferentially oxidize NADH (neutral nicotinamide structure) or reduce NAD+ (cationic nicotinamide structure). The ionized groups of protein structure may thus be important to redox specificity of the enzyme on basis of electrostatic preference for the oxidized (cationic nicotinamide) or reduced (neutral nicotinamide) coenzyme. Detailed comparisons of isozymes establish that the residues contributing in their redox specificity are scrambled in structure of the reductive enzyme.

High-throughput fluorescence assay for membrane-protein interaction

Membrane-protein interaction plays key roles in a wide variety of biological processes. Although various methods have been employed to measure membrane binding of soluble proteins, a robust high-throughput assay that is universally applicable to all proteins is lacking at present. Here we report a new fluorescence quenching assay utilizing enhanced green fluorescence protein (EGFP)-fusion proteins and a lipid containing a dark quencher, N-dimethylaminoazobenzenesulfonyl-phosphatidylethanolamine (dabsyl-PE). The EGFP fluorescence emission intensity showed a large decrease (i.e., >50%) when EGFP-fusion proteins bound the vesicles containing 5 mol% dabsyl-PE. This simple assay, which can be performed using either a cuvette-based spectrofluorometer or a fluorescence plate reader, allowed rapid, sensitive, and accurate determination of lipid specificity and affinity for various lipid binding domains, including two pleckstrin homology domains, an epsin N-terminal homology domain, and a phox homology domain. The assay can also be applied to high-throughput screening of small molecules that modulate membrane binding of proteins.

De Novo-Designed Metallopeptides with Type 2 Copper Centers: Modulation of Reduction Potentials and Nitrite Reductase Activities

Enzymatic reactions involving redox processes are highly sensitive to the local electrostatic environment. Despite considerable effort, the complex interactions among different influential factors in native proteins impede progress toward complete understanding of the structure-function relationship. Of particular interest is the type 2 copper center Cu(His)3, which may act as an electron transfer center in peptidylglycine α-hydroxylating monooxygenase (PHM) or a catalytic center in copper nitrite reductase (CuNiR). A de novo design strategy is used to probe the effect of modifying charged amino acid residues around, but not directly bound to, a Cu(His)3 center embedded in three-stranded coiled coils (TRI-H)3 [TRI-H = Ac-G WKALEEK LKALEEK LKALEEK HKALEEK G-NH2]. Specifically, the peptide TRI-EH (=TRI-HK22E) alters an important lysine to glutamate just above the copper binding center. With a series of TRI-EH peptides mutated below the metal center, we use a variety of spectroscopies (EPR, UV-vis, XAS) to show a direct impact on the protonation equilibria, copper binding affinities, reduction potentials, and nitrite reductase activities of these copper-peptide complexes. The potentials at a specific pH vary by 100 mV, and the nitrite reductase activities range over a factor of 4 in rates. We also observe that the affinities, potentials, and catalytic activities are strongly influenced by the pH conditions (pH 5.8-7.4). In general, Cu(II) affinities for the peptides are diminished at low pH values. The interplay among these factors can lead to a 200 mV shift in reduction potential across these peptides, which is determined by the pH-dependent affinities of copper in both oxidation states. This study illustrates the strength of de novo protein design in elucidating the influence of ionizable residues on a particular redox system, an important step toward understanding the factors that govern the properties of this metalloenzyme with a goal of eventually improving the catalytic activity.

Mutations of Asp540 and the domain-connecting residues synergistically enhance Pyrococcus furiosus DNA ligase activity

The structure of Pyrococcus furiosus DNA ligase (PfuLig), which architecturally resembles human DNA ligase I (hLigI), revealed that the C-terminal helix stabilizes the closed conformation through several ionic interactions between two domains (adenylylation domain (AdD) and C-terminal OB-fold domain (OBD)). This helix is oriented differently in DNA-bound hLigI, suggesting that the disruption of its interactions with AdD facilitates DNA binding. Previously, we demonstrated that the replacement of Asp540 with arginine improves the ligation activity. Here we report that the combination of the Asp540-replacement and the elimination of ionic residues in the helix, forming interactions with AdD, effectively enhanced the activity.

Kinetic role of a histidine residue in the T1 copper site of the laccase from Rigidoporus lignosus

Laccases (benzendiol:oxygen oxidoreductases; EC 1.10.3.2) catalyze the oxidation of a broad range of substrates, such as polyphenols, dyes and pollutants, and thus these enzymes are widely applied in industrial, biotechnological and environmental fields. In order to improve their biotechnological applications, a deep knowledge of structural factors involved in controlling their activity, in various experimental conditions and on different substrates, is required. In the present study, a laccase from the mushroom Rigidoporus lignosus was kinetically characterized. In particular, the stability, the effects of pH, ionic strength and fluoride ion concentration on the kinetic parameters were investigated, using three di-hydroxy-benzene isomers (1,2-dihydroxy-benzene, 1,3-dihydroxy-benzene and 1,4-dihydroxy-benzene) as substrates. The catalytic constant values of the laccase showed a bell-shaped pH profile, with the same optimum pH and pKa values for all tested substrates. This behavior appears to be due to the presence of an ionizable residue in the enzyme active site. To identify this residue, the enzyme was derivatized with diethyl-pyrocarbonate to modify accessible histidine residues, which, according to structural data, are present in the active site of this enzyme. The kinetic behavior of the derivatized laccase was compared with that of the native enzyme and the derivatized residues were identified by mass spectrometry. Mass spectrometry and kinetic results suggest the main role of His-457 in the control of the catalytic activity of laccase from R. lignosus.

PH modulates the binding of early growth response protein 1 transcription factor to DNA

The transcription factor early growth response protein (EGR)1 orchestrates a plethora of signaling cascades involved in cellular homeostasis, and its downregulation has been implicated in the development of prostate cancer. Herein, using a battery of biophysical tools, we show that the binding of EGR1 to DNA is tightly regulated by solution pH. Importantly, the binding affinity undergoes an enhancement of more than an order of magnitude with an increase in pH from 5 to 8, implying that the deprotonation of an ionizable residue accounts for such behavior. This ionizable residue is identified as His382 by virtue of the fact that its replacement by nonionizable residues abolishes the pH dependence of the binding of EGR1 to DNA. Notably, His382 inserts into the major groove of DNA, and stabilizes the EGR1-DNA interaction via both hydrogen bonding and van der Waals contacts. Remarkably, His382 is mainly conserved across other members of the EGR family, implying that histidine protonation-deprotonation may serve as a molecular switch for modulating the protein-DNA interactions that are central to this family of transcription factors. Collectively, our findings reveal an unexpected but a key step in the molecular recognition of the EGR family of transcription factors, and suggest that they may act as sensors of pH within the intracellular environment.

Depth: a web server to compute depth, cavity sizes, detect potential small-molecule ligand-binding cavities and predict the pKa of ionizable residues in proteins

Residue depth accurately measures burial and parameterizes local protein environment. Depth is the distance of any atom/residue to the closest bulk water. We consider the non-bulk waters to occupy cavities, whose volumes are determined using a Voronoi procedure. Our estimation of cavity sizes is statistically superior to estimates made by CASTp and VOIDOO, and on par with McVol over a data set of 40 cavities. Our calculated cavity volumes correlated best with the experimentally determined destabilization of 34 mutants from five proteins. Some of the cavities identified are capable of binding small molecule ligands. In this study, we have enhanced our depth-based predictions of binding sites by including evolutionary information. We have demonstrated that on a database (LigASite) of ∼200 proteins, we perform on par with ConCavity and better than MetaPocket 2.0. Our predictions, while less sensitive, are more specific and precise. Finally, we use depth (and other features) to predict pKas of GLU, ASP, LYS and HIS residues. Our results produce an average error of just <1 pH unit over 60 predictions. Our simple empirical method is statistically on par with two and superior to three other methods while inferior to only one. The DEPTH server (http://mspc.bii.a-star.edu.sg/depth/) is an ideal tool for rapid yet accurate structural analyses of protein structures.

A descoberta dos raios cósmicos ou o problema da ionização do ar atmosférico

A radiação cósmica foi descoberta por volta de 1911 por especialistas em condutividade no ar. Nenhum fenômeno teórico ou astronômico, ou ainda um fenômeno "mais próximo'', tinha sido previsto relativo à existência dos raios cósmicos. A sua descoberta resultou do estudo da causa geradora de uma grande quantidade de íons presentes na atmosfera. Neste artigo nós revisamos o contexto científico na época dessa descoberta e colocamos em perspectiva a participação dos diferentes atores. Nós esboçamos um inventário do problema da chamada ionização residual do ar tal como ele existia por volta de 1900.

Analysis of SARS-CoV E protein ion channel activity by tuning the protein and lipid charge

A partial characterization of the ion channels formed by the SARS coronavirus (CoV) envelope (E) protein was previously reported (C. Verdiá-Báguena et al., 2012 [12]). Here, we provide new significant insights on the involvement of lipids in the structure and function of the CoV E protein channel on the basis of three series of experiments. First, reversal potential measurements over a wide range of pH allow the dissection of the contributions to channel selectivity coming from ionizable residues of the protein transmembrane domain and also from the negatively charged groups of diphytanoyl phosphatidylserine (DPhPS) lipid. The corresponding effective pKas are consistent with the model pKas of the acidic residue candidates for titration. Second, the change of channel conductance with salt concentration reveals two distinct regimes (Donnan-controlled electrodiffusion and bulk-like electrodiffusion) fully compatible with the outcomes of selectivity experiments. Third, by measuring channel conductance in mixtures of neutral diphytanoyl phosphatidylcholine (DPhPC) lipids and negatively charged DPhPS lipids in low and high salt concentrations we conclude that the protein–lipid conformation in the channel is likely the same in charged and neutral lipids. Overall, the whole set of experiments supports the proteolipidic structure of SARS-CoV E channels and explains the large difference in channel conductance observed between neutral and charged membranes.

Effect of pH and ligand charge state on BACE-1 fragment docking performance

In this work we propose a protocol for estimating the effect of pH on the docking performance to BACE-1, which affords the charge state of the inhibitor as well as the protonation state of all ionisable residues in the protein at a given pH value. To the best of our knowledge, this is the first report of a protocol predicting the BACE-1 ligand docking poses not only at the neutral pH at which most crystallographic structures were obtained, but also at the optimal pH of the enzyme (in the acidic range), at which most of the BACE-1 binding affinity assays are performed. We have applied this protocol to a set of 23 fragment-like BACE-1 ligands that span four orders of magnitude in their binding affinities. The pK a values of the BACE-1 acidic residues deviate substantially from the estimates for model compounds in solution and display a ligand dependent variability, especially in the case of the catalytic Asp dyad residues. This outcome should have a strong bearing on the design of protocols for docking based BACE-1 screening campaigns. Finally, we were able to find an explanation for the poor docking success rate of some fragments based on the availability of anchoring points, a rationale that could help to improve hit rates in BACE-1 screening campaigns.

Disorder and order in unfolded and disordered peptides and proteins: A view derived from tripeptide conformational analysis. II. Tripeptides with short side chains populating asx and β‐type like turn conformations

In the preceding paper, we found that ensembles of tripeptides with long or bulky chains can include up to 20% of various turns. Here, we determine the structural and thermodynamic characteristics of GxG peptides with short polar and/or ionizable central residues (D, N, C), whose conformational distributions exhibit higher than average percentage (>20%) of turn conformations. To probe the side-chain conformations of these peptides, we determined the (3)J(H(α),H(β)) coupling constants and derived the population of three rotamers with χ1 -angles of -60°, 180° and 60°, which were correlated with residue propensities by DFT-calculations. For protonated GDG, the rotamer distribution provides additional evidence for asx-turns. A comparison of vibrational spectra and NMR coupling constants of protonated GDG, ionized GDG, and the protonated aspartic acid dipeptide revealed that side chain protonation increases the pPII content at the expense of turn populations. The charged terminal groups, however, have negligible influence on the conformational properties of the central residue. Like protonated GDG, cationic GCG samples asx-turns to a significant extent. The temperature dependence of the UVCD spectra and (3)J(H(N)H(α)) constants suggest that the turn populations of GDG and GNG are practically temperature-independent, indicating enthalpic and entropic stabilization. The temperature-independent J-coupling and UVCD spectra of GNG require a three-state model. Our results indicate that short side chains with hydrogen bonding capability in GxG segments of proteins may serve as hinge regions for establishing compact structures of unfolded proteins and peptides.

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.Electronic supplementary materialThe online version of this article (doi:10.1007/s10822-013-9643-9) contains supplementary material, which is available to authorized users.