Diffusivity of α-, β-, γ-cyclodextrin and the inclusion complex of β-cyclodextrin: Ibuprofen in aqueous solutions; A molecular dynamics simulation study

Enhanced stability and biocompatibility of HIPEs stabilized by cyclodextrin-metal organic frameworks with inclusion of resveratrol and soy protein isolate for β-carotene delivery

High internal phase Pickering emulsions (HIPEs) constitute a significant research domain within colloid interface chemistry, addressing the demand for robust emulsion systems across various applications. An innovative nanoparticle, synthesized from a cyclodextrin metal-organic framework encapsulated with a composite of resveratrol and soy isolate protein (RCS), was employed to fortify a high internal phase emulsion. The emulsion's three-dimensional printing capabilities, alongside the encapsulated delivery efficacy for β-carotene, were thoroughly examined. Cyclodextrin metal-organic frameworks (CD-MOFs), facilitated by cellulose nanofibrils, were synthesized to yield particles at the nanoscale, maintaining a remarkable 97.67 % cellular viability at an elevated concentration of 1000 μg/ml. The RCS nanoparticles demonstrated thermal stability and antioxidant capacities surpassing those of CD-MOF. The integration of soybean isolate protein augmented both the hydrophobicity (from 21.95 ± 0.64° to 59.15 ± 0.78°) and the interfacial tension (from 14.36 ± 0.46 mN/m to 5.34 ± 0.81 mN/m) of the CD-MOF encapsulated with resveratrol, thereby enhancing the RCS nanoparticles' adsorption at the oil-water interface with greater stability. The durability of the RCS-stabilized high internal phase emulsions was contingent upon the RCS concentration. Emulsions stabilized with 5 wt%-RCS exhibited optimal physical and chemical robustness, demonstrating superior performance in emulsion 3D printing and β-carotene encapsulation delivery. This investigation furnishes a novel perspective on the amalgamation of food customization and precision nutrition.

Interactions of a Novel Anthracycline with Oligonucleotide DNA and Cyclodextrins in an Aqueous Environment

Berubicin, a chemotherapy medication belonging to the class of anthracyclines, is simulated in double-stranded DNA sequences and cyclodextrins in an aqueous environment via full-atom molecular dynamics simulations on the time scale of microseconds. The drug is studied in both the neutral and protonated states so as to better comprehend the role of its charge in the formed complexes. The noncovalent berubicin-DNA and berubicin-cyclodextrin complexes are investigated in detail, paying special attention to their thermodynamic description by employing the double decoupling method, the solvent balance method, the weighted solvent accessible surface model, and the linear interaction energy method. A novel approach for extracting the desolvation thermodynamics of the binding process is also presented. Both the binding and desolvation Gibbs energies are decomposed into entropic and enthalpic contributions so as to elucidate the nature of complexation and its driving forces. Selected structural and geometrical properties of all the complexes, which are all stable, are analyzed. Both cyclodextrins under consideration are widely utilized for drug delivery purposes, and a comparative investigation between their bound states with berubicin is carried out.

Computing Accurate True Self-Diffusion Coefficients and Shear Viscosities Using the OrthoBoXY Approach

In a recent paper [Busch, J.; Paschek, D. J. Phys. Chem. B 2023, 127, 7983-7987], we have shown that for molecular dynamics (MD) simulations using orthorhombic periodic boundary conditions with "magic" box length ratios of Lz/Lx = Lz/Ly = 2.7933596497, the self-diffusion coefficients Dx and Dy in x- and y-directions are independent of the system size. They both represent the true self-diffusion coefficient D0 = (Dx + Dy)/2, while the shear viscosity can be calculated from diffusion coefficients in x-, y-, and z-directions, using η = kBT·8.1711245653/[3πLz(Dx + Dy - 2Dz)]. In this contribution, we test this "OrthoBoXY" approach by its application to a variety of different systems: liquid water, dimethyl ether, methanol, triglyme, water/methanol mixtures, water/triglyme mixtures, and imidazolium-based ionic liquids. The chosen systems range from small-sized molecular liquids to complex mixtures and ionic liquids, while spanning a viscosity range of almost 3 orders of magnitude. We assess the efficiency of the method for computing true self-diffusion and viscosity data and provide simple formulas for estimating the required MD simulation lengths and sizes for delivering reliable data with targeted uncertainty levels. Our analysis of the system size dependence of statistical uncertainties for both the viscosity and the self-diffusion coefficient leads us to the conclusion that it is preferable to extend the simulation length instead of increasing the system size. MD simulations consisting of 768 molecules or ion pairs seem to be perfectly adequate.

An OrthoBoXY-method for various alternative box geometries

We have shown in a recent contribution [Busch and Paschek, J. Phys. Chem. B, 2023 127, 7983-7987] that for molecular dynamics (MD) simulations of isotropic fluids based on orthorhombic periodic boundary conditions with "magical" box length ratios of Lz/Lx = Lz/Ly = 2.7933596497, the computed self-diffusion coefficients Dx and Dy in x- and y-direction become system size independent. They thus represent the true self-diffusion coefficient D0 = (Dx + Dy)/2, while the shear viscosity can be determined from diffusion coefficients in x-, y-, and z-direction, using the expression η = kBT·8.1711245653/[3πLz(Dx + Dy - 2Dz)]. Here we present a more generalized version of this "OrthoBoXY"-approach, which can be applied to any orthorhombic MD box of any shape. In particular, we would like to test, how the efficiency is affected by using a shape more akin to the cubic form, albeit with different box length ratios Lx/Lz ≠ Ly/Lz and Lx < Ly < Lz. We use NVT and NpT simulations of systems of 1536 TIP4P/2005 water molecules as a benchmark and explore different box geometries to determine the influence of the box shape on the computed statistical uncertainties for D0 and η. Moreover, another "magical" set of box length ratios is discovered with Ly/Lz = 0.57804765578 and Lx/Lz = 0.33413909235, where the self-diffusion coefficient in x-direction becomes system size independent, such that D0 = Dx.

β-Cyclodextrin functionalized agarose-based hydrogels for multiple controlled drug delivery of ibuprofen

Agarose hydrogels are three-dimensional hydrophilic polymeric frameworks characterised by high water content, viscoelastic properties, and excellent ability as cell and drug delivery systems. However, their hydrophilicity as gel systems makes loading of hydrophobic drugs difficult and often ineffective. The incorporation of amphiphilic molecules (e.g. cyclodextrins) into hydrogels as hosts able to form inclusion complexes with hydrophobic drugs could be a possible solution. However, if not properly confined, the host compounds can get out of the network resulting in uncontrolled release. Therefore, in this work, β-cyclodextrins-based host-guest supramolecular hydrogel systems were synthesised, with β-cyclodextrins (β-CD) covalently bound to the polymeric network, preventing leakage of the host molecules. Hydrogels were prepared at two different β-CD-functionalized polyvinyl alcohol (PVA)/agarose ratios, and characterised chemically and physically. Then ibuprofen, a drug often used as a gold standard in studies involving β-CD both in its hydrophilic and hydrophobic forms, was selected to investigate the release behavior of the synthesised hydrogels and the influence of β-CD on the release. The presence of β-CD linked to the polymeric 3D network ensured a higher and prolonged release profile for the hydrophobic drug and also seemed to have some influence on the hydrophilic one.

Theoretical investigations of free energy of binding and chiral recognition studies of (R)- and (S)-Noradrenaline towards β-cyclodextrin

Noradrenaline (NA), one of the important excitatory catecholamine neurotransmitters, is used as a medication for Parkinson's Disease (PD). The β-cyclodextrin (β-CD) is one of the most effective drug carrier & also used in chiral separation. So, in this theoretical investigation, the R/S-Noradrenaline (R/S-NA) forms binding & chiral recognition mechanisms and energies with β-CD were explored. Using the AutoDock, R/S forms were first docked into the cavity of β-CD giving host-guest complexes with the free energy of binding for S-NA (-4.81 kcal/mol) larger than R-NA (-4.53 kcal/mol). The host-guest inclusion 1:1 complexes between R/S-NA and β-CD have been also modeled and optimized with ONIOM2 (B3LYP/6-31g++DP: PM6) method by using the Gaussian software. Further, frequency calculations were carried out to obtain the free energies. In comparison to the R-NA (-54.59 kcal/mol), it was observed that the S-NA (-56.48 kcal/mol) with β-CD is more stable. Furthermore, the H-bond results from molecular dynamics simulation revealed that S-NA/β-CD was more stable than R-NA/β-CD. In addition, the thermodynamic properties, vibrational analysis (IR), HOMO-LUMO band gap energy, inter molecular hydrogen bond interactions, and conformational analysis were investigated for both the R/S forms to support & compare the stability of the inclusion complex. These inclusion & high stability of S-NA/β-CD and in turn its theoretical chiral recognition behavior observed agreeing well with the reported NMR experimental data have implications in drug delivery and chiral separation research.

A comparative study on inclusion complex formation between formononetin and β-cyclodextrin derivatives through multiscale classical and umbrella sampling simulations

Formononetin, a naturally occurring isoflavone exhibits a wide range of therapeutic applications including antioxidant, anti-tumor, antiviral, anti-diabetic and neuroprotective activities. However, the low hydro-solubility of formononetin has limited its prospective use in cosmetic, neutraceutical and pharmaceutical industries. Cyclodextrins (CDs), especially β-CD and its derivatives have emerged as promising agents to improve the water solubility of poorly hydrosoluble compounds by the formation of inclusion complexes. We employed multiscale (1000 ns) explicit solvent and umbrella sampling molecular dynamics (MD) simulations to study the interactions and thermodynamic parameters of inclusion complex formation between formononetin and five most commonly used β-CD derivatives. Classical MD simulations revealed two possible binding conformations of formononetin inside the central cavity of hydroxypropyl-β-CD (HP-β-CD), randomly methylated-β-CD (ME-β-CD), and sulfobutylether-β-CD (SBE-β-CD). The binding conformation with the benzopyrone ring of formononetin inside the central cavity of β-CD derivatives was more frequent than the phenyl group occupying the hydrophobic cavity. These interactions were supported by a variety of non-bonded contacts including hydrogen bonds, pi-lone pair, pi-sigma, and pi-alkyl interactions. Formononetin showed favorable end-state MD-driven thermodynamic binding free energies with all the selected β-CD derivatives, except succinyl-β-CD (S-β-CD). Furthermore, umbrella sampling simulations were used to investigate the interactions and thermodynamic parameters of the host-guest inclusion complexes. The SBE-β-CD/formononetin inclusion complex showed the lowest binding energy signifying the highest affinity among all the selected host-guest inclusion complexes. Our study could be used as a standard for analyzing and comparing the ability of different β-CD derivatives to enhance the hydro-solubility of poorly soluble molecules.

Effect of temperature on the structure and drug-release behaviour of inclusion complex of β-cyclodextrin with cyclophosphamide: a molecular dynamics study

Cyclodextrins (CDs) are suitable drug carriers because of their doughnut-shaped cavities with hydrophilic outer and hydrophobic inner surfaces. Temperature-responsive CD-based drug carriers are expected to be one of the most promising candidates for drug delivery systems. In this study, we performed molecular dynamics simulations of the inclusion complex of β-CD with cyclophosphamide (CP) at temperatures from 300 K to 400 K to investigate the temperature dependency of the release behaviour of CP and structural changes of β-CD in an aqueous solution. We analysed the distance between the centres of mass of β-CD and CP and the radius of gyration of β-CD. The CP molecule was released from the β-CD cavity at 400 K, whereas two different inclusion complexes, partially and completely, were observed at T < 400 K. β-CD encapsulating a CP molecule had a more spherical shape and rigidity than β-CD without a CP, and the rigidity of their inclusion complex decreased with increasing temperature. Our findings provide fundamental insights into the behaviours of the β-CD/CP complex and drug release at the molecular level and can facilitate the development of new temperature-responsive drug delivery systems with CD nanocarriers triggered by localised temperature increases using focused ultrasound.

Molecular effects of site-specific phosphate-methylated primer on the structure and motions of Taq DNA polymerase

Polymerase chain reaction (PCR) is a powerful molecular biology assay for gene detection and quantification. Conventional DNA primers for PCR often suffer from poor sensitivity in specific gene detection. Recently, oligonucleotides containing methyl phosphotriester (MPTE-DNA) have been developed with enhanced DNA hybridization and improved gene detection sensitivity. Yet, site-specific MPTE-modifications on DNA primers have been reported to affect PCR amplification efficiencies while the detailed mechanism remains elusive. Here, we utilized molecular dynamics (MD) simulation to examine the effects of site-specific MPTE-modified primers on the structure and motions of DNA/Taq polymerase complexes. All tested MPTE-DNA/Taq complexes exhibited RMSD values below 2 Å, indicating insignificant effects of all methylation sites on the complex stability. The energetic and hydrogen-bonding analyses suggest minor effects of methylation at t-3, t-4, t-6, and t-7 positions on the DNA-Taq interaction. Principal component analyses further reveal that only t-3, and t-7 methylations preserve the motions of the Taq thumb domain. The site-specific methylation affects the interaction between DNA and the surrounding protein residues, resulting in allosteric-like effects on the DNA/Taq complex. The MD data complement the best experimentally observed PCR efficacies for the t-3 and t-7 positions among all tested MPTE-primers. The unveiled molecular insights can benefit the design of novel PCR primers for improving genetic testing platforms.

Molecular simulations of fluoxetine in hydrated lipid bilayers, as well as in aqueous solutions containing β-cyclodextrin

Fluoxetine, which is a well-known antidepressant drug, is studied in hydrated cholesterol-free and cholesterol-containing lipid bilayers through unbiased and biased atomistic molecular dynamics simulations. The latter are conducted for the calculation of the potential of mean force (PMF) of fluoxetine along an axis perpendicular to the two leaflets of the bilayer. The PMF indicates that the drug prefers to reside inside the lipid phase and allows us to calculate important thermodynamic properties, such as the Gibbs energy difference of partitioning from the water to the lipid phase and the Gibbs energy barrier for hopping events between the two leaflets of the bilayer. The results from the biased simulations are in accord with the mass density profiles calculated from the unbiased simulations. Moreover, we estimate the effect of fluoxetine mole fraction on the order parameters of the lipid alkyl chains and on the area per lipid. It is also found that fluoxetine forms a hydrogen bond network with lipids and water molecules penetrating into the lipid phase. In addition, fluoxoetine is studied in detail in aqueous solutions containing β-cyclodextrin. It is observed from unbiased molecular dynamics simulations that the two aforementioned molecules form a noncovalent complex spontaneously and the calculated binding free energy is in agreement with the literature.

Solubilities and Transport Properties of CO2, Oxalic Acid, and Formic Acid in Mixed Solvents Composed of Deep Eutectic Solvents, Methanol, and Propylene Carbonate

Recently, deep eutectic solvents (DES) have been considered as possible electrolytes for the electrochemical reduction of CO₂ to value-added products such as formic and oxalic acids. The applicability of pure DES as electrolytes is hindered by high viscosities. Mixtures of DES with organic solvents can be a promising way of designing superior electrolytes by exploiting the advantages of each solvent type. In this study, densities, viscosities, diffusivities, and ionic conductivities of mixed solvents comprising DES (i.e., reline and ethaline), methanol, and propylene carbonate were computed using molecular simulations. To provide a quantitative assessment of the affinity and mass transport of CO₂ and oxalic and formic acids in the mixed solvents, the solubilities and self-diffusivities of these solutes were also computed. Our results show that the addition of DES to the organic solvents enhances the solubilities of oxalic and formic acids, while the solubility of CO₂ in the ethaline-containing mixtures are in the same order of magnitude with the respective pure organic components. A monotonic increase in the densities and viscosities of the mixed solvents is observed as the mole fraction of DES in the mixture increases, with the exception of the density of ethaline-propylene carbonate which shows the opposite behavior due to the high viscosity of the pure organic component. The self-diffusivities of all species in the mixtures significantly decrease as the mole fraction of DES approaches unity. Similarly, the self-diffusivities of the dissolved CO₂ and the oxalic and formic acids also decrease by at least 1 order of magnitude as the composition of the mixture shifts from the pure organic component to pure DES. The computed ionic conductivities of all mixed solvents show a maximum value for mole fractions of DES in the range from 0.2 to 0.6 and decrease as more DES is added to the mixtures. Since for most mixtures studied here no prior experimental measurements exist, our findings can serve as a first data set based on which further investigation of DES-containing electrolyte solutions can be performed for the electrochemical reduction of CO₂ to useful chemicals.

Local Ion Densities can Influence Transition Paths of Molecular Binding

Improper reaction coordinates can pose significant problems for path-based binding free energy calculations. Particularly, omission of long timescale motions can lead to over-estimation of the energetic barriers between the bound and unbound states. Many methods exist to construct the optimal reaction coordinate using a pre-defined basis set of features. Although simulations are typically conducted in explicit solvent, the solvent atoms are often excluded by these feature sets—resulting in little being known about their role in reaction coordinates, and ultimately, their role in determining (un)binding rates and free energies. In this work, analysis is done on an extensive set of host-guest unbinding trajectories, working to characterize differences between high and low probability unbinding trajectories with a focus on solvent-based features, including host-ion interactions, guest-ion interactions and location-dependent ion densities. We find that differences in ion densities as well as guest-ion interactions strongly correlate with differences in the probabilities of reactive paths that are used to determine free energies of (un)binding and play a significant role in the unbinding process.

Chiral induction in a novel self-assembled supramolecular system composed of α-cyclodextrin porous liquids, chiral silver nanoparticles and planar conjugated molecules

The newly developed porous liquids known as liquids with permanent microporosity, have considerable application potential in many unknown areas. Herein, a supramolecular system composed of α-cyclodextrin porous liquid, chiral silver nanoparticles and planar conjugated molecules (methylene blue and indigo carmine) was designed and the induced chirality of the system was observed. It was found that the induced chirality can be easily tuned by changing the pH value of the mixture solution. The induced chiral signal of methylene blue in the developed self-assembled supramolecular system occurred when the pH was between 8 and 10, and furthermore the induced chirality of indigo carmine was found when the pH was between 6.5 and 7.5. The intensity of induced chirality decreases upon increasing temperatures and ionic strength. This study may offer a new approach for the creation of a chiral supramolecular system based on host-guest and electrostatic interaction and make cyclodextrin porous liquids promising candidates for applications in chiral induction.

How sensitive are physical properties of choline chloride-urea mixtures to composition changes: Molecular dynamics simulations and Kirkwood-Buff theory

Deep eutectic solvents (DESs) have emerged as a cheaper and greener alternative to conventional organic solvents. Choline chloride (ChCl) mixed with urea at a molar ratio of 1:2 is one of the most common DESs for a wide range of applications such as electrochemistry, material science, and biochemistry. In this study, molecular dynamics simulations are performed to study the effect of urea content on the thermodynamic and transport properties of ChCl and urea mixtures. With increased mole fraction of urea, the number of hydrogen bonds (HBs) between cation-anion and ion-urea decreases, while the number of HBs between urea-urea increases. Radial distribution functions (RDFs) for ChCl-urea and ChCl-ChCl pairs shows a significant decrease as the mole fraction of urea increases. Using the computed RDFs, Kirkwood-Buff Integrals (KBIs) are computed. KBIs show that interactions of urea-urea become stronger, while interactions of urea-ChCl and ChCl-ChCl pairs become slightly weaker with increasing mole fraction of urea. All thermodynamic factors are found larger than one, indicating a non-ideal mixture. Our results also show that self- and collective diffusivities increase, while viscosities decrease with increasing urea content. This is mainly due to the weaker interactions between ions and urea, resulting in enhanced mobilities. Ionic conductivities exhibit a non-monotonic behavior. Up to a mole fraction of 0.5, the ionic conductivities increase with increasing urea content and then reach a plateau.

Computational Insights Into the Influence of Substitution Groups on the Inclusion Complexation of β-Cyclodextrin

Cyclodextrins (CDs) and their derivatives have good prospects in soil remediation application due to their ability to enhance the stability and solubility of low water-soluble compounds by inclusion performance. To investigate the effect of different structural properties of cyclodextrin and its derivatives on the inclusion complexation, molecular dynamic (MD) simulations were performed on the inclusion complexes formed by three kinds of CDs with polycyclic aromatic hydrocarbons (PAHs). Based on neutral β-CD, the other two CDs were modified by introducing substitutional groups, including 2-hydroxypropyl and sulfonated butyl (SBE) functional groups in the ring structure, called HP-CD and SBE-CD. MD results show that PAH can merely enter into the cavity of SBE–β-CD from its wide rim. The substitutional groups significantly affect the structure of CDs, which may also cause the flipping of the glucose units. However, the substitutional groups can also enlarge the volume of the hydrophobic cavity, resulting in a tight combination with the guest molecules.

Click preparation of multiple-thioether bridged cyclodextrin chiral materials for efficient enantioseparation in high-performance liquid chromatography

A highly efficient covalent immobilization procedure is considered as an essential tool for obtaining stable and reliable cyclodextrin (CD) chiral stationary phases (CSPs). This work reports the "thiolene" click immobilization of heptakis(6-mercapto-6-deoxy)-β-CD-CSP onto alkene functional silica to afford novel multiple-thioether bridged CD CSPs by controlling the surface CD concentration. Solid-state NMR, FTIR, TGA and X-ray photoelectron diffraction spectroscopy (XPS) results proved the successful preparation of the desired CSPs with different surface CD loadings. The surface CD concentrations were calculated to be 0.49 and 0.68 μmol m-2 according to the elemental analysis results. More than 60 chiral enantiomers including isoxazolines, chiral lactides, chiral ketones, dansyl amino acids, small molecule acids and alkalis as well as some flavonoids were resolved or partially separated in the reversed-phase HPLC mode. Compared with the previously prepared single thiolene bridged CD-CSP, the current multiple-thioether CD-CSP afforded much better enantioseparation ability due to the existence of the thiol moiety and a confined structure.

Evaluating the Performance of Water Models with Host-Guest Force Fields in Binding Enthalpy Calculations for Cucurbit[7]uril-Guest Systems

Computational prediction of thermodynamic components with computational methods has become increasingly routine in computer-aided drug design. Although there has been significant recent effort and improvements in the calculation of free energy, the prediction of enthalpy (and entropy) remains underexplored. Furthermore, there has been relatively little work reported so far that attempts to comparatively assess how well different force fields and water models perform in conjunction with each other. Here, we report a comprehensive assessment of force fields and water models using host-guest systems that mimic many features of protein-ligand systems. These systems are computationally inexpensive, possibly because of their small size compared to protein-ligand systems. We present absolute enthalpy calculations using the multibox approach on a set of 25 cucurbit[7]uril-guest pairs. Eight water models were considered (TIP3P, TIP4P, TIP4P-Ew, SPC, SPC/E, OPC, TIP5P, Bind3P), along with five force fields commonly used in the literature (GAFFv1, GAFFv2, CGenFF, Parsley, and SwissParam). We observe that host-guest binding enthalpies are strongly sensitive to the selection of force field and water model. In terms of water models, we find that TIP3P and its derivative Bind3P are the best performing models for this particular host-guest system. The performance is generally better for aliphatic compounds than for aromatic ones, suggesting that aromaticity remains a difficult property to include accurately in these simple force fields.