Studies on the interaction of extended terpyridyl and triazine metal complexes with DNA

The interaction with duplex DNA of a series of 12 previously synthesised ruthenium (II) complexes incorporating extended terpyridines or analogous ligands is reported. The interaction was investigated using absorption spectroscopy and also, when possible, luminescence spectroscopy. Binding parameters for the interactions were estimated using established models. It was found that binding affinities are relatively large and their magnitude is somewhat dependent on the metal centre attached to the ligand. Viscosity studies reveal that all the complexes interact with DNA through groove binding and not intercalation. The variation in binding affinities of these complexes is rationalised by a consideration of electrostatic and steric factors.     

Three-dimensional quantitative structure–farnesyltransferase inhibition analysis for some diaminobenzophenones

A 3D-QSAR investigation of 95 diaminobenzophenone yeast farnesyltransferase (FT) inhibitors selected from the work of Schlitzer et al. showed that steric, electrostatic, and hydrophobic properties play key roles in the bioactivity of the series. A CoMFA/CoMSIA combined model using the steric and electrostatic fields of CoMFA together with the hydrophobic field of CoMSIA showed significant improvement in prediction compared with the CoMFA steric and electrostatic fields model. The similarity of the 3D-QSAR field maps for yeast FT inhibition activity (from this work) and for antimalarial activity data (from previous work) and the correlation between those activities are discussed.     

Calculations of antibody-antigen interactions: microscopic and semi-microscopic evaluation of the free energies of binding of phosphorylcholine analogs to McPC603.

The study of antibody-antigen interactions should greatly benefit from the development of quantitative models for the evaluation of binding free energies in proteins. The present work addresses this challenge by considering the test case of the binding free energies of phosphorylcholine analogs to the murine myeloma protein McPC603. This includes the evaluation of the differential binding energy as well as the absolute binding energies and their corresponding electrostatic contributions. Four different approaches are examined: the Protein Dipoles Langevin Dipoles (PDLD) method, the semi-microscopic PDLD (PDLD/S) method, a free energy perturbation (FEP) method based on an adiabatic charging procedure and a linear response approximation that accelerates the FEP calculation. The PDLD electrostatic calculations are augmented by estimates of the relevant hydrophobic and steric contributions. The determination of the hydrophobic energy involves an approach which considers the modification of the effective surface area of the solute by local field effects. The steric contributions are analyzed in terms of the corresponding reorganization energies. This treatment, which considers the protein as a harmonic system, views the steric forces as the restoring forces for the electrostatic interactions. The FEP method is found to give unreliable results with regular cut-off radii and starts to give quantitative results only in very expensive treatment with very large cut-off radii. The PDLD and PDLD/S methods are much faster than the FEP approach and give reasonable results for both the relative and absolute binding energies. The speed and simplicity of the PDLD/S method make it an effective strategy for interactive docking studies and indeed such an option is incorporated in the program MOLARIS. A component analysis of the different energy contributions of the FEP treatment and a similar PDLD analysis indicate that electrostatic effects provide the largest contribution to the differential binding energy, while the hydrophobic and steric contributions are much smaller. This finding lends further support to the idea that electrostatic interactions play a major role in determining the antigen specificity of McPC603.     

High-throughput screening of ecdysone agonists using a reporter gene assay followed by 3-D QSAR analysis of the molting hormonal activity

In this study, 172 diacylhydrazine analogs were examined for their ability to activate an ecdysone (molting hormone)-dependent reporter gene in a silkworm (Bombyx mori) cell-based high-throughput screening assay. The measured EC(50) values (concentration required to cause an effect in 50% of the cells) were used to construct a 3-D QSAR model that describes the ecdysone agonist activities of the diacylhydrazine analogs. Of these compounds, 14 exhibited no activity and were excluded from the 3-D QSAR analysis. The resulting equation described approximately 74% of the activity for 158 compounds. The final equation consisted of 42% electrostatic and 58% steric effects (r(2) = 0.74 and q(2) = 0.45). Comparative molecular field analysis (CoMFA) was used to visualize the steric and electrostatic potential fields that were favorable and unfavorable for biological activity. Of particular interest was the observation that the hydrophobic parameter (logP) was not necessary for describing the observed activities, although previous studies have cited the importance of hydrophobic parameters in both classical and 3-D QSAR analyses of these compounds. Modeling studies of the B. mori ecdysone receptor supported the observed physicochemical parameters required for activity reported by the CoMFA models. Comparison of the present analysis with those performed using other lepidopteran assay systems evidenced a high degree of correlation (r(2) = 0.81 for a Sf-9 cell-based assay and r(2) = 0.89 for a Chilo suppressalis integument-based assay), indicating that it is valid to compare the results generated with the B. mori cell-based system to those generated with previous lepidopteran assays. This novel assay system is amendable to a high-throughput screening format and should greatly increase our ability to discover novel agonists of molting hormone (ecdysone) activity.     

Three-Dimension Quantitative Structure-Activity Relationship Analysis of Some Cinnamamides Using Comparative Molecular Similarity Indices Analysis (CoMSIA)

Comparative molecular similarity indices analysis (CoMSIA) has been applied to a data set of cinnamamides. Five different properties: steric, electrostatic, hydrophobic, H-bond donor and H-bond acceptor, assumed to cover the major contributions to ligand binding, were used to generate the 3D-QSAR model. A significant cross-validated correlation coefficient q ² (0.691) was obtained, indicating the statistical significance of the model for untested compounds of this class. The model was used to predict the anticonvulsant activities of five test-set compounds, and the predicted values were in good agreement with the experimental results. Moreover, from the contour maps, the key features vital to ligand binding have been identified, which are important for us to trace the important properties and gain insight into the potential mechanisms of intermolecular interactions between ligand and receptor.Electronic Supplementary Material available.     

Quantitative structure-activity relationships in a set of thiazolidin-4-ones acting as H1-histamine antagonists

Abstract

A series of 2-(3- and 4-substituted phenyl)-3-[3-(N,N-dimethylamino)propyl]-1,3-thiazolidin-4-ones acting as H₁-antihistaminics was investigated with a combined Hansch-CoMFA approach. The substituents at the 3- and 4-positions of the phenyl ring have been described through steric, electronic and hydrophobic parameters and correlated with pA₂ values. The obtained quantitative models suggest that affinity to the receptor is promoted by hydrophobic and small 4-substituents and by 3- and 4-substituents generating a positive electrostatic potential towards a complementary receptor region.     

Isolation,biological evaluation and 3D-QSAR studies of insecticidal/narcotic sesquiterpene polyol esters

For the first time, a set of (43) natural sesquiterpene polyol esters isolated from the root bark of Celastrus angulatus Maxim and Euonymus japonicus Thunb were subjected to 3D-QSAR comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) studies, with the aim of proposing novel sesquiterpene-based compounds with optimal narcotic or insecticidal activities. The established 3D-QSAR models exhibit reasonable statistical quality and prediction capabilities, with internal cross-validated Q ² values of ∼0.5 and external predicted R ² values of >0.9, respectively. The relative contributions of the steric/electrostatic fields of the 3D-QSAR models show that the electronic effect governs the narcotic activities of the molecules, but the hybrid effect of the electrostatic and hydrophobic interactions is more influential in the insecticidal activities of the compounds. These findings may have valuable implications for the development of novel natural insecticides.     

Pharmacophore modeling, virtual screening and 3D-QSAR studies of 5-tetrahydroquinolinylidine aminoguanidine derivatives as sodium hydrogen exchanger inhibitors

Sodium hydrogen exchanger (SHE) inhibitor is one of the most important targets in treatment of myocardial ischemia. In the course of our research into new types of non-acylguanidine, SHE inhibitory activities of 5-tetrahydroquinolinylidine aminoguanidine derivatives were used to build pharmacophore and 3D-QSAR models. Genetic Algorithm Similarity Program (GASP) was used to derive a 3D pharmacophore model which was used in effective alignment of data set. Eight molecules were selected on the basis of structure diversity to build 10 different pharmacophore models. Model 1 was considered as the best model as it has highest fitness score compared to other nine models. The obtained model contained two acceptor sites, two donor atoms and one hydrophobic region. Pharmacophore modeling was followed by substructure searching and virtual screening. The best CoMFA model, representing steric and electrostatic fields, obtained for 30 training set molecules was statistically significant with cross-validated coefficient (q(2)) of 0.673 and conventional coefficient (r(2)) of 0.988. In addition to steric and electrostatic fields observed in CoMFA, CoMSIA also represents hydrophobic, hydrogen bond donor and hydrogen bond acceptor fields. CoMSIA model was also significant with cross-validated coefficient (q(2)) and conventional coefficient (r(2)) of 0.636 and 0.986, respectively. Both models were validated by an external test set of eight compounds and gave satisfactory prediction (r(pred)(2)) of 0.772 and 0.701 for CoMFA and CoMSIA models, respectively. This pharmacophore based 3D-QSAR approach provides significant insights that can be used to design novel, potent and selective SHE inhibitors.     

Discriminating D1 and D2 agonists with a hydrophobic similarity index

Currently, methods for calculating molecular similarity indices have been developed for comparing steric, charge density, and molecular electrostatic potential (MEP) properties. Much of the existing technology may, however, be applied to the quantitative comparison of molecular hydrophobicities. In this article we present an empirical hydrophobic similarity index. We utilize atomic hydrophobic parameters derived from a quantum mechanical semiempirical wavefunction. Hydrophobicity at points on a grid is computed with a recently introduced “molecular lipophilicity potential”. The overlap of pairs of molecules is calculated with the metric introduced by Carbó. This approach is applied to a case in which steric and electrostatic criteria have already been shown to be inadequate in rationalizing selectivity, namely, requirements for recognition at the dopamine D1 and D2 receptors. We demonstrate that, for a set of dopamine agonists, D1 ligands show higher similarity in this property than D2 analogs. This indicator of similarity is more successful at accounting for D1 selectivity than previous methods.     

Docking by least-squares fitting of molecular surface patterns.

Molecular surfaces are fitted to each other by a new solution to the problem of docking a ligand into the active site of a protein molecule. The procedure constructs patterns of points on the surfaces and superimposes them upon each other using a least-squares best-fit algorithm. This brings the surfaces into contact and provides a direct measure of their local complementarity. The search over the ligand surface produces a large number of dockings, of which a small fraction having the best complementarity and the least steric hindrance are evaluated for electrostatic interaction energy. When applied to molecules taken from crystallographically observed complexes, this procedure consistently assigns the lowest electrostatic energies to correct dockings. On independently determined structures, the ability of the method to discern correct dockings depends on how much conformational difference there is between the free and complexed forms of the molecules. The procedure is found to be fast enough on contemporary workstation computers to permit many conformations to be considered, and tolerant enough to make rather coarse bond dihedral sampling a practicable way to overcome the problem of structural flexibility.     

Molecular modeling studies,synthesis and biological evaluation of dabigatran analogues as thrombin inhibitors

In this work, 48 thrombin inhibitors based on the structural scaffold of dabigatran were analyzed using a combination of molecular modeling techniques. We generated three-dimensional quantitative structure–activity relationship (3D-QSAR) models based on three alignments for both comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) to highlight the structural requirements for thrombin protein inhibition. In addition to the 3D-QSAR study, Topomer CoMFA model also was established with a higher leave-one-out cross-validation q² and a non-cross-validation r², which suggest that the three models have good predictive ability. The results indicated that the steric, hydrophobic and electrostatic fields play key roles in QSAR model. Furthermore, we employed molecular docking and re-docking simulation explored the binding relationship of the ligand and the receptor protein in detail. Molecular docking simulations identified several key interactions that were also indicated through 3D-QSAR analysis. On the basis of the obtained results, two compounds were designed and predicted by three models, the biological evaluation in vitro (IC₅₀) demonstrated that these molecular models were effective for the development of novel potent thrombin inhibitors.     

Complexes of 3,4-dihydroxyphenyl derivatives. VII. Mixed ligand complexes of L-dopa and related compounds

The stability constants of the mixed ligand complexes of manganese(II), cobalt(II), nickel(II), copper(II) and zinc(II) ions with L-dopa, dopamine, L-adrenaline and L-noradrenaline as ligand A and L-alanine, L-histidine, glycylglycine and ATP as ligand B were determined pH-metrically at 25 °C and an ionic strength of 0.2 mol/dm³ (KCl). From visible spectral studies, conclusions were drawn on the bonding properties of the ambidentate ligands in the complexes.It was found that, even in its mixed ligand complexes, the ambidentate nature of dopa permits its coordination via the amino acid side-chain (N, O) at lower pH, and via the ortho phenolic hydroxy groups (O, O) at higher pH. At the same time, the complex-forming tendency of the ethanolamine side-chain in the catecholamines plays a minor role in the mixed ligand complexes.The stabilities of the mixed ligand complexes were interpreted by consideration of the differences in the stepwise stability constants of the parent complexes, the back-coordination, the charge neutralization, and the electrostatic and hydrophobic interactions between the ligand pairs.     

Computer-aided design, synthesis and biological assay of p-methylsulfonamido phenylethylamine analogues

Class III antiarrhythmic agents selectively delay the effective refractory period (ERP) and increase the transmembrance action potential duration (APD). Based on our previous studies, a set of 17 methylsulfonamido phenylethylamine analogues were investigated by 3D-QSAR techniques of CoMFA and CoMSIA. The 3D-QSAR models proved a good predictive ability, and could describe the steric, electrostatic and hydrophobic requirements for recognition forces of the receptor site. According to the clues provided by this 3D-QSAR analysis, we designed and synthesized a series of new analogues of methanesulfonamido phenylethylamine (VIa-i). Pharmacological assay indicated that the effective concentrations of delaying the functional refractory period (FRP) 10ms of these new compounds have a good correlation with the 3D-QSAR predicted values. It is remarkable that the maximal percent change of delaying FRP in microM of compound VIc is much higher than that of dofetilide. The results showed that the 3D-QSAR models are reliable.     

Docking and 3D QSAR study of thiourea analogs as potent inhibitors of influenza virus neuraminidase

Surflex-Dock was applied to study interactions between 30 thiourea analogs and neuraminidase (NA). The docking results showed that hydrogen bonding and electrostatic interactions were highly correlated with the activities of neuraminidase inhibitors (NIs), followed by hydrophobic and steric factors. Moreover, there was a strong correlation between the predicted binding affinity (total score) and experimental pIC₅₀ (correlation coefficient r = 0.870; P < 0.0001). A three dimensional holographic vector of atomic interaction field (3D-HoVAIF) was employed to construct a QSAR model. The r ², q ² and r ² _test values of the optimal QSAR model were 0.849, 0.724 and 0.689, respectively. From the QSAR model, it could be seen that electrostatic, hydrophobic and steric interactions were closely related to inhibitory activity, which was consistent with the docking results. Based on the docking and QSAR results, five new compounds with high predicted activities were designed.     

3D-QSAR CoMFA studies on bis-coumarine analogues as urease inhibitors: a strategic design in anti-urease agents

A 3D-QSAR study has been performed on thirty (30) bis-coumarine derivatives to correlate their chemical structures with their observed urease inhibitory activity. Due to the absence of information on their active mechanism, comparative molecular field analysis (CoMFA) was used in the study. Two different properties: steric, electrostatic, assumed to cover the major contributions to ligand binding, were used to generate the 3D-QSAR model. Significant cross-validated correlation coefficients q(2) (0.558) and r(2) (0.992) for CoMFA were obtained, indicating the statistical significance of this class of compounds. The red electrostatic contour map highlighting those portion of compounds which may be interacting with nickel metal center in the active site of urease; while the blue contour map indicates positively charged groups in the ligands have improved biological activity and thus lower the IC(50)s. The steric contour map shows that bulkier substitutions at the 'R' position are detrimental to ligand receptor interaction. Actual urease inhibitory activities of this class and the predicted values were in good agreement with the experimental results. Moreover, from the contour maps, the key features vital to ligand binding have been identified, which are important for us to trace the important properties and gain insight into the potential mechanisms of intermolecular interactions between the ligand and receptor.     

CoMFA and CoMSIA 3D-QSAR analysis of diaryloxy-methano-phenanthrene derivatives as anti-tubercular agents

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) based on three-dimensional quantitative structure-activity relationship (3D-QSAR) studies were conducted on a series (44 compounds) of diaryloxy-methano-phenanthrene derivatives as potent antitubercular agents. The best predictions were obtained with a CoMFA standard model (q (2)=0.625, r (2)=0.994) and with CoMSIA combined steric, electrostatic, and hydrophobic fields (q (2)=0.486, r (2)=0.986). Both models were validated by a test set of seven compounds and gave satisfactory predictive r (2) values of 0.999 and 0.745, respectively. CoMFA and CoMSIA contour maps were used to analyze the structural features of the ligands to account for the activity in terms of positively contributing physicochemical properties: steric, electrostatic, and hydrophobic fields. The information obtained from CoMFA and CoMSIA 3-D contour maps can be used for further design of phenanthrene-based analogs as anti-TB agents. The resulting contour maps, produced by the best CoMFA and CoMSIA models, were used to identify the structural features relevant to the biological activity in this series of analogs. Further analysis of these interaction-field contour maps also showed a high level of internal consistency. This study suggests that introduction of bulky and highly electronegative groups on the basic amino side chain along with decreasing steric bulk and electronegativity on the phenanthrene ring might be suitable for designing better antitubercular agents.     

CoMFA and CoMSIA 3D-QSAR analysis of DMDP derivatives as anti-cancer agents

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) based on three dimensional quantitative structure-activity relationship (3D-QSAR) studies were conducted on a series (78 compounds) of 2, 4-diamino-5-methyl-5-deazapteridine (DMDP) derivatives as potent anticancer agents. The best prediction were obtained with a CoMFA standard model (q(2) = 0.530, r(2) = 0.903) and with CoMSIA combined steric, electrostatic, hydrophobic and hydrogen bond donor fields (q(2) = 0.548, r(2) = 0.909). Both models were validated by a test set of ten compounds producing very good predictive r(2) values of 0.935 and 0.842, respectively. CoMFA and CoMSIA contour maps were then used to analyze the structural features of ligands to account for the activity in terms of positively contributing physiochemical properties such as steric, electrostatic, hydrophobic and hydrogen bond donor fields. The resulting contour maps produced by the best CoMFA and CoMSIA models were used to identify the structural features relevant to the biological activity in this series of analogs. This study suggests that the highly electropositive substituents with low steric tolerance are required at 5 position of the pteridine ring and bulky electronegatve substituents are required at the meta-position of the phenyl ring. The information obtained from CoMFA and CoMSIA 3-D contour maps can be used for the design of deazapteridine-based analogs as anticancer agents.