Due to its protective properties of biological samples at low temperatures and under desiccation, dimethyl sulfoxide (DMSO) in aqueous solutions has been studied widely by many experimental approaches and molecular dynamics (MD) simulations. In the case of the latter, AMBER is among the most commonly used force fields for simulations of biomolecular systems; however, the parameters for DMSO published by Fox and Kollman in 1998 have only been tested for pure liquid DMSO. We have conducted an MD simulation study of DMSO in a water mixture and computed several structural and dynamical properties such as of the mean density, self-diffusion coefficient, hydrogen bonding and DMSO and water ordering. The AMBER force field of DMSO is seen to reproduce well most of the experimental properties of DMSO in water, with the mixture displaying strong and specific water ordering, as observed in experiments and multiple other MD simulations with other non-polarizable force fields.
In this article, we explore the capacity of formed Schiff base complexes to trap metal atoms or ions, using their aromatic ends. The intrinsic geometry of each complex defines the process of substitution. Two cases were studied; one involving a trans Schiff base complex and the other considering how a salen ligand, with nickel systems traps chromium. We also assessed the nature of the new bonds and the frontier molecular orbitals.
The aldol reaction in the presence of L-proline acting as an organocatalyst is a well-known example of asymmetric synthesis. Many theoretical and experimental studies have been carried out to probe the mechanism of this reaction. In this work, two levels of density functional theory in the gas phase and DMSO were used to elucidate the best pathways for this reaction, with the enamine and enol considered intermediates and L-proline considered either a reactant or a facilitator. The calculations indicated that both intermediates are formed simultaneously in the reaction medium. Interestingly, the formation of the enamine intermediate predominates in DMSO at room temperature, whereas the enol becomes the predominant intermediate upon the addition of water.
Unknown force-field parameters for metal organic beryllium complexes used in emitting and electron transporting layers of OLED structures are determined. These parameters can be used for the predictive atomistic simulations of the structure and properties of amorphous organic layers containing beryllium complexes. The parameters are found for the AMBER force field using a relaxed scan procedure and quantum-mechanical DFT calculations of potential energy curves for specific internal (angular) coordinates in a series of three Be complexes (Bebq2; Be(4-mpp)2; Bepp2). The obtained parameters are verified in calculations of some molecular and crystal structures available from either quantum-mechanical DFT calculations or experimental data.
The density functional theory method using the B3LYP/6-31G(d,p) level of theory was used to perform isoenergetic maps in order to determine the lower energy conformers of four disaccharides constituting alginic acids, which are based on β-D-mannuronic (M) and α-L-guluronic acid (G), called MM, GG, MG, and GM. The preferred structures are combined to monovalent (Li+, Na+, and K+) cations and further fully optimized, and an isoenergetic map corresponding to the complex (MG2?, 2Na+) was performed. Then, the reactivity of MG complexes with mono- and bivalent cations was studied using the global nucleophilic index. The position selectivity was also predicted using the local nucleophilic indices. It was demonstrated that experimental trends of relative reactivity and regioselectivity of the complexes are correctly predicted using these empirical indices of reactivity.
Magnetic shielding constants for an isolated fullerene C60, cucurbituril CB[9], and the host-guest complex C60@CB[9] were calculated as a function of separation of the monomers. Our results in the gas phase and water indicate a significant variation of the magnetic properties for all atoms of the monomers in the complex and after liberation of fullerene C60 from the interior of the CB[9] cavity. The interaction between the two monomers results in a charge transfer that collaborates with a redistribution of electron density to deshield the monomers.
A theoretical 1H NMR spectroscopy and thermodynamic analysis of the host–guest inclusion process involving the norfloxacin (NFX) into β-cyclodextrin (β-CD) was carried out. DFT structure and stabilization energies were obtained in both gas and aqueous phases. We could establish that the complex formation is enthalpy driven, and the hydrogen bonds established between NFX and β-CD play a major role in the complex stabilization. Besides, a theoretical 1H NMR analysis has shown to be a supplementary proceeding to predict appropriately the inclusion mode of norfloxacin molecule into the β-CD. In this work, a theoretical study of the NFX@β-CD complex is reported for the first time, seeking a deep understanding of topology and thermodynamics of the inclusion complex formation.
Mechanisms for the activation of water, ammonia, and other small molecules by the PCcarbeneP nickel pincer complex were studied computationally with the aid of density functional theory. The calculation results indicate that the strongly donating, nucleophilic carbene center can engage in a variety of heterolytic splitting of E?H (E=H, C, N, O) bonds, some of which are reversible. The cleavage of E?H bonds across the Ni=C bond represents a new mode of bond activation by ligand cooperativity in nickel pincer complex. On the basis of the calculations, we also demonstrate that reversible H2 activation across the Ir=C bond via the PCcarbeneP iridium pincer complex was observed in the experiments, while other E?H (E=C, N, O) bonds were not activated. Our calculations are in good agreement with experimental observations and could provide new insights into ligand cooperativity in nickel pincer complexes.
The absorption and emission spectra of dichlorvos and the dichlorvos-MAA complex in methanol, water, and chloroform in the molecularly imprinted recognition were investigated systematically. The M06-2X results revealed that: 1) the hydroxyl groups in polar solvents such as methanol and water may markedly influence the weak interactions, and then alter the adsorption and emission spectra; 2) the electronic excitation in absorption spectra of dichlorvos is dominated by the configuration HOMO?→?LUMO, but in the most stable dichlorvos-MAA it becomes the ππ* excitation of HOMO?→?LUMO?+?1; 3) Mulliken charges reveal that dichlorvos almost dissociates to Cl- and a cation in its S1 excitation state; 4) the phosphorescence spectra of dichlorvos-MAA are relatively weak.
Graphical Abstract The absorption and emission spectra of dichlorvos and the dichlorvos-MAA complex in the molecularly imprinted recognition of dichlorvos were investigated systematically in methanol, water, and chloroform as solvents.
Quantum chemical computations were used for prediction of the structure and color of alizarin complex with alkali metal hydroxides in methanolic solutions. The color prediction relying on the single Gaussian-like band once again proved the usefulness of the PBE0 density functional due to the observed smallest color difference between computed and experimentally derived values. It was found that the alkali metal hydroxide molecules can bind to the two oxygen atoms of both hydroxyl groups of alizarin or to one of these atoms and the oxygen atom from the keto group in a complex with three methanol molecules. This means that two electronic transitions need to be taken into account when considering the spectra of the studied complexes. The resulting bond lengths and angles are correlated with the properties of the alkali metal atoms. The molar mass, the atomic radius, and the Pauling electronegativity of studied metals are quite accurate predictors of the geometric properties of hydroxide complexes with alizarin in methanol solution.
Graphical abstract The spectra of the neutral and monoanionic form of alizarin together with color changes resulting from addition of different metal hydroxides and represented in CIE color space
The solubility advantage (SA) of meloxicam cocrystalized with mono- and dicarboxylic acids was expressed in terms of equilibrium constants involving active pharmaceutical ingredient and coformer in aqueous solutions. It is argued that SA can be quantified by concentration of pairs formed in water. The pH and concentration of dissolved components is included explicitly in the model. The alternative behavior of mono- and dicarboxylic acids was emphasized and addressed to different structural motifs. The structural and energetic properties of meloxicam and its complexes with carboxylic acids were characterized, including tautmerism and dissociation in aqueous media. In particular, performed in silico modeling confirmed experimental observation that meloxicam dissolved in water or modest acidic solutions is expected to be a mixture of anionic form in equilibrium with at least five neutral isomers. Tautomer-related diversity of pairs formation and the possibility of salt formation is also discussed.
In this study, the doped defects in nitromethane crystals were investigated using first-principles calculations for the first time. We introduce dopant atoms in the interstitial sites of the nitromethane lattice, aiming to study the effects of element-doping on the structural properties, electronic properties, and sensitivity characteristics. The obtained results show that doped defects obviously affect the neighboring nitromethane molecules. The modification of electronic properties shows that the band gaps are significantly influenced by doped defects. Partial density of states and population analysis further reveal the mechanism for sensitivity control of nitromethane. It is shown that the new electronic states were introduced in the forbidden bands and the doped defects resulted in charge redistributions in the systems.
Present molecular dynamics simulations indicate that the methanol component in a methanol/water mixture is more likely to be trapped in a cyclic peptide nanotube (CPNT), while water molecules tend to be present at the channel mouths as transient guests. Channel water resides mainly between methanol and the CPNT wall, resulting in a distinct decrease in the H-bond number per channel methanol. Six designed CPNTs with different channel diameters and outer surface characteristics all possess distinct selectivity to methanol over water. Of these, the amphipathic 8?×?(AQ)4-CPNT exhibits the best performance. Results in this study provide basic information for the application of a CPNT to enrich methanol from a methanol/water mixture.
The factors that explain the competition between intramolecular NO linkage photoisomerization and NO photorelease in five ruthenium nitrosyl complexes were investigated. By applying DFT-based methods, it was possible to characterize the ground states and lowest triplet potential energy surfaces of these species, and to establish that both photoisomerization and photorelease processes can occur in the lowest triplet state of each species. This work highlights the crucial role of the sideways-bonded isomer, a metastable state also known as the MS2 isomer, in the photochemical loss of NO, while the results obtained also indicate that the population of the triplet state of this isomer is compulsory for both processes and show how photoisomerization and photorelease interfere.
Bond critical points (BCPs) in the quantum theory of atoms in molecules (QTAIM) are shown to be a consequence of the molecular topology, symmetry, and the Poincaré-Hopf relationship, which defines the numbers of critical points of different types in a scalar field. BCPs can be induced by a polarizing field or by addition of a single non-bonded atom to a molecule. BCPs and their associated bond paths are therefore suggested not to be a suitable means of identifying chemical bonds, or even attractive intermolecular interactions.
The interaction of external water molecules with hydrated pyrrole-2-carboxaldehyde PCL/(H2O)n complexes was investigated. The work was supported by both theoretical [DFT/TD-DFT methods using 6-311G++(d,p) basis set in the ground (S0) and excited (S1, S2, S3)states] and experimental [UV-Vis, FTIR and Raman] verification. The focus of the present work was on the weak intermolecular O–H?O, N–H?O–H hydrogen bonded interaction (IerHB) between PCL and external water molecules, and the influence of increasing the number of water molecules to form hydrated PCL/(H2O)n complexes. Effects were observed on different vibrational normal modes and on electronic transition levels. A hydrogen-bonded network of water induces a shift to higher energy in certain normal modes of PCL to form stable PCL/(H2O)n complexes by lowering the barrier energy. Potential energy distribution (PED) analysis indicates a significant charge transfer from PCL to water by creating a water bridge. Hydrogen bonding effects account for the substantial red shift and broadness in νNH, νCO vibrational modes. Water rearrangement turns out to be the main driving force for hydrated complex formation.
The structure and stability of various ternary complexes in which an extended aromatic system such as coronene interacts with ions/atoms/molecules on opposite faces of the π-electron cloud were investigated using ab initio calculations. By characterizing the nature of the intermolecular interactions using an energy decomposition analysis, it was shown that there is an interplay between various types of interactions and that there are co-operativity effects, particularly when different types of interactions coexist in the same system.
Graphical abstract Weak OH-π, π-π and van der Waals-π ternary systems are stabilized through dispersion interactions. Cation-π ternary systems are stabilized by through-space electrostatic interactions.
The electronic structures and spectroscopic properties of two complexes [M(pic)3] (M = Ir, Rh) containing picolinate as bidentate ligands have been calculated by means density functional theory (DFT) and time-dependent DFT/TD-DFT using three hybrid functionals B3LYP, PBE0 and mPW1PW91. The PBE0 and mPW1PW91 functionals, which have the same HF exchange fraction (25%), give similar results and do not differ drastically from B3LYP results. Calculated geometric parameters of the complexes are in good agreement with the available experimental data. The UV absorptions observed in acetonitrile were assigned on the basis of singlet state transitions. The most intense band observed in the UV-C region corresponds to ligand-to-ligand charge transfer states (LLCT) in both complexes. The theoretical spectrum of the rhodium complex is characterized by a large degree of mixing between metal-to-ligand-charge-transfer (MLCT), LLCT and metal centered (MC) states in the UV-A region. The presence of low-lying excited states with MC character affects the absorption spectrum under spin-orbit coupling (SOC) effects and play important roles in the photochemical properties.
