Energy Decomposition Analysis Research Articles

The role of linkers and frustrated lewis pairs catalysts in the formation of zwitterionic 1,2-anti-addition product with non-conjugated terminal diacetylenes: A computational study

This study presents a computational investigation into the mechanistic pathway and the linker units involved in forming the zwitterionic 1,2-anti-addition product of non-conjugated diacetylenes, di(propargyl)ether (DPE), di(prop-2yn-1yl)sulfane (DPS) and 1,6-Heptadiyne (HD) catalyzed by the inter-molecular phosphine/borane frustrated Lewis pairs (FLPs), i.e., PPh2[C6H3(CF3)2](P-CF)/[B(C6F5)3]([B]) and P(o-tolyl)3(P-tol)/[B(C6F5)3]([B]). The potential energy surface (PES) calculations reveal that the anti-addition of P-CF to the internal C-atoms of acetylene units is energetically more favored than that of the addition of P-tol in DPE, DPS, and HD by ∼10.0, ∼9.2, and ∼6.0 kcal/mol, respectively. The calculations performed with DPE contain “—O—,” linker unit exhibits superior reactivity than DPS and HD, which suggests the electronegativity of linkers plays a significant role and facilitates the addition of Lewis bases. The higher electronegativity of linker units enables the 1,2-addition reaction by lowering the free energy activation barriers, as observed in the DFT calculations. The Molecular Electrostatic Potential (MESP) study shows that the electrostatic interactions favor the addition of P-CF to the active acetylene positions (C5/C4/C4) of [B]-DPE/DPS/HD-π complexes than the P-tol. The Distortion/Interaction (D/I) analysis reveals that transition states involving P-CF (TS1, TS3, and TS5) exhibit more interaction energy (ΔEInt) and less distortion energies (ΔEd) than that of the P-tol (TS2, TS4, and TS6). Further, the Energy Decomposition Analysis (EDA) also rationalizes the preferential approach of the electron-deficient Lewis base over the electron-rich one on the basis of the significant contribution of orbital interaction energies (ΔEorbital) in the cases of P-CF; TS1, TS3, and TS5. This study suggests that the electronic effects of substrates and the FLPs are crucial to facilitate the desired products formed with non-conjugated terminal alkynes.

Cuprophilic Interactions in Polymeric [Cu10O2(Mes)6]n.

The properties of cuprophilic compounds and the underlying fundamental principles responsible for the Cu(I)···Cu(I) interactions have been the subject of intense research as their diverse structural and physical attributes are being explored. In this light, we performed a new study of the compound [Cu10O2(Mes)6] reported by Haakansson et al. using state of the art experimental and theoretical analysis techniques. Doing this, we found the compound to be a polymer in the solid state, best written as [Cu10O2(Mes)6]n, with unsupported Cu(I)···Cu(I) contacts linking the monomers (2.776 Å). The monomeric unit also exhibits various cuprophilic contacts bridged by mesityl and/or oxo ligands. The compound was analyzed in its solid state, revealing luminescent properties resulting from two distinct fluorescent emissions, as well as in solution, in which its polymeric structure reversibly decomposes. A quantum theory of atoms in molecules (QTAIM) analysis based on density functional theory (DFT) calculations allows to characterize the various Cu(I)···Cu(I) contacts, in which only a few, and not necessarily the shortest, are associated with a bond critical point. Additionally, an energy decomposition analysis of the bonding between monomers indicates that it is dominated by dispersion forces in which the ligands play a dominant role, resulting in bonding energies significantly larger than found in previous DFT investigations based on less bulky models.

Synthetic and theoretical study on novel N-ylidic complexes of mercury as new antibacterial agents

The novel structure of Hg(II) complexes including the pyridinium ylide C5H5NCHC(O)C6H4-m-Br (Y) were synthesized and reported in this study. In the first step, the pyridinium salt C5H5NCH2C(O)C6H4-m-Br (S) was produced by reacting 2,3′-dibromoactophenone and pyridine. then, treatment of S with K2CO3 gave the related pyridinium ylide Y. Finally, the reaction of Y with HgX2 and Hg(NO3)2·H2O leads to the formation of novel binuclear [HgY2][HgX4] (X=Cl (1); X=Br (2); X=I (3)) and polymeric [HgY(NO3)2]n (4) complexes. The structure of complex 2 was also determined by X-ray diffraction analysis. The obtained analyses proved the coordination through the ylidic carbon to metallic center. Additionally, Natural Bond Orbital (NBO), Energy Decomposition Analysis (EDA), and EDA-NOCV studies are also used to investigate the nature of metal–ligand bonding in the complexes. Finally, the antibacterial activity of 1–4 was also examined against Gram positive and negative represented significant levels of inhibitory potency respected to used standards.

Copper-Catalyzed Asymmetric Hydrogenation of Unsymmetrical ortho-Br Substituted Benzophenones.

The asymmetric hydrogenation of benzophenones, catalyzed by low-activity earth-abundant metal copper, has hitherto remained a challenge due to the substrates equipped with two indistinguishably similar aryl groups. In this study, we demonstrated that the prochiral carbon of the ortho-bromine substrate exhibits the highest electrophilicity and high reactivity among the ortho-halogen substituted benzophenones, as determined by the Fukui function (f+) analysis and hydrogenation reaction. Considering that the enantiodirecting functional bromine group can be easily derivatized and removed in the products, we successfully achieved a green copper-catalyzed asymmetric hydrogenation of ortho-bromine substituted benzophenones. This method yielded a series of chiral benzhydrols with excellent results. The utility of this protocol has been validated through a gram-scale reaction and subsequent product transformations. Hirshfeld partition (IGMH) and energy decomposition analysis (EDA) indicate that the CH···HC multiple attractive dispersion interactions (MADI) effect between the catalyst and substrate enhances the catalyst's activity.

Prediction of Covalent Metal-Metal Bonding in Cp-M-M'-Nacnac Complexes of Group 2 and 12 Metals (Be, Mg, Ca, Zn, Cd, Hg).

Bimetallic CpMM'Nacnac molecules with group 2 and 12 metals (M=Be, Mg, Ca, Zn, Cd, Hg) that contain novel metal-metal bonding have been investigated in a theoretical study of their molecular and electronic structure, thermodynamic stability, and metal-metal bonding. In all cases the metal-metal bonds are characterized as electron-sharing covalent single bonds from natural bond orbital (NBO) and energy-decomposition analysis with natural orbitals of chemical valence (EDA-NOCV) analysis. The sum of [MM'] charges is relatively constant, with all complexes exhibiting a [MM']2+ core. Quantum theory of atoms in molecules (QTAIM) analysis indicates the presence of non-nuclear attractors (NNA) in the metal-metal bonds of the BeBe, MgMg, and CaCa complexes. There is substantial electron density (0.75-1.33 e) associated with the NNAs, which indicates that these metal-metal bonds, while classified as covalent electron-sharing bonds, retain significant metallic character that can be associated with reducing reactivity of the complex. The predicted stability of these complexes, combined with their novel covalent metal-metal bonding and potential as reducing agents, make them appealing targets for the synthesis of new metal-metal bonds.

A comprehensive analysis of the spectroscopic and drug-like properties of dimethyl 5-hydroxybenzene-1,3-dicarboxylate: insights from DFT and MD simulations

ABSTRACT An in-depth investigation was undertaken to explore the ground-state molecular and electronic structure of Dimethyl 5-hydroxybenzene-1,3-dicarboxylate (D5HD). The potential energy scan of D5HD at various dihedral angles identified the lowest energy conformer, which was re-optimised at the higher basis set. The recorded FT-Raman and FT-IR spectra were compared with corresponding theoretical spectra, demonstrating a high degree of coherence. A comprehensive Natural Bond Orbital analysis was done for D5HD to gain insights for the overall electronic distribution within the molecule. Additionally, nonlinear optical (NLO) properties along with various DFT-based descriptors were calculated. Molecular docking studies against human acetylcholinesterase enzyme (PDB 2JF0), D5HD displayed good binding affinity. The study revealed that hydrogen bonded as well as van der Waals interactions exist between D5HD and different residues in the binding cavity of the target protein. The 100 ns MD simulation confirms the conformational stability of protein (2JF0)–ligand (D5HD) complex. The binding energy and residue decomposition analysis performed using MM-GBSA approach showed positive results with binding energy of −17.76 kcal/mol. Furthermore, D5HD effectively passed all pharmacokinetic filters, establishing its potential as a promising candidate in the quest for novel acetylcholinesterase inhibitors. Highlights Systematic exploration of the structure and spectroscopic properties of D5HD. Molecular docking between D5HD and acetylcholinesterase (PDB 2JF0) reveals the establishment of hydrogen bonds and Van der Waals interactions. D5HD successfully complies with all relevant pharmacokinetic filters. MD simulation analysis on 2JF0–D5HD complex confirms the conformational stability. Energy and residue decomposition based on MMGBSA approach.

Theoretical Study of Metal-Ligand Interactions in Lead Complexes with Radiopharmaceutical Interest.

The 203Pb and 212Pb lead radioisotopes are attracting growing interest as they can aid in the development of personalized, targeted radionuclide treatment for advanced and currently untreatable cancers. In the present study, the bonding interactions of Pb2+ with twelve macrocyclic ligands, having an octa and nona coordination, were assessed using Density Functional Theory (DFT) calculations. The molecular structures in an aqueous solution were computed utilizing the polarized continuum model. The preference for the twisted square antiprismatic (TSAP) structure was confirmed for ten out of the eleven cyclen-based complexes. The characteristics of the bonding were assessed using a Natural Energy Decomposition Analysis (NEDA). The analysis revealed a strong electrostatic character of the bonding in the complexes, with minor variations in electrical terms. The charge transfer (CT) had a comparable energetic contribution only in the case of neutral ligands, while in general, it showed notable variations regarding the various donor groups. Our data confirmed the general superiority of the carboxylate O and aromatic N donors. The combination of the selected efficient pendant arms pointed out the superiority of the acetate pendant arms and the lack of significant cooperation between the different pendant arms in the probed ligands. Altogether, the combination led only to a marginal enhancement in the total CTs in the complexes.

"Hydridic Hydrogen-Bond Donors" Are Not Hydrogen-Bond Donors.

Herein, we dismiss a recent proposal by Civiš, Hobza, and co-workers to modify the IUPAC definition of hydrogen bonds in order to expand the scope from protonic Y-Hδ+ to hydridic Y-Hδ- hydrogen-bond donor fragments [J. Am. Chem. Soc. 2023, 145, 8550]. Based on accurate Kohn-Sham molecular orbital (KS-MO) analyses, we falsify the conclusion that interactions involving protonic and hydridic hydrogens are both hydrogen bonds; they are not. Instead, our quantitative KS-MO, energy decomposition, and Voronoi deformation density analyses reveal two fundamentally different bonding mechanisms for protonic Y-Hδ+ and hydridic Y-Hδ- fragments which go with charge transfer in opposite directions. On one hand, we confirm the IUPAC definition for regular hydrogen bonds in the case of protonic Y-Hδ+ fragments. On the other hand, complexes involving Y-Hδ- fragments are, in fact, acceptors in other well-known families of Lewis-acid/base interactions, such as halogen bonds, chalcogen bonds, and pnictogen bonds. These mechanisms lead to the same spectroscopic phenomenon in both the Y-Hδ+ and Y-Hδ- fragments, that is, the redshift in the Y-H stretching frequency, which is, thus, not an exclusive indicator for hydrogen bonding.

Self-assembled belt[12]pyridine nanotube for the adsorption and removal of anti-TB drugs from wastewater

Environmental safety is a major concern of today's world, and it has attracted the attention of the scientific research community. The environment gets contaminated due to various toxic drugs and pharmaceutical toxins coming from different sources. Therefore, removing these toxins is necessary for making the environment greener and safer. Nanomaterials are excellent candidates for adsorption and removal of harmful materials due to their unique properties. Herein we have studied self-assembled nanotube based on belt[12]pyridine (2-(12BPy)) as an adsorbent for the removal of antituberculosis drugs i.e. isoniazid (INZ) and pyrazinamide (PZA). The adsorption ability is studied by considering the adsorption energy (Eads), non-covalent interaction index analysis (NCI), quantum theory of atoms in molecules (QTAIM), energy decomposition analysis (EDA), charge transfer via natural bond orbital analysis (NBO) and electron localization function (ELF), and dipole moment. The results indicate strong adsorption of selected drugs with self-assembled nanotube without altering the electronic properties. Recovery time of the complexes is also calculated to understand the regeneration of the nanotube. Results also demonstrate that van der Waals forces are major contributor in the adsorption of drug molecules within the nanotube. Our results indicate that belt[12]pyridine (2-(12BPy)) nanotube can be used as an adsorbent for the removal of toxic drug molecules from the environment. This research will benefit the scientists working in the field of new sensors for biological components and materials.

Real-space energy decomposition analysis method for qualitative and quantitative interpretations.

In the work, a real-space energy decomposition analysis method, called DM-EDA(RS), is introduced based on our recently developed DM-EDA method [Zhang et al., J. Chem. Phys. 160, 174101 (2024)]. The EDA terms in DM-EDA(RS), including electrostatic, exchange, repulsion, polarization, and correlation, are expressed as the summations of grid-based energy density in real-space. This method is able to interpret intermolecular interactions in a unified qualitative and quantitative way. DM-EDA(RS) results provide not only comprehensive explanations for intermolecular interactions but also insights for sub-region interactions involving different functional groups.

Exploring the Impact of Protein Chain Selection in Binding Energy Calculations with DFT.

Calculation of binding free energies between a protein and a ligand are highly desired for computer-aided drug design. Here we approximate the binding energies of ABL1, an enzyme which is the target for drugs used in the treatment of chronic myeloid leukaemia, with minimal models and density functional theory (DFT). Starting from the crystal structures of protein-drug complexes, we estimated the binding free energies having used all available individual molecules (protein chains) within each structure, not only a single one as commonly used, in order to see if the choice of the protein chain is important in such calculations. Differences were observed between chains in the same file. Energy decomposition analysis (EDA) revealed that the most important factors for binding were exchange, repulsion and electrostatics. The desolvation term varied dramatically between the inhibitors (between 4.2 and 92.3 kcal/mol). All functionals showed similar patterns in the EDA and in discriminating between the ligands. Non-covalent interactions (NCI) analysis was used to further explain the differences between protein chains and functionals. Overall, it is shown that small minimal models of a drug binding site can be useful to infer on the suitability of an initial crystal structure for further analysis such as EDA.

Open-cage metallo-azafullerenes as efficient single-atom catalysts toward oxygen reduction reaction.

Very recently, open-cage metallo-azafullerenes PbC100N4H4 and Pb2C100N4H4 containing one Pb-N4-C moiety have been synthesized via the electron beam. Herein, we utilized density functional theory calculations in combination with abinitio molecular dynamics (AIMD) simulations to study the geometric and electronic structures, bonding properties, thermodynamic stability, and catalytic performance of MC100N4H4 and M2C100N4H4 (M = Ge, Sn, Pb). Metal-nitrogen distances and metal-metal distances increase along with the metal radius while the metal atom is positively charged. Energy decomposition analysis revealed that the bonding interactions between M and the C100N4H4 fragment could be described as the donor-acceptor interaction between M(ns0(n-1)d10np4) and C100N4H4 fragment, in which the orbital interactions terms contribute more than the electrostatic interactions. AIMD simulations demonstrate that those metallo-azafullerenes exhibit thermodynamic stability at room temperature. These metallo-azafullerenes, which could serve as typical carbon-supported single-atom catalysts, possess enhanced catalytic performance toward the oxygen reduction reaction (ORR) compared to the planar catalysts, which is attributed to the curvature of metallo-azafullerenes. GeC100N4H4 and SnC100N4H4 exhibit high catalytic performance in the 4e-ORR pathway to H2O, whereas only PbC100N4H4 is suitable for the 2e-ORR reaction pathway because of the difficulty in obtaining electrons. All M2C100N4H4 favors the 4e-reaction pathway due to the presence of the axial metal atom. Our finding of open-cage metallo-azafullerenes as efficient single-atom catalysts holds profound implications for both fundamental research in catalysis and practical applications in fuel cells and other electrochemical devices.

Regioselectivity and physical nature of the interactions between (methyl)guanine with HCl and CH3OH

A comprehensive study of the hydrogen bonding interactions between guanine (G) and methyl guanine derivatives (MGs) in the presence of HCl and MeOH is carried out at B3LYP, B3LYP-D3 and M062X/6–311 + + G(d.p) levels using molecular electrostatic potential, natural bond orbital, and symmetry adapted perturbation theory. Making use of these state-of-the-art techniques, this study attempts to elucidate the chemical bonding, regioselectivity, and physical nature of the interactions responsible for the stability of (M)G…L (L = HCl, MeOH) acid–base complexes. Our calculations reveal that 1-G, 3-MG, and 5-MG interact more strongly with MeOH than HCl due to the positive hydrogen bond cooperativity. Furthermore, the carbonyl site on G is found to be the most reactive site, while methyl substitution increases the basicity of the nucleobase, thus yieding more stable complexes. The strongest H-bond interaction in G-complexes is found when HCl and MeOH attack the carbonyl group in anti-position. Finally, energy decomposition analyses through the symmetry-adapted perturbation theory reveal that most complexes are mainly stabilized via electrostatic interactions. The energy difference between complex isomers shows a competition between 3-HCl-G (MG) and 4-HCl-G (MG) at ∆G level where thermal, BSSE and entropy terms are included.

Nature of metal–bis(β-diketonate) bonds in TM(acac’)2 (acac’ = acetylacetonate, hexafluoroacetylacetonate, hexamethylacetylacetonate; TM = Ni(II), Pd(II), Pt(II)) complexes: A density functional theory investigation

This computational study investigates the structural and bonding aspects of transition metal bis(β-diketonate) complexes, TM(acac’)2 (acac’ = acetylacetonate, hexafluoroacetylacetonate, hexamethylacetylacetonate; TM = Ni(II), Pd(II), Pt(II)), at BP86/def2-TZVP and M06/def2-TZVP levels of theory. Patterns observed in both computed interaction and stabilisation energies demonstrate an increasing in the order of hexafluoroacetylacetonate < hexamethylacetylacetonate < acetylacetonate. A comprehensive energy decomposition analysis (EDA) was performed to thoroughly examine the attributes of the bonds between the metal and ligands within complexes. According to the EDA results, it is evident that electrostatic interactions exert a greater influence on shaping the overall attractive interactions compared to orbital interactions. However, the significance of the b1g orbitals cannot be overstated, as they contribute more than 41% across all complexes. Additionally, the backdonation processes TM→Lσ and TM→Lπ exhibit notable weaknesses.

Exploring protein-berberine interactions via molecular dynamics and MM/GBSA calculations

Structural modification on berberine has garnered considerable attention among medicinal chemists as a means to enhance its therapeutic potential. The dissection of macrobiomolecule-berberine interactions is pivotal for guiding such modifications. In this study, four diverse crystal complexes of protein-berberine were selected, and simulated by molecular dynamics, in which each complex was conducted five times. Subsequently, the interactions of protein-berberine were collected, and analyzed by in-house scripts. Our findings indicated that the four oxygen atoms, namely O15, O16, O17, and O18, actively participated in the formation of hydrogen bonds and water bridges; the carbon atoms within the D and E rings were able to form hydrophobic interactions more than other carbon atoms; the positively charged N7 atom predominantly participated in the formation of cation-π interactions; π-π interactions were prevalent in the A, B, and D rings, with their frequency influenced by the presence of the residues Trp, Tyr, His, and Phe in the binding pocket. In the free energy decomposition analysis, the terms ΔG_Lipo and ΔG_vdW emerged as the most significant contributors to the overall energy, with ΔG_vdW accounting for over 50% of the energy contribution. Notably, the contribution of ΔG_Packing to the free energy was substantially greater than that of ΔG_Hbond, which contributed less than 1%. These insights into molecular interactions provide valuable guidance for the structural modification of berberine, such as site selection of structural modification, group optimization with the improvement of van der Waals and hydrophobic interactions, and the trade-off between hydrogen bonding and π-packing interactions in terms of interaction weights.

β-d-Glucopyranose – Silver+ (1:1) complex assisted aromaticity involving electron deficient X3 (X=Be, Mg) core

Seeing the frontier molecular orbitals (FMO) picture of β-d-glucopyranose–silver+ (1:1) complex ([Ag(C6H12O6)]+) reported in a recent article (Mondal, 2024) [53], we used the FMO based approach to obtain inverted sandwich like organometallic complex comprising X3 (X=Be, Mg) core, in silico. Thermodynamic feasibility of obtaining [(C6H12O6)Ag-X3-Ag(C6H12O6)] (X=Be, Mg) complexes are explored. Conceptual density functional theory (CDFT) based reactivity descriptors guide to ascertain the stability order of such complexes. By the help of natural bond orbital analysis, Wiberg bond indices (WBI), bond critical point analysis, electron density descriptors, and gradient isosurface between selected contacts, the nature of interaction between the X3 core and Ag metals are studied. FMO of the complexes further help in ascertaining the nature of binding. Energy decomposition analysis marks the interaction between X3 and Ag metals in [(C6H12O6)Ag-X3-Ag(C6H12O6)] (X=Be, Mg) as partially covalent. Negative values of nucleus independent chemical shift (NICS) at the triangular ring (X3) centres [NICS(0)] indicate aromatic property.

Trends in the Reactivity of Pentacyclic Ether Derivatives on Silicon and Germanium Surfaces Revealed by Energy Decomposition Analysis for Extended Systems

We analyze the adsorption and ring‐opening reaction of pentacyclic chalcogen alkyls with chalcogen atoms ranging from oxygen to tellurium using computational analysis. Thus extends our previous investigation of THF on silicon towards germanium surface and towards the heavier chalcogen homologues. We found an increasing dative bond strength of the precursor state with the period of the chalcogen atom on both surfaces. Using energy decomposition analysis for extended systems (pEDA), differences in the trends between silicon and germanium were revealed. Following the dative bound state, subsequent ring‐opening reactions were found to proceed via a nucleophilic back‐side attack, similar to a molecular SN2 reaction. Reaction energies and barrier heights show a maximum for the sulfur derivative. This trend was found to correlate with the position of the transition state in the reaction and agrees with the molecular reactivity. The uniquely low barrier of THF can be attributed to its ring strain. Our findings shed light on the fundamental aspects of surface chemistry and the benefits of using quantitative bonding analysis in this field. It also paves the way for future explorations into the functionalization of semiconductor surfaces using higher‐period ether derivatives.

4+2]‐Cycloaddition Reactions to Corannulene Accelerated by η6‐Coordination of Ruthenium Complexes

The influence of different (LRu)nII (L = Cp– (n = 1, 2); C6H6) moieties η6‐coordinated to corannulene (CpRu+, (Benzene)Ru2+, and (CpRu)22+ has been explored by using the Activation Strain Model (ASM) of reactivity together with the Energy Decomposition Analysis (EDA) scheme. The results obtained have been compared with those obtained for the non‐coordinated corannulene counterpart. It has been observed that the presence of the ruthenium moiety favors, both kinetically and thermodynamically, the [4+2]‐cycloaddition reaction with cyclopentadiene. The factors governing these differences depend on the nature of the [Ru] complex. Whereas the interaction energy is solely responsible for the lower activation energy for the monocationic system (CpRu+), not only the higher interaction but also the lower strain energy explain the much higher reactivity observed for the dicationic systems ((C6H6)Ru2+ and (CpRu)22+). In this latter case, the Pauli repulsion energy term is the factor controlling the differences in the interaction energies, following the so‐called Pauli repulsion‐lowering concept. Finally, when η6‐coordinating mono‐ or dicationic ruthenium complexes, although although the [4+2]‐cycloaddition reaction in the external rim bonds remains thermodynamically and kinetically more favourable, the cycloaddition becomes also thermodynamically feasible in some internal bonds.

Pharmacophore-Assisted Covalent Docking Identifies a Potential Covalent Inhibitor for Drug-Resistant Genotype 3 Variants of Hepatitis C Viral NS3/4A Serine Protease.

The emergence of drug-resistance-inducing mutations in Hepatitis C virus (HCV) coupled with genotypic heterogeneity has made targeting NS3/4A serine protease difficult. In this work, we investigated the mutagenic variations in the binding pocket of Genotype 3 (G3) HCV NS3/4A and evaluated ligands for efficacious inhibition. We report mutations at 14 positions within the ligand-binding residues of HCV NS3/4A, including H57R and S139P within the catalytic triad. We then modelled each mutational variant for pharmacophore-based virtual screening (PBVS) followed by covalent docking towards identifying a potential covalent inhibitor, i.e., cpd-217. The binding stability of cpd-217 was then supported by molecular dynamic simulation followed by MM/GBSA binding free energy calculation. The free energy decomposition analysis indicated that the resistant mutants alter the HCV NS3/4A-ligand interaction, resulting in unbalanced energy distribution within the binding site, leading to drug resistance. Cpd-217 was identified as interacting with all NS3/4A G3 variants with significant covalent docking scores. In conclusion, cpd-217 emerges as a potential inhibitor of HCV NS3/4A G3 variants that warrants further in vitro and in vivo studies. This study provides a theoretical foundation for drug design and development targeting HCV G3 NS3/4A.

The Chemical Bond at the Bottom of the Periodic Table: The Case of the Heavy Astatine and the Super-Heavy Tennessine.

In this work, we study the chemical bond in molecules containing heavy and super-heavy elements according to the current state-of-the-art bonding models. An Energy Decomposition Analysis in combination with Natural Orbital for Chemical Valence (EDA-NOCV) within the relativistic four-component Dirac-Kohn-Sham (DKS) framework is employed, which allows to successfully include the spin-orbit coupling (SOC) effects on the chemical bond description. Simple halogen-bonded adducts ClX⋯L (X=At, Ts; L=NH3, Br-, H2O, CO) of astatine and tennessine have been selected to assess a trend on descending along a group, while modulating the ClX⋯L bond features through the different electronic nature of the ligand L. Interesting effects caused by SOC have been revealed: i) a huge increase of the ClTs dipole moment (which is almost twice as that of ClAt), ii) a lowering of the ClX⋯L bonding energy arising from different contributions to the ClX…L interaction energy strongly depending on the nature of L, iii) a quenching of one of the π back-donation components to the bond. In the ClTs(CO) adduct, the back-donation from ClTs to CO becomes the most important component. The analysis of the electronic structure of the ClX dimers allows for a clear interpretation of the SOC effects in these systems.