A comparative molecular field analysis of the biotransformation of sulfides by Rhodococcus erythropolis

A comparative molecular field analysis (CoMFA) was used to model the efficacy with which the Rhodococcus erythropolis mono-oxygenase, DszC, catalyzes the enantioselective sulfoxidation of a broad range of substrates. Experimentally determined values of both the yield and enantiomeric excess for this reaction were employed to create these CoMFA models. A highly predictive CoMFA model was constructed for the prediction of enantiomeric excess of the sulfoxide product. The predictive ability of the model was demonstrated by both cross-validation of the training set (q² = 0.74) and for an external test set of substrates. The enantiomeric excesses of the members of the test set, which also included two amino acid sulfides that were structurally distinct from the membership of the training set, were predicted well by the CoMFA model. Product yield was not modelled well by any CoMFA model. Different models comparing the likely bioactive conformations of the substrates suggest that most compounds assume an ‘extended’ conformation upon binding. Contour diagrams illustrating significant substrate–enzyme interactions suggest that the model, which predicts the enantiomeric excess, is consistent with previous conclusions regarding the effect of various substrate substitutions on the enantiopurity of the product of the biotransformation.     

Synthesis and comparative molecular field analysis (CoMFA) of antitumor 3-arylisoquinoline derivatives

In this study a series of 3-arylisoquinoline derivatives were synthesized and cytotoxicity against human melanoma tumor cell evaluated, and a three dimensional quantitative structure—activity relationship was investigated using the comparative molecular field analysis (CoMFA). The results suggested that the electrostatic, steric and hydrophobic factors of 3-arylisoquinolines were strongly correlated with the antitumor activity. Considerable predictive ability (cross-validated r² as high as 0.721) was obtained through CoMFA.     

Three-dimensional structure-activity relationship modeling of cocaine binding to two monoclonal antibodies by comparative molecular field analysis

Successful immunotherapy of cocaine addiction and overdoses requires cocaine-binding antibodies with specific properties, such as high affinity and selectivity for cocaine. We have determined the affinities of two cocaine-binding murine monoclonal antibodies (mAb: clones 3P1A6 and MM0240PA) for cocaine and its metabolites by [3H]-radioligand binding assays. mAb 3P1A6 (K(d) = 0.22 nM) displayed a 50-fold higher affinity for cocaine than mAb MM0240PA (K(d) = 11 nM) and also had a greater specificity for cocaine. For the systematic exploration of both antibodies' binding specificities, we used a set of approximately 35 cocaine analogues as structural probes by determining their relative binding affinities (RBAs) using an enzyme-linked immunosorbent competition assay. Three-dimensional quantitative structure-activity relationship (3D-QSAR) models on the basis of comparative molecular field analysis (CoMFA) techniques correlated the binding data with structural features of the ligands. The analysis indicated that despite the mAbs' differing specificities for cocaine, the relative contributions of the steric (approximately 80%) and electrostatic (approximately 20%) field interactions to ligand-binding were similar. Generated three-dimensional CoMFA contour plots then located the specific regions about cocaine where the ligand/receptor interactions occurred. While the overall binding patterns of the two mAbs had many features in common, distinct differences were observed about the phenyl ring and the methylester group of cocaine. Furthermore, using previously published data, a 3D-QSAR model was developed for cocaine binding to the dopamine reuptake transporter (DAT) that was compared to the mAb models. Although the relative steric and electrostatic field contributions were similar to those of the mAbs, the DAT cocaine-binding site showed a preference for negatively charged ligands. Besides establishing molecular level insight into the interactions that govern cocaine binding specificity by biopolymers, the three-dimensional images obtained reflect the properties of the mAbs binding pockets and provide the initial information needed for the possible design of novel antibodies with properties optimized for immunotherapy.     

Comparative QSAR studies on peptide deformylase inhibitors

Comparative quantitative structure–activity relationship (QSAR) analyses of peptide deformylase (PDF) inhibitors were performed with a series of previously published (British Biotech Pharmaceuticals, Oxford, UK) reverse hydroxamate derivatives having antibacterial activity against Escherichia coli PDF, using 2D and 3D QSAR methods, comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), and hologram QSAR (HQSAR). Statistically reliable models with good predictive power were generated from all three methods (CoMFA r ² = 0.957, q ² = 0.569; CoMSIA r ² = 0.924, q ² = 0.520; HQSAR r ² = 0.860, q ² = 0.578). The predictive capability of these models was validated by a set of compounds that were not included in the training set. The models based on CoMFA and CoMSIA gave satisfactory predictive r ² values of 0.687 and 0.505, respectively. The model derived from the HQSAR method showed a low predictability of 0.178 for the test set. In this study, 3D prediction models showed better predictive power than 2D models for the test set. This might be because 3D information is more important in the case of datasets containing compounds with similar skeletons. Superimposition of CoMFA contour maps in the active site of the PDF crystal structure showed a meaningful correlation between receptor–ligand binding and biological activity. The final QSAR models, along with information gathered from 3D contour and 2D contribution maps, could be useful for the design of novel active inhibitors of PDF. Figure Superimposition of comparative molecular field analysis (CoMFA) contour plot in the active site of peptide deformylase (PDF)     

Convergent QSAR studies on a series of NK3 receptor antagonists for schizophrenia treatment

The dopamine hypothesis states that decreased dopaminergic neurotransmission reduces schizophrenia symptoms. Neurokinin-3 receptor (NK₃) antagonists reduce dopamine release and have shown positive effects in pre-clinical and clinical trials. We employed 2D and 3D-QSAR analysis on a series of 40 non-peptide NK₃ antagonists. Multivariate statistical analysis, PCA and HCA, were performed to rational training/test set splitting and PLS regression was employed to construct all QSAR models. We constructed one highly predictive CoMFA model (q^2?=?0.810 and r^2?=?0.929) and acceptable HQSAR and CoMSIA models (HQSAR q^2?=?0.644 and r^2?=?0.910; CoMSIA q^2?=?0.691, r^2?=?0.911). The three different techniques provided convergent physicochemical results. All models indicate cyclopropane, piperidine and di-chloro-phenyl ring attached to cyclopropane ring and also the amide group attached to the piperidine ring could play an important role in ligand–receptor interactions. These findings may contribute to develop potential NK₃ receptor antagonists for schizophrenia.     

Exploring the binding features of rimonabant analogues and acyclic CB<Subscript>1</Subscript> antagonists: docking studies and QSAR analysis

In order to elucidate the structural requirements for human CB₁ receptor antagonism, 78 antagonists belonging to five different chemical classes were selected from the literature and docked into the receptor binding site, built by homology modeling techniques. To further explore the structure-activity relationships within the considered chemical classes, a pharmacophore model and a QSAR analysis were developed. In a first step five alignments, one for each group of compounds were generated. All of them were then submitted to a MOE pharmacophore search in order to obtain a final pharmacophore model representative of the whole dataset which was used to elaborate the following 3D-QSAR analysis, by means of the CoMFA methodology. The results of these investigations are expected to be useful in the process of design and development of new potent CB₁ antagonists. Figure Compounds 1-78 are aligned into the putative CB1 receptor binding site. The three key features shared by all of them are reported in coloured spheres. The hydrophobic/aromatic ones are depicted in purple while the acceptor functions are coloured in blue.     

3D-QSAR studies on PU3 analogues by comparative molecular field analysis

A comparative molecular field analysis (CoMFA) of PU3 derivatives of Hsp90 (Heat shock protein 90) inhibitors has been performed to determine the factors contributing the corresponding activities. The energy minimized conformations were obtained by molecular mechanics using SYBYL package. The developed model, with r(2) value of 0.947, was verified by performing leave-one out (LOO) cross-validation, which showed q(2) value of 0.513. The calculated model not only elucidates the relationship between compound structures and biological activities but, more importantly, facilitates design of new Hsp90 inhibitors with calculated antiproliferative activity.     

Structure-activity relationship study of human liver microsomes-catalyzed hydrolysis rate of ester prodrugs of MENT by comparative molecular field analysis (CoMFA)

A series of MENT esters (3-71) was designed, prepared and tested to study the structure-activity relationship (SAR) of the hydrolysis rate with human liver microsomes of these prodrugs. Compounds were obtained covering a wide range of metabolic stability. The results are useful for the proper selection of prodrugs for different pharmaceutical formulations to deliver the potent and prostate-sparing androgen MENT. The MENT esters can especially be administered for male hormone replacement therapy and male contraception. Comparative molecular field analysis (CoMFA) was applied to a dataset of 28 esters, for which ED50 values could be obtained. The CoMFA model where the electrostatic and H-bond molecular fields were combined turned out to be most predictive. Despite the limited size of the dataset, CoMFA can help to rationalize the SAR of the ester hydrolysis rate of ester prodrugs of MENT.     

In silico modelling and molecular dynamics simulation studies of thiazolidine based PTP1B inhibitors

Protein tyrosine phosphatase 1B (PTP1B) has been identified as a negative regulator of insulin and leptin signalling pathway; hence, it can be considered as a new therapeutic target of intervention for the treatment of type 2 diabetes. Inhibition of this molecular target takes care of both diabetes and obesity, i.e. diabestiy. In order to get more information on identification and optimization of lead, pharmacophore modelling, atom-based 3D QSAR, docking and molecular dynamics studies were carried out on a set of ligands containing thiazolidine scaffold. A six-point pharmacophore model consisting of three hydrogen bond acceptor (A), one negative ionic (N) and two aromatic rings (R) with discrete geometries as pharmacophoric features were developed for a predictive 3D QSAR model. The probable binding conformation of the ligands within the active site was studied through molecular docking. The molecular interactions and the structural features responsible for PTP1B inhibition and selectivity were further supplemented by molecular dynamics simulation study for a time scale of 30 ns. The present investigation has identified some of the indispensible structural features of thiazolidine analogues which can further be explored to optimize PTP1B inhibitors.     

The conformation and activity relationship of benzofuran type of angiotensin II receptor antagonists

As a continuing effort to establish the structure and activity relationship in a benzofuran type of angiotensin II antagonist, we synthesized various regioisomers and performed a series of QSAR analyses. The conformational analyses of target isomers were carried out using molecular mechanics and fine-tuned using the information from the NMR NOE experiment. The conformations of compounds with a good binding activity are quite similar to that of DuP753, a prototype of AII antagonist, suggesting that these compounds also bind to the same site of AII receptor. We then studied the compounds with a varied length of the hydroxyl group bearing side chain to find out the optimum distance between the hydroxyl group and the imidazole ring. The CoMFA with these compounds gave acceptable statistical measures (cross-validated r2 and conventional r2 to be 0.881 and 0.974, respectively) and the map was well consistent with the previously proposed pharmacophore.     

Structure-based 3D-QSAR studies on thiazoles as 5-HT<Subscript>3</Subscript> receptor antagonists

Structure-based 3D-QSAR studies were performed on 20 thiazoles against their binding affinities to the 5-HT₃ receptor with comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The thiazoles were initially docked into the binding pocket of a human 5-HT_3A receptor homology model, constructed on the basis of the crystal structure of the snail acetylcholine binding protein (AChBP), using the GOLD program. The docked conformations were then extracted and used to build the 3D-QSAR models, with cross-validated values 0.785 and 0.744 for CoMFA and CoMSIA, respectively. An additional five molecules were used to validate the models further, giving satisfactory predictive values of 0.582 and 0.804 for CoMFA and CoMSIA, respectively. The results would be helpful for the discovery of new potent and selective 5-HT₃ receptor antagonists.      

Dissecting the molecular mechanism of drosophila odorant receptors through activity modeling and comparative analysis

To gain insight into the molecular mechanism of odorant receptors (ORs) in Drosophila species, we developed a Quantitative Structure Activity Relationship (QSAR) model that predicts experimentally measured electrophysiological activities between 24 D. melanogaster ORs and 108 odorants. Although the model is limited by the tested odorants,analyzing the model allowed dissection of specific topological and chemical properties necessary for an odorant to elicit excitatory or inhibitory receptor response. Linear odorants with five to eight nonhydrogen atoms at the main chain and hydrogen‐bond acceptor and/or hydrogen‐bond donor at its ends were found to stimulate strong excitatory response. A comparative sequence analysis of 90 ORs in 15 orthologous groups identified 15 putative specificity‐determining residues (SDRs) and 15 globally conserved residues that we postulate as functionally key residues. Mapping to a model of secondary structure resulted in 14 out of 30 key residues locating to the transmembrane (TM) domains. Twelve residues, including six SDRs and six conserved residues, are located at the extracellular halves of the TM domains. Combining the evidence from the QSAR modeling and the comparative sequence analysis, we hypothesize that the Drosophila ORs accept odorants into a binding pocket located on the extracellular halves of its TM domains. The QSAR modeling suggests that the binding pocket is around 15 Å in depth and about 6 Å in width. Twelve mainly polar or charged key residues, both SDRs and conserved, are located inthis pocket and postulated to distinguish docked odorants via primarily geometry fitting and hydrogen‐bond interaction. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Homology modelling and molecular dynamics simulations: comparative studies of human aquaporin-1

The structures of the mammalian water transport protein Aqp1 and of its bacterial homologue GlpF enables us to test whether homology models can be used to explore relationships between structure, dynamics and function in mammalian transport proteins. Molecular dynamics simulations (totalling almost 40 ns) were performed starting from: the X-ray structure of Aqp1; a homology model of Aqp1 based on the GlpF structure; and intermediate resolution structures of Aqp1 derived from electron microscopy. Comparisons of protein RMSDs vs. time suggest that the homology models are of comparable conformational stability to the X-ray structure, whereas the intermediate resolution structures exhibit significant conformation drift. For simulations based on the X-ray structure and on homology models, the flexibility profile vs. residue number correlates well with the crystallographic B-values for each residue. In the simulations based on intermediate resolution structures, mobility of the highly conserved NPA loops is substantially higher than in the simulations based on the X-ray structure or the homology models. Pore radius profiles remained relatively constant in the X-ray and homology model simulations but showed substantial fluctuations (reflecting the higher NPA loop mobility) in the intermediate resolution simulations. The orientation of the dipoles of water molecules within the pore is of key importance in maintaining low proton permeability through Aqp1. This property seems to be quite robust to the starting model used in the simulation. These simulations suggest that homology models based on bacterial homologues may be used to derive functionally relevant information on the structural dynamics of mammalian transport proteins.     

Models of Monoamine Oxidase A and B Active Sites Obtained by using 3D QSAR with CoMFA Analysis

Abstract

The monoamine oxidase catalyses the oxidative deamination of neuroactive amines. This enzyme exists in two forms A and B, which differ by substrates preference and inhibitors specificity. Investigation of the structures of these enzymes and design new selective inhibitors are of greatly interesting since MAO A inhibitors are used in therapeutic practice as antidepressants and MAO B inhibitors – in the treatment Parkinson's diseases. The three dimension structures of monoamine oxidases are still unknown. Therefore, one of the most perspective approach to define significant features of structure active site is method based on analysis of structure-activity relationship (3D QSAR) with comparison of molecular fields analysis (CoMFA) allowing to get the spatial distribution of important properties affecting the activity.

In present study we investigate the structures of active sites MAO A and B using 16 pyrazinocarbazole derivatives in variant conformation. Majority of pyrazinocarbazole derivatives have a rigit conformation, but three of those is sufficiently flexible. The latters can be in two conformation types: long molecules (substitution accommodate along axis of main structure) and short molecules (substitution accommodate at acute angle about of main structure). Several 3D QSAR and CoMFA models of MAO A and B active sites were design for data sets containing various types of flexible molecules conformation. All obtained models are statistical reliable and have sufficient predictive power for tested compound tetrindole. The best MAO A model that include two flexible molecules in long conformations was obtained, and the longest one of those in short conformation. In contrast, for MAO B model containing all flexible molecules in the short conformations is more preferred.

On the basis of obtained data the schematic models of MAO A and B active sites structures are proposed. According to these models MAO A active site have the narrow long cavity that accommodate long molecules, while MAO B active site is broader and shorter.     

Molecular modelling studies on d-annulated benzazepinones as VEGF-R2 kinase inhibitors using docking and 3D-QSAR

Chemotypes comprising the d-annulated 1,3-dihydro-2H-1-benzazepin-2-one scaffold derived from the paullone structure were found to be potent vascular endothelial growth factor receptor 2 (VEGF-R2) kinase inhibitors. Three-dimensional quantitative structure-activity relationship (3D-QSAR) and docking studies were performed on a series of d-annulated benzazepinones with VEGF-R2 kinase inhibition activities. The comparative molecular field analysis and comparative molecular similarity indices analysis models using 32 molecules in the training set gave r²_cv values of 0.811 and 0.769, r² values of 0.962 and 0.953, respectively. 3D contour maps generated from the two models revealed that the electron-withdrawing groups at R₁ and the bulky, electron-withdrawing as well as hydrogen bond donor groups at R₂ position are favourable; the bulky, hydrogen bond acceptor substituent at R₃ and the minor groups at R₄ position may benefit the potency. We have designed a series of novel VEGF-R2 inhibitors by utilizing the SAR results revealed in the present study, which were predicted with excellent potencies in the developed models. The results may aid in designing of potential VEGF-R2 inhibitors with better activities.     

A new approach to angiotensin antagonists: Methylation of the tyrosine hydroxyl in angiotensin II

[Sar¹, Tyr(Me)⁴] angiotensin II, synthesized by the solid phase method and purified by ion-exchange chromatography and reversed-phase HPLC, was found to inhibit the contractile response to angiotensin II in the rat isolated uterus and inhibit the pressor response to angiotensin II in the vagotomized ganglion-blocked rat. In the rat isolated uterus Schild plots gave a pA₂ of 8.1, and a slope of 0.9 indicative of competitive inhibition. In the rat pressor assay, infusion of the analogue at a rate of 500ng/kg/min caused a parallel displacement of the dose-response curve to ANG II to the right. In contrast, the classical angiotensin inhibitor [Sar¹, Ile⁸] ANG II appeared to demonstrate non-competitive inhibition in both the rat isolated uterus and the pressor assays. The phenolic hydroxyl or phenoxide anion of Tyr⁴ in angiotensin II appears to be critical for the activation of angiotensin receptors in smooth muscle. Alkylation of the tyrosine residue in angiotensin analogues provides a new route for the synthesis of potent competitive antagonists of angiotensin.