Comparative molecular field analysis (CoMFA) of new herbicidal sulfonylurea compounds

Ｓｕｌｆｏｎｙｌｕｒｅａｈｅｒｂｉｃｉｄｅｓｈａｖｅｂｅｅｎｕｓｅｄｆｏｒｓｅｖｅｒａｌｄｅｃａｄｅｓ,ｗｈｉｃｈｐｏｓｓｅｓｓｕｌｔｒａｌｏｗｄｏｓａｇｅ,ｌｏｗｍａｍｍａｌｉａｎｔｏｘｉｃｉｔｙａｎｄｄｅｓｉｒａｂｌｅｅｎｖｉｒｏｎｍｅｎｔａｌｐｒｏｐｅｒｔｉｅｓ.ＴｈｅｙｋｉｌｌｗｅｅｄｓｂｙｉｎｈｉｂｉｔｉｎｇｔｈｅａｃｔｉｖｉｔｙｏｆＡＬＳｅｎｚｙｍｅｏｆｔｈｅｗｅｅｄｓａｎｄｄｅｓｔｒｏｙｉｎｇｔｈｅｓｙｎｔｈｅｓｉｓｏｆｂｒａｎｃｈｃｈａｉｎａｍｉｎｏａｃｉｄ.Ｓｉｎｃｅｔｈｅｅｓｔａ…     

Comparative molecular field analysis (CoMFA) for phosphodiesterase (PDE) IV inhibitors

A comparative molecular field analysis (CoMFA) for phosphodiesterase (PDE) IV inhibitors has been performed to correlate their chemical structures with their observed biological activity. In this study, CoMFA model based on docking mode for active site of PDE IV can describe the relative change in magnitude of the steric and electrostatic fields as a function of the compounds. Pyridine N-oxide and pyridine group of each compound are aligned toward the metal ion in S2-sub pocket of PDE IV. The study provided a statistically valid model with good correlation and predictive power, and consequently we identified some key features that may be used to design new derivatives.     

Comparative molecular field analysis (CoMFA) model using a large diverse set of natural, synthetic and environmental chemicals for binding to the androgen receptor

A large number of natural, synthetic and environmental chemicals are capable of disrupting the endocrine systems of experimental animals, wildlife and humans. These so-called endocrine disrupting chemicals (EDCs), some mimic the functions of the endogenous androgens, have become a concern to the public health. Androgens play an important role in many physiological processes, including the development and maintenance of male sexual characteristics. A common mechanism for androgen to produce both normal and adverse effects is binding to the androgen receptor (AR). In this study, we used Comparative Molecular Field Analysis (CoMFA), a three-dimensional quantitative structure-activity relationship (3D-QSAR) technique, to examine AR-ligand binding affinities. A CoMFA model with r2 = 0.902 and q2 = 0.571 was developed using a large training data set containing 146 structurally diverse natural, synthetic, and environmental chemicals with a 10(6)-fold range of relative binding affinity (RBA). By comparing the binding characteristics derived from the CoMFA contour map with these observed in a human AR crystal structure, we found that the steric and electrostatic properties encoded in this training data set are necessary and sufficient to describe the RBA of AR ligands. Finally, the CoMFA model was challenged with an external test data set; the predicted results were close to the actual values with average difference of 0.637 logRBA. This study demonstrates the utility of this CoMFA model for real-world use in predicting the AR binding affinities of structurally diverse chemicals over a wide RBA range.     

Prediction of LUMO energy and rate constant by comparative molecular field analysis (CoMFA)

We used the comparative molecular field analysis (CoMFA) method to correlate the rate constant (log k) for the S_N2 reaction of benzyl benzenesulfonates and p-methoxybenzylamines. Molecular fields calculated with a C⁺ probe produced a good correlation with a small standard deviation and a high correlation coefficient with cross validation. This study demonstrated that CoMFA is an excellent method in predicting the physicochemical properties of the molecule such as LUMO energy and rate constants. © 1995 by John Wiley & Sons, Inc.     

Multilinear regression and comparative molecular field analysis (CoMFA) of azo dye-fiber affinities. 2. Inclusion of solution-phase molecular orbital descriptors

For a data set with 30 direct azo dyes taken from literature, quantitative structure-activity relationship (QSAR) analyses have been performed to model the affinity of the dye molecules for the cellulose fiber. The electronic structure of the compounds was characterized using quantum chemical gas-phase (AM1) and continuum-solvation molecular orbital parameters. As regards the solution phase, COSMO appears to be better suited than SM2 in quantifying relative trends of the aqueous solvation energy. For the dye-fiber affinity, the leave-one-out prediction capability of multilinear regression equations is superior to CoMFA, with predictive squared correlation coefficients ranging from 0.63 (pure CoMFA) to 0.89. At the same time, solution-phase CoMFA is superior to previously derived AM1-based CoMFA models. As a general trend, the dye-fiber affinity increases with increasing electron donor capacity that corresponds to an increasing hydrogen bond acceptor strength of the azo dyes. The discussion includes the consideration of structural features that are likely to be involved in dye-fiber and dye-dye hydrogen bonding interactions, and possible links between CoMFA electrostatic results and the atomic charge distribution of the compounds.     

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) of thiazolone derivatives as hepatitis C virus NS5B polymerase allosteric inhibitors

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models for a series of thiazolone derivatives as novel inhibitors bound to the allosteric site of hepatitis C virus (HCV) NS5B polymerase were developed based on CoMFA and CoMSIA analyses. Two different conformations of the template molecule and the combinations of different CoMSIA field/fields were considered to build predictive CoMFA and CoMSIA models. The CoMFA and CoMSIA models with best predictive ability were obtained by the use of the template conformation from X-ray crystal structures. The best CoMFA and CoMSIA models gave q (2) values of 0.621 and 0.685, and r (2) values of 0.950 and 0.940, respectively for the 51 compounds in the training set. The predictive ability of the two models was also validated by using a test set of 16 compounds which gave r (pred) (2) values of 0.685 and 0.822, respectively. The information obtained from the CoMFA and CoMSIA 3D contour maps enables the interpretation of their structure-activity relationship and was also used to the design of several new inhibitors with improved activity.     

Comparative molecular field analysis of non-steroidal aromatase inhibitors related to fadrozole

Summary A series of non-steroidal inhibitors of aromatase, structurally related to fadrozole (2), was investigated with the aim of developing a 3D QSAR model using the Comparative Molecular Field Analysis (CoMFA) technique. The alignment of the molecules was performed following two approaches (atom-by-atom and field fit), both starting from an initial hypothesis of superimposition of fadrozole to a steroidal inhibitor (3). From a number of CoMFA models built with different characteristics, one was recognized as the most statistically relevant; this one is discussed in detail. The features of the 3D QSAR model are consistent with those of other 3D and QSAR models of aromatase and its inhibitors.     

Use of the hydrogen bond potential function in a comparative molecular field analysis (CoMFA) on a set of benzodiazepines

Summary The results of the GRID-Comparative Molecular Field Analysis (CoMFA) were compared with those of the SYBYL-CoMFA in a study of benzodiazepines. The results demonstrate that the hydrogen bonding function using the GRID H₂O probe in a CoMFA can successfully describe the hydrophobic effects of substituents without any bias or preconcept of their effects in the development.     

GRID formalism for the comparative molecular surface analysis: application to the CoMFA benchmark steroids, azo dyes, and HEPT derivatives

Shape analysis is a powerful tool in chemistry and drug design, and molecular surface defines shape in the molecular scale. In the current publication we presented a novel formalism for the comparative molecular surface analysis (s-CoMSA). The method enables both quantitative modeling of 3D-QSAR and finding possible pharmacophoric sites. The method provides very predictive models for the CBG activity of the benchmark steroid series, tinctorial properties of the heterocyclic azo dyes and anti-HIV activity of the HEPT series.     

A comparative molecular field analysis (CoMFA) study using semiempirical, density functional, ab initio methods and pharmacophore derivation using DISCOtech on sigma 1 ligands

The Comparative Molecular Field Analysis (CoMFA) was developed to investigate a three-dimensional quantitative structure activity relationship (3D-QSAR) model of ligands for the sigma 1 receptor. The starting geometry of sigma-1 receptor ligands was obtained from the Tripos force field minimizations and conformations were decided from DISCOtech using the SYBYL 6.8. program. The structures of 48 molecules were fully optimized at the ab initio HF/3-21G* and semiempirical AM1 calculations using GAUSSIAN 98. The electrostatic charges were calculated using several methods such as semiempirical AM1, density functional B3LYP/3-21G*, and ab initio HF/3-21G*, MP2/3-21G* calculations within GAUSSIAN 98. Using the optimized geometries, the CoMFA results derived from the HF/3-21G method were better than those from AM1. The best CoMFA was obtained from HF/3-21G* optimized geometry and charges (R2 = 0.977). Using the optimized geometries, the CoMFA results derived from the HF/3-21G methods were better than those from AM1 calculations. The training set of 43 molecules gave higher R2 (0.989-0.977) from HF/3-21G* optimized geometries than R2 (0.966-0.911) values from AM1 optimized geometries. The test set of five molecules also suggested that HF/3-21G* optimized geometries produced good CoMFA models to predict bioactivity of sigma 1 receptor ligands but AM1 optimized geometries failed to predict reasonable bioactivity of sigma 1 receptor ligands using different calculations for atomic charges.     

Comparative molecular field analysis of protein tyrosine phosphatase low-molecular-weight substrates

A set of 25 low-molecular-weight substrates for protein tyrosine phosphatase PTP1 has been investigated by comparative molecular field analysis (CoMFA). The quality of the model was assessed by cross-validations, predictions for a test set of substrates, bootstrapping, randomization of biological data, grid characteristics tests, and comparison with the X-ray-determined binding site of a closely related enzyme. The CoMFA study provides new insight into the steric and electrostatic factors influencing binding around the active site in protein tyrosine phosphatase and complement existing information available from X-ray data.     

Comparative molecular field analysis and comparative molecular similarity index analysis studies on 1H NMR chemical shift of NH group of diaryl triazene derivatives

Comparative molecular field analysis (CoMFA), comparative molecular field analysis region focusing (CoMFA‐RF) for optimizing the region for the final partial least square analysis, and comparative molecular similarity indices analysis (CoMSIA) methods were employed to develop three‐dimensional quantitative structure–activity relationship (3D‐QSAR) models of ¹H NMR chemical shift of NH proton of diaryl triazene derivatives. The best orientation was searched by all‐orientation search (AOS) strategy to minimize the effect of the initial orientation of the structures. The predictive abilities of CoMFA‐RF and CoMSIA models were determined using a test set of ten compounds affording predictive correlation coefficients of 0.721 and 0.754, respectively, indicating good predictive power. For further model validation, cross validation (leave one out), progressive scrambling, and bootstrapping were also applied. The accuracy and speed of obtained 3D‐QSAR models for the prediction of ¹H NMR chemical shifts of NH group of diaryl triazene derivatives were greater compared to some computational well‐known procedures. Copyright © 2013 John Wiley & Sons, Ltd.     

Comparative analysis of binding energy of chymostatin with human cathepsin A and its homologous proteins by molecular orbital calculation

Cathepsin A is a mammalian lysosomal enzyme that catalyzes the hydrolysis of the carboxy-terminal amino acids of polypeptides and also regulates beta-galactosidase and neuraminidase-1 activities through the formation of a multienzymic complex in lysosomes. Human cathepsin A (hCathA), yeast carboxypeptidase (CPY), and wheat carboxypeptidase II (CPW) belong to the alpha/beta-hydrolase fold family. They have structurally similar active-site clefts, but there are small differences in the amino acid residues comprising their active sites that might determine the substrate specificity and sensitivity to microbial inhibitors including chymostatin. To examine the selectivity and binding mechanism of chymostatin as to hCathA, CPY, and CPW at the atomic level, we analyzed the interaction energy between chymostatin and each protein quantitatively by semiempirical molecular orbital calculation AM1 with the continuum solvent model. We predicted the electrostatic repulsion between the P3 cyclic arginine residue of the inhibitor and the Arg344 in the S3 active subsite of hCathA. Genetic conversion of Arg344 of the wild-type hCathA to Ile also caused an increase in its sensitivity to chymostatin, which was correlated with the decrease in the interaction energy calculated with the molecular orbital method. The present results suggest that such molecular calculation should be useful for evaluating the interactions between ligands, including inhibitors and homologous enzymes, in their docking models.     

Comparative molecular field analysis and energy interaction studies of thrombin-inhibitor complexes

A Comparative Molecular Field Analysis (CoMFA) and an interaction energy-based method were applied on a database holding the 3D structures of 29 thrombin-inhibitor complexes. Several parameters were optimized in both methods in order to obtain the best correlation between theoretical and experimentally determined binding (Ki) data. CoMFA, which only uses the information of the inhibitors, performed best (r = 0.99, q2 = 0.46, N = 29) when HF 6-31G charges were used in combination with a pharmacophore-based alignment. Inclusion of hydrophobic fields did not lead to improvements. The interaction energy-based approach uses the information of the whole thrombin-inhibitor complex. A statistically significant correlation (r = 0.74, N = 14) could only be obtained for a subset of the database containing the high resolution structures. Geometry optimization of the ligand only in combination with downscaled electrostatics performed best.     

Branched non-ionic oligo-oxyethylene Y-amphiphiles Effect of molecular geometry on the micellar shape 1

A homologous series of branched, non-ionic surfactants with the general formula C_nG(E_mM)₂, where C_n denotes an alkyl chain, G = glycerol and E_mM = oligo-oxyethylene mono-methyl ether, has been prepared from alkyl bromides (n = 10-16) and several monodisperse 1,3-di(methoxyoligo-oxyethylene) ethers of glycerol (m = 3-5). The branched hydrophilic chain is introduced to modify the interfacial area compared to corresponding linear oligo-oxyethylene surfactants (I-amphiphiles) without essentially changing the hydrophilic-lipophilic balance. The phase behaviour of these Y-surfactants in aqueous solution reveals that according to established packing models the branched hydrophilic group strongly stabilizes the cubic and hexagonal mesophases, while a lamellar phase is not observed.     

Structural analysis of the protein phosphatase 1 docking motif: molecular description of binding specificities identifies interacting proteins

The interplay between kinases and phosphatases represents a fundamental regulatory mechanism in biological systems. Being less numerous than kinases, phosphatases increase their diversity by the acquisition of a variety of binding partners, thereby forming a large number of holoenzymes. Proteins interacting with protein phosphatase 1 (PP1) often bind via a so-called docking motif to regulate its enzymatic activity, substrate specificity, and subcellular localization. Here, we systematically determined structural elements that mediate the binding specificity of PP1 interacting proteins, and propose a refined consensus sequence for high-affinity PP1 ligands. Applying this pattern to database searches, we predicted and experimentally confirmed several previously unknown PP1 interactors. Thus, the suggested PP1 docking motif enables a highly specific prediction of PP1 binding partners, thereby facilitating the genome-wide identification of PP1 interactors.     

Molecular modeling of the intestinal bile acid carrier: A comparative molecular field analysis study

A structure–binding activity relationship for the intestinal bile acidtransporter has been developed using data from a series of bile acid analogsin a comparative molecular field analysis (CoMFA). The studied compoundsconsisted of a series of bile acid–peptide conjugates, withmodifications at the 24 position of the cholic acid sterol nucleus, andcompounds with slight modifications at the 3, 7, and 12 positions. For theCoMFA study, these compounds were divided into a training set and a test set,comprising 25 and 5 molecules, respectively. The best three-dimensionalquantitative structure–activity relationship model found rationalizesthe steric and electrostatic factors which modulate affinity to the bile acidcarrier with a cross-validated, conventional and predictive r²of 0.63, 0.96, and 0.69, respectively, indicating a good predictive model forcarrier affinity. Binding is facilitated by positioning an electronegativemoiety at the 24–27 position, and also by steric bulk at the end of theside chain. The model suggests substitutions at positions 3, 7, 12, and 24that could lead to new substrates with reasonable affinity for the carrier.