The Stability of Hydrogen-Bonded Ion-Pair Complex Unexpectedly Increases with Increasing Solvent Polarity

The generally observed decrease of the electrostatic energy in the complex with increasing solvent polarity has led to the assumption that the stability of the complexes with ion-pair hydrogen bonds decreases with increasing solvent polarity. Besides, the smaller solvent-accessible surface area (SASA) of the complex in comparison with the isolated subsystems results in a smaller solvation energy of the latter, leading to a destabilization of the complex in the solvent compared to the gas phase. In our study, which combines Nuclear Magnetic Resonance, Infrared Spectroscopy experiments, quantum chemical calculations, and molecular dynamics (MD) simulations, we question the general validity of this statement. We demonstrate that the binding free energy of the ion-pair hydrogen-bonded complex between 2-fluoropropionic acid and n-butylamine (CH3CHFCOO-…NH3But+) increases with increased solvent polarity. This phenomenon is rationalized by a substantial charge transfer between the subsystems that constitute the ion-pair hydrogen-bonded complex. This unexpected finding introduces a new perspective to our understanding of solvation dynamics, emphasizing the interplay between solvent polarity and molecular stability within hydrogen-bonded systems.

Many-body Effects on Electronic Transport in Molecular Junctions: A Quantum Perspective

This concept delves into quantum particle transport at the nanoscale, with a particular focus on how electrons move through molecular circuits. The thriving field of single molecular electronics benefits from the unique electrical and other properties of nanostructures. It concentrates on single molecular junctions that serve as bridges between electrodes. In this context, the electronic correlation-induced many-body effect gives rise to resonant states. These states, along with conductance, depend on electron spin. Thus, the field acts as a bridge between quantum and macroscopic worlds, unveiling unique behaviors of electrons. Additionally, external factors, such as magnetic fields and voltages, offer means to control the electron correlation in these junctions.

Impact of dielectric constant of solvent on the formation of transition metal-ammine complexes

The DFT-level computational investigations into Gibbs free energies (ΔG) demonstrate that as the dielectric constant of the solvent increases, the stabilities of [M(NH3 )n ]2+/3+ (n = 4, 6; M = selected 3d transition metals) complexes decrease. However, there is no observed correlation between the stability of the complex and the solvent donor number. Analysis of the charge transfer and Wiberg bond indices indicates a dative-bond character in all the complexes. The solvent effect assessed through solvation energy is determined by the change in the solvent accessible surface area (SASA) and the change in the charge distribution that occurs during complex formation. It has been observed that the SASA and charge transfer are different in the different coordination numbers, resulting in a variation in the solvent effect on complex stability in different solvents. This ultimately leads to a change between the relative stability of complexes with different coordination numbers while increasing the solvent polarity for a few complexes. Moreover, the findings indicate a direct relationship between ΔΔG (∆Gsolvent -∆Ggas ) and ΔEsolv , which enables the computation of ΔG for the compounds in a particular solvent using only ΔGgas and ΔEsolv . This approach is less computationally expensive.

Trends in the stability of covalent dative bonds with variable solvent polarity depend on the charge transfer in the Lewis electron-pair system

In general, the stability of neutral complexes with dative bonds increases as the polarity of the solvent increases. This is based on the fact that the dipole moment of the complex increases as the charge transferred from the donor to the acceptor increases. As a result, the solvation energy of the complex becomes greater than that of subsystems, causing an increase in the stabilization energy with increasing solvent polarity. Our research confirms this assumption, but only when the charge transfer is sufficiently large. If it is below a certain threshold, the increase in the complex's dipole moment is insufficient to result in a higher solvation energy than subsystems. Thus, the magnitude of the charge transfer in the Lewis electron-pair system determines the stability trends of dative bonds with varying solvent polarity. We used molecular dynamics (MD) simulations based on an explicit solvent model, which is considered more reliable, to verify the results obtained with a continuous solvent model.

Strong Be-N Interaction Induced Complementary Chemical Tuning to Design a Dual-gated Single Molecule Junction

The interaction between pyridines and the π-hole of BeH2 leads to the formation of strong beryllium-bonded complexes. Theoretical investigations demonstrate that the Be-N bonding interaction can effectively regulate the electronic current through a molecular junction. The electronic conductance exhibits distinct switching behavior depending on the substituent groups at the para position of pyridine, highlighting the role of Be-N interaction as a potent chemical gate in the proposed device. The complexes exhibit short intermolecular distances ranging from 1.724 to 1.752 Å, emphasizing their strong binding. Detailed analysis of electronic rearrangements and geometric perturbations upon complex formation provides insights into the underlying reasons for the formation of such strong Be-N bonds, with bond strengths varying from -116.25 to -92.96 kJ/mol. Moreover, the influence of chemical substituents on the local electronic transmission of the beryllium-bonded complex offers valuable insights for the implementation of a secondary chemical gate in single-molecule devices. This study paves the way for the development of chemically gateable, functional single-molecule transistors, advancing the design and fabrication of multifunctional single-molecule devices in the nanoscale regime.

Visualization of π-hole in molecules by means of Kelvin probe force microscopy

Submolecular charge distribution significantly affects the physical-chemical properties of molecules and their mutual interaction. One example is the presence of a π-electron-deficient cavity in halogen-substituted polyaromatic hydrocarbon compounds, the so-called π-holes, the existence of which was predicted theoretically, but the direct experimental observation is still missing. Here we present the resolution of the π-hole on a single molecule using the Kelvin probe force microscopy, which supports the theoretical prediction of its existence. In addition, experimental measurements supported by theoretical calculations show the importance of π-holes in the process of adsorption of molecules on solid-state surfaces. This study expands our understanding of the π-hole systems and, at the same time, opens up possibilities for studying the influence of submolecular charge distribution on the chemical properties of molecules and their mutual interaction.

The Impact of the Solvent Dielectric Constant on A← NH3 Dative Bond Depends on the Nature of the Lewis Electron-Pair Systems

The present work aims to determine to what extent the value of the dielectric constant of the solvent can influence the dative bond in Lewis electron pair bonding systems. For this purpose, two different systems, namely H3B¬NH3 and {Zn¬(NH3)}2+, were studied in selected solvents with significantly different dielectric constants. Based on the results resulting from state-of-the-art computational methods using DFT, constrained DFT, energy decomposition analyses, solvent accessible surface area, and charge transfer calculations, we found that the stability of the neutral H3B¬NH3 system increases with increasing solvent polarity. In contrast, the opposite trend is observed for the positively charged {Zn¬(NH3)}2+. The observed changes are attributed to different charge redistributions in neutral and charged complexes, which are reflected by a different response to the solvent and are quantified by changes in solvation energies.

Hydrogen Bonding with Hydridic Hydrogen-Experimental Low-Temperature IR and Computational Study: Is a Revised Definition of Hydrogen Bonding Appropriate?

Spectroscopic characteristics of Me₃Si-H···Y complexes (Y = ICF₃, BrCN, and HCN) containing a hydridic hydrogen were determined experimentally by low-temperature IR experiments based on the direct spectral measurement of supersonically expanded intermediates on a cold substrate or by the technique of argon-matrix isolation as well as computationally at harmonic and one-dimensional anharmonic levels. The computations were based on DFT-D, MP2, MP2-F12, and CCSD(T)-F12 levels using various extended AO basis sets. The formation of all complexes related to the redshift of the Si-H stretching frequency upon complex formation was accompanied by an increase in its intensity. Similar results were obtained for another 10 electron acceptors of different types, positive σ-, π-, and p-holes and cations. The formation of HBe-H···Y complexes, studied only computationally and again containing a hydridic hydrogen, was characterized by the blueshift of the Be-H stretching frequency upon complexation accompanied by an increase in its intensity. The spectral shifts and stabilization energies obtained for all presently studied hydridic H-bonded complexes were comparable to those in protonic H-bonded complexes, which has prompted us to propose a modification of the existing IUPAC definition of H-bonding that covers, besides the classical protonic form, the non-classical hydridic and dihydrogen forms.

The unusual stability of H-bonded complexes in solvent caused by greater solvation energy of complex compared to those of isolated fragments

Here, the effect of solvent on the stability of non-covalent complexes, was studied. These complexes were from previously published S22, S66, and X40 datasets, which include hydrogen-, halogen- and dispersion-bonded complexes. It was shown that the charge transfer in the complex determines whether the complex is stabilized or destabilized in solvent.

Unexpected Strengthening of the H-Bond Complexes in a Polar Solvent Due to a More Efficient Solvation of the Complex Compared to Isolated Monomers

It is generally assumed that hydrogen-bonded complexes are less stable in solvents than in the gas phase and that their stability decreases with increasing solvent polarity. This assumption is based on the size of the area available to the solvent, which is always smaller in the complex compared to the subsystems, thereby reducing the solvation energy. This reduction prevails over the amplification of the electrostatic hydrogen bond by the polar solvent. In this work, we show, using experimental IR spectroscopy and DFT calculations, that there are hydrogen-bonded complexes whose stability becomes greater with increasing solvent polarity. The explanation for this surprising stabilization is based on the analysis of the charge redistribution in the complex leading to increase of its dipole moment and solvation energy. Constrained DFT calculations have shown a dominant role of charge transfer over polarization effects for dipole moment and solvation energy.

P-Doped graphene-C60 nanocomposite: a donor-acceptor complex with a P-C dative bond

Phosphorous-doped graphene can form a covalent dative bond with the electron acceptor, C₆₀ molecule. On the other hand, C₆₀ on graphene and N-doped graphene surfaces can only form vdW complexes. State-of-the-art quantum-chemical techniques have been used to characterise the nature of the P-C dative bond. A considerable amount of charge transfer from the P-Gr surface to C₆₀ has been observed. This complex formation may enable enhancement in the optical limiting response with potential application in energy harvesting. The stability of the P-C dative bond has been assessed using DFT-D molecular dynamics simulations at 300 K for 10 ps.

Real-space imaging of anisotropic charge of σ-hole by means of Kelvin probe force microscopy

[Figure: see text].

The Structure-Based Design of SARS-CoV-2 nsp14 Methyltransferase Ligands Yields Nanomolar Inhibitors

In this study, we have focused on the structure-based design of the inhibitors of one of the two SARS-CoV-2 methyltransferases (MTases), nsp14. This MTase catalyzes the transfer of the methyl group from S-adenosyl-l-methionine (SAM) to cap the guanosine triphosphate moiety of the newly synthesized viral RNA, yielding the methylated capped RNA and S-adenosyl-l-homocysteine (SAH). As the crystal structure of SARS-CoV-2 nsp14 is unknown, we have taken advantage of its high homology to SARS-CoV nsp14 and prepared its homology model, which has allowed us to identify novel SAH derivatives modified at the adenine nucleobase as inhibitors of this important viral target. We have synthesized and tested the designed compounds in vitro and shown that these derivatives exert unprecedented inhibitory activity against this crucial enzyme. The docking studies nicely explain the contribution of an aromatic part attached by a linker to the position 7 of the 7-deaza analogues of SAH.

Addition Reaction between Piperidine and C60 to Form 1,4-Disubstituted C60 Proceeds through van der Waals and Dative Bond Complexes: Theoretical and Experimental Study

A combined computational and experimental study reveals the character of the C₆₀ complexes with piperidine formed under different reaction conditions. The IR and NMR experiments detect the dative bond complex, which according to NMR, is stable in the oxygen-free environment and transforms to the adduct complex in the presence of O₂. Computational studies on the character of reaction channels rationalize the experimental observations. They show that the piperidine dimer rather than a single piperidine molecule is required for the complex formation. The calculations reveal significant differences in the dative bond and adduct complexes' character, suggesting a considerable versatility in their electronic properties modulated by the environment. This capability offers new application potential in several fields, such as in energy storage devices.

Cyclo[n]carbons Form Strong N → C Dative/Covalent Bonds with Piperidine

The newly synthesized C₁₈ ring is demonstrated as the smallest all-carbon acceptor that exhibits strong electron acceptance. This study provides a quantum-chemical investigation of the electron-acceptance behavior of monocyclic carbon rings with a particular emphasis on C₁₈ through the formation of a dative bond with piperidine. The results show that C_n rings form strong dative bonds with piperidine, whereas the respective van der Waals (vdW) complexes are higher in energy. The main driving force is the release of angle strain of cyclo[n]carbons caused by the change in hybridization from sp to sp² associated with the formation of the dative bond. On the contrary, other sp allotropes, diynes, favorably form vdW complexes. Molecular dynamics (MD) simulations support the stability of the dative bond throughout a simulation of 20 ps. This opens up the possibility of stabilizing highly reactive C₁₈ through dative/covalent functionalization.

Tuning the P-C dative/covalent bond formation in R3P-C60 complexes by changing the R group

The P-C dative/covalent bonds formed in R₃P-C₆₀ complexes (R = OCH₃, N(CH₃)₂, NC₄H₈) have been affected by the nature of the R group. The highest stabilisation (18.7 kcal mol^-1) has been found in the last system. The contribution of dispersion energies in the stabilisation also varies depending on the R group. The nature of the P→C bond has been characterised using state-of-the-art quantum-chemical techniques including NBO, AIM and ELF. The P→C dative bond is significantly different from the prototype dative bonds appearing in H₃N→BH₃ as well as in the fullerene - secondary-amine complexes previously studied by us. The findings obtained through electron structure theory have been supported by 10 ps DFT-D MD simulations.

Structure-directed formation of the dative/covalent bonds in complexes with C70piperidine

The combined experimental-computational study has been performed to investigate the complexes formed between C70 carbon allotrope and piperidine. The results of FT-IR, H-NMR, and C-NMR measurements, together with the calculations based on the DFT approach and molecular dynamics simulations, prove the existence of dative/covalent bonding in C70piperidine complexes. The dative bond forms not only at the region of five- and six-membered rings, observed previously with C60, but also at the region formed of six-membered rings. The structure, i.e., nonplanarity, explains the observed dative bond formation. New findings on the character of interaction of secondary amines with C70 bring new aspects for the rational design of modified fullerenes and their applications in electrocatalysis, spintronics, and energy storage.

Carboxylated Graphene for Radical-Assisted Ultra-Trace-Level Water Treatment and Noble Metal Recovery

Sorption technologies, enabling removal of heavy metals, play a pivotal role in meeting the global demands for unrestricted access to drinking water. Standard sorption technologies suffer from limited efficiency related to the weak sorbent-metal interaction. Further challenges include the development of technologies enabling smart metal recovery and sorbent regeneration. To this end, a densely functionalized graphene, with 33% by mass content of carboxyl groups, linked through direct C-C bonds (graphene acid, GA) represents a previously unexplored solution to this challenge. GA revealed excellent efficiency for removal of highly toxic metals, such as Cd²⁺ and Pb²⁺. Due to its selective chemistry, GA can bind heavy metals with high affinity, even at concentrations of 1 mg L^-1 and in the presence of competing ions of natural drinking water, and reduce them down to drinking water allowance levels of a few μg L^-1. This is not only due to carboxyl groups but also due to the stable radical centers of the GA structure, enabling metal ion-radical interactions, as proved by EPR, XPS, and density functional theory calculations. GA offers full structural integrity during the highly acidic and basic treatment, which is exploited for noble metal recovery (Ga³⁺, In³⁺, Pd²⁺) and sorbent regeneration. Owing to these attributes, GA represents a fully reusable metal sorbent, applicable also in electrochemical energy technologies, as illustrated with a GA/Pt catalyst derived from Pt⁴⁺-contaminated water.

Benchmark Data Sets of Boron Cluster Dihydrogen Bonding for the Validation of Approximate Computational Methods

The success of approximate computational methods, such as molecular mechanics, or dispersion-corrected density functional theory, in the description of non-covalent interactions relies on accurate parameterizations. Benchmark data sets are thus required. This area is well developed for organic molecules and biomolecules but practically non-existent for boron clusters, which have been gaining in importance in modern drug as well as material design. To fill this gap, we have introduced two data sets featuring the most common non-covalent interaction of boron clusters, the dihydrogen bond, and calculated reference interaction energies at the "golden standard" CCSD(T)/CBS level. The boron clusters studied interact with formamide, methanol, water and methane at various distances and in two geometrical arrangements. The performance of the tested approximate methods is variable and recommendations for further use are given.

SQM/COSMO Scoring Function: Reliable Quantum-Mechanical Tool for Sampling and Ranking in Structure-Based Drug Design

Invited for this month's cover is the group of Prof. Pavel Hobza, Czech Academy of Sciences, Prague. The cover picture shows a powerful automated quantum mechanics based SQM/COSMO approach to protein-ligand scoring. It comprises thorough preparation of ligand structures, extensive generation of binding complexes, fast geometry relaxation and reliable affinity prediction. Read the full text of the Minireview at 10.1002/cplu.202000120.

Chelating Polymers for Hereditary Hemochromatosis Treatment

Hemochromatosis (iron overload) encompasses a group of diseases that are characterized by a toxic hyperaccumulation of iron in parenchymal organs. Currently, only few treatments for this disease have been approved; however, all these treatments possess severe side effects. In this study, a paradigm for hemochromatosis maintenance/preventive therapy is investigated: polymers with negligible systemic biological availability form stable complexes with iron ions in the gastrointestinal tract, which reduces the biological availability of iron. Macroporous polymer beads are synthesized with three different iron-chelating moieties (benzene-1,2-diol, benzene-1,2,3-triol, and 1,10-phenanthroline). The polymers rapidly chelate iron ions from aqueous solutions in vitro in the course of minutes, and are noncytotoxic and nonprooxidant. Moreover, the in vivo biodistribution and pharmacokinetics show a negligible uptake from the gastrointestinal tract (using ¹²⁵ I-labeled polymer and single photon emission computed tomography/computed tomography), which generally prevents them from having systemic side effects. The therapeutic efficacy of the prepared polymers is successfully tested in vivo, and exhibits a significant inhibition of iron uptake from the gastrointestinal tract without any noticeable signs of toxicity. Furthermore, an in silico method is developed for the prediction of chelator selectivity. Therefore, this paradigm can be applied to the next-generation maintenance/preventive treatment for hemochromatosis and/or other diseases of similar pathophysiology.

Mechanical force-induced manipulation of electronic conductance in a spin-crossover complex: a simple approach to molecular electronics

The atomic-scale technological sophistication from the last half-decade provides new avenues for the atom-by-atom fabrication of nanostructures with extraordinary precision. This urges the appraisal of the fabrication scheme layout for a modular nanoelectronic device based on an individual molecular complex. The mechanical force-induced distortion to the metal coordination sphere triggers a low-spin (LS) to high-spin (HS) electronic transition in the complex. The controlled structural distortions (relative to a specific bond-angle) are deemed to be the switching parameter for the observed spin-transitions. Mechanical stretching is the key to engineering a spin-state switch in the proposed molecular device. The spin-dependent reversible variation in the electronic conductance concurrent to the unique spin-states can be understood from the state-of-the-art Nonequilibrium Green's Function (NEGF) calculations. Combined with NEGF calculations, the DFT study further provides a qualitative perception of the electronic conductance in the two-terminal device architecture. From the transport calculations, there is also evidence of considerable fluctuation in the spin-dependent electronic conductance at the molecular junction with relative variations in the scattering limit. Subsequently, the present study shows significant advances in the transmission probabilities for the high-spin state of the Fe(ii) complex. The results empower the progress of nanoelectronics at the single molecule level.

SQM/COSMO Scoring Function: Reliable Quantum-Mechanical Tool for Sampling and Ranking in Structure-Based Drug Design

Quantum mechanical (QM) methods have been gaining importance in structure-based drug design where a reliable description of protein-ligand interactions is of utmost significance. However, strategies i. e. QM/MM, fragmentation or semiempirical (SQM) methods had to be pursued to overcome the unfavorable scaling of QM methods. Various SQM-based approaches have significantly contributed to the accuracy of docking and improvement of lead compounds. Parametrizations of SQM and implicit solvent methods in our laboratory have been instrumental to obtain a reliable SQM-based scoring function. The experience gained in its application for activity ranking of ligands binding to tens of protein targets resulted in setting up a faster SQM/COSMO scoring approach, which outperforms standard scoring methods in native pose identification for two dozen protein targets with ten thousand poses. Recently, SQM/COSMO was effectively applied in a proof-of-concept study of enrichment in virtual screening. Due to its superior performance, feasibility and chemical generality, we propose the SQM/COSMO approach as an efficient tool in structure-based drug design.

Ground state of the Fe(ii)-porphyrin model system corresponds to quintet: a DFT and DMRG-based tailored CC study

Fe(ii)-porphyrins play an important role in many reactions, due to their closely lying spin states. We present a thorough study of a Fe(ii)-porphyrin model system, in which we examine how the geometrical parameters influence its spin state ordering.

Spin modification of iron(ii) complexes via covalent (dative) and dispersion guided non-covalent bonding with N-heterocyclic carbenes: DFT, DLPNO-CCSD(T) and MCSCF studies

Spin crossover from high spin Fe(ii)-phthalocyanine to low or intermediate spin via either dative covalent or non-covalent interaction by just varying the substituent using the same core ligand.

Mutations at hypothetical binding site 2 in insulin and insulin-like growth factors 1 and 2 result in receptor- and hormone-specific responses

Information on how insulin and insulin-like growth factors 1 and 2 (IGF-1 and -2) activate insulin receptors (IR-A and -B) and the IGF-1 receptor (IGF-1R) is crucial for understanding the difference in the biological activities of these peptide hormones. Cryo-EM studies have revealed that insulin uses its binding sites 1 and 2 to interact with IR-A and have identified several critical residues in binding site 2. However, mutagenesis studies suggest that Ile-A10, Ser-A12, Leu-A13, and Glu-A17 also belong to insulin's site 2. Here, to resolve this discrepancy, we mutated these insulin residues and the equivalent residues in IGFs. Our findings revealed that equivalent mutations in the hormones can result in differential biological effects and that these effects can be receptor-specific. We noted that the insulin positions A10 and A17 are important for its binding to IR-A and IR-B and IGF-1R and that A13 is important only for IR-A and IR-B binding. The IGF-1/IGF-2 positions 51/50 and 54/53 did not appear to play critical roles in receptor binding, but mutations at IGF-1 position 58 and IGF-2 position 57 affected the binding. We propose that IGF-1 Glu-58 interacts with IGF-1R Arg-704 and belongs to IGF-1 site 1, a finding supported by the NMR structure of the less active Asp-58–IGF-1 variant. Computational analyses indicated that the aforementioned mutations can affect internal insulin dynamics and inhibit adoption of a receptor-bound conformation, important for binding to receptor site 1. We provide a molecular model and alternative hypotheses for how the mutated insulin residues affect activity.

Impressive Enrichment of Semiempirical Quantum Mechanics-Based Scoring Function: HSP90 Protein with 4541 Inhibitors and Decoys

This paper describes the excellent performance of a newly developed scoring function (SF), based on the semiempirical QM (SQM) PM6-D3H4X method combined with the conductor-like screening implicit solvent model (COSMO). The SQM/COSMO, Amber/GB and nine widely used SFs have been evaluated in terms of ranking power on the HSP90 protein with 72 biologically active compounds and 4469 structurally similar decoys. Among conventional SFs, the highest early and overall enrichment measured by EF₁ and AUC% obtained using single-scoring-function ranking has been found for Glide SP and Gold-ASP SFs, respectively (7, 75 % and 3, 76 %). The performance of other standard SFs has not been satisfactory, mostly even decreasing below random values. The SQM/COSMO SF, where P-L structures were optimised at the advanced Amber level, has resulted in a dramatic enrichment increase (47, 98 %), almost reaching the best possible receiver operator characteristic (ROC) curve. The best SQM frame thus inserts about seven times more active compounds into the selected dataset than the best standard SF.

Non-covalent control of spin-state in metal-organic complex by positioning on N-doped graphene

Nitrogen doping of graphene significantly affects its chemical properties, which is particularly important in molecular sensing and electrocatalysis applications. However, detailed insight into interaction between N-dopant and molecules at the atomic scale is currently lacking. Here we demonstrate control over the spin state of a single iron(II) phthalocyanine molecule by its positioning on N-doped graphene. The spin transition was driven by weak intermixing between orbitals with z-component of N-dopant (pz of N-dopant) and molecule (dxz, dyz, dz2) with subsequent reordering of the Fe d-orbitals. The transition was accompanied by an electron density redistribution within the molecule, sensed by atomic force microscopy with CO-functionalized tip. This demonstrates the unique capability of the high-resolution imaging technique to discriminate between different spin states of single molecules. Moreover, we present a method for triggering spin state transitions and tuning the electronic properties of molecules through weak non-covalent interaction with suitably functionalized graphene.

Understanding the non-covalent interaction mediated modulations on the electronic structure of quasi-zero-dimensional graphene nanoflakes

In recent years, magnetic or electric field induced modulations on the electronic environment of single molecular systems are common practice. In this particular study, we have instigated the possibility of controlling the electronic and spin-dependent properties of hydrogen-terminated graphene fragments, so-called graphene nanoflakes (GNF), using weak non-covalent interactions as the external stimuli. The topological frustration in the graphene fragment appreciated the compelling electronic behavior of the system. This leads to some unorthodox spin-distribution in the system and it is possible to synchronize this electronic perturbation switching through a non-covalent interaction. These findings institute a new avenue for sculpting such donor-acceptor composites as self-regulated spintronic devices in next generation electronics.

Comparison of the DFT-SAPT and Canonical EDA Schemes for the Energy Decomposition of Various Types of Noncovalent Interactions

Interaction energies computed with density functional theory can be divided into physically meaningful components by symmetry-adapted perturbation theory (DFT-SAPT) or the canonical energy decomposition analysis (EDA). In this work, the decomposition results obtained by these schemes were compared for more than 200 hydrogen-, halogen-, and pnicogen-bonded, dispersion-bound, and mixed complexes to investigate their similarity in the evaluation of the nature of noncovalent interactions. BLYP functional with D3(BJ) correction was used for the EDA scheme, whereas asymptotically corrected PBE0 functional for DFT-SAPT provided some of the best combinations for description of noncovalent interactions. Both schemes provide similar results concerning total interaction energies and insight into the individual energy components. For most complexes, the dominant energetic term was identified equally by both decomposition schemes. Because the canonical EDA is computationally less demanding than the DFT-SAPT, the former can be especially used in cases where the systems investigated are very large.

Ranking Power of the SQM/COSMO Scoring Function on Carbonic Anhydrase II-Inhibitor Complexes

Accurate prediction of protein-ligand binding affinities is essential for hit-to-lead optimization and virtual screening. The reliability of scoring functions can be improved by including quantum effects. Here, we demonstrate the ranking power of the semiempirical quantum mechanics (SQM)/implicit solvent (COSMO) scoring function by using a challenging set of 10 inhibitors binding to carbonic anhydrase II through Zn²⁺ in the active site. This new dataset consists of the high-resolution (1.1-1.4 Å) crystal structures and experimentally determined inhibitory constant (K_i ) values. It allows for evaluation of the common approximations, such as representing the solvent implicitly or by using a single target conformation combined with a set of ligand docking poses. SQM/COSMO attained a good correlation of R² of 0.56-0.77 with the experimental inhibitory activities, benefiting from careful handling of both noncovalent interactions (e.g. charge transfer) and solvation. This proof-of-concept study of SQM/COSMO ranking for metalloprotein-ligand systems demonstrates its potential for hit-to-lead applications.

Sequential BN-doping induced tuning of electronic properties in zigzag-edged graphene nanoribbons: a computational approach

We employed first-principles methods to elaborate doping induced electronic and magnetic perturbations in one-dimensional zigzag graphene nanoribbon (ZGNR) superlattices. Consequently, the incorporation of alternate boron and nitrogen (hole–electron) centers into the hexagonal network instituted substantial modulations to electronic and magnetic properties of ZGNR. Our theoretical analysis manifested some controlled changes to electronic and magnetic properties of the ZGNR by tuning the positions (array) of impurity centers in the carbon network. Subsequent DFT based calculations also suggested that the site-specific alternate electron–hole (B/N) doping could regulate the band-gaps of the superlattices within a broad range of energy. The consequence of variation in the width of ZGNR in the electronic environment of the system was also tested. The systematic analysis of various parameters such as the structural orientations, spin-arrangements, the density of states (DOS), band structures, and local density of states envisioned a basis for the band-gap engineering in ZGNR and attributed to its feasible applications in next generation electronic device fabrication.

Incorporation of an alternate impurity array in the ZGNR don't break the spin-degeneracy, providing the freedom to tune the electronic behaviors without affecting the spin-dependent properties.

The role of the σ-holes in stability of non-bonded chalcogenide⋯benzene interactions: the ground and excited states

Different characters of molecular electrostatic potential (MEP) in the ground and excited states of chalcogenides are responsible for changes in conformer stability of T-shape and stacked non-bonded chalcogenide–benzene complexes upon electronic excitation.

Various types of non-covalent interactions contributing towards crystal packing of halogenated diphospha-dicarbaborane with an open pentagonal belt

We have synthesized and crystalized 3-Cl-10-I-nido-7,8,9,11-P₂C₂B₇H₇. Quantum chemical calculations have demonstrated that the obtained crystal structure is stabilized by hydrogen, dihydrogen and pnictogen bonds.

Selective binding of choline by a phosphate-coordination-based triple helicate featuring an aromatic box

In nature, proteins have evolved sophisticated cavities tailored for capturing target guests selectively among competitors of similar size, shape, and charge. The fundamental principles guiding the molecular recognition, such as self-assembly and complementarity, have inspired the development of biomimetic receptors. In the current work, we report a self-assembled triple anion helicate (host 2) featuring a cavity resembling that of the choline-binding protein ChoX, as revealed by crystal and density functional theory (DFT)-optimized structures, which binds choline in a unique dual-site-binding mode. This similarity in structure leads to a similarly high selectivity of host 2 for choline over its derivatives, as demonstrated by the NMR and fluorescence competition experiments. Furthermore, host 2 is able to act as a fluorescence displacement sensor for discriminating choline, acetylcholine, l-carnitine, and glycine betaine effectively.

The choline-binding protein ChoX exhibits a synergistic dual-site binding mode that allows it to discriminate choline over structural analogues. Here, the authors design a biomimetic triple anion helicate receptor whose selectivity for choline arises from a similar binding mechanism.

Superior Performance of the SQM/COSMO Scoring Functions in Native Pose Recognition of Diverse Protein-Ligand Complexes in Cognate Docking

General and reliable description of structures and energetics in protein-ligand (PL) binding using the docking/scoring methodology has until now been elusive. We address this urgent deficiency of scoring functions (SFs) by the systematic development of corrected semiempirical quantum mechanical (SQM) methods, which correctly describe all types of noncovalent interactions and are fast enough to treat systems of thousands of atoms. Two most accurate SQM methods, PM6-D3H4X and SCC-DFTB3-D3H4X, are coupled with the conductor-like screening model (COSMO) implicit solvation model in so-called "SQM/COSMO" SFs and have shown unique recognition of native ligand poses in cognate docking in four challenging PL systems, including metalloprotein. Here, we apply the two SQM/COSMO SFs to 17 diverse PL complexes and compare their performance with four widely used classical SFs (Glide XP, AutoDock4, AutoDock Vina, and UCSF Dock). We observe superior performance of the SQM/COSMO SFs and identify challenging systems. This method, due to its generality, comparability across the chemical space, and lack of need for any system-specific parameters, gives promise of becoming, after comprehensive large-scale testing in the near future, a useful computational tool in structure-based drug design and serving as a reference method for the development of other SFs.

Adsorption of Organic Molecules to van der Waals Materials: Comparison of Fluorographene and Fluorographite with Graphene and Graphite

Understanding strength and nature of noncovalent binding to surfaces imposes significant challenge both for computations and experiments. We explored the adsorption of five small nonpolar organic molecules (acetone, acetonitrile, dichloromethane, ethanol, ethyl acetate) to fluorographene and fluorographite using inverse gas chromatography and theoretical calculations, providing new insights into the strength and nature of adsorption of small organic molecules on these surfaces. The measured adsorption enthalpies on fluorographite range from -7 to -13 kcal/mol and are by 1-2 kcal/mol lower than those measured on graphene/graphite, which indicates higher affinity of organic adsorbates to fluorographene than to graphene. The dispersion-corrected functionals performed well, and the nonlocal vdW DFT functionals (particularly optB86b-vdW) achieved the best agreement with the experimental data. Computations show that the adsorption enthalpies are controlled by the interaction energy, which is dominated by London dispersion forces (∼70%). The calculations also show that bonding to structural features, like edges and steps, as well as defects does not significantly increase the adsorption enthalpies, which explains a low sensitivity of measured adsorption enthalpies to coverage. The adopted Langmuir model for fitting experimental data enabled determination of adsorption entropies. The adsorption on the fluorographene/fluorographite surface resulted in an entropy loss equal to approximately 40% of the gas phase entropy.