The presented paper is focused on the calculation of hyperfine coupling constants (HFCC) of Cu2+ ion in water environment. To simulate the conditions of the electron paramagnetic resonance (EPR) experiment in aqueous phase, molecular dynamics using the density functional theory (DFT) was employed. In total three different functionals (BLYP, B3LYP, M06) were employed for studying their suitability in describing coordination of Cu2+ by water molecules. The system of our interest was composed of one Cu2+ cation surrounded by a selected number (between thirty and fifty) of water molecules. Besides the non-relativistic HFCCs (Fermi contact terms) of Cu2+ also the four-component relativistic HFCC calculations are presented. The importance of the proper evaluation of HFCCs, the inclusion of spin-orbit term, for Cu2+ containing systems (Neese, J. Chem. Phys. 118, 3939 2003; Almeida et al., Chem. Phys. 332, 176 2007) is confirmed at the relativistic four-component level of theory.
Graphical Abstract Five and six coordinated copper dication is solvated by adding extra water molecules to simulate conditions in aqueous solution. Molecular dynamics study is performed and nonrelativistic and relativistic hyperfine coupling constants are calculated subsequently.
Ionic liquids (ILs) constitute a fast growing class of compounds finding multiple applications in science and technology. Morpholinium-based ILs (MBILs) and their mixtures with polar molecular co-solvents are interesting as sustainable electrolyte systems for electrochemistry. We investigate local structures of protic and apropic morpholinium cations in acetonitrile (ACN) using semi-empirical molecular dynamics (MD) simulations. An impact of an anion (acetate) on the cation solvation regularities is discussed. Unlike oxygen, nitrogen of the morpholine ring is a strong electrophilic binding center. This site is responsible for the interactions of the cation with the solvent and with the anion. In protic MBILs, the role of nitrogen is delegated to the proton, which is linked to nitrogen. The acetate anion weakens solvation of the cation due to occupation of space near nitrogen or proton. The analysis reveals a favorable solvation of MBILs in ACN, which is a prerequisite for a new high-performance electrolyte system. The reported structural data were validated through point-to-point comparison with the MP2 post-Hartree-Fock theory and density functional theory.
This work investigated interactions between calcium cations (Ca2+) and three common types of oxygen-based functional groups of concrete superplasticizers using density functional theory (DFT) calculations and all-atom molecular dynamics (MD) simulations. The three common types of oxygen-based functional groups were modeled as three hypothetical, low-molecular-weight organic molecules, each containing a methyl-terminated oxyethylene dimer and an adsorbing head of two oxygen-based functional groups, and are referred to as carboxylate, sulfonate, and phosphate groups, respectively, following the usual terminology in the field of concrete admixtures. Our DFT results show that the binding strength of the three groups with calcium cations follows (from high to low) phosphate>carboxylate>sulfonate, and both the electrophilic attack and the chemical reactivity of the three groups contribute significantly to the binding strength. The MD simulation results indicate that the adsorption of the three small molecules on the calcite (1 0 4) surface in aqueous solution shares a similar pattern in the sense that just two oxygen atoms of two adjacent anchor groups adsorb on the calcium atoms on the top layer of the crystal. The adsorption strength among the three types of functional groups follows the same order as the binding strength obtained from DFT calculations; both results corroborate a similar rule-of-thumb established by experiments. Furthermore, interactions of the three types of groups with water molecules suggest that strong hydrogen-bonding interactions exist in those systems.
Efficient design of ionic compounds requires a systematic understanding of cation–anion interactions. Weakening of electrostatic attraction is essential to increase the liquid range of the ionic compound and decrease its melting point. Here, we report simulations of the closest-approach cation–anion distances in a variety of ion pairs containing the tetrakis(pentafluorophenyl)borate (TFPB–) anion. Small alkali cations (Li+, Na+) penetrate the TFPB– core, whereas K+ and larger organic cations do not. In the latter case, the shortest possible distance from the cations to the boron atom of TFPB– ranges from 0.50 nm to 0.63 nm. TFPB– was shown to be substantially rigid, providing a steric hindrance to thermodynamically efficient cation–anion coordination. Our results prove that TFPB– is more efficient for electrostatic charge confinement than the tetraoctylammonium cation, whereas the perfluorophenyl group is more efficient than linear alkyl chains. These simulations will motivate development of TFPB–-based ionic liquids with low phase transition points.
Recently, a series of xanthone analogues has been identified as α-glucosidase inhibitors. To provide deeper insight into the three-dimensional (3D) structural requirements for the activities of these molecules, CoMFA and CoMSIA approaches were employed on 54 xanthones to construct 3D-QSAR models. Their bioactive conformations were first investigated by docking studies and optimized by subsequent molecular dynamics (MD) simulations using the homology modeled structure of the target protein. Based on the docking/MD-determined conformers, 3D-QSAR studies generated several significant models in terms of 47 molecules as the training set. The best model (CoMSIA-SHA) yielded q2 of 0.713, r2 of 0.967 and F of 140.250. The robustness of the model was further externally confirmed by a test set of the remaining molecules (q2 = 0.793, r2 = 0.902, and k = 0.905). Contour maps provided much information for future design and optimization of new compounds with high inhibitory activities towards α-glucosidase.
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.
The number of hydrogen bonds and detailed information on the interlayer spacing of graphene oxide (GO) confined water molecules were calculated through experiments and molecular dynamics simulations. Experiments play a crucial role in the modeling strategy and verification of the simulation results. The binding of GO and water molecules is essentially controlled by hydrogen bond networks involving functional groups and water molecules confined in the GO layers. With the increase in the water content, the clusters of water molecules are more evident. The water molecules bounding to GO layers are transformed to a free state, making the removal of water molecules from the system difficult at low water contents. The diffuse behaviors of the water molecules are more evident at high water contents. With an increase in the water content, the functional groups are surrounded by fewer water molecules, and the distance between the functional groups and water molecules increases. As a result, the water molecules adsorbed into the GO interlamination will enlarge the interlayer spacing. The interlayer spacing is also affected by the number of GO layers. These results were confirmed by the calculations of number of hydrogen bonds, water state, mean square displacement, radial distribution function, and interlayer spacing of hydrated GO.
Graphical Abstract This work research the interaction between GO functional groups and confined water molecules. The state of water molecules and interlayer spacing of graphene oxide were proved to be related to the number of hydrogen bonds.
Development of new energetic salts is the key factor in replacing low performance compounds in conventional formulations of high explosives as well as propellants. Ten salts based on the nitroformate anion and various nitrogen-rich cations were designed and their geometric optimizations carried out using the density functional method. With reasonable oxygen balance (from ?36 % to 0 %), heats of formation (47–624 kJ mol?1) and high densities (1.81–1.89 g cm?3), the detonation velocity (D) and pressure (P) values of salts were calculated as 8.62–9.36 km s?1 and 33.10–40.01 GPa, respectively. Lastly, the nitroformate salts studied in this work are of prospective interest as high performance explosives.
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.
Complexes of the dipeptide phenylalanine–phenylalanine (Phe–Phe) with divalent metal cations (Cu2+, Zn2+, Ca2+ and Ba2+) were studied at the B3LYP and MP2 levels of theory with the basis sets 6-311++G(d,p) and 6-31 + G(d) in the gas phase. The relative energies of these complexes indicated that cation–π bidentate/tridentate conformations are more favourable than other conformations with uncoordinated rings. These findings were confirmed by the calculated values of thermodynamic parameters such as the Gibbs free energy. Natural bond orbital (NBO) analysis was carried out to explore the metal–ligand coordination in Phe–Phe–Cu2+/Zn2+ complexes. Possible orbital transitions, types of orbitals and their occupancies were determined for a range of Phe–Phe–Cu2+/Zn2+ complexes. The charge transfer involved in various orbital transitions was explored by considering the second-order perturbation energy. NBO analysis revealed that the change transfer is stronger when the metal cation uses both the 4s + 4p subshells rather than just its 4p subshell. We also performed molecular dynamics (MD) simulations to check the stability and consistency of the most favourable binding motifs of Cu2+, Zn2+, Ca2+ and Ba2+ with Phe–Phe over time. The structures of the Phe–Phe–Cu2+/Zn2+/Ca2+/Ba2+ complexes obtained using MD simulation were found to be in good agreement with those obtained in the DFT-based calculations.
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).
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.
A mechanistic investigation using Becke3LYP density functional theory (DFT) was carried out on the palladium-catalyzed amidition of bromobenzene and tBu-isocyanide. The whole catalytic cycle consists of five steps: oxidative addition, migratory insertion, anion exchange, reductive elimination, and hydrogen migration. The rate-determining step is oxidative addition, with a small Gibbs free energy of 14.6 kcal mol?1. In the migratory insertion step, tBu-isocyanide provides an important source of carboxy and amino groups to establish the amide group. For anion exchange, path 1a is suggested as the most favorable pathway with the help of the base, and water provides a source of oxygen which is perfectly in line with experimental observations. Finally, in the hydrogen migration step, we illustrate that the six-membered ring path is energetically favored due to the assisting influence of water. In addition, our calculations indicate that using dimethyl sulfoxide as a solvent does not change the rate-determining step.
The interactions of the drugs amlodipine and paroxetine, which are prescribed respectively for treatment of hypertension and depression, with the metabolizing enzyme cytochrome CYP2B4 as the drug target, have been studied by molecular dynamics (MD) simulation. Poly ethylene glycol was used to control the drugs’ interactions with each other and with the target CYP2B4. Thirteen simulation systems were carefully designed, and the results obtained from MD simulations indicated that amlodipine in the PEGylated form prescribed with paroxetine in the nonPEGylated form promotes higher cytochrome stability and causes fewer fluctuations as the drugs approach the target CYP2B4 and interact with it. The simulation results led us to hypothesize that the combination of the drugs with a specific drug ratio, as proposed in this work, manifests more effective diffusivity and less instability while metabolizing with enzyme CYP2B4. Also, the active residues in the CYP2B4 enzyme that interact with the drugs were determined by MD simulation, which were consistent with the reported experimental results.
This work is focused on the donor properties of cobalt-exchanged cationic sites in zeolites. It is based on cluster and periodic density functional theory modeling for relevant {[Co(II)(NH3)n]–NO} adducts, where Co(II) means a cobalt cation embedded either in a periodic model of chabasite (CHA) zeolite or in model clusters. NO stretching frequencies were derived from MD trajectories and compared to harmonic values from cluster calculations. By relating calculated NO frequencies to experimental FTIR spectra, it was shown that the forms of {Co(II)-NO} adducts comprising three or four ammonia co-ligands dominate the spectrum taken in ammonia-saturation conditions while forms with two NH3 ligands prevail under intermediate ammonia saturation. Finally, this work confirms the critical dependence of Co(II) activation ability towards NO upon the center donor properties, reinforced by ligation of strong donor ammonia ligands. However, strongly bound ligands appear also to compete with interaction of the center with the electron-rich framework, and a balance must be observed to maintain optimal activation ability.
Graphical abstract A snapshot from MD trajectory showing a fragment of periodic framework with twoCo(II)–NO centers, bound to one framework oxygen and strongly coordinating three ammonia ligands with four others forming the second coordination sphere
Beryllium telluride (BeTe) with cubic zinc-blende (ZB) structure was studied using ab initio constant pressure method under high pressure. The ab initio molecular dynamics (MD) approach for constant pressure was studied and it was found that the first order phase transition occurs from the ZB structure to the nickel arsenide (NiAs) structure. It has been shown that the MD simulation predicts the transition pressure PT more than the value obtained by the static enthalpy and experimental data. The structural pathway reveals MD simulation such as cubic → tetragonal → orthorhombic → monoclinic → orthorhombic → hexagonal, leading the ZB to NiAs phase. The phase transformation is accompanied by a 10% volume drop and at 80 GPa is likely to be around 35 GPa in the experiment. In the present study, our obtained values can be compared with the experimental and theoretical results.
In this work we used a combination of classical molecular dynamics and simulated annealing techniques to shed more light on the conformational flexibility of 12 adenosine triphosphate (ATP) analogues in a water environment. We present simulations in AMBER force field for ATP and 12 published analogues [Shah et al. (1997) Proc Natl Acad Sci USA 94: 3565–3570]. The calculations were carried out using the generalized Born (GB) solvation model in the presence of the cation Mg2+. The ion was placed at a close distance (2 Å) from the charged oxygen atoms of the beta and gamma phosphate groups of the ?3 negatively charged ATP analogue molecules. Analysis of the results revealed the distribution of inter-proton distances H8–H1′ and H8–H2′ versus the torsion angle ψ (C4–N9-C1′–O4′) for all conformations of ATP analogues. There are two gaps in the distribution of torsion angle ψ values: the first is between ?30 and 30 degrees and is described by cis-conformation; and the second is between 90 and 175 degrees, which mostly covers a region of anti conformation. Our results compare favorably with results obtained in experimental assays [Jiang and Mao (2002) Polyhedron 21:435–438].
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.
Human matrix metalloproteinase (MMP)-1 or collagenase–1 plays a significant role in embryonic development, tissue remodeling, and is also involved in several diseases like arthritis, metastasis, etc. Molecular dynamics simulation studies on hMMP-1 X-ray structures (PDB Id. 1CGE, 1CGF, 1CGL, 1HFC, and 2TCL) suggest that the three conserved water molecules (WH/1, WI, WS) are coordinated with catalytic zinc (ZnC), and one water molecule (W) is associated at structural zinc ion (ZnS). Transition of the coordination geometry around ZnC from tetrahedral to octahedral and tetrahedral to trigonal bipyramidal at ZnS are also observed during the dynamics. Recognition of two zinc ions through water mediated bridges (ZnC – WH (W1)…W2….H183 – ZnS) and stabilization of secondary coordination zone around the metal ions indicates the possibility of ZnC…ZnS coupled catalytic mechanism in hMMP-I. This study not only reveals a functionally important role of conserved water molecules in hMMP-I but also highlights the involvement of other non catalytic residues, such as S172 and D170 in the catalytic mechanism. The results obtained in this study could be relevant for importance of conserved water mediated recognition site of the sequence residue id. 202(RWTNNFREY)210, interaction of W(tryptophan)203 to zinc bound histidine, their influence on the water molecules that are involved in bridging between ZnC and ZnS, and structure-based design of specific hMMP inhibitors.
Multiple molecular dynamics simulations have been performed to explore the transport properties of single methane, methanol, and ethanol molecules through the water-filled transmembrane cyclic peptide nanotube (CPNT) of \( 8\times {\left(\mathrm{W}\underline {\mathrm{L}}\right)}_4-\mathrm{POPE} \), as well as the potential application of this CPNT in the separation of an alcohol/water mixture. Molecular size and hydrophilicity/hydrophobicity were found to significantly influence molecular diffusion behavior in the channel. Methane and ethanol display more explicit distributions in midplane regions, while methanol mainly occurs in α-plane zones. Methane and ethanol drift faster near an α-plane zone, whereas methanol diffuses uniformly throughout the whole transmembrane region. The dipole orientation of channel methanol is significantly affected by the bare carbonyl groups at the tube mouths and flips mainly in gap 4, whereas the rotation of ethanol is blocked. Ball-shaped hydrophobic methane experiences more flips in gap 4. The PMF (potential of mean force) profiles of the three organic molecules disclose their different diffusion behaviors in the CPNT. Amphiphilic alcohols are able to form direct H-bonds with channel water and the tube. Both single and double water bridges with the tube were observed in the methanol and ethanol systems. The different adsorption behaviors of the alcohols and water in the dehydrated CPNT may lead to the potential application of the CPNT as a means of separating alcohols from water.
