High-Accuracy Quantum Mechanical Studies of π−π Interactions in Benzene Dimers

MP2-Based Composite Extrapolation Schemes Can Predict Core-Ionization Energies for First-Row Elements with Coupled-Cluster Level Accuracy

X-ray photoelectron spectroscopy (XPS) measures core-electron binding energies (CEBEs) to reveal element-specific insights into the chemical environment and bonding. Accurate theoretical CEBE prediction aids XPS interpretation but requires proper modeling of orbital relaxation and electron correlation upon core-ionization. This work systematically investigates basis set selection for extrapolation to the complete basis set limit of CEBEs from ΔMP2 and ΔCC energies across 94 K-edges in diverse organic molecules. We demonstrate that an alternative composite scheme using ΔMP2 in a large basis corrected by ΔCC-ΔMP2 difference in a small basis can quantitatively recover optimally extrapolated ΔCC CEBEs within 0.02 eV. Unlike ΔCC, MP2 calculations do not suffer from convergence issues and are computationally cheaper, and thus, the composite ΔMP2/ΔCC scheme balances accuracy and cost, overcoming limitations of solely using either method. We conclude by providing a comprehensive analysis of the choice of small and large basis sets for the composite schemes and provide practical recommendations for highly accurate (within 0.10-0.15 eV MAE) ab initio prediction of XPS data.

Exploring the persistence of the fluorinated thiolate 2,3,5,6-S(C6F4H-4) motif to establish πF-πF stacking in metal complexes: a crystal engineering perspective

π-π stacking interactions are versatile because they are involved in many processes, such as protein folding, DNA stacking, and drug recognition. However, from the point of view of crystal engineering, there is an incipient knowledge of its exploitation. A comparison of these interactions with hydrogen bonds shows a huge difference in their employment as a reliable non-covalent interaction. And different reasons can be listed to explain why hydrogen bonding can be considered a more robust interaction than π-π stacking. For instance, hydrogen bonds encompass a wide energy range (25-40 kJ mol-1). From this, these interactions can be classified as strong, moderate, and weak. Hence, the first two can be considered highly to moderately directional to be exploited in crystal engineering. This aspect is relevant for them to be used in a relatively reliably way in this area of supramolecular chemistry. On the other hand, in the case of π-π stacking, the energy range is 0-10 kJ mol-1, thus implying that hydrogen bonds or any other energetically more robust contact would predominate in the competition for establishing packing interactions in a given arrangement. In this sense, if stacking is pretended to be exploited from the point of view of crystal engineering, one of the points that must be ensured is that this interaction will be the one energetically predominant. However, although there are other factors to consider, it seems that energetics is the dominant one. In this line, our research group has obtained and studied many single-crystalline structures of coordination and organometallic compounds containing fluorinated thiolates. This being particularly true in the case of the thiolate 2,3,5,6-S(C6F4H-4) bound to different metals, where it has been observed that they preferentially tend to establish πF-πF stacking interactions, results that have been reported in several papers. Thus, from this perspective, we have explored, using ConQuest (CCDC) a number of structures to observe how feasible is to find stacking in coordination and organometallic compounds containing the thiolate 2,3,5,6-S(C6F4H-4).

Computational investigation on the binding modes of PET polymer to PETase

Poly(ethylene terephthalate) (PET) has been widely utilized in daily life, but its non-degradability has induced severe environmental and health problems. Recently, PETase, which has been isolated from bacterium Ideonella sakaiensisis, was reported to have the highest PET degradation activity and specificity under room temperature, but no crystal structure for PET in complex with PETase has been reported. To provide deep insight into the binding mode of PET polymer on PETase and the binding interactions, we employed molecular docking and molecular dynamics simulations to study the substrate binding at the atomic level. Different PET oligomers have been studied with chain lengths varying from 2 to 8. In addition, the binding energies and hot-spot residues were analyzed to gain better insights into the binding mechanism by MM/GBSA approach. The PET oligomers adopt stable and reactive conformations in a shallow cleft on a flat surface of PETase. The binding cleft can only accommodate four moieties, and others beyond the region will be stabilized by the π-stacking interactions with Trp156 at the terephthalic acid terminal. Our studies provide a clear picture of how the binding mode of PET polymer and its interactions with PETase change with the chain length. Those studies would provide useful information for the rational design of catalytically more efficient PETase variants toward plastic degradation.Communicated by Ramaswamy H. Sarma.

Quantum Mechanics-Based Fast and Reliable Prediction of Binding Pose Structures

Predicting the binding poses of docking with an accurate estimation of binding energies is highly important but very challenging in computational drug design. A quantum mechanics (QM) calculation-based docking approach considering multiple conformations and orientations of the ligand is introduced here to tackle the problem. This QM docking consists of three steps: generating an ensemble of binding poses with a conventional docking simulation, computing the binding energies with self-consistent charge density functional theory tightly binding with dispersion correction (DFTB-D) to selecting the 10 top binding modes, and optimizing the selected binding mode structures using the ONIOM(DFTB:PM7) technique to determine the binding poses. The ONIOM(DFTB-D:PM6) docking approach is tested on 121 ligand-receptor biocomplexes with the crystal structures obtained from the Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB). The result shows that the new method is highly satisfactory for the accurate prediction of the binding poses. The new docking method should be beneficial to structure-based drug design.

Managing Pb-Related Imperfections via Rationally Designed Aniline Derivative with Bilateral Cyano and Acetyl Groups as Lewis Base for High-Efficiency Perovskite Solar Cells Exceeding 24

Pb-related imperfections (surface or halide vacancy induced uncoordinated Pb2+, Pb-I antisite, and Pb2+ vacancy defects) of the ionic crystal perovskite film seriously restrict the photovoltaic performance of perovskite solar cells (PSCs). Here, an aniline derivative N-(4-cyanophenyl)acetamide (CAL) is rationally designed, incorporating bilateral functional sites of cyano and acetyl groups, acting as Lewis base molecule for managing the Pb-related imperfections in perovskite surface through post-treatment. Theoretical calculation and experimental verification together proved the reduced defect density, improved crystallinity, and inhibited ion migration in the CAL-modified perovskite. Precisely, cyano as a side group and acetyl as another side group can both coordinate with Pb2+ for its low electrostatic potential energy. Further, the aniline core and the π-π conjugate structure in the benzene ring of the ligand tend to form a dimer to improve the mobility for carrier transportation and collection. The strategy demonstrates a champion PCE of 24.35% for the air-processed PSCs with over 1200 hours of maximum power point tracking (MPPT) stability. This study presents a comprehensive approach to overcoming the current Pb-related imperfections induced limitations in PSCs, paving the way for their integration into mainstream solar technologies.

Closed-loop recyclability of a biomass-derived epoxy-amine thermoset by methanolysis

Epoxy resin thermosets (ERTs) are an important class of polymeric materials. However, owing to their highly cross-linked nature, they suffer from poor recyclability, which contributes to an unacceptable level of environmental pollution. There is a clear need for the design of inherently recyclable ERTs that are based on renewable resources. We present the synthesis and closed-loop recycling of a fully lignocellulose-derivable epoxy resin (DGF/MBCA), prepared from dimethyl ester of 2,5-furandicarboxylic acid (DMFD), 4,4'-methylenebis(cyclohexylamine) (MBCA), and glycidol, which displays excellent thermomechanical properties (a glass transition temperature of 170°C, and a storage modulus at 25°C of 1.2 gigapascals). Notably, the material undergoes methanolysis in the absence of any catalyst, regenerating 90% of the original DMFD. The diamine MBCA and glycidol can subsequently be reformed by acetolysis. Application and recycling of DGF/MBCA in glass and plant fiber composites are demonstrated.

Cation desolvation-induced capacitance enhancement in reduced graphene oxide (rGO)

Understanding the local electrochemical processes is of key importance for efficient energy storage applications, including electrochemical double layer capacitors. In this work, we studied the charge storage mechanism of a model material - reduced graphene oxide (rGO) - in aqueous electrolyte using the combination of cavity micro-electrode, operando electrochemical quartz crystal microbalance (EQCM) and operando electrochemical dilatometry (ECD) tools. We evidence two regions with different charge storage mechanisms, depending on the cation-carbon interaction. Notably, under high cathodic polarization (region II), we report an important capacitance increase in Zn2+ containing electrolyte with minimum volume expansion, which is associated with Zn2+ desolvation resulting from strong electrostatic Zn2+-rGO interactions. These results highlight the significant role of ion-electrode interaction strength and cation desolvation in modulating the charging mechanisms, offering potential pathways for optimized capacitive energy storage. As a broader perspective, understanding confined electrochemical systems and the coupling between chemical, electrochemical and transport processes in confinement may open tremendous opportunities for energy, catalysis or water treatment applications in the future.

The effect of weak π-π interactions on single-molecule electron transport properties of the tetraphenylethene molecule and its derivatives: a first-principles study

Intramolecular π-π interactions are a significant research focus in fields such as chemistry, biology, and materials science. Different configurations of benzene-benzene moieties within a molecule can affect the magnitude of their π-π interactions, consequently influencing the electronic transport capabilities of the molecule. In this study, we designed three π-conjugated molecules, TPEM, TPEEM, and TEEPM, based on tetraphenylethene (TPE). These three molecules exhibit three distinct π-conjugated structures: linear cis-π-conjugation, linear trans-π-conjugation, and cross-π-conjugation. Thereinto, TPEM and TPEEM molecules share the same TPE core, with identical π-π interaction distances, while the TEEPM molecule has acetylene groups between the TPE units, thereby increasing the π-π interaction distances between the benzene moieties. Using density functional theory calculations combined with non-equilibrium Green's function (DFT+NEGF), our results reveal that the conductance order of different π-conjugated structures in TPEM and TPEEM molecules is as follows: cis > cross ≈ trans. Through analysis of transmission spectra, transmission pathways, and the innermost π orbitals, we find that in TPEM and TPEEM molecules, the cis- and cross-π-conjugated structures exhibit π-π interactions between benzene moieties and provide special through-space electron transport pathways, enhancing their electronic transport capabilities in coordination with the bonded molecular framework, whereas their trans-conjugated structures only allow electron transport along the molecular backbone. In contrast, in TEEPM molecule, due to the absence of π-π interactions, the conductance of different π-conjugated structures is primarily determined by the molecular backbone and follows the order: trans > cis > cross. These findings provide a theoretical basis for designing single-molecule electronic devices with multiple electron channels based on intramolecular π-π interactions.

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

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

Excitonic Coupling in Fluorene-Based Bichromophoric Systems: Vibrational Quenching and the Transition from Weak to Intermediate Coupling

Excimeric systems (i.e., excited dimers) have well served as model compounds for the study of the delocalization of electronic energy over weakly interacting chromophores. However, there remain relatively few isolated systems in which such interactions can be studied experimentally at a level to afford detailed comparisons with theory. In this Article, we examine a series of covalently and noncovalently linked dimers of fluorene, as a model aromatic chromophore, where the formation of excimers requires a π-stacked, cofacial orientation at van der Waals contact. Building upon a series of seminal prior studies that examined vibronic quenching of the excitation interaction in van der Waals dimers, the key question that we sought to address here is whether a single quenching factor could reproduce experimental excitonic splittings across a series of covalently and noncovalently linked bichromophoric systems built from the same chromophore. In comparing experimentally measured excitonic splittings with calculated static splittings using time-dependent density functional methods, we find that all systems save one fall on a line with a slope of 0.080(8), reflecting a vibrational quenching of roughly 1 order of magnitude. The outlier, which shows a significantly reduced quenching factor, represents a cyclophane-linked system where the fluorene moieties are constrained in a cofacial arrangement. We argue that this system evidences the transition from the weak to intermediate coupling regime.

Adsorption of glyphosate and metabolite aminomethylphosphonic acid (AMPA) from water by polymer-based spherical activated carbon (PBSAC)

Glyphosate (GLY) is the most commonly used herbicide worldwide, and aminomethylphosphonic acid (AMPA) is its main metabolite. Their occurrence in ground and surface waters causes diseases in humans, while complex physico-chemical properties hinder detection and effective removal. Polymer-based spherical activated carbon (PBSAC) can adsorb many micropollutants efficiently and, hence, overcome the shortfalls of conventional treatment methods. The static adsorption of a mixture of GLY and AMPA by PBSAC was investigated with varying PBSAC properties and relevant solution chemistry. The results show that PBSAC can remove 95% GLY and 57% AMPA from an initial concentration of 1 µg/L at pH 8.2. PBSAC properties (size, activation level, and surface charge) have a strong influence on herbicide removal, where surface area plays a key role. Low to neutral pH favors non-charge interactions and results in good adsorption, while higher temperatures equally enhance GLY/AMPA adsorption by PBSAC. The work demonstrated the effective removal of GLY to meet the European guideline concentration (0.1 µg/L), while AMPA could not be removed to the required level.

Systematic elucidation of the second coordination sphere effect on the structure and properties of a blue copper protein, pseudoazurin

The rational structural and computational studies of a blue copper protein, pseudoazurin (PAz), and its Met16X (X = Phe, Leu, Val, Ile) variants gave clear functional meanings of the noncovalent interaction (NCI) through the second coordination sphere. The high-resolution X-ray crystal structures of Met16X PAz demonstrated that the active site geometry is significantly affected by the substitution of Met16, which is located within the NCI distance from the His81 imidazole ring at the copper active site. The computational chemistry calculations based on the crystal structure analyses confirmed that the NCI of S-π/CH-π (wild-type), π-π (Met16Phe), double CH-π (Met16Leu), and single CH-π (Met16Val and Met16Ile). The estimated interaction energies for the NCI demonstrated that the fine-tuning of the protein stability and Cu site properties form the second coordination sphere of PAz.

How great is the stabilization of crowded polyphenylbiphenyls by London dispersion?

Decaphenylbiphenyl (1) and 2,2',4,4',6,6'-hexaphenylbiphenyl (2) are bulky molecules expected to be greatly destabilized by steric crowding. Herein, through a combined experimental and computational approach, we evaluate the molecular energetics of crowded biphenyls. This is complemented by the study of phase equilibria for 1 and 2. Compound 1 shows a rich phase behavior, displaying an unusual interconversion between two polymorphs. Surprisingly, the polymorph with distorted molecules of C₁ symmetry is found to have the highest melting point and to be the one that is preferentially formed. The thermodynamic results also indicate that the polymorph displaying the more regular D₂ molecular geometry has larger heat capacity and is probably the more stable at lower temperatures. The melting and sublimation data clearly reveal the weakening of cohesive forces in crowded biphenyls due to the lower molecular surface area. The experimental quantification of the intramolecular interactions in 1 and 2 indicated, using homodesmotic reactions, a molecular stabilization of about 30 kJ mol^-1. We attribute the origin of this stabilization in both compounds to the existence of two parallel-displaced π⋯π interactions between the ortho-phenyl substituents on each side of the central biphenyl. Computational calculations with dispersion-corrected DFT methods underestimate the stabilization in 1, unless the steric crowding is well balanced in a homodesmotic scheme. This work demonstrates that London dispersion forces are important in crowded aromatic systems, making these molecules considerably more stable than previously thought.

Adsorption behavior of aniline pollutant on polystyrene microplastics

Microplastic contamination is ubiquitous in aquatic environments. As global plastic production increases, the abundance of microplastic contaminants released into the environment has also continued to soar. The hydrophobic surfaces of plastic particles can adsorb a variety of chemical pollutants, and could therefore facilitate toxin accumulation through the food chain. In this study, the adsorption behavior of aniline, a priority environmental pollutant from industrial production, on the surface of polystyrene microplastics (mPS) was investigated. The results showed that the maximum adsorption capacity of mPS was 0.060 mg/g. Adsorption equilibrium was reached after 24 h, and the pseudo-second-order model was employed to explain the adsorption kinetics of aniline on the mPS particles. The Freundlich models could describe the adsorption isotherms. The potential adsorption mechanisms may include π-π interactions and hydrophobic interactions. pH, ionic strength, and ambient temperature of the solution played important roles in the adsorption process.

High-Order Quantum-Mechanical Analysis of Hydrogen Bonding in Hachimoji and Natural DNA Base Pairs

High-order quantum chemistry is applied to hydrogen-bonded natural DNA nucleobase pairs [adenine:thymine (A:T) and guanine:cytosine (G:C)] and non-natural Hachimoji nucleobase pairs [isoguanine:1-methylcytosine (B:S) and 2-aminoimidazo[1,2a][1,3,5]triazin-4(1H)-one:6-amino-5-nitropyridin-2-one (P:Z)] to see how the intermolecular interaction energies and their energetic components (electrostatics, exchange-repulsion, induction/polarization, and London dispersion interactions) vary among the base pairs. We examined the Hoogsteen (HG) geometries in addition to the traditional Watson-Crick (WC) geometries. Coupled-cluster theory through perturbative triples [CCSD(T)] extrapolated to the complete basis set (CBS) limit and high-order symmetry-adapted perturbation theory (SAPT) at the SAPT2+(3)(CCD)δMP2/aug-cc-pVTZ level are used to estimate highly accurate noncovalent interaction energies. Electrostatic interactions are the most attractive component of the interaction energies, but the sum of induction/polarization and London dispersion is nearly as large, for all base pairs and geometries considered. Interestingly, the non-natural Hachimoji base pairs interact more strongly than the corresponding natural base pairs, by -21.8 (B:S) and -0.3 (P:Z) kcal mol^-1 in the WC geometries, according to CCSD(T)/CBS. This is consistent with the H-bond distances being generally shorter in the non-natural base pairs. The natural base pairs are energetically more stabilized in their Hoogsteen geometries than in their WC geometries. The Hoogsteen geometry makes the A:T base pair slightly more stable, by -0.8 kcal mol^-1, and it greatly stabilizes the G:C⁺ base pair, by -15.3 kcal mol^-1. The G:C⁺ stabilization is mainly due to the fact that C has typically added a proton when found in Hoogsteen geometries. By contrast, Hoogsteen geometries are substantially less favorable than WC geometries for non-natural Hachimoji base pairs, by 17.3 (B:S) and 13.8 (P:Z) kcal mol^-1.

Effect of stacking interactions on charge transfer states in photoswitches interacting with ion channels

The activity of ion channels can be reversibly photo-controlled via the binding of molecular photoswitches, often based on an azobenzene scaffold. Those azobenzene derivatives interact with aromatic residues of the protein via stacking interactions. In the present work, the effect of face-to-face and t-shaped stacking interactions on the excited state electronic structure of azobenzene and p-diaminoazobenzene integrated into the NaV1.4 channel is computationally investigated. The formation of a charge transfer state, caused by electron transfer from the protein to the photoswitches, is observed. This state is strongly red shifted when the interaction takes place in a face-to-face orientation and electron donating groups are present on the aromatic ring of the amino acids. The low-energy charge transfer state can interfere with the photoisomerization process after excitation to the bright state by leading to the formation of radical species.

Using a diphenyl-bi-(1,2,4-triazole) tricarbonylrhenium(I) complex with intramolecular π-π stacking interaction for efficient solid-state luminescence enhancement

Since intramolecular π-π stacking interactions can modify the geometry, crystal packing mode, or even the electronic properties of transition metal complexes, they are also likely to influence the solid-state luminescence properties. Following this concept, a new tricarbonylrhenium(I) complex (Re-BPTA) was designed, based on a simple symmetrical 5,5'-dimethyl-4,4'-diphenyl-3,3'-bi-(1,2,4-triazole) organic ligand. The complex was prepared in good yield using a three-step procedure. The crystallographic study revealed that both phenyl rings are located on the same side of the molecule, and twisted by 71° and 62°, respectively, with respect to the bi-(1,2,4-triazole) unit. They overlap significantly, although they are slipped parallel to each other to minimize the intramolecular interaction energy. The π-π stacking interaction was also revealed by ¹H NMR spectroscopy, in good agreement with the results of theoretical calculations. In organic solutions, a peculiar electrochemical signature was observed compared to closely-related pyridyl-triazole (pyta)-based complexes. With regard to the optical properties, the stiffness of the Re-BPTA complex led to the stabilization of the ³MLCT state, and thus to an enhancement of the red phosphorescence emission compared to the more flexible pyta complexes. However, an increased sensitivity to quenching by oxygen appeared. In the microcrystalline phase, the Re-BPTA complex showed strong photoluminescence (PL) emission in the green-yellow wavelength range (λ_PL = 548 nm, Φ_PL = 0.52, 〈τ_PL〉 = 713 ns), and thus a dramatic solid-state luminescence enhancement (SLE) effect. These attractive emission properties can be attributed to the fact that the molecule undergoes little distortion between the ground state and the triplet excited state, as well as to a favorable intermolecular arrangement that minimizes detrimental interactions in the crystal lattice. The aggregation-induced phosphorescence emission (AIPE) effect was clear, with a 7-fold increase in emission intensity at 546 nm, although the aggregates formed in aqueous medium were much less emissive than the native microcrystalline powder. In this work, the rigidity of the Re-BPTA complex is reinforced by the intramolecular π-π stacking interaction of the phenyl rings. This original concept provides a rhenium tricarbonyl compound with very good SLE properties, and could be used more widely to successfully develop this area of research.

Structures and Supramolecular Properties of Inclusion Complexes of Anthracene-Triptycene Nanocages with Fullerene Guests and Their Dynamic Motion as Molecular Gyroscopes

Three derivatives of macrocyclic cage compounds consisting of diarylanthracene and triptycene units were synthesized. These nanocages formed host-guest complexes with C₆₀ and other fullerene guests as confirmed by ¹ H NMR and fluorescence spectroscopy. The association constant of the mesityl and 2,4,6-tributoxyphenyl derivatives with C₆₀ was determined to be 2.2 × 10⁴  L mol^-1 , which was larger than that of the pentafluorophenyl derivative. Direct experimental evidence of the complexation was obtained by X-ray diffraction analysis: the guest C₆₀ molecule was included in the cavity via multipoint CH⋅⋅⋅π interactions. Dynamic disorders of the included C₆₀ molecule in variable-temperature X-ray analysis indicated uniaxial motion, such as gyroscopic motion. The unique dynamic behavior of the spherical C₆₀ rotor anchored by the cage stator via CH⋅⋅⋅π interactions in the crystal, as well as substituent effects on the association properties, are discussed with the aid of DFT calculations.

Rationalizing the Binding Modes of PET Radiotracers Targeting the Norepinephrine Transporter

PURPOSE

A new PET radiotracer ¹⁸F-AF78 showing great potential for clinical application has been reported recently. It belongs to a new generation of phenethylguanidine-based norepinephrine transporter (NET)-targeting radiotracers. Although many efforts have been made to develop NET inhibitors as antidepressants, systemic investigations of the structure-activity relationships (SARs) of NET-targeting radiotracers have rarely been performed.

METHODS

Without changing the phenethylguanidine pharmacophore and 3-fluoropropyl moiety that is crucial for easy labeling, six new analogs of ¹⁸F-AF78 with different meta-substituents on the benzene-ring were synthesized and evaluated in a competitive cellular uptake assay and in in vivo animal experiments in rats. Computational modeling of these tracers was established to quantitatively rationalize the interaction between the radiotracers and NET.

RESULTS

Using non-radiolabeled reference compounds, a competitive cellular uptake assay showed a decrease in NET-transporting affinity from meta-fluorine to iodine (0.42 and 6.51 µM, respectively), with meta-OH being the least active (22.67 µM). Furthermore, in vivo animal studies with radioisotopes showed that heart-to-blood ratios agreed with the cellular experiments, with AF78(F) exhibiting the highest cardiac uptake. This result correlates positively with the electronegativity rather than the atomic radius of the meta-substituent. Computational modeling studies revealed a crucial influence of halogen substituents on the radiotracer-NET interaction, whereby a T-shaped π-π stacking interaction between the benzene-ring of the tracer and the amino acid residues surrounding the NET binding site made major contributions to the different affinities, in accordance with the pharmacological data.

CONCLUSION

The SARs were characterized by in vitro and in vivo evaluation, and computational modeling quantitatively rationalized the interaction between radiotracers and the NET binding site. These findings pave the way for further evaluation in different species and underline the potential of AF78(F) for clinical application, e.g., cardiac innervation imaging or molecular imaging of neuroendocrine tumors.

The Auger spectrum of benzene

We present an ab initio computational study of the Auger electron spectrum of benzene. Auger electron spectroscopy exploits the Auger-Meitner effect, and although it is established as an analytic technique, the theoretical modeling of molecular Auger spectra from first principles remains challenging. Here, we use coupled-cluster theory and equation-of-motion coupled-cluster theory combined with two approaches to describe the decaying nature of core-ionized states: (i) Feshbach-Fano resonance theory and (ii) the method of complex basis functions. The spectra computed with these two approaches are in excellent agreement with each other and also agree well with experimental Auger spectra of benzene. The Auger spectrum of benzene features two well-resolved peaks at Auger electron energies above 260 eV, which correspond to final states with two electrons removed from the 1e_1g and 3e_2g highest occupied molecular orbitals. At lower Auger electron energies, the spectrum is less well resolved, and the peaks comprise multiple final states of the benzene dication. In line with theoretical considerations, singlet decay channels contribute more to the total Auger intensity than the corresponding triplet decay channels.

Mussel-Based Biomimetic Strategies in Musculoskeletal Disorder Treatment: From Synthesis Principles to Diverse Applications

Musculoskeletal disorders are the second leading cause of disability worldwide, posing a huge global burden to the public sanitation system. Currently, tissue engineering-based approaches act as effective strategies, which are, however, challenging in limited application scenarios. Mussel-based biomimetic materials, exhibit numerous unique properties such as intense adhesion, biocompatibility, moisture resistance, and injectability, to name only a few, and have attracted extensive research interest. In particular, featuring state-of-the-art properties, mussel-inspired biomaterials have been widely explored in innumerable musculoskeletal disorder treatments including osteochondral defects, osteosarcoma, osteoarthritis, ligament rupture, and osteoporosis. Nevertheless, a comprehensive and timely discussion of their applications in musculoskeletal disorders is insufficient. In this review, we emphasize on (1) the main categories and characteristics of mussel foot proteins and their fundamental mechanisms for the spectacular adhesion in mussels; (2) the diverse synthetic methods and modification of various polymers; and (3) the emerging applications of mussel-biomimetic materials, the future perspectives, and challenges, especially in the area of musculoskeletal disorder. We envision that this review will provide a unique and insightful perspective to improve the development of a new generation of mussel biomimetic strategies.

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.

Thermodynamics of π–π Interactions of Benzene and Phenol in Water

The π–π interaction is a major driving force that stabilizes protein assemblies during protein folding. Recent studies have additionally demonstrated its involvement in the liquid–liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs). As the participating residues in IDPs are exposed to water, π–π interactions for LLPS must be modeled in water, as opposed to the interactions that are often established at the hydrophobic domains of folded proteins. Thus, we investigated the association of free energies of benzene and phenol dimers in water by integrating van der Waals (vdW)-corrected density functional theory (DFT) and DFT in classical explicit solvents (DFT-CES). By comparing the vdW-corrected DFT and DFT-CES results with high-level wavefunction calculations and experimental solvation free energies, respectively, we established the quantitative credibility of these approaches, enabling a reliable prediction of the benzene and phenol dimer association free energies in water. We discovered that solvation influences dimer association free energies, but not significantly when no direct hydrogen-bond-type interaction exists between two monomeric units, which can be explained by the enthalpy–entropy compensation. Our comprehensive computational study of the solvation effect on π–π interactions in water could help us understand the molecular-level driving mechanism underlying the IDP phase behaviors.

σ-σ Stacked supramolecular junctions

Intermolecular charge transport plays an essential role in organic electronic materials and biological systems. To date, experimental investigations of intermolecular charge transport in molecular materials and electronic devices have been restricted to conjugated systems in which π-π stacking interactions are involved. Herein we demonstrate that the σ-σ stacking interactions between neighbouring non-conjugated molecules offer an efficient pathway for charge transport through supramolecular junctions. The conductance of σ-σ stacked molecular junctions formed between two non-conjugated cyclohexanethiol or single-anchored adamantane molecules is comparable to that of π-π stacked molecular junctions formed between π-conjugated benzene rings. The current-voltage characteristics and flicker noise analysis demonstrate the existence of stacked molecular junctions formed between the electrode pairs and exhibit the characteristics of through-space charge transport. Density functional theory calculations combined with the non-equilibrium Green's function method reveal that efficient charge transport occurs between two molecules configured with σ-σ stacking interactions.

Comparison of approximate intermolecular potentials for ab initio fragment calculations on medium sized N-heterocycles

The ground state intermolecular potential of bimolecular complexes of N-heterocycles is analyzed for the impact of individual terms in the interaction energy as provided by various, conceptually different theories. Novel combinations with several formulations of the electrostatic, Pauli repulsion, and dispersion contributions are tested at both short- and long-distance sides of the potential energy surface, for various alignments of the pyrrole dimer as well as the cytosine-uracil complex. The integration of a DFT/CCSD density embedding scheme, with dispersion terms from the effective fragment potential (EFP) method is found to provide good agreement with a reference CCSD(T) potential overall; simultaneously, a quantum mechanics/molecular mechanics approach using CHELPG atomic point charges for the electrostatic interaction, augmented by EFP dispersion and Pauli repulsion, comes also close to the reference result. Both schemes have the advantage of not relying on predefined force fields; rather, the interaction parameters can be determined for the system under study, thus being excellent candidates for ab initio modeling.

Evaluating Geometric Definitions of Stacking for RNA Dinucleoside Monophosphates Using Molecular Mechanics Calculations

RNA modulation via small molecules is a novel approach in pharmacotherapies, where the determination of the structural properties of RNA motifs is considered a promising way to develop drugs capable of targeting RNA structures to control diseases. However, due to the complexity and dynamic nature of RNA molecules, the determination of RNA structures using experimental approaches is not always feasible, and computational models employing force fields can provide important insight. The quality of the force field will determine how well the predictions are compared to experimental observables. Stacking in nucleic acids is one such structural property, originating mainly from London dispersion forces, which are quantum mechanical and are included in molecular mechanics force fields through nonbonded interactions. Geometric descriptions are utilized to decide if two residues are stacked and hence to calculate the stacking free energies for RNA dinucleoside monophosphates (DNMPs) through statistical mechanics for comparison with experimental thermodynamics data. Here, we benchmark four different stacking definitions using molecular dynamics (MD) trajectories for 16 RNA DNMPs produced by two different force fields (RNA-IL and ff99OL3) and show that our stacking definition better correlates with the experimental thermodynamics data. While predictions within an accuracy of 0.2 kcal/mol at 300 K were observed in RNA CC, CU, UC, AG, GA, and GG, stacked states of purine-pyrimidine and pyrimidine-purine DNMPs, respectively, were typically underpredicted and overpredicted. Additionally, population distributions of RNA UU DNMPs were poorly predicted by both force fields, implying a requirement for further force field revisions. We further discuss the differences predicted by each RNA force field. Finally, we show that discrete path sampling (DPS) calculations can provide valuable information and complement the MD simulations. We propose the use of experimental thermodynamics data for RNA DNMPs as benchmarks for testing RNA force fields.

The influence of a solvent environment on direct non-covalent interactions between two molecules: A symmetry-adapted perturbation theory study of polarization tuning of π-π interactions by water

High-level quantum chemical computations have provided significant insight into the fundamental physical nature of non-covalent interactions. These studies have focused primarily on gas-phase computations of small van der Waals dimers; however, these interactions frequently take place in complex chemical environments, such as proteins, solutions, or solids. To better understand how the chemical environment affects non-covalent interactions, we have undertaken a quantum chemical study of π-π interactions in an aqueous solution, as exemplified by T-shaped benzene dimers surrounded by 28 or 50 explicit water molecules. We report interaction energies (IEs) using second-order Møller-Plesset perturbation theory, and we apply the intramolecular and functional-group partitioning extensions of symmetry-adapted perturbation theory (ISAPT and F-SAPT, respectively) to analyze how the solvent molecules tune the π-π interactions of the solute. For complexes containing neutral monomers, even 50 explicit waters (constituting a first and partial second solvation shell) change total SAPT IEs between the two solute molecules by only tenths of a kcal mol^-1, while significant changes of up to 3 kcal mol^-1 of the electrostatic component are seen for the cationic pyridinium-benzene dimer. This difference between charged and neutral solutes is attributed to large non-additive three-body interactions within solvated ion-containing complexes. Overall, except for charged solutes, our quantum computations indicate that nearby solvent molecules cause very little "tuning" of the direct solute-solute interactions. This indicates that differences in binding energies between the gas phase and solution phase are primarily indirect effects of the competition between solute-solute and solute-solvent interactions.

Novel Thermostable Heparinase Based on the Genome of Bacteroides Isolated from Human Gut Microbiota

Among the nutrients available to the human gut microbiota, the complex carbohydrates and glycosaminoglycans are important sources of carbon for some of the species of human gut microbiota. Glycosaminoglycan (heparin) from the host is a highly preferred carbohydrate for Bacteroides. To explore how gut microbiota can effectively use heparin as a carbon source for growth, we conducted a screening of the Carbohydrate-Active enzymes (CAZymes) database for lytic enzymes of the PL13 family and Research Center of Food Biotechnology at School of Food Science and Technology of Jiangnan University database of Bacteroides to identify novel glycosaminoglycan-degrading bacterial strains. Four Bacteroides species (Bacteroides eggerthii, Bacteroides clarus, Bacteroides nordii, and Bacteroides finegoldii) that degraded heparin were selected for further studies. Analysis of the polysaccharide utilization sites of the four strains revealed that all of them harbored enzyme encoding genes of the PL13 family. Functional analysis revealed the activity of CAZymes in a medium containing heparin as the sole carbon source, suggesting their potential to degrade heparin and support growth. The four enzymes were heterologous expressed, and their enzymatic properties, kinetics, and thermal stability were determined. The lytic enzyme of B. nordii had high enzymatic activity and thermal stability. The features that cause this high thermal stability were elucidated based on an examination of the three-dimensional structure of the protein. Our findings provide an important theoretical basis for the application of glycosaminoglycans and glycosaminoglycan-degrading enzymes in the medical and biotechnology industries, and an important scientific basis for precision nutrition and medical intervention studies using gut microbiota or enzymes as targets.

Ultrafast energy transfer between π-stacked aromatic rings upon inner-valence ionization

Non-covalently bound aromatic systems are ubiquitous and govern the physicochemical properties of various organic materials. They are important to many phenomena of biological and technological relevance, such as protein folding, base-pair stacking in nucleic acids, molecular recognition and self-assembly, DNA-drug interactions, crystal engineering and organic electronics. Nevertheless, their molecular dynamics and chemical reactivity, particularly in electronic excited states, are not fully understood. Here, we observe intermolecular Coulombic decay in benzene dimers, (C₆H₆)₂-the simplest prototypes of noncovalent π-π interactions between aromatic systems. Intermolecular Coulombic decay is initiated by a carbon 2s vacancy state produced by electron-impact ionization and proceeds through ultrafast energy transfer between the benzene molecules. As a result, the dimer relaxes with the emission of a further low-energy electron (<10 eV) and a pair of C₆H₆⁺ cations undergoing Coulomb explosion. Coincident fragment-ion and electron momentum spectroscopy, accompanied by ab initio calculations, enables us to elucidate the dynamical details of this ultrafast relaxation process.

Unraveling the Role of π-Stacking Interactions in Ligand Binding to the Thiamine Pyrophosphate Riboswitch with High-level Quantum Chemical Calculations and Docking Study

The thiamine pyrophosphate (TPP) riboswitch has emerged as the new target for designing new ligands for antibiotic purpose. Binding of the natural ligand TPP to the TPP riboswitch causes downregulation of the genes responsible for its biosynthesis. We have reported the role of π-stacking energy contributions to ligand binding with a TPP riboswitch. In conjunction with the docking study, the higher-level quantum chemical calculations performed with the wB97XD and Def2TZVPP basis set in the aqueous phase revealed that the optimum ring size is crucial to attain the effective binding efficiency of ligands with a TPP riboswitch. The π-stacking energy contributions observed for the ligands studied are largely similar; however, the cases studied with higher π-stacking energies with larger rings have a weaker ability to displace the radiolabeled thiamine from the riboswitch. The EDA and NCI analyses suggest the role of larger dispersive interactions in stabilizing the π-stacking rings. The contribution from hydrogen-bonding interactions of the hydrogen-bond donor groups on the A ring augments the binding affinity of the ligand. This study sheds light on various factors that contribute to the design of new ligands for efficient binding with a TPP riboswitch and inhibition of gene expression.

Dynamical Behavior of Aromatic Trimer Complexes in Unimolecular Dissociation Reaction at High Temperatures. Case Studies on C6H6-C6F6-C6H6 and C6H6 Trimer Complexes

The intramolecular vibrational energy redistribution (IVR) dynamics during unimolecular dissociation of aromatic trimers at high temperatures is the primary interest of this study. Chemical dynamics simulations are performed for the unimolecular dissociation of benzene-hexafluorobenzene-benzene (Bz-HFB-Bz) and benzene trimer (Bz-trimer) complexes at a temperature range of 1000-2000 K. Partial dissociation of both the complexes is observed, which leads to a dimer and a monomer in the dynamics. However, the probability of such dissociation was found much lower in the case of the Bz-trimer, which further decreases with the increase of temperature. The rate of partial dissociation of Bz-HFB-Bz is faster at 1500, 1800, and 2000 K, whereas the rate of complete dissociation of the Bz-trimer is significantly faster than Bz-HFB-Bz at all temperatures. This is just the opposite of the corresponding dimer's dissociation, where benzene-hexafluorobenzene (Bz-HFB) dissociates at a faster rate than the benzene dimer (Bz-dimer). Thus, the dissociation dynamics of the trimer is different than that of the dimer. Simulations with excited intramolecular and intermolecular modes of the trimer complexes reveal that energy flows from intermolecular to intramolecular modes of Bz-HFB-Bz more freely than the Bz-trimer, and the dissociation process becomes slower for the former. Calculated activation energies for both types of dynamics are much lower than the corresponding binding energies, which may be due to the anharmonicity. The Arrhenius equation with an anharmonic correction factor is considered to recalculate the activation energy and pre-exponential factor.

Stacked but not Stuck: Unveiling the Role of π→π* Interactions with the Help of the Benzofuran-Formaldehyde Complex

The 1:1 benzofuran-formaldehyde complex has been chosen as model system for analyzing π→π* interactions in supramolecular organizations involving heteroaromatic rings and carbonyl groups. A joint "rotational spectroscopy-quantum chemistry" strategy unveiled the dominant role of π→π* interactions in tuning the intermolecular interactions of such adduct. The exploration of the intermolecular potential energy surface led to the identification of 14 low-energy minima, with 4 stacked isomers being more stable than those linked by hydrogen bond or lone-pair→π interactions. All energy minima are separated by loose transition states, thus suggesting an effective relaxation to the global minimum under the experimental conditions. This expectation has been confirmed by the experimental detection of only one species, which was unambiguously assigned owing to the computation of accurate spectroscopic parameters and the characterization of 11 isotopologues. The large number of isotopic species opened the way to the determination of the first semi-experimental equilibrium structure for a molecular complex of such a dimension.

Energetic and Topological Supramolecular Study and Nucleation Proposal of Halogenated Aryl Amides

A supramolecular study of 12 halogenated aryl amides (F, Cl, Br, I) in the ortho, meta, or para positions was realized. From mono crystal X-ray diffraction data and theoretical calculations...

Parallel-stacked aromatic molecules in hydrogen-bonded inorganic frameworks

By precisely constructing molecules and assembling these into well-defined supramolecular structures, novel physical properties and functionalities can be realized, and new areas of the chemical space can be accessed. In both materials science and biology, a deeper understanding of the properties and exploitation of the reversible character of weak bonds and interactions, such as hydrogen bonds and π–π interactions, is anticipated to lead to the development of materials with novel properties and functionalities. We apply the hydrogen-bonded organic frameworks (HOFs) strategy to inorganic materials science using the cubic octamer of orthosilicic acid, [Si₈O₁₂][OH]₈, as a building block, and find that various types of hydrogen-bonded inorganic frameworks (HIFs). We succeed in parallel π-stacking pure benzene, thiophene, selenophene, p-benzoquinone, thiophene·p-benzoquinone, and benzene·p-benzoquinone polymers infinitely. These polymers interact via their π-systems by taking advantage of the flexible pores of the three-dimensional nano-honeycomb HIFs, which consist of periodic wide and narrow segments.

Hydrogen-bonded inorganic frameworks are porous structures that may lead to novel materials with unprecedented properties and functionalities. Here the authors report the solid-state structures of orthosilicic acid-based hydrogen-bonded inorganic frameworks that can encapsulate small unsaturated cyclic molecules such as benzene, which are found stacked in parallel.

Structures, binding energies and non-covalent interactions of furan clusters

Understanding of the furan solvent is subjected to the knowledge of the structures of the furan clusters and interactions taking place therein. Although, furan clusters can be very important to determine the dynamics and the properties of the furan solvent, there has been only a few investigations reported on furan dimer. In this work, we have explored the potential energy surfaces (PESs) of the furan clusters using two incremental levels of theory. Structures have been initially generated using classical molecular dynamics followed by full optimization at the MP2/aug-cc-pVDZ level of theory. The results show that the most stable structure of the furan dimer has a stacking configuration while that of the trimer has a cyclic configuration. We have noted that the structures of the furan tetramer have no definite configurations. In addition, we have performed a quantum theory of atoms in molecule (QTAIM) analysis to identify all possible non-covalent interactions of the furan clusters. The results show that six different types of non-covalent interactions can be identified in furan clusters. We have noted that the CH⋯C and CH⋯O hydrogen bondings are the strongest non-covalent interactions while the H⋯H bonding interaction is found to be the weakest. Furthermore, we have assessed the performance of ten DFT functionals in calculating the binding energies of the furan clusters. The ten DFT functionals (M05, M05-2X, M06, M06-2X, M08HX, PBE0, ωB97XD, PW6B95D3, APFD and MN15) have been benchmarked to DLPNO-CCSD(T)/CBS. The functionals M05-2X and M06 are recommended for further affordable investigations of the furan clusters.

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.