The effect of Li+ on GSK-3 inhibition: molecular dynamics simulation.

Glycogen synthase kinase-3 (GSK-3) is a kind of serine-threonine protein kinase. It places important roles in several signaling pathways and it is a key therapeutic target for a number of diseases, such as diabetes, cancer, Alzheimer's disease and chronic inflammation. Mg(2+) ions which interact with ATP are conserved in GSK. They are important in phosphoryl transfer. Li(+) is an inhibitor for GSK-3. It is used to treat bipolar mood disorder. This paper illustrates the effect of Li(+) on GSK-3. When Mg(I)(2+) is replaced by Li(+), the atom fluctuation of GSK-3 will rise, and the in-line phosphoryl transfer mechanism is probably demolished and the binding of pre-phosphorylated substrates may be disturbed. All the results we obtained clearly suggest that inhibition to GSK-3 is caused by the Mg(I)(2+) replacement with Li(+).

Molecular docking and molecular dynamics simulation studies of GPR40 receptor-agonist interactions.

In order to explore the agonistic activity of small-molecule agonists to GPR40, AutoDock and GROMACS software were used for docking and molecular dynamics studies. A molecular docking of eight structurally diverse agonists (six carboxylic acids (CAs) agonist, and two non-carboxylic acids (non-CAs) agonist) was performed and the differences in their binding modes were investigated. Moreover, a good linear relationship based on the predicted binding affinities (pK(i)) determined by using AutoDock and experimental activity values (pEC50) was obtained. Then, the 10ns molecular dynamics (MD) simulations of three obtained ligand-receptor complexes embedded into the phospholipid bilayer were carried out. The position fluctuations of the ligands located inside the transmembrane domain were explored, and the stable binding modes of the three studied agonists were determined. Furthermore, the residue-based decomposition of interaction energies in three systems identified several critical residues for ligand binding.

Ab initio calculations on halogen-bonded complexes and comparison with density functional methods.

A systematic theoretical investigation on a series of dimeric complexes formed between some halocarbon molecules and electron donors has been carried out by employing both ab initio and density functional methods. Full geometry optimizations are performed at the Moller-Plesset second-order perturbation (MP2) level of theory with the Dunning's correlation-consistent basis set, aug-cc-pVDZ. Binding energies are extrapolated to the complete basis set (CBS) limit by means of two most commonly used extrapolation methods and the aug-cc-pVXZ (X = D, T, Q) basis sets series. The coupled cluster with single, double, and noniterative triple excitations [CCSD(T)] correction term, determined as a difference between CCSD(T) and MP2 binding energies, is estimated with the aug-cc-pVDZ basis set. In general, the inclusion of higher-order electron correlation effects leads to a repulsive correction with respect to those predicted at the MP2 level. The calculations described herein have shown that the CCSD(T) CBS limits yield binding energies with a range of -0.89 to -4.38 kcal/mol for the halogen-bonded complexes under study. The performance of several density functional theory (DFT) methods has been evaluated comparing the results with those obtained from MP2 and CCSD(T). It is shown that PBEKCIS, B97-1, and MPWLYP functionals provide accuracies close to the computationally very expensive ab initio methods.

Effect of mutation K85R on GSK-3beta: Molecular dynamics simulation.

As a serine-threonine protein kinase, glycogen synthase kinase-3 (GSK-3) regulates the synthesis of glycogen and plays important roles in several signaling pathways. It is a key therapeutic target for a number of diseases, such as diabetes, cancer, Alzheimer's disease and chronic inflammation. The conserved Lys85 is important to GSK-3beta activity and in this paper we illustrate the significant role of Lys85 using dynamic simulation. We find that when Lys85 is mutated to Arg, one of the two conserved hydrogen bonds between Lys85 and ATP disappears, the salt bridge between Lys85 and Glu97 cannot form, and conformational changes of Phe93, Arg96 and Glu211 occur. These will cause conformational changes of the substrate binding groove that would inhibit the activity of GSK-3beta. MM-GBSA calculations reveal that the K85R mutation could lead to a less energy-favorable complex, which is consistent with the structural analysis.

Quantitative structure-chromatographic retention relationship for polycyclic aromatic sulfur heterocycles.

Polycyclic aromatic sulfur heterocycles (PASHs) are of concern in petroleum geochemistry and environmental chemistry. In the present study, geometrical optimization and electrostatic potential calculations have been performed for 114 PASHs reported previously at the HF/6-31G* level of theory. A group of 25 statistically based parameters have been extracted. Linear relationships between gas-chromatographic retention index (RI) and the structural descriptors have been established by stepwise linear regression analysis. The result shows that two quantities derived from positive electrostatic potential on molecular surface, V(s)(+) (the average value of the positive electrostatic potentials on molecular surface) and sigma(+)(2) (a measure of dispersion tendency of positive electrostatic potential), together with V(mc) (the molecular volume) and E(HOMO) (the energy of the highest occupied molecular orbital) can be well used to express the quantitative structure-retention relationship (QSRR) of PASHs. Predictive capability of the model has been demonstrated by leave-one-out cross-validation with the cross-validated correlation coefficient (R(CV)) of 0.992. Furthermore, when splitting the 114 PASH samples into calibration and test sets in the ratio of 2:1, a similar treatment yields an equation of almost equal statistical quality and very similar regression coefficients, validating the robustness of our model. Predictions for six PASHs from other source have also been made. The QSRR model established may provide a new powerful method for predicting chromatographic properties of aromatic organosulfur compounds.

Identification of ligand features essential for HDACs inhibitors by pharmacophore modeling.

Histone deacetylases (HDACs) enzyme plays a significant role in transcriptional regulation by modifying the core histones of the nucleosome. It has emerged as an important therapeutic target for the treatment of cancer and other diseases. Inhibitors of HDACs become a new class of anticancer agents and have provoked much interest amongst pharmacologists and cancer researchers. To facilitate the discovery of specific HDACs inhibitors, a 3D chemical-feature-based QSAR pharmacophore model was developed and was well consistent with the structure-functional requirements for the binding of the HDAC inhibitors. Using this model, the interactions between the benzamide MS-275 and HDAC were explored. The result showed that the type and spatial location of chemical features encoded in the pharmacophore are in full agreement with the enzyme-inhibitor interaction pattern identified from molecular docking.

Ab initio investigation of the complexes between bromobenzene and several electron donors: some insights into the magnitude and nature of halogen bonding interactions.

Halogen bonding, a specific intermolecular noncovalent interaction, plays crucial roles in fields as diverse as molecular recognition, crystal engineering, and biological systems. This paper presents an ab initio investigation of a series of dimeric complexes formed between bromobenzene and several electron donors. Such small model systems are selected to mimic halogen bonding interactions found within crystal structures as well as within biological molecules. In all cases, the intermolecular distances are shown to be equal to or below sums of van der Waals radii of the atoms involved. Halogen bonding energies, calculated at the MP2/aug-cc-pVDZ level, span over a wide range, from -1.52 to -15.53 kcal/mol. The interactions become comparable to, or even prevail over, classical hydrogen bonding. For charge-assisted halogen bonds, calculations have shown that the strength decreases in the order OH- > F- > HCO2- > Cl- > Br-, while for neutral systems, their relative strengths attenuate in the order H2CS > H2CO > NH3 > H2S > H2O. These results agree with those of the quantum theory of atoms in molecules (QTAIM) since bond critical points (BCPs) are identified for these halogen bonds. The QTAIM analysis also suggests that strong halogen bonds are more covalent in nature, while weak ones are mostly electrostatic interactions. The electron densities at the BCPs are recommended as a good measure of the halogen bond strength. Finally, natural bond orbital (NBO) analysis has been applied to gain more insights into the origin of halogen bonding interactions.

Three-dimensional QSAR of HPPD inhibitors, PSA inhibitors, and anxiolytic agents: effect of tautomerism on the CoMFA models.

The present study was design to examine the effect of tautomerism upon the CoMFA results. Three selected data sets involving protropic tautomerism, which are 21 p-hydroxyphenylpyruvate dioxygenase (HPPD) inhibitors, 35 inhibitors of puromycin-sensitive aminopeptidase (PSA), and 67 anxiolytic agents, were used for this purpose. Atom-by-atom alignment technique was adopted to superimpose the molecules in the data sets onto a template. The structural alignments using different tautomeric forms had no significant difference except the atoms involved in tautomerism, which ensures, to a great extent, that the differences of the CoMFA results result primarily from the tautomerism. All-orientation and all-placement search (AOS-APS) based CoMFA models, in addition to the conventional ones, were derived for each system and proved to be capable of yielding much improved statistical results. In the cases of the data sets of HPPD inhibitors and PSA inhibitors, excellent AOS-APS CoMFA models (q2>0.8 with four components for the former and q2>0.7 with seven components for the latter) were obtained, and almost no significant difference in statistical quality was observed when using different tautomeric forms to derive the models. However, it was not the case when treating the data set of anxiolytic agents. The keto tautomer, which was the active form of the PBI type inhibitors, produced measurably better results (q2=0.54 with eight components) than that the enol one (q2=0.37 with five components), indicating the importance of selecting proper tautomer in the CoMFA studies. Furthermore, there existed some substantial differences of the electrostatic field contours between the two different tautomeric forms for all of the three systems considered, whereas the differences in the steric field contour maps were limited. This implies that the resulting new potent ligands may be quite different if one utilizes the CoMFA models of different tautomeric forms for guiding further structural refinements.

QSPR/QSAR models for prediction of the physicochemical properties and biological activity of polybrominated diphenyl ethers.

Polybrominated diphenyl ethers (PBDEs) are a group of important persistent organic pollutants. In the present study, geometrical optimization and electrostatic potential calculations have been performed for all 209 PBDE congeners at the HF/6-31G level of theory. A number of statistically-based parameters have been obtained. Linear relationships between gas-chromatographic relative retention time (RRT), n-octanol/air partition coefficient (lgK(OA)), 298 K supercooled liquid vapour pressures (lgp(L)), Henry's law constant (lgH) and Ah receptor binding affinity (-lgRBA) of PBDEs and the structural descriptors have been established by multiple regression method. The result shows that the quantities derived from electrostatic potential V(s,max),V(s,min),Pi, Sigma V+(S), V-(S) , nu, sigma 2(tot), and N-(v), together with the molecular volume (Vmc) can be well used to express the quantitative structure-property relationships of PBDEs, which proves the general applicability of this parameter set to a great extent. Good predictive capabilities have also been demonstrated. Based on these equations, the predicted values have been presented for those PBDE congeners whose experimentally determined physicochemical properties are unavailable. The QSAR model for the Ah receptor binding affinity is relatively poor, which can be ascribed to the complexity of factors which affect biological activity and the limitations of the present parameter set in describing steric characters of the molecule.

Proton exchanges between phenols and ammonia or amines: a computational study.

Density functional theory calculations at the B3LYP/6-31+G(d,p) level of theory have been performed to explore proton exchanges between phenols and ammonia or amines, which can be used to account for previous NMR experiments. For the parent phenol-NH(3) system, a transition state with a symmetric phenolate-NH(4)(+)-like structure, which lies about 35 kcal mol(-1) in energy above the hydrogen-bonded complex, has been successfully located. An intrinsic reaction coordinate (IRC) analysis indicates that the proton exchange is a concerted process, which can be roughly divided into four continuous subprocesses. A series of para-substituted phenol-NH(3) systems have been considered to investigate the substituent effect. Whereas introduction of an electron-withdrawing group on the phenol appreciably reduces the barrier, an opposite effect is observed for an electron-donating group. Moreover, it has been disclosed that there exists a good linear correlation between the activation barriers and the interaction energies between the phenols and NH(3), indicating the important role of proton transfer (or hydrogen bonding) in determining the proton exchange. Also considered are the proton exchanges between phenol and amines and those for some sterically hindered systems. The results show that the phenol tends to exchange hydrogen with the amines, preferably the secondary amines, and that the steric effect is favorable for the proton exchange, which imply that, as the IRC analysis suggested, besides the proton transfer, the flip of the ammonium-like moiety may play a significant role in the course of proton exchange. For all of these systems, we investigated the solvent effects and found that the barrier heights of proton exchange decrease remarkably as compared to those in a vacuum due to the ion pair feature of the transition state. Finally, we explored the phenol radical cation-NH(3) system; the barrierless proton transfer and remarkably low barrier (5.2 kcal mol(-1)) of proton exchange provide further evidence for the importance of proton transfer in the proton exchange.

Electrochemical characteristics of gatifloxacin and its interaction with DNA.

The electrochemical behavior of gatifloxacin (GTFX) and its interaction with natural calf thymus DNA (ctDNA) is investigated by differential pulse voltammetry (DPV) on a carbon paraffined electrode. According to the suggested electrochemical equation, a binding constant of 1.7058 x 10(5) (mol L(-1))(-1) and binding sizes s = 3.09 (base pairs) of GTFX with ctDNA are obtained by nonlinear fit analysis of electrochemical data. The results demonstrate that GTFX has the properties of an intercalative binder.

Effects of the dietary supplementation with fructooligosaccharides on the excretion of nitrogen and phosphorus in Miichthys miiuy fries.

Effects of dietary supplementation with fructooligosaccharides on the excretion of nitrogen and phosphorus in Miichthys miiuy fries were investigated. Nine hundred Miichthys miiuy fries were divided into 3 groups, each with triplicates. The basal diet and the basal diet supplemented with carnitine groups were considered as the negative and positive controls respectively. Results showed that the nitrogen concentration in excreted feces decreased significantly in fries fed the diet supplementation with 1000 x 10(-6) fructooligosaccharides and 200 x 10(-6) carnitine (P<0.05). The ammonic-nitrogen concentration decreased significantly in the carnitine group only (P<0.05), indicating the decreasing tendency caused by the supplementation with fructooligosaccharides. Supplementation with both did not have significant effects on the concentration of phosphorus in feces of Miichthys miiuy fries.

[Exclusion of the association of five known mutations with congenital stationary nyctalopia in a large Chinese family].

OBJECTIVE: To detect gene mutations associated with autosomal dominant congenital stationary night blindness(ADCSNB) in a large Chinese family. METHODS: Genomic DNAs were extracted from peripheral blood samples of 16 affected and 14 unaffected family members. According to 5 missense mutations in 3 genes reported previously, 4 pairs of primers were designed and corresponding exons containing the five mutation sites were amplified by polymerase chain reaction. Amplified products were purified and sequenced by MegaBACE1000 capillary array electrophoresis DNA sequencer. Full field electroretinogram (ERG, ISCEV) of patients was recorded and analyzed by Roland Consult System. RESULTS: Dark-adapted ERG showed a-wave was normal, but b-wave of the patients was markedly decreased. None of the five missense mutations were detected in 16 affected and 14 unaffected family members. CONCLUSION: The molecular pathogenesis of ADCSNB in this family does not involve point mutations or deletions of these five sites, which indicates the heterogeneity of ADCSNB.

Ab initio study of the complexes of halogen-containing molecules RX (X=Cl, Br, and I) and NH3: towards understanding the nature of halogen bonding and the electron-accepting propensities of covalently bonded halogen atoms.

Ab initio calculations have been performed on a series of complexes formed between halogen-containing molecules and ammonia to gain a deeper insight into the nature of halogen bonding. It appears that the dihalogen molecules form the strongest halogen-bonded complexes with ammonia, followed by HOX; the charge-transfer-type contribution has been demonstrated to dominate the halogen bonding in these complexes. For the complexes involving carbon-bound halogen molecules, our calculations clearly indicate that electrostatic interactions are mainly responsible for their binding energies. Whereas the halogen-bond strength is significantly enhanced by progressive fluorine substitution, the substitution of a hydrogen atom by a methyl group in the CH(3)X...NH(3) complex weakened the halogen bonding. Moreover, remote substituent effects have also been noted in the complexes of halobenzenes with different para substituents. The influence of the hybridization state of the carbon atom bonded to the halogen atom has also been examined and the results reveal that halogen-bond strengths decrease in the order HC triple bond CX > H(2)C=CHX approximately O=CHX approximately C(6)H(5)X > CH(3)X. In addition, several excellent linear correlations have been established between the interaction energies and both the amount of charge transfer and the electrostatic potentials corresponding to an electron density of 0.002 au along the R-X axis; these correlations provide good models with which to evaluate the electron-accepting abilities of the covalently bonded halogen atoms. Finally, some positively charged halogen-bonded systems have been investigated and the effect of the charge has been discussed.

Binding interaction of gatifloxacin with bovine serum albumin.

The binding of gatifloxacin to bovine serum albumin (BSA) in aqueous solution was studied using fluorescence spectroscopy and absorbance spectra, Further, the interactions influenced by Fe3+ and Cu2+ were also explored in this work. Based on Scatchard's site-binding model and florescence quenching, practical formulas for small molecule ligands to bio-macromolecules have been proposed. The binding parameters were measured according to suggested models, and the binding distance and the transfer efficiency of energy between gatifloxacin and BSA were also obtained in view of the Förster theory of non-radiation energy transfer. The effect of gatifloxacin on the conformation of BSA has also been analyzed using synchronous fluorescence spectroscopy.