The occurrence of CH…O hydrogen bonds in α-helices and helix termini in globular proteins

Differing Effects of Nonlinearity around the Proton Acceptor on CH··O and NH··O H-Bond Strength within Proteins

The effects of deviations from nonlinearity around the carbonyl proton acceptor of an amide group are assessed by DFT quantum chemical calculations for both CH··O and NH··O H-bonds. The proton donors are the imidazole functional group of His and the indole of Trp, which are paired respectively with N-methylacetamide and acetamide. The displacement of either CH or NH group toward the carbonyl O sp2 lone pairs stabilizes the system and strengthens the H-bond. But the two donor groups differ in their response to a shift out of the amide plane. While the NH··O H-bond is weakened by this displacement, a substantial strengthening is observed when the CH donor is moved out of this plane, in one direction versus the other. This pattern is explained on the basis of simple Coulombic considerations.

Intramolecular CH-Hydrogen Bonding During the Dissociation of the Oxaphosphetane Intermediate Facilitates Z/E-Selectivity in Wittig Olefination

Herein, DFT studies corroborating experimental results revealed that the shortest intramolecular hydrogen bonding distance of cis/trans-oxaphosphetane (OPA) oxygen with the CH-hydrogen of a triphenylphosphine phenyl ring provides good evidence for the attained olefin Z/E-selectivity in Wittig olefination of the studied examples. 2-Nitrobenzaldehyde, 3-nitrobenzaldehyde, 2-nitro-3-bromobenzaldehyde, 2-nitro-5-bromobenzaldehyde and 2-nitro-5-arylbenzaldehydes provided Z-nitrostilbenes with (2-chloro-4-hydroxy-3-methoxy-5-(methoxycarbonyl)benzyl) triphenylphosphonium chloride as the major products. However, 4-nitrobenzaldehyde and 2-nitro-6-bromobenzaldehydes furnished E-nitrostilbenes as the major products in high yields. Furthermore, the DFT computed intramolecular CH1/CH2-hydrogen bond distances with Cl/NO2 of selected stilbene derivatives were in good agreement with intramolecular hydrogen bond distances measured from single-crystal X-ray diffraction measurements.

The active site of the SGNH hydrolase-like fold proteins: Nucleophile-oxyanion (Nuc-Oxy) and Acid-Base zones

SGNH hydrolase-like fold proteins are serine proteases with the default Asp-His-Ser catalytic triad. Here, we show that these proteins share two unique conserved structural organizations around the active site: (1) the Nuc-Oxy Zone around the catalytic nucleophile and the oxyanion hole, and (2) the Acid-Base Zone around the catalytic acid and base. The Nuc-Oxy Zone consists of 14 amino acids cross-linked with eight conserved intra- and inter-block hydrogen bonds. The Acid-Base Zone is constructed from a single fragment of the polypeptide chain, which incorporates both the catalytic acid and base, and whose N- and C-terminal residues are linked together by a conserved hydrogen bond. The Nuc-Oxy and Acid-Base Zones are connected by an SHLink, a two-bond conserved interaction from amino acids, adjacent to the catalytic nucleophile and base.

A Highly Ordered Nitroxide Side Chain for Distance Mapping and Monitoring Slow Structural Fluctuations in Proteins

Site-directed spin labeling electron paramagnetic resonance (SDSL-EPR) is an established tool for exploring protein structure and dynamics. Although nitroxide side chains attached to a single cysteine via a disulfide linkage are commonly employed in SDSL-EPR, their internal flexibility complicates applications to monitor slow internal motions in proteins and to structure determination by distance mapping. Moreover, the labile disulfide linkage prohibits the use of reducing agents often needed for protein stability. To enable the application of SDSL-EPR to the measurement of slow internal dynamics, new spin labels with hindered internal motion are desired. Here, we introduce a highly ordered nitroxide side chain, designated R9, attached at a single cysteine residue via a non-reducible thioether linkage. The reaction to introduce R9 is highly selective for solvent-exposed cysteine residues. Structures of R9 at two helical sites in T4 Lysozyme were determined by X-ray crystallography and the mobility in helical sequences was characterized by EPR spectral lineshape analysis, Saturation Transfer EPR, and Saturation Recovery EPR. In addition, interspin distance measurements between pairs of R9 residues are reported. Collectively, all data indicate that R9 will be useful for monitoring slow internal structural fluctuations, and applications to distance mapping via dipolar spectroscopy and relaxation enhancement methods are anticipated.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00723-023-01618-8.

C-H Groups as Donors in Hydrogen Bonds: A Historical Overview and Occurrence in Proteins and Nucleic Acids

Hydrogen bonds constitute a unique type of non-covalent interaction, with a critical role in biology. Until fairly recently, the canonical view held that these bonds occur between electronegative atoms, typically O and N, and that they are mostly electrostatic in nature. However, it is now understood that polarized C-H groups may also act as hydrogen bond donors in many systems, including biological macromolecules. First recognized from physical chemistry studies, C-H…X bonds were visualized with X-ray crystallography sixty years ago, although their true significance has only been recognized in the last few decades. This review traces the origins of the field and describes the occurrence and significance of the most important C-H…O bonds in proteins and nucleic acids.

Sulfur-mediated chalcogen versus hydrogen bonds in proteins: a see-saw effect in the conformational space

Divalent sulfur (S) forms a chalcogen bond (Ch-bond) via its σ-holes and a hydrogen bond (H-bond) via its lone pairs. The relevance of these interactions and their interplay for protein structure and function is unclear. Based on the analyses of the crystal structures of small organic/organometallic molecules and proteins and their molecular electrostatic surface potential, we show that the reciprocity of the substituent-dependent strength of the σ-holes and lone pairs correlates with the formation of either Ch-bond or H-bond. In proteins, cystines preferentially form Ch-bonds, metal-chelated cysteines form H-bonds, while methionines form either of them with comparable frequencies. This has implications for the positioning of these residues and their role in protein structure and function. Computational analyses reveal that the S-mediated interactions stabilise protein secondary structures by mechanisms such as helix capping and protecting free β-sheet edges by negative design. The study highlights the importance of S-mediated Ch-bond and H-bond for understanding protein folding and function, the development of improved strategies for protein/peptide structure prediction and design and structure-based drug discovery.

Design and Synthesis of a Ferrocene-Based Diol Library and Application in the Hetero-Diels-Alder Reaction

Diol scaffolds have been utilized as highly effective catalysts and ligands in a wide range of catalytic asymmetric transformations. For scaffolds to be successful as broadly used motifs, they should be prepared cheaply through facile routes and be easily handled. Herein, the synthesis of a family of planar chiral diols based on a modular and robust three-step route is described, which yields catalytically active diols in >99 % de and >99 % ee, with up to seven different chiral elements. These diols have been characterized by X-ray crystallographic analysis, which provides clear evidence for the likely transition state when applied in the asymmetric hetero-Diels-Alder reaction. Without altering the stereochemistry of the catalyst backbone, it is possible to access both enantiomers of the product by varying the substitution of the catalyst at the α-position.

Proline C-H Bonds as Loci for Proline Assembly via C-H/O Interactions

Proline residues within proteins lack a traditional hydrogen bond donor. However, the hydrogens of the proline ring are all sterically accessible, with polarized C-H bonds at Hα and Hδ that exhibit greater partial positive character and can be utilized as alternative sites for molecular recognition. C-H/O interactions, between proline C-H bonds and oxygen lone pairs, have been previously identified as modes of recognition within protein structures and for higher-order assembly of protein structures. In order to better understand intermolecular recognition of proline residues, a series of proline derivatives was synthesized, including 4R-hydroxyproline nitrobenzoate methyl ester, acylated on the proline nitrogen with bromoacetyl and glycolyl groups, and Boc-4S-(4-iodophenyl)hydroxyproline methyl amide. All three derivatives exhibited multiple close intermolecular C-H/O interactions in the crystallographic state, with H⋅⋅⋅O distances as close as 2.3 Å. These observed distances are well below the 2.72 Å sum of the van der Waals radii of H and O, and suggest that these interactions are particularly favorable. In order to generalize these results, we further analyzed the role of C-H/O interactions in all previously crystallized derivatives of these amino acids, and found that all 26 structures exhibited close intermolecular C-H/O interactions. Finally, we analyzed all proline residues in the Cambridge Structural Database of small-molecule crystal structures. We found that the majority of these structures exhibited intermolecular C-H/O interactions at proline C-H bonds, suggesting that C-H/O interactions are an inherent and important mode for recognition of and higher-order assembly at proline residues. Due to steric accessibility and multiple polarized C-H bonds, proline residues are uniquely positioned as sites for binding and recognition via C-H/O interactions.

Molecular duplexes featuring NH···N, CH···O and CH···π interactions in solid-state self-assembly of triazine-based compounds

Synthetic supramolecular structures constructed through the cooperative action of numerous non-covalent forces are highly desirable as models to unravel and understand the complexity of systems created in nature via self-assembly. Taking advantage of the low cost of 2,4,6-trichloro-1,3,5-triazine (cyanuric chloride) and the sequential nucleophilic substitution reactions with almost all types of nucleophiles, a series of six structurally related novel s-triazine derivatives 1-6 were synthesized and structurally characterized based on their physical, spectral and crystallographic data. The solid-state structures of all the six compounds showed intriguing and unique molecular duplexes featuring NH···N, CH···O and CH···π interactions. Careful analysis of different geometric parameters of the involved H-bonds indicates that they are linear, significant and are therefore responsible for guiding the three-dimensional structure of these compounds in the solid state. The prevalence of sextuple hydrogen bond array-driven molecular duplexes and the possibility of structural modifications on the s-triazine ring render these novel triazine derivatives 1-6 attractive as a platform to create heteroduplex constructs and their subsequent utility in the field of supramolecular chemistry and crystal engineering.

Supramolecular Self-Assembly Mediated by Multiple Hydrogen Bonds and the Importance of C-S···N Chalcogen Bonds in N'-(Adamantan-2-ylidene)hydrazide Derivatives

The present article comprehensively examines six N′-(adamantan-2-ylidene)hydrazide derivatives using the Hirshfeld surface analysis, PIXEL energy for molecular dimers, lattice energies for crystal packing, and topological analysis for intramolecular and intermolecular interactions. The crystal structure of one of the N′-(adamantan-2-ylidene)hydrazide derivatives, namely, N′-(adamantan-2-ylidene)-5-bromothiophene-2-carbohydrazide 1, C₁₅H₁₇N₂OSBr, has been determined and analyzed in detail along with five closely related structures. The molecular conformation of 1 is locked by an intramolecular C–S···N chalcogen bond as found in one of its closely related structure, namely, N′-(adamantan-2-ylidene)thiophene-2-carbohydrazide. Furthermore, a detailed potential energy surface scan analysis has been performed to highlight the importance of a chalcogen bond. Two of these compounds possess syn-orientation for amide units, whereas the corresponding moiety exhibits anti-conformations in the remaining four structures. The Hirshfeld surface and its decomposed fingerprint plots provide a qualitative picture of acyl substituent effects on the intermolecular interactions toward crystal packing of these six structures. Intermolecular interaction energies for dimers observed in these structures calculated by density functional theory (B97D3/def2-TZVP) and PIXEL (MP2/6-31G**) methods are comparable. This study also identifies that multiple hydrogen bonds, including N/C–H···O/N and C–H···π interactions, are collectively responsible for a self-assembled synthon. The nature and strength of these interactions have been studied using atoms in molecule topological analysis. The in vitro antiproliferative activity of compound 1 was assessed against five human tumor cell lines and showed marked antiproliferative activity.

Phosphine oxides as NMR and IR spectroscopic probes for geometry and energy of PO···H–A hydrogen bonds

In this work we evaluate the possibility to use the NMR and IR spectral properties of P=O group to estimate the geometry and strength of hydrogen bonds which it forms...

O-H stretching frequency red shifts do not correlate with the dissociation energies in the dimethylether and dimethylsulfide complexes of phenol derivatives

In this perspective, we present a comprehensive report on the spectroscopic and computational investigations of the hydrogen bonded (H-bonded) complexes of Me2O and Me2S with seven para-substituted H-bond donor phenols. The salient finding was that although the dissociation energies, D0, of the Me2O complexes were consistently higher than those of the analogous Me2S complexes, the red-shifts in phenolic O-H frequencies, Δν(O-H), showed the exactly opposite trend. This is in contravention of the general perception that the red shift in the X-H stretching frequency in the X-HY hydrogen bonded complexes is a reliable indicator of H-bond strength (D0), a concept popularly known as the Badger-Bauer rule. This is also in contrast to the trend reported for the H-bonded complexes of H2S/H2O with several para substituted phenols of different pKa values wherein the oxygen centered hydrogen bonded (OCHB) complexes consistently showed higher Δν(O-H) and D0 compared to those of the analogous sulfur centered hydrogen bonded (SCHB) complexes. Our effort was to understand these intriguing observations based on the spectroscopic investigations of 1 : 1 complexes in combination with a variety of high level quantum chemical calculations. Ab initio calculations at the MP2 level and the DFT calculations using various dispersion corrected density functionals (including DFT-D3) were performed on counterpoise corrected surfaces to compute the dissociation energy, D0, of the H-bonded complexes. The importance of anharmonic frequency computations is underscored as they were able to correctly reproduce the observed trend in the relative OH frequency shifts unlike the harmonic frequency computations. We have attempted to find a unified correlation that would globally fit the observed red shifts in the O-H frequency with the H-bonding strength for the four bases, namely, H2S, H2O, Me2O, and Me2S, in this set of H-bond donors. It was found that the proton affinity normalized Δν(O-H) values scale very well with the H-bond strength.

Molecular Interactions of Pyrazine-Based Compounds to Proteins

Pyrazine-based compounds are of great importance in medicinal chemistry. Due to their heteroaromatic nature, they uniquely combine properties of heteroatoms (polar interactions) with the properties of aromatic moieties (nonpolar interactions). This review summarizes results of a systematic analysis of RCSB PDB database focused on important binding interactions of pyrazine-based ligands cocrystallized in protein targets. The most frequent interaction of pyrazine was hydrogen bond to pyrazine nitrogen atom as an acceptor, followed by weak hydrogen bond with pyrazine hydrogen as donor. We also identified intramolecular hydrogen bonds within pyrazine ligands, π-interactions, coordination to metal ions, and few halogen bonds in chloropyrazines. In many cases the binding mode of the pyrazine fragment was complex, involving a combination of several interactions. We conclude that pyrazine as a molecular fragment should not be perceived as a simple aromatic isostere but rather as a readily interacting moiety of drug-like molecules with high potential for interactions to proteins.

The CH-π Interaction in Protein-Carbohydrate Binding: Bioinformatics and In Vitro Quantification

The molecular recognition of carbohydrates by proteins plays a key role in many biological processes including immune response, pathogen entry into a cell, and cell-cell adhesion (e.g., in cancer metastasis). Carbohydrates interact with proteins mainly through hydrogen bonding, metal-ion-mediated interaction, and non-polar dispersion interactions. The role of dispersion-driven CH-π interactions (stacking) in protein-carbohydrate recognition has been underestimated for a long time considering the polar interactions to be the main forces for saccharide interactions. However, over the last few years it turns out that non-polar interactions are equally important. In this study, we analyzed the CH-π interactions employing bioinformatics (data mining, structural analysis), several experimental (isothermal titration calorimetry (ITC), X-ray crystallography), and computational techniques. The Protein Data Bank (PDB) has been used as a source of structural data. The PDB contains over 12 000 protein complexes with carbohydrates. Stacking interactions are very frequently present in such complexes (about 39 % of identified structures). The calculations and the ITC measurement results suggest that the CH-π stacking contribution to the overall binding energy ranges from 4 up to 8 kcal mol^-1 . All the results show that the stacking CH-π interactions in protein-carbohydrate complexes can be considered to be a driving force of the binding in such complexes.

Cooperativity in alkane-1,2- and 1,3-polyols: NMR, QTAIM, and IQA study of O─H… OH and C─H… OH bonding interactions

Proton nuclear magnetic resonance chemical shifts and atom-atom interaction energies for alkanepolyols with 1,2-diol and 1,3-diol repeat units, and for their 1:1 pyridine complexes, are computed by density functional theory calculations. In the 1,3-polyols, based on a tG'Gg' repeat unit, the only important intramolecular hydrogen bonding interactions are O─H^… OH. By quantum theory of atoms in molecules analysis of the electron density, unstable bond and ring critical points are found for such interactions in 1,2-polyols with tG'g repeat units, from butane-1,2,3,4-tetrol onwards and in their pyridine complexes from propane-1,2,3-triol onwards. Several features (OH proton shifts and charges, and interaction energies computed by the interacting quantum atoms approach) are used to monitor the dependence of cooperativity on chain length: This is much less regular in 1,2-polyols than in 1,3-polyols and by most criteria has a higher damping factor. Well defined C─H^… OH interactions are found in butane-1,2,3,4-tetrol and higher members of the 1,2-polyol series, as well as in their pyridine complexes: There is no evidence for cooperativity with O─H^… OH bonding. For the 1,2-polyols, there is a tenuous empirical relationship between the existence of a bond critical point for O─H^… OH hydrogen bonding and the interaction energies of competing exchange channels, but the primary/secondary ratio is always less than unity.

Strategies in Translating the Therapeutic Potentials of Host Defense Peptides

The golden era of antibiotics, heralded by the discovery of penicillin, has long been challenged by the emergence of antimicrobial resistance (AMR). Host defense peptides (HDPs), previously known as antimicrobial peptides, are emerging as a group of promising antimicrobial candidates for combatting AMR due to their rapid and unique antimicrobial action. Decades of research have advanced our understanding of the relationship between the physicochemical properties of HDPs and their underlying antimicrobial and non-antimicrobial functions, including immunomodulatory, anti-biofilm, and wound healing properties. However, the mission of translating novel HDP-derived molecules from bench to bedside has yet to be fully accomplished, primarily attributed to their intricate structure-activity relationship, toxicity, instability in host and microbial environment, lack of correlation between in vitro and in vivo efficacies, and dwindling interest from large pharmaceutical companies. Based on our previous experience and the expanding knowledge gleaned from the literature, this review aims to summarize the novel strategies that have been employed to enhance the antimicrobial efficacy, proteolytic stability, and cell selectivity, which are all crucial factors for bench-to-bedside translation of HDP-based treatment. Strategies such as residues substitution with natural and/or unnatural amino acids, hybridization, L-to-D heterochiral isomerization, C- and N-terminal modification, cyclization, incorporation with nanoparticles, and "smart design" using artificial intelligence technology, will be discussed. We also provide an overview of HDP-based treatment that are currently in the development pipeline.

Unusual behaviour of the spin-spin coupling constant 1JCH upon formation of CHX hydrogen bond

One-bond coupling constants ¹J_XY are usually used as a measure of the corresponding XY interatomic distances. However, the physical nature of this correlation is not well understood and, in some cases, a counterintuitive behaviour of ¹J_XY upon hydrogen bonded complex formation has been reported. In this work, the behavior of ¹J_CH upon formation and strengthening of complexes with CHX hydrogen bonds and upon a proton transfer process is investigated by means of ¹H NMR spectroscopy and quantum chemical calculations. ¹H NMR spectra of 1,1-dinitroethane solution at room temperature in various solvents (carbon tetrachloride, chloroform, dichloromethane, acetone, dimethylformamide and dimethyl sulfoxide) illustrate the increase of ¹J_CH by several Hz upon an increase of the complex strength. Computational results (MP2/aug-cc-pVDZ) reproduce this observation and allow one to conclude that the increase of ¹J_CH is mainly caused by the change of the carbon hybridization (an increase of s-character), rather than by the change in interatomic distance r_CH. The behavior of ¹J_CH was also examined computationally for a wide range of CHX hydrogen bond energies and geometries. For this purpose, quantum-chemical modeling of the partial proton transfer process for complexes formed by 1,1-dinitroethane and trinitromethane as hydrogen bond donors with acetone, pyridine and fluoride anion as hydrogen bond acceptors was performed. The obtained results have confirmed the above-mentioned idea - for rather weak complexes, the dominant impact on the change of ¹J_CH magnitude is the increase of the s-character of carbon atom hybridization, while for complexes with a significantly transferred proton, the exponential decrease of the Fermi-contact term dominates.

Cooperative nature of the sulfur centered hydrogen bond: investigation of (H2S)n (n = 2-4) clusters using an affordable yet accurate level of theory

Existing studies have shown that appreciably high level quantum chemical calculations are required to reproduce experimental energetic and geometric features of a H₂S dimer. This condition severely restricts any practical possibility of obtaining reliable results for higher order H₂S clusters. We have shown here that the binding energies calculated at the CCSD(T)/CBS level with counterpoise corrected geometries calculated at the MP2/aug-cc-pV(Q+d)Z level of theory excellently match with the experimental results for the H₂S dimer. Subsequently, the above mentioned levels of theory were used for trimers and tetramers. (H₂S)_n (n = 2-4) clusters were found to show cooperative strengthening of S-HS hydrogen bonds, which is clearly evident from the evolution of binding energies and hydrogen bond lengths, with increasing cluster size. Localized molecular orbital energy decomposition analyses have been carried out to understand how the contributions of various energy components modulate with the size of the clusters and what are their relative contributions towards the overall stabilization of the clusters. Natural bond orbital and atoms in molecules analyses were also carried out in order to look into the evolution of the electronic charge transfer and electron density topology with cluster size.

Sub-Atomic Resolution Crystal Structures Reveal Conserved Geometric Outliers at Functional Sites

Myelin protein 2 (P2) is a peripheral membrane protein of the vertebrate nervous system myelin sheath, having possible roles in both lipid transport and 3D molecular organization of the multilayered myelin membrane. We extended our earlier crystallographic studies on human P2 and refined its crystal structure at an ultrahigh resolution of 0.72 Å in perdeuterated form and 0.86 Å in hydrogenated form. Characteristic differences in C–H…O hydrogen bond patterns were observed between extended β strands, kinked or ending strands, and helices. Often, side-chain C–H groups engage in hydrogen bonding with backbone carbonyl moieties. The data highlight several amino acid residues with unconventional conformations, including both bent aromatic rings and twisted guanidinium groups on arginine side chains, as well as non-planar peptide bonds. In two locations, such non-ideal conformations cluster, providing proof of local functional strain. Other ultrahigh-resolution protein structures similarly contain chemical groups, which break planarity rules. For example, in Src homology 3 (SH3) domains, a conserved bent aromatic residue is observed near the ligand binding site. Fatty acid binding protein (FABP) 3, belonging to the same family as P2, has several side chains and peptide bonds bent exactly as those in P2. We provide a high-resolution snapshot on non-ideal conformations of amino acid residues under local strain, possibly relevant to biological function. Geometric outliers observed in ultrahigh-resolution protein structures are real and likely relevant for ligand binding and conformational changes. Furthermore, the deuteration of protein and/or solvent are promising variables in protein crystal optimization.

Secondary Forces in Protein Folding

A complete inventory of the forces governing protein folding is critical for productive protein modeling, including structure prediction and de novo design, as well as understanding protein misfolding diseases of clinical significance. The dominant contributors to protein folding include the hydrophobic effect and conventional hydrogen bonding, along with Coulombic and van der Waals interactions. Over the past few decades, important additional contributors have been identified, including C-H···O hydrogen bonding, n→π* interactions, C5 hydrogen bonding, chalcogen bonding, and interactions involving aromatic rings (cation-π, X-H···π, π-π, anion-π, and sulfur-arene). These secondary contributions fall into two general classes: (1) weak but abundant interactions of the protein main chain and (2) strong but less frequent interactions involving protein side chains. Though interactions with high individual energies play important roles in specifying nonlocal molecular contacts and ligand binding, we estimate that weak but abundant interactions are likely to make greater overall contributions to protein folding, particularly at the level of secondary structure. Further research is likely to illuminate additional roles of these noncanonical interactions and could also reveal contributions yet unknown.

Exploring the Role of Consecutive Addition of Nitrogen Atoms on Stability and Reactivity of Hydrogen-Bonded Azine-Water Complexes

The second-order Møller-Plesset perturbation theory (MP2) and density functional theory with dispersion function calculations have been applied to investigate the hydrogen-bonding interaction between azines and water. The study suggests that the ability of nitrogen present in azine to act as a hydrogen-bond acceptor decreases in the order of pyridine (PY) > diazine (DZ) > triazine (TZ) > tetrazine (TTZ) > pentazine (PZ) > hexazine (HZ). Natural bond orbital (NBO) analysis, atoms in molecules, symmetry-adapted perturbation theory (SAPT), and molecular electrostatic potential studies reflect the factors important for hydrogen-bond strength as well as for the structural, electronic, and vibrational changes occurring during complexation. NBO analysis reflects that upon gradual addition of nitrogen atoms, hyperconjugation leads to an increase in the population of antibonding O-H bond, thus causing elongation and weakening of O-H bond in complexes incorporating N···H-O_W interaction, whereas rehybridization leads to an increase in the s character of the carbon hybrid orbital in C-H bond, thus causing contraction and shortening of C-H bond in complexes having C-H···O_W interactions. From the topological analysis, an excellent linear correlation is found to exist between stabilization energy (ΔE _BSSE), electron density (ρ_c), and its Laplacian (∇²ρ_c) at the bond critical points.

Intramolecular O-H⋯O and C-H⋯O hydrogen bond cooperativity in D-glucopyranose and D-galactopyranose-A DFT/GIAO, QTAIM/IQA, and NCI approach

Density functional theory calculations are used to compute proton nuclear magnetic resonance (NMR) chemical shifts, interatomic distances, atom-atom interaction energies, and atomic charges for partial structures and conformers of α-D-glucopyranose, β-D-glucopyranose, and α-D-galactopyranose built up by introducing OH groups into 2-methyltetrahydropyran stepwisely. For the counterclockwise conformers, the most marked effects on the NMR shift and the charge on the OH1 proton are produced by OH2, those of OH3 and OH4 being somewhat smaller. This argues for a diminishing cooperative effect. The effect of OH6 depends on the configuration of the hydroxymethyl group and the position, axial or equatorial, of OH4, which controls hydrogen bonding in the 1,3-diol motif. Variations in the interaction energies reveal that a "new" hydrogen bond is sometimes formed at the expense of a preexisting one, probably due to geometrical constraints. Whereas previous work showed that complexing a conformer with pyridine affects only the nearest neighbour, successive OH groups increase the interaction energy of the N⋯H1 hydrogen bond and reduce its length. Analogous results are obtained for the clockwise conformers. The interaction energies for C-H⋯OH hydrogen bonding between axial CH protons and OH groups in certain conformers are much smaller than for O-H⋯OH bonds but they are largely covalent, whereas those of the latter are predominantly coulombic. These interactions are modified by complexation with pyridine in the same way as O-H⋯OH interactions: the computed NMR shifts of the CH protons increase, the atom-atom distances are shorter, and interaction energies are enhanced.

C-H···O Hydrogen Bonding. The Prototypical Methane-Formaldehyde System: A Critical Assessment

Distinguishing the functionality of C-H···O hydrogen bonds (HBs) remains challenging, because their properties are difficult to quantify reliably. Herein, we present a study of the model methane-formaldehyde complex (MFC). Six stationary points on the MFC potential energy surface (PES) were obtained at the CCSD(T)/ANO2 level. The CCSDT(Q)/CBS interaction energies of the conformers range from only -1.12 kcal mol^-1 to -0.33 kcal mol^-1, denoting a very flat PES. Notably, only the lowest energy stationary point (MFC1) corresponds to a genuine minimum, whereas all other stationary points-including the previously studied ideal case of a_e(C-H···O) = 180°-exhibit some degree of freedom that leads to MFC1. Despite the flat PES, we clearly see that the HB properties of MFC1 align with those of the prototypical water dimer O-H···O HB. Each HB property generally becomes less prominent in the higher-energy conformers. Only the MFC1 conformer prominently exhibits (1) elongated C-H donor bonds, (2) attractive C-H···O═C interactions, (3) n(O) → σ*(C-H) hyperconjugation, (4) critical points in the electron density from Bader's method and from the noncovalent interactions method, (5) positively charged donor hydrogen, and (6) downfield NMR chemical shifts and nonzero ²J(C_M-H_M···O_F) coupling constants. Based on this research, some issues merit further study. The flat PES hinders reliable determinations of the HB-induced shifts of the C-H stretches; a similarly difficult challenge is observed for the experiment. The role of charge transfer in HBs remains an intriguing open question, although our BLW and NBO computations suggest that it is relevant to the C-H···O HB geometries. These issues notwithstanding, the prominence of the HB properties in MFC1 serves as clear evidence that the MFC is predominantly bound by a C-H···O HB.

The effect of oxidation on the structure of G-actin and its binding ability with aroma compounds in carp grass skeletal muscle

To investigate the influence of oxidative modifications of G-actin on its binding ability with aroma compounds, the influence of H₂O₂ treatments on G-actin structure and the absorption for alcohols and aldehydes was investigated. Raman spectroscopy and scanning electron microscopy were used to evaluate structural changes of G-actin; GC-MS was used to analyze the binding with alcohols and aldehydes. Results showed that 0-5mM H₂O₂ enhanced the absorption of G-actin toward alcohols involved in the formation of hydrogen bonds by increasing α-helix and carbonyl values. 0-1mM H₂O₂ caused the release of aldehydes with decreased sulfhydryl sites. 1-20mM H₂O₂ increased the retention of aldehydes, due to the increased hydrophobic sites by G-actin rebuilding and aggregating. The aggregated G-actin favoured the hydrophobic interactions with aroma compounds, forming the protein-aroma compound complex, thus enhancing the resultant binding ability, as evidenced by scanning electron microscopy and GC/MS analysis.

An optimal hydrogen-bond surrogate for α-helices

Substitution of a main chain i → i + 4 hydrogen bond with a covalent bond can nucleate and stabilize the α-helical conformation in peptides. Herein we describe the potential of different alkene isosteres to mimic intramolecular hydrogen bonds and stabilize α-helices in diverse peptide sequences.

Helix-helix interactions in membrane domains of bitopic proteins: Specificity and role of lipid environment

Interaction between transmembrane helices often determines biological activity of membrane proteins. Bitopic proteins, a broad subclass of membrane proteins, form dimers containing two membrane-spanning helices. Some aspects of their structure-function relationship cannot be fully understood without considering the protein-lipid interaction, which can determine the protein conformational ensemble. Experimental and computer modeling data concerning transmembrane parts of bitopic proteins are reviewed in the present paper. They highlight the importance of lipid-protein interactions and resolve certain paradoxes in the behavior of such proteins. Besides, some properties of membrane organization provided a clue to understanding of allosteric interactions between distant parts of proteins. Interactions of these kinds appear to underlie a signaling mechanism, which could be widely employed in the functioning of many membrane proteins. Treatment of membrane proteins as parts of integrated fine-tuned proteolipid system promises new insights into biological function mechanisms and approaches to drug design. This article is part of a Special Issue entitled: Lipid order/lipid defects and lipid-control of protein activity edited by Dirk Schneider.

Data on diverse roles of helix perturbations in membrane proteins

The various structural variations observed in TM helices of membrane proteins have been deconstructed into 9 distinct types of helix perturbations. These perturbations are defined by the deviation of TM helices from the predominantly observed linear α-helical conformation, to form 3₁₀- and π-helices, as well as adopting curved and kinked geometries. The data presented here supplements the article ‘Helix perturbations in Membrane Proteins Assist in Inter-helical Interactions and Optimal Helix Positioning in the Bilayer’ (A. Shelar, M. Bansal, 2016) [1]. This data provides strong evidence for the role of various helix perturbations in influencing backbone torsion angles of helices, mediating inter-helical interactions, oligomer formation and accommodation of hydrophobic residues within the bilayer. The methodology used for creation of various datasets of membrane protein families (Sodium/Calcium exchanger and Heme Copper Oxidase) has also been mentioned.

SpeeDB: fast structural protein searches

MOTIVATION

Interactions between amino acids are important determinants of the structure, stability and function of proteins. Several tools have been developed for the identification and analysis of such interactions in proteins based on the extensive studies carried out on high-resolution structures from Protein Data Bank (PDB). Although these tools allow users to identify and analyze interactions, analysis can only be performed on one structure at a time. This makes it difficult and time consuming to study the significance of these interactions on a large scale.

RESULTS

SpeeDB is a web-based tool for the identification of protein structures based on structural properties. SpeeDB queries are executed on all structures in the PDB at once, quickly enough for interactive use. SpeeDB includes standard queries based on published criteria for identifying various structures: disulphide bonds, catalytic triads and aromatic-aromatic, sulphur-aromatic, cation-π and ionic interactions. Users can also construct custom queries in the user interface without any programming. Results can be downloaded in a Comma Separated Value (CSV) format for further analysis with other tools. Case studies presented in this article demonstrate how SpeeDB can be used to answer various biological questions. Analysis of human proteases revealed that disulphide bonds are the predominant type of interaction and are located close to the active site, where they promote substrate specificity. When comparing the two homologous G protein-coupled receptors and the two protein kinase paralogs analyzed, the differences in the types of interactions responsible for stability accounts for the differences in specificity and functionality of the structures.

AVAILABILITY AND IMPLEMENTATION

SpeeDB is available at http://www.parallelcomputing.ca as a web service.

CONTACT

d@drobilla.net

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

An approach to creating a more realistic working model from a protein data bank entry

An accurate model of three-dimensional protein structure is important in a variety of fields such as structure-based drug design and mechanistic studies of enzymatic reactions. While the entries in the Protein Data Bank ( http://www.pdb.org ) provide valuable information about protein structures, a small fraction of the PDB structures were found to contain anomalies not reported in the PDB file. The semiempirical PM7 method in MOPAC2012 was used for identifying anomalously short hydrogen bonds, C-H⋯O/C-H⋯N interactions, non-bonding close contacts, and unrealistic covalent bond lengths in recently published Protein Data Bank files. It was also used to generate new structures with these faults removed. When the semiempirical models were compared to those of PDB_REDO (http://www.cmbi.ru.nl/pdb_redo/), the clashscores, as defined by MolProbity ( http://molprobity.biochem.duke.edu/), were better in about 50% of the structures. The semiempirical models also had a lower root-mean-square-deviation value in nearly all cases than those from PDB_REDO, indicative of a better conservation of the tertiary structure. Finally, the semiempirical models were found to have lower clashscores than the initial PDB file in all but one case. Because this approach maintains as much of the original tertiary structure as possible while improving anomalous interactions, it should be useful to theoreticians, experimentalists, and crystallographers investigating the structure and function of proteins.

Mechanism of stabilization of helix secondary structure by constrained Cα-tetrasubstituted α-amino acids

The theoretical basis behind the ability of constrained Cα-tetrasubstituted amino acids (CTAAs) to induce stable helical conformations has been studied through Replica Exchange Molecular Dynamics Potential of Mean Force Quantum Theory of Atoms In Molecules calculations on Ac-l-Ala-CTAA-l-Ala-Aib-l-Ala-NHMe peptide models. We found that the origin of helix stabilization by CTAAs can be ascribed to at least two complementary mechanisms limiting the backbone conformational freedom: steric hindrance predominantly in the (+x,+y,-z) sector of a right-handed 3D Cartesian space, where the z axis coincides with the helical axis and the Cα of the CTAA lies on the +y axis (0,+y,0), and the establishment of additional and relatively strong C-H···O interactions involving the CTAA.

ω-Turn: a novel β-turn mimic in globular proteins stabilized by main-chain to side-chain C−H···O interaction

Mimicry of structural motifs is a common feature in proteins. The 10-membered hydrogen-bonded ring involving the main-chain C − O in a β-turn can be formed using a side-chain carbonyl group leading to Asx-turn. We show that the N − H component of hydrogen bond can be replaced by a C(γ) -H group in the side chain, culminating in a nonconventional C − H···O interaction. Because of its shape this β-turn mimic is designated as ω-turn, which is found to occur ∼ three times per 100 residues. Three residues (i to i + 2) constitute the turn with the C − H···O interaction occurring between the terminal residues, constraining the torsion angles ϕi + 1, ψi + 1, ϕi + 2 and χ'1(i + 2) (using the interacting C(γ) atom). Based on these angles there are two types of ω-turns, each of which can be further divided into two groups. C(β) -branched side-chains, and Met and Gln have high propensities to occur at i + 2; for the last two residues the carbonyl oxygen may participate in an additional interaction involving the S and amino group, respectively. With Cys occupying the i + 1 position, such turns are found in the metal-binding sites. N-linked glycosylation occurs at the consensus pattern Asn-Xaa-Ser/Thr; with Thr at i + 2, the sequence can adopt the secondary structure of a ω-turn, which may be the recognition site for protein modification. Location between two β-strands is the most common occurrence in protein tertiary structure, and being generally exposed ω-turn may constitute the antigenic determinant site. It is a stable scaffold and may be used in protein engineering and peptide design.

Strong ionic hydrogen bonding causes a spectral isotope effect in photoactive yellow protein

Standard hydrogen bonds are of great importance for protein structure and function. Ionic hydrogen bonds often are significantly stronger than standard hydrogen bonds and exhibit unique properties, but their role in proteins is not well understood. We report that hydrogen/deuterium exchange causes a redshift in the visible absorbance spectrum of photoactive yellow protein (PYP). We expand the range of interpretable isotope effects by assigning this spectral isotope effect (SIE) to a functionally important hydrogen bond at the active site of PYP. The inverted sign and extent of this SIE is explained by the ionic nature and strength of this hydrogen bond. These results show the relevance of ionic hydrogen bonding for protein active sites, and reveal that the inverted SIE is a novel, to our knowledge, tool to probe ionic hydrogen bonds. Our results support a classification of hydrogen bonds that distinguishes the properties of ionic hydrogen bonds from those of both standard and low barrier hydrogen bonds, and show how this classification helps resolve a recent debate regarding active site hydrogen bonding in PYP.

C–H/O interactions of nucleic bases with a water molecule: a crystallographic and quantum chemical study

The C–H/O interactions of nucleic bases are substantially stronger than the C–H/O interactions of benzene and pyridine. These results can be very important for molecular recognition of DNA and RNA.

First steps in growth of a polypeptide toward β-sheet structure

The full conformational energy surface is examined for a molecule in which a dipeptide is attached to the same spacer group as another peptide chain, so as to model the seminal steps of β-sheet formation. This surface is compared with the geometrical preferences of the isolated dipeptide to extract the perturbations induced by interactions with the second peptide strand. These interpeptide interactions remove any tendency of the dipeptide to form a C5 ring structure, one of its two normally stable geometries. A C7 structure, the preferred conformation of the isolated dipeptide, remains as the global minimum in the full molecule. However, the stability of this structure is highly dependent upon interpeptide H-bonds with the second chain. The latter forces include not only the usual NH···O interaction, but also a pair of CH···O H-bonds. The secondary minimum is also of C7 type and likewise depends in part upon CH···O H-bonds for its stability. The latter interactions also play a part in the tertiary minimum. A two-strand β-sheet structure is not yet in evidence for this small model system, requiring additional peptide units to be added to each chain.

Stacking interactions between carbohydrate and protein quantified by combination of theoretical and experimental methods

Carbohydrate – receptor interactions are an integral part of biological events. They play an important role in many cellular processes, such as cell-cell adhesion, cell differentiation and in-cell signaling. Carbohydrates can interact with a receptor by using several types of intermolecular interactions. One of the most important is the interaction of a carbohydrate's apolar part with aromatic amino acid residues, known as dispersion interaction or CH/π interaction. In the study presented here, we attempted for the first time to quantify how the CH/π interaction contributes to a more general carbohydrate - protein interaction. We used a combined experimental approach, creating single and double point mutants with high level computational methods, and applied both to Ralstonia solanacearum (RSL) lectin complexes with α-l-Me-fucoside. Experimentally measured binding affinities were compared with computed carbohydrate-aromatic amino acid residue interaction energies. Experimental binding affinities for the RSL wild type, phenylalanine and alanine mutants were −8.5, −7.1 and −4.1 kcal.mol⁻¹, respectively. These affinities agree with the computed dispersion interaction energy between carbohydrate and aromatic amino acid residues for RSL wild type and phenylalanine, with values −8.8, −7.9 kcal.mol⁻¹, excluding the alanine mutant where the interaction energy was −0.9 kcal.mol⁻¹. Molecular dynamics simulations show that discrepancy can be caused by creation of a new hydrogen bond between the α-l-Me-fucoside and RSL. Observed results suggest that in this and similar cases the carbohydrate-receptor interaction can be driven mainly by a dispersion interaction.