Theoretical Investigation of the Binding of Nucleobases to Cucurbiturils by Dispersion Corrected DFT Approaches

Driving forces and large scale affinity calculations for piperine/γ-cyclodxetrin complexes: Mechanistic insights from umbrella sampling simulation and DFT calculations

Piperine (PiP), a bioactive molecule, exhibits numerous health benefits and is frequently employed as a co-delivery agent with various phytomedicines (e.g., curcumin) to enhance their bioavailability. This is attributed to PiP's inhibitory activity against drug-metabolizing proteins, notably CYP3A4. Nevertheless, PiP encounters solubility challenges addressed in this study using cyclodextrins (CDs). Specifically, γ-CD and its derivatives, Hydroxypropyl-γ-CD (HP-γ-CD), and Octakis (6-O-sulfo)-γ-CD (Octakis-S-γ-CD), were employed to form supramolecular complexes with PiP. The conformational space of the complexes was assessed through 1 μs molecular dynamics simulations and umbrella sampling. Additionally, quantum mechanical calculations using wB97X-D dispersion-corrected DFT functional and 6-311 + G(d,p) basis set were conducted on the complexes to examine the thermodynamics and kinetic stability. Results indicated that Octakis-S-γ-CD exhibits superior host capabilities for PiP, with the most favorable complexation energy (-457.05 kJ/mol), followed by HP-γ-CD (-249.16 kJ/mol). Furthermore, two conformations of the Octakis-S-γ-CD/PiP complex were explored to elucidate the optimal binding orientation of PiP within the binding pocket of Octakis-S-γ-CD. Supramolecular chemistry relies significantly on non-covalent interactions. Therefore, our investigation extensively explores the critical atoms involved in these interactions, elucidating the influence of substituted groups on the stability of inclusion complexes. This comprehensive analysis contributes to emphasizing the γ-CD derivatives with improved host capacity.

Host-Guest Interactions of Ruthenium(II) Arene Complexes with Cucurbit[7/8]uril

Cucurbit[n]urils (CB[n]s) have been recognized for their chemical and thermal stability, and their ability to bind many neutral and cationic guest molecules makes them excellent hosts in a range of supramolecular applications. In drug delivery, CB[n]s can enhance drug solubility, improve chemical and physical drug stability, and allow for triggered and controlled release. This study aimed to investigate the ability of CB[7] and CB[8] as molecular hosts to bind ruthenium(II) arene complexes that are current anticancer lead structures in the area of metallodrugs. Both, experimental and computational methods, led to insights into the binding preferences and geometries of [RuII(cym)Cl2]2 (1; cym = η6-p-cymene), [RuII(cym)(dmb)Cl2]) (2; cym = η6-p-cymene; dmb = 1,3-dimethylbenzimidazol-2-ylidene), and [RuII(cym)(pta)Cl2] (3, RAPTA-C; cym = η6-p-cymene; pta = 1,3,5-triaza-7-phospha-adamantane) with CB[7] and CB[8]. Competition experiments by mass spectrometry revealed clear preferences of 2 for CB[8] and 3 for CB[7]. Based on a comparison of the associated interaction energies from quantum chemical calculations as well as experimental data, 3@CB[7] clearly prefers a binding mode, where the pta ligand is located inside the cavity of the host, and the metal ion interacts with two of the carbonyl groups on the rim of CB[7]. In contrast, complex 2 binds in two different orientations with interaction energies similar to those of both CB[n]s, with the cym ligand being either inside or outside of the cavity. These findings suggest that ruthenium(II) arene complexes are able to form stable host-guest interactions with CB[n]s, which can be exploited as drug delivery vehicles in further metallodrug development to improve their chemical stability.

Dispersion-corrected DFT calculations and umbrella sampling simulations to investigate stability of Chrysin-cyclodextrin inclusion complexes

The study of inclusion complexes of Chrysin (ChR) with three forms of cyclodextrins (CDs) α-, β-, and γ-CD was accomplished to examine the stability of ChR inside the central cavities of CDs. The aim of study was to identify the most suitable form of CD to improve the hydro-solubility of poorly soluble ChR bioactive molecule. Microsecond timescale molecular dynamics (MD) simulations were performed on four inclusion complexes (α-CD/ChR, β-CD/ChR, and two conformations of γ-CD/ChR) to examine the dynamics of ChR inside the cavity of CDs. The first conformation of γ-CD/ChR inclusion complex (γ-CD1/ChR) was identified to possess the highest affinity between host and guest molecule on the basis of binding energy calculated by employing Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) and umbrella sampling simulations. To further strengthen the claims of classical and biased MD studies, Our own N-layered Integrated molecular Orbital and Molecular mechanics (ONIOM) (wB97XD/6-311+g(d,p):pm7) calculations were performed on the selected inclusion complexes. The ONIOM based complexation energy reaffirmed that ChR had highest affinity for the γ-CD1 host molecule. Further, the non-covalent interaction analysis was conducted using Multiwfn software on QM-optimized inclusion complexes with wB97XD/6-311+G(d,p) model chemistry, revealing non-covalent interactions between ChR and CDs. This atomic level information helped us to gain better insights into critical atoms of ChR and CD that participated in intermolecular interactions and identify γ-CD as a suitable host molecule for improving the hydro-solubulity of ChR. The structural insights would help to derive new derivatives of γ-CD with better host capacity.

A perception into binary and ternary copper (II) complexes: synthesis, characterization, DFT modeling, antimicrobial activity, protein binding screen, and amino acid interaction

Ensuring healthy lives and promoting well-being for all at all ages is the third goal of the sustainable development plan, so it was necessary to identify the most important problems that threaten health in our world. The World Health Organization declared that antibiotic resistance is one of the uppermost global public health threats facing humanity and searching for new antibiotics is slow. This problem can be approached by improving available drugs to combat various bacterial threats. To circumvent bacterial resistance, three copper(II) complexes based on the pefloxacin drug were prepared and characterized using analytical, spectroscopic, and thermal techniques. The resulting data suggested the formation of one octahedral binary and two distorted square pyramidal ternary complexes. Fluorescence spectra results revealed the formation of a turn-on fluorophore for amino acid detection. Computational calculations investigated quantum and reactivity parameters. Molecular electrostatic potential profiles and noncovalent bond interaction-reduced density gradient analysis indicated the active sites on the complex surface. The complexes were subjected to six microbial species, where the octahedral binary complex provoked its antimicrobial potency in comparison with ternary complexes. The enhanced antimicrobial activity against gram-negative bacterium E-coli compared to gentamicin was exhibited by the three complexes. Docking simulation was performed based on the crystal structure of E. coli and S. pneumoniae receptors using 5I2D and 6O15 codes. The binary complex exhibited a potent fitness score with 5I2D (TBE = - 107 kcal/mol) while ternary complexes displayed the highest docked score of fitness with 6O15.

In acid-aminopyrimidine continuum: experimental and computational studies of furan tetracarboxylate-2-aminopyrimidinium salt

Salts and cocrystals are the two important solid forms when a carboxylic acid crystallizes with an aminopyrimidine base such that the extent of proton transfer distinguishes between them. The ΔpK_a value (pK_a(base) − pK_a(acid)) predicts whether the proton transfer will occur or not. However, the ΔpK_a range, 0 < ΔpK_a < 3, is elusive where the formation of cocrystal or salt cannot be predicted. The current study has been done to obtain a generalization in this elusive range with the Cambridge Structural Database (CSD). Based on the generalization, a novel salt (FTCA)⁻(2-AP)⁺ of furantetracarboxylic acid (FTCA) with 2-aminopyrimidine (2-AP) is obtained. The structural confirmation was done by single-crystal X-ray diffraction (SCXRD). Density functional theory (DFT) calculations were performed at the IEF-PCM-B3LYP-D3/6-311G(d,p) level to optimize the geometrical coordinates of salt for frontier molecular orbitals (FMOs) and molecular electrostatic potential (MESP). The geometrical parameters of most of the atoms of the optimized salt structure were comparable with SCXRD data. Additionally, results of other computational methods such as ab initio (Hartree–Fock; HF and second-order-Møller–Plesset perturbation; MP2) and semi-empirical were also compared with experimental results of the salt. Quantum theory of atoms in molecules (QTAIM), reduced density gradient (RDG), and natural bond orbital (NBO) analyses were done to calculate the strength and nature of non-covalent interactions present in the salt. Furthermore, Hirshfeld surface analysis, interaction energy calculations, and total energy frameworks were performed for qualitative and quantitative estimations of strong and weak intermolecular interactions.

Generalization in the elusive ΔpK_a range, experimental and computational studies of furan tetracarboxylate-2-aminopyrimidinium salt.

Host-Guest Complexation of Cucurbit[7]Uril and Cucurbit[8]Uril with the Antineoplastic and Multiple Sclerosis Agent Mitoxantrone (Novantrone)

The nature of interactions between the neutral/protonated mitoxantrone and the cucurbit[n]uril (n = 7, 8) host system was analyzed by employing density functional theory calculations. A comparison between the inclusion complexes of CB[7] and CB[8] shows various subtle differences in the complexation thermodynamics, given as changes in the Gibbs energy. Doubly and quadruply charged mitoxantrone (MX) molecules spontaneously form complexes in a water solvent, which are modeled using the polarizable continuum model approach. Both CB[7] and CB[8] complexes are stable as the geometry of the cavity allows for electrostatic interactions between the charged MX arms and the rim of the CB cavity. CB[8] also forms a stable complex with two mitoxantrone molecules with their aromatic rings stacked inside the cavity. Both CB[7] and CB[8] show properties that can be utilized in drug delivery.

Theoretical Insight into the Interaction between Chloramphenicol and Functional Monomer (Methacrylic Acid) in Molecularly Imprinted Polymers

Molecular imprinting technology is a promising method for detecting chloramphenicol (CAP), a broad-spectrum antibiotic with potential toxicity to humans, in animal-derived foods. This work aimed to investigate the interactions between the CAP as a template and functional monomers required for synthesizing efficient molecularly imprinted polymers for recognition and isolation of CAP based on density functional theory. The most suitable monomer, methacrylic acid (MAA), was determined based on interaction energies and Gibbs free energy changes. Further, the reaction sites of CAP and MAA was predicted through the frontier molecular orbitals and molecular electrostatic potentials. Atoms in molecules topology analysis and non-covalent interactions reduced density gradient were applied to investigate different types of non-covalent and inter-atomic interactions. The simulation results showed that CAP was the main electron donor, while MAA was the main electron acceptor. Moreover, the CAP–MAA complex simultaneously involved N-H···O and C=O···H double hydrogen bonds, where the strength of the latter was greater than that of the former. The existence of hydrogen bonds was also confirmed by theoretical and experimental hydrogen nuclear magnetic resonance and Fourier transform infrared spectroscopic analyses. This research can act as an important reference for intermolecular interactions and provide strong theoretical guidance regarding CAP in the synthesis of molecularly imprinted polymers.

Encapsulation of the Sulfur Compounds by Cucurbit[7]uril: A Quantum Chemistry Study

Benzothiophene (BT) and dibenzothiophene (DT) are the most important contaminants in the petroleum derivatives responsible for serious environmental and health problems. Therefore, we have investigated the absorption of these compounds for the first time by considering cucurbit[7]uril (CB[7]) as the host molecule and using the theoretical levels of density functional theory//B3LYP-D3/6-31G(d). BT and DT absorbed into CB[7] do not undergo a significant structural change in the CB[7] structure. The energy gap of the S-compounds@CB[7] in water and hexane solvents was approximately 5 eV, and this large value implies that the complexes have high chemical stability. Moreover, the absorption of the BT and DT into CB[7] in the water and hexane solvents is a favorable process, whereas the lowest binding energy was observed between the dibenzothiophene and CB[7] in the DT@CB[7] complex. The solvation enthalpy shows a preferential solvation of the complexes in water than in hexane solvent. This trend is confirmed by the AIM analysis that shows the highest stability for the DT@CB[7] system with the contribution of cooperative hydrogen bonding. The transfer free energy of S-compounds@CB[7] complexes from hexane to water are -66.12 and -59.56 kcal/mol for BT@CB[7] and DT@CB[7], respectively, implying the spontaneous transference of these complexes from hexane to water solvent. Overall, our results show that the cucurbiturils can be a new class of host molecules to be used in the removal of S-compounds from petroleum derivatives. Finally, a schematic flow diagram of the desulfurization process by cucurbiturils was proposed.

Controlling anticancer drug mediated G-quadruplex formation and stabilization by a molecular container

Controlling of ligand mediated G-quadruplex DNA (GQ-DNA) formation and stabilization is an important and challenging aspect due to its active involvement in many biologically important processes such as DNA replication, transcription, etc. Here, we have demonstrated that topotecan (TPT), a potential anticancer drug, can instigate the formation and stabilization of GQ-DNA (H24 → GQ-DNA) in the absence of Na⁺/K⁺ ions via circular dichroism, fluorescence, NMR, UV melting and molecular dynamics (MD) simulation studies. The primary binding mode of TPT to GQ was found to be stacking at the terminal rather than binding to the groove. We have also reverted this conformational transition (GQ-DNA → H24) using a molecular container, cucurbit[7]uril (CB7), by means of the translocation of the drug (TPT) from GQ-DNA to its nanocavity. Importantly, we have carried out the detection of these conformational transitions using the fluorescence color switch of the drug, which is more direct and simple than some of the other methods that involve sophisticated and complex detection techniques.

Are cucurbiturils better drug carriers for bent metallocenes? Insights from theory

Bent metallocenes (BM) have anti-tumor properties but they face a serious drug efficacy problem due to poor aqueous solubility and rapid hydrolysis under physiological conditions. These two problems can be fixed by encapsulating them in host molecules such as cyclodextrin (CD), cucurbituril (CB) etc. Experimentally, CD-BM, CB-BM host-guest complexes have been investigated to check the efficiency of the drug delivery and efficiency of the encapsulated drug. CB has been reported to be a better host than CD but the reasons for this has not been figured out. This can be done by finding out the mechanism of binding and the nature of the binding forces in both the inclusion complexes. This is exactly done here by performing a DFT study at BP86/TZP level on CB-BM host-guest systems. For comparison CD-BM with β-cyclodextrin as host have been studied. Four BMs (Cp₂MCl₂, M=Ti, V, Nb, Mo) and their corresponding cations (Cp₂MCl⁺, Cp₂M²⁺) are chosen as guests and they are encapsulated into cucurbit-[6]-uril (CB[6]) and cucurbit-[7]-uril(CB[7]) host systems. Computations reveal that CB[7] accommodates well the BMs over CB[6] due to their larger cavity size and also CB[7] is found to be a better host than β-cyclodextrin. BMs enter vertically rather than horizontally into the CB cavity. The reversible binding of BMs within CB[7] is controlled by various non-bonding interactions and mainly by hydrogen bonding between the portal oxygen atoms and Cp protons as revealed by QTAIM analysis. On the other hand, the interaction between the wall nitrogen atoms in CB[7] and chlorine atoms attached to the metal in BM strengthens the M-Cl bonds that prevents rapid hydrolysis of M-Cl and M-Cp bonds saving the drug. Comparatively, BMs experience less electrostatic attraction and more Pauli repulsion within β-cyclodextrin cavity and this affects the drug binding with CD. This makes β-cyclodextrin a less suitable drug carrier for BMs than CBs. Among the four BMs, niobocene binds strongly and titanocene binds weakly with CBs. EDA clearly shows that all the interactions between the guest and host are non-covalent in nature and electrostatic interactions outperform high-repulsion resulting in stable complexes. Cations form stronger complexes than neutral BMs. FMO analysis reveals that neutral BMs are less reactive compared to their cations and complexes are more reactive in CB[6] environment due to excess strain. QTAIM analysis helps to bring out the newer insights in these types of host-guest systems.

Nature of bonding and cooperativity in linear DMSO clusters: A DFT, AIM and NCI analysis

This study aims to cast light on the nature of interactions and cooperativity that exists in linear dimethyl sulfoxide (DMSO) clusters using dispersion corrected density functional theory. In the linear DMSO, DMSO molecules in the middle of the clusters are bound strongly than at the terminal. The plot of the total binding energy of the clusters vs the cluster size and mean polarizabilities vs cluster size shows an excellent linearity demonstrating the presence of cooperativity effect. The computed incremental binding energy of the clusters remains nearly constant, implying that DMSO addition at the terminal site can happen to form an infinite chain. In the linear clusters, two σ-hole at the terminal DMSO molecules were found and the value on it was found to increase with the increase in cluster size. The quantum theory of atoms in molecules topography shows the existence of hydrogen and SO⋯S type in linear tetramer and larger clusters. In the dimer and trimer SO⋯OS type of interaction exists. In 2D non-covalent interactions plot, additional peaks in the regions which contribute to the stabilization of the clusters were observed and it splits in the trimer and intensifies in the larger clusters. In the trimer and larger clusters in addition to the blue patches due to hydrogen bonds, additional, light blue patches were seen between the hydrogen atom of the methyl groups and the sulphur atom of the nearby DMSO molecule. Thus, in addition to the strong H-bonds, strong electrostatic interactions between the sulphur atom and methyl hydrogens exists in the linear clusters.

Surface modification induced enhanced CO2 sorption in cucurbit[6]uril, an organic porous material

Theoretical and experimental studies on CO₂ sorption of the organic macrocycle cucurbit[6]uril (CB[6]).