"Strong" hydrogen bonds in chemistry and biology

Direct evidence for a deprotonated lysine serving as a H-bond "acceptor" in a photoreceptor protein

Deprotonation or suppression of the pKa of the amino group of a lysine sidechain is a widely recognized phenomenon whereby the sidechain amino group transiently can act as a nucleophile at the active site of enzymatic reactions. However, a deprotonated lysine and its molecular interactions have not been directly experimentally detected. Here, we demonstrate a deprotonated lysine stably serving as an "acceptor" in a H-bond between the photosensor protein RcaE and its chromophore. Signal splitting and trans-H-bond J coupling observed by NMR spectroscopy provide direct evidence that Lys261 is deprotonated and serves as a H-bond acceptor for the chromophore NH group. Quantum mechanical/molecular mechanical calculations also indicate that this H-bond exists stably. Interestingly, the sidechain amino group of the lysine can act as both donor and acceptor. The remarkable shift in the H-bond characteristics arises from a decrease in solvation, triggered by photoisomerization. Our results provide insights into the dual role of this lysine. This mechanism has broad implications for other biological reactions in which lysine plays a role.

From pixels to druggable leads: A CADD strategy for the design and synthesis of potent DDR1 inhibitors

BACKGROUND AND OBJECTIVE

While numerous in silico tools exist for target-based drug discovery, the inconsistent integration of in vitro data with predictive models hinders research and development productivity. This is particularly apparent during the Hit-to-Lead stage, where unreliable in-silico tools often lead to suboptimal lead selection. Herein, we address this challenge by presenting a CADD-guided pipeline that successfully integrates rational drug design with in-silico hits to identify a promising DDR1 lead.

METHODS

2 × 1000 ns MD simulations along with their respective FEL and MMPBSA analyses were employed to guide the rational design and synthesis of 12 novel compounds which were evaluated for their DDR inhibition.

RESULTS

The molecular dynamics investigation of the initial hit led to the identification of key structural features within the DDR1 binding pocket. The identified key features were used to guide the rational design and synthesis of twelve novel derivatives. SAR analysis, biological evaluation, molecular dynamics, and free energy calculations were carried out for the synthesized derivatives to understand their mechanism of action. Compound 4c exhibited the strongest inhibition and selectivity for DDR1, with an IC50 of 0.11 µM.

CONCLUSIONS

The MD simulations led to the identification of a key hydrophobic groove in the DDR1 binding pocket. The integrated approach of SAR analysis with molecular dynamics led to the identification of compound 4c as a promising lead for further development of potent and selective DDR1 inhibitors. Moreover, this work establishes a protocol for translating in silico hits to real world bioactive druggable leads.

Unveiling the Reaction Mystery Between Lithium Polysulfides and Lithium Metal Anode in Lithium-Sulfur Batteries

Lithium-sulfur (Li-S) batteries are widely regarded as one of the most promising next-generation high-energy-density energy storage devices. However, soluble lithium polysulfides (LiPSs) corrode Li metal and deteriorate the cycling stability of Li-S batteries. Understanding the reaction mechanism between LiPSs and Li metal anode is imperative. Herein, the reaction rate and products of LiPSs with Li metal anode, the composition and structure of the as-generated solid electrolyte interphase (SEI), and the mechanism of lithium nitrate (LiNO3) additives for inhibiting the corrosion reactions are systematically unveiled. Concretely, LiPSs react with Li metal anode more rapidly than Li salt and generate a Li2S-rich SEI. The Li2S-rich SEI is highly reactive with LiPSs, which exacerbates the formation of dendritic Li and the continuous corrosion of active Li. LiNO3 functions dominantly by modulating the solvation structure of LiPSs and inherently reducing the reactivity of LiPSs, rather than the conventional understanding of LiNO3 participating in the formation of SEI. This work reveals the reaction mechanism between LiPSs and Li metal anode and inspires rational regulating of the solvation structure of LiPSs for stabilizing Li metal anode in Li-S batteries.

A low-barrier proton shared between two aspartates acts as a conformational switch that changes the substrate specificity of the β-lactamase BlaC

Serine β-lactamases inactivate β-lactam antibiotics in a two-step mechanism comprising acylation and deacylation. For the deacylation step, a water molecule is activated by a conserved glutamate residue to release the adduct from the enzyme. The third-generation cephalosporin ceftazidime is a poor substrate for the class A β-lactamase BlaC from Mycobacterium tuberculosis but it can be hydrolyzed faster when the active site pocket is enlarged, as was reported for mutant BlaC P167S. The conformational change in the Ω-loop of the P167S mutant displaces the conserved glutamate (Glu166), suggesting it is not required for deacylation of the ceftazidime adduct. Here, we report the characterization of wild type BlaC and BlaC E166A at various pH values. The presence of Glu166 strongly enhances activity against nitrocefin but not ceftazidime, indicating it is indeed not required for deacylation of the adduct of the latter substrate. At high pH wild type BlaC was found to exist in two states, one of which converts ceftazidime much faster, resembling the open state previously reported for the BlaC mutant P167S. The pH-dependent switch between the closed and open states is caused by the loss at high pH of a low-barrier hydrogen bond, a proton shared between Asp172 and Asp179. These results illustrate how readily shifts in substrate specificity can occur as a consequence of subtle changes in protein structure.

Novel Inhibitory Site Revealed by Mode of Action on the Metabotropic Glutamate 7 Receptor Venus Flytrap Domain

Metabotropic glutamate (mGlu) receptors play a key role in modulating most synapses in the brain. The mGlu7 receptors inhibit presynaptic neurotransmitter release and offer therapeutic possibilities for post-traumatic stress disorders or epilepsy. Screening campaigns provided mGlu7-specific allosteric modulators as the inhibitor XAP044 (Gee et al. J. Biol. Chem. 2014). In contrast to other mGlu receptor allosteric modulators, XAP044 does not bind in the transmembrane domain but to the extracellular domain of the mGlu7 receptor and not at the orthosteric site. Here, we identified the mode of action of XAP044, combining synthesis of derivatives, modeling and docking experiments, and mutagenesis. We propose a unique mode of action of these inhibitors, preventing the closure of the Venus flytrap agonist binding domain. While acting as a noncompetitive antagonist of L-AP4, XAP044 and derivatives act as apparent competitive antagonists of LSP4-2022. These data revealed more potent XAP044 analogues and new possibilities to target mGluRs.

Melanin: Nature's 4th bioorganic polymer

The pigments known as the melanins are widely recognized for their responsibility in the coloration of human skin, eyes, hair, and minimising the harmful effects of solar ultraviolet radiation. But specialists are aware that the melanins are present in all living kingdoms, barring viruses, and have functionality that extends beyond neutralizing ionising radiation. The ubiquitous presence of melanin in almost all human organs, recognized in recent years, as well as the presence of melanin in organisms that are evolutionarily distant from each other, indicate the fundamental importance of this class of material for all life forms. In this review, we argue for the need to accept melanins as the fourth primordial class of biological polymers, along with nucleic acids, proteins and polysaccharides. We consistently compare the properties of these canonical biological polymers with the properties of melanin and highlight key features that fundamentally distinguish melanins, their function and its mysteries.

Computational study of the dimerization of glyphosate: mechanism and effect of solvent

A computational study on the structure and stability of different series of glyphosate (Glyph) dimers comprising nonionized (N) and zwitterionic structures (Z) for neutral monomers, followed by an analysis of energetics of Glyph dimerization process have been performed by means of quantum chemical calculations in different media. Optimized geometries for energy minima, as well as relative potential and free energies of the possible various conformers of each series of Glyph dimer were computed as a function of the medium at B3LYP-D3/6-311++G(2d,2p) level. The solvation model based on density (SMD) is employed for all solution phase computations. Non-ionized dimers (DN), anion-cation (AC) and either zwitterion-zwitterion (DZP and DZC) or non-ionized-phosphonate zwitterion (NZP) ionized neutral forms of Glyph dimer are predicted to exist in the gas phase and in solution in large contrast to Glyph monomers. The DZC dimer form exhibiting a centrosymmetric arrangement of two carboxylate zwitterion units was found to be the most stable dimer structure in all media. In aqueous solution, the DZP and AC dimer type structures are significantly stabilized by hydration. The tautomerisms between DZC, DZP and AC dimer type structures have been investigated in the gas phase and in solution. The DZC type structures are more prone to experience proton transfer in water than in the gas phase and in cyclohexane. The mechanism for the tautomerization process in neutral ionized Glyph dimers proceeds via two direct proton transfer paths: DZP ⇋ AC ⇋ DZC. Results show that solvents play a key role in modulating the energetics of the dimerization process of Glyph. Solvation in cyclohexane, favors the dimerization process however, hydration opposed it. In aqueous solution, the mechanism of the dimerization of Glyph from its phosphonate zwitterionic monomer form (ZP) could be described by a set of equilibria including direct proton transfer paths as follows: 2ZP ⇋ DZP ⇋ AC ⇋ DZC. According to our results, in aqueous solution, DZC Glyph dimers and their corresponding DZP and AC tautomers should be present in higher concentration than phosphonate zwitterionic Glyph monomers for high Glyph concentration, a fact that seems controversial in the literature.

Phosphorylation of collagen fibrils enhances intrafibrillar mineralization and dentin remineralization

The hierarchical assembly of nanoapatite within a type I collagen matrix was achieved through biomimetic mineralization in vitro, cooperatively regulated by non-collagenous proteins and small biomolecules. Here, we demonstrated that IP6 could significantly promote intrafibrillar mineralization in two- and three-dimensional collagen models through binding to collagen fibrils via hydrogen bonds (the interaction energy ∼10.21 kJ mol-1), as confirmed by the FTIR spectra and isothermal experimental results. In addition, we find that IP6 associated with dental collagen fibrils can also enhance the remineralization of calcium-depleted dentin and restore its mechanical properties similar to the natural dentin within 4 days. The promoting effect is mainly due to the chemical modification of IP6, which alters the interfacial physicochemical properties of collagen fibrils, strengthening the interaction of calcium phosphate minerals and mineral ions with collagen fibrils. This strategy of interfacial regulation to accelerate the mineralization of collagen fibrils is essential for dental repair and the development of a clinical product for the remineralization of hard tissue.

Unique hydrogen-bonding network in a viral channelrhodopsin

Channelrhodopsins (CRs) are used as key tools in optogenetics, and novel CRs, either found from nature or engineered by mutation, have greatly contributed to the development of optogenetics. Recently CRs were discovered from viruses, and crystal structure of a viral CR, OLPVR1, reported a very similar water-containing hydrogen-bonding network near the retinal Schiff base to that of a light-driven proton-pump bacteriorhodopsin (BR). In both OLPVR1 and BR, nearly planar pentagonal cluster structures are comprised of five oxygen atoms, three oxygens from water molecules and two oxygens from the Schiff base counterions. The planar pentagonal cluster stabilizes a quadrupole, two positive charges at the Schiff base and an arginine, and two negative charges at the counterions, and thus plays important roles in light-gated channel function of OLPVR1 and light-driven proton pump function of BR. Despite similar pentagonal cluster structures, present FTIR analysis revealed different hydrogen-bonding networks between OLPVR1 and BR. The hydrogen bond between the protonated Schiff base and a water is stronger in OLPVR1 than in BR, and internal water molecules donate hydrogen bonds much weaker in OLPVR1 than in BR. In OLPVR1, the bridged water molecule between the Schiff base and counterions forms hydrogen bonds to D76 and D200 equally, while the hydrogen-bonding interaction is much stronger to D85 than to D212 in BR. The present interpretation is supported by the mutation results, where D76 and D200 equally work as the Schiff base counterions in OLPVR1, but D85 is the primary counterion in BR. This work reports highly sensitive hydrogen-bonding network in the Schiff base region, which would be closely related to each function through light-induced alterations of the network.

Energetics of the H-Bond Network in Exiguobacterium sibiricum Rhodopsin

Exiguobacterium sibiricum rhodopsin (ESR) functions as a light-driven proton pump utilizing Lys96 for proton uptake and maintaining its activity over a wide pH range. Using a combination of methodologies including the linear Poisson-Boltzmann equation and a quantum mechanical/molecular mechanical approach with a polarizable continuum model, we explore the microscopic mechanisms underlying its pumping activity. Lys96, the primary proton uptake site, remains deprotonated owing to the loss of solvation in the ESR protein environment. Asp85, serving as a proton acceptor group for Lys96, does not form a low-barrier H-bond with His57. Instead, deprotonated Asp85 forms a salt-bridge with protonated His57, and the proton is predominantly located at the His57 moiety. Glu214, the only acidic residue at the end of the H-bond network exhibits a pKa value of ∼6, slightly elevated due to solvation loss. It seems likely that the H-bond network [Asp85···His57···H2O···Glu214] serves as a proton-conducting pathway toward the protein bulk surface.

Study on 2D Molecular Networks of Flexible Pentacarboxylic Acid Ligands Induced by Ether Bonds in Response to Selective Guest Inclusion

The flexibility of ligands allows for their bending, twisting, or rotation to adopt various conformations, leading to distinct symmetries during the self-assembled process. Flexible aromatic acid ligands modified by ether bonds are a promising type of self-assembled module when it comes to surfaces. Here, two pentacarboxylic acid ligands (H5L1 and H5L2) with minor skeleton differences have successfully self-assembled into disparate porous networks on the graphite surface and demonstrated excellent potential for the inclusion of guest molecules. The H5L1 molecule's network structure only accommodates coronene (COR) molecules. With fewer COR molecules, H5L1 molecules act as a host template to accommodate the COR molecules. When there are too many COR molecules, COR molecules will induce H5L1 molecules to transform into a new host-guest nanostructure. Additionally, H5L2 molecules showed the ability to capture C70 molecules and exhibited cavity selectivity. However, the assembled network of H5L2 was slightly deformed in attempts to trap the COR molecules. To understand these phenomena more deeply, various assembled mechanisms were analyzed in combination with building theoretical models and energy analysis. These results reveal the great potential of flexible aromatic acid ligands in two-dimensional self-assembly and host-guest systems for their application in related fields.

Redox-Mediated Amination of Pyrogallol-Based Polyphenols

Flavonoids are known to covalently modify amyloidogenic peptides by amination reactions. The underlying coupling process between polyphenols and N-nucleophiles is assessed by several in vitro and in silico approaches. The coupling reaction involves a sequence of oxidative dearomatization, amination, and reductive amination (ODARA) reaction steps. The C6-regioselectivity of the product is confirmed by crystallographic analysis. Under aqueous conditions, the reaction of baicalein with lysine derivatives yields C-N coupling as well as hydrolysis products of transient imine intermediates. The observed C-N coupling reactions work best for flavonoids combining a pyrogallol substructure with an electron-withdrawing group attached to the C4a-position. Thermodynamic properties such as bond dissociation energies also highlight the key role of pyrogallol units for the antioxidant ability. Combining the computed electronic properties and in vitro antioxidant assays suggests that the studied pyrogallol-containing flavonoids act by various radical-scavenging mechanisms working in synergy. Multivariate analysis indicates that a small number of descriptors for transient intermediates of the ODARA process generates a model with excellent performance (r=0.93) for the prediction of cross-coupling yields. The same model has been employed to predict novel antioxidant flavonoid-based molecules as potential covalent inhibitors, opening a new avenue to the design of therapeutically relevant anti-amyloid compounds.

Design of Highly Efficient Electronic Energy Transfer in Functionalized Quantum Dots Driven Specifically by Ethylenediamine

The exploration of emerging functionalized quantum dots (QDs) through modulating the effective interaction between the sensing element and target analyte is of great significance for high-performance trace sensing. Here, the chromone-based ligand grafted QDs (QDs-Chromone) were initiated to realize the electronic energy transfer (EET) driven specifically by ethylenediamine (EDA) in the absence of spectral overlap. The fluorescent and colorimetric dual-mode responses (from red to blue and from colorless to yellow, respectively) resulting from the expanded conjugated ligands reinforced the analytical selectivity, endowing an ultrasensitive and specific response to submicromolar-liquid of EDA. In addition, a QDs-Chromone-based sensing chip was constructed to achieve the ultrasensitive recognition of EDA vapor with a naked-eye observed response at a concentration as low as 10 ppm, as well as a robust anti-interfering ability in complicated scenarios monitoring. We expect the proposed EET strategy in shaping functionalized QDs for high-performance sensing will shine light on both rational probe design methodology and deep sensing mechanism exploration.

Deliver on a Promise: Hydrogen-Bonded Polymer Nanomedicine with a Precise Ratio of Chemodrug and Photosensitizer for Intelligent Cancer Therapy

The outcomes of combined cancer therapy are largely related to loading content and contribution of each therapeutic agent; however, fine-tuning the ratio of two coloaded components toward precise cancer therapy is a great challenge and still remains in its infancy. We herein develop a supramolecular polymer scaffold to optimize the coloading ratio of chemotherapeutic agent and photosensitizer through hydrogen-bonding (H-bonding) interaction, for maximizing the efficacy of intelligent cancer chemo/photodynamic therapies (CT/PDT). To do so, we first synthesize a thymine (THY)-functionalized tetraphenylporphyrin photosensitizer (i.e., TTPP), featuring the same molecular configuration of H-bonding array with chemotherapeutic carmofur (e.g., 1-hexylcarbamoyl-5-fluorouracil, HCFU). Meanwhile, a six-arm star-shaped amphiphilic polymer vehicle P(DAPA-co-DPMA-co-OEGMA)6 (poly(diaminopyridine acrylamide-co-2-(diisopropylamino)ethyl methacrylate-co-oligo(ethylene glycol) monomethyl ether methacrylate)6) is prepared, bearing hydrophilic and biocompatible POEGMA segment, along with hydrophobic PDAPA and PDPMA segments, characterizing the randomly dispersed dual functionalities, i.e., heterocomplementary H-bonding DAP motifs and pH-responsive protonation DPMA content. Thanks to the identical DAP/HCFU and DAP/TTPP H-bonding association capability, the incorporation of both HCFU and TTPP to six-arm star-shaped P(DAPA-co-DPMA-co-OEGMA)6 vehicle, with an optimized coloading ratio, can be straightforwardly realized by adjusting the feeding concentrations, thus yielding the hydrogen-bonded supramolecular nanoparticles (i.e., HCFU-TTPP-SPNs), demonstrating the codelivery of two components with the promise to optimize the combined CT/PDT efficacy.

Computational and crystallographic study of hydrogen bonds in the second coordination sphere of chelated amino acids with a free water molecule: Influence of complex charge and metal ion

Hydrogen bonds of glycine complexes were calculated using quantum chemistry calculations at M06L-GD3/def2-TZVPP level and by analyzing the crystal structures from the Cambridge Structural Database (CSD). One hydrogen bond where amino acid plays the role of the H-donor (NH/O), and two where it plays the role of the H-acceptor (O1/HO, O1 is a coordinated oxygen atom, and, O2/HO, O2 is a non-coordinated oxygen atom) were investigated. The calculations were done on octahedral nickel(II), square pyramidal copper(II), square planar copper(II), palladium(II), and platinum(II) glycine complexes with different charges adjusted using water(s) and/or chlorine ion(s) as the remaining ligands. For NH/O hydrogen bond, interaction energies of neutral complexes are the weakest, from -5.2 to -7.2 kcal/mol for neutral, stronger for singly positive, from -8.3 to -12.1 kcal/mol, and the strongest for doubly positive complex, -16.9 kcal/mol. For O1/HO and O2/HO interactions, neutral complexes have weaker interaction energies (from -2.2 to -5.1 kcal/mol for O1/HO, and from -3.7 to -5.0 kcal/mol for O2/HO), than for singly negative complexes (from -6.9 to -8.2 kcal/mol for O1/HO, and from -8.0 to -9.0 kcal/mol for O2/HO). Additionally to the complex charge, metal oxidation number, coordination number, and metal atomic number also influence the hydrogen bond strength, however, the influence is smaller.

Combined FTIR/Raman spectroscopic studies and ab initio electronic structure calculations of Dithiothreitol

A total of seven minimum energy geometries were obtained on exploring the conformational landscape of dithiothreitol (DTT) by varying the prominent dihedral angles in the molecule through a relaxed scan with a step size of 5° at B3LYP/cc-pVTZ with further geometry optimization at CCSD/cc-pVDZ level of theory. Single point energies were calculated for all the conformers at CCSD(T)/CBS limit with cc-pVNZ (N = T, Q) level of theory and revealed the similar energy pattern. The two conformers, namely G'TG'1 and G'TT, were found iso-energic even though they differed in their structure significantly and were of the lowest energy compared to others. Energies corresponding to the cyclic as well as other configurational counterpart of the global minimum were found much higher in energy compared to the global minimum structure. Intramolecular sulfur centered hydrogen bond was seen to stabilize the global minimum structure of DTT as revealed by AIM, NBO, FMO and ESP charge analysis. Computed NMR of DTT matched well with the experimental data gleaned from the literature. Vibrational spectra (Raman and IR) were measured and compared with computed normal modes of DTT, which were found in good agreement.

Stretching vibrational frequencies and pKa differences in H-bond networks of protein environments

The experimentally measured stretching vibrational frequencies of O-D [νO-D(donor)] and C=O [νC=O(donor)] H-bond donor groups can provide valuable information about the H-bonds in proteins. Here, using a quantum mechanical/molecular mechanical approach, the relationship between these vibrational frequencies and the difference in pKa values between H-bond donor and acceptor groups [ΔpKa(donor … acceptor)] in bacteriorhodopsin and photoactive yellow protein environments was investigated. The results show that νO-D(donor) is correlated with ΔpKa(donor … acceptor), regardless of the specific protein environment. νC=O(donor) is also correlated with ΔpKa(donor … acceptor), although the correlation is weak because the C=O bond does not have a proton. Importantly, the shifts in νO-D(donor) and νC=O(donor) are not caused by changes in pKa(donor) alone, but rather by changes in ΔpKa(donor … acceptor). Specifically, a decrease in ΔpKa(donor … acceptor) can lead to proton release from the H-bond donor group toward the acceptor group, resulting in shifts in the vibrational frequencies of the protein environment. These findings suggest that changes in the stretching vibrational frequencies, in particular νO-D(donor), can be used to monitor proton transfer in protein environments.

Structure of a ribonucleotide reductase R2 protein radical

Aerobic ribonucleotide reductases (RNRs) initiate synthesis of DNA building blocks by generating a free radical within the R2 subunit; the radical is subsequently shuttled to the catalytic R1 subunit through proton-coupled electron transfer (PCET). We present a high-resolution room temperature structure of the class Ie R2 protein radical captured by x-ray free electron laser serial femtosecond crystallography. The structure reveals conformational reorganization to shield the radical and connect it to the translocation path, with structural changes propagating to the surface where the protein interacts with the catalytic R1 subunit. Restructuring of the hydrogen bond network, including a notably short O-O interaction of 2.41 angstroms, likely tunes and gates the radical during PCET. These structural results help explain radical handling and mobilization in RNR and have general implications for radical transfer in proteins.

Sulfur Centered Hydrogen Bonding in Thioglycolic Acid and Its Clusters: A Computational Exploration

The conformational landscape of thioglycolic acid (TGA) was investigated by using the CCSD/cc-pVTZ level of theory. The GGC conformer was identified as the global minimum, followed by the GAC conformer. The calculated rotational constant for the GGC conformer exhibited good agreement with the previously reported experimental results. Subsequently, the study delved into the exploration of sulfur-centered hydrogen bonding in TGA's dimer and trimer clusters, employing the CCSD/cc-pVDZ level of theory. These clusters revealed the participation of both oxygen and sulfur atoms in noncovalent H-bonding, contributing to their stability. The presence of these noncovalent interactions in TGA clusters was elucidated through Atoms in Molecule (AIM), reduced density gradient (RDG), and natural bond order (NBO) analysis, while electrostatic potential (ESP) charge and vibrational mode analysis further supported these findings.

Discovery of a novel series of medium-sized cyclic enteropeptidase inhibitors

Enteropeptidase is located in the duodenum that involved in intestinal protein digestion. We have reported enteropeptidase inhibitors with low systemic exposure. The aim of this study was to discover novel enteropeptidase inhibitors showing more potent in vivo efficacy while retaining low systemic exposure. Inhibitory mechanism-based drug design led us to cyclize ester 2 to medium-sized lactones, showing potent enteropeptidase inhibitory activity and improving the ester stability, thus increasing fecal protein output in vivo. Optimization on the linker between two benzene rings resulted in discovery of ether lactone 6b, exhibiting further enhanced enteropeptidase inhibitory activity and long duration of inhibitory state. Oral administration of 6b in mice significantly elevated fecal protein output compared with the lead 2. In addition, 6b showed low systemic exposure along with low intestinal absorption. Furthermore, we identified the 10-membered lactonization method for scale-up synthesis of 6b, which does not require high-dilution conditions.

New insights into the proton pumping mechanism of ba3 cytochrome c oxidase: the functions of key residues and water

As the terminal oxidase of cell respiration in mitochondria and aerobic bacteria, the proton pumping mechanism of ba3-type cytochrome c oxidase (CcO) of Thermus thermophiles is still not fully understood. Especially, the functions of key residues which were considered as the possible proton loading sites (PLSs) above the catalytic center, as well as water located above and within the catalytic center, remain unclear. In this work, molecular dynamic simulations were performed on a set of designed mutants of key residues (Asp287, Asp372, His376, and Glu126II). The results showed that Asp287 may not be a PLS, but it could modulate the ability of the proton transfer pathway to transfer protons through its salt bridge with Arg225. Maintaining the closed state of the water pool above the catalytic center is necessary for the participation of inside water molecules in proton transfer. Water molecules inside the water pool can form hydrogen bond chains with PLS to facilitate proton transfer. Additional quantum cluster models of the Fe-Cu metal catalytic center are established, indicating that when the proton is transferred from Tyr237, it is more likely to reach the OCu atom directly through only one water molecule. This work provides a more profound understanding of the functions of important residues and specific water molecules in the proton pumping mechanism of CcO.

Human 2'-Deoxynucleoside 5'-Phosphate N-Hydrolase 1: Mechanism of 2'-Deoxyuridine 5'-Monophosphate Hydrolysis

The enzyme 2'-deoxynucleoside 5'-phosphate N-hydrolase 1 (DNPH1) catalyzes the N-ribosidic bond cleavage of 5-hydroxymethyl-2'-deoxyuridine 5'-monophosphate to generate 2-deoxyribose 5-phosphate and 5-hydroxymethyluracil. DNPH1 accepts other 2'-deoxynucleoside 5'-monophosphates as slow-reacting substrates. DNPH1 inhibition is a promising strategy to overcome resistance to and potentiate anticancer poly(ADP-ribose) polymerase inhibitors. We solved the crystal structure of unliganded human DNPH1 and took advantage of the slow reactivity of 2'-deoxyuridine 5'-monophosphate (dUMP) as a substrate to obtain a crystal structure of the DNPH1:dUMP Michaelis complex. In both structures, the carboxylate group of the catalytic Glu residue, proposed to act as a nucleophile in covalent catalysis, forms an apparent low-barrier hydrogen bond with the hydroxyl group of a conserved Tyr residue. The crystal structures are supported by functional data, with liquid chromatography-mass spectrometry analysis showing that DNPH1 incubation with dUMP leads to slow yet complete hydrolysis of the substrate. A direct UV-vis absorbance-based assay allowed characterization of DNPH1 kinetics at low dUMP concentrations. A bell-shaped pH-rate profile indicated that acid-base catalysis is operational and that for maximum kcat/KM, two groups with an average pKa of 6.4 must be deprotonated, while two groups with an average pKa of 8.2 must be protonated. A modestly inverse solvent viscosity effect rules out diffusional processes involved in dUMP binding to and possibly uracil release from the enzyme as rate limiting to kcat/KM. Solvent deuterium isotope effects on kcat/KM and kcat were inverse and unity, respectively. A reaction mechanism for dUMP hydrolysis is proposed.

Deciphering the driving forces in crystal packing by analysis of electrostatic energies and contact enrichment ratios

Hirshfeld surface analysis is a widely used tool for identifying the types of intermolecular contacts that contribute most significantly to crystal packing stabilization. One useful metric for analyzing these contacts is the contact enrichment descriptor, which indicates the types of contacts that are over- or under-represented. In this statistical study, enrichment ratios were combined with electrostatic energy (Eelec) data for a variety of compound families. To compute the electrostatic interaction energy between atoms, charge density models from the ELMAM2 database of multipolar atoms were used. As expected, strong hydrogen bonds such as O/N-H...N and O/N-H...O typically display large enrichment values and have the most negative (i.e. favorable) electrostatic energies. Conversely, contacts that are repulsive from an electrostatic perspective are usually the most under-represented. Analyzing the enrichment ratio and electrostatic energy indicators was shown to help identify which favorable contacts are the most competitive with each other. For weaker interactions, such as hydrophobic contacts, the behavior is less clear cut and can depend on other factors such as the chemical content of the molecule. The anticorrelation between contact enrichment and Eelec is generally lost for weaker contacts. However, we observed that C...C contacts are often enriched in crystal structures containing heterocycles, despite the low electrostatic attraction. For molecules with only weak hydrogen bond donors/acceptors and hydrophobic groups, the correlation between contact enrichment and Eelec is still evident for the strongest of these interactions. However, there are some exceptions where the most favorable contacts from an electrostatic perspective are not the most over-represented. This can occur in cases where the shape of the molecule is complex or elongated, favoring dispersion forces and shape complementarity in the packing.

Chemical features and machine learning assisted predictions of protein-ligand short hydrogen bonds

There are continuous efforts to elucidate the structure and biological functions of short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms reside more than 0.3 Å closer than the sum of their van der Waals radii. In this work, we evaluate 1070 atomic-resolution protein structures and characterize the common chemical features of SHBs formed between the side chains of amino acids and small molecule ligands. We then develop a machine learning assisted prediction of protein-ligand SHBs (MAPSHB-Ligand) model and reveal that the types of amino acids and ligand functional groups as well as the sequence of neighboring residues are essential factors that determine the class of protein-ligand hydrogen bonds. The MAPSHB-Ligand model and its implementation on our web server enable the effective identification of protein-ligand SHBs in proteins, which will facilitate the design of biomolecules and ligands that exploit these close contacts for enhanced functions.

Proton-Tunneling Dynamics along Low-Barrier Hydrogen Bonds: A Full-Dimensional Instanton Study of 6-Hydroxy-2-formylfulvene

Understanding the dynamics of proton transfer along low-barrier hydrogen bonds remains an outstanding challenge of great fundamental and practical interest, reflecting the central role of quantum effects in reactions of chemical and biological importance. Here, we combine ab initio calculations with the semiclassical ring-polymer instanton method to investigate tunneling processes on the ground electronic state of 6-hydroxy-2-formylfulvene (HFF), a prototypical neutral molecule supporting low-barrier hydrogen-bonding. The results emerging from a full-dimensional ab initio instanton analysis reveal that the tunneling path does not pass through the instantaneous transition-state geometry. Instead, the tunneling process involves a multidimensional reaction coordinate with concerted reorganization of the heavy-atom skeletal framework to substantially reduce the donor-acceptor distance and drive the ensuing intramolecular proton-transfer event. The predicted tunneling-induced splittings for HFF isotopologues are in good agreement with experimental findings, leading to percentage deviations of only 20-40%. Our full-dimensional results allow us to characterize vibrational contributions along the tunneling path, highlighting the intrinsically multidimensional nature of the attendant hydron-migration dynamics.

Length-Dependent Transition from Extended to Folded Shapes in Short Oligomers of an Azetidine-Based α-Amino Acid: The Critical Role of NH···N H-Bonds

Hydrogen bonds (H-bonds) are ubiquitous in peptides and proteins and are central to the stabilization of their structures. Inter-residue H-bonds between non-adjacent backbone amide NH and C=O motifs lead to the well-known secondary structures of helices, turns and sheets, but it is recognized that other H-bonding modes may be significant, including the weak intra-residue H-bond (called a C5 H-bond) that implicates the NH and C=O motifs of the same amino acid residue. Peptide model compounds that adopt stable C5 H-bonds are not readily available and the so-called 2.0₅-helix, formed by successive C5 H-bonds, is an elusive secondary structure. Using a combination of theoretical chemistry and spectroscopic studies in both the gas phase and solution phase, we have demonstrated that derivatives of 3-amino-1-methylazetidine-3-carboxylic acid, Aatc(Me) can form sidechain-backbone N-H···N C6γ H-bonds that accompany-and thereby stabilize-C5 H-bonds. In the capped trimer of Aatc(Me), extended C5/C6γ motifs are sufficiently robust to challenge classical 3₁₀-helix formation in solution and the fully-extended 2.0₅-helix conformer has been characterized in the gas phase. Concurrent H-bonding support for successive C5 motifs is a new axiom for stabilizing the extended backbone secondary structure in short peptides.

First-principles density functional theoretical study on the structures, reactivity and spectroscopic properties of (NH) and (OH) Tautomer's of 4-(methylsulfanyl)-3[(1Z)-1-(2-phenylhydrazinylidene) ethyl] quinoline-2(1H)-one

The tautomerizations mechanism of 4-(methylsulfanyl)-3[(1Z)-1-(2-phenylhydrazinylidene) ethyl] quinoline-2(1H)-one were inspected in the gas phase and ethanol using density function theory (DFT) M06-2X and B3LYP methods. Thermo-kinetic features of different conversion processes were estimated in temperature range 273-333 K using the Transition state theory (TST) accompanied with one dimensional Eckert tunneling correction (1D-Eck). Acidity and basicity were computed as well, and the computational results were compared against the experimental ones. Additionally, NMR, global descriptors, Fukui functions, NBO charges, and electrostatic potential (ESP) were discussed. From thermodynamics analysis, the keto form of 4-(methylsulfanyl)-3-[(1Z)-1-(2 phenylhydrazinylidene) quinoline-2(1H)-one is the most stable form in the gas phase and ethanol and the barrier heights required for tautomerization process were found to be high in the gas phase and ethanol ~ 38.80 and 37.35 kcal/mol, respectively. DFT methods were used for UV-Vis electronic spectra simulation and the time-dependent density functional theory solvation model (TDDFT-SMD) in acetonitrile compounds.

Biomolecular condensates: Formation mechanisms, biological functions, and therapeutic targets

Biomolecular condensates are cellular structures composed of membraneless assemblies comprising proteins or nucleic acids. The formation of these condensates requires components to change from a state of solubility separation from the surrounding environment by undergoing phase transition and condensation. Over the past decade, it has become widely appreciated that biomolecular condensates are ubiquitous in eukaryotic cells and play a vital role in physiological and pathological processes. These condensates may provide promising targets for the clinic research. Recently, a series of pathological and physiological processes have been found associated with the dysfunction of condensates, and a range of targets and methods have been demonstrated to modulate the formation of these condensates. A more extensive description of biomolecular condensates is urgently needed for the development of novel therapies. In this review, we summarized the current understanding of biomolecular condensates and the molecular mechanisms of their formation. Moreover, we reviewed the functions of condensates and therapeutic targets for diseases. We further highlighted the available regulatory targets and methods, discussed the significance and challenges of targeting these condensates. Reviewing the latest developments in biomolecular condensate research could be essential in translating our current knowledge on the use of condensates for clinical therapeutic strategies.

Brønsted Acids Promote Olefin Oxidations by Bioinspired Nonheme CoIII(PhIO)(OH) Complexes: A Role for Low-Barrier Hydrogen Bonds

Introduction of Brønsted acids into biomimetic nonheme reactions promotes the oxidative ability of metal-oxygen complexes significantly. However, the molecular machinery of the promoted effects is missing. Herein, a comprehensive investigation of styrene oxidation by a cobalt(III)-iodosylbenzene complex, [(TQA)Co^III(OIPh)(OH)]²⁺ (1, TQA = tris(2-quinolylmethyl)amine), in the presence and absence of triflic acid (HOTf) was performed using density functional theory calculations. Results revealed for the first time that there is a low-barrier hydrogen bond (LBHB) between HOTf and the hydroxyl ligand of 1, which forms two valence-resonance structures [(TQA)Co^III(OIPh)(HO^---HOTf)]²⁺ (1_LBHB) and [(TQA)Co^III(OIPh)(H₂O--OTf^-)]²⁺ (1'_LBHB). Due to the oxo-wall, these complexes (1_LBHB and 1'_LBHB) cannot convert to high-valent cobalt-oxyl species. Instead, styrene oxidation by these oxidants (1_LBHB and 1'_LBHB) shows novel spin-state selectivity, i.e., on the ground closed-shell singlet state, styrene is oxidized to an epoxide, whereas on the excited triplet and quintet states, an aldehyde product, phenylacetaldehyde, is formed. The preferred pathway is styrene oxidation by 1'_LBHB, which is initiated by a rate-limiting bond-formation-coupled electron transfer process with an energy barrier of 12.2 kcal mol^-1. The nascent PhIO-styrene-radical-cation intermediate undergoes an intramolecular rearrangement to produce an aldehyde. The halogen bond between the OH^-/H₂O ligand and the iodine of PhIO modulates the activity of the cobalt-iodosylarene complexes 1_LBHB and 1'_LBHB. These new mechanistic findings enrich our knowledge of nonheme chemistry and hypervalent iodine chemistry and will play a positive role in the rational design of new catalysts.

Exploring the Non-Covalent Bonding in Water Clusters

QTAIM and source function analysis were used to explore the non-covalent bonding in twelve different water clusters (H2O)n obtained by considering n = 2–7 and various geometrical arrangements. A total of seventy-seven O−H⋯O hydrogen bonds (HBs) were identified in the systems under consideration, and the examination of the electron density at the bond critical point (BCP) of these HBs revealed the existence of a great diversity of O−H⋯O interactions. Furthermore, the analysis of quantities, such as |V(r)|/G(r) and H(r), allowed a further description of the nature of analogous O−H⋯O interactions within each cluster. In the case of 2-D cyclic clusters, the HBs are nearly equivalent between them. However, significant differences among the O−H⋯O interactions were observed in 3-D clusters. The assessment of the source function (SF) confirmed these findings. Finally, the ability of SF to decompose the electron density (ρ) into atomic contributions allowed the evaluation of the localized or delocalized character of these contributions to ρ at the BCP associated to the different HBs, revealing that weak O−H⋯O interactions have a significant spread of the atomic contributions, whereas strong interactions have more localized atomic contributions. These observations suggest that the nature of the O−H⋯O hydrogen bond in water clusters is determined by the inductive effects originated by the different spatial arrangements of the water molecules in the studied clusters.

Intermolecular hydrogen bonding of alcohols with dinitrobenzene radical anion and dianion: A combined electrochemical and DFT study

Hydrogen bonding is one of the most important inter-molecular interactions in the field of biochemistry and medicinal chemistry. Such non-covalent interactions play a vital role in self-assembly phenomena, chemical structures, material properties and enzymatic catalysis. Herein, we present hydrogen bonding phenomenon in alcohols-dinitrobenzene (DNB) radical anion/dianion systems using electrochemical and computational approaches. First, 1,3-DNB radical anion and dianion were generated through electrochemical method and then hydrogen bonding interactions with aliphatic alcohols in DMSO are studied through cyclic voltammetry (CV). CV results show that the cathodic peak potential of 1,3-Dinitrobenzene in Dimethyl sulfoxide is shifted catholically upon addition of alcohols which represent hydrogen bonding phenomenon. Theoretical investigations are performed to gain deep insight on hydrogen bonding interaction strength in DNB-alcohol systems. H-bonding interaction of all isomers of DNB (1,2-, 1,3-, 1,4-), its corresponding radical anion, and dianion with aliphatic alcohols is studied using density functional calculations. The strength of H-bonding is evaluated both qualitatively and quantitatively using interaction energies, vibrational and electronic spectroscopic analysis. Understanding of these interactions will be helpful in gaining insight into biological systems where these interactions play significant role.

The effect of spin exchange interaction on protein structural stability

Partially charged chiral molecules act as spin filters, with preference for electron transport toward one type of spin ("up" or "down"), depending on their handedness. This effect is named the chiral induced spin selectivity (CISS) effect. A consequence of this phenomenon is spin polarization concomitant with electric polarization in chiral molecules. These findings were shown by adsorbing chiral molecules on magnetic surfaces and investigating the spin-exchange interaction between the surface and the chiral molecule. This field of study was developed using artificial chiral molecules. Here we used such magnetic surfaces to explore the importance of the intrinsic chiral properties of proteins in determining their stability. First, proteins were adsorbed on paramagnetic and ferromagnetic nanoparticles in a solution, and subsequently urea was gradually added to induce unfolding. The structural stability of proteins was assessed using two methods: bioluminescence measurements used to monitor the activity of the Luciferase enzyme, and fast spectroscopy detecting the distance between two chromophores implanted at the termini of a Barnase core. We found that interactions with magnetic materials altered the structural and functional resilience of the natively folded proteins, affecting their behavior under varying mild denaturing conditions. Minor structural disturbances at low urea concentrations were impeded in association with paramagnetic nanoparticles, whereas at higher urea concentrations, major structural deformation was hindered in association with ferromagnetic nanoparticles. These effects were attributed to spin exchange interactions due to differences in the magnetic imprinting properties of each type of nanoparticle. Additional measurements of proteins on macroscopic magnetic surfaces support this conclusion. The results imply a link between internal spin exchange interactions in a folded protein and its structural and functional integrity on magnetic surfaces. Together with the accumulating knowledge on CISS, our findings suggest that chirality and spin exchange interactions should be considered as additional factors governing protein structures.

Promoting Co-Crystallization in Poly(butylene succinate) and Poly(butylene fumarate) Blends via End-Group Functionalization

Co-crystallization plays a crucial role in the integration and regulation of thermal and mechanical properties in polymer blends, but the poor compatibility of the components in the crystal phase has always been a major obstacle to co-crystallization, which puts forward stricter requests for linkage and interaction between different entities. On the basis of the hydrogen-bonding interaction that can promote chain stacking and thus improve miscibility, we propose that crystalline/crystalline blends of 2-ureido-4[1H]-pyrimidinone (UPy)-functionalized poly(butylene succinate) and poly(butylene fumarate) (PBS-UPy/PBF-UPy) where UPy groups with quadruple hydrogen-bonding interaction are employed to connect different chain ends, could inhibit phase separation and improve co-crystallization. PBS-UPy/PBF-UPy blends exhibit complex component-dependent and cooling-rate-dependent co-crystallization behavior. A high level of co-crystallization occurs in the range of PBS-UPy-rich fractions, and the proportion could approach over 98% under optimized conditions with the aid of UPy quadruple hydrogen bonds interaction. This work enriches the understanding of co-crystallization in crystalline/crystalline polymer blends and provides more possibility for the design of structures and properties of polymer materials.

Ensemble-function relationships to dissect mechanisms of enzyme catalysis

Decades of structure-function studies have established our current extensive understanding of enzymes. However, traditional structural models are snapshots of broader conformational ensembles of interchanging states. We demonstrate the need for conformational ensembles to understand function, using the enzyme ketosteroid isomerase (KSI) as an example. Comparison of prior KSI cryogenic x-ray structures suggested deleterious mutational effects from a misaligned oxyanion hole catalytic residue. However, ensemble information from room-temperature x-ray crystallography, combined with functional studies, excluded this model. Ensemble-function analyses can deconvolute effects from altering the probability of occupying a state (P-effects) and changing the reactivity of each state (k-effects); our ensemble-function analyses revealed functional effects arising from weakened oxyanion hole hydrogen bonding and substrate repositioning within the active site. Ensemble-function studies will have an integral role in understanding enzymes and in meeting the future goals of a predictive understanding of enzyme catalysis and engineering new enzymes.

Theoretical study on L-H+-L with identical donors: short strong hydrogen bond or not?

Short strong hydrogen bonds (SSHBs) play crucial role in many chemical processes. Recently, as the representative of SSHBs, [F-H-F]- was experimentally observed. [F-H-F]- has a symmetric structure, which can be described as a H+ acid shared by two terminal F- donors (F--H+-F-). To explore whether two identical donors are bound to result in SSHBs, we performed theoretical studies on a series of compounds (L-H+-L) with two identical electron donors (L corresponds to donors containing group 14, 15, 16 and 17 elements). The results show that identical donors do not definitely lead to SSHBs. Instead, typical hydrogen bonds also exist. We found that both electronegativity and basicity contribute to the patterns of hydrogen bonds, where more electronegative and weaker donors benefit to SSHBs. Besides, it was found that zero-point energies also respond to the hydrogen bonding systems. This systemic work is expected to provide more insights into SSHBs.

Promotion of collagen mineralization and dentin repair by succinates

Biomineralization of collagen fibers is regulated by non-collagenous proteins and small biomolecules, which are essential in bone and teeth formation. In particular, small biomolecules such as succinic acid (SA) exist at a high level in hard tissues, but their role is yet unclear. Here, our work demonstrated that SA could significantly promote intrafibrillar mineralization in two- and three-dimensional collagen models, where the relative mineralization rate was 16 times faster than the control group. Furthermore, the FTIR spectra and isothermal experimental results showed that collagen molecules could interact with SA via a hydrogen bond and that the interaction energy was about 4.35 kJ mol^-1. As expected, the SA-pretreated demineralized dentin obtained full remineralization within two days, whereas it took more than four days in the control group, and their mechanical properties were considerably enhanced compared with those of the demineralized one. The possible mechanism of the promotion effect of SA was ultimately illustrated, with SA modification strengthening the capacity of the collagen matrix to attract more calcium ions, which might create a higher local concentration that could accelerate the mineralization of collagen fibers. These findings not only advance the understanding of the vital role of small biomolecules in collagen biomineralization but also facilitate the development of an effective strategy to repair hard tissues.

Small Molecules Promote Selective Denaturation and Degradation of Tubulin Heterodimers through a Low-Barrier Hydrogen Bond

Here, we report a novel mechanism to selectively degrade target proteins. 3-(3-Phenoxybenzyl)amino-β-carboline (PAC), a tubulin inhibitor, promotes selective degradation of αβ-tubulin heterodimers. Biochemical studies have revealed that PAC specifically denatures tubulin, making it prone to aggregation that predisposes it to ubiquitinylation and then degradation. The degradation is mediated by a single hydrogen bond formed between the pyridine nitrogen of PAC and βGlu198, which is identified as a low-barrier hydrogen bond (LBHB). In contrast, another two tubulin inhibitors that only form normal hydrogen bonds with βGlu198 exhibit no degradation effect. Thus, the LBHB accounts for the degradation. We then screened for compounds capable of forming an LBHB with βGlu198 and demonstrated that BML284, a Wnt signaling activator, also promotes tubulin heterodimer degradation through the LBHB. Our study provided a unique example of LBHB function and identified a novel approach to obtain tubulin degraders.

Release of a Proton and Formation of a Low-Barrier Hydrogen Bond between Tyrosine D and D2-His189 in Photosystem II

In photosystem II (PSII), the second-lowest oxidation state (S₁) of the oxygen-evolving Mn₄CaO₅ cluster is the most stable, as the radical form of the redox-active D2-Tyr160 is considered to be a candidate that accepts an electron from the lowest oxidation state (S₀) in the dark. Using quantum mechanical/molecular mechanical calculations, we investigated the redox potential (E _m) of TyrD and its H-bond partner, D2-His189. The potential energy profile indicates that the release of a proton from the TyrD...D2-His189 pair leads to the formation of a low-barrier H-bond. The E _m depends on the H⁺ position along the low-barrier H-bond, e.g., 680 mV when the H⁺ is at the D2-His189 moiety and 800 mV when the H⁺ is at the TyrD moiety, which can explain why TyrD mediates both the S₀ to S₁ oxidation and the S₂ to S₁ reduction.

Self-Organization Dynamics of Collagen-like Peptides Crosslinking Is Driven by Rose-Bengal-Mediated Electrostatic Bridges

The present work focuses on the computational study of the structural micro-organization of hydrogels based on collagen-like peptides (CLPs) in complex with Rose Bengal (RB). In previous studies, these hydrogels computationally and experimentally demonstrated that when RB was activated by green light, it could generate forms of stable crosslinked structures capable of regenerating biological tissues such as the skin and cornea. Here, we focus on the structural and atomic interactions of two collagen-like peptides (collagen-like peptide I (CLPI), and collagen-like peptide II, (CLPII)) in the presence and absence of RB, highlighting the acquired three-dimensional organization and going deep into the stabilization effect caused by the dye. Our results suggest that the dye could generate a ternary ground-state complex between collagen-like peptide fibers, specifically with positively charged amino acids (Lys in CLPI and Arg in CLPII), thus stabilizing ordered three-dimensional structures. The discoveries generated in this study provide the structural and atomic bases for the subsequent rational development of new synthetic peptides with improved characteristics for applications in the regeneration of biological tissues during photochemical tissue bonding therapies.

Hydrogen-Bonds Mediated Nanomedicine: Design, Synthesis and Applications

Among the various challenges in medicine, diagnosis, complete cure and healing of cancers remain difficult given the heterogeneity and complexity of such disease. Differing from conventional platforms with often unsatisfactory theranostic capabilities, the contribution of supramolecular interactions, such as hydrogen-bonds (H-bonds), to cancer nanotheranostics opens new perspectives for the design of biomedical materials, exhibiting remarkable properties and easier processability. Thanks to their dynamic characteristics, a feature generally observed for non-covalent interactions, H-bonding (macro)molecules can be used as supramolecular motifs for yielding drug- and diagnostic carriers that possess attractive features, arising from the combination of assembled nanoplatforms and the responsiveness of H-bonds. Thus H-bonded nanomedicine provides a rich toolbox that is useful to fulfill biomedical needs with unique advantages in early-stage diagnosis and therapy, demonstrating the promising potential in clinical translations and applications. We here summarize the design and synthetic routes towards H-bonded nanomedicines, focus on the growing understanding of the structure-function relationship for efficient cancer treatment. We propose a guidance for designing new H-bonded intelligent theranostic agents, to inspire more successful explorations of cancer nanotheranostics and finally to promote potential clinical translations. This article is protected by copyright. All rights reserved.

Absorption wavelength along chromophore low-barrier hydrogen bonds

In low-barrier hydrogen bonds (H-bonds), the pK_a values for the H-bond donor and acceptor moieties are nearly equal, whereas the redox potential values depend on the H⁺ position. Spectroscopic details of low-barrier H-bonds remain unclear. Here, we report the absorption wavelength along low-barrier H-bonds in protein environments, using a quantum mechanical/molecular mechanical approach. Low-barrier H-bonds form between Glu46 and p-coumaric acid (pCA) in the intermediate pR_CW state of photoactive yellow protein and between Asp116 and the retinal Schiff base in the intermediate M-state of the sodium-pumping rhodopsin KR2. The H⁺ displacement of only ∼0.4 Å, which does not easily occur without low-barrier H-bonds, is responsible for the ∼50 nm-shift in the absorption wavelength. This may be a basis of how photoreceptor proteins have evolved to proceed photocycles using abundant protons.

•

The low-barrier H-bond formation is a prerequisite for proton transfer

•

How the absorption wavelength changes as H⁺ moves is an open question

•

The H⁺ displacement of ∼0.4 Å leads to the absorption wavelength shift of ∼50 nm

•

The localization of the molecular orbitals plays a key role in the wavelength shift