Cation–π interactions: Complexes of guanidinium and simple aromatic systems

Fast single metal cation conduction in ion-water aggregated aqueous battery electrolytes

Metal ion transport in solution is closely linked to its interactions with counter anions and solvent molecules. This interplay creates a longstanding trade-off between transference number (tMn⁺), ionic conductivity (δ), and solvation process. Advanced aqueous batteries with metal negative electrode require electrolytes with unity tMn⁺, high δ and low solvation energy. Here we introduce guanidinium sulfate (Gdm2SO4) into metal sulfate aqueous solutions to construct the ion-water aggregated electrolytes. These electrolytes exhibit fast single ion conduction, approaching unity tMn⁺ and high δ over 50 mS cm-1 for various metal cations (M= Zn, Cu, Fe, Sn and Li). The ion-water aggregates, dynamically formed by strong hydrogen bonding between sulfate anions, guanidinium cations and water, featuring an unfrustrated topological structure to suppress both anion mobility and water activity. This general configuration decouples the metal charge carrier from its coordination sheath, resulting in decreased solvation energy. These merits lead to homogeneous metal plating/stripping behavior with high coulombic efficiency of 99.9%. Moreover, the ion-water aggregates with reinforced kosmotropic characteristics significantly decrease the freezing point of the sulfate-based electrolytes to -28 oC, making them widely applicable in aqueous metal batteries for both intercalation and conversion positive electrodes.

Cryo-EM reconstruction of oleate hydratase bound to a phospholipid membrane bilayer

Oleate hydratase (OhyA) is a bacterial peripheral membrane protein that catalyzes FAD-dependent water addition to membrane bilayer-embedded unsaturated fatty acids. The opportunistic pathogen Staphylococcus aureus uses OhyA to counteract the innate immune system and support colonization. Many Gram-positive and Gram-negative bacteria in the microbiome also encode OhyA. OhyA is a dimeric flavoenzyme whose carboxy terminus is identified as the membrane binding domain; however, understanding how OhyA binds to cellular membranes is not complete until the membrane-bound structure has been elucidated. All available OhyA structures depict the solution state of the protein outside its functional environment. Here, we employ liposomes to solve the cryo-electron microscopy structure of the functional unit: the OhyA•membrane complex. The protein maintains its structure upon membrane binding and slightly alters the curvature of the liposome surface. OhyA preferentially associates with 20-30 nm liposomes with multiple copies of OhyA dimers assembling on the liposome surface resulting in the formation of higher-order oligomers. Dimer assembly is cooperative and extends along a formed ridge of the liposome. We also solved an OhyA dimer of dimers structure that recapitulates the intermolecular interactions that stabilize the dimer assembly on the membrane bilayer as well as the crystal contacts in the lattice of the OhyA crystal structure. Our work enables visualization of the molecular trajectory of membrane binding for this important interfacial enzyme.

Excited States of apo-Guanidine-III Riboswitch Contribute to Guanidinium Binding through Both Conformational and Induced-Fit Mechanisms

Riboswitches are mRNA segments that regulate gene expression through conformational changes driven by their cognate ligand binding. The ykkC motif forms a riboswitch class that selectively senses a guanidinium ion (Gdm+) and regulates the downstream expression of proteins which aid in the efflux of excess Gdm+ from the cells. The aptamer domain (AD) of the guanidine-III riboswitch forms an H-type pseudoknot with a triple helical domain that binds a Gdm+. We studied the binding of Gdm+ to the AD of the guanidine (ykkC)-III riboswitch using computer simulations to probe the specificity of the riboswitch to Gdm+ binding. We show that Gdm+ binding is a fast process occurring on the nanosecond time scale, with minimal conformational changes to the AD. Using machine learning and Markov-state models, we identified the excited conformational states of the AD, which have a high Gdm+ binding propensity, making the Gdm+ binding landscape complex exhibiting both conformational selection and induced-fit mechanisms. The proposed apo-AD excited states and their role in the ligand-sensing mechanism are amenable to experimental verification. Further, targeting these excited-state conformations in discovering new antibiotics can be explored.

Synthesis and Exploitation of the Biological Profile of Novel Guanidino Xylofuranose Derivatives

The synthesis and biological evaluation of novel guanidino sugars as isonucleoside analogs is described. 5-Guanidino xylofuranoses containing 3- O -saturated/unsaturated hydrocarbon or aromatic-containing moieties were accessed from 5-azido xylofuranoses via reduction followed by guanidinylation with N , N '-bis( tert -butoxycarbonyl)- N ''-triflylguanidine. Molecules comprising novel types of isonucleosidic structures including 5-guanidino 3- O -methyl-branched N -benzyltriazole isonucleosides and a guanidinomethyltriazole 3'- O -dodecyl xylofuranos-5'-yl isonucleoside were accessed. The guanidinomethyltriazole derivative and a 3- O -dodecyl ( N -Boc)guanidino xylofuranose were revealed as selective inhibitors of acetylcholinesterase ( K i = 22.87 and 7.49 µM, respectively). The latter also showed moderate antiproliferative effects in chronic myeloid leukemia (K562) and breast cancer (MCF-7) cells. An aminomethyltriazole 5'-isonucleoside was the most potent molecule with low micromolar GI 50 values in both cells (GI 50 = 6.33 μM, 8.45 μM), similar to that of the drug 5-fluorouracil in MCF-7 cells. Moreover, the most bioactive compounds showed low toxicity in human fibroblasts, further indicating their interest as promising lead molecules.

Theoretical Insights into the Role of Water Molecules in the Guanidinium-Based Protein Denaturation Process in Specific to Aromatic Amino Acids

Noncovalent interactions between the guanidinium cation (Gdm⁺) and aromatic amino acids (AAs) in the water molecules have been studied using quantum chemical calculation and molecular dynamics (MD) simulations. Our studies show that there are two different modes of interactions between Gdm⁺ and AAs with and without water molecules. It is observed that nonhydrated Gdm⁺ interacts with AAs through N-H···π interactions, whereas hydrated clusters of Gdm⁺ are stabilized by stacking interactions with the help of the water-mediated hydrogen bond. Thus, different hydration patterns have significant effects on the predominant cation···π interactions in AAs-Gdm⁺ complexes. Findings from MD simulation elicit that the interaction pattern of Gdm⁺ with AAs varies as Phe < Tyr < Trp. Both the QM and MD calculations show a similar trend in the interaction of AAs with Gdm⁺. Moreover, the interaction of AAs with Gdm⁺ depends on the spatial orientation of AAs in the protein and the concomitant local structure, that is, the AAs present in the unstructured region of protein such as coils and bends exhibit higher binding for Gdm⁺ when compared to the AAs present in the structured region of the protein such as the α-helix and the β-sheet. Our study clearly reveals that H-bonded water molecules and the hydration pattern of Gdm⁺ as well as the positional presence of these AAs in the protein structure context play determining roles in the denaturation of protein by the Gdm⁺ cation.

Parametrization of halogen bonds in the CHARMM general force field: Improved treatment of ligand-protein interactions

A halogen bond is a highly directional, non-covalent interaction between a halogen atom and another electronegative atom. It arises due to the formation of a small region of positive electrostatic potential opposite the covalent bond to the halogen, called the 'sigma hole.' Empirical force fields in which the electrostatic interactions are represented by atom-centered point charges cannot capture this effect because halogen atoms usually carry a negative charge and therefore interact unfavorably with other electronegative atoms. A strategy to overcome this problem is to attach a positively charged virtual particle to the halogen. In this work, we extend the additive CHARMM General Force Field (CGenFF) to include such interactions in model systems of phenyl-X, with X being Cl, Br or I including di- and trihalogenated species. The charges, Lennard-Jones parameters, and halogen-virtual particle distances were optimized to reproduce the orientation dependence of quantum mechanical interaction energies with water, acetone, and N-methylacetamide as well as experimental pure liquid properties and relative hydration free energies with respect to benzene. The resulting parameters were validated in molecular dynamics simulations on small-molecule crystals and on solvated protein-ligand complexes containing halogenated compounds. The inclusion of positive virtual sites leads to better agreement across experimental observables, including preservation of ligand binding poses as a direct result of the improved representation of halogen bonding.

A Study of π-π Stacking Interactions and Aromaticity in Polycyclic Aromatic Hydrocarbon/Nucleobase Complexes

We analysed the interactions and aromaticity electron-density delocalisation observed in π-π complexes between the phenalenyl radical and acenaphthylene, and the DNA and RNA nucleobases (adenine, guanine, cytosine, thymine and uracil). Interaction energies are obtained at the M06-2X/6-311++G(2df,p) computational level for gas phase and PCM-water conditions. For both the phenalenyl radical and acenaphthylene, the complexes formed with guanine are the most stable ones. Atoms in molecules and natural bond orbital results reveal weak π-π interactions between both interacting moieties, characterized by bond critical points between C⋅⋅⋅C and C⋅⋅⋅N atoms. Nucleus independent chemical shifts (NICS) indicate the retention of the aromatic character of the monomers in the outer region of the complex. The fluctuation indexes reveal a loss of electron delocalisation upon complexation for all cases except guanine complexes. Nevertheless, the interface region shows large negative NICS values, which is not associated with an increase of the aromaticity or electron-density delocalisation, but with magnetic couplings of both molecules, leading to an unrealistic description of the aromatic behaviour in that region.

A theoretical study of ternary indole-cation-anion complexes

The simultaneous interactions of an anion and a cation with a π system were investigated by MP2 and M06-2X theoretical calculations. Indole was chosen as a model π system for its relevance in biological environments. Two different orientations of the anion, interacting with the N-H and with the C-H groups of indole, were considered. The four cations (Na(+), NH4(+), C(NH2)3(+) and N(CH3)4(+)) and the four anions (Cl(-), NO3(-), HCOO(-) and BF4(-)) included in the study are of biological interest. The total interaction energy of the ternary complexes was calculated and separated into its two- and three-body components and all of them are further divided into their electrostatic, exchange, repulsion, polarization and dispersion contributions using the local molecular orbital-energy decomposition analysis (LMO-EDA) methodology. The binding energy of the indole-cation-anion complexes depends on both ions, with the cation having the strongest effect. The intense cation-anion attraction determines the geometric and energetic features in all ternary complexes. These structures, with both ions on the same side of the π system, show an anti-cooperative interaction. However, the interaction is not only determined by electrostatics, but also the polarization contribution is important. Specific interactions like the one established between the anion and the N-H group of indole or the proton transfer between an acidic cation and a basic anion play a significant role in the energetics and the structure of particular complexes. The presence of the polar solvent as modelled with the polarizable continuum model (PCM) does not seem to have a significant effect on the geometry of the ternary complexes, but drastically weakens the interaction energy. Also, the strength of the interaction is reduced at a faster rate when the anion is pushed away, compared to the results obtained in the gas phase. The combination of PCM with the addition of one water molecule indicates that the PCM method properly reproduces the main energetic and geometrical changes, even at the quantitative level, but the explicit hydration allows refining the solvent effect and detecting cases that do not follow the general trend.

Aromatic Amino Acids-Guanidinium Complexes through Cation-π Interactions

Continuing with our interest in the guanidinium group and the different interactions than can establish, we have carried out a theoretical study of the complexes formed by this cation and the aromatic amino acids (phenylalanine, histidine, tryptophan and tyrosine) using DFT methods and PCM-water solvation. Both hydrogen bonds and cation-π interactions have been found upon complexation. These interactions have been characterized by means of the analysis of the molecular electron density using the Atoms-in-Molecules approach as well as the orbital interactions using the Natural Bond Orbital methodology. Finally, the effect that the cation-π and hydrogen bond interactions exert on the aromaticity of the corresponding amino acids has been evaluated by calculating the theoretical NICS values, finding that the aromatic character was not heavily modified upon complexation.

Syntheses of Cyclic Guanidine-Containing Natural Products

Naturally occurring guanidine derivatives frequently display medicinally useful properties. Among them, the higher order pyrrole-imidazole alkaloids, the dragmacidins, the crambescidins/batzelladines, and the saxitoxins/tetradotoxins have stimulated the development of many new synthetic methods over the past decades. We provide here an overview of the syntheses of these cyclic guanidine-containing natural products.

Hydrogen bonding and π–π interactions in imidazolium-chloride ionic liquid clusters

The importance of 1° and 2° hydrogen-bonding and anion–π⁺ interactions for ionic liquid structuring.

On the interaction between the imidazolium cation and aromatic amino acids. A computational study

Phe, Tyr and Trp form parallel complexes with cation⋯π interactions. His complexes are the strongest, but without making contact with the aromatic cloud.

Significant strength of charged DNA-protein π-π interactions: a preliminary study of cytosine

The present work characterized the preferred gas-phase structure and optimum interaction energy of both parallel stacked and perpendicular T-shaped dimers between cytosine (C), as a representative nucleobase, and aspartic/glutamic acid (DE), aspartate/glutamate (DE(-)) or arginine (R(+)), using detailed M06-2X/6-31+G(d,p) potential energy surface scans as a function of the relative monomer orientation. Through comparison to previous literature on the π-π interactions between the DNA nucleobases and the aromatic amino acid residues, this work will allow for comparisons between DNA-protein interactions involving aromatic and acyclic R-side chains, as well as comparisons of the relative geometric dependence and magnitude of π-π (C:DE), πcation-π (C:R(+)), and πanion-π (C:DE(-)) interactions. Our results show that the preferred relative monomer orientation is highly dependent on the monomer composition and charge, and is dictated by electrostatic-driven interactions. More importantly, for the first time, we report that the π-π interactions between cytosine and (neutral) aspartic/glutamic acid are up to approximately -40 kJ mol(-1), while the πcation-π or πanion-π interactions between cytosine and arginine or aspartate/glutamate are up to approximately -90 and -99 kJ mol(-1), respectively. An extensive investigation of the effects of the computational methodology implemented, including comparisons to detailed CCSD(T)/CBS potential energy surfaces and interaction energies, supports the use of M06-2X, as well as ωB97X-D, to study DNA-protein π-π interactions of varying composition and charge. Most importantly, the CCSD(T)/CBS results verify the strong nature of these DNA-protein π-π interactions, as well as the unique nature of the πcation-π and πanion-π counterparts. Therefore, our results emphasize that a wide variety of different types of noncovalent interactions between both cyclic and acyclic π-containing components can significantly contribute to the stability of DNA-protein complexes and likely play a larger role in biology than currently accepted.

Pyridin-2-yl guanidine derivatives: conformational control induced by intramolecular hydrogen-bonding interactions

The synthesis and conformational analysis of a series of pyridin-2-yl guanidine derivatives using NMR, X-ray crystallography, and B3LYP/6-31+G** theoretical studies are reported. A remarkable difference was observed in the (1)H NMR spectra of the guanidinium salts as compared with their N,N'-di-Boc protected and neutral analogues. This difference corresponds to a 180° change in the dihedral angle between the guanidine/ium moiety and the pyridine ring in the salts as compared to the Boc-protected derivatives, a conclusion that was supported by theoretical studies, X-ray data, and NMR analysis. Moreover, our data sustain the existence of two intramolecular hydrogen-bonding systems: (i) between the pyridine N1 atom and the guanidinium protons in the salts and (ii) within the tert-butyl carbamate groups of the Boc-protected derivatives. To verify that the observed conformational control arises from these intramolecular interactions, a new series of N-Boc-N'-propyl-substituted pyridin-2-yl guanidines were also prepared and studied.