Aliphatic C−H/π Interactions: Methane−Benzene, Methane−Phenol, and Methane−Indole Complexes

Origins of the unphysical noncovalent interaction energy curves obtained with the 2011 and 2012 Minnesota density functionals

With the noncovalent interaction energy curves of the methane dimer [(CH4)2], we have clarified two different origins of the unphysical noncovalent interaction energy curves obtained with the Minnesota density functionals of M11-L, MN12-L, and MN12-SX. For the M11-L functional, the unphysical inflection point on the (CH4)2 interaction energy curve originates from the inclusion of the long-range exchange. As to the MN12-L and MN12-SX functionals, the lack of smoothness restraints results in unphysical inflection points on the corresponding (CH4)2 interaction energy curves. As a result, exchange functionals are as important as dispersion corrections for density functionals to map noncovalent interaction energy surfaces reasonably. Moreover, very highly parameterized functionals with smoothness restraints are suggested for investigating noncovalent interaction energy surfaces.

Analysing the effect caused by increasing the molecular volume in M1-AChR receptor agonists and antagonists: a structural and computational study

M1 muscarinic acetylcholine receptor (M1-AChR), a member of the G protein-coupled receptors (GPCR) family, plays a crucial role in learning and memory, making it an important drug target for Alzheimer's disease (AD) and schizophrenia. M1-AChR activation and deactivation have shown modifying effects in AD and PD preclinical models, respectively. However, understanding the pharmacology associated with M1-AChR activation or deactivation is complex, because of the low selectivity among muscarinic subtypes, hampering their therapeutic applications. In this regard, we constructed two quantitative structure-activity relationship (QSAR) models, one for M1-AChR agonists (total and partial), and the other for the antagonists. The binding mode of 59 structurally different compounds, including agonists and antagonists with experimental binding affinity values (pKi), were analyzed employing computational molecular docking over different structures of M1-AChR. Furthermore, we considered the interaction energy (Einter), the number of rotatable bonds (NRB), and lipophilicity (ilogP) for the construction of the QSAR model for agonists (R2 = 89.64, QLMO2 = 78, and Qext2 = 79.1). For the QSAR model of antagonists (R2 = 88.44, QLMO2 = 82, and Qext2 = 78.1) we considered the Einter, the fraction of sp3 carbons fCsp3, and lipophilicity (MlogP). Our results suggest that the ligand volume is a determinant to establish its biological activity (agonist or antagonist), causing changes in binding energy, and determining the affinity for M1-AChR.

Protein-Ligand CH-π Interactions: Structural Informatics, Energy Function Development, and Docking Implementation

Here, we develop an empirical energy function based on quantum mechanical data for the interaction between methane and benzene that captures the contribution from CH-π interactions. Such interactions are frequently observed in protein-ligand crystal structures, particularly for carbohydrate ligands, but have been hard to quantify due to the absence of a model for CH-π interactions in typical molecular mechanical force fields or docking scoring functions. The CH-π term was added to the AutoDock Vina (AD VINA) scoring function enabling its performance to be evaluated against a cohort of more than 1600 occurrences in 496 experimental structures of protein-ligand complexes. By employing a conformational grid search algorithm, inclusion of the CH-π term was shown to improve the prediction of the preferred orientation of flexible ligands in protein-binding sites and to enhance the detection of carbohydrate-binding sites that display CH-π interactions. Last but not least, this term was also shown to improve docking performance for the CASF-2016 benchmark set and a carbohydrate set.

Computation of CCSD(T)-Quality NMR Chemical Shifts via Δ-Machine Learning from DFT

NMR spectroscopy undoubtedly plays a central role in determining molecular structures across different chemical disciplines, and the accurate computational prediction of NMR parameters is highly desirable. In this work, a new Δ-machine learning approach is presented to correct DFT-computed NMR chemical shifts using input features from the calculation and in addition highly accurate reference data at the CCSD(T)/pcSseg-2 level of theory with a basis set extrapolation scheme. The model is trained on a data set containing 1000 optimized and geometrically distorted structures of small organic molecules comprising most elements of the first three periods and containing data for 7090 ¹H and 4230 ¹³C NMR chemical shifts. Applied to the PBE0/pcSseg-2 method, the mean absolute deviation (MAD) on the internal NMR shift test set is reduced by 81% for ¹H and 92% for ¹³C at virtually no additional computational cost. For 12 different DFT functional and basis set combinations, the MAD of the ML-corrected NMR shifts ranges from 0.021 to 0.039 ppm (¹H) and from 0.38 to 1.07 ppm (¹³C). Importantly, the new method consistently outperforms the simple and widely used linear regression correction technique. This behavior is reproduced on three different external benchmark sets, confirming the generality and robustness of the correction scheme, which can easily be applied in DFT-based spectral simulations.

Barrier properties of Nup98 FG phases ruled by FG motif identity and inter-FG spacer length

Nup98 FG repeat domains comprise hydrophobic FG motifs linked through uncharged spacers. FG motifs capture nuclear transport receptors (NTRs) during nuclear pore complex (NPC) passage, confer inter-repeat cohesion, and condense the domains into a selective phase with NPC-typical barrier properties. We show that shortening inter-FG spacers enhances cohesion, increases phase density, and tightens such barrier - all consistent with a sieve-like phase. Phase separation tolerates mutating the Nup98-typical GLFG motifs, provided domain-hydrophobicity remains preserved. NTR-entry, however, is sensitive to (certain) deviations from canonical FG motifs, suggesting co-evolutionary adaptation. Unexpectedly, we observed that arginines promote FG-phase-entry apparently also by hydrophobic interactions/ hydrogen-bonding and not just through cation-π interactions. Although incompatible with NTR·cargo complexes, a YG phase displays remarkable transport selectivity, particularly for engineered GFP^NTR-variants. GLFG to FSFG mutations make the FG phase hypercohesive, precluding NTR-entry. Extending spacers relaxes this hypercohesion. Thus, antagonism between cohesion and NTR·FG interactions is key to transport selectivity.

Cumulative Millisecond-Long Sampling for a Comprehensive Energetic Evaluation of Aqueous Ionic Liquid Effects on Amino Acid Interactions

The interactions of amino acid side-chains confer diverse energetic contributions and physical properties to a protein's stability and function. Various computational tools estimate the effect of changing a given amino acid on the protein's stability based on parametrized (free) energy functions. When parametrized for the prediction of protein stability in water, such energy functions can lead to suboptimal results for other solvents, such as ionic liquids (IL), aqueous ionic liquids (aIL), or salt solutions. However, to our knowledge, no comprehensive data are available describing the energetic effects of aIL on intramolecular protein interactions. Here, we present the most comprehensive set of potential of mean force (PMF) profiles of pairwise protein-residue interactions to date, covering 50 relevant interactions in water, the two biotechnologically relevant aIL [BMIM/Cl] and [BMIM/TfO], and [Na/Cl]. These results are based on a cumulated simulation time of >1 ms. aIL and salt ions can weaken, but also strengthen, specific residue interactions by more than 3 kcal mol^-1, depending on the residue pair, residue-residue configuration, participating ions, and concentration, necessitating considering such interactions specifically. These changes originate from a complex interplay of competitive or cooperative noncovalent ion-residue interactions, changes in solvent structural dynamics, or unspecific charge screening effects and occur at the contact distance but also at larger, solvent-separated distances. This data provide explanations at the atomistic and energetic levels for complex IL effects on protein stability and should help improve the prediction accuracies of computational tools that estimate protein stability based on (free) energy functions.

Dynamic Long-Range Interactions Influence Substrate Binding and Catalysis by Human Histidine Triad Nucleotide-Binding Proteins (HINTs), Key Regulators of Multiple Cellular Processes and Activators of Antiviral ProTides

Human histidine triad nucleotide-binding (hHINT) proteins catalyze nucleotide phosphoramidase and acyl-phosphatase reactions that are essential for the activation of antiviral proTides, such as Sofosbuvir and Remdesivir. hHINT1 and hHINT2 are highly homologous but exhibit disparate roles as regulators of opioid tolerance (hHINT1) and mitochondrial activity (hHINT2). NMR studies of hHINT1 reveal a pair of dynamic surface residues (Q62, E100), which gate a conserved water channel leading to the active site 13 Å away. hHINT2 crystal structures identify analogous residues (R99, D137) and water channel. hHINT1 Q62 variants significantly alter the steady-state k_cat and K_m for turnover of the fluorescent substrate (TpAd), while stopped-flow kinetics indicate that K_D also changes. hHINT2, like hHINT1, exhibits a burst phase of adenylation, monitored by fluorescent tryptamine release, prior to rate-limiting hydrolysis and nucleotide release. hHINT2 exhibits a much smaller burst-phase amplitude than hHINT1, which is further diminished in hHINT2 R99Q. Kinetic simulations suggest that amplitude variations can be accounted for by a variable fluorescent yield of the E·S complex from changes in the environment of bound TpAd. Isothermal titration calorimetry measurements of inhibitor binding show that these hHINT variants also alter the thermodynamic binding profile. We propose that these altered surface residues engender long-range dynamic changes that affect the orientation of bound ligands, altering the thermodynamic and kinetic characteristics of hHINT active site function. Thus, studies of the cellular roles and proTide activation potential by hHINTs should consider the importance of long-range interactions and possible protein binding surfaces far from the active site.

Enantioselective Synthesis of Quaternary Oxindoles: Desymmetrizing Staudinger-Aza-Wittig Reaction Enabled by a Bespoke HypPhos Oxide Catalyst

This paper describes a catalytic asymmetric Staudinger-aza-Wittig reaction of (o-azidoaryl)malonates, allowing access to chiral quaternary oxindoles through phosphine oxide catalysis. We designed a novel HypPhos oxide catalyst to enable the desymmetrizing Staudinger-aza-Wittig reaction through the PIII/PV═O redox cycle in the presence of a silane reductant and an IrI-based Lewis acid. The reaction occurs under mild conditions, with good functional group tolerance, a wide substrate scope, and excellent enantioselectivity. Density functional theory revealed that the enantioselectivity in the desymmetrizing reaction arose from the cooperative effects of the IrI species and the HypPhos catalyst. The utility of this methodology is demonstrated by the (formal) syntheses of seven alkaloid targets: (-)-gliocladin C, (-)-coerulescine, (-)-horsfiline, (+)-deoxyeseroline, (+)-esermethole, (+)-physostigmine, and (+)-physovenine.

Benzene, an Unexpected Binding Unit in Anion–π Recognition: The Critical Role of CH/π Interactions

We report high-level ab initio calculations (CCSD(T)(full)/CBS//SCS-RI-MP2(full)/aug-cc-pwCVTZ) that demonstrate the importance of cooperativity effects when Anion–π and CH/π interactions are simultaneously established with benzene as the π-system. In fact, most of the complexes exhibit high cooperativity energies that range from 17% to 25.3% of the total interaction energy, which is indicative of the strong influence of the CH/π on the Anion–π interaction and vice versa. Moreover, the symmetry-adapted perturbation theory (SAPT) partition scheme was used to study the different energy contributions to the interaction energies and to investigate the physical nature of the interplay between both interactions. Furthermore, the Atoms in Molecules (AIM) theory and the Non-Covalent Interaction (NCI) approach were used to analyze the two interactions further. Finally, a few examples from the Protein Data Bank (PDB) are shown. All results stress that the concurrent formation of both interactions may play an important role in biological systems due to the ubiquity of CH bonds, phenyl rings, and anions in biomolecules.

From groove binding to intercalation: unravelling the weak interactions and other factors modulating the modes of interaction between methylated phenanthroline-based drugs and duplex DNA

Several antitumor drugs base their cytotoxicity on their capacity to intercalate between base pairs of DNA. Nevertheless, it has been established that the mechanism of intercalation of drugs in DNA starts with the prior groove binding mode of interaction of the drug with DNA. Sometimes, for some kind of flat small molecules, groove binding does not produce any cytotoxic effect and the fast transition of such flat small molecules to the cytotoxic intercalation mode is desirable. This is the case of methylated phenanthroline (phen) derivatives, where, changes in the substitution in the position and number of methyl groups determine their capability as cytotoxic compounds and, therefore, it is a way for the modulation of cytotoxic effects. In this work, we studied this modulation by means of the interaction of the [Pt(en)(phen)]²⁺ complex and several derivatives by methylation of phen in different number and position and the d(GTCGAC)₂ DNA hexamer via groove binding using PM6-DH2 and DFT-D methods. The analysis of the geometries, electronic structure and energetics of the studied systems was compared to experimental works to gain insight into the relation structure-interaction for the studied systems with cytotoxicity. The trends are explained by means of the Non-Covalent Interaction (NCI) index, the Energy Decomposition Analysis (EDA) and solvation contributions. Our results are in agreement with the experiments, in which the methylation of position 4 of phen seems to favour the interaction via groove binding thus making the transition to the intercalation cytotoxic mode difficult. Looking at the NCI results, these interactions come not only from the CH/π and CH/n interactions of the methyl group in position 4 but also from the ethylenediamine (en) ligand, whose orientation in the Pt complex was found in such a way that it produces a high number of weak interactions with DNA, especially with the sugar and phosphate backbone.

Biotin's Lessons in Drug Design

At the heart of drug design is the discovery of molecules that bind with high affinity to their drug targets. Biotin forms the strongest known noncovalent ligand-protein interactions with avidin and streptavidin, achieving femtomolar and picomolar affinities, respectively. This is made even more exceptional because biotin achieves this with a meagre molecular weight of 240 Da. Surprisingly, the approaches by which biotin achieves this are not in the standard repertoire of current medicinal chemistry practice. Biotin's biggest lesson is the importance of nonclassical H-bonds in protein-ligand complexes. Most of biotin's affinity stems from its flexible valeric acid side chain that forms CH-π, CH-O, and classical H-bonds with the lipophilic region of the binding pocket. Biotin also utilizes an oxyanion hole, a sulfur-centered H-bond, and water solvation in the bound state to achieve its potency. The facets and advantages of biotin's approach to binding should be more widely adopted in drug design.

NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts

In this work, we present NENCI-2021, a benchmark database of ∼8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation-π and anion-π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion-π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (E_int) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The E_int values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) E_int value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent-these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)-based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation-π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the E_int values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive E_int values, and maximum errors of ±0.2-0.3 kcal/mol (i.e., ∼±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.

Local energy decomposition analysis and molecular properties of encapsulated methane in fullerene (CH4@C60)

Methane has been successfully encapsulated within cages of C₆₀ fullerene, which is an appropriate model system to study confinement effects. Its chemistry and physics are also relevant for theoretical model descriptions. Here we provide insights into intermolecular interactions and predicted spectroscopic responses of the CH₄@C₆₀ complex and compared them with results from other methods and with data from the literature. Local energy decomposition analysis (LED) within the domain-based local pair natural orbital coupled cluster singles, doubles, and perturbative triples (DLPNO-CCSD(T)) framework was used, and an efficient protocol for studies of endohedral complexes of fullerenes is proposed. This approach allowed us to assess energies in relation to electronic and geometric preparation, electrostatics, exchange, and London dispersion for the CH₄@C₆₀ endohedral complex. The calculated stabilization energy of CH₄ inside the C₆₀ fullerene was -13.5 kcal mol^-1 and its magnitude was significantly larger than the latent heat of evaporation of CH₄. Evaluation of vibrational frequencies and polarizabilities of the CH₄@C₆₀ complex revealed that the infrared (IR) and Raman bands of the endohedral CH₄ were essentially "silent" due to the dielectric screening effect of C₆₀, which acted as a molecular Faraday cage. Absorption spectra in the UV-vis domain and ionization potentials of C₆₀ and CH₄@C₆₀ were predicted. They were almost identical. The calculated ¹H/¹³C NMR shifts and spin-spin coupling constants were in very good agreement with experimental data. In addition, reference DLPNO-CCSD(T) interaction energies for complexes with noble gases (Ng@C₆₀; Ng = He, Ne, Ar, Kr) were calculated. The values were compared with those derived from supramolecular MP2/SCS-MP2 calculations and estimates with London-type formulas by Pyykkö and coworkers [Phys. Chem. Chem. Phys., 2010, 12, 6187-6203], and with values derived from DFT-based symmetry-adapted perturbation theory (DFT-SAPT) by Hesselmann & Korona [Phys. Chem. Chem. Phys., 2011, 13, 732-743]. Selected points at the potential energy surface of the endohedral He₂@C₆₀ trimer were considered. In contrast to previous theoretical attempts with the DFT/MP2/SCS-MP2/DFT-SAPT methods, our calculations at the DLPNO-CCSD(T) level of theory predicted the He₂@C₆₀ trimer to be thermodynamically stable, which is in agreement with experimental observations.

S. Azevedo D, Guzmán-Vargas A, Felipe C. Effect of Calcination Temperature and Chemical Composition of PAN-Derived Carbon Microfibers on N2, CO2, and CH4 Adsorption

This work investigates the interplay of carbonization temperature and the chemical composition of carbon microfibers (CMFs), and their impact on the equilibration time and adsorption of three molecules (N₂, CO₂, and CH₄). PAN derived CMFs were synthesized by electrospinning and calcined at three distinct temperatures (600, 700 and 800 °C), which led to samples with different textural and chemical properties assessed by FTIR, TGA/DTA, XRD, Raman, TEM, XPS, and N₂ adsorption. We examine why samples calcined at low/moderate temperatures (600 and 700 °C) show an open hysteresis loop in nitrogen adsorption/desorption isotherms at -196.15 °C. The equilibrium time in adsorption measurements is nearly the same for these samples, despite their distinct chemical compositions. Increasing the equilibrium time did not allow for the closure of the hysteresis loop, but by rising the analysis temperature this was achieved. By means of the isosteric enthalpy of adsorption measurements and ab initio calculations, adsorbent/adsorbate interactions for CO₂, CH₄ and N₂ were found to be inversely proportional to the temperature of carbonization of the samples (CMF-600 > CMF-700 > CMF-800). The enhancement of adsorbent/adsorbate interaction at lower carbonization temperatures is directly related to the presence of nitrogen and oxygen functional groups on the surface of CMFs. Nonetheless, a higher concentration of heteroatoms also causes: (i) a reduction in the adsorption capacity of CO₂ and CH₄ and (ii) open hysteresis loops in N₂ adsorption at cryogenic temperatures. Therefore, the calcination of PAN derived microfibers at temperatures above 800 °C is recommended, which results in materials with suitable micropore volume and a low content of surface heteroatoms, leading to high CO₂ uptake while keeping acceptable selectivity with regards to CH₄ and moderate adsorption enthalpies.

Switching Aromatic Character by Complexation: π to π* Change Seen in Molecular Rotation Spectra

A dominating F···π*_aromatic interaction is found to govern the benzaldehyde···tetrafluoromethane complex, as revealed by this rotational spectroscopic investigation. Secondary F···π*_-C=O- and C σ*_CF4···π_aromatic interactions also contribute to the stability of the observed isomer. Narrow splittings have been observed in the rotational spectrum originating from a 3-fold internal rotation of CF₄ above the aromatic moiety, and a corresponding V₃ barrier was determined to be 1.572(14) kJ mol^-1. This is the first rotational spectroscopic evidence in the literature implying that the aromatic π* antibonding orbital can be activated not only by electron-withdrawing substituents but also by complexation partners containing atoms with high electronegativity, like CF₄. The results emphasize the partner molecules' role to modulate the π electron structure and show a change in the orbital character (π or π*) when participating in the formation of noncovalent interactions.

Theoretical investigation of the effect of hydrogen bonding on the stereoselectivity of the Diels–Alder reaction

Here, we report the computational prediction of high exo selectivities in a series of Diels–Alder reactions with H-bonding interaction.

Self-adjusting binding pockets enhance H2 and CH4 adsorption in a uranium-based metal-organic framework

A new, air-stable, permanently porous uranium(iv) metal-organic framework U(bdc)2 (1, bdc2- = 1,4-benzenedicarboxylate) was synthesized and its H2 and CH4 adsorption properties were investigated. Low temperature adsorption isotherms confirm strong adsorption of both gases in the framework at low pressures. In situ gas-dosed neutron diffraction experiments with different D2 loadings revealed a rare example of cooperative framework contraction (ΔV = -7.8%), triggered by D2 adsorption at low pressures. This deformation creates two optimized binding pockets for hydrogen (Q st = -8.6 kJ mol-1) per pore, in agreement with H2 adsorption data. Analogous experiments with CD4 (Q st = -24.8 kJ mol-1) and N,N-dimethylformamide as guests revealed that the binding pockets in 1 adjust by selective framework contractions that are unique for each adsorbent, augmenting individual host-guest interactions. Our results suggest that the strategic combination of binding pockets and structural flexibility in metal-organic frameworks holds great potential for the development of new adsorbents with an enhanced substrate affinity.

Highly selective Diels-Alder and Heck arylation reactions in a divergent synthesis of isoindolo- and pyrrolo-fused polycyclic indoles from 2-formylpyrrole

A highly regio-, chemo- and stereoselective divergent synthesis of isoindolo- and pyrrolo-fused polycyclic indoles is herein described, starting from 2-formylpyrrole and employing Diels-Alder and Heck arylation reactions. 3-(N-Benzyl-2-pyrrolyl)acrylates and 4-(pyrrol-2-yl)butenones underwent a highly endo-Diels-Alder cycloaddition with maleimides to furnish octahydropyrrolo[3,4-e]indoles, which served as precursors in the regioselective synthesis of aza-polycyclic skeletons via an intramolecular Heck arylation reaction. Through the latter reaction, the 3-(N-benzyl-2-pyrrolyl)acrylates give rise to 3-(pyrrolo[2,1-a]isoindol-3-yl)acrylates. A further oxidative aromatization of the polycyclic intermediates provides the corresponding polycyclic pyrrolo-isoindoles and isoindolo-pyrrolo-indoles. A theoretical study on the stereoselective Diels-Alder reactions, carried out by calculating the endo/exo transition states, revealed the assistance of non-covalent interactions in governing the endo stereocontrol.

Evidence of Facilitated Transport in Crowded Nanopores

Fluid transport in nature often occurs through crowded nanopores, where a number of phenomena can affect it, because of fluid-fluid and fluid-solid interactions, as well as the presence of organic compounds filling the pores and their structural fluctuations. Employing molecular dynamics, we probe here the transport of fluid mixtures (CO₂-CH₄ and H₂S-CH₄) through silica nanopores filled with benzene. Both CO₂ and H₂S are strongly adsorbed within the organic-filled pore, partially displacing benzene. Unexpectedly, CO₂/H₂S adsorption facilitates CH₄ transport. Analysis of the trajectories suggests that both CO₂ and H₂S act as vehicle-like carriers and might swell benzene, generating preferential transport pathways within the crowded pore. The results are useful for identifying unexpected transport mechanisms and for developing engineering approaches that could lead to storage of CO₂ in caprocks.

exo/endo Selectivity Control in Diels-Alder Reactions of Geminal Bis(silyl) Dienes: Theoretical and Experimental Studies

The factors controlling the reactivity and exo/ endo selectivity of Diels-Alder reactions of geminal bis(silyl) dienes catalyzed by AlEt₂Cl are studied at the B3LYP-D3(BJ)(SMD,CH₂Cl₂)/6-31++G**//B3LYP-D3(BJ)(SMD,CH₂Cl₂)/6-31+G* theoretical level. The reaction proceeds via a two-stage one-step mechanism, and the AlEt₂Cl as a Lewis acid catalyst enhances the electrophilicity of the carbonyl compound by coordination, consequently accelerating a cycloaddition reaction with a low energy barrier. A geminal bis(silyl) group of the diene and an α-substituent in α,β-unsaturated carbonyl compounds adjust the interaction energy (Δ E_int) as well as the deformation energy (Δ E_strain) of the diene and dienophile, affecting the barrier height and the diastereochemical outcome accordingly. The steric repulsion between the geminal bis(silyl) group and Al(III) catalyst increases the Pauli repulsion energy (Δ E_Pauli) and strain energy of dienophile fragment (Δ E_{strain(dienophile-LA)}) in the endo pathway, ensuring the exo selectivity. The introduction of a halogen atom (Cl or Br) or methyl group at the α-position of α,β-unsaturated carbonyl compounds increases the deformation energy of the diene fragment. Meanwhile, the noncovalent interactions (that is, dispersion and electrostatic interaction) stabilize the endo transition state, leading to predominant endo products. The theoretical predictions of the exo/ endo selectivity for Diels-Alder reactions of the substituted α,β-unsaturated carbonyl compound with Cl or Br atoms by the DFT method are also well confirmed by experiment.

C-H/π and C-H-O Interactions in Concert: A Study of the Anisole-Methane Complex using Resonant Ionization and Velocity Mapped Ion Imaging

Noncovalent forces such as hydrogen bonding, halogen bonding, π-π stacking, and C-H/π and C-H/O interactions hold the key to such chemical processes as protein folding, molecular self-assembly, and drug-substrate interactions. Invaluable insight into the nature and strength of these forces continues to come from the study of isolated molecular clusters. In this work, we report on a study of the isolated anisole-methane complex, where both C-H/π and C-H/O interactions are possible, using a combination of theory and experiments that include mass-selected two-color resonant two-photon ionization spectroscopy, two-color appearance potential (2CAP) measurements, and velocity mapped ion imaging (VMI). Using 2CAP and VMI, we derive the binding energies of the complex in ground, excited, and cation radical states. The experimental values from the two methods are in excellent agreement, and they are compared with selected theoretical values calculated using density functional theory and ab initio methods. The optimized ground-state cluster geometry, which is consistent with the experimental observations, shows methane sitting above the ring, interacting with anisole via both C-H/π and C-H/O interactions, and this dual mode of interaction is reflected in a larger ground-state binding energy as compared with the prototypical benzene-methane system.

Molecular design principles towards exo-exclusive Diels–Alder reactions

The exo selective Diels–Alder reactions, reported as special cases, usually involve catalytic reaction conditions and specific cyclic structural motifs on the diene and/or the dienophile. Here we report a systematic computational investigation on the substituent effect for simple, linear dienes and dienophiles towards exo control in Diels–Alder reactions under thermal conditions. Through detailed characterization of reaction pathways for Diels–Alder cycloadditions between linear dienes and dienophiles with various substituents, we summarize a set of design principles aiming for an optimal and nearly-exclusive exo selectivity. These results shall lead to valuable guidelines and more versatile strategies in organic synthesis that are in accordance with the principles of green chemistry.

Exo-exclusive stereoselectivity for simple, terminal-substituted dienes and dienophiles may be achieved under thermal conditions through a delicate control of substituent identities.

Two Highly Stable Luminescent Lead Phosphonates Based on Mixed Ligands: Highly Selective and Sensitive Sensing for Thymine Molecule and VO3 - Anion

Two luminescent lead phosphonates with two-dimensional (2D) layer and three-dimensional (3D) framework structure, namely, Pb₃[(L₁)₂(Hssc)(H₂O)₂] (1) and [Pb₂(L₂)_0.5(bts)(H₂O)₂]·H₂O (2) (H₂L₁ = O(CH₂CH₂)₂NCH₂PO₃H₂, H₄L₂ = H₂PO₃CH₂NH(C₂H₄)₂NHCH₂PO₃H₂, H₃ssc = 5-sulfosalicylic acid, NaH₂bts = 5-sulfoisophthalic acid sodium) have been prepared via hydrothermal techniques. The two compounds not only show excellent thermal stability but also remain intact in aqueous solution within an extensive pH range. Moreover, the atomic absorption spectroscopy analysis experiment indicates that there does not exist the leaching of Pb²⁺ ions from the lead phosphonates, which show they are nontoxic in aqueous solution. In compound 1, the Pb(1)O₄, Pb(2)O₇, Pb(3)O₄, and CPO₃ polyhedra are interlinked into a one-dimensional chain, which is further connected to adjacent chain by sharing the Hssc^2- to form a 2D layer. Interestingly, compound 1 as a highly selective and sensitive luminescent material can be used to detect the thymine molecule with a very low detection limit of 8.26 × 10^-7 M. In compound 2, the Pb(1)O₆ and Pb(2)O₅ polyhedra are interlinked into a dimer via edge sharing, which is further connected to adjacent dimer to form a tetramer via corner sharing, and such a tetramer is then interlinked into a 2D layer through bts^3- ligands; the adjacent 2D layers are finally constructed to a 3D structure by sharing the L₂ ^4- ligand. Compound 2 can be applied as an excellent luminescent sensor for sensing of VO₃ ^- anion. Furthermore, the probable fluorescent quenching mechanisms of the two compounds have also been studied.

Weak interactions in furan dimers

Dimers of furan, 2,3-dihydrofuran, 2,5-dihydrofuran and tetrahydrofuran were investigated with the use of theoretical methods to determine the interactions that keep the molecules together. The QTAIM and NCI methods confirmed that for furan dimers the C-H⋯O hydrogen bond and stacking interactions can form the dimers with similar energy. For 2,3-dihydrofuran, 2,5-dihydrofuran and tetrahydrofuran, the decisive mechanism of dimer formation is the stacking interaction between the furan rings.

A quantification scheme for non-covalent interactions in the enantio-controlling transition states in asymmetric catalysis

The origin of enantioselectivity in asymmetric catalysis is attributed to the energy difference between lower and higher energy diastereomeric transition states, which are respectively responsible for the formation of major and minor enantiomers. Although the increase in the number of transition state models emphasizes the role of weak non-covalent interactions in asymmetric induction, the strength of such interactions is seldom quantified. Through this article, we propose a simple and effective method of quantifying the total non-covalent interaction in stereocontrolling transition states belonging to a group of three representative asymmetric catalytic reactions involving chiral phosphoric acids. Our method relies on rational partitioning of a given transition state into two (or three) sub-units, such that the complex network of intramolecular interactions can be ameliorated to a set of intermolecular interactions between two sub-units. The computed strength of interaction obtained using the counterpoise (CP) method on suitably partitioned transition states provides improved estimates of non-covalent interactions, which are also devoid of basis set superposition error (BSSE). It has been noted that catalysts decorated with larger aromatic arms provide cumulative non-covalent interactions (C-Hπ, N-Hπ and ππ) to the tune of 10 to 15 kcal mol-1. Fine-tuning of the magnitude and nature of these interactions can provide valuable avenues in the design of asymmetric catalysts.

Site specificity of halogen bonding involving aromatic acceptors

Based on Cambridge structural database survey and quantum chemical studies, aromatic halogen bond (XB) acceptors are found to have unique pattern of XB binding sites and rim specificity.

Intrinsically disordered sequences enable modulation of protein phase separation through distributed tyrosine motifs

Liquid–liquid phase separation (LLPS) is thought to contribute to the establishment of many biomolecular condensates, eukaryotic cell structures that concentrate diverse macromolecules but lack a bounding membrane. RNA granules control RNA metabolism and comprise a large class of condensates that are enriched in RNA-binding proteins and RNA molecules. Many RNA granule proteins are composed of both modular domains and intrinsically disordered regions (IDRs) having low amino acid sequence complexity. Phase separation of these molecules likely plays an important role in the generation and stability of RNA granules. To understand how folded domains and IDRs can cooperate to modulate LLPS, we generated a series of engineered proteins. These were based on fusions of an IDR derived from the RNA granule protein FUS (fused in sarcoma) to a multivalent poly-Src homology 3 (SH3) domain protein that phase-separates when mixed with a poly-proline–rich-motif (polyPRM) ligand. We found that the wild-type IDR promotes LLPS of the polySH3–polyPRM system, decreasing the phase separation threshold concentration by 8-fold. Systematic mutation of tyrosine residues in Gly/Ser-Tyr-Gly/Ser motifs of the IDR reduced this effect, depending on the number but not on the position of these substitutions. Mutating all tyrosines to non-aromatic residues or phosphorylating the IDR raised the phase separation threshold above that of the unmodified polySH3–polyPRM pair. These results show that low-complexity IDRs can modulate LLPS both positively and negatively, depending on the degree of aromaticity and phosphorylation status. Our findings provide plausible mechanisms by which these sequences could alter RNA granule properties on evolutionary and cellular timescales.

Crystallization-induced emission enhancement: A novel fluorescent Au-Ag bimetallic nanocluster with precise atomic structure

We report the first noble metal nanocluster with a formula of Au4Ag13(DPPM)3(SR)9 exhibiting crystallization-induced emission enhancement (CIEE), where DPPM denotes bis(diphenylphosphino)methane and HSR denotes 2,5-dimethylbenzenethiol. The precise atomic structure is determined by x-ray crystallography. The crystalline state of Au4Ag13 shows strong luminescence at 695 nm, in striking contrast to the weak emission of the amorphous state and hardly any emission in solution phase. The structural analysis and the density functional theory calculations imply that the compact C-H⋯π interactions significantly restrict the intramolecular rotations and vibrations and thus considerably enhance the radiative transitions in the crystalline state. Because the noncovalent interactions can be easily modulated via varying the chemical environments, the CIEE phenomenon might represent a general strategy to amplify the fluorescence from weakly (or even non-) emissive nanoclusters.

Effect of Lewis acid bulkiness on the stereoselectivity of Diels–Alder reactions between acyclic dienes and α,β-enals

The influence of Lewis acid bulkiness on the stereoselectivity of Diels–Alder reactions is analysed computationally in detail.

in/out Isomerism of cyclophanes: a theoretical account of 2,6,15-trithia-[34,10][7]metacyclophane and [34,10][7]metacyclophane as well as their halogen substituted analogues

What are the effects responsible for cyclophane isomer stability?

Torsional and Electronic Factors Control the C-H⋅⋅⋅O Interaction

The precise role of non‐conventional hydrogen bonds such as the C−H⋅⋅⋅O interaction in influencing the conformation of small molecules remains unresolved. Here we survey a series of β‐turn mimetics using X‐ray crystallography and NMR spectroscopy in conjunction with quantum calculation, and conclude that favourable torsional and electronic effects are important for the population of states with conformationally influential C−H⋅⋅⋅O interactions. Our results also highlight the challenge in attempting to deconvolute a myriad of interdependent noncovalent interactions in order to focus on the contribution of a single one. Within a small molecule that is designed to resemble the complexity of the environment within peptides and proteins, the interplay of different steric burdens, hydrogen‐acceptor/‐donor properties and rotational profiles illustrate why unambiguous conclusions based solely on NMR chemical shift data are extremely challenging to rationalize.

Analysis of the interactions of sulfur-containing amino acids in membrane proteins

The interactions of Met and Cys with other amino acid side chains have received little attention, in contrast to aromatic-aromatic, aromatic-aliphatic or/and aliphatic-aliphatic interactions. Precisely, these are the only amino acids that contain a sulfur atom, which is highly polarizable and, thus, likely to participate in strong Van der Waals interactions. Analysis of the interactions present in membrane protein crystal structures, together with the characterization of their strength in small-molecule model systems at the ab-initio level, predicts that Met-Met interactions are stronger than Met-Cys ≈ Met-Phe ≈ Cys-Phe interactions, stronger than Phe-Phe ≈ Phe-Leu interactions, stronger than the Met-Leu interaction, and stronger than Leu-Leu ≈ Cys-Leu interactions. These results show that sulfur-containing amino acids form stronger interactions than aromatic or aliphatic amino acids. Thus, these amino acids may provide additional driving forces for maintaining the 3D structure of membrane proteins and may provide functional specificity.

Analysis of transition state stabilization by non-covalent interactions in the Houk-List model of organocatalyzed intermolecular Aldol additions using functional-group symmetry-adapted perturbation theory

Rational design of catalysts would be aided by a better understanding of how non-covalent interactions stabilize transition states. Here, we apply the newly-developed Functional-Group Symmetry-Adapted Perturbation Theory (F-SAPT) to quantify non-covalent interactions in transition states of the proline-catalyzed intermolecular aldol reaction between benzaldehyde and cyclohexanone, according to the Houk-List mechanism [Bahmanyar et al., J. Am. Chem. Soc., 2003, 125, 2475]. A recent re-examination of this organocatalytic reaction by Rzepa and co-workers [Armstrong et al., Chem. Sci., 2014, 5, 2057] used electron density analysis to identify three key non-covalent interactions thought to influence stereoselectivity: (1) a favorable electrostatic interaction (originally identified by Houk and List) between the NCH(δ+) group of the enamine intermediate and the (δ-)O[double bond, length as m-dash]C of benzaldehyde; (2) a C-H/π interaction between the cyclohexene group of the enamine intermediate and the benzaldehyde phenyl ring; (3) a stabilizing contact between an ortho-hydrogen of the phenyl and an oxygen of the carboxylic acid group of the enamine. These three interactions have been directly computed using F-SAPT, which confirms the stabilizing interaction between an ortho-hydrogen and the carboxylic acid in the (S,S) and (R,S) transition state stereoisomers. F-SAPT analysis also finds stabilizing dispersion and electrostatic interactions due to a C-H/π interaction between the cyclohexene and phenyl groups in the (S,S) and (R,R) transition states. However, unfavorable exchange-repulsion cancels the attractive terms that favor these stereoisomers. Surprisingly, the interaction thought to be most important for stereoselectivity, the NCH(δ+)(δ-)O[double bond, length as m-dash]C interaction, is actually found to be repulsive due to the negative charge on the nitrogen. Hence, our results indicate that geometric analysis and/or density-based analysis does not necessarily produce a reliable picture of non-covalent stabilization. As confirmed by high-level coupled-cluster computations, intermolecular interaction energies are strongest for the (R,R) transition states, which are not the experimentally favored products. This suggests that at least for this reaction, stereoselectivity is also strongly dependent on the energy required to distort the reacting molecules into the transition state geometry.

An Isochemogenic Set of Inhibitors To Define the Therapeutic Potential of Histone Deacetylases in β-Cell Protection

Modulation of histone deacetylase (HDAC) activity has been implicated as a potential therapeutic strategy for multiple diseases. However, it has been difficult to dissect the role of individual HDACs due to a lack of selective small-molecule inhibitors. Here, we report the synthesis of a series of highly potent and isoform-selective class I HDAC inhibitors, rationally designed by exploiting minimal structural changes to the clinically experienced HDAC inhibitor CI-994. We used this toolkit of isochemogenic or chemically matched inhibitors to probe the role of class I HDACs in β-cell pathobiology and demonstrate for the first time that selective inhibition of an individual HDAC isoform retains beneficial biological activity and mitigates mechanism-based toxicities. The highly selective HDAC3 inhibitor BRD3308 suppressed pancreatic β-cell apoptosis induced by inflammatory cytokines, as expected, or now glucolipotoxic stress, and increased functional insulin release. In addition, BRD3308 had no effect on human megakaryocyte differentiation, while inhibitors of HDAC1 and 2 were toxic. Our findings demonstrate that the selective inhibition of HDAC3 represents a potential path forward as a therapy to protect pancreatic β-cells from inflammatory cytokines and nutrient overload in diabetes.