Protonation and geometry of histidine rings

A physicochemical rationale for the varied catalytic efficiency in RNase J paralogues

Paralogs of the bifunctional nuclease, Ribonuclease J (RNase J), demonstrate varied catalytic efficiencies despite extensive sequence and structural similarity. Of the two Staphylococcus aureus RNase J paralogues, RNase J1 is substantially more active than RNase J2. Mutational analysis of active site residues revealed that only H80 and E166 were critical for nuclease activity. Electronic properties of active site residues were further evaluated using density functional theory in conjunction with molecular mechanics. This analysis suggested that multiple residues at the active site can function as Lewis bases or acids in RNase J2. The bond dissociation energy, on the other hand, suggested that the Mn ion in RNase J2, located at a structurally identical location to that in RNase J1, is crucial for overall structural integrity. Structures of mutant enzymes lacking the metal ion were seen to adopt a different orientation between the substrate binding and catalytic domain than wild-type RNase J2. A surprising finding was that the RNase J2 H78 A mutant was five-fold more active than the wild-type enzyme. Structural and biochemical experiments performed in light of this observation revealed that the RNase J2 catalytic mechanism is distinct from both two-metal ion and one-metal ion reaction mechanisms proposed for RNase J nucleases. Different activity levels in RNase J paralogues can thus be ascribed to the diversity in catalytic mechanisms.

Protonation of histidine rings using quantum-mechanical methods

Histidine can be protonated on either or both of the two N atoms of the imidazole moiety. Each of the three possible forms occurs as a result of the stereochemical environment of the histidine side chain. In an atomic model, comparing the possible protonation states in situ, looking at possible hydrogen bonding and metal coordination, it is possible to predict which is most likely to be correct. A more direct method is described that uses quantum-mechanical methods to calculate, also in situ, the minimum geometry and energy for comparison, and therefore to more accurately identify the most likely protonation state.

Effect of histidine protonation state on ligand binding at the ATP-binding site of human protein kinase CK2

Histidine residues contribute to numerous molecular interactions, owing to their structure with the ionizable aromatic side chain with pKa close to the physiological pH. Herein, we studied how the two histidine residues, His115 and His160 of the catalytic subunit of human protein kinase CK2, affect the binding of the halogenated heterocyclic ligands at the ATP-binding site. Thermodynamic studies on the interaction between five variants of hCK2α (WT protein and four histidine mutants) and three ionizable bromo-benzotriazoles and their conditionally non-ionizable benzimidazole counterparts were performed with nanoDSF, MST, and ITC. The results allowed us to identify the contribution of interactions involving the particular histidine residues to ligand binding. We showed that despite the well-documented hydrogen bonding/salt bridge formation dragging the anionic ligands towards Lys68, the protonated His160 also contributes to the binding of such ligands by long-range electrostatic interactions. Simultaneously, His 115 indirectly affects ligand binding, placing the hinge region in open/closed conformations.

Phospholipid binding of the dengue virus envelope E protein segment containing the conserved His residue

Flaviviruses encompass many important human pathogens, including Dengue, Zika, West Nile, Yellow fever, Japanese encephalitis, and Tick-borne encephalitis viruses as well as several emerging viruses that affect millions of people worldwide. They enter cells by endocytosis, fusing their membrane with the late endosomal one in a pH-dependent manner, so membrane fusion is one of the main targets for obtaining new antiviral inhibitors. The envelope E protein, a class II membrane fusion protein, is responsible for fusion and contains different domains involved in the fusion mechanism, including the fusion peptide. However, other segments, apart from the fusion peptide, have been implicated in the mechanism of membrane fusion, in particular a segment containing a His residue supposed to act as a specific pH sensor. We have used atomistic molecular dynamics to study the binding of the envelope E protein segment containing the conserved His residue in its three different tautomer forms with a complex membrane mimicking the late-endosomal one. We show that this His-containing segment is capable of spontaneous membrane binding, preferentially binds electronegatively charged phospholipids and does not bind cholesterol. Since Flaviviruses have caused epidemics in the past, continue to do so and will undoubtedly continue to do so, this specific segment could characterise a new target that would allow finding effective antiviral molecules against DENV virus in particular and Flaviviruses in general.

Protonation State of an Important Histidine from High Resolution Structures of Lytic Polysaccharide Monooxygenases

Lytic Polysaccharide Monooxygenases (LPMOs) oxidatively cleave recalcitrant polysaccharides. The mechanism involves (i) reduction of the Cu, (ii) polysaccharide binding, (iii) binding of different oxygen species, and (iv) glycosidic bond cleavage. However, the complete mechanism is poorly understood and may vary across different families and even within the same family. Here, we have investigated the protonation state of a secondary co-ordination sphere histidine, conserved across AA9 family LPMOs that has previously been proposed to be a potential proton donor. Partial unrestrained refinement of newly obtained higher resolution data for two AA9 LPMOs and re-refinement of four additional data sets deposited in the PDB were carried out, where the His was refined without restraints, followed by measurements of the His ring geometrical parameters. This allowed reliable assignment of the protonation state, as also validated by following the same procedure for the His brace, for which the protonation state is predictable. The study shows that this histidine is generally singly protonated at the Nε2 atom, which is close to the oxygen species binding site. Our results indicate robustness of the method. In view of this and other emerging evidence, a role as proton donor during catalysis is unlikely for this His.

Docking studies and molecular dynamics simulation of triazole benzene sulfonamide derivatives with human carbonic anhydrase IX inhibition activity

Carbonic anhydrase IX has been used as a hypoxia endogenous marker in a range of solid tumors including renal cell, lung, bladder and tumors of the head and neck. α-CA IX isozyme is over-expressive in hypoxic environment which becomes an attractive target for the design of inhibitors' targeting cancer particularly, tumor progression and invasion. In the process of designing new leads for the inhibition of tumor-associated hCA IX, the best triazole benzene sulfonamide derivatives were obtained from the QSAR model published in the research paper as cited. The statistically validated QSAR model was utilized for bioactivity prediction of novel leads. Further the designed molecules having good scores were subjected to molecular docking studies and molecular dynamic simulation studies. Designed compounds 1, 2, 20, 24 and 27 have shown predicted bioactivity of 9.13, 9.65, 10.05, 10.03 and 10.104 logarithmic units respectively using QSAR model 2. The low energy conformations of the above compounds exhibited good Autodock binding energy scores (-8.1, -8.2, -8.1, -8.3 and -9.2 K cal mol-1) and interactions with Gln92, Thr200, Asn66 and His68. Desmond's molecular dynamics simulations studies for 100 ns of compound 27 compared to reference SLC0111 provided useful structural insights of human carbonic anhydrase IX inhibition. Compound 27 with new chemical structure displayed both hydrophobic and hydrophilic stable interactions in the active site. RMSD, RMSF, RoG, H-bond and SASA analysis confirmed the stable binding of compound 27 with 5FL4 structure. In addition, MM-PBSA and MM-GBSA also affirm the docking results. We propose the designed compound 27 (predicted Ki = ∼0.07 nM) as the best theoretical lead which may further be experimentally studied for selective inhibition.

Crystallographic models of SARS-CoV-2 3CLpro: in-depth assessment of structure quality and validation

The appearance at the end of 2019 of the new SARS-CoV-2 coronavirus led to an unprecedented response by the structural biology community, resulting in the rapid determination of many hundreds of structures of proteins encoded by the virus. As part of an effort to analyze and, if necessary, remediate these structures as deposited in the Protein Data Bank (PDB), this work presents a detailed analysis of 81 crystal structures of the main protease 3CLpro, an important target for the design of drugs against COVID-19. The structures of the unliganded enzyme and its complexes with a number of inhibitors were determined by multiple research groups using different experimental approaches and conditions; the resulting structures span 13 different polymorphs representing seven space groups. The structures of the enzyme itself, all determined by molecular replacement, are highly similar, with the exception of one polymorph with a different inter-domain orientation. However, a number of complexes with bound inhibitors were found to pose significant problems. Some of these could be traced to faulty definitions of geometrical restraints for ligands and to the general problem of a lack of such information in the PDB depositions. Several problems with ligand definition in the PDB itself were also noted. In several cases extensive corrections to the models were necessary to adhere to the evidence of the electron-density maps. Taken together, this analysis of a large number of structures of a single, medically important protein, all determined within less than a year using modern experimental tools, should be useful in future studies of other systems of high interest to the biomedical community.

An unprecedented insight into the catalytic mechanism of copper nitrite reductase from atomic-resolution and damage-free structures

Copper-containing nitrite reductases (CuNiRs), encoded by nirK gene, are found in all kingdoms of life with only 5% of CuNiR denitrifiers having two or more copies of nirK Recently, we have identified two copies of nirK genes in several α-proteobacteria of the order Rhizobiales including Bradyrhizobium sp. ORS 375, encoding a four-domain heme-CuNiR and the usual two-domain CuNiR (Br 2DNiR). Compared with two of the best-studied two-domain CuNiRs represented by the blue (AxNiR) and green (AcNiR) subclasses, Br 2DNiR, a blue CuNiR, shows a substantially lower catalytic efficiency despite a sequence identity of ~70%. Advanced synchrotron radiation and x-ray free-electron laser are used to obtain the most accurate (atomic resolution with unrestrained SHELX refinement) and damage-free (free from radiation-induced chemistry) structures, in as-isolated, substrate-bound, and product-bound states. This combination has shed light on the protonation states of essential catalytic residues, additional reaction intermediates, and how catalytic efficiency is modulated.

Hydrogen-Mediated Noncovalent Interactions in Solids: What Can NMR Crystallography Tell About?

Hydrogen atoms play a crucial role in the aggregation of organic (bio)molecules through diverse number of noncovalent interactions that they mediate, such as electrostatic in proton transfer systems, hydrogen bonding, and CH-π interactions, to mention only the most prominent. To identify and adequately describe such low-energy interactions, increasingly sensitive methods have been developed over time, among which quantum chemical computations have witnessed impressive advances in recent years. For reaching the present state-of-the-art, computations had to rely on a pool of relevant experimental data, needed at least for validation, if not also for other purposes. In the case of molecular crystals, the best illustration for the synergy between computations and experiment is given by the so-called NMR crystallography approach. Originally designed to increase the confidence level in crystal structure determination of organic compounds from powders, NMR crystallography is able now to offer also a wealth of information regarding the noncovalent interactions that drive molecules to pack in a given crystalline pattern or another. This is particularly true for the noncovalent interactions which depend on the exact location of labile hydrogen atoms in the system: in such cases, NMR crystallography represents a valuable characterization tool, in some cases complementing even the standard single-crystal X-ray diffraction technique. A concise introduction in the field is made in this mini-review, which is aimed at providing a comprehensive picture with respect to the current accuracy level reached by NMR crystallography in the characterization of hydrogen-mediated noncovalent interactions in organic solids. Different types of practical applications are illustrated with the example of molecular crystals studied by our research group, but references to other representative developments reported in the literature are also made. By summarizing the major concepts and methodological progresses, the present work is also intended to be a guide to the practical potential of this relatively recent analytical tool for the scientists working in areas where crystal engineering represents the main approach for rational design of novel materials.

A secreted fungal histidine- and alanine-rich protein regulates metal ion homeostasis and oxidative stress

Like pathogens, beneficial endophytic fungi secrete effector proteins to promote plant colonization, for example, through perturbation of host immunity. The genome of the root endophyte Serendipita indica encodes a novel family of highly similar, small alanine- and histidine-rich proteins, whose functions remain unknown. Members of this protein family carry an N-terminal signal peptide and a conserved C-terminal DELD motif. Here we report on the functional characterization of the plant-responsive DELD family protein Dld1 using a combination of structural, biochemical, biophysical and cytological analyses. The crystal structure of Dld1 shows an unusual, monomeric histidine zipper consisting of two antiparallel coiled-coil helices. Similar to other histidine-rich proteins, Dld1 displays varying affinity to different transition metal ions and undergoes metal ion- and pH-dependent unfolding. Transient expression of mCherry-tagged Dld1 in barley leaf and root tissue suggests that Dld1 localizes to the plant cell wall and accumulates at cell wall appositions during fungal penetration. Moreover, recombinant Dld1 enhances barley root colonization by S. indica, and inhibits H₂ O₂ -mediated radical polymerization of 3,3'-diaminobenzidine. Our data suggest that Dld1 has the potential to enhance micronutrient accessibility for the fungus and to interfere with oxidative stress and reactive oxygen species homeostasis to facilitate host colonization.

Deuteron Solid-State NMR Relaxation Measurements Reveal Two Distinct Conformational Exchange Processes in the Disordered N-Terminal Domain of Amyloid-β Fibrils

We employed deuterium solid-state NMR techniques under static conditions to discern the details of the μs-ms timescale motions in the flexible N-terminal subdomain of Aβ1-40 amyloid fibrils, which spans residues 1-16. In particular, we utilized a rotating frame (R1ρ ) and the newly developed time domain quadrupolar Carr-Purcell-Meiboom-Gill (QCPMG) relaxation measurements at the selectively deuterated side chains of A2, H6, and G9. The two experiments are complementary in terms of probing somewhat different timescales of motions, governed by the tensor parameters and the sampling window of the magnetization decay curves. The results indicated two mobile "free" states of the N-terminal domain undergoing global diffusive motions, with isotropic diffusion coefficients of 0.7-1 ⋅ 108 and 0.3-3 ⋅ 106 ad2  s-1 . The free states are also involved in the conformational exchange with a single bound state, in which the diffusive motions are quenched, likely due to transient interactions with the structured hydrophobic core. The conformational exchange rate constants are 2-3 ⋅ 105  s-1 and 2-3 ⋅ 104  s-1 for the fast and slow diffusion free states, respectively.

Discrimination of Protein Amino Acid or Its Protonated State at Single-Residue Resolution by Graphene Nanopores

The function of a protein is determined by the composition of amino acids and is essential to proteomics. However, protein sequencing remains challenging due to the protein's irregular charge state and its high-order structure. Here, a proof of principle study on the capability of protein sequencing by graphene nanopores integrated with atomic force microscopy is performed using molecular dynamics simulations. It is found that nanopores can discriminate a protein sequence and even its protonation state at single-residue resolution. Both the pulling forces and current blockades induced by the permeation of protein residues are found to be highly correlated with the type of amino acids, which makes the residues identifiable. It is also found that aside from the dimension, both the conformation and charge state of the residue can significantly influence the force and current signal during its permeation through the nanopore. In particular, due to the electro-osmotic flow effect, the blockade current for the double-protonated histidine is slightly smaller than that for single-protonated histidine, which makes it possible for discrimination of different protonation states of amino acids. The results reported here present a novel protein sequencing scheme using graphene nanopores combined with nanomanipulation technology.

Solid-state NMR reveals a comprehensive view of the dynamics of the flexible, disordered N-terminal domain of amyloid-β fibrils

Amyloid fibril deposits observed in Alzheimer's disease comprise amyloid-β (Aβ) protein possessing a structured hydrophobic core and a disordered N-terminal domain (residues 1-16). The internal flexibility of the disordered domain is likely essential for Aβ aggregation. Here, we used ²H static solid-state NMR methods to probe the dynamics of selected side chains of the N-terminal domain of Aβ_1-40 fibrils. Line shape and relaxation data suggested a two-state model in which the domain's free state undergoes a diffusive motion that is quenched in the bound state, likely because of transient interactions with the structured C-terminal domain. At 37 °C, we observed freezing of the dynamics progressively along the Aβ sequence, with the fraction of the bound state increasing and the rate of diffusion decreasing. We also found that without solvation, the diffusive motion is quenched. The solvent acted as a plasticizer reminiscent of its role in the onset of global dynamics in globular proteins. As the temperature was lowered, the fraction of the bound state exhibited sigmoidal behavior. The midpoint of the freezing curve coincided with the bulk solvent freezing for the N-terminal residues and increased further along the sequence. Using ²H R _1ρ measurements, we determined the conformational exchange rate constant between the free and bound states under physiological conditions. Zinc-induced aggregation leads to the enhancement of the dynamics, manifested by the faster conformational exchange, faster diffusion, and lower freezing-curve midpoints.

Structure Refinement at Atomic Resolution

X-Ray diffraction data at atomic resolution, i.e., beyond 1.2 Å, provide the most detailed and reliable information we have about the structure of macromolecules, which is especially important for validating new discoveries and resolving subtle issues of molecular mechanisms. Refinement at atomic resolution allows reliable interpretation of static disorder and solvent structure, as well as modeling of anisotropic atomic vibrations and even of H atoms. Stereochemical restraints can be relaxed or removed, providing unbiased information about macromolecular stereochemistry, which in turn can be used to define improved conformation-dependent libraries, and the surplus of data allows estimation of least-squares uncertainties in the derived parameters. At ultrahigh resolution it is possible to study charge density distribution by multipolar refinement of electrons in non-spherical orbitals.

Neutron structure of the T26H mutant of T4 phage lysozyme provides insight into the catalytic activity of the mutant enzyme and how it differs from that of wild type

T4 phage lysozyme is an inverting glycoside hydrolase that degrades the murein of bacterial cell walls by cleaving the β‐1,4‐glycosidic bond. The substitution of the catalytic Thr26 residue to a histidine converts the wild type from an inverting to a retaining enzyme, which implies that the original general acid Glu11 can also act as an acid/base catalyst in the hydrolysis. Here, we have determined the neutron structure of the perdeuterated T26H mutant to clarify the protonation states of Glu11 and the substituted His26, which are key in the retaining reaction. The 2.09‐Å resolution structure shows that the imidazole group of His26 is in its singly protonated form in the active site, suggesting that the deprotonated Nɛ2 atom of His26 can attack the anomeric carbon of bound substrate as a nucleophile. The carboxyl group of Glu11 is partially protonated and interacts with the unusual neutral state of the guanidine moiety of Arg145, as well as two heavy water molecules. Considering that one of the water‐binding sites has the potential to be occupied by a hydronium ion, the bulk solvent could be the source for the protonation of Glu11. The respective protonation states of Glu11 and His26 are consistent with the bond lengths determined by an unrestrained refinement of the high‐resolution X‐ray structure of T26H at 1.04‐Å resolution. The detail structural information, including the coordinates of the deuterium atoms in the active site, provides insight into the distinctively different catalytic activities of the mutant and wild type enzymes.

PDB Code(s): 5XPE; 5XPF

Stereochemistry and Validation of Macromolecular Structures

Macromolecular structure is governed by the strict rules of stereochemistry. Several approaches to the validation of the correctness of the interpretation of crystallographic and NMR data that underlie the models deposited in the PDB are utilized in practice. The stereochemical rules applicable to macromolecular structures are discussed in this chapter. Practical, computer-based methods and tools of verification of how well the models adhere to those established structural principles to assure their quality are summarized.

Unifying the microscopic picture of His-containing turns: from gas phase model peptides to crystallized proteins

The presence in crystallized proteins of a local anchoring between the side chain of a His residue, located in the central position of a γ- or β-turn, and its local main chain environment, is assessed by the comparison of protein structures with relevant isolated model peptides.

Geometry of guanidinium groups in arginines

The restraints in common usage today have been obtained based on small molecule X-ray crystal structures available 25 years ago and recent reports have shown that the values of bond lengths and valence angles can be, in fact, significantly different from those stored in libraries, for example for the peptide bond or the histidine ring geometry. We showed that almost 50% of outliers found in protein validation reports released in the Protein Data Bank on 23 March 2016 come from geometry of guanidine groups in arginines. Therefore, structures of small molecules and atomic resolution protein crystal structures have been used to derive new target values for the geometry of this group. The most significant difference was found for NE-CZ-NH1 and NE-CZ-NH2 angles, showing that the guanidinium group is not symmetric. The NE-CZ-NH1 angle is larger, 121.5(10)˚, than NE-CZ-NH2, 119.2(10)˚, due to the repulsive interaction between NH1 and CD1 atom.

Transferable aspherical atom model refinement of protein and DNA structures against ultrahigh-resolution X-ray data

In contrast to the independent-atom model (IAM), in which all atoms are assumed to be spherical and neutral, the transferable aspherical atom model (TAAM) takes into account the deformed valence charge density resulting from chemical bond formation and the presence of lone electron pairs. Both models can be used to refine small and large molecules, e.g. proteins and nucleic acids, against ultrahigh-resolution X-ray diffraction data. The University at Buffalo theoretical databank of aspherical pseudo-atoms has been used in the refinement of an oligopeptide, of Z-DNA hexamer and dodecamer duplexes, and of bovine trypsin. The application of the TAAM to these data improves the quality of the electron-density maps and the visibility of H atoms. It also lowers the conventional R factors and improves the atomic displacement parameters and the results of the Hirshfeld rigid-bond test. An additional advantage is that the transferred charge density allows the estimation of Coulombic interaction energy and electrostatic potential.