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.
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.
We present a theoretical study on the detailed mechanism and kinetics of the H+HCN →H+HNC process. The potential energy surface was calculated at the complete basis set quantum chemical method, CBS-QB3. The vibrational frequencies and geometries for four isomers (H2CN, cis-HCNH, trans-HCNH, CNH2), and seven saddle points (TSn where n = 1 ? 7) are very important and must be considered during the process of formation of the HNC in the reaction were calculated at the B3LYP/6-311G(2d,d,p) level, within CBS-QB3 method. Three different pathways (PW1, PW2, and PW3) were analyzed and the results from the potential energy surface calculations were used to solve the master equation. The results were employed to calculate the thermal rate constant and pathways branching ratio of the title reaction over the temperature range of 300 up to 3000 K. The rate constants for reaction H + HCN → H + HNC were fitted by the modified Arrhenius expressions. Our calculations indicate that the formation of the HNC preferentially occurs via formation of cis–HCNH, the fitted expression is kPW2(T) = 9.98 × 10?22T2.41 exp(?7.62 kcal.mol?1/RT) while the predicted overall rate constant kOverall(T) = 9.45 × 10?21T2.15 exp(?8.56 kcal.mol?1/RT) in cm3 molecule?1s?1.
Graphical Abstract (a) Potential energy surface, (b) thermal rate constants as a function of temperature and (c) the branching ratios (%) of PW1, PW2, PW3 pathways involved in rm H + HCN → H + HNC process.
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).
Density functional theory (DFT) was used to study the stereochemistry, thermodynamic stability, and spectra of recently synthesized bis-Betti base derivatives of 2,3-dihydroxynaphthalene obtained through multicomponent reactions of 2,3-dihydroxynaphthalene with aminoisoxazole and benzaldehyde derivatives. The stereochemistry of the products was investigated by theoretically calculating the infrared (IR) and proton nuclear magnetic resonance (1H NMR) spectra of the diastereomers and comparing them to the corresponding experimental data. The thermochemical properties of the reactions, including the enthalpy, internal energy, entropy, and Gibbs free energy, were also calculated. The diastereoselectivity of the reactions was estimated from the equilibrium distribution of diastereomers. According to the results, the synthesis of bis-Betti bases is exothermic and accompanied by a decrease in entropy. The energy difference between the diastereomers is quite small, but the Gibbs free energy change for the equilibrium syn\( \rightleftharpoons \)anti favors the anti over syn configuration. These results are in good agreement with experimental observations.
Spin–spin coupling constants in water monomer and dimer have been calculated using several wave function and density functional-based methods. CCSD, MCSCF, and SOPPA wave functions methods yield similar results, specially when an additive approach is used with the MCSCF. Several functionals have been used to analyze their performance with the Jacob’s ladder and a set of functionals with different HF exchange were tested. Functionals with large HF exchange appropriately predict 1JOH, 2JHH and 2hJOO couplings, while 1hJOH is better calculated with functionals that include a reduced fraction of HF exchange. Accurate functionals for 1JOH and 2JHH have been tested in a tetramer water model. The hydrogen bond effects on these intramolecular couplings are additive when they are calculated by SOPPA(CCSD) wave function and DFT methods.
The first cytogenetic comparisons of five snapper species from Thailand were presented here. Renal cell samples were taken from blacktail snapper (Lutjanus fulvus), five lined snapper (L. quinquelineatus), dory snapper (L. fulviflamma), brownstripe red snapper (L. vitta), and mangrove red snapper (L. argentimaculatus). The mitotic chromosome preparation was prepared directly from kidney cells. Conventional staining and Ag-NOR banding techniques were applied to stain the chromosomes. The results exhibited that all five snapper species have the diploid chromosome numbers of 2n = 48 and the fundamental numbers (NF) of 48. The presences of large, medium, and small telocentric chromosomes were 22-24-2, 24-20-4, 36-10-2, 28-16-4 and 36-10-2, respectively. The Ag- NORs banding technique provides the pair of nucleolar organizer regions (NORs) at subcentromeric region of the long arm of the respective telocentric chromosome pairs 9, 1, 3, 4 and 9. Their karyotype formulas is as follows: L. fulvus (2n = 48): L22t + M24t + S2t, L. quinquelineatus (2n = 48): L24t + M20t + S4t, L. fulviflamma (2n = 48): Lt36 + Mt10 + St2, L. vitta (2n = 48): L28t + M16t + S4t, and L. argentimaculatus (2n = 48): L36t + M10t + S2t. 相似文献
Based on dissipative particle dynamics (DPD) methods and experimental data, we used an empirical relationship between the DPD temperature and the real temperature to build a model that describes the viscosity of molten TNT fluids. The errors in the predicted viscosity based on this model were no more than 2.3 %. We also studied the steady-state shear rheological behavior of molten TNT fluids containing nanoparticles (“nanofluids”). The dependence of the nanofluid viscosity on the temperature was found to satisfy an Arrhenius-type equation, η?=?AeB/T, where B, the flow activation energy, depends on particle content, size, and shape. We modified the Einstein-type viscosity model to account for the effects of nanoparticle solvation in TNT nanofluids. The resulting model was able to correctly predict the viscosities of suspensions containing nano- to microsized particles, and did not require any changes to the physical background of Einstein’s viscosity theory.
A series of conjugated multi-structured fluorescent probe molecules based on a salen ligand were designed and investigated in dimethyl sulfoxide solvent using a quantum-chemical method. The results indicate that the one-photon absorption and fluorescence emission spectra (λO and λEM) of these molecules generally show redshifts (of 23.1–74.5 and 22.7–116.6 nm, respectively) upon the coordination of the molecules to Zn2+. Large Stokes shifts (1511.2–11744.1 cm?1) were found for the molecules, meaning that interference between λO and λEM can be avoided for these molecules. The two-photon absorption spectra of the molecules usually present blueshifts, but the two-photon absorption cross-section (δ) greatly increases (by 221.5–868.0 GM) upon the coordination of the molecules with Zn2+. Most of the molecules show strong two-photon absorption peaks in the range 678.2–824.4 nm, i.e., in the near-infrared region. In a word, the expanded π-conjugated frameworks of these molecules lead to redshifted λO and λEM and enhanced δ values. Moreover, (L-phenyl)?2 and (L-phenyl-ethynyl)2 are the most suitable of the multi-structured molecules examined in this work for use as two-photon fluorescent probes for zinc ion detection in vivo.
Graphical Abstract Scheme of the calculated transition energies (E0k and E0n) and the transition dipole moments (M0k and Mkn). NTO 109, NTO 197 and NTO 228 of Zn(L-phenyl-ethynyl), Zn2(L-phenyl-ethynyl)2 and Zn3(L-phenyl)3 for one-photon absorption, respectively.
Designing and synthesizing novel electron-donor polymers with the high photovoltaic performances has remained a major challenge and hot issue in organic electronics. In this work, the exciton-dissociation (kdis) and charge-recombination (krec) rates for the PC61BM-PTDPPSe system as a promising polymer-based solar cell candidate have been theoretically investigated by means of density functional theory (DFT) calculations coupled with the non-adiabatic Marcus charge transfer model. Moreover, a series of regression analysis has been carried out to explore the rational structure–property relationship. Results reveal that the PC61BM-PTDPPSe system possesses the large open-circuit voltage (0.77 V), middle-sized exiton binding energy (0.457 eV), and relatively small reorganization energies in exciton-dissociation (0.273 eV) and charge-recombination (0.530 eV) processes. With the Marcus model, the kdis, krec, and the radiative decay rate (ks), are estimated to be 3.167×1011 s?1, 3.767×1010 s?1, and 7.930×108 s?1 respectively in the PC61BM-PTDPPSe interface. Comparably, the kdis is as 1～3 orders of magnitude larger than the krec and the ks, which indicates a fast and efficient photoinduced exciton-dissociation process in the PC61BM-PTDPPSe interface.
Graphical Abstract PTDPPSe is predicted to be a promising electron donor polymer, and the PC61BM-PTDPPSe system is worthy of further device research by experiments.
The mechanism of phenanthridines synthesis from the nitrogenation of 2-acetylbiphenyls (1) by TMSN3 in TFA has been studied by DFT calculations. Results at the B3LYP/6-311G(d) level showed that: 1) reaction of TMSN3/HN3 with the protonated form of 1 (1H+), which generates the key intermediate Cx+ by removal of TMSOH/H2O, is the rate determining step, and TMSN3 as the nitrogen source is certainly preferred over HN3. 2) from Cx+, the two pathways leading to 2xH+ and 3xH+ are both thermodynamically and kinetically feasible and competitive to each other. 3) The high barriers of the reverse reactions suggest that the ratio of the final products 2x:3x is determined by the branching ratio of reaction rates of Cx+ to intermediates Dx+ in pass_I and Ex+ in pass_II.
Graphical Abstract DFT results indicate that the replacement of -OH by -N3 which generates Cx+ controls the consumption rate of 1x H+, and the ratio of Cx+ transforms to Dx+ and Cx+ transforms to Ex+ (k:k') determines the final ratio of products 2x:3x.
The present study reports the geometries, electronic structures, growth behavior, and stabilities of neutral and ionized copper-doped germanium clusters containing 1–20 Ge atoms within the framework of linear combination of atomic orbitals density functional theory (DFT) under the spin-polarized generalized gradient approximation. It was found that Cu-capped Gen (or Cu-substituted Gen+1) and Cu-encapsulated Gen clusters mostly occur in the ground state at a particular cluster size (n). In order to explain the relative stabilities of the ground-state clusters, parameters such as the average binding energy per atom (BE), the embedding energy (EE), and the fragmentation energy (FE) of the clusters were calculated, and the resulting values are discussed. To explain the chemical stabilities of the clusters, parameters such as the energy gap between the highest occupied and the lowest unoccupied molecular orbitals (the HOMO–LUMO gap), the ionization energy (IP), the electron affinity (EA), the chemical potential (μ), the chemical hardness (η), and the polarizability were calculated, and the resulting values are also discussed. Natural atomic orbital (NAO) and natural bond orbital (NBO) analyses were also used to determine the electron-counting rule that should be applied to the most stable Ge10Cu cluster. Finally, the relevance of the calculated results to the design of Ge-based superatoms is discussed.
Figure Contributions of the valance orbitals of the Ge and Cu atom(s) to the HOMO of the ground-state icosahedral Ge10Cu cluster obtained from NBO analysis. The numbers below the clusters represent the occupancies of the HOMO orbitals
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.
Density functional theory (B3LYP, B3LYP-D2 and wB97XD functionals) was used in finite models of zigzag carbon nanotubes (CNT), (n,0)×k with n?=?6–9 and k?=?2–4, to systematically investigate the effects of size on their structural and electronic properties. We found that the ratio between the length (Lt) and the diameter (dt) of the pristine CNT has to be larger than 2, i.e., Lt/dt?>?2, in order to provide the observed experimental trends of C=C bond distances, as well as to maintain the atomic charges nearly constant and zero around the center of the tube. Therefore, the concepts of useful length and volume were developed and tested for the encapsulation process of HCN and C2H2 into CNTs. The energies involved in these processes, as well as the changes in molecular structure and electronic properties of the dopants and the CNTs are discussed and rationalized by the amount of charge transferred between dopant and CNT.
Density functional theory (DFT) was used to study the cobalt(I)-catalyzed enantioselective intramolecular hydroacylation of ketones and alkenes. All intermediates and transition states were fully optimized at the M06/6-31G(d,p) level (LANL2DZ(f) for Co). The results demonstrated that the ketone and alkene present different reactivities in the enantioselective hydroacylation. In ketone hydroacylation catalyzed by the cobalt(I)–(R,R)-Ph-BPE complex, reaction channel “a” to (R)-phthalide was more favorable than channel “b” to (S)-phthalide. Hydrogen migration was both the rate-determining and chirality-limiting step, and this step was endothermic. In alkene hydroacylation catalyzed by the cobalt(I)–(R,R)-BDPP complex, reaction channel “c” leading to the formation of (S)-indanone was the most favorable, both thermodynamically and kinetically. Reductive elimination was the rate-determining step, but the chirality-limiting step was hydrogen migration, which occurred easily. The results also indicated that the alkene hydroacylation leading to (S)-indanone formation was more energetically favorable than the ketone hydroacylation that gave (R)-phthalide, both thermodynamically and kinetically.
Graphical abstract A DFT study demonstrated that the ketone and alkene in the cobalt(I)-catalyzed enantioselective intramolecular hydroacylation showed different reactivities
Molecular dynamics simulations were performed to investigate the separation of trihalomethanes (THMs) from water using boron nitride nanosheets (BNNSs). The studied systems included THM molecules and a functionalized BNNS membrane immersed in an aqueous solution. An external pressure was applied to the z axis of the systems. Two functionalized BNNSs with large fluorinated-hydrogenated pore (F-H-pores) and small hydrogen-hydroxyl pore (H-OH-pores) were used. The pores of the BNNS membrane were obtained by passivating each nitrogen and boron atoms at the pore edges with fluorine and hydrogen atoms in the large pore or with hydroxyl and hydrogen atoms in the small pore. The results show that the BNNS with a small functionalized pore was impermeable to THM molecules, in contrast to the BNNS with a large functionalized pore. Using these membranes, water contaminants can be removed at lower cost.
Graphical Abstract A snapshot of the simulation system. The BNNS membrane with the large functionalized pore is located in the middle of the box. The size of the box is 3 × 3 × 5 nm3. Green chlorine, cyan carbon, red oxygen, white hydrogen
A grid-based variational method was proposed and applied to the ground state energies of atoms from the first to the third period of the periodic table. The nonuniform grid in the radial coordinate was defined by a q-exponential sequence. Some unusual properties between the optimum q-parameters and the electronic energies or atomic numbers are described. The behavior of the electronic energy, with respect to the q-parameter, yields near Hartree-Fock accuracy with a relatively small number of integration points. A simple relationship between the optimum q-parameters and the atomic numbers was found, which allowed the determination of the optimum q-parameters for atoms of the same period from two results. The remarkable results provide a simple alternative route to reach accurate results. The consistent results also suggest that this is not a random or accidental effect, but some optimum condition achieved by using a q-exponential mesh grid.
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
