Strength of a bifurcated H bond

Hydrogen Bond Strengthens Acceptor Group: The Curious Case of the C-H···O=C Bond

An H-bond involves the sharing of a hydrogen atom between an electronegative atom to which it is covalently bound (the donor) and another electronegative atom serving as an acceptor. Such bonds represent a critically important geometrical force in biological macromolecules and, as such, have been characterized extensively. H-bond formation invariably leads to a weakening within the acceptor moiety due to the pulling exerted by the donor hydrogen. This phenomenon can be compared to a spring connecting two masses; pulling one mass stretches the spring, similarly affecting the bond between the two masses. Herein, we describe the opposite phenomenon when investigating the energetics of the C-H···O=C bond. This bond underpins the most prevalent protein transmembrane dimerization motif (GxxxG) in which a glycine Cα-H on one helix forms a hydrogen bond with a carbonyl in a nearby helix. We use isotope-edited FT-IR spectroscopy and corroborating computational approaches to demonstrate a surprising strengthening of the acceptor C=O bond upon binding with the glycine Cα-H. We show that electronic factors associated with the Cα-H bond strengthen the C=O oscillator by increasing the s-character of the σ-bond, lowering the hyperconjugative disruption of the π-bond. In addition, a reduction of the acceptor C=O bond's polarity is observed upon the formation of the C-H···O=C bond. Our findings challenge the conventional understanding of H-bond dynamics and provide new insights into the structural stability of inter-helical protein interactions.

Stabilizing Bifurcated Hydrogen Bond in 8-Aminoquinoline Appended Peptides

The hydrogen bonding interaction between an amide N-H and the amide N of the preceding residue is prevalent in proline-containing proteins and peptides. However, the N-H⋅⋅⋅N hydrogen bonding interaction is rare in non-prolyl natural peptides due to restricted dihedral angles. Herein, we stabilize this type of interaction in 8-aminoquinoline appended non-prolyl peptides through bifurcated N⋅⋅⋅H⋅⋅⋅N hydrogen bond. The 8-aminoquinoline-incorporated model peptides 2 a-i were designed, synthesized, and the crystal structures of 2 a-c and 2 i were solved. Analysis of crystal data reveals that the amide N-H of aminoquinoline is involved in bifurcated hydrogen bonding interaction with the nitrogen of the preceding amino acid residue and the nitrogen in quinoline. Analysis of crystal packing, Hirshfeld surface and fingerprint plots confirms that the intermolecular O⋅⋅⋅H contacts significantly contribute to stabilizing bifurcated N⋅⋅⋅H⋅⋅⋅N hydrogen bonding interaction. Furthermore, NMR experiments and CD spectroscopy were conducted to examine the preferred conformation in solution, and the data corroborate with the crystal structure conformation.

Cholesterol and Lipid Rafts in the Biogenesis of Amyloid-β Protein and Alzheimer's Disease

Cholesterol has been conjectured to be a modulator of the amyloid cascade, the mechanism that produces the amyloid-β (Aβ) peptides implicated in the onset of Alzheimer's disease. We propose that cholesterol impacts the genesis of Aβ not through direct interaction with proteins in the bilayer, but indirectly by inducing the liquid-ordered phase and accompanying liquid-liquid phase separations, which partition proteins in the amyloid cascade to different lipid domains and ultimately to different endocytotic pathways. We explore the full process of Aβ genesis in the context of liquid-ordered phases induced by cholesterol, including protein partitioning into lipid domains, mechanisms of endocytosis experienced by lipid domains and secretases, and pH-controlled activation of amyloid precursor protein secretases in specific endocytotic environments. Outstanding questions on the essential role of cholesterol in the amyloid cascade are identified for future studies.

Fluoride-Based Deep Eutectic Solvents with Amide Dual-Hydrogen-Bond Donors

Developing F--containing electrolytes is crucial for electrochemical and chemical fluorination. However, balancing the F- concentration and electrochemical stability of the electrolytes remains a challenge. In this study, fluoride-based deep eutectic solvents (F-DESs) were obtained by using amide hydrogen-bond donors (HBDs) containing dual N-H bonds. The obtained F-DES, [TMA]F·3.5[1,3-DMU], was prepared by facilely mixing solid compounds of tetramethylammonium fluoride ([TMA]F) and 1,3-dimethylurea (1,3-DMU), resulting in a high F- concentration (2.6 mol dm-3) and a wide electrochemical window (3.1 V) at room temperature. The electrochemical window was much wider than that of [TMA]F·3.5[EG] (EG, ethylene glycol) as another F-DES with an alcohol HBD (1.9 V). Moreover, [TMA]F·3.5[1,3-DMU] exhibited an ionic conductivity that was 2 orders of magnitude higher than that of [TMA]F·3.5[1,3-DMTU] (1,3-DMTU, 1,3-dimethylthiourea) around room temperature because of the bifurcated hydrogen bonds between the dual N-H bonds of 1,3-DMU and one F-. Thus, [TMA]F·3.5[1,3-DMU] was demonstrated to be applicable to electrochemical fluorination.

Molecular dynamics and NMR reveal the coexistence of H-bond-assisted and through-space JFH coupling in fluorinated amino alcohols

The JFH coupling constants in fluorinated amino alcohols were investigated through experimental and theoretical approaches. The experimental JFH couplings were only reproduced theoretically when explicit solvation through molecular dynamics (MD) simulations was conducted in DMSO as the solvent. The combination of MD conformation sampling and DFT NMR spin-spin coupling calculations for these compounds reveals the simultaneous presence of through-space (TS) and hydrogen bond (H-bond) assisted JFH coupling between fluorine and hydrogen of the NH group. Furthermore, MD simulations indicate that the hydrogen in the amino group participates in both an intermolecular bifurcated H-bond with DMSO and in transmitting the observed JFH coupling. The contribution of TS to the JFH coupling is due to the spatial proximity of the fluorine and the NH group, aided by a combination of the non-bonding transmission pathway and the hydrogen bonding pathway. The experimental JFH coupling observed for the molecules studied should be represented as 4TS/1hJFH coupling.

N3-Methyluridine and 2'-O-Alkyl/2'-Fluoro-N3-methyluridine functionalized nucleic acids improve nuclease resistance while maintaining duplex geometry

Herein, we report the synthesis of 2'-O-alkyl/2'-fluoro-N3-methyluridine (2'-O-alkyl/2'-F-m3U) phosphoramidites and their incorporation in DNA and RNA oligonucleotides. The duplex binding affinity and base discrimination studies showed that all 2'-O-alkyl/2'-F-m3U modifications significantly decreased the thermal stability and base-pairing discrimination ability. Serum stability study of dT20 with 2'-O-alkyl-m3U modification exhibited excellent nuclease resistance when incubated with 3'-exonucleases (SVPD) or 5'-exonucleases (PDE-II) as compared to m3U, 2'-F, 2'-OMe modified oligonucleotides. MD simulation studies with RNA tetradecamer duplexes illustrated that the m3U and 2'-O-methyl-m3U modifications reduce the duplex stabilities by disrupting the Watson-Crick hydrogen bonding and base-stacking interactions. Further molecular modelling investigations demonstrated that the 2'-O-propyl-m3U modification exhibits steric interactions with amino acid residues in the active site of 3'- and 5'-exonuclease, leading to enhanced stability. These combined data indicate that the 2'-modified-m3U nucleotides can be used as a promising tool to enhance the stability, silencing efficiency, and drug-like properties of antisense/siRNA-based therapeutics.

Dimer Stabilization by SpyTag/SpyCatcher Coupling of the Reductase Domains of a Chimeric P450 BM3 Monooxygenase from Bacillus spp. Improves its Stability, Activity, and Purification

The vast majority of known enzymes exist as oligomers, which often gives them high catalytic performance but at the same time imposes constraints on structural conformations and environmental conditions. An example of an enzyme with a complex architecture is the P450 BM3 monooxygenase CYP102A1 from Bacillus megaterium. Only active as a dimer, it is highly sensitive to dilution or common immobilization techniques. In this study, we engineered a thermostable P450BM3 chimera consisting of the heme domain of a CYP102A1 variant and the reductase domain of the homologous CYP102A3. The dimerization of the hybrid was even weaker compared to the corresponding CYP102A1 variant. To create a stable dimer, we covalently coupled the C-termini of two monomers of the chimera via SpyTag003/SpyCatcher003 interaction. As a result, purification, thermostability, pH stability, and catalytic activity were improved. Via a bioorthogonal two-step affinity purification, we obtained high purity (94 %) of the dimer-stabilized variant being robust against heme depletion. Long-term stability was increased with a half-life of over 2 months at 20 °C and 80-90 % residual activity after 2 months at 5 °C. Most catalytic features were retained with even an enhancement of the overall activity by ~2-fold compared to the P450BM3 chimera without SpyTag003/SpyCatcher003.

Photoinduced radical formation in hydrogen-bonded organic frameworks

Hydrogen-bonded organic frameworks (HOFs) constructed from naphthalene-diimide bearing tectons undergo photochromic changes whilst forming radical bearing species within the framework structure.

Cucurbit[7]uril Enhances Distance Measurements of Spin-Labeled Proteins

We report complex formation between the chloroacetamide 2,6-diazaadamantane nitroxide radical (ClA-DZD) and cucurbit[7]uril (CB-7), for which the association constant in water, Ka = 1.9 × 106 M-1, is at least 1 order of magnitude higher than the previously studied organic radicals. The radical is highly immobilized by CB-7, as indicated by the increase in the rotational correlation time, τrot, by a factor of 36, relative to that in the buffer solution. The X-ray structure of ClA-DZD@CB-7 shows the encapsulated DZD guest inside the undistorted CB-7 host, with the pendant group protruding outside. Upon addition of CB-7 to T4 Lysozyme (T4L) doubly spin-labeled with the iodoacetamide derivative of DZD, we observe the increase in τrot and electron spin coherence time, Tm, along with the narrowing of interspin distance distributions. Sensitivity of the DEER measurements at 83 K increases by a factor 4-9, compared to the common spin label such as MTSL, which is not affected by CB-7. Interspin distances of 3 nm could be reliably measured in water/glycerol up to temperatures near the glass transition/melting temperature of the matrix at 200 K, thus bringing us closer to the goal of supramolecular recognition-enabled long-distance DEER measurements at near physiological temperatures. The X-ray structure of DZD-T4L 65 at 1.12 Å resolution allows for unambiguous modeling of the DZD label (0.88 occupancy), indicating an undisturbed structure and conformation of the protein.

n → π* Interaction Enabling Transient Inversion of Chirality

This study reports the observation and characterization of two isomers of the acrolein dimer by using high-resolution rotational spectroscopy in pulsed jets. The first isomer is stabilized by two hydrogen bonds, adopting a planar configuration, and is energetically favored over the second isomer, which exhibits a dominant n → π* interaction in a nearly orthogonal arrangement. Surprisingly, the n → π* interaction was revealed to enable a concerted tunneling motion of two moieties along the carbonyl group. This motion leads to the inversion of transient chirality associated with the exchange of donor-acceptor roles, as revealed by the spectral feature of quadruplets. Inversion of transient chirality is a fundamental phenomenon in quantum mechanics and commonly observed for only inversional motions of protons. It is the first discovery, to the best of our knowledge, that such heavy moieties can also undergo chirality inversion.

Theoretical investigation of tube-like supramolecular structures formed through bifurcated lithium bonds

The stability of three supramolecular naostructures, which are formed through the aggregation of identical belts of [12] arene containing p-nitrophenyllithium, 1,4-dilithiatedbenzene and 1,4-dinitrobenzene units, is investigated by density functional theory. The electrostatic potential calculations indicate the ability of these belts in forming bifurcated lithium bonds (BLBs) between the Li atoms of one belt and the oxygen atoms of the NO2 groups in the other belt, which is also confirmed by deformation density maps and quantum theory of atoms in molecules (QTAIM) analysis. Topological analysis and natural bond analysis (NBO) imply to ionic character for these BLBs with binding energies up to approximately - 60 kcal mol-1. The many-body interaction energy analysis shows the strong cooperativity belongs to the configuration with the highest symmetry (C4v) containing p-nitrophenyllithium fragments as the building unit. Therefore, it seems that this configuration could be a good candidate for designing a BLB-based supramolecular nanotube with infinite size in this study.

High-resolution crystal structure of the Mu8.1 conotoxin from Conus mucronatus

Marine cone snails produce a wealth of peptide toxins (conotoxins) that bind their molecular targets with high selectivity and potency. Therefore, conotoxins constitute valuable biomolecular tools with a variety of biomedical purposes. The Mu8.1 conotoxin from Conus mucronatus is the founding member of the newly identified saposin-like conotoxin class of conotoxins and has been shown to target Cav2.3, a voltage-gated calcium channel. Two crystal structures have recently been determined of Mu8.1 at 2.3 and 2.1 Å resolution. Here, a high-resolution crystal structure of Mu8.1 was determined at 1.67 Å resolution in the high-symmetry space group I4122. The asymmetric unit contained one molecule, with a symmetry-related molecule generating a dimer equivalent to that observed in the two previously determined structures. The high resolution allows a detailed atomic analysis of a water-filled cavity buried at the dimer interface, revealing a tightly coordinated network of waters that shield a lysine residue (Lys55) with a predicted unusually low side-chain pKa value. These findings are discussed in terms of a potential functional role of Lys55 in target interaction.

Cucurbit[7]uril Enhances Distance Measurements of Spin-Labeled Proteins

We report complex formation between the chloroacetamide 2,6-diazaadamantane nitroxide radical (ClA-DZD) and cucurbit[7]uril (CB-7), for which the association constant in water, Ka = 1.9 × 106 M-1, is at least one order of magnitude higher than the previously studied organic radicals. The radical is highly immobilized by CB-7, as indicated by the increase of the rotational correlation time, τrot, by a factor of 36, relative to that in the buffer solution. The X-ray structure of ClA-DZD@CB-7 shows the encapsulated DZD guest inside the undistorted CB-7 host, with the pendant group protruding outside. Upon addition of CB-7 to T4 Lysozyme (T4L) doubly spin-labeled with the iodoacetamide derivative of DZD, we observe the increase in τrot and electron spin coherence time, Tm, along with the narrowing of inter-spin distance distributions. Sensitivity of the DEER measurements at 83 K increases by a factor 4 - 9, compared to the common spin label such as MTSL, which is not affected by CB-7. Inter-spin distances of 3-nm could be reliably measured in water/glycerol up to temperatures near the glass transition/melting temperature of the matrix at 200 K, thus bringing us closer to the goal of supramolecular recognition-enabled long-distance DEER measurements at near physiological temperatures. The X-ray structure of DZD-T4L 65 at 1.12 Å resolution allows for unambiguous modeling of the DZD label (0.88 occupancy), indicating undisturbed structure and conformation of the protein.

Synergistic Spectroscopic and Computational Characterization Evidencing the Preservation or Flipping of the Hydroxyl Group of 2-Phenylethyl Alcohol upon Single and Double Hydration

Even apparently simple, obtaining and analyzing observations on molecules and clusters and unambiguously assigning their structures is challenging. We report here the first ionization-loss Raman spectra compared to quantum chemical predictions for establishing the structural preferences of hydrates of the neurotransmitters hydroxy analogue, 2-phenylethyl alcohol (PEAL). The spectra encode two monohydrates and two previously unnoticed dihydrates, consequences of water insertion and sidewise attachment to the O-H group of gauche PEALs, in PEAL-H2O and PEAL-(H2O)2, or the higher-energy gauche-trans PEAL in the latter. The electronic structures retain the stable PEAL or flip its O-H to convert the gauche-trans PEAL conformer to the global minimum-energy dihydrate. We disclose conventional and bifurcated hydrogen bonds and electron steric repulsions by noncovalent interaction analysis and correlations between the experimental O-H stretching vibrational frequencies and the O-H and H···X bond lengths and electron densities, pointing to implications on hydrate forms and our approach virtue.

Molecular insights into the inhibitory potential of anthocyanidins on glucokinase regulatory protein

Computational methods were used to investigate six anthocyanidins exhibiting antidiabetic activity by inhibiting glucokinase regulatory protein (GKRP) activity. Density functional theory was used to optimise the geometry of anthocyanidins and calculate their quantum chemical properties. A blind docking method was employed to conduct a molecular docking study, which revealed that delphinidin (Del), cyanidin (Cya), and pelargonidin (Pel) as potential GKRP inhibitors with the lowest binding free energy of -8.7, -8.6, and -8.6 kcal/mol, corresponding to high binding affinity. The molecular dynamics study further verified the blind docking results by showing high GKRP-F1P complex stability and high binding affinity calculated through the MM/GBSA method, upon the binding of pelargonidin. The lower RMSF values of pivotal GK-interacting residues for GKRP-F1P-Pel compared to GKRP-F1P, as a positive control, indicating pelargonidin ability to maintain the inactive conformation of GKRP through the inhibition of GK binding. The key residues that control the binding of the F1P to GKRP and anthocyanidin to GKRP-F1P were also identified in this study. Altogether, pelargonidin is anthocyanidins-derived natural products that have the most potential to act as inhibitors of GKRP and as antidiabetic nutraceuticals.

Transition State Stabilization by SCF2 -H⋅⋅⋅O Bifurcated Hydrogen Bond

The difluoromethylthio group (SCF₂ H), which is generally considered a highly lipophilic weak hydrogen bonding donor, has attracted special interest from the pharmaceutical and agrochemical industry. Remarkably, there have been relatively few literature investigations of SCF₂ H hydrogen bonding interactions. Here, we report the determination of the hydrogen bond acidity parameter A of the SCF₂ H in the most popularly used electrophilic difluoromethylthiolating reagent. We present kinetic and computational evidence of the RSCF₂ -H⋅⋅⋅O2 bifurcated hydrogen bond for stabilizing the SCF₂ H-transferring transition state, which could cause a reversal of apparent electrophilic reactivity of difluoromethylthiolating and trifluoromethylthiolating reagents. Solvent effects on the RSCF₂ -H⋅⋅⋅O2 bifurcated hydrogen bonds will also be discussed.

Impacts of noncovalent interactions involving sulfur atoms on protein stability, structure, folding, and bioactivity

This review discusses the various types of noncovalent interactions in which sulfur atoms participate and their effects on protein stability, structure, folding and bioactivity. Current approaches and recommendations for modelling these noncovalent interactions (in terms of both geometries and interaction energies) are highlighted.

A proof-of-concept study of the secondary structure of influenza A, B M2 and MERS- and SARS-CoV E transmembrane peptides using folding molecular dynamics simulations in a membrane mimetic solvent

Here, we have carried out a proof-of-concept molecular dynamics (MD) simulation with adaptive tempering in a membrane mimetic environment to study the folding of single-pass membrane peptides. We tested the influenza A M2 viroporin, influenza B M2 viroporin, and protein E from coronaviruses MERS-Cov-2 and SARS-CoV-2 peptides with known experimental secondary structures in membrane bilayers. The two influenza-derived peptides are significantly different in the peptide sequence and secondary structure and more polar than the two coronavirus-derived peptides. Through a total of more than 50 μs of simulation time that could be accomplished in trifluoroethanol (TFE), as a membrane model, we characterized comparatively the folding behavior, helical stability, and helical propensity of these transmembrane peptides that match perfectly their experimental secondary structures, and we identified common motifs that reflect their quaternary organization and known (or not) biochemical function. We showed that BM2 is organized into two structurally distinct parts: a significantly more stable N-terminal half, and a fast-converting C-terminal half that continuously folds and unfolds between α-helical structures and non-canonical structures, which are mostly turns. In AM2, both the N-terminal half and C-terminal half are very flexible. In contrast, the two coronavirus-derived transmembrane peptides are much more stable and fast helix-formers when compared with the influenza ones. In particular, the SARS-derived peptide E appears to be the fastest and most stable helix-former of all the four viral peptides studied, with a helical structure that persists almost without disruption for the whole of its 10 μs simulation. By comparing the results with experimental observations, we benchmarked TFE in studying the conformation of membrane and hydrophobic peptides. This work provided accurate results suggesting a methodology to run long MD simulations and predict structural properties of biologically important membrane peptides.

The role of cholesterol binding in the control of cholesterol by the Scap-Insig system

Scap and Insig, two proteins embedded in the membrane of the endoplasmic reticulum (ER), regulate the synthesis of cholesterol in animal cells by forming a dimer in the presence of high concentrations of cholesterol. Cryo-electron microscopic structures for the Scap-Insig dimer show a sterol-binding site at the dimer interface, but none of the structures include cholesterol itself. Here, a molecular docking approach developed to characterise cholesterol binding to the transmembrane (TM) regions of membrane proteins is used to characterise cholesterol binding to sites on the TM surface of the dimer and to the interfacial binding site. Binding of cholesterol is also observed at sites on the extra-membranous luminal domains of Scap, but the properties of these sites suggest that they will be unoccupied in vivo. Comparing the structure of Scap in the dimer with that predicted by AlphaFold for monomeric Scap suggests that dimer formation could result in relocation of TM helix 7 of Scap and of the loop between TM6 and 7, and that this could be the key change on Scap that signals that there is a high concentration of cholesterol in the ER.

Understanding the role of non-Watson-Crick base pairs in DNA-protein recognition: Structural and energetic aspects using crystallographic database analysis and quantum chemical calculation

Specific recognition of DNA base sequences by proteins is vital for life-cycles of all organisms. In a large number of crystal structures of protein-DNA complexes, DNA conformation significantly deviates from the canonical B-DNA structure. A key question is whether such alternate conformations exist prior to protein binding and one is selected for complexation or the structure observed is induced by protein binding. Non-canonical base pairs, such as Hoogsteen base pairs, are often observed in crystal structures of protein-DNA complexes. We decided to explore whether the occurrence of such non-canonical base pairs in protein-DNA complexes is induced by the protein or is selected from pre-existing conformations. Detailed quantum chemical calculations with dispersion-corrected density functional theory (DFT-D) indicated that most of the non-canonical base pairs with DNA bases are stable even in the absence of the interacting amino acids. However, the G:G Hoogsteen base pair, which also appears in the telomere structure, appears to be unstable in the absence of other stabilizing agents, such as positively charged amino acids. Thus, the stability of many of the non-canonical base pair containing duplexes may be close to the canonical B-DNA structure and hence energetically accessible in the ground state; suggesting that the selection from pre-existing conformations may be an important mechanism for observed non-canonical base pairs in protein-DNA complexes.

Revisiting the concept of peptide bond planarity in an iron-sulfur protein by neutron structure analysis

The planarity of the peptide bond is important for the stability and structure formation of proteins. However, substantial distortion of peptide bonds has been reported in several high-resolution structures and computational analyses. To investigate the peptide bond planarity, including hydrogen atoms, we report a 1.2-Å resolution neutron structure of the oxidized form of high-potential iron-sulfur protein. This high-resolution neutron structure shows that the nucleus positions of the amide protons deviate from the peptide plane and shift toward the acceptors. The planarity of the H─N─C═O plane depends strongly on the pyramidalization of the nitrogen atom. Moreover, the orientation of the amide proton of Cys⁷⁵ is different in the reduced and oxidized states, possibly because of the electron storage capacity of the iron-sulfur cluster.

Carbamoylated Guanidine-Containing Polymers for Non-Covalent Functional Protein Delivery in Serum-Containing Media

Despite the high potential of controlling cellular processes and treating various diseases by intracellularly delivered proteins, current delivery systems exhibit poor efficiency due to poor serum stability, cellular entry, and cytosolic availability of proteins. Here, we report a novel functional group, phenyl carbamoylated guanidine (Ph-CG), that greatly enhances the delivery efficiency to various types of cells. Owing to the substantially lowered pK_a , the hydrophobic Ph-CG offers optimized inter-macromolecular interactions via enhanced hydrogen-bonding and hydrophobic interactions. The coplanarity of Ph-CG also leads to the better intracellular entry of protein complexes. Intracellularly delivered apoptosis-inducing enzymes and antibodies significantly induce cell viability inhibitions in a serum-containing medium. The newly developed Ph-CG can be introduced to various existing carriers, leading to the realization of future therapeutic protein delivery.

Exciton interactions in helical crystals of a hydrogen-bonded eumelanin monomer

Eumelanin, a naturally occurring group of heterogeneous polymers/aggregates providing photoprotection to living organisms, consist of 5,6-dihydroxyindole (DHI) and 5,6-dihydroxyindole-2-carboxylic acid (DHICA) building blocks. Despite their prevalence in the animal world, the structure and therefore the mechanism behind the photoprotective broadband absorption and non-radiative decay of eumelanin remain largely unknown. As a small step towards solving the incessant mystery, DHI is crystallized in a non-protic solvent environment to obtain DHI crystals having a helical packing motif. The present approach reflects the solitary directional effect of hydrogen bonds between the DHI chromophores for generating the crystalline assembly and filters out any involvement of the surrounding solvent environment. The DHI single crystals having an atypical chiral packing motif (P212121 Sohncke space group) incorporate enantiomeric zig-zag helical stacks arranged in a herringbone fashion with respect to each other. Each of the zig-zag helical stacks originates from a bifurcated hydrogen bonding interaction between the hydroxyl substituents in adjacent DHI chromophores which act as the backbone structure for the helical assembly. Fragment-based excited state analysis performed on the DHI crystalline assembly demonstrates exciton delocalization along the DHI units that connect each enantiomeric helical stack while, within each stack, the excitons remain localized. Fascinatingly, over the time evolution for generation of single-crystals of the DHI-monomer, mesoscopic double-helical crystals are formed, possibly attributed to the presence of covalently connected DHI trimers in chloroform solution. The oligomeric DHI (in line with the chemical disorder model) along with the characteristic crystalline packing observed for DHI provides insights into the broadband absorption feature exhibited by the chromophore.

Folded network and structural transition in molten tin

The fundamental relationships between the structure and properties of liquids are far from being well understood. For instance, the structural origins of many liquid anomalies still remain unclear, but liquid-liquid transitions (LLT) are believed to hold a key. However, experimental demonstrations of LLTs have been rather challenging. Here, we report experimental and theoretical evidence of a second-order-like LLT in molten tin, one which favors a percolating covalent bond network at high temperatures. The observed structural transition originates from the fluctuating metallic/covalent behavior of atomic bonding, and consequently a new paradigm of liquid structure emerges. The liquid structure, described in the form of a folded network, bridges two well-established structural models for disordered systems, i.e., the random packing of hard-spheres and a continuous random network, offering a large structural midground for liquids and glasses. Our findings provide an unparalleled physical picture of the atomic arrangement for a plethora of liquids, shedding light on the thermodynamic and dynamic anomalies of liquids but also entailing far-reaching implications for studying liquid polyamorphism and dynamical transitions in liquids.

Unraveling the structural origin of liquid anomalies remains a challenging topic. Xu et al. propose a folded-network structural model for molten tin and provide insights into the observed second-order-like structural transition.

On cellulose spatial organization and interactions as unraveled by diffraction and spectroscopic methods throughout the 20th century

This contribution attempts to describe the path towards determination of cellulose crystal structure down to atomic coordinates, towards the determination of its molecular conformation, as well as towards the details of the intricate pattern of hydrogen bonds and their dynamics. This path started at the beginning of the 20th century with X-ray diffraction, continued with electron diffraction, infrared and Raman spectroscopy, and significant knowledge was gained by methods of NMR spectroscopy. Towards the end of the 20th century and at the beginning of the 21st century, X-ray diffraction in conjunction with neutron diffraction provided the position of hydrogens, which led to detailed description of the geometry of hydrogen bonding network in cellulose. Quantum chemical and molecular dynamics calculations, polarized infrared spectroscopy and sum frequency generation vibrational spectroscopy were used to identify the origins of the vibrational modes in cellulose and to describe their extensive coupling mediated by hydrogen bonds. The role of amphiphilic character of cellulose macromolecule (and consequent hydrophobic interactions) in cellulose properties and behavior has been gaining more recognition in the 21st century.

An empirical model for solvation based on surface site interaction points

Surface site interaction points (SSIP) provide a quantitative description of the non-covalent interactions a molecule makes with the environment based on specific intermolecular contacts, such as H-bonds. Summation of the free energy of interaction of each SSIP across the surface of a molecule allows calculation of solvation energies and partition coefficients. A rule-based approach to the assignment of SSIPs based on chemical structure has been developed, and a combination of experimental data on the formation of 1 : 1 H-bonded complexes in non-polar solvents and partition of solutes between different solvents was used to parameterise the method. The resulting model is simple to implement using just a spreadsheet and accurately describes the transfer of a wide range of different solutes from water to a wide range of different organic solvents (overall rmsd is 1.4 kJ mol-1 for 1713 data points). The hydrophobic effect as well as the properties of perfluorocarbon solvents are described well by the model, and new descriptors have been determined for range of organic solvents that were not accessible by direct investigation of H-bond formation in non-polar solvents.

Charge-assisted bond and molecular self-assembly drive the gelation of lenvatinib mesylate

Lenvatinib mesylate (LM) is a first-line anticancer agent for the treatment of unresectable hepatocellular carcinoma, while it formed viscoelastic hydrogel when contacting with aqueous medium, which would significantly hinder its in vitro dissolution. The aim of this study was to systematicly explore the gelation mechanism and gel properties via thermal analysis, rheology, morphology and spectroscopy studies. The formed hydrogel was found to be composed of a new polymorph of crystalline LM, and its mechanical strength depended on the cross-linking degree of the fibrillar network structure. Spectroscopy analyses revealed that the intermolecular hydrogen bonds (the bifurcated hydrogen bond between the adjacent urea groups and the NH⋯OC hydrogen bond between the primary amide groups) as well as π-π stacking interactions (between the benzene ring and the quinoline ring) were suggested to be the driving forces for the self-assembly of LM during gelation process. Additionally, no gelation phenomenon was observed when suspending the base form lenvatinib in water, while it could form gel in various acidic solutions (e.g. hydrochloric acid, phosphoric acid and methanesulfonic acid) because the regenerated N⁺-H group increased the solubility of lenvatinib and promoted the balance between the dissolution or aggregation of LX (X: acid radical ion) molecules in solutions. In conclusion, the charge-assisted bond N⁺-H in LM molecule and intermolecular non-covalent interactions drived the hydrogel formation of LM in aqueous media. This study elucidates the gelation mechanism and gel properties of LM hydrogel, which would be helpful to figure out strategy to eliminate its gelation fundamentally and pave the way for its further formulation development in future.

Biomedical Applications of Supramolecular Materials in the Controllable Delivery of Steroids

Glucocorticoids are a class of steroid hormones secreted from the adrenal glands. The strong anti-inflammatory effects make it be one of the most popular and versatile drugs available to treat chronic inflammatory diseases. Additionally, supramolecular materials have been widely exploited in drug delivery, due to their biocompatibility, tunability, and predictability. Thus, steroid-based supramolecular materials and the release of steroids have been applied in the treatment of inflammatory diseases. This mini-review summarized recent advances in supramolecular materials loaded with glucocorticoid drugs in terms of hydrophobic interactions, electrostatic interactions, hydrogen bonding, and π-π stackings. We also discussed and prospected the application of the glucocorticoid drugs-based supramolecular system on chronic rhinosinusitis, multifactorial inflammatory disease of the nasal and paranasal sinuses mucosal membranes. Overall, supramolecular materials can provide an alternative to traditional materials as a novel delivery platform in clinical practice.

A Synthetic Galectin Mimic

Galectins are a galactoside specific subclass of carbohydrate binding proteins (lectins) involved in various cellular activities, certain cancers, infections, inflammations, and many other biological processes. The molecular basis for the selectivity of galectins is well-documented and revolves around appropriate interaction complementarity: an aromatic residue for C-H⋅⋅⋅π interactions and polar residues for (charge assisted) hydrogen bonds with the axial hydroxyl group of a galactoside. However, no synthetic mimics are currently available. We now report on the design and synthesis of the first galectin mimic (6), and show that it has a higher than 65-fold preference for n-octyl-β-galactoside (8) over n-octyl-β-glucoside (7) in CD2 Cl2 containing 5 % [D6 ]DMSO (with Ka ≥4500 M-1 for 6:8). Molecular modeling informed by nOe studies reveal a high degree of interaction complementarity between 6 and galactoside 8, which is very similar to the interaction complementarity found in natural galectins.

Interfacial binding sites for cholesterol on GABAA receptors and competition with neurosteroids

γ-Aminobutyric acid type A (GABA_A) receptors in the brain are located in the outer membranes of brain cells where the concentration of cholesterol is high. Of the 25 available high-resolution structures available for GABA_A receptors, none were determined in the presence of cholesterol, but four include resolved molecules of cholesterol hemisuccinate (CHS). Here, a molecular docking procedure is used to sweep the transmembrane (TM) surfaces of the receptors for cholesterol binding sites. Cholesterol docking poses determined in this way match 89% of the resolved CHS when CHS molecules deemed unlikely to represent typical bound cholesterols are excluded. The receptors are pentameric, and their TM surfaces consist of a set of five facets, each including pairs of TM helices from two adjacent subunits. Each facet contains hydrophobic hollows running from one side of the membrane to the other, within which are six potential binding sites for cholesterol, three on each side of the membrane. High-resolution structures of GABA_A receptors with bound neurosteroids show that neurosteroids bind in these cholesterol binding sites, so the binding of neurosteroids and cholesterol will be competitive.

The First Insight Into the Supramolecular System of D,L-α-Difluoromethylornithine: A New Antiviral Perspective

Targeting the polyamine biosynthetic pathway by inhibiting ornithine decarboxylase (ODC) is a powerful approach in the fight against diverse viruses, including SARS-CoV-2. Difluoromethylornithine (DFMO, eflornithine) is the best-known inhibitor of ODC and a broad-spectrum, unique therapeutical agent. Nevertheless, its pharmacokinetic profile is not perfect, especially when large doses are required in antiviral treatment. This article presents a holistic study focusing on the molecular and supramolecular structure of DFMO and the design of its analogues toward the development of safer and more effective formulations. In this context, we provide the first deep insight into the supramolecular system of DFMO supplemented by a comprehensive, qualitative and quantitative survey of non-covalent interactions via Hirshfeld surface, molecular electrostatic potential, enrichment ratio and energy frameworks analysis visualizing 3-D topology of interactions in order to understand the differences in the cooperativity of interactions involved in the formation of either basic or large synthons (Long-range Synthon Aufbau Modules, LSAM) at the subsequent levels of well-organized supramolecular self-assembly, in comparison with the ornithine structure. In the light of the drug discovery, supramolecular studies of amino acids, essential constituents of proteins, are of prime importance. In brief, the same amino-carboxy synthons are observed in the bio-system containing DFMO. DFT calculations revealed that the biological environment changes the molecular structure of DFMO only slightly. The ADMET profile of structural modifications of DFMO and optimization of its analogue as a new promising drug via molecular docking are discussed in detail.

Through space JFH spin-spin coupling constant transmission pathways in 2-(trifluoromethyl)thiophenol: formation of unusual stabilizing bifurcated CFHS and CFSH interactions

Given its importance and the possibility of organic F to participate in hydrogen bonds (H-bonds), the understanding of its behavior as a H-bond acceptor with different donors is crucial. The interest in organofluorine compounds and the works related to the study of the participation of this atom in non-covalent interactions is constantly growing. Following recent studies in this subject, we evaluated the existence of two bifurcated intramolecular interactions, a bifurcated CFHS H-bond in the cis conformer of 2-trifluoromethylthiophenol and an unusual, bifurcated CFSH interaction in the trans conformer. The JFH spin-spin coupling constant (SSCC) was evaluated for 2-trifluoromethylthiophenol both experimentally by 1H and 19F NMR and theoretically using the natural bond orbitals (NBO), the quantum theory of atoms in molecules (QTAIM) and the non-covalent interactions (NCI) framework. Although both interactions are crucial for the stabilization of the conformer geometries, the observed positive JFH spin-spin coupling constant (SSCC) is mainly resultant from the trans conformer, which has a large calculated positive SSCC, and is transmitted through steric interactions involving the F lone pairs and the σSH bonding orbital.

Competing HB acceptors: an extensive NMR investigations corroborated by single crystal XRD and DFT calculations

A series of N-benzoylanthranilamide derivatives have been synthesized with the substitution of competitive HB acceptors and investigated by NMR spectroscopy and single crystal XRD. The interesting rivalry for HB acceptance between [double bond splayed left]C[double bond, length as m-dash]O and X (F or OMe) is observed in the investigated molecules which leads to an unusual increase in the electron density at the site of one of the NH protons, reflecting in the high field resonance in the 1H NMR spectrum. The NMR experimental findings and single crystal XRD are further reinforced by the DFT studies.

Base Pairing and Functional Insights into N3-Methylcytidine (m3C) in RNA

N³-methylcytidine (m³C) is present in both eukaryotic tRNA and mRNA and plays critical roles in many biological processes. We report the synthesis of the m³C phosphoramidite building block and its containing RNA oligonucleotides. The base-pairing stability and specificity studies show that the m³C modification significantly disrupts the stability of the Watson-Crick C:G pair. Further m³C decreases the base pairing discrimination between C:G and the other mismatched C:A, C:U, and C:C pairs. Our molecular dynamic simulation study further reveals the detailed structural insights into the m³C:G base pairing pattern in an RNA duplex. More importantly, the biochemical investigation of m³C using reverse transcription in vitro shows that N³-methylation specifies the C:A pair and induces a G to A change using HIV-1-RT, MMLV-RT, and MutiScribe-RT enzymes, all with relatively low replication fidelity. For other reverse transcriptases with higher fidelity like AMV-RT, the methylation could completely shut down DNA synthesis. Our work provides detailed insights into the thermostability of m³C in RNA and a foundation for developing new molecular tools for mapping m³C in different RNA contexts and exploring the biochemical and biomedical potentials of m³C in the design and development of RNA based therapeutics.

(Me2NH2)10[H2-Dodecatungstate] polymorphs: dodecatungstate cages embedded in a variable dimethylammonium cation + water of crystallization matrix

Two polymorphs and a solvatomorph of a new dimethylammonium polytungstate-decakis(dimethylammonium) dihydrogendodecatungstate, (Me₂NH₂)₁₀(W₁₂O₄₂)·nH₂O (n = 10 or 11)-have been synthesized. Their structures were characterized by single-crystal X-ray diffraction and solid-phase NMR methods. The shape of the dodecatungstate anions is essentially the same in all three structures, their interaction with the cations and water of crystallization, however, is remarkably variable, because the latter forms different hydrogen-bonded networks, and provides a highly versatile matrix. Accordingly, the N-H⋯O and C-H⋯O hydrogen bonds are positioned in each crystal lattice in a variety of environments, characteristic to the structure, which can be distinguished by solid-state ¹H-CRAMPS, ¹³C, ¹⁵N CP MAS and ¹H-¹³C heteronuclear correlation NMR. Thermogravimetry of the solvatomorphs also reflect the difference and multiformity of the environment of the water molecules in the different crystal lattices. The major factors behind the variability of the matrix are the ability of ammonium cations to form two hydrogen bonds and the rigidity of the polyoxometalate anion cage. The positions of the oxygen atoms in the latter are favourable for the formation of bifurcated and trifurcated cation-anion hydrogen bonds, some which are so durable that they persist after the crystals are dissolved in water, forming ion associates even in dilute solutions. The H atom involved in furcated hydrogen bonds cannot be exchanged by deuterium when the compound is dissolved in D₂O. An obvious consequence of the versatility of the matrix is the propensity of these compounds to form multiple polymorphs.

Membrane-Mediated Lateral Interactions Regulate the Lifetime of Gramicidin Channels

The lipid matrix of cellular membranes is an elastic liquid crystalline medium. Its deformations regulate the functionality and interactions of membrane proteins,f membrane-bound peptides, lipid and protein-lipid domains. Gramicidin A (gA) is a peptide, which incorporates into membrane leaflets as a monomer and may form a transmembrane dimer. In both configurations, gA deforms the membrane. The transmembrane dimer of gA is a cation-selective ion channel. Its electrical response strongly depends on the elastic properties of the membrane. The gA monomer and dimer deform the membrane differently; therefore, the elastic energy contributes to the activation barriers of the dimerization and dissociation of the conducting state. It is shown experimentally that channel characteristics alter if gA molecules have been located in the vicinity of the conducting dimer. Here, based on the theory of elasticity of lipid membranes, we developed a quantitative theoretical model which allows explaining experimentally observed phenomena under conditions of high surface density of gA or its analogues, i.e., in the regime of strong lateral interactions of gA molecules, mediated by elastic deformations of the membrane. The model would be useful for the analysis and prediction of the gA electrical response in various experimental conditions. This potentially widens the possible applications of gA as a convenient molecular sensor of membrane elasticity.

Transmembrane Polar Relay Drives the Allosteric Regulation for ABCG5/G8 Sterol Transporter

The heterodimeric ATP-binding cassette (ABC) sterol transporter, ABCG5/G8, is responsible for the biliary and transintestinal secretion of cholesterol and dietary plant sterols. Missense mutations of ABCG5/G8 can cause sitosterolemia, a loss-of-function disorder characterized by plant sterol accumulation and premature atherosclerosis. A new molecular framework was recently established by a crystal structure of human ABCG5/G8 and reveals a network of polar and charged amino acids in the core of the transmembrane domains, namely, a polar relay. In this study, we utilize genetic variants to dissect the mechanistic role of this transmembrane polar relay in controlling ABCG5/G8 function. We demonstrated a sterol-coupled ATPase activity of ABCG5/G8 by cholesteryl hemisuccinate (CHS), a relatively water-soluble cholesterol memetic, and characterized CHS-coupled ATPase activity of three loss-of-function missense variants, R543S, E146Q, and A540F, which are respectively within, in contact with, and distant from the polar relay. The results established an in vitro phenotype of the loss-of-function and missense mutations of ABCG5/G8, showing significantly impaired ATPase activity and loss of energy sufficient to weaken the signal transmission from the transmembrane domains. Our data provide a biochemical evidence underlying the importance of the polar relay and its network in regulating the catalytic activity of ABCG5/G8 sterol transporter.

Impact of Cholesterol Concentration and Lipid Phase on Structure and Fluctuation of Amyloid Precursor Protein

Elevated levels of cellular cholesterol have been identified as one factor contributing to the onset of Alzheimer's disease (AD). Specific interaction between cholesterol and the amyloid precursor protein (APP), investigated via NMR experiments and computational studies, has been proposed to play a critical role in the processing of APP by secretases and the biogenesis of amyloid-β (Aβ) protein. We present all-atom molecular dynamics simulations of the 40-residue congener of the C-terminal domain of APP, C9916-55 (C99), in cholesterol-enriched DMPC lipid bilayers. We investigated the effect of cholesterol concentration on the conformational ensemble of wild-type C99 and C99-cholesterol associations at the low pH of endosomal environments, at which residues E22 and D23 are neutral. C99 was also characterized in liquid ordered domains for Dutch (E22Q) and Iowa (D23N) Familial AD mutants at low pH and for the wild-type sequence using protonation states characteristic of neutral pH. Our results reproduce the equilibrium constant of past NMR characterizations of the C99-cholesterol interaction but are not consistent with the C99-cholesterol binding hypothesis. We find that the lifetimes of both DMPC and cholesterol complexed with C99 display a power-law distribution of residence lifetimes. Longer-lived C99-DMPC and C99-cholesterol complexes are primarily stabilized by salt bridges and hydrogen bonds of lysine amines to phosphate and hydroxyl groups. Nevertheless, specific interfaces for C99-cholesterol association which are not present for DMPC can be identified. Changes to C99-cholesterol interfaces are found to depend on C99 tilt angle and orientation of the juxtamembrane domain of C99 containing residues E22 and D23. These observations support a more nuanced view of the C99-cholesterol interaction than has previously been suggested. We propose that cholesterol modulates the conformation and activity of C99 and other small transmembrane proteins indirectly through induction of the liquid ordered phase and directly through hydrogen bonding. This suggests a critical role for membrane heterogeneity introduced by cholesterol in modulating the structural ensemble of C99 and the production of Aβ.

The balance between side-chain and backbone-driven association in folding of the α-helical influenza A transmembrane peptide

The correct balance between attractive, repulsive and peptide hydrogen bonding interactions must be attained for proteins to fold correctly. To investigate these important contributors, we sought a comparison of the folding between two 25-residues peptides, the influenza A M2 protein transmembrane domain (M2TM) and the 25-Ala (Ala₂₅ ). M2TM forms a stable α-helix as is shown by circular dichroism (CD) experiments. Molecular dynamics (MD) simulations with adaptive tempering show that M2TM monomer is more dynamic in nature and quickly interconverts between an ensemble of various α-helical structures, and less frequently turns and coils, compared to one α-helix for Ala₂₅ . DFT calculations suggest that folding from the extended structure to the α-helical structure is favored for M2TM compared with Ala₂₅ . This is due to CH⋯O attractive interactions which favor folding to the M2TM α-helix, and cannot be described accurately with a force field. Using natural bond orbital (NBO) analysis and quantum theory atoms in molecules (QTAIM) calculations, 26 CH⋯O interactions and 22 NH⋯O hydrogen bonds are calculated for M2TM. The calculations show that CH⋯O hydrogen bonds, although individually weaker, have a cumulative effect that cannot be ignored and may contribute as much as half of the total hydrogen bonding energy, when compared to NH⋯O, to the stabilization of the α-helix in M2TM. Further, a strengthening of NH⋯O hydrogen bonding interactions is calculated for M2TM compared to Ala₂₅ . Additionally, these weak CH⋯O interactions can dissociate and associate easily leading to the ensemble of folded structures for M2TM observed in folding MD simulations.

Structural modelling of the lumenal domain of human GPAA1, the metallo-peptide synthetase subunit of the transamidase complex, reveals zinc-binding mode and two flaps surrounding the active site

Background

The transamidase complex is a molecular machine in the endoplasmic reticulum of eukaryotes that attaches a glycosylphosphatidylinositol (GPI) lipid anchor to substrate proteins after cleaving a C-terminal propeptide with a defined sequence signal. Its five subunits are very hydrophobic; thus, solubility, heterologous expression and complex reconstruction are difficult. Therefore, theoretical approaches are currently the main source of insight into details of 3D structure and of the catalytic process.

Results

In this work, we generated model 3D structures of the lumenal domain of human GPAA1, the M28-type metallo-peptide-synthetase subunit of the transamidase, including zinc ion and model substrate positions. In comparative molecular dynamics (MD) simulations of M28-type structures and our GPAA1 models, we estimated the metal ion binding energies with evolutionary conserved amino acid residues in the catalytic cleft. We find that canonical zinc binding sites 2 and 3 are strongest binders for Zn1 and, where a second zinc is available, sites 2 and 4 for Zn2. Zinc interaction of site 5 with Zn1 enhances upon substrate binding in structures with only one zinc. Whereas a previously studied glutaminyl cyclase structure, the best known homologue to GPAA1, binds only one zinc ion at the catalytic site, GPAA1 can sterically accommodate two. The M28-type metallopeptidases segregate into two independent branches with regard to one/two zinc ion binding modality in a phylogenetic tree where the GPAA1 family is closer to the joint origin of both groups. For GPAA1 models, MD studies revealed two large loops (flaps) surrounding the active site being involved in an anti-correlated, breathing-like dynamics.

Conclusions

In the light of combined sequence-analytic and phylogenetic arguments as well as 3D structural modelling results, GPAA1 is most likely a single zinc ion metallopeptidase. Two large flaps environ the catalytic site restricting access to large substrates.

Reviewers

This article was reviewed by Thomas Dandekar (MD) and Michael Gromiha.

Quantitative Analysis of Multiplex H-Bonds

H-bonding is the predominant geometrical determinant of biomolecular structure and interactions. As such, considerable analyses have been undertaken to study its detailed energetics. The focus, however, has been mostly reserved for H-bonds comprising a single donor and a single acceptor. Herein, we measure the prevalence and energetics of multiplex H-bonds that are formed between three or more groups. We show that 92% of all transmembrane helices have at least one non-canonical H-bond formed by a serine or threonine residue whose hydroxyl side chain H-bonds to an over-coordinated carbonyl oxygen at position i-4, i-3, or i in the sequence. Isotope-edited FTIR spectroscopy, coupled with DFT calculations, enables us to determine the bond enthalpies, pointing to values that are up to 127% higher than that of a single canonical H-bond. We propose that these strong H-bonds serve to stabilize serine and threonine residues in hydrophobic environments while concomitantly providing them flexibility between different configurations, which may be necessary for function.

A detailed look at the bonding interactions in the microsolvation of monoatomic cations

Global and local descriptors of the properties of intermolecular bonding, formally derived from independent methodologies (QTAIM, NCI, NBO, density differences) afford a highly complex picture of the bonding interactions responsible for microsolvation of monoatomic cations. In all cases, the dominant factor dictating geometries and interaction strengths is the electrophilic power of the metal cation. The formal charge disrupts the hydrogen bonding network otherwise present in pristine water clusters, making the hydrogen bonds considerably stronger, even inducing some degree of covalency. All MO interactions are highly ionic, with strengths than in some cases approach that of the reference LiCl bond. Accumulation of electron density in the region connecting MO is observed, thus, ionic bonding in the microsolvation of monoatomic cations is not as simple as an electrostatic interaction between opposing charges.