Determination Of Protonation States Research Articles

Determination of protonation state in molecular salt of minoxidil and 2,4-dihydroxybenzoic acid through a combined experimental and theoretical study: influence of proton transfer on biological activities

Topical minoxidil (MXD) is frequently used as therapy for alopecia. In this study, a new MXD molecular salt formed by MXD and the 2,4-dihydroxybenzoic acid (24HA) were prepared and characterized. The nature of the MXD/24HA salt was studied using infrared spectrometry, single-crystal and powder X-ray diffraction, X-ray photoelectron spectroscopy, Hirshfeld surface analysis, and configuration interaction calculations. The theoretical results and the experimental data show that the proton transfer accompanies the formation of the MXD/24HA salt. The resulted MXD/24HA salt crystallizes in themonoclinic space group P21/c. The in vivo assessment with the hair growth model in mice reveals that the MXD/24HA salt is superior to the raw MXD. Thus, the MXD/24HA salt can potentially be a candidate drug for hair loss. The combination of theory and experiment enables an innovative methodology to understand the differences in structure and properties of pharmaceutical salts.

What are the current limits on determination of protonation state using neutron MX?

What are the current limits on determination of protonation state using neutron MX?

Role of water in the determination of protonation states of titratable residues.

Water is the fundamental unit for living being, and its contribution in variety of crucial cellular functions is widely accepted. The presence of water molecules in protein's environment also accounts for structural optimization, in which highly conserved water molecules ensure structural stability of the biomolecule by providing protein-water (solute-solvent) hydrogen-bond interaction networks. Similarly, protonation states and pKa values of individual amino acid residues are also influenced by neighboring water molecules present in the protein's vicinity. In the present study, we have highlighted the role of water molecules in hydrogen-bond optimization, in determining pKa values and protonation states of titratable residues in JH2 domain of JAK2 apo protein. We found that inclusion or exclusion of water molecules while calculating pKa and assigning protonation states to amino acid residues during the molecular system build-up step resulted in slight differences in pKa values of few titratable residues and alternative protonation states of a certain residue. Accordingly, different protonation states of ionizable residues offer differing interaction patterns. Thus, we inferred that the presence of water optimizes the hydrogen-bond interactions by forming direct protein-water interactions and by linking via protein-protein bridging interactions. However, in the absence of water, the interaction pattern is somewhat disrupted. We assume that water molecules could modulate the plausibility of a particular protonation state of titratable residues on the basis of its fit with the local environment, by utilizing some particular hydrogen-bond contacts that would remain unexploited in the absence of water.

What are the current limits on determination of protonation state using neutron macromolecular crystallography?

The rate of deposition of models determined by neutron diffraction, or a hybrid approach that combines X-ray and neutron diffraction, has increased in recent years. The benefit of neutron diffraction is that hydrogen atom (H) positions are detectable, allowing for the determination of protonation state and water molecule orientation. This study analyses all neutron models deposited in the Protein Data Bank to date, focusing on protonation state and properties of H (or deuterium, D) atoms as well as the details of water molecules. In particular, clashes and hydrogen bonds involving H or D atoms are investigated. As water molecules are typically the least reproducible part of a structural model, their positions in neutron models were compared to those in homologous high-resolution X-ray structures. For models determined by joint refinement against X-ray and neutron data, the water structure comparison was also carried out for models re-refined against the X-ray data alone. The homologues have generally fewer conserved water molecules where X-ray only was used and the positions of equivalent waters vary more than in the case of the hybrid X-ray model. As neutron diffraction data are generally less complete than X-ray data, the influence of neutron data completeness on nuclear density maps was also analyzed. We observe and discuss systematic map quality deterioration as result of data incompleteness.

Evaluation of models determined by neutron diffraction and proposed improvements to their validation and deposition.

The Protein Data Bank (PDB) contains a growing number of models that have been determined using neutron diffraction or a hybrid method that combines X-ray and neutron diffraction. The advantage of neutron diffraction experiments is that the positions of all atoms can be determined, including H atoms, which are hardly detectable by X-ray diffraction. This allows the determination of protonation states and the assignment of H atoms to water molecules. Because neutrons are scattered differently by hydrogen and its isotope deuterium, neutron diffraction in combination with H/D exchange can provide information on accessibility, dynamics and chemical lability. In this study, the deposited data, models and model-to-data fit for all PDB entries that used neutron diffraction as the source of experimental data have been analysed. In many cases, the reported Rwork and Rfree values were not reproducible. In such cases, the model and data files were analysed to identify the reasons for this mismatch. The issues responsible for the discrepancies are summarized and explained. The analysis unveiled limitations to the annotation, deposition and validation of models and data, and a lack of community-wide accepted standards for the description of neutron models and data, as well as deficiencies in current model refinement tools. Most of the issues identified concern the handling of H atoms. Since the primary use of neutron macromolecular crystallography is to locate and directly visualize H atoms, it is important to address these issues, so that the deposited neutron models allow the retrieval of the maximum amount of information with the smallest effort of manual intervention. A path forward to improving the annotation, validation and deposition of neutron models and hybrid X-ray and neutron models is suggested.

Red/Green Color Tuning of Visual Rhodopsins: Electrostatic Theory Provides a Quantitative Explanation

We present a structure-based theory of the long-wavelength (red/green) color tuning in visual rhodopsins and its application to the analysis of site-directed mutagenesis experiments. Using a combination of electrostatic and molecular-mechanics methods, we explain the measured mutant-minus-wild-type absorption shifts and conclude that the dominant mechanism of the color tuning in these systems is electrostatic pigment-protein coupling. An important element of our analysis is the independent determination of protonation states of titratable residues in the wild type and the mutant protein as well as the self-consistent reoptimization of hydrogen atom positions, which includes the relaxation of the hydrogen bonding network and the reorientation of water molecules. On the basis of this analysis, we propose a "dipole-orientation rule" according to which both the position and the orientation of a polar group introduced in the protein environment determine the direction of the transition energy shift of the retinal chromophore.

Determination of protonation states of iminosugar-enzyme complexes using photoinduced electron transfer.

A series of N-alkylated analogues of 1-deoxynojirimycin containing a fluorescent 10-chloro-9-anthracene group in the N-alkyl substituent were prepared. The anthracene group acted as a reporting group for protonation at the nitrogen in the iminosugar because an unprotonated amine was found to quench fluorescence by photoinduced electron transfer. The new compounds were found to inhibit β-glucosidase from Phanerochaete chrysosporium and α-glucosidase from Aspergillus niger, with Ki values in the low micro- to nanomolar range. Fluorescence and inhibition versus pH studies of the β-glucosidase-iminosugar complexes revealed that the amino group in the inhibitor is unprotonated when bound, while one of the active site carboxylates is protonated.

Direct determination of protonation states and visualization of hydrogen bonding in a glycoside hydrolase with neutron crystallography

Glycoside hydrolase (GH) enzymes apply acid/base chemistry to catalyze the decomposition of complex carbohydrates. These ubiquitous enzymes accept protons from solvent and donate them to substrates at close to neutral pH by modulating the pKa values of key side chains during catalysis. However, it is not known how the catalytic acid residue acquires a proton and transfers it efficiently to the substrate. To better understand GH chemistry, we used macromolecular neutron crystallography to directly determine protonation and ionization states of the active site residues of a family 11 GH at multiple pD (pD=pH+0.4) values. The general acid glutamate (Glu) cycles between two conformations, upward and downward, but is protonated only in the downward orientation. We performed continuum electrostatics calculations to estimate the pKa values of the catalytic Glu residues in both the apo- and substrate-bound states of the enzyme. The calculated pKa of the Glu increases substantially when the side chain moves down. The energy barrier required to rotate the catalytic Glu residue back to the upward conformation, where it can protonate the glycosidic oxygen of the substrate, is 4.3 kcal/mol according to free energy simulations. These findings shed light on the initial stage of the glycoside hydrolysis reaction in which molecular motion enables the general acid catalyst to obtain a proton from the bulk solvent and deliver it to the glycosidic oxygen.

BioDiff - a neutron diffractometer for protein crystallography

The research reactor Heinz Maier-Leibnitz (FRM II) is a modern high flux neutron source which feeds some 30 state of the art neutron beam instruments. Currently 24 are operational, others in commissioning or under construction. The newly built neutron single crystal diffractometer BIODIFF is especially designed to collect data from crystals with large unit cells. The main field of application is the structural analysis of proteins, especially the determination of hydrogen atom positions. BIODIFF is a joint project of the Jülich Centre for Neutron Science (JCNS) and the FRM II. Typical scientific questions addressed are the determination of protonation states of amino acid side chains (see e. g. [1,2]) and the characterization of the hydrogen bonding network between the protein active centre and an inhibitor or substrate. BIODIFF is designed as a monochromatic instrument. By using a highly orientated pyrolytic graphite monochromator (PG002) the diffractometer is able to operate in the wavelength range of 2.4 Å to about 5.6 Å. Contaminations of higher order wavelengths are removed by a neutron velocity selector. To cover a large solid angle the main detector of BIODIFF consists of a neutron imaging plate in a cylindrical geometry with online read-out capability. A fast Li/ZnS scintillator CCD camera is available for additional detection abilities. An optical CCD-camera pointing at the sample position is used to quickly align the crystal with respect to the neutron beam. The main advantage of BIODIFF is the possibility to adapt the wavelength to the size of the unit cell of the sample crystal while operating with a clean monochromatic beam that keeps the background level low. BIODIFF is equipped with a standard Oxford Cryosystem "Cryostream 700+" which allows measurements in the temperature regime from 90 K up to 500 K (see Figure underneath).

Looking for Hydrogen Atoms: Neutron Crystallography Provides Novel Insights Into Protein Structure and Function

Neutron crystallography allows direct localization of hydrogen positions in biological macromolecules. Within enzymes, hydrogen atoms play a pivotal role in catalysis. Recent advances in instrumentation and sample preparation have helped to overcome the difficulties of performing neutron diffraction experiments on protein crystals. The application of neutron macromolecular crystallography to a growing number of proteins has yielded novel structural insights. The ability to accurately position water molecules, hydronium ions, and hydrogen atoms within protein structures has helped in the study of low-barrier hydrogen bonds and hydrogen-bonding networks. The determination of protonation states of protein side chains, substrates, and inhibitors in the context of the macromolecule has provided important insights into enzyme chemistry and ligand binding affinities, which can assist in the design of potent therapeutic agents. In this review, we give an overview of the method and highlight advances in knowledge attained through the application of neutron protein crystallography.

The protonation states of oxo-bridged Mn(IV) dimers resolved by experimental and computational Mn K pre-edge X-ray absorption spectroscopy.

In nature, the protonation of oxo bridges is a commonly encountered mechanism for fine-tuning chemical properties and reaction pathways. Often, however, the protonation states are difficult to establish experimentally. This is of particular importance in the oxygen evolving complex of photosystem II, where identification of the bridging oxo protonation states is one of the essential requirements toward unraveling the mechanism. In order to establish a combined experimental and theoretical protocol for the determination of protonation states, we have systematically investigated a series of Mn model complexes by Mn K pre-edge X-ray absorption spectroscopy. An ideal test case for selective bis-μ-oxo-bridge protonation in a Mn dimer is represented by the system [Mn(IV)2(salpn)2(μ-OHn)2](n+). Although the three species [Mn(IV)2(salpn)2(μ-O)2], [Mn(IV)2(salpn)2(μ-O)(μ-OH)](+) and [Mn(IV)2(salpn)2(μ-OH)2](2+) differ only in the protonation of the oxo bridges, they exhibit distinct differences in the pre-edge region while maintaining the same edge energy. The experimental spectra are correlated in detail to theoretically calculated spectra. A time-dependent density functional theory approach for calculating the pre-edge spectra of molecules with multiple metal centers is presented, using both high spin (HS) and broken symmetry (BS) electronic structure solutions. The most intense pre-edge transitions correspond to an excitation of the Mn 1s core electrons into the unoccupied orbitals of local e(g) character (d(z)(2) and d(xy) based in the chosen coordinate system). The lowest energy experimental feature is dominated by excitations of 1s-α electrons, and the second observed feature is primarily attributed to 1s-β electron excitations. The observed energetic separation is due to spin polarization effects in spin-unrestricted density functional theory and models final state multiplet effects. The effects of spin polarization on the calculated Mn K pre-edge spectra, in both the HS and BS solutions, are discussed in terms of the strength of the antiferromagnetic coupling and associated changes in the covalency of Mn-O bonds. The information presented in this paper is complemented with the X-ray emission spectra of the same compounds published in an accompanying paper. Taken together, the two studies provide the foundation for a better understanding of the X-ray spectroscopic data of the oxygen evolving complex (OEC) in photosystem II.

A QM/MM study on the catalytic mechanism of pyruvate decarboxylase

Pyruvate decarboxylase (PDC) is a typical thiamin diphosphate (ThDP)-dependent enzyme with widespread applications in industry. Though studies regarding the reaction mechanism of PDC have been reported, they are mainly focused on the formation of ThDP ylide and some elementary steps in the catalytic cycle, studies about the whole catalytic cycle of PDC are still not completed. In these previous studies, a major controversy is whether the key active residues (Glu473, Glu50′, Asp27′, His113′, His114′) are protonated or ionized during the reaction. To explore the catalytic mechanism and the role of key residues in the active site, three whole-enzyme models were considered, and the combined QM/MM calculations on the nonoxidative decarboxylation of pyruvate to acetaldehyde catalyzed by PDC were performed. According to our computational results, the fundamental reaction pathways, the complete energy profiles of the whole catalytic cycle, and the specific role of key residues in the common steps were obtained. It is also found that the same residue with different protonation states will lead to different reaction pathways and energy profiles. The mechanism derived from the model in which the residues (Glu473, Glu50′, Asp27′, His113′, His114′) are in their protonated states is most consistent with experimental observations. Therefore, extreme care must be taken when assigning the protonation states in the mechanism study. Because the experimental determination of protonation state is currently difficult, the combined QM/MM method provides an indirect means for determining the active-site protonation state.

Thermodynamic integration to predict host-guest binding affinities

An alchemical free energy method with explicit solvent molecular dynamics simulations was applied as part of the blind prediction contest SAMPL3 to calculate binding free energies for seven guests to an acyclic cucurbit-[n]uril host. The predictions included determination of protonation states for both host and guests, docking pose generation, and binding free energy calculations using thermodynamic integration. We found a root mean square error (RMSE) of 3.6 kcal mol(-1) from the reference experimental results, with an R(2) correlation of 0.51. The agreement with experiment for the largest contributor to this error, guest 6, is improved by 1.7 kcal mol(-1) when a periodicity-induced free energy correction is applied. The corrections for the other ligands were significantly smaller, and altogether the RMSE was reduced by 0.4 kcal mol(-1). We link properties of the host-guest systems during simulation to errors in the computed free energies. Overall, we show that charged host-guest systems studied here, initialized in unconfirmed docking poses, present a challenge to accurate alchemical simulation methods.

Salt or Co-Crystal? Determination of Protonation State by X-Ray Photoelectron Spectroscopy (XPS)

Combined (15)N ssNMR and X-ray photoelectron spectroscopy (XPS) investigations for theophylline, a theophylline co-crystal, and a theophyllinium salt demonstrate that XPS allows direct observation of the degree of proton transfer, and thus identification of whether a salt or a co-crystal has been formed. The presence of a strongly binding-energy-shifted N 1s XPS peak with protonation indicates a salt (C==NH(+)), while this peak is unmistakably absent in the co-crystal. XPS should be considered as an alternative and complementary technique to single crystal X-ray diffraction and solid-state nuclear magnetic resonance spectroscopy (ssNMR).

Direct Determination of Protonation States of Histidine Residues in a 2 Å Neutron Structure of Deoxy-Human Normal Adult Hemoglobin and Implications for the Bohr Effect

We have investigated the protonation states of histidine residues (potential Bohr groups) in the deoxy form (T state) of human hemoglobin by direct determination of hydrogen (deuterium) positions with the neutron protein crystallography technique. The reversible binding of protons is key to the allosteric regulation of human hemoglobin. The protonation states of 35 of the 38 His residues were directly determined from neutron scattering omit maps, with 3 of the remaining residues being disordered. Protonation states of 5 equivalent His residues—αHis20, αHis50, αHis89, βHis143, and βHis146—differ between the symmetry-related globin subunits. The distal His residues, αHis58 and βHis63, are protonated in the α1β1 heterodimer and are neutral in α2β2. Buried residue αHis103 is found to be protonated in both subunits. These distal and buried residues have the potential to act as Bohr groups. The observed protonation states of His residues are compared to changes in their p K a values during the transition from the T to the R state and the results provide some new insights into our understanding of the molecular mechanism of the Bohr effect.

The determination of protonation states in proteins

The determination of protonation states in proteins

The determination of protonation states in proteins

The protonation states of aspartic acids and glutamic acids as well as histidine are investigated in four X-ray cases: Ni,Ca concanavalin A at 0.94 A, a thrombin-hirugen binary complex at 1.26 A resolution and two thrombin-hirugen-inhibitor ternary complexes at 1.32 and 1.39 A resolution. The truncation of the Ni,Ca concanavalin A data at various test resolutions between 0.94 and 1.50 A provided a test comparator for the ;unknown' thrombin-hirugen carboxylate bond lengths. The protonation states of aspartic acids and glutamic acids can be determined (on the basis of convincing evidence) even to the modest resolution of 1.20 A as exemplified by our X-ray crystal structure refinements of Ni and Mn concanavalin A and also as indicated in the 1.26 A structure of thrombin, both of which are reported here. The protonation-state indication of an Asp or a Glu is valid provided that the following criteria are met (in order of importance). (i) The acidic residue must have a single occupancy. (ii) Anisotropic refinement at a minimum diffraction resolution of 1.20 A (X-ray data-to-parameter ratio of approximately 3.5:1) is required. (iii) Both of the bond lengths must agree with the expectation (i.e. dictionary values), thus allowing some relaxation of the bond-distance standard uncertainties required to approximately 0.025 A for a '3sigma' determination or approximately 0.04 A for a '2sigma' determination, although some variation of the expected bond-distance values must be allowed according to the microenvironment of the hydrogen of interest. (iv) Although the F(o) - F(c) map peaks are most likely to be unreliable at the resolution range around 1.20 A, if admitted as evidence the peak at the hydrogen position must be greater than or equal to 2.5 sigma and in the correct geometry. (v) The atomic B factors need to be less than 10 A(2) for bond-length differentiation; furthermore, the C=O bond can also be expected to be observed with continuous 2F(o) - F(c) electron density and the C-OH bond with discontinuous electron density provided that the atomic B factors are less than approximately 20 A(2) and the contour level is increased. The final decisive option is to carry out more than one experiment, e.g. multiple X-ray crystallography experiments and ideally neutron crystallography. The complementary technique of neutron protein crystallography has provided evidence of the protonation states of histidine and acidic residues in concanavalin A and also the correct orientations of asparagine and glutamine side chains. Again, the truncation of the neutron data at various test resolutions between 2.5 and 3.0 A, even 3.25 and 3.75 A resolution, examines the limits of the neutron probe. These various studies indicate a widening of the scope of both X-ray and neutron probes in certain circumstances to elucidate the protonation states in proteins.

Comparison of hydrogen determination with X-ray and neutron crystallography in a human aldose reductase–inhibitor complex

Protonation states determination by neutron (2.2 A at room temperature) and X-ray (0.66 A at 100 K) crystallographic studies were compared for a medium size enzyme, human aldose reductase (MW=36 kDa), complexed with its NADP+ coenzyme and a selected inhibitor of therapeutic interest. The neutron resolution could be achieved only with the ab initio fully deuterated protein and the subsequent crystallization in D2O of the complex. We used the largest good-quality crystal (1.00x0.67x0.23 mm, i.e. volume of 0.15 mm3) that we were able to grow so far. Both studies enable the determination of protonation states, with a clear advantage for neutrons in the case of less-ordered atoms (B>5 A2). Hydrogen atoms are best determined by a complementary analysis of the Fourier maps obtained from both methods.

Quantum chemical analysis of the interactions of transition state analogs with leucine aminopeptidase

AbstractThe physical nature of the intermolecular interactions between several leucine aminopeptidase inhibitors, transition state analogs differing in functional groups, and various constituents of the enzyme active site was analyzed using the hybrid variation–perturbation decomposition of self‐consistent field and second‐order Møller–Plesset perturbation theory interaction energies. The electrostatic term constitutes the dominant contribution in the total interaction energy, although the magnitude of the remaining terms—exchange, delocalization, and correlation—seems to be non‐negligible. The total MP2 interaction energy and its dominant electrostatic term correlate reasonably well with the experimentally measured activities of the inhibitors. The application of this method for activity prediction of leucine aminopeptidase inhibitors resulted in very good agreement between calculated and measured inhibition constant values. Results confirm that the applied approach can be a valuable tool for structure‐based drug design, prediction of binding affinities, determination of protonation state and binding mode in ligand–receptor systems. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem 84: 302–310, 2001