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.
We investigated the interaction (hyper)polarizability of neon–dihydrogen pairs by performing high-level ab initio calculations with atom/molecule-specific, purpose-oriented Gaussian basis sets. We obtained interaction-induced electric properties at the SCF, MP2, and CCSD levels of theory. At the CCSD level, for the T-shaped configuration, around the respective potential minimum of 6.437 a0, the interaction-induced mean first hyperpolarizability varies for 5?<? R/a0?<?10 as
Again, at the CCSD level, but for the L-shaped configuration around the respective potential minimum of 6.572 a0, this property varies for 5?<? R/a0?<?10 as
Graphical Abstract Interaction-induced mean dipole polarizability (\( \overline{a} \)) for the T-shaped configuration of H2–Ne calculated at the SCF, MP2, and CCSD levels of theory
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.
In this article, we explore, both theoretically and experimentally, the general reactivity of alkyl hydrogeno-phenylphosphinates with alcohols. We show that alcohol molecules act exclusively as nucleophilic species, and add to alkyl hydrogeno-phenylphosphinates, leading to pentacoordinated intermediates. These intermediates are shown to subsequently competitively undergo alcohol eliminations and/or Berry pseudorotations. This offers several possible routes for racemizations and/or alcohol exchange reactions. Transition standard Gibbs free energies predicted from DFT calculations for the overall alcohol exchange mechanism are shown to be compatible with those experimentally measured in case ethanol reacts with ethyl hydrogeno-phenylphosphinate (134.5~136.0 kJ mol?1 at 78 °C).
In this work we studied the structural and electronic properties of the metal–Schiff base complexes Ni\( {\mathrm{L}}_2^2 \) (1), Pd\( {\mathrm{L}}_2^1 \) (2), Zn\( {\mathrm{L}}_2^2 \) (3), and Ni\( {\mathrm{L}}_2^1 \)(4), where L1 and L2 are Schiff bases synthesized from salicylaldehyde and 2-hydroxy-5-methylbenzaldehyde, respectively. Natural bond analysis showed that in complexes 1 and 2, the metal ion coordinates to the ligands through electron donation from lone pairs on ligand nitrogen and oxygen atoms to s and d orbitals on the metal ion. In complex 3, metal–N and metal–O bonds are formed through charge transfer from the lone pairs on nitrogen and oxygen atoms to an s orbital of Zn. Dimethylation of the phenolate rings in the ligands decreases the energy gap and redshifts the spectrum of the nickel complex. The main absorptions observed were assigned on the basis of singlet-state transitions. The simulated spectra of the two complexes 1 and 2 are characterized by excited states with ligand-to-ligand charge-transfer (LLCT), metal-to-ligand charge-transfer (MLCT), ligand-to-metal charge-transfer (LMCT), and metal-centered (MC) character.
Bioethanol is one of the world’s most extensively produced biofuels. However, it is difficult to purify due to the formation of the ethanol–water azeotrope. Knowledge of the azeotrope structure at the molecular level can help to improve existing purification methods. In order to achieve a better understanding of this azeotrope structure, the characterization of (ethanol)5–water heterohexamers was carried out by analyzing the results of electronic structure calculations performed at the B3LYP/6-31+G(d) level. Hexamerization energies were found to range between ?36.8 and ?25.8 kcal/mol. Topological analysis of the electron density confirmed the existence of primary (OH…O) hydrogen bonds (HBs), secondary (CH…O) HBs, and H…H interactions in these clusters. Comparison with three different solvated alcohol systems featuring the same types of atom–atom interactions permitted the following order of stability to be determined: (methanol)5–water > (methanol)6 > (ethanol)5–water > (ethanol)6. These findings, together with accompanying geometric and spectroscopic analyses, show that similar cooperative effects exist among the primary HBs for structures with the same arrangement of primary HBs, regardless of the nature of the molecules involved. This result provides an indication that the molecular ratio can be considered to determine the unusual behavior of the ethanol–water system. The investigation also highlights the presence of several types of weak interaction in addition to primary HBs.
Graphical Abstract Water-ethanol clusters exhibit a variety of interaction types between their atoms, such as primary OH...O (blue), secondary CH...O (green) and H...H (yellow) interactions as revealed by Quantum Chemical Topology
The chelating properties of diethanoldithiocarbamate (DEDC) and π-electron flow from the nitrogen atom to the sulfur atom via a plane-delocalized π-orbital system (quasi ring) was studied using a density functional theory method. The molecular structure of DEDC and its complexes with Zn(II), Cd(II), and Hg(II) were also considered. First, the geometries of this ligand and DEDC-Zn(II), DEDC-Cd(II), and DEDC-Hg(II) were optimized, and the formation energies of these complexes were then calculated based on the electronic energy, or sum of electronic energies, with the zero point energy of each species. Formation energies indicated the DEDC-Zn(II) complex as the most stable complex, and DEDC-Cd(II) as the least stable. Structural data showed that the N1–C2 π-bond was localized in the complexes rather than the ligand, and a delocalized π-bond over S7–C2–S8 was also present. The stability of DEDC-Zn(II), DEDC-Cd(II), and DEDC-Hg(II) complexes increased in the presence of the non-specific effects of the solvent (PCM model), and their relative stability did not change. There was π-electron flow or resonance along N1–C2–S7 and along S7–C2–S8 in the ligand. The π-electron flow or resonance along N1–C2–S7 was abolished when the metal interacted with sulfur atoms. Energy belonging to van der Waals interactions and non-covalent delocalization effects between the metal and sulfur atoms of the ligand was calculated for each complex. The results of nucleus-independent chemical shift (NICS) indicated a decreasing trend as Zn(II)?<?Cd(II)?<?Hg(II) for the aromaticity of the quasi-rings. Finally, by ignoring van der Waals interactions and non-covalent delocalization effects between the metal and sulfur atoms of the ligand, the relative stability of the complexes was changed as follows:
In respective water or ethanol polarizable continuum cavity environments, simultaneous aldol condensation was performed using density functional theory (DFT) computational method to model the synthesis of optically active (RS)-1,2,4-butanetriol trinitrate (BTTN). The results of reaction energy barrier analysis suggested feasible routes with lower activation energies to obtain either the (R)- or (S)-configuration product in ethanolic solution. In addition, local analysis of average inter-particulate distances of reaction species revealed that a stronger inter-particulate interaction accompanied a shorter average distance in the ethanol system. The stabilization effect also indicated that related syntheses would be able to proceed in ethanol. Furthermore, relative to the production of (R)-BTTN, a lower overall energy of 425.3 kJ/mol was required for the synthesis of (S)-BTTN. Through analysis of the effects of temperature on the reaction rates of individual parallel stages of (R)- and (S)-species synthesis, it was simple to adjust the reaction temperature accordingly to differentiate between relative rates in order to obtain a product of a specific configuration.
The evolution of structural properties, thermodynamics and averaged (dynamic) total hardness values as a function of the composition of binary water–organic solvents, was rationalized in view of the intermolecular interactions. The organic solvents considered were ethanol, acetonitrile, and isopropanol at 0.25, 0.5, 0.75, and 1 mass fractions, and the results were obtained using molecular dynamics simulations. The site-to-site radial distribution functions reveal a well-defined peak for the first coordination shell in all solvents. A characteristic peak of the second coordination shell exists in aqueous mixtures of acetonitrile, whereas in the water–alcohol solvents, a second peak develops with the increase in alcohol content. From the computed coordination numbers, averaged hydrogen bonds and their lifetimes, we found that water mixed with acetonitrile largely preserves its structural features and promotes the acetonitrile structuring. Both the water and alcohol structures in their mixtures are disturbed and form hydrogen bonds between molecules of different kinds. The dynamic hardness values are obtained as the average over the total hardness values of 1200 snapshots per solvent type, extracted from the equilibrium dynamics. The dynamic hardness profile has a non-linear evolution with the liquid compositions, similarly to the thermodynamic properties of these non-ideal solvents.
Graphical abstract Computed dynamic total hardness, as a function of the cosolvent mass fraction for water–ethanol (EtOH), water–isopropanol (2PrOH) and water–acetonitrile (AN)
The diffusion coefficients of 14 n-alkanes (ranging from methane to n-tetradecane) in liquid and supercritical methanol at infinite dilution (at a pressure of 10.5 MPa and at temperatures of 299 K and 515 K) were deduced via molecular dynamics simulations. Values for the radial distribution function, coordination number, and number of hydrogen bonds were then calculated to explore the local structure of each fluid. The flexibility of the n-alkane (as characterized by the computed dihedral distribution, end-to-end distance, and radius of gyration) was found to be a major influence and hydrogen bonding to be a minor influence on the local structure. Hydrogen bonding reduces the flexibility of the n-alkane, whereas increasing the temperature enhances its flexibility, with temperature having a greater effect than hydrogen bonding on flexibility.
Graphical abstract The flexibility of the alkane is a major influence and the hydrogen bonding is a minor influence on the first solvation shell; the coordination numbers of long-chain n-alkanes in the first solvation shell are rather low
Protonation in the two-electron/two-proton reduction processes of 2,6-dichlorophenolindophenolate (DCIP) is investigated combining density functional theory (DFT) and molecular dynamics (MD) methods. DCIP (anion), DCIP?– (radical anion), and DCIP2? (dianion) are considered, including the electronic structure analysis from the prospective of quantum theory of atoms and molecules (QTAIM). It is shown that oxygen on the indophenolate moiety and nitrogen are the first and/or the second proton acceptor sites and their energetic order depends on the total charge of the system. MD simulations of differently charged species interacting with the solvent molecules have been performed for methanol, water, and oxonium cation (H3O+). Methanol and water molecules are found to form only hydrogen bonds with the solute irrespective of its charge. The calculated pKa values show that the imino group of DCIPH? is a weaker acid than water. While in the case of DCIP (and DCIP?–) plus oxonium cation, proton transfer from the solvent to the solute was evidenced for both aforementioned acceptor sites. In addition, MD simulations of bulks containing 15 and 43 molecules of water around the DCIP molecule have been performed, revealing the formation of 2–4 hydrogen bonds.
Graphical Abstract 2,6-Dichlorophenolindophenolate interacts with solvent molecules (water, oxonium cation and methanol). Hydrogen transfer and electronic structure are studied by DFT and molecular dynamics methods
We use the TDDFT/PBE0/6-31+G* method to determine the electronic absorption and emission energies, in different media, of the four forms of tocopherol, which differ by the number and the position of methyl groups on the chromanol. Geometries of the ground state S0 and the first singlet excited state S1 were optimized in the gas phase, and various solvents. The solvent effect is evaluated using an implicit solvation model (IEF-PCM). Our results are compared to the experimental ones obtained for the vitamin E content in several vegetable oils. For all forms of tocopherols, the HOMO–LUMO first vertical excitation is a π–π* transition. Gas phase and non-polar solvents (benzene and toluene) give higher absorption wavelengths than polar solvents (acetone, ethanol, methanol, DMSO, and water); this can be interpreted by a coplanarity between the O-H group and the chroman, allowing a better electronic resonance of the oxygen lone pairs and the aromatic ring, and therefore giving an important absorption wavelength, whereas the polar solvents give high emission wavelengths comparatively to gas phase and non-polar solvents. Fluorescence spectra permit the determination, the separation, and the identification of the four forms of tocopherols by a large difference in emission wavelength values.
Antiamoebin I (AAM-I) and zervamicin II (Zrv-IIB) are peptaibols that exert antibiotic activity through the insertion/disruption of cell membranes. In this study, we investigated how the folding of these peptaibols are affected when some of their native residues are replaced with proline analogues and asymmetrical D-α,α-dialkyl glycines (two classes of noncanonical amino acids). Systematic substitutions of native Aib, Pro, Hyp, and Iva residues were performed to elucidate the folding properties of the modified peptaibols incorporating noncanonical residues. The secondary structure of a peptaibol influences its ability to incorporate into membranes and therefore its function. Our findings reveal that native Zrv-IIB unfolds considerably in water. The presence of Iva and the noncanonical proline analogue cis-3-amino-L-proline (ALP) in both peptaibols induces helical structures. Inserting asymmetric glycines such as α-methyl-D-leucine (MDL) and α-methyl-D-phenylalanine (MDP) into the peptaibols induces folding. This preorganization in water may help to overcome the energy barrier required for peptaibol insertion into the membrane, as well as to facilitate the formation of transmembrane channels.
Graphical abstract AAM-I and Zrv-IIB peptidomimetics carrying MDL and ALP noncanonical amino acids, exhibiting improved helical secondary structure in water
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
Coarse-grained dynamical simulations have been performed to investigate the behavior of a surfactant micelle in the presence of six different alcohols: hexanol, octanol, decanol, dodecanol, tetradecanol, and hexadecanol. The self-assembly of sodium dodecyl sulfate (SDS) is modified by the alcohol molecules into cylindrical and bilayer micelles as a function of the alcohol/SDS mass ratio. Therefore, in order to understand, from a molecular point of view, how SDS and alcohol molecules self-organize to form the new micelles, different studies were carried out. Analysis of micelle structures, density profiles, and parameters of order were conducted to characterize the shape and size of those micelles. The density profiles revealed that the alcohol molecules were located at the water–micelle interface next to the SDS molecules at low alcohol/SDS mass ratio. At high alcohol/SDS mass ratios, alcohol molecules moved to the middle of the micelle by increasing their size and by producing a structural change. Moreover, micelle structures and sizes were influenced not only by the alcohol/SDS mass ratio but also by the order of the SDS and alcohol tails. Finally, the size of the micelles and enthalpy calculations were used as order parameters to determine a structural phase diagram of alcohol/SDS mixtures in water.
The human respiratory syncytial virus (hRSV) is the major cause of lower respiratory tract infection in children and elderly people worldwide. Its genome encodes 11 proteins including SH protein, whose functions are not well known. Studies show that SH protein increases RSV virulence degree and permeability to small compounds, suggesting it is involved in the formation of ion channels. The knowledge of SH structure and function is fundamental for a better understanding of its infection mechanism. The aim of this study was to model, characterize, and analyze the structural behavior of SH protein in the phospholipids bilayer environment. Molecular modeling of SH pentameric structure was performed, followed by traditional molecular dynamics (MD) simulations of the protein immersed in the lipid bilayer. Molecular dynamics with excited normal modes (MDeNM) was applied in the resulting system in order to investigate long time scale pore dynamics. MD simulations support that SH protein is stable in its pentameric form. Simulations also showed the presence of water molecules within the bilayer by density distribution, thus confirming that SH protein is a viroporin. This water transport was also observed in MDeNM studies with histidine residues of five chains (His22 and His51), playing a key role in pore permeability. The combination of traditional MD and MDeNM was a very efficient protocol to investigate functional conformational changes of transmembrane proteins that act as molecular channels. This protocol can support future investigations of drug candidates by acting on SH protein to inhibit viral infection.
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.
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.
Catalytic fields illustrate topology of the optimal charge distribution of a molecular environment reducing the activation energy for any process involving barrier crossing, like chemical reaction, bond rotation etc. Until now, this technique has been successfully applied to predict catalytic effects resulting from intermolecular interactions with individual water molecules constituting the first hydration shell, aminoacid mutations in enzymes or Si→Al substitutions in zeolites. In this contribution, hydrogen to fluorine (H→F) substitution effects for two model reactions have been examined indicating qualitative applicability of the catalytic field concept in the case of systems involving intramolecular interactions.
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.