The occurrence of C-H...O hydrogen bonds in proteins

Ultrapotent influenza hemagglutinin fusion inhibitors developed through SuFEx-enabled high-throughput medicinal chemistry

Seasonal and pandemic-associated influenza strains cause highly contagious viral respiratory infections that can lead to severe illness and excess mortality. Here, we report on the optimization of our small-molecule inhibitor F0045(S) targeting the influenza hemagglutinin (HA) stem with our Sulfur-Fluoride Exchange (SuFEx) click chemistry-based high-throughput medicinal chemistry (HTMC) strategy. A combination of SuFEx- and amide-based lead molecule diversification and structure-guided design led to identification and validation of ultrapotent influenza fusion inhibitors with subnanomolar EC50 cellular antiviral activity against several influenza A group 1 strains. X-ray structures of six of these compounds with HA indicate that the appended moieties occupy additional pockets on the HA surface and increase the binding interaction, where the accumulation of several polar interactions also contributes to the improved affinity. The compounds here represent the most potent HA small-molecule inhibitors to date. Our divergent HTMC platform is therefore a powerful, rapid, and cost-effective approach to develop bioactive chemical probes and drug-like candidates against viral targets.

Structure, Self-Assembly, and Phase Behavior of Neuroactive N-Acyl GABAs: Doxorubicin Encapsulation in NPGABA/DPPC Liposomes and Release Studies

N-Acylated amino acids and neurotransmitters in mammals exert significant biological effects on the nervous system, immune responses, and vasculature. N-Acyl derivatives of γ-aminobutyric acid (N-acyl GABA), which belong to both classes mentioned above, are prominent among them. In this work, a homologous series of N-acyl GABAs bearing saturated N-acyl chains (C8-C18) have been synthesized and characterized with respect to self-assembly, thermotropic phase behavior, and supramolecular organization. Differential scanning calorimetric studies revealed that the transition enthalpies and entropies of N-acyl GABAs are linearly dependent on the acyl chain length. The crystal structure of N-tridecanoyl GABA showed that the molecules are packed in bilayers with the acyl chains aligned parallel to the bilayer normal and that the carboxyl groups from opposite layers associate to form dimeric structures involving strong O-H···O hydrogen bonds. In addition, N-H···O and C-H···O hydrogen bonds between amide moieties of adjacent molecules within each layer stabilize the molecular packing. Powder X-ray diffraction studies showed odd-even alternation in the d spacings, suggesting that the odd chain and even chain compounds pack differently. Equimolar mixtures of N-palmitoyl GABA and dipalmitoylphosphatidylcholine (DPPC) were found to form stable unilamellar vesicles with diameters of ∼300-340 nm, which could encapsulate doxorubicin, an anticancer drug, with higher efficiency and better release characteristics than DPPC liposomes at physiologically relevant pH. These liposomes exhibit faster release of doxorubicin at acidic pH (<7.0), indicating their potential utility as drug carriers in cancer chemotherapy.

De novo antioxidant peptide design via machine learning and DFT studies

Antioxidant peptides (AOPs) are highly valued in food and pharmaceutical industries due to their significant role in human function. This study introduces a novel approach to identifying robust AOPs using a deep generative model based on sequence representation. Through filtration with a deep-learning classification model and subsequent clustering via the Butina cluster algorithm, twelve peptides (GP1-GP12) with potential antioxidant capacity were predicted. Density functional theory (DFT) calculations guided the selection of six peptides for synthesis and biological experiments. Molecular orbital representations revealed that the HOMO for these peptides is primarily localized on the indole segment, underscoring its pivotal role in antioxidant activity. All six synthesized peptides exhibited antioxidant activity in the DPPH assay, while the hydroxyl radical test showed suboptimal results. A hemolysis assay confirmed the non-hemolytic nature of the generated peptides. Additionally, an in silico investigation explored the potential inhibitory interaction between the peptides and the Keap1 protein. Analysis revealed that ligands GP3, GP4, and GP12 induced significant structural changes in proteins, affecting their stability and flexibility. These findings highlight the capability of machine learning approaches in generating novel antioxidant peptides.

Effects of filler type and content on mechanical, thermal, and physical properties of carrageenan biocomposite films

This study investigates the influence of various fillers on the properties of carrageenan, a natural polymer derived from red seaweed. Despite its potential for enhanced biocomposite film development, carrageenan faces challenges related to strength. The incorporation of nanoclay into the carrageenan film resulted in a significant increase in film thickness from 0.026 to 0.068 mm. The UV light transmission value for the carrageenan film alone was measured at 30.9 %, whereas films containing 5 wt% of Tetraethyl orthosilicate (TEOS), 3-Aminopropyltriethoxysilane (APTES), and nanoclay exhibited reduced transmission values of 23 %, 18 %, and 1 %, respectively. Notably, the tensile strength of the unfilled carrageenan film was 38.4 MPa, which increased to 38.6, 57, and 60 MPa upon the addition of 3 wt% of nanoclay, APTES, and TEOS fillers, respectively. All fillers contributed to improved tensile strength, with TEOS demonstrating the highest enhancement. The optimal filler content was determined to be 3 wt%. Regarding thermal properties, films containing TEOS displayed higher thermal stability compared to those with APTES, while films incorporating nanoclay exhibited the lowest stability. Findings provide insights into the effects of different fillers on the mechanical, physical and thermal properties of carrageenan films, supporting the development of improved biocomposite materials suitable for application in food packaging.

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.

AQcalc: A web server that identifies weak molecular interactions in protein structures

Weak molecular interactions play an important role in protein structure and function. Computational tools that identify weak molecular interactions are, therefore, valuable for the study of proteins. Here, we present AQcalc, a web server (https://aqcalcbiocomputing.com/) that can be used to identify anion-quadrupole (AQ) interactions, which are weak interactions involving aromatic residue (Trp, Tyr, and Phe) ring edges and anions (Asp, Glu, and phosphate ion) both within proteins and at their interfaces (protein-protein, protein-nucleic acids, and protein-lipid bilayer). AQcalc identifies AQ interactions as well as clusters involving AQ, cation-π, and salt bridges, among others. Utilizing AQcalc we analyzed weak interactions in protein models, even in the absence of experimental structures, to understand the contributions of weak interactions to deleterious structural changes, including those associated with oncogenic and germline disease variants. We identified several deleterious variants with disrupted AQ interactions (comparable in frequency to cation-π disruptions). Amyloid fibrils utilize AQ to bury anions at frequencies that far exceed those observed for globular proteins. AQ interactions were detected three and five times more frequently than the hydrogen-bonded AQ (HBAQ) in fibril structures and protein-lipid bilayer interfaces, respectively. By contrast, AQ and HBAQ interactions were detected with similar frequencies in globular proteins. Collectively, these findings suggest AQcalc will be effective in facilitating fine structural analysis. As other web utilities designed to identify protein residue interaction networks do not report AQ interactions, wide use of AQcalc will enrich our understanding of residue interaction networks and facilitate hypothesis testing by identifying and experimentally characterizing these comparably weak but important interactions.

Porous and Close Packed Supramolecular Assemblies from 2,4-Difluoronitrobenzene with Three Different Linkers and an n-Butylamine Cap

A porous structure formed from sheets with cavities and two close packed structures were crystallised from building blocks prepared from 2,4-difluoronitrobenzene, a diamine linker and n-butylamine. The porous structure crystallised from a flexible building block prepared using 1,4-diaminobutane as linker. The close packed structures were prepared using either piperazine or 1,4-bis(aminomethyl)benzene as a linker and have less conformational freedom.

AlphaFold2 models of the active form of all 437 catalytically competent human protein kinase domains

Humans have 437 catalytically competent protein kinase domains with the typical kinase fold, similar to the structure of Protein Kinase A (PKA). Only 155 of these kinases are in the Protein Data Bank in their active form. The active form of a kinase must satisfy requirements for binding ATP, magnesium, and substrate. From structural bioinformatics analysis of 40 unique substrate-bound kinases, we derived several criteria for the active form of protein kinases. We include requirements on the DFG motif of the activation loop but also on the positions of the N-terminal and C-terminal segments of the activation loop that must be placed appropriately to bind substrate. Because the active form of catalytic kinases is needed for understanding substrate specificity and the effects of mutations on catalytic activity in cancer and other diseases, we used AlphaFold2 to produce models of all 437 human protein kinases in the active form. This was accomplished with templates in the active form from the PDB and shallow multiple sequence alignments of orthologs and close homologs of the query protein. We selected models for each kinase based on the pLDDT scores of the activation loop residues, demonstrating that the highest scoring models have the lowest or close to the lowest RMSD to 22 non-redundant substrate-bound structures in the PDB. A larger benchmark of all 130 active kinase structures with complete activation loops in the PDB shows that 80% of the highest-scoring AlphaFold2 models have RMSD < 1.0 Å and 90% have RMSD < 2.0 Å over the activation loop backbone atoms. Models for all 437 catalytic kinases are available at http://dunbrack.fccc.edu/kincore/activemodels. We believe they may be useful for interpreting mutations leading to constitutive catalytic activity in cancer as well as for templates for modeling substrate and inhibitor binding for molecules which bind to the active state.

Human Endocrine-Disrupting Effects of Phthalate Esters through Adverse Outcome Pathways: A Comprehensive Mechanism Analysis

Phthalate esters (PAEs) are widely exposed in the environment as plasticizers in plastics, and they have been found to cause significant environmental and health hazards, especially in terms of endocrine disruption in humans. In order to investigate the processes underlying the endocrine disruption effects of PAEs, three machine learning techniques were used in this study to build an adverse outcome pathway (AOP) for those effects on people. According to the results of the three machine learning techniques, the random forest and XGBoost models performed well in terms of prediction. Subsequently, sensitivity analysis was conducted to identify the initial events, key events, and key features influencing the endocrine disruption effects of PAEs on humans. Key features, such as Mol.Wt, Q+, QH+, ELUMO, minHCsats, MEDC-33, and EG, were found to be closely related to the molecular structure. Therefore, a 3D-QSAR model for PAEs was constructed, and, based on the three-dimensional potential energy surface information, it was discovered that the hydrophobic, steric, and electrostatic fields of PAEs significantly influence their endocrine disruption effects on humans. Lastly, an analysis of the contributions of amino acid residues and binding energy (BE) was performed, identifying and confirming that hydrogen bonding, hydrophobic interactions, and van der Waals forces are important factors affecting the AOP of PAEs' molecular endocrine disruption effects. This study defined and constructed a comprehensive AOP for the endocrine disruption effects of PAEs on humans and developed a method based on theoretical simulation to characterize the AOP, providing theoretical guidance for studying the mechanisms of toxicity caused by other pollutants.

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.

Kinetic versus thermodynamic polymorph stabilization of a tri-carboxylic acid derivative at the solid-liquid interface

The carboxylic acid moiety gives rise to structural variability in surface-supported self-assembly due to the common expression of various H-bonding motifs. Self-assembly of 3-fold symmetric tricarboxylic acid derivatives on surfaces typically results in monolayer structures that feature the common 2-fold cyclic R22(8) H-bond motif for at least one of the carboxylic acid groups. Polymorphs that are exclusively based on 3-fold cyclic R33(12) H-bonds were predicted but remained elusive. Here, we show the emergence of such a superflower (SF) structure purely based on R33(12) H-bonds for L-benzene-1,3,5-tricarbonyl phenylalanine (L-BTA), a molecule derived from the well-studied trimesic acid (TMA). In contrast to TMA, L-BTA is not completely planar and is also equipped with additional functional groups for the formation of secondary intermolecular bonds. At the heptanoic acid-graphite interface we transiently observe a SF structure, which is dynamically converted into a chicken-wire structure that only exhibits R22(8) H-bonds. Interestingly, when using nonanoic acid as a solvent the initially formed SF structure remained stable. This unexpected behaviour is rationalized by accompanying force field simulations and experimental determination of solvent-dependent L-BTA solubility.

Norethindrone suppress the germ cell development via androgen receptor resulting in male bias

Progestins are widely used and detected in surface waters, and can affect gonad development and sexual differentiation in fish. However, the toxicological mechanisms of sexual differentiation induced by progestins are not well understood. Here, we investigated the effects of norethindrone (NET) and androgen receptor (AR) antagonist flutamide (FLU) on gonadal differentiation in zebrafish from 21 dpf (days post-fertilization) to 49 dpf. The results showed that NET caused male bias, while FLU resulted in female bias at 49 dpf. The NET and FLU mixtures significantly decreased the percentage of males compared to the NET single exposure. Molecular docking analysis showed that FLU and NET had similar docking pocket and docking posture with AR resulting in competitively forming the hydrogen bond with Thr334 of AR. These results suggested that binding to AR was the molecular initiating event of sex differentiation induced by NET. Moreover, NET strongly decreased transcription of biomarker genes (dnd1, ddx4, dazl, piwil1 and nanos1) involved in germ cell development, while FLU significantly increased transcription of these target genes. There was an increase in the number of juvenile oocytes, which was consistent with the female bias in the combined groups. The bliss independence model analysis further showed that NET and FLU had antagonistic effect on transcription and histology during gonadal differentiation. Thus, NET suppressed the germ cell development via AR, resulting in male bias. Understanding the molecular initiation of sex differentiation in progestins is essential to provide a comprehensive biological basis for ecological risk assessment.

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.

Vibrational Solvatochromism Study of the C-H···O Improper Hydrogen Bond

The improper C-H···O hydrogen bond is an important weak interaction, with broad implications for protein and nucleic acid structure, molecular recognition, enzyme catalysis, and drug interaction. Despite its wide identification in crystal structures, the general existence of C-H···O hydrogen bonds remains elusive especially for natural C-H groups in bulk aqueous solutions at room temperature. Vibrational spectroscopy is a promising methodology to tackle this challenge, as formation of C-H···O hydrogen bonds usually causes shifts of the C-H stretch frequency. Yet, prior observations are inconclusive, as they are all based on a simple blue-shift in aqueous solution and cannot distinguish if it is an effect caused by solvent reorganization or a specific hydrogen-bonding interaction. In this work, we used vibrational solvatochromism as a calibration of the solvent reorganization effect and identified a specific H-bonding interaction. We performed vibrational solvatochromism study of C-H(D) of multiple alcohol molecules including the CH mode of CD₃CH(OH)CD₃ and the CD₃ modes of CD₃OH, CD₃CH₂OH, and CD₃CH(OH)CD₃ in a series of solvents. We found an abnormal blue-shift of the Raman frequency of the C-H and C-D bonds at both the C_α and C_β positions of alcohols in water, which lies in an opposite direction to the expected trend due to vibrational solvatochromism. This experimental evidence supports that the improper C-H···O hydrogen bonds might generally exist between nonpolarized C-H and water in liquid solutions at room temperature.

Non-covalent interactions reveal the protein chain δ conformation in a flexible single-residue model

The δ conformation is a local secondary structure in proteins that implicates a π_amide N-H⋯N interaction between a backbone N atom and the NH of the following residue. Small-molecule models thereof have been limited so far to rigid proline-type compounds. We show here that in derivatives of a cyclic amino acid with a sulphur atom in the γ-position, specific side-chain/backbone N-H⋯S interactions stabilize the δ conformation sufficiently to allow it to compete with classical C5 and C7 H-bonded conformers.

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.

(-)-Epicatechin Is a Biased Ligand of Apelin Receptor

(-)-Epicatechin (EC) is part of a large family of biomolecules called flavonoids and is widely distributed in the plant kingdom. Several studies have shown the beneficial effects of EC consumption. Many of these reported effects are exerted by activating the signaling pathways associated with the activation of two specific receptors: the G protein-coupled estrogen receptor (GPER), a transmembrane receptor, and the pregnane X receptor (PXR), which is a nuclear receptor. However, the effects of EC are so diverse that these two receptors cannot describe the complete phenomenon. The apelin receptor or APLNR is classified within the G protein-coupled receptor (GPCR) family, and is capable of activating the G protein canonical pathways and the β-arrestin transducer, which participates in the phenomenon of receptor desensitization and internalization. β-arrestin gained interest in selective pharmacology and mediators of the so-called "biased agonism". With molecular dynamics (MD) and in vitro assays, we demonstrate how EC can recruit the β-arrestin in the active conformation of the APLN receptor acting as a biased agonist.

Dydrogesterone and levonorgestrel at environmentally relevant concentrations have antagonist effects with rhythmic oscillation in brain and eyes of zebrafish

Synthetic progestins levonorgestrel (LNG) and dydrogesterone (DDG) are frequency detected in surface water. Combined effects of LNG and DDG on gonad differentiation are similar to LNG single exposure in juvenile zebrafish. However, LNG and DDG mixtures have stronger effects on spermatogenesis in testes of adult zebrafish, which show variable at different life stage. Effects of LNG and DDG mixtures on eyes and brain remain unknown. Here we investigated effects of LNG, DDG and their mixtures on eyes and brain. Zebrafish were exposed to LNG, DDG and their mixtures from 2 hpf to 144 dpf. Rhythm and vision related biological processes were enriched in eyes and brain in LNG and DDG treatments, which indicated rhythmic oscillation in eyes and brain. The qPCR data revealed that both LNG and DDG decreased transcription of arntl2 and clocka, while increased transcription of per1a, per1b, rpe65a and tefa in eyes and brain. However, DDG and LNG mixtures had slight effect on transcription of genes related to rhythm and vision. In addition, LNG and DDG reduced the thickness of inner nuclear layer in the eyes. Bliss independent model revealed that LNG and DDG had antagonist effects on transcription and histology in eyes and brain. Moreover, LNG and DDG formed the same hydrogen bonds with green-sensitive opsin-4 and rhodopsin kinase GRK7a. Taken together, LNG and DDG competed with each other for the same binding residues resulting in antagonist effect in their mixtures treatments, and have significant ecological implications to assess combined effects of progestins mixtures on fish in different organs.

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.

Deciphering C–H⋯O/X weak hydrogen bonding and halogen bonding interactions in aromatic peptoids

We deciphered weak interactions in aromatic peptoids, such as C–H⋯O/X, and simultaneously identified strong interactions, including N–H⋯N and N–H⋯O, in this class of foldamer.

Newly identified C–H⋯O hydrogen bond in histidine

Histidine C–H bonds observed in protein structures include (clockwise from top left): myoglobin, β-lactamase, and photoactive yellow protein; calculations indicate that tautomeric/protonation state influences H-bonding ability (bottom left).

Bispidine as a β-strand nucleator: from a β-arch to self-assembled cages and vesicles

The development of synthetic scaffolds that nucleate well-folded secondary structures is highly challenging. Herein, we designed and synthesized a series of core-modified peptides (F1, F2, F3, and F4) that fold into β-strand structures. These bispidine-scaffolded peptides were studied by CD, IR, NMR, single crystal XRD, and Molecular Dynamics (MD) simulations to investigate their conformational preferences. Solid-state and solution studies revealed that bispidine is a versatile scaffold that could be placed either at the terminal or at the middle of the peptide strand for nucleating the β-strand structure. Scaffolds that nucleate an isolated β-strand conformation are rare. Bispidine placed at the C-terminus of the peptide chain could nucleate a β-strand conformation, while bispidine placed at the middle resulted in a β-arch conformation. This nucleation activity stems from the ability to restrict the psi torsion angle (ψ) through intramolecular C5 hydrogen bonding between the equatorial hydrogen(s) of bispidine and the carbonyl oxygen(s) of the amino acid close to the scaffold. Furthermore, the bispidine peptidomimetic with a super secondary structure, namely β-arch, assembled into single-hole submicron cages and spherical vesicles as evident from microscopic studies. The design logic defined here will be a significant strategy for the development of β-strand mimetics and super secondary structures.

Characterization of a triazole scaffold compound as an inhibitor of Mycobacterium tuberculosis imidazoleglycerol-phosphate dehydratase

Mycobacterium tuberculosis (Mtb), the causative agent of human tuberculosis (TB), employs ten enzymes including imidazoleglycerol-phosphate dehydratase (IGPD) for de novo biosynthesis of histidine. The absence of histidine-biosynthesis in humans combined with its essentiality for Mtb makes the enzymes of this pathway major anti-TB drug targets. We explored the inhibitory potential of a small molecule β-(1,2,4-Triazole-3-yl)-DL-alanine (DLA) against Mtb IGPD. DLA exhibits an in vitro inhibitory efficacy in the lower micromolar range. Higher-resolution crystal structures of native and substrate-bound Mtb IGPD provided additional structural features of this important drug target. Crystal structure of IGPD-DLA complex at a resolution of 1.75 Å, confirmed that DLA locks down the function of the enzyme by binding in the active site pocket of the IGPD mimicking the substrate-binding mode to a high degree. In our biochemical study, DLA showed an efficient inhibition of Mtb IGPD. Furthermore, DLA also showed bactericidal activity against Mtb and Mycobacterium smegmatis and inhibited their growth in respective culture medium. Importantly, owing to the favorable ADME and physicochemical properties, it serves as an important lead molecule for further derivatizations.

CH⋯O H-bond mediated tautomerization of 2-methyl-1,3-cyclohexanedione: A combined IR spectroscopic and theoretical study

Molecular association and its impact on the keto-enol tautomerization of 2-methyl-1,3-cyclohexanedione (MCHD) have been investigated in low temperature argon matrix and thin solid film. The system exists exclusively in diketo tautomeric form in argon matrix. The CH⋯O H-bonded homodimers of the diketo tautomer are produced by annealing the matrix at 28 K. No trace of the keto-enol tautomer is observed in matrix isolated homodimers in the temperature range of 8-28 K. However, tautomeric conversion initiates in a thin film of pure diketo tautomer when the temperature of the film is raised to ~170 K. Transition state calculations on the monomeric and dimeric MCHD demonstrate that CH⋯O H-bond formations between diketo tautomers play a vital role in lowering the tautomerization barrier. However, the extent of CH⋯O H-bonded dimer formation in matrix isolation, as well as extent of tautomerization in the neat sample are found to be smaller than that for the previously reported 1,3-cyclohexanedione (CHD) under similar experimental conditions (J. Phys. Chem. A 2012, 116, 3836-3845). Electronic structure calculations suggest that formation of the CH⋯O H-bonded dimer is less feasible in presence of the bulky 2-methyl groups of MCHD, as compared to CHD. Additionally, the transition state geometry of the dimeric keto-enol product of MCHD, as compared to the same for CHD, is more strained and offers a weaker CH---O H-bond that contributes to lesser tautomeric conversion in the former.

In silico approach of modified melanoma peptides and their immunotherapeutic potential

Melanoma is a type of skin cancer with increasing incidence worldwide and high lethality. Conventional forms of treatment are not effective in advanced cancer stages. Hence, immunotherapeutic approaches have been tested to modulate immune response against tumor cells. Some vaccine models using tumor-associated antigens (TAAs) such as glycoprotein 100 (gp100) have been studied, but their expected effectiveness has not been shown until now. Antigen immunogenicity is a crucial point to improve the immune response, and therefore mutations are inserted in peptide sequences. It is possible to understand the interactions which occur between peptides and immune system molecules through computer simulation, and this is essential in order to guide efficient vaccine models. In this work, we have calculated the interaction binding energies of crystallographic data based on modified gp100 peptides and HLA-A*0201 using density functional theory (DFT) and the molecular fractionation with conjugated caps (MFCC) approach. Our results show the most relevant residue-residue interactions, the impact of three mutations in their binding sites, and the main HLA-A*0201 amino acids for peptide-HLA binding.

When are two hydrogen bonds better than one? Accurate first-principles models explain the balance of hydrogen bond donors and acceptors found in proteins

Hydrogen bonds (HBs) play an essential role in the structure and catalytic action of enzymes, but a complete understanding of HBs in proteins challenges the resolution of modern structural (i.e., X-ray diffraction) techniques and mandates computationally demanding electronic structure methods from correlated wavefunction theory for predictive accuracy. Numerous amino acid sidechains contain functional groups (e.g., hydroxyls in Ser/Thr or Tyr and amides in Asn/Gln) that can act as either HB acceptors or donors (HBA/HBD) and even form simultaneous, ambifunctional HB interactions. To understand the relative energetic benefit of each interaction, we characterize the potential energy surfaces of representative model systems with accurate coupled cluster theory calculations. To reveal the relationship of these energetics to the balance of these interactions in proteins, we curate a set of 4000 HBs, of which >500 are ambifunctional HBs, in high-resolution protein structures. We show that our model systems accurately predict the favored HB structural properties. Differences are apparent in HBA/HBD preference for aromatic Tyr versus aliphatic Ser/Thr hydroxyls because Tyr forms significantly stronger O–H⋯O HBs than N–H⋯O HBs in contrast to comparable strengths of the two for Ser/Thr. Despite this residue-specific distinction, all models of residue pairs indicate an energetic benefit for simultaneous HBA and HBD interactions in an ambifunctional HB. Although the stabilization is less than the additive maximum due both to geometric constraints and many-body electronic effects, a wide range of ambifunctional HB geometries are more favorable than any single HB interaction.

Correlated wavefunction theory predicts and high-resolution crystal structure analysis confirms the important, stabilizing effect of simultaneous hydrogen bond donor and acceptor interactions in proteins.

A Geometric Definition of Short to Medium Range Hydrogen-Mediated Interactions in Proteins

We present a method to rapidly identify hydrogen-mediated interactions in proteins (e.g., hydrogen bonds, hydrogen bonds, water-mediated hydrogen bonds, salt bridges, and aromatic π-hydrogen interactions) through heavy atom geometry alone, that is, without needing to explicitly determine hydrogen atom positions using either experimental or theoretical methods. By including specific real (or virtual) partner atoms as defined by the atom type of both the donor and acceptor heavy atoms, a set of unique angles can be rapidly calculated. By comparing the distance between the donor and the acceptor and these unique angles to the statistical preferences observed in the Protein Data Bank (PDB), we were able to identify a set of conserved geometries (15 for donor atoms and 7 for acceptor atoms) for hydrogen-mediated interactions in proteins. This set of identified interactions includes every polar atom type present in the Protein Data Bank except OE1 (glutamate/glutamine sidechain) and a clear geometric preference for the methionine sulfur atom (SD) to act as a hydrogen bond acceptor. This method could be readily applied to protein design efforts.

Pinacolone-Alcohol Gas-Phase Solvation Balances as Experimental Dispersion Benchmarks

The influence of distant London dispersion forces on the docking preference of alcohols of different size between the two lone electron pairs of the carbonyl group in pinacolone was explored by infrared spectroscopy of the OH stretching fundamental in supersonic jet expansions of 1:1 solvate complexes. Experimentally, no pronounced tendency of the alcohol to switch from the methyl to the bulkier tert-butyl side with increasing size was found. In all cases, methyl docking dominates by at least a factor of two, whereas DFT-optimized structures suggest a very close balance for the larger alcohols, once corrected by CCSD(T) relative electronic energies. Together with inconsistencies when switching from a C4 to a C5 alcohol, this points at deficiencies of the investigated B3LYP and in particular TPSS functionals even after dispersion correction, which cannot be blamed on zero point energy effects. The search for density functionals which describe the harmonic frequency shift, the structural change and the energy difference between the docking isomers of larger alcohols to unsymmetric ketones in a satisfactory way is open.

Papain-like cysteine proteinase zone (PCP-zone) and PCP structural catalytic core (PCP-SCC) of enzymes with cysteine proteinase fold

There are several families of cysteine proteinases with different folds – for example the (chymo)trypsin fold family and papain-like fold family – but in both families the hydrolase activity of cysteine proteinases requires a cysteine residue as the catalytic nucleophile. In this work, we have analyzed the topology of the active site regions in 146 three-dimensional structures of proteins belonging to the Papain-like Cysteine Proteinase (PCP) superfamily, which includes papain as a typical representative of this protein superfamily. All analyzed enzymes contain a unique structurally closed conformation – a “PCP-Zone” – which can be divided into two groups, Class A and Class B. Eight structurally conserved amino acids of the PCP-Zone form a common Structural Core. The Structural Core, catalytic nucleophile, catalytic base and residue Xaa – which stabilizes the side-chain conformation of the catalytic base – make up a PCP Structural Catalytic Core (PCP-SCC). The PCP-SCC of Class A and Class B are divided into 5 and 2 types, respectively. Seven variants of the mutual arrangement of the amino-acid side chains of the catalytic triad – nucleophile, base and residue Xaa – within the same fold clearly demonstrate how enzymes with the papain-like fold adapt to the need to perform diverse functions in spite of their limited structural diversity. The roles of both the PCP-Zone of SARS-CoV-2-PLpro described in this study and the NBCZone of SARS-CoV-2-3CLpro presented in our earlier article (Denesyuk AI, Johnson MS, Salo-Ahen OMH, Uversky VN, Denessiouk K. Int J Biol Macromol. 2020;153:399-411) that are in contacts with inhibitors are discussed.

Computational and In Vitro Investigation of (-)-Epicatechin and Proanthocyanidin B2 as Inhibitors of Human Matrix Metalloproteinase 1

Matrix metalloproteinases 1 (MMP-1) energetically triggers the enzymatic proteolysis of extracellular matrix collagenase (ECM), resulting in progressive skin aging. Natural flavonoids are well known for their antioxidant properties and have been evaluated for inhibition of matrix metalloproteins in human. Recently, (-)-epicatechin and proanthocyanidin B2 were reported as essential flavanols from various natural reservoirs as potential anti-inflammatory and free radical scavengers. However, their molecular interactions and inhibitory potential against MMP-1 are not yet well studied. In this study, sequential absorption, distribution, metabolism, and excretion (ADME) profiling, quantum mechanics calculations, and molecular docking simulations by extra precision Glide protocol predicted the drug-likeness of (-)-epicatechin (−7.862 kcal/mol) and proanthocyanidin B2 (−8.145 kcal/mol) with the least reactivity and substantial binding affinity in the catalytic pocket of human MMP-1 by comparison to reference bioactive compound epigallocatechin gallate (−6.488 kcal/mol). These flavanols in docked complexes with MMP-1 were further studied by 500 ns molecular dynamics simulations that revealed substantial stability and intermolecular interactions, viz. hydrogen and ionic interactions, with essential residues, i.e., His218, Glu219, His222, and His228, in the active pocket of MMP-1. In addition, binding free energy calculations using the Molecular Mechanics Generalized Born Surface Area (MM/GBSA) method suggested the significant role of Coulomb interactions and van der Waals forces in the stability of respective docked MMP-1-flavonol complexes by comparison to MMP-1-epigallocatechin gallate; these observations were further supported by MMP-1 inhibition assay using zymography. Altogether with computational and MMP-1–zymography results, our findings support (-)-epicatechin as a comparatively strong inhibitor of human MMP-1 with considerable drug-likeness against proanthocyanidin B2 in reference to epigallocatechin gallate.

Proving and Probing the Presence of the Elusive C-H⋅⋅⋅O Hydrogen Bond in Liquid Solutions at Room Temperature

Hydrogen bonds (H bonds) play a major role in defining the structure and properties of many substances, as well as phenomena and processes. Traditional H bonds are ubiquitous in nature, yet the demonstration of weak H bonds that occur between a highly polarized C-H group and an electron-rich oxygen atom, has proven elusive. Detailed here are linear and nonlinear IR spectroscopy experiments that reveal the presence of H bonds between the chloroform C-H group and an amide carbonyl oxygen atom in solution at room temperature. Evidence is provided for an amide solvation shell featuring two clearly distinguishable chloroform arrangements that undergo chemical exchange with a time scale of about 2 ps. Furthermore, the enthalpy of breaking the hydrogen bond is found to be 6-20 kJ mol^-1 . Ab-initio computations support the findings of two distinct solvation shells formed by three chloroform molecules, where one thermally undergoes hydrogen-bond making and breaking.

Conformation control through concurrent N-H⋯S and N-H⋯O[double bond, length as m-dash]C hydrogen bonding and hyperconjugation effects

In addition to the classical N-H⋯O[double bond, length as m-dash]C non-covalent interaction, less conventional types of hydrogen bonding, such as N-H⋯S, may play a key role in determining the molecular structure. In this work, using theoretical calculations in combination with spectroscopic analysis in both gas phase and solution phase, we demonstrate that both these H-bonding modes exist simultaneously in low-energy conformers of capped derivatives of Attc, a thietane α-amino acid. 6-Membered ring inter-residue N-H⋯S interactions (C6^γ), assisted by hyperconjugation between the thietane ring and the backbone, combine with 5-membered ring intra-residue backbone N-H⋯O[double bond, length as m-dash]C interactions (C5) to provide a C5-C6^γ feature that stabilizes a planar geometry in the monomer unit. Two contiguous C5-C6^γ features in the planar dimer implicate an unprecedented three-centre H-bond of the type C[double bond, length as m-dash]O⋯H(N)⋯SR₂, while the trimer adopts two C5-C6^γ features separated by a Ramachandran α-type backbone configuration. These low-energy conformers are fully characterized in the gas phase and support is presented for their existence in solution state.

Increasing protein stability by engineering the n → π* interaction at the β-turn

Abundant n → π* interactions between adjacent backbone carbonyl groups, identified by statistical analysis of protein structures, are predicted to play an important role in dictating the structure of proteins. However, experimentally testing the prediction in proteins has been challenging due to the weak nature of this interaction. By amplifying the strength of the n → π* interaction via amino acid substitution and thioamide incorporation at a solvent exposed β-turn within the GB1 proteins and Pin 1 WW domain, we demonstrate that an n → π* interaction increases the structural stability of proteins by restricting the ϕ torsion angle. Our results also suggest that amino acid side-chain identity and its rotameric conformation play an important and decisive role in dictating the strength of an n → π* interaction.

Amino acid residues adopt a right-handed α-helical conformation with increasing strength of the n → π* interaction. We also demonstrate a direct consequence of n → π* interactions on enhancing the structural stability of proteins.

Tuning Supramolecular Selectivity for Hydrosulfide: Linear Free Energy Relationships Reveal Preferential C-H Hydrogen Bond Interactions

Supramolecular anion receptors can be used to study the molecular recognition properties of the reactive yet biologically critical hydrochalcogenide anions (HCh^-). Achieving selectivity for HCh^- over the halides is challenging but necessary for not only developing future supramolecular probes for HCh^- binding and detection, but also for understanding the fundamental properties that govern these binding and recognition events. Here we demonstrate that linear free energy relationships (LFERs)-including Hammett and Swain-Lupton plots-reveal a clear difference in sensitivity to the polarity of an aryl C-H hydrogen bond (HB) donor for HS^- over other HCh^- and halides. Analysis using electrostatic potential maps highlights that this difference in sensitivity results from a preference of the aryl C-H HB donor for HS^- in this host scaffold. From this study, we demonstrate that LFERs are a powerful tool to gain interpretative insight into motif design for future anion-selective supramolecular receptors and highlight the importance of C-H HB donors for HS^- recognition. From our results, we suggest that aryl C-H HB donors should be investigated in the next generation of HS^- selective receptors based on the enhanced HS^- selectivity over other competing anions in this system.

C─H⋯O hydrogen bonds in kinase-inhibitor interfaces

C─H⋯O hydrogen bonds constitute a unique class of cohesive interactions. Their properties are similar to those of canonical H-bonds, although their energy is significantly lower, typically in the 0.5-2.5 kcal/mol range. Polarised C─H groups, such as those adjacent to electronegative groups, or within aromatic moieties, are particularly strong donors. C─H⋯O bonds are ubiquitous in nucleic acids and in proteins, notably stabilizing the β-sheet secondary structure. They have also been observed in numerous protein-ligand interactions. Here, we analysed crystal structures, deposited in the Protein Data Bank, of complexes of FDA-approved protein kinase inhibitors with cognate kinases, to assess the possible role of C─H^inhibitor ⋯O^protein hydrogen bonds. The conserved hinge motif of protein kinases with two solvent-exposed carbonyl groups and one exposed backbone amide, is well known to be involved in canonical H-bonding with inhibitors. We now find that in virtually all complexes where the inhibitor interacts with the hinge backbone, at least one of the hinge carbonyl groups accepts an H-bond from a C─H inhibitor group, which is either aromatic or adjacent to an electronegative group. These observations are important for design of hinge-binding scaffolds of novel kinase inhibitors for therapeutic use.

The acid-base-nucleophile catalytic triad in ABH-fold enzymes is coordinated by a set of structural elements

The alpha/beta-Hydrolases (ABH) are a structural class of proteins that are found widespread in nature and includes enzymes that can catalyze various reactions in different substrates. The catalytic versatility of the ABH fold enzymes, which has been a valuable property in protein engineering applications, is based on a similar acid-base-nucleophile catalytic mechanism. In our research, we are concerned with the structure that surrounds the key units of the catalytic machinery, and we have previously found conserved structural organizations that coordinate the catalytic acid, the catalytic nucleophile and the residues of the oxyanion hole. Here, we explore the architecture that surrounds the catalytic histidine at the active sites of enzymes from 40 ABH fold families, where we have identified six conserved interactions that coordinate the catalytic histidine next to the catalytic acid and the catalytic nucleophile. Specifically, the catalytic nucleophile is coordinated next to the catalytic histidine by two weak hydrogen bonds, while the catalytic acid is directly involved in the coordination of the catalytic histidine through by two weak hydrogen bonds. The imidazole ring of the catalytic histidine is coordinated by a CH-π contact and a hydrophobic interaction. Moreover, the catalytic triad residues are connected with a residue that is located at the core of the active site of ABH fold, which is suggested to be the fourth member of a “structural catalytic tetrad”. Besides their role in the stability of the catalytic mechanism, the conserved elements of the catalytic site are actively involved in ligand binding and affect other properties of the catalytic activity, such as substrate specificity, enantioselectivity, pH optimum and thermostability of ABH fold enzymes. These properties are regularly targeted in protein engineering applications, and thus, the identified conserved structural elements can serve as potential modification sites in order to develop ABH fold enzymes with altered activities.

Per-Residue Program of Multiple Backbone Dihedral Angles of β-Peptoids via Backbone Substitutions

Unique folded structures of natural and synthetic oligomers are the most fundamental basis for their unique functions. N-Substituted β-peptides, or β-peptoids, are synthetic oligomers with great potential to fold into diverse three-dimensional structures because of the existence of four rotatable bonds in a monomer with highly modular synthetic accessibility. However, the existence of the four rotatable bonds poses a challenge for conformational control of β-peptoids. Here, we report a strategy for per-residue programming of two dihedral angles of β-peptoids, which is useful for restricting the conformational space of the oligomers. The oligomer was found to form a unique loop conformation that is stabilized by the backbone rotational restrictions. Circular dichroism and NMR spectroscopic analyses and X-ray crystallographic analysis of the oligomer are presented. The strategy would significantly facilitate the discovery of many more unique folded structures of β-peptoids.

Sagrillo F, R. Freitas MC, A. Alves M, Thota S, Tinoco LW, Magalhães A, Sant’Anna CMR, Fraga CAM. Structure–property relationship studies of 3-acyl-substituted furans: the serendipitous identification and characterization of a new non-classical hydrogen bond donor moiety

A serendipitous identification and characterization of a new non-classical hydrogen bond donor moiety found in N-acylhydrazones containing 3-acyl-substituted furan subunit is presented.

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.

Ribosomal RNA-Aminoglycoside Hygromycin B Interaction Energy Calculation within a Density Functional Theory Framework

We intend to investigate the drug-binding energy of each nucleotide inside the aminoglycoside hygromycin B (hygB) binding site of 30S ribosomal RNA (rRNA) subunit by using the molecular fractionation with conjugate caps (MFCC) strategy based on the density functional theory (DFT), considering the functional LDA/PWC, OBS, and the dielectric constant parametrization. Aminoglycosides are bactericidal antibiotics that have high affinity to the prokaryotic rRNA, inhibiting the synthesis of proteins by acting on the main stages of the translation mechanism, whereas binding to rRNA 16S, a component of the 30S ribosomal subunit in prokaryotes. The identification of the nucleotides presenting the most negative binding energies allows us to stabilize hygB in a suitable binding pocket of the 30S ribosomal subunit. In addition, it should be highlighted that mutations in these residues may probably lead to resistance to ribosome-targeting antibiotics. Quantum calculations of aminoglycoside hygromycin B-ribosome complex might contribute to further quantum studies with antibiotics like macrolides and other aminoglycosides.