Insights into the interactions and dynamics of a DES formed by phenyl propionic acid and choline chloride

Pharmaceutical applications of therapeutic deep eutectic systems (THEDES) in maximising drug delivery

The issue of poor solubility of active pharmaceutical ingredients (APIs) has been a salient area of investigation and novel drug delivery systems are being developed to improve the solubility of drugs, enhance their permeability and thereby their efficacy. Several techniques for solubilization enhancement of poorly soluble drugs are often employed at various stages of pharmaceutical drug product development. One such delivery system is the therapeutic deep eutectic system (THEDES), which showed great potential in the enhancement of solubility and permeability of drugs and ultimately augmenting their bioavailability. THEDES are made by mixing drugs with deep eutectic solvents (DESs) in a definite molar ratio by the hit and trial method. The DESs are a new class of green solvents which are non-toxic, cheap, easy to prepare, biodegradable and have multiple applications in the pharmaceutical industry. The terminologies such as ionic liquids (ILs), DES, THEDES, and therapeutic liquid eutectic systems (THELES) have been very much in use recently, and it is important to highlight the pharmaceutical applications of these unexplored reservoirs in drug solubilization enhancement, drug delivery routes, and in the management of various diseases. This review is aimed at discussing the components, formulation strategies, and routes of administration of THEDES that are used in developing the formulation. Also, the major pharmaceutical applications of THEDES in the treatment of various metabolic and non-metabolic diseases are reviewed.

Biophysical assessment of amantadine and SDS surfactant mixture onto boron nitride nanotube: a molecular dynamics investigation

CONTEXT

The aggregation and adsorption of amantadine and sodium dodecyl sulfate on the boron nitride nanotubes in aqueous system were investigated, employing the classical molecular dynamics method. The aqueous system containing amantadine and sodium dodecyl sulfate was investigated by self-diffusion of the molecules, with particular emphasis on their center of mass mean square displacement. The radial distribution function and dipole moment measurement were used to analyze the water effect on the aggregation and molecular interaction between amantadine and sodium dodecyl sulfate in the complex. In turn, the amantadine molecules are able to develop water monolayers at the aqueous solution but show no noticeable aggregated adsorption in the fluid. In contrast, sodium dodecyl sulfate self-aggregation on the boron nitride nanotube efficiently occurs, and so the effect of sodium dodecyl sulfate on aggregated adsorption of amantadine was studied. Our results show that in the presence of sodium dodecyl sulfate at critical micelle concentration, its effect on the aggregation of amantadine was enhanced onto boron nitride nanotube compared to just pure sodium dodecyl sulfate or amantadine on boron nitride nanotube. Our simulations show a remarkable restructuring and enhanced orientation of the water molecules around sodium dodecyl sulfate and amantadine adsorbed on boron nitride nanotube.

METHOD

All molecular dynamics simulations were done using NAMD-2.9 package with CHARMM-36 force field. The nanotube PDB file was made by VMD with (12, 12) indexes. The nanotube atoms were considered to be fixed and periodic during the simulation process, and the ends did not cap with hydrogens. The initial PDB file of components was saved from ChemOffice software, and the webservers to obtain the charge and topology files were Reddb and Swiss Param webservers. The VMD software was used for structural visualizations and graphical analysis.

Mechanistic insights into the lignin dissolution behavior in amino acid based deep eutectic solvents

Deep eutectic solvents (DESs) composed by amino acids (L-arginine, L-proline, L-alanine) as the hydrogen bond acceptors (HBAs) and carboxylic acids (formic acid, acetic acid, lactic acid, levulinic acid) as hydrogen bond donors (HBDs) were prepared and used for the dissolution of dealkaline lignin (DAL). The mechanism of lignin dissolution in DESs was explored at molecular level by combining the analysis of Kamlet-Taft (K-T) solvatochromic parameters, FTIR spectrum and density functional theory (DFT) calculations of DESs. Firstly, it was found that the formation of new hydrogen bonds between lignin and DESs mainly drove the dissolution of lignin, which were accompanied by the erosion of hydrogen bond networks in both lignin and DESs. The nature of hydrogen bond network within DESs was fundamentally determined by the type and number of functional groups in both HBA and HBD, which affected its ability to form hydrogen bond with lignin. One hydroxyl group and carboxyl group in HBDs provided active protons, which facilitated proton-catalyzed cleavage of β-O-4, thus enhancing the dissolution of DESs. The superfluous functional group resulted in more extensive and stronger hydrogen bond network in the DESs, thus decreasing the lignin dissolving ability. Moreover, it was found that lignin solubility had a closed positive correlation with the subtraction value of α and β (net hydrogen donating ability) of DESs. Among all the investigated DESs, L-alanine/formic acid (1:3) with the strong hydrogen-bond donating ability (acidity), weak hydrogen-bond accepting ability (basicity) and small steric-hindrance effect showed the best lignin dissolving ability (23.99 wt%, 60 °C). On top of that, the value of α and β of L-proline/carboxylic acids DESs showed some positive correlation with the global electrostatic potential (ESP) maxima and minima of the corresponding DESs respectively, indicating the analysis of ESP quantitative distributions of DESs could be an effective tool for DESs screening and design for lignin dissolution as well as other applications.

Liquid Structures and Diffusion Dynamics of Alkyl-Pyrene Liquids Studied by Molecular Dynamics Simulations

Functional molecular liquids (FMLs) based on alkylated π-conjugated molecules have attracted attention as solvent-free and nonvolatile liquid materials with prominent optoelectronic features. Recently, novel FML compounds containing pyrene as the functional core were synthesized, and their rheological and photochemical properties were investigated. Although the molecules differ only in the number of alkyl chain substituents and their substitution positions, their viscosity coefficients are largely different beyond the Stokes-Einstein relation on the assumption of identical microscopic friction, indicating that local microscopic molecular interactions are crucial for the macroscopic rheological properties. Here, we report a theoretical study on the rheological properties of the alkyl-pyrene liquids by means of atomistic molecular dynamics (MD) simulations. We performed long-time MD simulations for tens of microseconds to obtain ample statistical samples of the alkyl-pyrene liquids and analyzed their liquid structures and diffusion dynamics based on spatiotemporal correlation functions. We found the formation of characteristic local liquid structures of π-π stacking of the pyrene moieties and locally anisotropic and anomalous diffusion dynamics, which remarkably vary depending on the alkyl substituent patterns. The present results provide an atomistic insight into the macroscopic rheological properties of alkyl-π FMLs and molecular design strategy for them.

Analysis of the microscopic interactions between processed Polygonatum cyrtonema polysaccharides and water

The dissolution and microscopic interactions of processed Polygonatum cyrtonema polysaccharides in water are extremely important because they strongly influence the process to extract these polysaccharides from water. In this paper, molecular dynamics simulation methods were used to analyse the influence of extraction temperature, concentration and molecular weight on the radial distribution function (RDF), mean square displacement (MSD), diffusion coefficient (D), radius of gyration (R_g), and microstructure of processed Polygonatum cyrtonema polysaccharides in water as well as the intrinsic viscosity (η), hydrogen bond characteristics and microscopic interactions in the solutions. The research results showed that the extraction temperature, concentration and molecular weight of the polysaccharides had important effects on the RDF, MSD, D, R_g, η, hydrogen bond characteristics and the microstructure of the polysaccharide molecules, but there were some major differences. The order of the influence of the factors affecting the dissolution of polysaccharides in water was temperature > concentration > molecular weight. Extraction temperature, material fluid ratio and molecular weight had greater influence on the fluidity and dissolution state of the polysaccharides in water. As the solute concentration and molecular weight increased, hydrogen bonds between polysaccharides and water inhibited their dissolution and diffusion. Properly increasing the temperature, reducing the concentration and selecting low molecular weight polysaccharides enhanced the dissolution and diffusion of the polysaccharides in the solution system. Molecular weight was the main factor affecting the η of the polysaccharide solutions. These results can provide theoretical guidance for the extraction or tea brewing process of Polygonatum cyrtonema polysaccharides in future work.

Ternary CE2Ba2 (E = As, Sb) Clusters: New Pentaatomic Planar Tetracoordinate Carbon Species with 18 Valence Electrons

18-valence-electron (ve) rule is one important guide for us to design planar tetracoordinate carbon (ptC) species. Using the "polarization of ligands" strategy, the new pentaatomic ptC species CE₂Ba₂ (E = As, Sb) with 18 ve are designed in this work. Computer structural searches and high-level calculations reveal that the ptC CE₂Ba₂ (E = As, Sb) species are global minima (GMs) on the potential energy surfaces, whose C center is coordinated by the interspaced E and Ba atoms. CE₂Ba₂ (E = As, Sb) are also kinetically stable. Chemical bonding analyses reveal that the ptC core is stabilized by two localized C-E σ bonds, one delocalized five-center two-electron (5c-2e) σ bond and one delocalized 5c-2e π bond. One π and three σ bonds collectively conform to the 8-electron counting, which determines the stability of ptC CE₂Ba₂ (E = As, Sb) species. Interestingly, the delocalized 2π and 2σ electrons render the ptC systems π/σ double aromaticity. Additional 10 electrons contribute to peripheral lone pairs of E and E-Ba bonding.

Computer Simulations of Deep Eutectic Solvents: Challenges, Solutions, and Perspectives

Deep eutectic solvents (DESs) are one of the most rapidly evolving types of solvents, appearing in a broad range of applications, such as nanotechnology, electrochemistry, biomass transformation, pharmaceuticals, membrane technology, biocomposite development, modern 3D-printing, and many others. The range of their applicability continues to expand, which demands the development of new DESs with improved properties. To do so requires an understanding of the fundamental relationship between the structure and properties of DESs. Computer simulation and machine learning techniques provide a fruitful approach as they can predict and reveal physical mechanisms and readily be linked to experiments. This review is devoted to the computational research of DESs and describes technical features of DES simulations and the corresponding perspectives on various DES applications. The aim is to demonstrate the current frontiers of computational research of DESs and discuss future perspectives.

Fine-Tuning the Polarizable CL&Pol Force Field for the Deep Eutectic Solvent Ethaline

Polarizable force fields are gradually becoming a common choice for ionic soft matter, in particular, for molecular dynamics (MD) simulations of ionic liquids (ILs) and deep eutectic solvents (DESs). The CL&Pol force field introduced in 2019 is the first general, transferable, and polarizable force field for MD simulations of different types of DESs. The original formulation contains, however, some problems that appear in simulations of ethaline and may also have a broader impact. First, the originally proposed atomic diameter parameters are unbalanced, resulting in too weak interactions between the chlorides and the hydroxyl groups of the ethylene glycol molecules. This, in turn, causes an artificial phase separation in long simulations. Second, there is an overpolarization of chlorides due to strong induced dipoles that give rise to the presence of peaks and antipeaks at very low q-vector values (2.4 nm^-1) in the partial components of the structure factors. In physical terms, this is equivalent to overestimated spatial nanoscale heterogeneity. To correct these problems, we adjusted the chloride-hydroxyl radial distribution functions against ab initio data and then extended the use of the Tang-Toennis damping function for the chlorides' induced dipoles. These adjustments correct the problems without losing the robustness of the CL&Pol force field. The results were also compared with the nonpolarizable version, the CL&P force field. We expect that the corrections will facilitate reliable use of the CL&Pol force field for other types of DESs.

Molecular and Spectroscopic Insights into a Metal Salt-Based Deep Eutectic Solvent: A Combined Quantum Theory of Atoms in Molecules, Noncovalent Interaction, and Density Functional Theory Study

Deep eutectic solvents (DESs) based on metal halide salts are highly catalytic, low toxic, reusable, cost-effective, and have higher thermal stability than their analogue ionic liquids (ILs). In this work, we have reported the formation mechanism of metal salt-based DESs at the molecular level along with their charge-transfer analysis and thermodynamics associated with their formation using density functional theory. The DES systems analyzed in the present work were choline chloride and tin(II)chloride (DES1) and choline chloride and zinc(II)chloride (DES2), both in a molar ratio of 1:2, respectively. An excellent correlation is obtained between the theoretically calculated IR spectra of the DES systems and the previously reported experimental findings for the formation of the complex systems. The DESs were found to be stable systems due to traditional hydrogen bonding and electrostatic interactions resulting in the ionic species [Sn₂Cl₅]^- and [Zn₂Cl₅]^- and are elucidated with the help of electronic structure calculations. CHELPG partial charge analysis and natural bond orbital analysis suggest a charge transfer from Cl^- (chloride) to Ch⁺ (choline) and metal salts in the DES structures. The atom-in-molecules and noncovalent interaction (NCI) analysis suggest a strong electrostatic interaction within the DES2 system as compared to DES1. Higher stability and reactivity are observed in the DES2 system based on the frontier molecular orbital analysis. Our analysis offers important insights into the formation mechanism of these economic IL analogues.

Formulation design and mechanism study of hydrogel based on computational pharmaceutics theories

Formulation design and mechanism study of the drug delivery system (DDS) is an important but difficult subject in pharmaceutical research. The study of formulation factors is the most time- and labor-consuming work of formulation design. In this paper, a multiscale computational pharmaceutics strategy was developed to guide the systematic study of formulation factors of a typical polymer-based DDS, hydrogel, and further to guide the formulation design. According to the strategy, the combination of solubility parameter (δ) and diffusion coefficient (D) calculated by the AA-MD simulation was suggested as the general evaluation method for the matrix screening of the hydrogels at the pre-formulation stage. At the formulation design stage, the CG-MD simulation method was suggested to predict the morphology and drug-releasing behavior of the hydrogels under different formulation factors. The influence mechanism can be explained by the combination of multiple parameters, such as the microstructure diagram, the radius of gyration (Rg), the radial distribution function (RDF), and the free diffusion volume (V_diffusion). The simulation results are in good agreement with the in vitro release experiment, indicating that the strategy has good applicability.