Journal of Computer Aided Molecular Design

Influence of metal cofactors and water on the catalytic mechanism of creatininase-creatinine in aqueous solution from molecular dynamics simulation and quantum study.

Publication Date: 2010 Aug 28 PMID: 20803053
Authors: Lee, V. S. - Kodchakorn, K. - Jitonnom, J. - Nimmanpipug, P. - Kongtawelert, P. - Premanode, B.
Journal: J Comput Aided Mol Des

The reaction mechanism of creatinine-creatininase binding to form creatine as a final product has been investigated by using a combined ab initio quantum mechanical/molecular mechanical approach and classical molecular dynamics (MD) simulations. In MD simulations, an X-ray crystal structure of the creatininase/creatinine was modified for creatininase/creatinine complexes and the MD simulations were run for free creatininase and creatinine in water. MD results reveal that two X-ray water molecules can be retained in the active site as catalytic water. The binding free energy from Molecular Mechanics Poisson-Boltzmann Surface Area calculation predicted the strong binding of creatinine with Zn(2+), Asp45 and Glu183. Two step mechanisms via Mn(2+)/Zn(2+) (as in X-ray structure) and Zn(2+)/Zn(2+) were proposed for water adding step and ring opening step with two catalytic waters. The pathway using synchronous transit methods with local density approximations with PWC functional for the fragment in the active region were obtained. Preferable pathway Zn(2+)/Zn(2+) was observed due to lower activation energy in water adding step. The calculated energy in the second step for both systems were comparable with the barrier of 26.03 and 24.44 kcal/mol for Mn(2+)/Zn(2+) and Zn(2+)/Zn(2+), respectively.

post to: CiteULike

Rationalizing fragment based drug discovery for BACE1: insights from FB-QSAR, FB-QSSR, multi objective (MO-QSPR) and MIF studies.

Publication Date: 2010 Aug 26 PMID: 20740315
Authors: Manoharan, P. - Vijayan, R. S. - Ghoshal, N.
Journal: J Comput Aided Mol Des

The ability to identify fragments that interact with a biological target is a key step in FBDD. To date, the concept of fragment based drug design (FBDD) is increasingly driven by bio-physical methods. To expand the boundaries of QSAR paradigm, and to rationalize FBDD using In silico approach, we propose a fragment based QSAR methodology referred here in as FB-QSAR. The FB-QSAR methodology was validated on a dataset consisting of 52 Hydroxy ethylamine (HEA) inhibitors, disclosed by GlaxoSmithKline Pharmaceuticals as potential anti-Alzheimer agents. To address the issue of target selectivity, a major confounding factor in the development of selective BACE1 inhibitors, FB-QSSR models were developed using the reported off target activity values. A heat map constructed, based on the activity and selectivity profile of the individual R-group fragments, and was in turn used to identify superior R-group fragments. Further, simultaneous optimization of multiple properties, an issue encountered in real-world drug discovery scenario, and often overlooked in QSAR approaches, was addressed using a Multi Objective (MO-QSPR) method that balances properties, based on the defined objectives. MO-QSPR was implemented using Derringer and Suich desirability algorithm to identify the optimal level of independent variables (X) that could confer a trade-off between selectivity and activity. The results obtained from FB-QSAR were further substantiated using MIF (Molecular Interaction Fields) studies. To exemplify the potentials of FB-QSAR and MO-QSPR in a pragmatic fashion, the insights gleaned from the MO-QSPR study was reverse engineered using Inverse-QSAR in a combinatorial fashion to enumerate some prospective novel, potent and selective BACE1 inhibitors.

post to: CiteULike

QMOD: physically meaningful QSAR.

Publication Date: 2010 Aug 19 PMID: 20721601
Authors: Jain, A. N.
Journal: J Comput Aided Mol Des

Computational methods for predicting ligand affinity where no protein structure is known generally take the form of regression analysis based on molecular features that have only a tangential relationship to a protein/ligand binding event. Such methods have utility in retrospective rationalization of activity patterns of substituents on a common scaffold, but are limited when either multiple scaffolds are present or when ligand alignment varies significantly based on structural changes. In addition, such methods generally assume independence and additivity of effect from scaffold substituents. Collectively, these non-physical modeling assumptions sharply limit the utility of widely used QSAR approaches for prospective prediction of ligand activity. The recently introduced Surflex-QMOD approach, by virtue of constructing physical models of binding sites, comes closer to a modeling approach that is congruent with protein ligand binding events. A set of congeneric CDK2 inhibitors showed that induced binding pockets can be quite congruent with the enzyme's active site but that model predictivity within a chemical series does not necessarily depend on congruence. Muscarinic antagonists were used to show that the QMOD approach is capable of making accurate predictions in cases where highly non-additive structure activity effects exist. The QMOD method offers a means to go beyond non-causative correlations in QSAR analysis.

post to: CiteULike

Mechanisms of amphipathic helical peptide denaturation by guanidinium chloride and urea: a molecular dynamics simulation study.

Publication Date: 2010 Aug 10 PMID: 20697778
Authors: Mehrnejad, F. - Khadem-Maaref, M. - Ghahremanpour, M. M. - Doustdar, F.
Journal: J Comput Aided Mol Des

Urea and GdmCl are widely used to denature proteins at high concentrations. Here, we used MD simulations to study the denaturation mechanisms of helical peptide in different concentrations of GdmCl and urea. It was found that the helical structure of the peptide in water simulation is disappeared after 5 ns while the helicity of the peptide is disappeared after 70 ns in 2 M urea and 25 ns in 1 M GdmCl. Surprisingly, this result shows that the helical structure in low concentration of denaturants is remained more with respect to that solvated in water. The present work strongly suggests that urea interact more preferentially to non-polar and aromatic side chains in 2 M urea; therefore, hydrophobic residues are in more favorable environment in 2 M urea. Our results also reveal that the hydrogen bonds between urea and the backbone is the dominant mechanism by which the peptide is destabilized in high concentration of urea. In 1 M and 2 M GdmCl, GdmCl molecules tend to engage in transient stacking interactions with aromatics and hydrophobic planar side chains that lead to displacement of water from the hydration surface, providing more favorable environment for them. This shows that accumulation of GdmCl around hydrophobic surfaces in 1 M and 2 M GdmCl solutions prevents proper solvation of the peptide at the beginning. In high GdmCl concentrations, water solvate the peptide better than 1 M and 2 M GdmCl. Therefore, our results strongly suggest that hydrogen bonds between water and the peptide are important factors in the destabilization of peptide in GdmCl solutions.

post to: CiteULike

T-Analyst: a program for efficient analysis of protein conformational changes by torsion angles.

Publication Date: 2010 Aug 6 PMID: 20689979
Authors: Ai, R. - Qaiser Fatmi, M. - Chang, C. E.
Journal: J Comput Aided Mol Des

T-Analyst is a user-friendly computer program for analyzing trajectories from molecular modeling. Instead of using Cartesian coordinates for protein conformational analysis, T-Analyst is based on internal bond-angle-torsion coordinates in which internal torsion angle movements, such as side-chain rotations, can be easily detected. The program computes entropy and automatically detects and corrects angle periodicity to produce accurate rotameric states of dihedrals. It also clusters multiple conformations and detects dihedral rotations that contribute hinge-like motions. Correlated motions between selected dihedrals can also be observed from the correlation map. T-Analyst focuses on showing changes in protein flexibility between different states and selecting representative protein conformations for molecular docking studies. The program is provided with instructions and full source code in Perl.

post to: CiteULike

Homology modeling and metabolism prediction of human carboxylesterase-2 using docking analyses by GriDock: a parallelized tool based on AutoDock 4.0.

Publication Date: 2010 Sep PMID: 20623318
Authors: Vistoli, G. - Pedretti, A. - Mazzolari, A. - Testa, B.
Journal: J Comput Aided Mol Des

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted to developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas less has been done to predict the activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES2. The study involved first a homology modeling of the hCES2 protein based on the model of hCES1 since the two proteins share a high degree of homology (congruent with 73%). A set of 40 known substrates of hCES2 was taken from the literature; the ligands were docked in both their neutral and ionized forms using GriDock, a parallel tool based on the AutoDock4.0 engine which can perform efficient and easy virtual screening analyses of large molecular databases exploiting multi-core architectures. Useful statistical models (e.g., r (2) = 0.91 for substrates in their unprotonated state) were calculated by correlating experimental pK(m) values with distance between the carbon atom of the substrate's ester group and the hydroxy function of Ser228. Additional parameters in the equations accounted for hydrophobic and electrostatic interactions between substrates and contributing residues. The negatively charged residues in the hCES2 cavity explained the preference of the enzyme for neutral substrates and, more generally, suggested that ligands which interact too strongly by ionic bonds (e.g., ACE inhibitors) cannot be good CES2 substrates because they are trapped in the cavity in unproductive modes and behave as inhibitors. The effects of protonation on substrate recognition and the contrasting behavior of substrates and products were finally investigated by MD simulations of some CES2 complexes.

post to: CiteULike

An efficient synthesis of a rationally designed 1,5 disubstituted imidazole AT(1) angiotensin II receptor antagonist: reorientation of imidazole pharmacophore groups in losartan reserves high receptor affinity and confirms docking studies.

Publication Date: 2010 Sep PMID: 20623162
Authors: Agelis, G. - Roumelioti, P. - Resvani, A. - Durdagi, S. - Androutsou, M. E. - Kelaidonis, K. - Vlahakos, D. - Mavromoustakos, T. - Matsoukas, J.
Journal: J Comput Aided Mol Des

A new 1,5 disubstituted imidazole AT(1) Angiotensin II (AII) receptor antagonist related to losartan with reversion of butyl and hydroxymethyl groups at the 2-, 5-positions of the imidazole ring was synthesized and evaluated for its antagonist activity (V8). In vitro results indicated that the reorientation of butyl and hydroxymethyl groups on the imidazole template of losartan retained high binding affinity to the AT(1) receptor concluding that the spacing of the substituents at the 2,5- positions is of primary importance. The docking studies are confirmed by binding assay results which clearly show a comparable binding score of the designed compound V8 with that of the prototype losartan. An efficient, regioselective and cost effective synthesis renders the new compound as an attractive candidate for advanced toxicological evaluation and a drug against hypertension.

post to: CiteULike

ElectroShape: fast molecular similarity calculations incorporating shape, chirality and electrostatics.

Publication Date: 2010 Sep PMID: 20614163
Authors: Armstrong, M. S. - Morris, G. M. - Finn, P. W. - Sharma, R. - Moretti, L. - Cooper, R. I. - Richards, W. G.
Journal: J Comput Aided Mol Des

We present ElectroShape, a novel ligand-based virtual screening method, that combines shape and electrostatic information into a single, unified framework. Building on the ultra-fast shape recognition (USR) approach for fast non-superpositional shape-based virtual screening, it extends the method by representing partial charge information as a fourth dimension. It also incorporates the chiral shape recognition (CSR) method, which distinguishes enantiomers. It has been validated using release 2 of the Directory of useful decoys (DUD), and shows a near doubling in enrichment ratio at 1% over USR and CSR, and improvements as measured by Receiver Operating Characteristic curves. These improvements persisted even after taking into account the chemotype redundancy in the sets of active ligands in DUD. During the course of its development, ElectroShape revealed a difference in the charge allocation of the DUD ligand and decoy sets, leading to several new versions of DUD being generated as a result. ElectroShape provides a significant addition to the family of ultra-fast ligand-based virtual screening methods, and its higher-dimensional shape recognition approach has great potential for extension and generalisation.

post to: CiteULike

In silico analysis of the histaprodifen induced activation pathway of the guinea-pig histamine H(1)-receptor.

Publication Date: 2010 Sep PMID: 20607589
Authors: Strasser, A. - Wittmann, H. J.
Journal: J Comput Aided Mol Des

The binding of (partial) agonists in the binding pocket of biogenic amine receptors induces a conformational change from the inactive to the active state of the receptors. There is only little knowledge about the binding pathways of ligands into binding pocket on molecular level. So far, it was not possible with molecular dynamic simulations to observe the ligand binding and receptor activation. Furthermore, there is nearly nothing known, in which state of ligand binding, the receptor gets activated. The aim of this study was to get more detailed insight into the process of ligand binding and receptor activation. With the recently developed LigPath algorithm, we scanned the potential energy surface of the binding process of dimeric histaprodifen, a partial agonist at the histamine H(1)-receptor, into the guinea pig histamine H(1)-receptor, taking also into account the receptor activation. The calculations exhibited large conformational changes of Trp(6.48) and Phe(6.55) during ligand binding and receptor activation. Additionally, conformational changes were also observed for Phe(6.52), Tyr(6.51) and Phe(6.44). Conformational changes of Trp(6.48) and Phe(6.52) are discussed in literature as rotamer toggle switch in context with receptor activation. Additionally, the calculations indicate that the binding of dimeric histaprodifen, accompanied by receptor activation is energetically preferred. In general, this study gives new, theoretical insights onto ligand binding and receptor activation on molecular level.

post to: CiteULike

Phosphorylation and ATP-binding induced conformational changes in the PrkC, Ser/Thr kinase from B. subtilis.

Publication Date: 2010 Sep PMID: 20563625
Authors: Gruszczynski, P. - Obuchowski, M. - Kazmierkiewicz, R.
Journal: J Comput Aided Mol Des

Recent studies on the PrkC, serine-threonine kinase show that that the enzyme is located at the inner membrane of endospores and is responsible for triggering spore germination. The activity of the protein increases considerably after phosphorylation of four threonine residues placed on the activation loop and one serine placed in the C-terminal lobe of the PrkC. The molecular relationship between phosphorylation of these residues and enzyme activity is not known. In this work molecular dynamics simulation is performed on four forms of the protein kinase PrkC from B. subtilis-phosphorylated or unphosphorylated; with or without ATP bound-in order to gain insight into phosphorylation and ATP binding on the conformational changes and functions of the protein kinase. Our results show how phosphorylation, as well as the presence of ATP, is important for the activity of the enzyme through its molecular interaction with the catalytic core residues. Three of four threonine residues were found to be involved in the interactions with conservative motifs important for the enzyme activity. Two of the threonine residues (T167 and T165) are involved in ionic interactions with an arginine cluster from alphaC-helix. The third residue (T163) plays a crucial role, interacting with His-Arg-Asp triad (HRD). Last of the threonine residues (T162), as well as the serine (S214), were indicated to play a role in the substrate recognition or dimerization of the enzyme. The presence of ATP in the unphosphorylated model induced conformational instability of the activation loop and Asp-Phe-Gly motif (DFG). Based on our calculations we put forward a hypothesis suggesting that the ATP binds after phosphorylation of the activation loop to create a fully active conformation in the closed position.

post to: CiteULike

Flexibility of a biotinylated ligand in artificial metalloenzymes based on streptavidin--an insight from molecular dynamics simulations with classical and ab initio force fields.

Publication Date: 2010 Sep PMID: 20526651
Authors: Panek, J. J. - Ward, T. R. - Jezierska-Mazzarello, A. - Novic, M.
Journal: J Comput Aided Mol Des

In the field of enzymatic catalysis, creating activity from a non catalytic scaffold is a daunting task. Introduction of a catalytically active moiety within a protein scaffold offers an attractive means for the creation of artificial metalloenzymes. With this goal in mind, introduction of a biotinylated d(6)-piano-stool complex within streptavidin (SAV) affords enantioselective artificial transfer-hydrogenases for the reduction of prochiral ketones. Based on an X-ray crystal structure of a highly selective hybrid catalyst, displaying significant disorder around the biotinylated catalyst [eta(6)-(p-cymene)Ru(Biot-p-L)Cl], we report on molecular dynamics simulations to shed light on the protein-cofactor interactions and contacts. The results of these simulations with classical force field indicate that the SAV-biotin and SAV-catalyst complexes are more stable than ligand-free SAV. The point mutations introduced did not affect significantly the overall behavior of SAV and, unexpectedly, the P64G substitution did not provide additional flexibility to the protein scaffold. The metal-cofactor proved to be conformationally flexible, and the S112K or P64G mutants proved to enhance this effect in the most pronounced way. The network of intermolecular hydrogen bonds is efficient at stabilizing the position of biotin, but much less at fixing the conformation of an extended biotinylated ligand. This leads to a relative conformational freedom of the metal-cofactor, and a poorly localized catalytic metal moiety. MD calculations with ab initio potential function suggest that the hydrogen bonds alone are not sufficient factors for full stabilization of the biotin. The hydrophobic biotin-binding pocket (and generally protein scaffold) maintains the hydrogen bonds between biotin and protein.

post to: CiteULike