Two new rotenoid glycosides from the rhizomes of Stemona curtisii Hook. f.

Two new rotenoid glycosides named stemonal 11-O-ß-D-glucopyranoside and 6-O-methylstemonal 11-O-ß-D-glucopyranoside together with ten known metabolites were isolated from the rhizomes of Stemona curtisii. The chemical structures of the new compounds were elucidated based on the analysis of their 1D and 2D NMR and HRESIMS, while the sugar unit and absolute configuration were determined by chemical hydrolysis and ECD analysis. Among the tested compounds for anti-α-glucosidase assay, stemonal showed an inhibitory effect (IC50 = 38.67 µM), which is 2.4-fold more potent than acarbose. Cytotoxic evaluation against the lung adenocarcinoma A549 cell line indicated that none of the compounds were strongly active to suppress the cancer cell growth at 100 µM. This work describes the occurrence of rotenoids bearing a sugar moiety, which are reported for the first time in the genus Stemona. The isolated compound's α-glucosidase inhibitory potential provides insight for further investigation of natural rotenoids as anti-diabetic agents.

DFT insight into metals and ligands substitution effects on reactivity of phenoxy-imine catalysts for ethylene polymerization.

The reaction mechanism of ethylene (ET) polymerization catalyzed by the phenoxy-imine (FI) ligands using DFT calculations was studied. Among five possible isomers, isomer A which has an octahedral geometry and a (cis-N/trans-O/cis-Cl) arrangement is the most stable pre-reaction Ti-FI dichloride complex. The isomer A can be activated by MAO to form the active catalyst and the active form was used for the study of the mechanism for Ti-FI. The second ethylene insertion was found to be the rate-determining step of the catalyzed ethylene polymerization. To examine the effect of group IVB transition metals (M = Ti, Zr, Hf) substitutions, calculated activation energies at the rate-determining step (EaRDS) were compared, where values of EaRDS of Zr < Hf < Ti agree with experiments. Moreover, we examined the effect of substitution on (O, X) ligands of the Ti-phenoxy-imine (Ti-1) based catalyst. The results revealed that EaRDS of (O, N) > (O, O) > (O, P) > O, S). Hence, the (O, S) ligand has the highest potential to improve the catalytic activity of the Ti-FI catalyst. We also found the activation energy to be related to the Ti-X distance. In addition, a novel Ni-based FI catalyst was investigated. The results indicated that the nickel (II) complex based on the phenoxy-imine (O, N) ligand in the square-planar geometry is more active than in the octahedral geometry. This work provides fundamental insights into the reaction mechanism of M - FI catalysts which can be used for the design and development of M - FI catalysts for ET polymerization.

Schomburginones A‒J, geranylated benzophenones from the leaves of Garcinia schomburgkiana and their cytotoxic and anti-inflammatory activities.

Ten undescribed benzophenones, schomburginones A-J, together with 14 known analogs were isolated from the leaves of Garcinia schomburgkiana, an edible plant native to the Indochina region. The structures of the undescribed compounds were elucidated by NMR combined with HRMS spectroscopy, while their absolute configurations were determined using ECD and single-crystal X-ray diffraction analysis. The isolated metabolites represent benzophenone derivatives containing a modified monoterpene unit, including tri- and tetracyclic skeletons, which are rarely found in genus Garcinia. The cytotoxic evaluation on three cancerous cell lines demonstrated that schomburginone G, schomburginone H, and 3-geranyl-2,4,6-trihydroxybenzophenone were active against HeLa cells with IC50 values in the range of 12.2-15.7 µM, respectively, and selective compared to the non-cancerous L929 cells (SI > 3.5). In addition, the three cytotoxic compounds together with clusiacyclol A showed significant NO inhibitory activity in RAW 264.7 macrophage cells over 85% inhibition without obvious cytotoxicity at a final concentration of 100 µM. The promising activities of these compounds in cytotoxic and anti-inflammatory assays make them attractive for further study in the development of anticancer drugs.

Garcowacinols A-J, cytotoxic polyprenylated benzoylphloroglucinol derivatives from the twigs of Garcinia cowa.

Ten undescribed polyprenylated benzoylphloroglucinol derivatives named garcowacinols A‒J (1-10) and four known analogues (11-14) were isolated from the twigs of Garcinia cowa. Their structures were determined by spectroscopic data analysis (1D and 2D NMR and HRESIMS), and their absolute configurations were established based on NOESY and ECD data. All isolated compounds were evaluated for their cytotoxicity against five types of human cancer cells (KB, HeLa S3, MCF-7, Hep G2, and HT-29) as well as Vero cells by MTT colorimetric assay. Garcowacinol C was significantly active against all the five cancer cells with IC50 values in the range of 0.61-9.50 µM. Selective proliferative inhibitions were observed on garcowacinol F and 7-epiclusianone against KB cells, and guttiferone Q toward MCF-7 cells with IC50 values less than 10 µM.

Roles of Salicylate Donors in Enhancement of Productivity and Isotacticity of Ziegler-Natta Catalyzed Propylene Polymerization.

Roles of internal salicylate donors (SID) in enhancing activity and stereoselectivity of Ziegler-Natta catalyzed propylene (PP) polymerization were examined using DFT calculations. Five salicylate donors were studied. The chelate mode is the preferred adsorption mode. The linear relationship (R2 = 0.96) between calculated adsorption energies (Eads) of five SIDs and the experimental PP activities was observed. Thus, the SID with the strongest adsorption energy will provide the highest activity in agreement with our previous studies. Compared with diisobutyl phthalate (DIBP), which is the industrial electron donor, SID has stronger Eads. The insertion step, which involves the π-complex formation (∆Eπ and the insertion activation or intrinsic activation energy (Ea) for PP polymerization was also examined. The relation between ln(activity) and apparent activation energy (Ea(app)), which is ∆Eπ + Ea for the primary(1,2)-re insertion with R2 = 0.99, was observed. The salicylate donor also has a lower Ea(app) than that of DIBP. This explains the better catalytic performance of SID. Our results also demonstrated that the size and the type of hydrocarbon substituents play a key role in controlling stereoselectivity and activity. In addition, we found a good relationship between Eads and both intrinsic (Ea) and apparent (Ea(app)) activation energies of five salicylate donors with R2 of 0.90 and 0.97, respectively.

Effects of shape, size, and pyrene doping on electronic properties of graphene nanoflakes.

Effects of size, shape, and pyrene doping on electronic properties of graphene nanoflakes (GNFs) were theoretically investigated using density functional theory method with PBE, B3PW91, and M06-2X functionals and cc-pVDZ basis set. Two shapes of zigzag GNFs, hexagonal (HGN) and rhomboidal (RGN), were considered. The energy band gap of GNF depends on shape and decreases with size. The HGN has larger band gap energy (1.23-3.96 eV) than the RGN (0.13-2.12 eV). The doping of pyrene and pyrene derivatives on both HGN and RGN was also studied. The adsorption energy of pyrene and pyrene derivatives on GNF does not depend on the shape of GNFs with energies between 21 and 27 kcal mol-1. The substituent on pyrene enhances the binding to GNF but the strength does not depend on electron withdrawing or donating capability. The doping by pyrene and pyrene derivatives also shifts the HOMO and LUMO energies of GNFs. Both positive (destabilizing) and negative (stabilizing) shifts on HOMO and LUMO of GNFs were seen. The direction and magnitude of the shift do not follow the electron withdrawing and donating capability of pyrene substituents. However, only a slight shift was observed for doped RGN. A shift of 0.19 eV was noticed for HOMO of HGN doped with 1-aminopyrene (pyNH2) and of 0.04 eV for LUMO of HGN doped with 1-pyrenecarboxylic acid (pyCOOH). Graphical Abstract HOMO and LUMO Energies of pyrene/pyrene derivatives doped Graphene Nanoflakes.

How strong is the edge effect in the adsorption of anticancer drugs on a graphene cluster?

The adsorption of nucleobase-analog anticancer drugs (fluorouracil, thioguanine, and mercaptopurine) on a graphene flake (C54H18) was investigated by shifting the site at which adsorption occurs from one end of the sheet to the other end. The counterpoise-corrected M06-2X/cc-pVDZ binding energies revealed that the binding stability decreases in the sequence thioguanine > mercaptopurine > fluorouracil. We found that adsorption near the middle of the sheet is more favorable than adsorption near the edge due to the edge effect. This edge effect is stronger for the adsorption of thioguanine or mercaptopurine than for fluorouracil adsorption. However, the edge effect reduces the binding energy of the drug to the flake by only a small amount, <5 kcal/mol, depending on the adsorption site and the alignment of the drug at this site.

Glucose transformation to 5-hydroxymethylfurfural in acidic ionic liquid: A quantum mechanical study.

Isomerization and transformation of glucose and fructose to 5-hydroxymethylfurfural (HMF) in both ionic liquids (ILs) and water has been studied by the reference interaction site model self-consistent field spatial electron density distribution (RISM-SCF-SEDD) method coupled with ab initio electronic structure theory, namely coupled cluster single, double, and perturbative triple excitation (CCSD(T)). Glucose isomerization to fructose has been investigated via cyclic and open chain mechanisms. In water, the calculations support the cyclic mechanism of glucose isomerization; with the predicted activation free energy is 23.8 kcal mol(-1) at experimental condition. Conversely, open ring mechanism is more favorable in ILs with the energy barrier is 32.4 kcal mol(-1) . Moreover, the transformation of fructose into HMF via cyclic mechanism is reasonable; the calculated activation barriers are 16.0 and 21.5 kcal mol(-1) in aqueous and ILs solutions, respectively. The solvent effects of ILs could be explained by the decomposition of free energies and radial distribution functions of solute-solvent that are produced by RISM-SCF-SEDD.

The Biginelli Reaction Is a Urea-Catalyzed Organocatalytic Multicomponent Reaction.

The recently developed artificial force induced reaction (AFIR) method was applied to search systematically all possible multicomponent pathways for the Biginelli reaction mechanism. The most favorable pathway starts with the condensation of the urea and benzaldehyde, followed by the addition of ethyl acetoacetate. Remarkably, a second urea molecule catalyzes nearly every step of the reaction. Thus, the Biginelli reaction is a urea-catalyzed multicomponent reaction. The reaction mechanism was found to be identical in both protic and aprotic solvents.

An n-channel two-dimensional covalent organic framework.

Co-condensation of metallophthalocyanine with an electron-deficient benzothiadiazole (BTDA) block leads to the formation of a two-dimensional covalent organic framework (2D-NiPc-BTDA COF) that assumes a belt shape and consists of AA stacking of 2D polymer sheets. Integration of BTDA blocks at the edges of a tetragonal metallophthalocyanine COF causes drastic changes in the carrier-transport mode and a switch from a hole-transporting skeleton to an electron-transporting framework. 2D-NiPc-BTDA COF exhibits broad and enhanced absorbance up to 1000 nm, shows panchromatic photoconductivity, is highly sensitive to near-infrared photons, and has excellent electron mobility as high as 0.6 cm(2) V(-1) s(-1).

Theoretical studies on the molecular basis of HIV-1RT/NNRTIs interactions.

Molecular dynamics simulations (MD) of the human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) complexed with the four non-nucleoside reverse transcriptase inhibitors (NNRTIs): efavirenz (EFV), emivirine (EMV), etravirine (ETV) and nevirapine (NVP), were performed to examine the structures, binding free energies and the importance of water molecules in the binding site. The binding free energy, calculated using molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), was found to decrease in the following order: EFV ~ ETV > EMV > NVP. The decrease in stability of the HIV-1 RT/NNRTI complexes is in good agreement with the experimentally derived half maximal inhibitory concentration (IC(50)) values. The interaction energy of the protein-inhibitor complexes was found to be essentially associated with the cluster of seven hydrophobic residues, L100, V106, Y181, Y188, F227, W229 and P236, and two basic residues, K101 and K103. Moreover, these residues are considered to be the most frequently detected mutated amino acids during treatment by various NNRTIs and therefore, those most likely to have been selected in the population for resistance.

Theoretical investigations on the stereoselectivity of the proline catalyzed Mannich reaction in DMSO.

The stereocontrol steps of the (S)-proline catalyzed Mannich reaction of cyclohexanone, formaldehyde, and aniline were theoretically investigated. The geometries of reactants, products, and transition states were optimized using density functional theory using the B3LYP functional with the 6-31++G(d,p) basis set. The energies of these compounds were then more accurately determined at the MP2 level, and the effect of DMSO as the solvent was included using a polarizable continuum model (PCM). The reaction was modeled from the previously proposed mechanism that cyclohexanone reacts with (S)-proline to generate an enamine, while formaldehyde reacts with aniline to produce an imine, and that the conformation around the C-N bond of the enamine 1 is crucial for the further enantioselective step. The formation of two conformations of the enamine via a proton transfer process was examined, revealing activation barriers for syn- and anti-enamine proton transfer of 10.2 and 17.9 kcal/mol, respectively. The transformation of syn- to anti-enamine through C-N bond rotation, however, was predicted to require only 4.2 kcal/mol, while the (S)- and (R)-intermediates could be obtained from subsequent reactions between enamine and imine with energy barriers of 8.5 and 12.4 kcal/mol, respectively. The difference between these barriers, but not the C-N rotation energy, becomes larger at the MP2 level and when DMSO as a solvent is included. This predicted enantioselective reaction, through the kinetic and thermodynamic favoring of the (S)-pathway, is in agreement with experimental results, which have reported the (S)-configuration as the major product.

Source of high pathogenicity of an avian influenza virus H5N1: why H5 is better cleaved by furin.

The origin of the high pathogenicity of an emerging avian influenza H5N1 due to the -RRRKK- insertion at the cleavage loop of the hemagglutinin H5, was studied using the molecular dynamics technique, in comparison with those of the noninserted H5 and H3 bound to the furin (FR) active site. The cleavage loop of the highly pathogenic H5 was found to bind strongly to the FR cavity, serving as a conformation suitable for the proteolytic reaction. With this configuration, the appropriate interatomic distances were found for all three reaction centers of the enzyme-substrate complex: the arrangement of the catalytic triad, attachment of the catalytic Ser(368) to the reactive S1-Arg, and formation of the oxyanion hole. Experimentally, the--RRRKK--insertion was also found to increase in cleavage of hemagglutinin by FR. The simulated data provide a clear answer to the question of why inserted H5 is better cleaved by FR than the other subtypes, explaining the high pathogenicity of avian influenza H5N1.

Hybrid quantum mechanical/molecular mechanical molecular dynamics simulations of HIV-1 integrase/inhibitor complexes.

Human immunodeficiency virus (HIV)-1 integrase (IN) is an attractive target for development of acquired immunodeficiency syndrome chemotherapy. In this study, conventional and coupled quantum mechanical and molecular mechanical (QM/MM) molecular dynamics (MD) simulations of HIV-1 IN complexed with 5CITEP (IN-5CITEP) were carried out. In addition to differences in the bound position of 5CITEP, significant differences at the two levels of theory were observed in the metal coordination geometry and the areas involving residues 116-119 and 140-166. In the conventional MD simulation, the coordination of Mg(2+) was found to be a near-perfect octahedral geometry whereas a distorted octahedral complex was observed in QM/MM. All of the above reasons lead to a different pattern of protein-ligand salt link formation that was not observed in the classical MD simulation. Furthermore to provide a theoretical understanding of inhibition mechanisms of 5CITEP and its derivative (DKA), hybrid QM/MM MD simulations of the two complexes (IN-5CITEP and IN-DKA) have been performed. The results reveal that areas involving residues 60-68, 116-119, and 140-149 were substantially different among the two systems. The two systems show similar pattern of metal coordination geometry, i.e., a distorted octahedron. In IN-DKA, both OD1 and OD2 of Asp-64 coordinate the Mg(2+) in a monodentate fashion whereas only OD1 is chelated to the metal as observed in IN-5CITEP. The high potency of DKA as compared to 5CITEP is supported by a strong salt link formed between its carboxylate moiety and the ammonium group of Lys-159. Detailed comparisons between HIV-1 IN complexed with DKA and with 5CITEP provide information about ligand structure effects on protein-ligand interactions in particular with the Lys-159. This is useful for the design of new selective HIV-1 IN inhibitors.

On the lower susceptibility of oseltamivir to influenza neuraminidase subtype N1 than those in N2 and N9.

Aiming to understand, at the molecular level, why oseltamivir (OTV) cannot be used for inhibition of human influenza neuraminidase subtype N1 as effectively as for subtypes N2 and N9, molecular dynamics simulations were carried out for the three complexes, OTV-N1, OTV-N2, and OTV-N9. The three-dimensional OTV-N2 and OTV-N9 initial structures were represented by the x-ray structures, whereas that of OTV-N1, whose x-ray structure is not yet solved, was built up using the aligned sequence of H5N1 isolated from humans in Thailand with the x-ray structure of the N2-substrate as the template. In comparison to the OTV-N2 and OTV-N9 complexes, dramatic changes were observed in the OTV conformation in the OTV-N1 complex in which two of its bulky side chains, N-acethyl (-NHAc) and 1-ethylproxy group (-OCHEt2), were rotated to adjust the size to fit into the N1 catalytic site. This change leads directly to the rearrangements of the OTV's environment, which are i), distances to its neighbors, W-178 and E-227, are shorter whereas those to residues R-224, E-276, and E-292 are longer; ii), hydrogen bonds to the two nearest neighbors, R-224 and E-276, are still conserved in distance and number as well as percentage occupation; iii), the calculated ligand/enzyme binding free energies of -7.20, -13.44, and -13.29 kcal/mol agree with their inhibitory activities in terms of the experimental IC50 of 36.1-53.2 nM, 1.9-2.7 nM, and 9.5-17.7 nM for the OTV-N1, OTV-N2, and OTV-N9 complexes, respectively; and iv), hydrogen-bond breaking and creation between the OTV and neighborhood residues are accordingly in agreement with the ligand solvation/desolvation taking place in the catalytic site.

Three-dimensional quantitative structure: activity relationship studies on diverse structural classes of HIV-1 integrase inhibitors using CoMFA and CoMSIA.

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA), three-dimensional quantitative structure-activity relationship (3D-QSAR) techniques, were applied to a set of 89 HIV-1 integrase (IN) inhibitors (training set=61, test set=28), belonging to 11 structurally different classes. The biological data for 3' processing mechanism were used. For CoMFA calculations, three different fitting methods for alignment process were investigated. The best CoMFA model yielded the cross-validated r(2) r(2)(cv) =0.698 and the non-cross-validated r(2) (r(2))=0.947. The derived model indicated the importance of steric (60.8%) as well as electrostatic (39.2%) contributions. For CoMSIA calculations, different combinations of the fields were tested. The best CoMSIA model gave r(2)(cv) =0.724 and r(2)=0.864. This model showed that steric (30.3%), hydrogen bond donor (43.4%) and hydrogen bond acceptor (26.3%) properties played major roles in HIV-1 IN inhibition. The mapping of hydrogen bond interaction fields with the HIV-1 IN active site gave details on hydrogen bond forming between ligands and enzyme. These obtained results agree well with the experimental observations that there should be hydrogen bond interactions between ligands and Glu152, Lys156 and Lys159 residues. The results not only lead to a better understanding of structural requirements of HIV-1 IN inhibitors but also can help in the design of new IN inhibitors.

Effective discrimination of antimalarial potency of artemisinin compounds based on quantum chemical calculations of their reaction mechanism.

The reaction mechanism of 12 antimalarial artemisinin compounds with two competitive pathways was studied by means of quantum chemical calculations using the IMOMO(B3LYP/6-31(d,p):HF/3-21G) method. The oxygen-centered radicals, carbon-centered radicals, and transition states (TS) in both pathways were geometrically optimized. The obtained kinetic and thermodynamic energy profiles show that homolytic C-C cleavage reaction (pathway 2) is energetically more preferable than an intramolecular 1,5-hydrogen shift process (pathway 1), which is consistent with the docking calculations. However, compounds that can easily proceed along the pathway 1 have high activity. Therefore, both pathways are important for antimalarial activity. Moreover, effective discrimination between high and low activity compounds using EA1, deltaE1, and deltaE(1A-2A) was accomplished.

Structure, dynamics and solvation of HIV-1 protease/saquinavir complex in aqueous solution and their contributions to drug resistance: molecular dynamic simulations.

As it is known that the understanding of the basic properties of the enzyme/inhibitor complex leads directly to enhancing the capability in drug designing and drug discovery. Molecular dynamics simulations have been performed to examine detailed information on the structure and dynamical properties of the HIV-1 PR complexed with saquinavir in the three protonated states, monoprotonates at Asp25 (Mono-25) and Asp25'(Mono-25') and diprotonate (Di-Pro) at both Asp25 and Asp25'. The obtained results support clinical data which reveal that Ile84 and Gly48 are two of the most frequent residues where mutation toward a protease inhibitor takes place. In contrast to the Ile84 mutation due to high displacement of Ile84 in the presence of saquinavir, source of the Gly48 mutation was observed to be due to the limited space in the HIV-1 PR pocket. The Gly48 was, on one side, found to form strong hydrogen bonds with saquinavir, while on the other side this residue was repelled by the hydrophobic Phe53 residue. In terms of inhibitor/enzyme binding, interactions between saquinavir and a catalytic triad of the HIV-1 PR were calculated using the ab initio method. The results show an order of the binding energy of Mono25

Insights into saquinavir resistance in the G48V HIV-1 protease: quantum calculations and molecular dynamic simulations.

The spread of acquired immune deficiency syndrome has increasingly become a great concern owing largely to the failure of chemotherapies. The G48V is considered the key signature residue mutation of HIV-1 protease developing with saquinavir therapy. Molecular dynamics simulations of the wild-type and the G48V HIV-1 protease complexed with saquinavir were carried out to explore structure and interactions of the drug resistance. The molecular dynamics results combined with the quantum-based and molecular mechanics Poisson-Boltzmann surface area calculations indicated a monoprotonation took place on D25, one of the triad active site residues. The inhibitor binding of the triad residues and its interaction energy in the mutant were similar to those in the wild-type. The overall structure of both complexes is almost identical. However, the steric conflict of the substituted valine results in the conformational change of the P2 subsite and the disruption of hydrogen bonding between the -NH of the P2 subsite and the backbone -CO of the mutated residue. The magnitude of interaction energy changes was comparable to the experimental K(i) data. The designing for a new drug should consider a reduction of steric repulsion on P2 to enhance the activity toward this mutant strain.