The advanced treatment of hydrogen bonding in quantum crystallography.
Ontology highlight
ABSTRACT:
SUBMITTER: Malaspina LA
PROVIDER: S-EPMC8202034 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Approximate wavefunctions can be improved by constraining them to reproduce observations derived from diffraction and scattering experiments. Conversely, charge density models, incorporating electron-density distributions, atomic positions and atomic motion, can be improved by supplementing diffraction experiments with quantum chemically calculated, tailor-made electron densities (form factors). In both cases quantum chemistry and diffraction/scattering experiments are combined into a single, integrated tool. The development of quantum crystallographic research is reviewed. Some results obtained by quantum crystallography illustrate the potential and limitations of this field.
| S-EPMC5576428 | biostudies-literature
Project description: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.
| S-EPMC4603456 | biostudies-literature
Project description:The structure of the O-methyl glycoside of the naturally occurring 6-O-[(R)-1-carboxyethyl]-α-D-galactopyranose, C10H18O8, has been determined by X-ray crystallography at 100 K, supplementing the previously determined structure obtained at 293 K (Acta Crystallogr.1996, C52, 2285-2287). Molecular dynamics simulations of this glycoside were performed in the crystal environment with different numbers of units cells included in the primary simulation system at both 100 and 293 K. The calculated unit cell parameters and the intramolecular geometries (bonds, angles, and dihedrals) agree well with experimental results. Atomic fluctuations, including B-factors and anisotropies, are in good agreement with respect to the relative values on an atom-by-atom basis. In addition, the fluctuations increase with increasing simulation system size, with the simulated values converging to values lower than those observed experimentally indicating that the simulation model is not accounting for all possible contributions to the experimentally observed B-factors, which may be related to either the simulation time scale or size. In the simulations, the hydroxyl group of O7 is found to form bifurcated hydrogen bonds with O6 and O8 of an adjacent molecule, with the interactions dominated by the HO7-O6 interaction. Quantum mechanical calculations support this observation.
| S-EPMC3771504 | biostudies-literature
Project description:Hirshfeld atom refinement (HAR) is a method which enables the user to obtain more accurate positions of hydrogen atoms bonded to light chemical elements using X-ray data. When data quality permits, this method can be extended to hydrogen-bonded transition metals (TMs), as in hydride complexes. However, addressing hydrogen thermal motions with HAR, particularly in TM hydrides, presents a challenge. At the same time, proper description of thermal vibrations can be vital for determining hydrogen positions correctly. In this study, we employ tools such as SHADE3 and Normal Mode Refinement (NoMoRe) to estimate anisotropic displacement parameters (ADPs) for hydrogen atoms during HAR and IAM refinements performed for seven structures of TM (Fe, Ni, Cr, Nb, Rh and Os) and metalloid (Sb) hydride complexes for which both the neutron and the X-ray structures have been determined. A direct comparison between neutron and HAR/SHADE3/NoMoRe ADPs reveals that the similarity between neutron hydrogen ADPs and those estimated with NoMoRe or SHADE3 is significantly higher than when hydrogen ADPs are refined with HAR. Regarding TM-H bond lengths, traditional HAR exhibits a slight advantage over the other methods. However, combining NoMoRe/SHADE3 with HAR results in a minor decrease in agreement with neutron TM-H bond lengths. For the Cr complex, for which high-resolution X-ray data were collected, an investigation of resolution-related effects was possible.
| S-EPMC10833390 | biostudies-literature
Project description:The strong force that originates from breaking covalent bonds can be easily quantified through various testing platforms, while weak interfacial sliding resistance (ISR), originating from hydrogen bonding or van der Waals (vdW) forces, is very challenging to measure. Facilitated by an in-house nanomechanical testing system, we are able to precisely quantify and clearly distinguish the interfacial interactions between individual carbon fibers and several substrates governed by either hydrogen bonding or vdW forces. The specific ISR of the interface dominated by vdW forces is 3.55 ± 0.50 μN mm-1 and it surprisingly increases to 157.86 ± 44.18 μN mm-1 if the interface is bridged by hydrogen bonding. The ad hoc studies demonstrate that hydrogen bonding rather than vdW forces has great potential in sewing the interface if both surfaces are supportive of the formation of hydrogen bonds. The findings will enlighten the engineering of interfacial interactions and further mediate the entire mechanical performance of structures.
| S-EPMC9053468 | biostudies-literature
Project description:Quantum crystallography is a fast-developing multidisciplinary area of crystallography. In this work, we analyse the influence of different charge density models (i.e., the multipole model (MM), Hirshfeld atom refinement (HAR), and the transferable aspherical atom model (TAAM)), modelling of the thermal motion of hydrogen atoms (anisotropic, isotropic, and with the aid of SHADE or NoMoRe), and the type of radiation used (Mo Kα and Cu Kα) on the final results. To achieve this aim, we performed a series of refinements against X-ray diffraction data for three model compounds and compared their final structures, geometries, shapes of ADPs, and charge density distributions. Our results were also supported by theoretical calculations that enabled comparisons of the lattice energies of these structures. It appears that geometrical parameters are better described (closer to the neutron values) when HAR is used; however, bonds to H atoms more closely match neutron values after MM or TAAM refinement. Our analysis shows the superiority of the NoMoRe method in the description of H-atom ADPs. Moreover, the shapes of the ADPs of H atoms, as well as their electron density distributions, were better described with low-resolution Cu Kα data in comparison to low-resolution Mo Kα data.
| S-EPMC8233966 | biostudies-literature
Project description:We have carried out extensive computational analyses of the structure and bonding mechanism in trihalides DX⋅⋅⋅A(-) and the analogous hydrogen-bonded complexes DH⋅⋅⋅A(-) (D, X, A=F, Cl, Br, I) using relativistic density functional theory (DFT) at zeroth-order regular approximation ZORA-BP86/TZ2P. One purpose was to obtain a set of consistent data from which reliable trends in structure and stability can be inferred over a large range of systems. The main objective was to achieve a detailed understanding of the nature of halogen bonds, how they resemble, and also how they differ from, the better understood hydrogen bonds. Thus, we present an accurate physical model of the halogen bond based on quantitative Kohn-Sham molecular orbital (MO) theory, energy decomposition analyses (EDA) and Voronoi deformation density (VDD) analyses of the charge distribution. It appears that the halogen bond in DX⋅⋅⋅A(-) arises not only from classical electrostatic attraction but also receives substantial stabilization from HOMO-LUMO interactions between the lone pair of A(-) and the σ* orbital of D-X.
| S-EPMC3922460 | biostudies-literature
Project description:The pH-dependence of the degree of hydrogen-bonding between a base and its conjugate acid is considered. When only a small proportion of the total base is complexed, the amount complexed is proportional to (1+coshp)(-1) where p=2.303 (pK(a)-pH), pK(a) being the dissociation constant of the conjugate acid. This represents sharp pH-dependence. As the proportion complexed increases, the curve broadens, eventually becoming flat-topped, with more than half the base complexed over the range of pH values pK(a)+/-logKC, approximately. (K is the complex association constant and C is the formal base concentration, including all forms.) There are similarities to the extent of mono-protonation of a dibasic acid.
| S-EPMC1168448 | biostudies-other
Project description:Plant lectins are carbohydrate-binding proteins with various biomedical applications. Concanavalin A (Con A) holds promise in treating cancerous tumors. To better understand the Con A carbohydrate binding specificity, we obtained a room-temperature neutron structure of this legume lectin in complex with a disaccharide Man?1-2Man, mannobiose. The neutron structure afforded direct visualization of the hydrogen bonding between the protein and ligand, showing that the ligand is able to alter both protonation states and interactions for residues located close to and distant from the binding site. An unprecedented low-barrier hydrogen bond was observed forming between the carboxylic side chains of Asp28 and Glu8, with the D atom positioned equidistant from the oxygen atoms having an O···D···O angle of 101.5°.
| S-EPMC5754002 | biostudies-literature
Project description:Archaeal family-B DNA polymerases stall replication on encountering the pro-mutagenic bases uracil and hypoxanthine. This publication describes an X-ray crystal structure of Thermococcus gorgonarius polymerase in complex with a DNA containing hypoxanthine in the single-stranded region of the template, two bases ahead of the primer-template junction. Full details of the specific recognition of hypoxanthine are revealed, allowing a comparison with published data that describe uracil binding. The two bases are recognized by the same pocket, in the N-terminal domain, and make very similar protein-DNA interactions. Specificity for hypoxanthine (and uracil) arises from a combination of polymerase-base hydrogen bonds and shape fit between the deaminated bases and the pocket. The structure with hypoxanthine at position 2 explains the stimulation of the polymerase 3'-5' proofreading exonuclease, observed with deaminated bases at this location. A beta-hairpin element, involved in partitioning the primer strand between the polymerase and exonuclease active sites, inserts between the two template bases at the extreme end of the double-stranded DNA. This denatures the two complementary primer bases and directs the resulting 3' single-stranded extension toward the exonuclease active site. Finally, the relative importance of hydrogen bonding and shape fit in determining selectivity for deaminated bases has been examined using nonpolar isosteres. Affinity for both 2,4-difluorobenzene and fluorobenzimidazole, non-hydrogen bonding shape mimics of uracil and hypoxanthine, respectively, is strongly diminished, suggesting polar protein-base contacts are important. However, residual interaction with 2,4-difluorobenzene is seen, confirming a role for shape recognition.
| S-EPMC2901926 | biostudies-literature