Molecular Dynamics Study of β-Cyclodextrin–Phenylalanine (1:1) Inclusion Complex in Aqueous Medium

Evaluating solubility, stability, and inclusion complexation of oxyresveratrol with various β-cyclodextrin derivatives using advanced computational techniques and experimental validation

Oxyresveratrol (OXY), a natural stilbenoid in mulberry fruits, is known for its diverse pharmacological properties. However, its clinical use is hindered by low water solubility and limited bioavailability. In the present study, the inclusion complexes of OXY with β-cyclodextrin (βCD) and its three analogs, dimethyl-β-cyclodextrin (DMβCD), hydroxypropyl-β-cyclodextrin (HPβCD) and sulfobutylether-β-cyclodextrin (SBEβCD), were investigated using in silico and in vitro studies. Molecular docking revealed two binding orientations of OXY, namely, 4',6'-dihydroxyphenyl (A-form) and 5,7-benzenediol ring (B-form). Molecular Dynamics simulations suggested the formation of inclusion complexes with βCDs through two distinct orientations, with OXY/SBEβCD exhibiting maximum atom contacts and the lowest solvent-exposed area in the hydrophobic cavity. These results corresponded well with the highest binding affinity observed in OXY/SBEβCD when assessed using the MM/GBSA method. Beyond traditional simulation methods, Ligand-binding Parallel Cascade Selection Molecular Dynamics method was employed to investigate how the drug enters and accommodates within the hydrophobic cavity. The in silico results aligned with stability constants: SBEβCD (2060 M-1), HPβCD (1860 M-1), DMβCD (1700 M-1), and βCD (1420 M-1). All complexes exhibited a 1:1 binding mode (AL type), with SBEβCD enhancing OXY solubility (25-fold). SEM micrographs, DSC thermograms, FT-IR and 1H NMR spectra confirm the inclusion complex formation, revealing novel surface morphologies, distinctive thermal behaviors, and new peaks. Notably, the inhibitory impact on the proliferation of breast cancer cell lines, MCF-7, exhibited by inclusion complexes particularly OXY/DMβCD, OXY/HPβCD, and OXY/SBEβCD were markedly superior compared to that of OXY alone.

Update on chiral recognition mechanisms in separation science

The stereospecific analysis of chiral molecules is an important issue in many scientific fields. In separation sciences, this is achieved via the formation of transient diastereomeric complexes between a chiral selector and the selectand enantiomers driven by molecular interactions including electrostatic, ion-dipole, dipole-dipole, van der Waals or π-π interactions as well as hydrogen or halogen bonds depending on the nature of selector and selectand. Nuclear magnetic resonance spectroscopy and molecular modeling methods are currently the most frequently applied techniques to understand the selector-selectand interactions at a molecular level and to draw conclusions on the chiral separation mechanism. The present short review summarizes some of the recent achievements for the understanding of the chiral recognition of the most important chiral selectors combining separation techniques with molecular modeling and/or spectroscopic techniques dating between 2020 and early 2024. The selectors include polysaccharide derivatives, cyclodextrins, macrocyclic glycopeptides, proteins, donor-acceptor type selectors, ion-exchangers, crown ethers, and molecular micelles. The application of chiral ionic liquids and chiral deep eutectic solvents, as well as further selectors, are also briefly addressed. A compilation of all published literature on chiral selectors has not been attempted.

Molecular level investigation of interactions between pharmaceuticals and β-cyclodextrin (β-CD) functionalized adsorption sites for removal of pharmaceutical contaminants from water

This study describes a novel application of the use of molecular modeling tools for investigating the adsorption of organic micropollutants (OMPs) from water by nanocomposites. The partitioning of pharmaceuticals onto β-Cyclodextrin (β-CD) functionalized adsorbents was investigated at the molecular level to explore the atomistic interactions of pharmaceutical contaminants in water systems with β-CD and to provide insight into possible approaches for removal of pharmaceuticals from water. Molecular electrostatic surface potential mapping of β-CD derivatives was employed to examine the impact of substitution degree of β-CD and type of grafting agent on host-guest complexation. The stability of the complexes of selected pharmaceuticals and β-CD derivatives were assessed via molecular dynamics simulations to evaluate competitive adsorption between organic micropollutants (OMPs) and between OMPs and fulvic acid, a representative natural organic material (NOM) component found in water systems. Molecular electrostatic surface potential maps showed that grafting agents with aromatic and amine functional groups were found to provide attractive interactions for negatively charged OMPs. In addition, optimization of substitution degree of β-CD is necessary to enhance adsorption of target OMPs. Furthermore, it was found that aromatic ring bearing grafting agents can provide additional electrostatic attractions by π-π interactions with the aromatic ring of the OMPs. The impact of common water quality characteristics on adsorption was assessed and it was revealed that the effect of pH and calcium on adsorption depends on the ionizable functional groups present on the grafting agent. Molecular dynamics simulations showed that adsorption of target OMPs does not solely depend on hydrophobicity but is affected by electrostatic interactions. The simulations revealed that fulvic acid which is commonly present in environmental waters can be a competitor with ibuprofen for the β-CD cavity. Ultimately, this study showed that molecular level simulation can be effectively employed to investigate adsorption of OMPs by β-CD functionalized adsorbents and could be employed to enhance their design and subsequent environmental applications.

Long-Lived States Provide Insights from NMR into the β-Cyclodextrin Drug Assemblies

In the last two decades, extending spin memory in NMR has been used for several purposes. Long-lived states (LLS) or singlet states are one of the first spin memory enhancement techniques used. LLS have the potential to extract structural information and intra- and intermolecular interactions of complex systems other than studying slow phenomenon. The motional regime of β-cyclodextrin (β-CD) drug inclusion complexes generally lies in the intermediate region, where ωτ_c ≈ 1, and the standard methods of studying these interactions, i.e., NOE and chemical shift monitoring, suffer from insufficient output information. The sensitivity of LLS toward the environmental changes is utilized here to gain insights into the drug assemblies formed by β-CD. One can use change in relaxation of LLS to study the structural changes during complexation. The examples of β-CD with the drugs indomethacin, paracetamol, gliclazide, and CI-933 (a precursor 4-methoxybenzamide) were studied. Indomethacin, paracetamol, and 4-methoxybenzamide show strong interaction through the para-substituted benzene ring, unlike gliclazide. Relaxation of LLS in β-CD-drug complexes is modeled using standard Redfield Relaxation Theory. Computational studies performed support the experimental observations. Docking and molecular dynamics simulation provided the explanation of the relaxation properties of these drug molecules.

Dynamic gelatin-based hydrogels promote the proliferation and self-renewal of embryonic stem cells in long-term 3D culture

Long-term maintenance of embryonic stem cells (ESCs) in the undifferentiated state is still challenging. Compared with traditional 2D culture methods, 3D culture in biomaterials such as hydrogels is expected to better support the long-term self-renewal of ESCs by emulating the biophysical and biochemical properties of the extracellular matrix (ECM). Although prior studies showed that soft and degradable hydrogels favor the 3D growth of ESCs, few studies have examined the impact of the structural dynamics of the hydrogel matrix on ESC behaviors. Herein, we report a gelatin-based structurally dynamic hydrogel (GelCD hydrogel) that emulates the intrinsic structural dynamics of the ECM. Compared with covalently crosslinked gelatin hydrogels (GelMA hydrogels) with similar stiffness and biodegradability, GelCD hydrogels significantly promote the clonal expansion and viability of encapsulated mouse ESCs (mESCs) independent of MMP-mediated hydrogel degradation. Furthermore, GelCD hydrogels better maintain the pluripotency of encapsulated mESCs than do traditional 2D culture methods that use MEF feeder cells or medium supplementation with GSK3β and MEK 1/2 inhibitors (2i). When cultured in GelCD hydrogels for an extended period (over 2 months) with cell passaging every 7 days, mESCs preserve their normal morphology and maintain their pluripotency and full differentiation capability. Our findings highlight the critical role of the structural dynamics of the hydrogel matrix in accommodating the volume expansion that occurs during clonal ESC growth, and we believe that our dynamic hydrogels represent a valuable tool to support the long-term 3D culture of ESCs.

Formation and molecular dynamics simulation of inclusion complex of large-ring cyclodextrin and 4-terpineol

In the present study, the formation and structure of the inclusion compound of large-ring cyclodextrin and 4-terpineol were obtained through different experiments and molecular dynamics (MD) simulation. The analysis of FTIR, ¹ H-NMR, and thermodynamic results confirmed the formation of clathrates. Analysis of molecular structure (root-mean-square deviation and radius of gyration), solubility, and interaction energy (Coul, H bond) based on MD simulations further clarified the nature of the clathrate and the conformational changes caused by guest molecules as well as inclusion complexes process trends. The inclusion complex reportedly has a new crystal structure with improved thermal stability. PRACTICAL APPLICATION: This is the first work to demonstrate the complex formation between 4-terpineol and large-ring cyclodextrin by molecular dynamics simulation. Molecular dynamics simulation confirmed the formation of inclusion complexes theoretically. Conformational changes of the molecules and the formation of complexes with improved thermal stability were observed. Complexing with large-ring cyclodextrin can be used as an effective means to encapsulate the aroma/flavor compounds.

In Search of the Most Stable Molecular Configuration of Heptakis(2,6-O-dimethyl)-β-cyclodextrin and Its Complex with Mianserin: A Comparison of the B3LYP-GD2 and M062X-GD3 Results

Cyclodextrins are well known for their ability to form stable, highly soluble complexes with various substances, which makes them widely used as excipients in food, cosmetics, and pharmaceuticals. In this work, properties of heptakis(2,6-O-dimethyl)-β-cyclodextrin (DM-β-CD) in vacuo and in water, as well as its ability to bind the antidepressant drug mianserin (MIA) in aqueous solution, are investigated computationally. The results are shown to depend strongly on the density functional theory (DFT) applied. The most stable conformers of DM-β-CD found with the B3LYP-GD2 method differ from these indicated by M062X-GD3 and other functionals. According to the latter, two crystal structures, ZULQAY and BOYFOK03, optimized in vacuo and in water, respectively, have the lowest energy. Both the B3LYP-GD2 and M062X-GD3 results show that all tested inclusion and noninclusion complexes of MIA:DM-β-CD in stoichiometry 1:1 are stable in water. However, the structures and their energetic properties obtained with each method differ: in the most stable configurations, different aromatic rings of MIA are embedded inside DM-β-CD, and the corresponding complexation energies (calculated with the 6-31++G(d,p) basis set and corrected for the basis set superposition error) are -29.6 (B3LYP-GD2) and -23.9 (M062X-GD3) kcal/mol. The NMR spectra of DM-β-CD and MIA:DM-β-CD are also compared.

Application of Molecular Dynamics Simulations in the Analysis of Cyclodextrin Complexes

Cyclodextrins (CDs) are highly respected for their ability to form inclusion complexes via host-guest noncovalent interactions and, thus, ensofance other molecular properties. Various molecular modeling methods have found their applications in the analysis of those complexes. However, as showed in this review, molecular dynamics (MD) simulations could provide the information unobtainable by any other means. It is therefore not surprising that published works on MD simulations used in this field have rapidly increased since the early 2010s. This review provides an overview of the successful applications of MD simulations in the studies on CD complexes. Information that is crucial for MD simulations, such as application of force fields, the length of the simulation, or solvent treatment method, are thoroughly discussed. Therefore, this work can serve as a guide to properly set up such calculations and analyze their results.

Mimicking the Endometrial Cancer Tumor Microenvironment to Reprogram Tumor-Associated Macrophages in Disintegrable Supramolecular Gelatin Hydrogel

Purpose

Besides the tumor cells themselves, solid tumors are comprised of numerous cell types including infiltrating immune cells such as tumor-associated macrophages (TAMs). TAMs are vital stromal components of host immune system and play a critical role in the development of cancer. TAMs can be divided into two subtypes: M1 tumor-suppressive macrophage and M2 tumor-supportive macrophage. To better address the observations of TAMs functional performance, we describe an in vitro system that mimics the populations of TAMs infiltrated into the tumor mass by using our disintegrable supramolecular gelatin (DSG) hydrogels, which are physically crosslinked by host-guest complexations.

Materials and Methods

The host–guest interaction was adopted between the aromatic groups of gelatin and the photocrosslinkable acrylated β-cyclodextrins (Ac-β-CDs) to form the DSG hydrogels. The convenient macrophage/endometrial cancer cells heterospheroid 3D model was set up by DSG hydrogels. RT-PCR and Western blot assays were developed to evaluate the efficiencies of inducers on the macrophages. The ELISA and oxygen saturation assays were performed to measure the secretion of VEGF and consumption of oxygen of tumor and/or macrophages, respectively. To determine the antitumor effects of M2 reprogrammed macrophages in vitro and in vivo, migration assay and tumor xenograft model were used, respectively.

Results

The host-guest complexations of DSG hydrogels were controllably broken efficiently by soaking into the solution of competitive guest monomers 1-adamantanamine hydrochloride. The DSG hydrogels help IFN-γ reprogram the M2 to M1 and then decrease the tumor/M2 reprogrammed macrophage cells heterospheroid secretion of VEGF and increase the relative oxygen saturation. Significantly, the co-cultural tumor/M2 reprogrammed group from the disintegrated DSG hydrogels reduced the migration of cancer cells in vitro and the tumor growth in vivo.

Conclusion

We obtain a TAMs/tumor microenvironment-responsive 3D model based on the novel DSG hydrogels, and will be of utility in cancer therapy and drug discovery.

Effect of Photoacid Strength on Fluorescence Modulation of 2-Naphthol Derivatives inside β-Cyclodextrin Cavity: Insights from Fluorescence, Isothermal Calorimetry, and Molecular Dynamics Simulations

Fluorescence response of a photoacid inside a confined environment often differs markedly from the bulk response. Is there any correlation between the extent of fluorescence modulation and the strength of the photoacid? Here, we used three photoacids: 2-naphthol (2OH, pK_a* = 3.3), 6-sulfonate-2-naphthol (6SO₃-2OH, pK_a* = 3.06), and 6-cyano-2-naphthol (6CN-2OH, pK_a* = 0.6) with remarkably different excited-state acidities to investigate fluorescence modulation inside the nanocavity of β-cyclodextrin (β-CD). Interestingly, we found strong fluorescence modulation for 2OH and 6SO₃-2OH but almost none for 6CN-2OH. Isothermal calorimetry measurements showed that all three fluorophores form 1:1 inclusion complex with comparable binding constants (285, 420, and 580 M^-1 for 2OH, 6SO₃-2OH, and 6CN-2OH, respectively). Molecular dynamics simulation further revealed that binding modes are quite similar, and the distribution of water molecules around the proton-donating hydroxyl group of the photoacids are also comparable. Consequently, the difference in the fluorescence response should be accounted solely to the difference in the photoacidity strengths.

Elucidating the H-Bonding Environment of Coumarin 102 in a Phenol-Cyclohexane Mixture by Molecular Dynamics Simulation: Implications for H-Bond-Guided Photoinduced Electron Transfer

Recently, we have experimentally demonstrated that the fluorescence intensity of coumarin 102 (C102) modulates anomalously upon hydrogen bonding to phenol in a nonpolar solvent: cyclohexane. The fluorescence intensity is first quenched gradually up to a particular mole fraction (X_PH ≈ 0.013) but thereafter increases with further increases in the phenol mole fraction. These studies speculate about the importance of C102-phenol H-bonding to induce photoinduced electron transfer (PET) and propose a competition between the C102-phenol and phenol-phenol H-bonding to account for the anomalous fluorescence modulation. In this work, we investigate the exact H-bonding environment around the acceptor C102 at various compositions by molecular dynamics simulation and correlate the H-bonding environment to the observed fluorescence quenching. In addition to the 1:1 C102-phenol complex, 1:2 C102-(phenol)₂ complexes with two different types of geometries were also found. Furthermore, density functional theory (DFT) and time-dependent density functional theory (TDDFT) calculations were carried out to understand the H-bonding in these complexes in the ground state and in the excited state and their possible contribution to the observed fluorescence quenching.

Factors Affecting the Formation of 2:1 Host:Guest Inclusion Complexes of 2-[(R-Phenyl)amine]-1,4-naphthalenediones (PAN) in β- and γ-Cyclodextrins

The molecular hosts cyclodextrins form inclusion complexes with a wide variety of guests, resulting in complexes with various host:guest stoichiometries. In the case of a series of 19 1,4-naphthoquinolines as guests with either β- or γ-cyclodextrin studied using electrospray mass spectroscopy, in most cases only 1:1 complexes were observed, with 2:1 host:guest complexes observed in just 6 out of 38 host:guest combinations. It is shown that these higher-order complexes were observed only in the case of small (or no) electronically withdrawing substituents, and were much less likely in the case of the larger γ-cyclodextrin host. The size and electronic properties of the substituents involved shows that both steric and electronic factors must be taken into account in predicting which cyclodextrin host:guest stoichiometries will be stable enough to form (or once formed, be robust enough to be observed in the ESI-MS experiments). It is clear that the prediction of host-guest stoichiometry for a specific host-guest pair is complicated, and involves a subtle interplay of both electronic and steric factors. However, there are definite trends, which can be used to help predict host:guest stoichiometry for a given host-guest pair.

Conformational Change in the Mechanism of Inclusion of Ketoprofen in β-Cyclodextrin: NMR Spectroscopy, Ab Initio Calculations, Molecular Dynamics Simulations, and Photoreactivity

Inclusion of drugs in cyclodextrins (CDs) is a recognized tool for modifying several properties such as solubility, stability, bioavailability, and so on. The photoreactive behavior of the β-CD/ketoprofen (KP) complex upon UV exposure showed a significant increase in photodecarboxylation, whereas the secondary degradation products by hydroxylation of the benzophenone moiety were inhibited. The results may account for an improvement of KP photophysical properties upon inclusion, thus better fostering its topical use. To correlate the structural details of the inclusion with these results, an NMR spectroscopic study of KP upon inclusion in β-CD was performed. Effects of the magnetically anisotropic centers of KP, changing their orientations upon inclusion and giving chemical shift variations, were specifically correlated with the results of the molecular dynamic simulations and ab initio calculations. In the large variety of papers focusing on the structural analysis of β-CD complexes, this work represents one of the few examples in which a detailed analysis of these simultaneous upfield-downfield NMR shifts of the same aromatic molecule upon inclusion is reported. Interestingly, the results demonstrate that the observed upfield and downfield shifts upon inclusion are not related to any direct magnetic role of β-CD. The conformational change of KP upon the inclusion process consists of a slight reduction in the angle between the two phenyl rings and in a remarkable reduction in the mobility of the carboxyl group, the latter being one of the main contributions to the NMR resonance shifts. These structural details help in understanding the features of the inclusion complex and, eventually, the driving force for its formation.

β-Cyclodextrin at the Water/1-Bromobutane Interface: Molecular Insight into Reverse Phase Transfer Catalysis

Molecular insight into the role of β-cyclodextrin (βCD) as a phase transfer catalyst at the liquid/liquid interface is obtained by molecular dynamics simulations of the structure and dynamics of βCD adsorbed at the interface between water and 1-bromobutane. In particular, we consider the structure and dynamics of the water and bromobutane molecules inside the βCD cavity and compare them with the behavior when βCD is dissolved in bulk water. βCD is preferentially oriented at the interface, with the cavity opening along the interface normal. While in bulk water the cavity includes 6-8 water molecules that are relatively mobile with short residence time, at the interface the cavity is mostly dehydrated and includes a single bromobutane molecule. This inclusion complex is stable in bulk water. The implication of this behavior for reverse phase transfer catalysis is discussed.

In silico understanding of the cyclodextrin-phenanthrene hybrid assemblies in both aqueous medium and bacterial membranes

The explicit-solvent molecular dynamic (MD) simulation and adaptive biased forces (ABF) methods were employed to systemically study the structural and thermodynamic nature of the β-cyclodextrin (βCD) monomer, phenanthrene (Phe) monomer, and their inclusion complexes in both the aqueous and membrane environments, aiming at clarifying the atomic-level mechanisms underlying in the CD-enhanced degradation of polycyclic aromatic hydrocarbons (PAHs) by bacteria. Simulations showed that βCD and Phe monomers could associate together to construct two distinctive assemblies, i.e, βCD1-Phe1 and βCD2-Phe1, respectively. The membrane-involved equilibrium simulations and the data of potential of mean forces (PMFs) further confirmed that Phe monomer was capable of penetrating through the membranes without confronting any large energy barrier, whereas, the single βCD and βCD-involved assemblies were unable to pass across the membranes. These observations clearly suggested that βCD only served as the carrier to enhance the bioavailability of Phe rather than the co-substrate in the Phe biodegradation process. The Phe-separation PMF profiles indicated that the maximum of the Phe uptake by bacteria would be achieved by the "optimal" βCD:Phe molar ratio, which facilitated the maximal formation of βCD1-Phe1 inclusion and the minimal construction of βCD2-Phe1 complex.

Resolving Internal Motional Correlations to Complete the Conformational Entropy Meter

Conformational entropy (SΩ) has long been used to theoretically characterize the dynamics of proteins, DNA, and other polymers. Though recent advances enabled its calculation also from simulations and nuclear magnetic resonance (NMR) relaxation experiments, correlated molecular motion has hitherto greatly hindered both numerical and experimental determination, requiring demanding empirical and computational calibrations. Herein, we show that these motional correlations can be estimated directly from the temperature-dependent SΩ series that reveal effective persistence lengths of the polymers, which we demonstrate by measuring SΩ of amphiphilic molecules in model lipid systems by spin-labeling electron paramagnetic resonance (EPR) spectroscopy. We validate our correlation-corrected SΩ meter against the basic biophysical interactions underlying biomembrane formation and stability, against the changes in enthalpy and diffusion coefficients upon phase transitions, and against the energetics of fatty acid dissociation. As the method can be directly applied to conformational analysis of proteins and other polymers, as well as adapted to NMR or polarized fluorescence techniques, we believe that the approach can greatly enrich the scope of experimentally available statistical thermodynamics, offering new physical insights into the behavior of biomolecules.

Sensitivity of local hydration behaviour and conformational preferences of peptides to choice of water model

Secondary structural preferences of the beta-hairpin of the 2GB1 protein in the folded and unfolded ensembles are shown to be sensitive to the choice of water model.