Impact of 2-hydroxypropyl-ß-cyclodextrin inclusion complex formation on dopamine receptor-ligand interaction - A case study.

The octanol-water distribution coefficient (logP), used as a measure of lipophilicity, plays a major role in the drug design and discovery processes. While average logP values remain unchanged in approved oral drugs since 1983, current medicinal chemistry trends towards increasingly lipophilic compounds that require adapted analytical workflows and drug delivery systems. Solubility enhancers like cyclodextrins (CDs), especially 2-hydroxypropyl-ß-CD (2-HP-ß-CD), have been studied in vitro and in vivo investigating their ADMET (adsorption, distribution, metabolism, excretion and toxicity)-related properties. However, data is scarce regarding the applicability of CD inclusion complexes (ICs) in vitro compared to pure compounds. In this study, dopamine receptor (DR) ligands were used as a case study, utilizing a combined in silico/in vitro workflow. Media-dependent solubility and IC stoichiometry were investigated using HPLC. NMR was used to observe IC formation-caused chemical shift deviations while in silico approaches utilizing basin hopping global minimization were used to propose putative IC binding modes. A cell-based in vitro homogeneous time-resolved fluorescence (HTRF) assay was used to quantify ligand binding affinity at the DR subtype 2 (D2R). While all ligands showed increased solubility using 2-HP-ß-CD, they differed regarding IC stoichiometry and receptor binding affinity. This case study shows that IC-formation was ligand-dependent and sometimes altering in vitro binding. Therefore, IC complex formation can't be recommended as a general means of improving compound solubility for in vitro studies as they may alter ligand binding.

Polarity profiling of porous architectures: solvatochromic dye encapsulation in metal-organic frameworks.

Metal-organic frameworks (MOFs) have gathered significant interest due to their tunable porosity leading to diverse potential applications. In this study, we investigate the incorporation of the fluorosolvatochromic dye 2-butyl-5,6-dimethoxyisoindoline-1,3-dione ([double bond, length as m-dash]Phth) into various MOF structures as a means to assess the polarity of these porous materials. As a purely inorganic compound, zeolite Y was tested for comparison. The fluorosolvatochromic behavior of Phth, which manifests as changes in its emission spectra in response to solvent polarity, provides a sensitive probe for characterizing the local environment within the MOF pores. Through systematic variation of the MOF frameworks, we demonstrate the feasibility of using (fluoro-)solvatochromic dyes as probes for assessing the polarity gradients within MOF structures. Additionally, the fluorosolvatochromic response was studied as a function of loading amount. Our findings not only offer insights into the interplay between MOF architecture and guest molecule interactions but also present a promising approach for the rational design and classification of porous materials based on their polarity properties.

The impact of alcohol and ammonium fluoride on pressure-induced amorphization of cubic structure I clathrate hydrates.

We have investigated pressure-induced amorphization (PIA) of an alcohol clathrate hydrate (CH) of cubic structure type I (sI) in the presence of NH4F utilizing dilatometry and x-ray powder diffraction. PIA occurs at 0.98 GPa at 77 K, which is at a much lower pressure than for other CHs of the same structure type. The amorphized CH also shows remarkable resistance against crystallization upon decompression. While amorphized sI CHs could not be recovered previously at all, this is possible in the present case. By contrast to other CHs, the recovery of the amorphized CHs to ambient pressure does not even require a high-pressure annealing step, where recovery without any loss of amorphicity is possible at 120 K and below. Furthermore, PIA is accessible upon compression at unusually high temperatures of up to 140 K, where it reaches the highest degree of amorphicity. Molecular dynamics simulations confirm that polar alcoholic guests, as opposed to non-polar guests, induce cage deformation at lower pressure. The substitution of NH4F into the host-lattice stabilizes the collapsed state more than the crystalline state, thereby enhancing the collapse kinetics and lowering the pressure of collapse.

Probing the range of applicability of structure- and energy-adjusted QM/MM link bonds III: QM/MM MD simulations of solid-state systems at the example of layered carbon structures.

The previously introduced workflow to achieve an energetically and structurally optimized description of frontier bonds in quantum mechanical/molecular mechanics (QM/MM)-type applications was extended into the regime of computational material sciences at the example of a layered carbon model systems. Optimized QM/MM link bond parameters at HSEsol/6-311G(d,p) and self-consistent density functional tight binding (SCC-DFTB) were derived for graphitic systems, enabling detailed investigation of specific structure motifs occurring in graphene-derived structures v i a $$ via $$ quantum-chemical calculations. Exemplary molecular dynamics (MD) simulations in the isochoric-isothermic (NVT) ensemble were carried out to study the intercalation of lithium and the properties of the Stone-Thrower-Wales defect. The diffusivity of lithium as well as hydrogen and proton adsorption on a defective graphene surface served as additional example. The results of the QM/MM MD simulations provide detailed insight into the applicability of the employed link-bond strategy when studying intercalation and adsorption properties of graphitic materials.

Characterization of the Coordination and Solvation Dynamics of Solvated Systems─Implications for the Analysis of Molecular Interactions in Solutions and Pure H2O.

The characterization of solvation shells of atoms, ions, and molecules in solution is essential to relate solvation properties to chemical phenomena such as complex formation and reactivity. Different definitions of the first-shell coordination sphere from simulation data can lead to potentially conflicting data on the structural properties and associated ligand exchange dynamics. The definition of a solvation shell is typically based on a given threshold distance determined from the respective solute-solvent pair distribution function g(r) (i.e., GC). Alternatively, a nearest neighbor (NN) assignment based on geometric properties of the coordination complex without the need for a predetermined cutoff criterion, such as the relative angular distance (RAD) or the modified Voronoi (MV) tessellation, can be applied. In this study, the effect of different NN algorithms on the coordination number and ligand exchange dynamics evaluated for a series of monatomic ions in aqueous solution, carbon dioxide in aqueous and dichloromethane solutions, and pure liquid water has been investigated. In the case of the monatomic ions, the RAD approach is superior in achieving a well separated definition of the first solvation layer. In contrast, the MV algorithm provides a better separation of the NNs from a molecular point of view, leading to better results in the case of solvated CO2. When analyzing the coordination environment in pure water, the cutoff-based GC framework was found to be the most reliable approach. By comparison of the number of ligand exchange reactions and the associated mean ligand residence times (MRTs) with the properties of the coordination number autocorrelation functions, it is shown that although the average coordination numbers are sensitive to the different definitions of the first solvation shell, highly consistent estimates for the associated MRT of the solvated system are obtained in the majority of cases.

Modelling bulk and surface characteristics of cubic CeO2, Gd2O3, and gadolinium-doped ceria using a partial charge framework.

The development and characterization of materials for solid oxide fuel cells (SOFC) is an important step towards sustainable energy technologies. This present study models cubic CeO2, Gd2O3, and gadolinium-doped ceria (GDC) using newly constructed interaction potentials based on a partial atom charge framework. The interaction model was validated by comparing the structural properties with experimental reference data, which were found to be in good agreement. Validation of the potential model was conducted considering the surface stability of CeO2 and Gd2O3. Additionally, the accuracy of the novel potential model was assessed by comparing the oxygen diffusion coefficient in GDCn (n = 4-15) and the associated activation energy. The results demonstrate that the novel potential model is capable of describing the oxygen diffusion in GDC. In addition, this study compares the vibrational properties of the bulk with density functional theory (DFT) calculations, using a harmonic frequency analysis that avoids the need for computationally expensive quantum mechanical molecular dynamics (QM MD) simulations. The potential is compatible with a reactive water model, thus providing a framework for the simulation of solid-liquid interfaces.

Perylenetetracarboxylic Diimide Composite Electrodes as Organic Cathode Materials for Rechargeable Sodium-Ion Batteries: A Joint Experimental and Theoretical Study.

The organic semiconductor 3,4,9,10-perylenetetracarboxylic diimide (PTCDI), a widely used industrial pigment, has been identified as a diffusion-less Na-ion storage material, allowing for exceptionally fast charging/discharging rates. The elimination of diffusion effects in electrochemical measurements enables the assessment of interaction energies from simple cyclic voltammetry experiments through the theoretical work of Laviron and Tokuda. In this work, the two N-substituted perylenes, N,N'-dimethyl-3,4,9,10-perylenetetracarboxylic diimide (Me2PTCDI) and N,N'-diphenyl-3,4,9,10-perylenetetracarboxylic diimide (Ph2PTCDI), as well as the parent molecule 3,4,9,10-perylenetetracarboxylic diimide (H2PTCDI) are investigated as thin-film composite electrodes on carbon fibers for sodium-ion batteries. The composite electrodes are analyzed with Raman spectroscopy. Interaction parameters are extracted from cyclic voltammetry measurements. The stability and rate capability of the three PTCDI derivatives are examined through galvanostatic measurements in sodium-ion half-cell batteries and the influence of the interactions on those parameters is evaluated. In addition, self-consistent charge density function tight binding calculations of the different PTCDI systems interacting with graphite have been carried out. The results show that the binding motif displays notable deviations from an ideal ABA stacking, especially for the neutral state. In addition, data obtained for the electron-transfer integrals show that the difference in performance between different PTCDI thin-film batteries cannot be solely explained by the electron-transfer properties and other factors such as H-bonding have to be considered.

Towards hybrid quantum mechanical/molecular mechanical simulations of Li and Na intercalation in graphite - force field development and DFTB parametrisation.

In this work a previously established QM/MM simulation protocol for the treatment of solid-state interfaces was extended towards the treatment of layered bulk materials enabling for instance investigation of metal intercalation in graphitic carbon materials. In order to study the intercalation of Li in graphite, new density functional tight binding (DFTB) parameters for Li have been created. Molecular dynamics (MD) simulations at constant temperatures (273.15, 298.15 and 323.15 K) have been carried out to assess the performance of the presented DFTB MD simulation approach. The intercalation of variable lithium and sodium content was investigated via z-distribution functions and analysis of the diffusivity in the direction parallel to the graphene plane. Both the calculated diffusion coefficients and the activation energy in case of lithium are in good agreement with experimental data. The comparison of the QM/MM MD simulation results provide detailed insights into the structural and dynamical properties of intercalated metal ions.

Storage and Diffusion of Carbon Dioxide in the Metal Organic Framework MOF-5─A Semi-empirical Molecular Dynamics Study.

Metal-organic frameworks (MOFs) have attracted increasing attention due to their high porosity for exceptional gas storage applications. MOF-5 belongs to the family of isoreticular MOFs (IRMOFs) and consists of Zn4O6+ clusters linked by 1,4-benzenedicarboxylate. Due to the large number of atoms in the unit cell, molecular dynamics simulation based on density functional theory has proved to be too demanding, while force field models are often inadequate to model complex host-guest interactions. To overcome this limitation, an alternative semi-empirical approach using a set of approximations and extensive parametrization of interactions called density functional tight binding (DFTB) was applied in this work to study CO2 in the MOF-5 host. Calculations of pristine MOF-5 yield very good agreement with experimental data in terms of X-ray diffraction patterns as well as mechanical properties, such as the negative thermal expansion coefficient and the bulk modulus. In addition, different loadings of CO2 were introduced, and the associated self-diffusion coefficients and activation energies were investigated. The results show very good agreement with those of other experimental and theoretical investigations. This study provides detailed insights into the capability of semi-empirical DFTB-based molecular dynamics simulations of these challenging guest@host systems. Based on the comparison of the guest-guest pair distributions observed inside the MOF host and the corresponding gas-phase reference, a liquid-like structure of CO2 can be deduced upon storage in the host material.

Accessing Position Space Wave Functions in Band Structure Calculations of Periodic Systems─A Generalized, Adapted Numerov Implementation for One-, Two-, and Three-Dimensional Quantum Problems.

In this work, a generalized, adapted Numerov implementation capable of determining band structures of periodic quantum systems is outlined. Based on the input potential, the presented approach numerically solves the Schrödinger equation in position space at each momentum space point. Thus, in addition to the band structure, the method inherently provides information about the state functions and probability densities in position space at each momentum space point considered. The generalized, adapted Numerov framework provided reliable estimates for a variety of increasingly complex test suites in one, two, and three dimensions. The accuracy of the proposed methodology was benchmarked against results obtained for the analytically solvable Kronig-Penney model. Furthermore, the presented numerical solver was applied to a model potential representing a 2D optical lattice being a challenging application relevant, for example, in the field of quantum computing.

Terminal methylene phosphonium ions: precursors for transient monosubstituted phosphinocarbenes.

Isolable singlet carbenes are among the most important tools in chemistry, but generally require the interaction of two substituents with the electron deficient carbon atom. We herein report a synthetic approach to monosubstituted phosphinocarbenes via deprotonation of hitherto unknown diprotic terminal methylene phosphonium ions. Two methylene phosphonium salts bearing bulky N-heterocyclic imine substituents at the phosphorus atom were isolated and fully characterized. Deprotonation studies indicate the formation of transient monosubstituted carbenes that undergo intermolecular cycloadditions or intramolecular Buchner ring expansion to afford a cycloheptatriene derivative. The reaction mechanism of the latter transformation was elucidated using DFT calculations, which reveal the ambiphilic nature of the phosphinocarbene enabling the insertion into the aromatic C-C bond. Additional computational studies on the role of substituent effects are presented.

Beyond the Status Quo: Density Functional Tight Binding and Neural Network Potentials as a Versatile Simulation Strategy to Characterize Host-Guest Interactions in Metal- and Covalent Organic Frameworks.

In recent years, research focused on synthesis, characterization, and application of metal-organic frameworks (MOFs) has attracted increased interest, from both an experimental as well as a theoretical perspective. Self-consistent charge density functional tight binding (SCC DFTB) in conjunction with a suitable constrained molecular dynamics (MD) simulation protocol provides a versatile and flexible platform for the study of pristine MOFs as well as guest@MOF systems. Although being a semi-empirical quantum mechanical method, SCC DFTB inherently accounts for polarization and many-body contributions, which may become a limiting factor in purely force field-based simulation studies. A number of examples such as CO2, indigo, and drug molecules embedded in various MOF hosts are discussed to highlight the capabilities of the presented simulation approach. Furthermore, a promising extension of the outlined simulation strategy toward the treatment of covalent organic frameworks utilizing state-of-the-art neural network potentials providing a description at DFT accuracy and force field cost is outlined.

Deciphering Structural and Dynamical Properties of Hydrated Cobalt Porphyrins via Ab Initio Quantum Mechanical Charge Field Molecular Dynamics Simulation.

The present study successfully implemented the ab initio quantum mechanical charge field molecular dynamics (QMCF MD) formalism for the investigation of structural and dynamical properties of hydrated cobalt-porphyrin complexes. Considering the significance of cobalt ions in biological systems (for instance, vitamin B12), which reportedly incorporate cobalt ions in a d6, low spin, +3 state chelated in the corrin ring, an analog of porphyrin, the current study is focused on cobalt in the oxidation states +2 and +3 bound to the parent porphyrin lead structures embedded in an aqueous solution. These cobalt-porphyrin complexes were investigated in terms of their structural and dynamical properties at the quantum chemical level. The structural attributes of these hydrated complexes revealed the contrasting features of the water binding to these solutes, including a detailed evaluation of the associated dynamics. The study also yielded notable findings in regard to the respective electronic configurations vs coordination, which suggested that Co(II)-POR possesses a 5-fold square pyramidal coordination geometry in an aqueous solution containing the metal ion coordinating to four nitrogen atoms of the porphyrin ring and one axial water as the fifth ligand. On the other hand, high-spin Co(III)-POR was hypothesized to be more stable due to the smaller size-to-charge ratio of the cobalt ion, but the high-spin complex demonstrated unstable structural and dynamical behavior. However, the corresponding properties of the hydrated Co(III)LS-POR revealed a stable structure in an aqueous solution, thus suggesting the Co(III) ion to be in a low-spin state when bound to the porphyrin ring. Moreover, the structural and dynamical data were augmented by computing the free energy of water binding to the cobalt ions and the solvent-accessible surface area, which provide further information on thermochemical properties of the metal-water interaction and the hydrogen bonding potential of the porphyrin ring in these hydrated systems.

Storage and diffusion of CO2 in covalent organic frameworks-A neural network-based molecular dynamics simulation approach.

As a consequence of the accelerated climate change, solutions to capture, store and potentially activate carbon dioxide received increased interest in recent years. Herein, it is demonstrated, that the neural network potential ANI-2x is able to describe nanoporous organic materials at approx. density functional theory accuracy and force field cost, using the example of the recently published two- and three-dimensional covalent organic frameworks HEX-COF1 and 3D-HNU5 and their interaction with CO2 guest molecules. Along with the investigation of the diffusion behaviour, a wide range of properties of interest is analyzed, such as the structure, pore size distribution and host-guest distribution functions. The workflow developed herein facilitates the estimation of the maximum CO2 adsorption capacity and is easily generalizable to other systems. Additionally, this work illustrates, that minimum distance distribution functions can be a highly useful tool in understanding the nature of interactions in host-gas systems at the atomic level.

Structural Properties of Metal-Organic Frameworks at Elevated Thermal Conditions via a Combined Density Functional Tight Binding Molecular Dynamics (DFTB MD) Approach.

The performance of different density functional tight binding (DFTB) methods for the description of six increasingly complex metal-organic framework (MOF) compounds have been assessed. In particular the self-consistent charge density functional tight binding (SCC DFTB) approach utilizing the 3ob and matsci parameter sets have been considered for a set of four Zn-based and two Al-based MOF systems. Moreover, the extended tight binding for geometries, frequencies, and noncovalent interactions (GFN2-xTB) approach has been considered as well. In addition to the application of energy minimizations of the respective unit cells, molecular dynamics (MD) simulations at constant temperature and pressure conditions (298.15 K, 1.013 bar) have been carried out to assess the performance of the different DFTB methods at nonzero thermal conditions. In order to obtain the XRD patterns from the MD simulations, a flexible workflow to obtain time-averaged XRD patterns from (in this study 5000) individual snapshots taken at regular intervals over the simulation trajectory has been applied. In addition, the comparison of pair-distribution functions (PDFs) directly accessible from the simulation data shows very good agreement with experimental reference data obtained via measurements employing synchrotron radiation in case of MOF-5. The comparison of the lattice constants and the associated X-ray diffraction (XRD) patterns with the experimental reference data demonstrate, that the SCC DFTB approach provides a highly efficient and accurate description of the target systems.

First ternary tungsten tellurate(IV) WTe2O7 with unique crystal structure type.

At multianvil high-pressure/high-temperature conditions of 10 GPa and 1273 K, the first ternary tungsten tellurate WTe2O7 is formed, starting from a stoichiometric mixture of WO3 and TeO2. The compound crystallizes triclinic in a hitherto unknown crystal structure type with the space group P1̄; (no. 2), and was refined from single-crystal X-ray diffractometer data: a = 538.3(1), b = 687.5(1), c = 802.3(1) pm, α = 72.4(1)°, ß = 85.7(1)°, γ = 68.1(1)°, wR2 = 0.0323, GooF = 1.048, 3157 F2 values, and 106 variables. The main motifs of the crystal structure are pairs of edge-linked [WO6]6- octahedra and fourfold oxygen-coordinated Te4+ atoms. The oxidation state of W6+ and Te4+ was further verified by measuring the characteristic binding energy values for the W 4f and the Te 3d core levels via X-ray photoelectron spectroscopy (XPS). In addition, DFT calculations of the structure, the associated electron localisation functions (ELF) and vibrational spectra have been carried out. The theoretical data clearly demonstrates the impact of the residual electron density located at the Te4+ ions, which can be directly interpreted as the presence of lone electron pairs within the solid structure.

From vibrational spectroscopy and quantum tunnelling to periodic band structures - a self-supervised, all-purpose neural network approach to general quantum problems.

In this work, a feed-forward artificial neural network (FF-ANN) design capable of locating eigensolutions to Schrödinger's equation via self-supervised learning is outlined. Based on the input potential determining the nature of the quantum problem, the presented FF-ANN strategy identifies valid solutions solely by minimizing Schrödinger's equation encoded in a suitably designed global loss function. In addition to benchmark calculations of prototype systems with known analytical solutions, the outlined methodology was also applied to experimentally accessible quantum systems, such as the vibrational states of molecular hydrogen H2 and its isotopologues HD and D2 as well as the torsional tunnel splitting in the phenol molecule. It is shown that in conjunction with the use of SIREN activation functions a high accuracy in the energy eigenvalues and wavefunctions is achieved without the requirement to adjust the implementation to the vastly different range of input potentials, thereby even considering problems under periodic boundary conditions.

Direct Electrochemical CO2 Capture Using Substituted Anthraquinones in Homogeneous Solutions: A Joint Experimental and Theoretical Study.

Electrochemical capture of carbon dioxide (CO2) using organic quinones is a promising and intensively studied alternative to the industrially established scrubbing processes. While recent studies focused only on the influence of substituents having a simple mesomeric or nucleophilicity effect, we have systematically selected six anthraquinone (AQ) derivatives (X-AQ) with amino and hydroxy substituents in order to thoroughly study the influence thereof on the properties of electrochemical CO2 capture. Experimental data from cyclic voltammetry (CV) and UV-Vis spectroelectrochemistry of solutions in acetonitrile were analyzed and compared with innovative density functional tight binding computational results. Our experimental and theoretical results provide a coherent explanation of the influence of CO2 on the CV data in terms of weak and strong binding nomenclature of the dianions. In addition to this terminology, we have identified the dihydroxy substituted AQ as a new class of molecules forming rather unstable [X-AQ-(CO2) n ]2- adducts. In contrast to the commonly used dianion consideration, the results presented herein reveal opposite trends in stability for the X-AQ-CO2 •- radical species for the first time. To the best of our knowledge, this study presents theoretically calculated UV-Vis spectra for the various CO2-AQ reduction products for the first time, enabling a detailed decomposition of the spectroelectrochemical data. Thus, this work provides an extension of the existing classification with proof of the existence of X-AQ-CO2 species, which will be the basis of future studies focusing on improved materials for electrochemical CO2 capture.

Fused Bis(hemi-indigo): Broad-Range Wavelength-Independent Photoswitches.

Wavelength-independent conversion of organic photoswitches in the photostationary state is a rare phenomenon that opens up a way for many practical applications. In this work, three fused bis(hemi-indigo) derivatives with different substitution patterns were synthesized and their photoswitching was investigated by optical spectroscopy, real-time NMR spectroscopy and TD-DFT calculations. We disclosed that the Z-E photoisomerization of the meta-bis(hemi-indigo) derivative was remarkably independent of the irradiation wavelength from UV up to yellow light. The wavelength-independent forward photoswitching together with the inhibited backward photoisomerization, high thermal stability of the photoinduced isomers as well as significant overlap between the photoswitch absorption and the solar spectrum allows to suggest bis(hemi-indigo) derivatives as promising candidates for molecular solar thermal energy storage (MOST) systems.

Anthraquinone and its derivatives as sustainable materials for electrochemical applications - a joint experimental and theoretical investigation of the redox potential in solution.

Anthraquinone (AQ) has long been identified as a highly promising lead structure for various applications in organic electronics. Considering the enormous number of possible substitution patterns of the AQ lead structure, with only a minority being commercially available, a systematic experimental screening of the associated electrochemical potentials represents a highly challenging and time consuming task, which can be greatly enhanced via suitable virtual pre-screening techniques. In this work the calculated electrochemical reduction potentials of pristine AQ and 12 hydroxy- or/and amino-substituted AQ derivatives in N,N-dimethylformamide have been correlated against newly measured experimental data. In addition to the calculations performed using density functional theory (DFT), the performance of different semi-empirical density functional tight binding (DFTB) approaches has been critically assessed. It was shown that the SCC DFTB/3ob parametrization in conjunction with the COSMO solvation model provides a highly adequate description of the electrochemical potentials also in the case of the two-fold reduced species. While the quality in the correlation against the experimental data proved to be slightly inferior compared to the employed DFT approach, the highly advantageous cost-accuracy ratio of the SCC DFTB/3ob/COSMO framework has important implications in the formulation of hierarchical screening strategies for materials associated with organic electronics. Based on the observed performance, the low-cost method provides sufficiently accurate results to execute efficient pre-screening protocols, which may then be followed by a DFT-based refinement of the best candidate structures to facilitate a systematic search for new, high-performance organic electronic materials.