Exploration of various electronic properties along the reaction coordinate for hydration of Pt(II) and Ru(II) complexes; the CCSD, MPx, and DFT computational study

In the study behavior of molecular electrostatic potential, averaged local ionization energy, and reaction electronic flux along the reaction coordinate of hydration process of three representative Ru(II) and Pt(II) complexes were explored using both post-HF and DFT quantum chemical approximations. Previously determined reaction mechanisms were explored by more detailed insight into changes of electronic properties using ωB97XD functional and MP2 method with 6–311++G(2df,2pd) basis set and CCSD/6–31(+)G(d,p) approach. The dependences of all examined properties on reaction coordinate give more detailed understanding of the hydration process.

Theoretical study of the decomposition of ethyl and ethyl 3-phenyl glycidate

The mechanism of the decomposition of ethyl and ethyl 3-phenyl glycidate in gas phase was studied by density functional theory (DFT) and MP2 methods. A proposed mechanism for the reaction indicates that the ethyl side of the ester is eliminated as ethylene through a concerted six-membered cyclic transition state, and the unstable intermediate glycidic acid decarboxylates rapidly to give the corresponding aldehyde. Two possible pathways for glycidic acid decarboxylation were studied: one via a five-membered cyclic transition state, and the other via a four-membered cyclic transition state. The results of the calculations indicate that the decarboxylation reaction occurs via a mechanism with five-membered cyclic transition state.

DFT studies of the conversion of four mesylate esters during reaction with ammonia

The energetics of the Menshutkin-like reaction between four mesylate derivatives and ammonia have been computed using B3LYP functional with the 6-31+G** basis set. Additionally, MPW1K/6-31+G** level calculations were carried out to estimate activation barrier heights in the gas phase. Solvent effect corrections were computed using PCM/B3LYP/6-31+G** level. The conversion of the reactant complexes into ion pairs is accompanied by a strong energy decrease in the gas phase and in all solvents. The ion pairs are stabilized with two strong hydrogen bonds in the gas phase. The bifurcation at C2 causes a significant activation barrier increase. Also, bifurcation at C5 leads to noticeable barrier height differentiation. Both B3LYP/6-31+G** and MPW1K/6-31+G** activation barriers suggest the reaction 2 (2a?+?NH₃) to be the fastest in the gas phase. The reaction 4 is the slowest one in all environments.

Perspectives on the reaction force constant

A synchronous, concerted chemical process is rigorously divided by the reaction force F(R), the negative gradient of V(R), into “reactant” and “product” regions which are dominated by structural changes and an intervening “transition” region which is electronically intensive. The reaction force constant κ(R), the second derivative of V(R), is negative throughout the transition region, not just at the nominal transition state, at which κ(R) has a minimum. This is consistent with experimental evidence that there is a transition region, not simply a specific point. We show graphically that significant nonsynchronicity in the process is associated with the development of a maximum of κ(R) in the transition region, which increases as the process becomes more nonsynchronous. (We speculate that for a nonconcerted process this maximum is actually positive.) Thus, κ(R) can serve as an indicator of the level of nonsynchronicity.

Computationally designed prodrugs of statins based on Kirby’s enzyme model

DFT calculations at B3LYP/6-31G(d,p) for intramolecular proton transfer in Kirby’s enzyme models 1–7 demonstrated that the reaction rate is dependent on the distance between the two reacting centers, r_GM, and the hydrogen bonding angle, α, and the rate of the reaction is linearly correlated with r_GM and α. Based on these calculation results three simvastatin prodrugs were designed with the potential to provide simvastatin with higher bioavailability. For example, based on the calculated log EM for the three proposed prodrugs, the interconversion of simvastatin prodrug ProD 3 to simvastatin is predicted to be about 10 times faster than that of either simvastatin prodrug ProD 1 or simvastatin ProD 2. Hence, the rate by which the prodrug releases the statin drug can be determined according to the structural features of the promoiety (Kirby’s enzyme model).

Variation of the electronic dipole polarizability on the reaction path

The reaction force and the electronic flux, first proposed by Toro-Labbé et al. (J Phys Chem A 103:4398, 1999) have been expressed by the existing conceptual DFT apparatus. The critical points (extremes) of the chemical potential, global hardness and softness have been identified by means of the existing and computable energy derivatives: the Hellman-Feynman force, nuclear reactivity and nuclear stiffness. Specific role of atoms at the reaction center has been unveiled by indicating an alternative method of calculation of the reaction force and the reaction electronic flux. The electron dipole polarizability on the IRC has been analyzed for the model reaction HF + CO→HCOF. The electron polarizability determined on the IRC α _e (ξ) was found to be reasonably parallel to the global softness curve S(ξ). The softest state on the IRC (not TS) coincides with zero electronic flux.

Investigation of the binding network of IGF-I on the cavity surface of IGFBP4

Insulin-like growth factor-binding proteins (IGFBPs) control bioactivity and distribution of insulin-like growth factors (IGFs) through high-affinity complex of IGFBP and IGF. To get more insight into the binding interaction of IGF system, the site-directed mutagenesis and force-driving desorption methods were employed to study the interaction mechanism of IGFBP4 and IGF-I by molecular dynamics (MD) simulation. In IGF-I, residues Gly7 to Asp12 were found to be the hot spots and they mainly anchored on the N-domain of IGFBP4. The contact area, the shape and size of protein, the surroundings of the binding site, the hydrophobic and electrostatic interaction between the two proteins worked as a complex network to regulate the protein-protein interaction. It was also found that the unfolding of the helix was not inevitable in the mutant, and it could be regulated by careful selection of the substituted amino acid.

Can Satraplatin be hydrated before the reduction process occurs? The DFT computational study

Hydration reactions of two anticancer Pt(IV) complexes JM149 and JM216 (Satraplatin) were studied computationally together with the hydration of the Pt(II) complex JM118, which is a product of the Satraplatin reduction. Thermodynamic and kinetic parameters of the reactions were determined at the B3LYP/6-311++G(2df.2pd)//B3LYP/6-31 + G(d)) level of theory. The water solution was modeled using the COSMO implicit solvation model, with cavities constructed using Klamt’s atomic radii. It was found that hydration of the Pt(IV) complexes is an endergonic/endothermic reaction. It follows the (pseudo)associative mechanism is substantially slower (k?≈?10^-11 s^?1) than the corresponding reaction of Pt(II) analogues ((k?≈?10^-5 s^?1). Such a low value of the reaction constant signifies that the hydration of JM149 and Satraplatin is with high probability a kinetically forbidden reaction. Similarly to JM149 and Satraplatin, the hydration of JM118 is an endothermic/endoergic reaction. On the other hand, the kinetic parameters are similar to those of cisplatin Zimmermann et al. (J Mol Model 17:2385–2393, 2011), allowing the hydration reaction to occur at physiological conditions. These results suggest that in order to become active Satraplatin has to be first reduced to JM118, which may be subsequently hydrated to yield the active species.

Nickel/zinc-catalyzed decarbonylative addition of anhydrides to alkynes: A DFT study

Density functional theory (DFT) was used to investigate the nickel- or nickel(0)/zinc- catalyzed decarbonylative addition of phthalic anhydrides to alkynes. All intermediates and transition states were optimized completely at the B3LYP/6-31+G(d,p) level. Calculated results indicated that the decarbonylative addition of phthalic anhydrides to alkynes was exergonic, and the total free energy released was ?87.6 kJ mol^?1. In the five-coordinated complexes M4a and M4b, the insertion reaction of alkynes into the Ni-C bond occurred prior to that into the Ni-O bond. The nickel(0)/zinc-catalyzed decarbonylative addition was much more dominant than the nickel-catalyzed one in whole catalytic decarbonylative addition. The reaction channel CA→M1'→T1'→M2'→T2'→M3a'→M4a'→T3a1'→M5a1' →T4a1'→M6a'→P was the most favorable among all reaction pathways of the nickel- or nickel(0)/zinc- catalyzed decarbonylative addition of phthalic anhydrides to alkynes. And the alkyne insertion reaction was the rate-determining step for this channel. The additive ZnCl₂ had a significant effect, and it might change greatly the electron and geometry structures of those intermediates and transition states. On the whole, the solvent effect decreased the free energy barriers.

Zinc-, cadmium-, and mercury-containing one-dimensional tetraphenylporphyrin arrays: a DFT study

This work describes theoretical and experimental studies on glycerol esterification to obtain acetins focusing on the obtained isomers. The reaction of glycerol with acetic acid was carried out on Amberlyst 36 wet. Density functional theory calculations on the level of M06-2X functional and 6-311+G(d,p) basis set are carried out and the most stable structures of the reactants and products are located by considering a large number of conformers. The thermodynamics is discussed in terms of the calculated reaction Gibbs free energy. The AIM theory was used to characterize reactants and products. The glycerol esterification with acetic acid is found to be thermodynamically favored, with exothermal property. These agree well with experiments and allow us to explain the relative selectivity of products.

Rational design of the survivin/CDK4 complex by combining protein–protein docking and molecular dynamics simulations

Survivin, the smallest inhibitor of apoptosis protein (IAP), is a valid target for cancer research. It mediates both the apoptosis pathway and the cell cycle and has been proposed to form a complex with the cyclin-dependent kinase protein CDK4. The resulting complex transports CDK4 from the cytosol to the nucleus, where CDK4 participates in cell division. Survivin has been recognized as a node protein that interacts with several partners; disruption of the formed complexes can lead to new anticancer compounds. We propose a rational model of the survivin/CDK4 complex that fulfills the experimental evidence and that can be used for structure-based design of inhibitors modifying its interface recognition. In particular, the suggested complex involves the alpha helical domain of survivin and resembles the mode of binding of survivin in the survivin/borealin X-ray structure. The proposed model has been obtained by combining protein–protein docking, fractal-based shape complementarity, electrostatics studies and extensive molecular dynamics simulations.

Conformational entropy of a polymer chain grafted to rough surfaces

A polymer molecule (represented by a statistical chain) end-grafted to a topologically rough surface was studied by static MC simulations. A modified self-avoiding walk on a cubic lattice was used to model the polymer in an athermal solution. Different statistical models of surface roughness were applied. Conformational entropies of chains attached to uncorrelated Gaussian, Brownian, and fractional Brownian surfaces were calculated. Results were compared with the predictions of a simple analytical model of a macromolecule end-grafted to a fractal surface.

Insight into substituent effects in Cal-B catalyzed transesterification by combining experimental and theoretical approaches

Candida antarctica lipase B (Cal-B) is one of the most recognized biocatalysts because of its high degree of selectivity in a broad range of synthetic applications of industrial importance. Herein, the substituent effects involved in transesterification catalyzed by Cal-B are explored in detail using a combination of experimental analysis and theoretical modeling. The transesterification ability of Cal-B was experimentally determined with 22 vinyl ester analogs and ribavirin as substrates and, on this basis, a series of quantitative structure-activity relationship (QSAR) models are developed using various structural parameters characterizing the variation in substituent groups of the substrate molecules. The resulting models exhibit a good stability and predictive power, from which five most important properties are highlighted and engaged to ascertain the structural basis and reaction mechanism underlying the transesterification. From the modeling analysis it is seen that the size, geometry, and charge distributions of substrate exert a significant effect on reaction yield, where, the size of the substituent group was the most significant impact factor on the reaction yield, the charge distribution was the second, and then the topological structure of the substrate.

Theoretical evaluation of flotation performance of carboxyl hydroxamic acids with different number of polar groups on the surfaces of diaspore (010) and kaolinite (001)

The adsorption behaviors of three carboxyl hydroxamic acids on diaspore (010) and kaolinite (001) have been studied by density functional theory (DFT) and molecular dynamics (MD) method. The results indicated that carboxyl hydroxamic acids could adsorb on diaspore surface by ionic bonds and hydrogen bonds, and adsorb on kaolinite surface by hydrogen bonds. The models of carboxyl hydroxamic acids adsorbed on diaspore and kaolinite surfaces are proposed.

Heavy periodane

The potential energy surface of the hypothetical NaMgAlSiPSCl system (heavy periodane) is exhaustively analyzed via the gradient embedded genetic algorithm (GEGA) in combination with density functional theory (DFT) computations. The electronegativity differences among the elements in both the second and third rows of the periodic table indicate that low-energy heavy periodane structures are obtained when highly electronegative and electropositive elements are bound together, but the global minimum of the heavy periodane system is completely different to its second-row analog (LiBeBCNOF).

Charge sensitivity approach to mutual polarization of reactants: molecular mechanics perspective

Charge sensitivity analysis (CSA) in force-field atoms resolution was applied to describe the mutual polarization of reactants as well as charge-transfer (CT) effects. An inclusion complex of β-cyclodextrin with salicylic acid was used as a model system. Three CSA models were taken into account and verified on a Born–Oppenheimer molecular dynamics (BOMD) trajectory. The models differed in terms of the equilibrium conditions imposed on the system. It was demonstrated that mutual polarization is an important source of stabilization, in contrast to the results obtained from static charge calculations. The energy lowering induced by CT was small and comparable to the CT stabilization that occurs in hydrogen-bonded systems. All models correctly described the main topological features of the BOMD energy surface. CSA in force-field atoms resolution qualitatively reproduced the charge reorganization accompanying hydrogen-bond formation. It was shown that CSA parameters are very sensitive to the bond formation process, which suggests that they could be applied in reactive force fields as detectors of newly formed chemical bonds.

Prodrugs for masking bitter taste of antibacterial drugs—a computational approach

DFT calculations for the acid-catalyzed hydrolysis of several maleamic acid amide derivatives revealed that the reaction rate-limiting step is determined on the nature of the amine leaving group. Further, it was established that when the amine leaving group was a secondary amine, acyclovir or cefuroxime moiety the tetrahedral intermediate formation was the rate-limiting step such as in the cases of acyclovir ProD 1- ProD 4 and cefuroxime ProD 1- ProD 4. In addition, the linear correlation between the calculated and experimental rates provided a credible basis for designing prodrugs for masking bitter taste of the corresponding parental drugs which have the potential to release the parent drug in a sustained release fashion. For example, based on the DFT calculated rates the predicted t_1/2 (a time needed for 50 % of the reactant to be hydrolyzed to products) for cefuroxime prodrugs, cefuroxime ProD 1- ProD 4, were 12 min, 18 min, 200 min and 123 min, respectively.