Hydrogen bonding patterns and C-H...π interactions in the structure of the antiparkinsonian drug (R)-rasagiline mesylate determined using laboratory and synchrotron X-ray powder diffraction data.

The structure of (R)-rasagiline mesylate [(R)-RasH+·Mes-], an active pharmaceutical ingredient used to treat Parkinson's disease, is presented. The structure was determined from laboratory and synchrotron powder diffraction data, refined using the Rietveld method, and validated and optimized using dispersion-corrected DFT calculations. The unit-cell parameters obtained in both experiments are in good agreement and the refinement with both datasets converged to good agreement factors. The final parameters obtained from laboratory data were a = 5.4905 (8), b = 6.536 (2), c = 38.953 (3) Å, V = 1398.0 (4) Å3 and from synchrotron powder data were a = 5.487530 (10) Å, b = 6.528939 (12) Å, c = 38.94313 (9) Å, V = 1395.245 (5) Å3 with Z = 4 and space group P212121. Preferred orientation was properly accounted for using the synchrotron radiation data, leading to a March-Dollase parameter of 1.140 (1) instead of the 0.642 (1) value obtained from laboratory data. In the structure, (R)-RasH+ moieties form layers parallel to the ab plane connected by mesylate ions through N-H...O and C-H...O hydrogen bonds. These layers stack along the c axis and are further connected by C-H...π interactions. Hirshfeld surface analysis and fingerprint plot calculations indicate that the main interactions are: H...H (50.9%), H...C/C...H (27.1%) and H...O/O...H (21.1%).

Leucopterin, the white pigment in butterfly wings: structural analysis by PDF fit, FIDEL fit, Rietveld refinement, solid-state NMR and DFT-D.

Leucopterin (C6H5N5O3) is the white pigment in the wings of Pieris brassicae butterflies, and other butterflies; it can also be found in wasps and other insects. Its crystal structure and its tautomeric form in the solid state were hitherto unknown. Leucopterin turned out to be a variable hydrate, with 0.5 to about 0.1 molecules of water per leucopterin molecule. Under ambient conditions, the preferred state is the hemihydrate. Initially, all attempts to grow single crystals suitable for X-ray diffraction were to no avail. Attempts to determine the crystal structure by powder diffraction using the direct-space method failed, because the trials did not include the correct, but rare, space group P2/c. Attempts were made to solve the crystal structure by a global fit to the pair distribution function (PDF-Global-Fit), as described by Prill and co-workers [Schlesinger et al. (2021). J. Appl. Cryst. 54, 776-786]. The approach worked well, but the correct structure was not found, because again the correct space group was not included. Finally, tiny single crystals of the hemihydrate could be obtained, which allowed at least the determination of the crystal symmetry and the positions of the C, N and O atoms. The tautomeric state of the hemihydrate was assessed by multinuclear solid-state NMR spectroscopy. 15N CPMAS spectra showed the presence of one NH2 and three NH groups, and one unprotonated N atom, which agreed with the 1H MAS and 13C CPMAS spectra. Independently, the tautomeric state was investigated by lattice-energy minimizations with dispersion-corrected density functional theory (DFT-D) on 17 different possible tautomers, which also included the prediction of the corresponding 1H, 13C and 15N chemical shifts in the solid. All methods showed the presence of the 2-amino-3,5,8-H tautomer. The DFT-D calculations also confirmed the crystal structure. Heating of the hemihydrate results in a slow release of water between 130 and 250 °C, as shown by differential thermal analysis and thermogravimetry (DTA-TG). Temperature-dependent powder X-ray diffraction (PXRD) showed an irreversible continuous shift of the reflections upon heating, which reveals that leucopterin is a variable hydrate. This observation was also confirmed by PXRD of samples obtained under various synthetic and drying conditions. The crystal structure of a sample with about 0.2 molecules of water per leucopterin was solved by a fit with deviating lattice parameters (FIDEL), as described by Habermehl et al. [Acta Cryst. (2022), B78, 195-213]. A local fit, starting from the structure of the hemihydrate, as well as a global fit, starting from random structures, were performed, followed by Rietveld refinements. Despite dehydration, the space group remains P2/c. In both structures (hemihydrate and variable hydrate), the leucopterin molecules are connected by 2-4 hydrogen bonds into chains, which are connected by further hydrogen bonds to neighbouring chains. The molecular packing is very efficient. The density of leucopterin hemihydrate is as high as 1.909 kg dm-3, which is one of the highest densities for organic compounds consisting of C, H, N and O only. The high density might explain the good light-scattering and opacity properties of the wings of Pieris brassicae and other butterflies.

How Accurate Can Crystal Structure Predictions Be for High-Energy Molecular Crystals?

Molecular crystals have shallow potential energy landscapes, with multiple local minima separated by very small differences in total energy. Predicting molecular packing and molecular conformation in the crystal generally requires ab initio methods of high accuracy, especially when polymorphs are involved. We used dispersion-corrected density functional theory (DFT-D) to assess the capabilities of an evolutionary algorithm (EA) for the crystal structure prediction (CSP) of well-known but challenging high-energy molecular crystals (HMX, RDX, CL-20, and FOX-7). While providing the EA with the experimental conformation of the molecule quickly re-discovers the experimental packing, it is more realistic to start instead from a naïve, flat, or neutral initial conformation, which reflects the limited experimental knowledge we generally have in the computational design of molecular crystals. By doing so, and using fully flexible molecules in fully variable unit cells, we show that the experimental structures can be predicted in fewer than 20 generations. Nonetheless, one must be aware that some molecular crystals have naturally hindered evolutions, requiring as many attempts as there are space groups of interest to predict their structures, and some may require the accuracy of all-electron calculations to discriminate between closely ranked structures. To save resources in this computationally demanding process, we showed that a hybrid xTB/DFT-D approach could be considered in a subsequent study to push the limits of CSP beyond 200+ atoms and for cocrystals.

Bracelet-like Complexes of Lithium Fluoride with Aromatic Tetraamides, and Their Potential for LiF-Mediated Self-Assembly: A DFT Study.

Geometries and binding energies of complexes between a LiF molecule and a model aromatic tetraamide are obtained using various DFT methods. The tetraamide consists of a benzene ring and four amides positioned so that the LiF molecule can bind via Li⋯O=C or N-H⋯F interactions. The complex with both interactions is the most stable one, followed by the complex with only N-H⋯F interactions. Doubling the size of the former resulted in a complex with a LiF dimer sandwiched between the model tetraamides. In turn, doubling the size of the latter resulted in a more stable tetramer with bracelet-like geometry having the two LiF molecules also sandwiched but far apart from each other. Additionally, all methods show that the energy barrier to transition to the more stable tetramer is small. The self-assembly of the bracelet-like complex mediated by the interactions of adjacent LiF molecules is demonstrated by all computational methods employed.

Combined DFT-D3 Computational and Experimental Studies on g-C3N4: New Insight into Structure, Optical, and Vibrational Properties.

Graphitic carbon nitride (g-C3N4) has emerged as one of the most promising solar-light-activated polymeric metal-free semiconductor photocatalysts due to its thermal physicochemical stability but also its characteristics of environmentally friendly and sustainable material. Despite the challenging properties of g-C3N4, its photocatalytic performance is still limited by the low surface area, together with the fast charge recombination phenomena. Hence, many efforts have been focused on overcoming these drawbacks by controlling and improving the synthesis methods. With regard to this, many structures including strands of linearly condensed melamine monomers, which are interconnected by hydrogen bonds, or highly condensed systems, have been proposed. Nevertheless, complete and consistent knowledge of the pristine material has not yet been achieved. Thus, to shed light on the nature of polymerised carbon nitride structures, which are obtained from the well-known direct heating of melamine under mild conditions, we combined the results obtained from XRD analysis, SEM and AFM microscopies, and UV-visible and FTIR spectroscopies with the data from the Density Functional Theory method (DFT). An indirect band gap and the vibrational peaks have been calculated without uncertainty, thus highlighting a mixture of highly condensed g-C3N4 domains embedded in a less condensed "melon-like" framework.

Fine Tuning the Intermolecular Interactions of Water Clusters Using the Dispersion-Corrected Density Functional Theory.

Dispersion-inclusive density functional theory (DFT) methods have unequivocally demonstrated improved performances with respect to standard DFT approximations for modeling large and extended molecular systems at the quantum mechanical level. Yet, in some cases, disagreements with highly accurate reference calculations, such as CCSD(T) and quantum Monte Carlo (MC) calculations, still remain. Furthermore, the application of general-purpose corrections, such as the popular Grimme's semi-classical models (DFT-D), to different Kohn-Sham exchange-correlation functionals sometimes leads to variable and inconsistent results, which recommend a careful prior evaluation. In a recent study, we proposed a simple optimization protocol for enhancing the accuracy of these DFT-D methods by following an alternative and system-specific approach. Here, adopting the same computational strategy, we show how the accurate MC intermolecular interactions of a large set of water clusters of variable sizes (i.e., 300 (H2O)n structures, n = 9, 15, 27) can be reproduced remarkably well by dispersion-corrected DFT models (i.e., B3LYP-D4, PBE-D4, revPBE(0)-D4) upon re-optimization, reaching a mean absolute error per monomer of ~0.1 kcal/mol. Hence, the obtained results support the use of this procedure for fine-tuning tailored DFT-D models for the accurate description of targeted molecular systems.

Complementary base pair interactions between different rare tautomers of the second-generation artificial genetic alphabets.

The functionality of a semisynthetic DNA in the biological environment will depend on the base pair nature of its complementary base pairs. To understand this, base pair interactions between complementary bases of recently proposed eight second-generation artificial nucleobases are studied herein by considering their rare tautomeric conformations and a dispersion-corrected density functional theoretic method. It is found that the binding energies of two hydrogen-bonded complementary base pairs are more negative than those of the three hydrogen-bonded base pairs. However, as the former base pairs are endothermic, the semisynthetic duplex DNA would involve the latter base pairs.

Tuning the catalytic activity of Ag7Au6 cluster for oxygen reduction reaction via support interactions.

We investigate the suitability as oxygen reduction reaction (ORR) catalysts of Ag7Au6 clusters supported by nitrogen-doped graphene (N-graphene) sheets with two different nitrogen contents (3.61% and 8.43%), as well as TiO2(101) monolayers, using post-DFT methods. We employ two energetic descriptors (namely, the O2 and OH binding energies and O-O bond length), to compare the effect of two different supports on the catalytic activity of the Ag7Au6 cluster. We find that low adsorption energies for O2 (approximately -1 eV) and OH (below the one of its reaction on Pt), and O-O bond elongation (which reaches up to 1.34 Å in our study, corresponding to the superoxide bond length) cause the high ORR activity of our materials of interest. Our results show that the supported clusters have higher catalytic performance compared to the unsupported ones due to the lower binding energies of the O2 and OH molecules. Moreover, the TiO2(101) monolayer is a better support for the Ag7Au6 cluster, due to the O-O bond distance increasing up to 1.33 Å, facilitating O2 dissociation. In addition, the low (-0.35 eV) OH adsorption energy decreases the difficulty of OH removal. Through free energy calculations, we prove the ORR mechanism in acidic conditions to display an over-potential of -1.76 V, competitive with to the ones featured by currently available catalysts.

Stability, spectroscopic, electrochemistry and QTAIM analysis of Cu-Znn-1On clusters for glucose sensing application: A study on theoretical and experimental insights.

Clusters of (ZnO)n (n = 2-4) have been shown to play a central role in the detection of glucose entity based on the existence of photo-induced electrons (PE), which facilitates the interaction between (ZnO)n clusters and glucose entity guests. The electrochemistry experiment has confirmed the detection of glucose by the title clusters. The optimization, energetic parameters, and vibrational frequency calculations have indicated that the Cu-Znn-1On-glucose are more stable than the (ZnO)n-glucose complexes. It has been demonstrated that the Cu doping enhanced the chemical behavior of the clusters and formed a high intramolecular charge transfer (ICT) in the system. The glucose sensing by all the forms of Cu-Znn-1On clusters showed that the Cu-Zn3O4, Cu-Wurtzite, and Cu-Rocksalt clusters are the most suitable for adsorbing the glucose guest. The HOMO/LUMO iso-surfaces of the complexes showed that the electron concentrations are localized in the d orbitals and mainly in the form of the d10 orbitals around Zn atoms. The molecular electrostatic potential (MEP) has clearly indicated that a high charge transfer occurs between the copper and the oxygen atoms, which facilitate the adsorption of glucose. The reactivity parameters also indicated that the Wurtzite-glucose complex has a high electrophilicity index (ω), which means a good acceptor behavior to interact with glucose. Additionally, the bond between the (ZnO)n clusters and the glucose polar element has been studied in detail by using QTAIM theory. Finally, the theoretical and experimental studies prove that the Cu-Znn-1On clusters are very suitable and competent compounds for detecting glucose.

Density Functional Theory Study on NiNx (x = 1, 2, 3, 4) Catalytic Hydrogenation of Acetylene.

In this study, using the application of density functional theory, the mechanism of graphene-NiNx (x = 1, 2, 3, 4) series non-noble metal catalysts in acetylene hydrogenation was examined under the B3LYP/6-31G** approach. With the DFT-D3 density functional dispersion correction, the effective core pseudopotential basis set of LANL2DZ was applied to metallic Ni atoms. The reaction energy barriers of NiNx catalysts are different from the co-adsorption structure during the catalytic hydrogenation of graphene-NiNx (x = 1, 2, 3, 4). The calculated results showed that the energy barrier and selectivity of graphene-NiN4 for ethylene production were 25.24 kcal/mol and 26.35 kcal/mol, respectively. The low energy barrier and high activity characteristics showed excellent catalytic performance of the catalyst. Therefore, graphene-NiN4 provides an idea for the direction of catalytic hydrogenation.

New salts of teriflunomide (TFM) - Single crystal X-ray and solid state NMR investigation.

New salts of teriflunomide TFM (drug approved for Multiple Sclerosis treatment) with inorganic counterions: lithium (TFM_Li), sodium (TFM_Na), potassium (TFM_K), rubidium (TFM_Rb), caesium (TFM_Cs) and ammonium (TFM_NH4) were prepared and investigated employing solid state NMR Spectroscopy, Powder X-ray Diffraction PXRD and Single Crystal X-ray Diffraction (SC XRD). Crystal and molecular structures of three salts: TFM_Na (CCDC: 2173257), TFM_Cs (CCDC: 2165288) and TFM_NH4 (CCDC: 2165281) were determined and deposited. Compared to the native TFM, for all crystalline salt structures, a conformational change of the teriflunomide molecule involving about 180-degree rotation of the end group, forming an intramolecular hydrogen bond N-H⋯O is observed. By applying a complementary multi-technique approach, employing 1D and 2D solid state MAS NMR techniques, single and powder X-ray diffraction measurements, as well as the DFT-based GIPAW calculations of NMR chemical shifts for TFM_Na and TFM_Cs allowed to propose structural features of TFM_Li for which it was not possible to obtain adequate material for single crystal X-Ray measurement.

In Search of Preferential Macrocyclic Hosts for Sulfur Mustard Sensing and Recognition: A Computational Investigation through the New Composite Method r2SCAN-3c of the Key Factors Influencing the Host-Guest Interactions.

Sulfur mustard (SM) is a harmful warfare agent that poses a serious threat to human health and the environment. Thus, the design of porous materials capable of sensing and/or capturing SM is of utmost importance. In this paper, the interactions of SM and its derivatives with ethylpillar[5]arene (EtP[5]) and the interactions between SM and a variety of host macrocycles were investigated through molecular docking calculations and non-covalent interaction (NCI) analysis. The electronic quantum parameters were computed to assess the chemical sensing properties of the studied hosts toward SM. It was found that dispersion interactions contributed significantly to the overall complexation energy, leading to the stabilization of the investigated systems. DFT energy computations showed that SM was more efficiently complexed with DCMP[5] than the other hosts studied here. Furthermore, the studied macrocyclic containers could be used as host-based chemical sensors or receptors for SM. These findings could motivate experimenters to design efficient sensing and capturing materials for the detection of SM and its derivatives.

RNABPDB: Molecular Modeling of RNA Structure-From Base Pair Analysis in Crystals to Structure Prediction.

The stable three-dimensional structure of RNA is known to play several important biochemical roles, from post-transcriptional gene regulation to enzymatic action. These structures contain double-helical regions, which often have different types of non-canonical base pairs in addition to Watson-Crick base pairs. Hence, it is important to study their structures from experimentally obtained or even predicted ones, to understand their role, or to develop a drug against the potential targets. Molecular Modeling of RNA double helices containing non-canonical base pairs is a difficult process, particularly due to the unavailability of structural features of non-Watson-Crick base pairs. Here we show a composite web-server with an associated database that allows one to generate the structure of RNA double helix containing non-canonical base pairs using consensus parameters obtained from the database. The database classification is followed by an evaluation of the central tendency of the structural parameters as well as a quantitative estimation of interaction strengths. These parameters are used to construct three-dimensional structures of double helices composed of Watson-Crick and/or non-canonical base pairs. Our benchmark study to regenerate double-helical fragments of many experimentally derived RNA structures indicate very high accuracy. This composite server is expected to be highly useful in understanding functions of various pre-miRNA by modeling structures of the molecules and estimating binding efficiency. The database can be accessed from http://hdrnas.saha.ac.in/rnabpdb .

Influence of conventional hydrogen bonds in the intercalation of phenanthroline derivatives with DNA: The important role of the sugar and phosphate backbone.

The influence of hydrogen bonds in model intercalated systems between guanine-cytosine and adenine-thymine DNA base pairs (bps) was analyzed with the popular intercalator 1,10-phenanthroline (phen) and derivatives obtained by substitution with OH and NH2 groups in positions 4 and 7. Semiempirical and Density Functional Theory (DFT) methods were used both including dispersion effects: PM6-DH2, M06-2X and B3LYP-D3 along with the recently developed near linear-scaling coupled cluster method DLPNO-CCSD(T) for benchmark calculations. Our results given by QTAIM and non-covalent interaction analysis confirmed the existence of hydrogen bonds created by OH and NH2 . The trends in the energy decomposition analysis for the interaction energy, ΔEint , showed that the ΔEelstat contributions are equal or even a little bit higher than the values for ΔEdisp . Such important ΔEelstat attractive contribution comes mainly from the conventional hydrogen bonds formed by OH and NH2 functional groups with DNA not only with bps but specially with the sugar and phosphate backbone. This behavior is very different from that of phen and other classical intercalators that cannot form conventional hydrogen bonds, where the ΔEdisp is the most important attractive contribution to the ΔEint . The inclusion of explicit water molecules in molecular dynamics simulations showed, as a general trend, that the hydrogen bonds with the bps disappear during the simulations but those with the sugar and phosphate backbone remain in time, which highlights the important role of the sugar and phosphate backbone in the stabilization of these systems.

Synthesis and molecular modelling of thiadizole based hydrazone derivatives as acetylcholinesterase and butyrylcholinesterase inhibitory activities.

Some novel substituted thiazolylhydrazine derivatives were designed, synthesized and their inhibitory effects on acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE) enzymes and antioxidant activities were investigated. The structures of the synthesized compounds were determined using different spectroscopic techniques such as 1H-NMR, 13C-NMR, and HRMS. According to the enzyme inhibition results, the synthesized compounds showed selectivity against BuChE enzyme inhibition. Compounds 5e, 5g, 5i and 5j displayed significant BuChE inhibition potencies. Among them, compound 5i was found to be the most effective derivative with an IC50 value of 56.01 ± 0.054 µM. In addition, their antioxidant properties were evaluated in vitro through the 2,2-diphenyl-1-picrylhydrazyl (DPPH) assay. For compounds 5e, 5g, 5i and 5j in silico molecular docking and 100 ns molecular dynamics simulations studies against the BuChE enzyme were performed to determine possible protein-ligand interactions and stability. DFT-D3 study was performed to stabilize of compounds 5e, 5g, 5i and 5j both in gas and solvent medium and investigated their electronic properties. Of all geometries, that of DMSO is the lowest one.

Crystal structure from X-ray powder diffraction data, DFT-D calculation, Hirshfeld surface analysis, and energy frameworks of (RS)-trichlorme-thia-zide.

The structure of racemic (RS)-trichlorme-thia-zide [systematic name: (RS)-6-chloro-3-(di-chloro-meth-yl)-1,1-dioxo-3,4-di-hydro-2H-1λ6,2,4-benzo-thia-di-azine-7-sulfonamide], C8H8Cl3N3O4S2 (RS-TCMZ), a diuretic drug used in the treatment of oedema and hypertension, was determined from laboratory X-ray powder diffraction data using DASH [David et al. (2006 ▸). J. Appl. Cryst. 39, 910-915.], refined by the Rietveld method with TOPAS-Academic [Coelho (2018 ▸). J. Appl. Cryst. 51, 210-218], and optimized using DFT-D calculations. The extended structure consists of head-to-tail dimers connected by π-π inter-actions which, in turn, are connected by C-Cl⋯π inter-actions. They form chains propagating along [101], further connected by N-H⋯O hydrogen bonds to produce layers parallel to the ac plane that stack along the b-axis direction, connected by additional N-H⋯O hydrogen bonds. The Hirshfeld surface analysis indicates a major contribution of H⋯O and H⋯Cl inter-actions (32.2 and 21.7%, respectively). Energy framework calculations confirm the major contribution of electrostatic inter-actions (E elec) to the total energy (E tot). A comparison with the structure of S-TCMZ is also presented.

Removal of methylene blue from wastewater using ternary nanocomposite aerogel systems: Carboxymethyl cellulose grafted by polyacrylic acid and decorated with graphene oxide.

In this study, environmentally-friendly nanocomposite hydrogels were fabricated. These hydrogels consisted of semi-interpenetrating networks of carboxymethyl cellulose (CMC) molecules grafted to polyacrylic acid (PAA), as an eco-friendly and non-toxic polymer with numerous carboxyl and hydroxyl functional groups, which were reinforced with different levels of graphene oxide particles (0.5, 1.5 or 3% wt). Field-emission electron scanning microscopy (FESEM) images indicated that the pore size of the nanocomposites decreased with increasing graphic oxide concentration. The presence of the graphic oxide increased the storage modulus and thermal stability of the nanocomposite hydrogels. The hydrogels had an adsorption capacity of 138 mg/g of a model cationic dye pollutant (methylene blue) after 250 min. Moreover, a reusability test showed that the adsorption capacity remained at around 90% after 9 cycles. Density functional theory (DFT) simulations suggested that the adsorption of methylene blue was mainly a result of π-π bonds, hydrogen bonds, and electrostatic interactions with graphene oxide. Our results indicated that the nanocomposite hydrogels fabricated in this study may be eco-friendly, stable, efficient, and reusable adsorbents for ionic pollutants in wastewater treatment.

Solid-State Study of the Structure, Dynamics, and Thermal Processes of Safinamide Mesylate─A New Generation Drug for the Treatment of Neurodegenerative Diseases.

Safinamide mesylate (SM), the pure active pharmaceutical ingredient (API) recently used in Parkinson disease treatment, recrystallized employing water-ethanol mixture of solvents (vol/vol 1:9) gives a different crystallographic form compared to SM in Xadago tablets. Pure SM crystallizes as a hemihydrate in the monoclinic system with the P21 space group. Its crystal and molecular structure were determined by means of cryo X-ray crystallography at 100 K. SM in the Xadago tablet exists in anhydrous form in the orthorhombic crystallographic system with the P212121 space group. The water migration and thermal processes in the crystal lattice were monitored by solid-state NMR spectroscopy, differential scanning calorimetry, and thermogravimetric analysis. SM in Xadago in the high-humidity environment undergoes phase transformation to the P21 form which can be easily reversed just by heating up to 80 °C. For the commercial form of the API, there is also a reversible thermal transformation observed between Z' = 1 ↔ Z' = 3 crystallographic forms in the 0-20 °C temperature range. Analysis of molecular motion in the crystal lattice proves that the observed conformational polymorphism is forced by intramolecular dynamics. All above-mentioned processes were analyzed and described employing the NMR crystallography approach with the support of advanced theoretical calculations.

DFT-D4 Insight into the Inclusion of Amphetamine and Methamphetamine in Cucurbit[7]uril: Energetic, Structural and Biosensing Properties.

The host-guest interactions of cucurbit[7]uril (CB[7]) as host and amphetamine (AMP), methamphetamine (MET) and their enantiomeric forms (S-form and R-form) as guests were computationally investigated using density functional theory calculations with the recent D4 atomic-charge dependent dispersion corrections. The analysis of energetic, structural and electronic properties with the aid of frontier molecular orbital analysis, charge decomposition analysis (CDA), extended charge decomposition analysis (ECDA) and independent gradient model (IGM) approach allowed to characterize the host-guest interactions in the studied systems. Energetic results indicate the formation of stable non-covalent complexes where R-AMP@CB[7] and S-AMP@CB[7] are more stable thermodynamically than R-MET@CB[7] and S-MET@CB[7] in gas phase while the reverse is true in water solvent. Based on structural analysis, a recognition mechanism is proposed, which suggests that the synergistic effect of van der Waals forces, ion-dipole interactions, intermolecular charge transfer interactions and intermolecular hydrogen bonding is responsible for the stabilization of the complexes. The geometries of the complexes obtained theoretically are in good agreement with the X-ray experimental structures and indicate that the phenyl ring of amphetamine and methamphetamine is deeply buried into the cavity of CB[7] through hydrophobic interactions while the ammonium group remains outside the cavity to establish hydrogen bonds with the portal oxygen atoms of CB[7].

A Dispersion Corrected DFT Investigation of the Inclusion Complexation of Dexamethasone with ß-Cyclodextrin and Molecular Docking Study of Its Potential Activity against COVID-19.

The encapsulation mode of dexamethasone (Dex) into the cavity of ß-cyclodextrin (ß-CD), as well as its potential as an inhibitor of the COVID-19 main protease, were investigated using density functional theory with the recent dispersion corrections D4 and molecular docking calculations. Independent gradient model and natural bond orbital approaches allowed for the characterization of the host-guest interactions in the studied systems. Structural and energetic computation results revealed that hydrogen bonds and van der Waals interactions played significant roles in the stabilization of the formed Dex@ß-CD complex. The complexation energy significantly decreased from -179.50 kJ/mol in the gas phase to -74.14 kJ/mol in the aqueous phase. A molecular docking study was performed to investigate the inhibitory activity of dexamethasone against the COVID-19 target protein (PDB ID: 6LU7). The dexamethasone showed potential therapeutic activity as a SARS CoV-2 main protease inhibitor due to its strong binding to the active sites of the protein target, with predicted free energy of binding values of -29.97 and -32.19 kJ/mol as calculated from AutoDock4 and AutoDock Vina, respectively. This study was intended to explore the potential use of the Dex@ß-CD complex in drug delivery to enhance dexamethasone dissolution, thus improving its bioavailability and reducing its side effects.