Molecular modeling and 3D-QSAR studies in 2-aziridinyl-and 2,3-bis(aziridinyl)-1,4-naphthoquinonyl sulfonate and acylate derivatives as potential antimalarial agents

Malaria is still continuing to be one of the most dreadful diseases of the tropical countries particularly due to the development of resistance to the existing antimalarials. From observed, antimalarial activity of 2-aziridinyl- and 2,3-bis(aziridinyl)-1,4-naphthoquinonyl sulfonate and acylate derivatives acting through redox cycling mechanism, molecular modeling and three dimensional-quantitative structure activity relationship (3D-QSAR) studies have been carried out on a set of 63 compounds to identify important pharmacophors. Among several 3D-QSAR models generated, three models with correlation coefficient r > 0.82, match > 0.60 and chance = 0.00 have shown two common biophoric sites: one being the oxygen atom at position 1 of the naphthoquinone ring in terms of pi-population, charge and electron donating ability while the second being the center of the phenyl ring in terms of its 6pi-electrons. In addition to these sites, the models also share two common secondary sites: one positively contributing H-acceptor site while the second site contributing negatively in terms of steric refractivity. All these models showed good agreement between the experimental, calculated and predicted antimalarial activities.     

Comparison of CoMFA Models for Salmonella Typhimurium TA98, TA100+TA98+S9 and TA100+S9 Mutagenicity of Nitroaromatics

Abstract

Comparative Molecular Field Analysis (CoMFA) was applied to a comprehensive data set of heterogeneous nitroaromatics tested in Salmonella typhimurium TA98 and TA100 with and without S9 microsomal activation. The four CoMFA models developed agree with postulated mechanisms of mutagenicity, and explain over 70% of the corresponding mutagenic variance. The standard deviation coefficient contours common in the four models included high electronic density regions equivalent to C4-C5 in the pyrene ring, and an electron deficient site equivalent to C6. These areas are associated with high mutagenicity. Electron deficient areas may be related with the nitroreductive bioactivation of nitroaromatics. Electron rich sites may be involved with oxidative mechanisms analogous to the bioactivation pathway of polycyclic aromatic hydrocarbons. The contribution of steric factors to mutagenicity follows the order TA98 + S9 > TA98 > TA100 + S9 > TA100. The models indicated that increasing bulk perpendicular to the aromatic plane would decrease mutagenicity, but increasing the aromatic ring system along a region corresponding to C6-C7 in 1-nitropyrene would increase mutagenicity.     

Molecular modelling studies on adamantane-based Ebola virus GP-1 inhibitors using docking,pharmacophore and 3D-QSAR

The pathogenic Ebola virus (EBOV) causes a potential health risk and global spread. To date, few drugs are available for the treatment of Ebola virus disease (EVD) that allow researchers to use computational methods for designing potential drugs. The developed PHASE-based common six-point pharmacophore hypothesis (AADHPR_1) showed the necessity of two hydrogen bond acceptor features, one hydrogen bond donor feature, one hydrophobic group feature, one positively ionizable and one aromatic ring feature for further designing. We developed best 3D-QSAR models with high regression coefficients for the training (r²>0.82) and test (Q²>0.5) sets for both atoms-based and field-based 3D-QSAR models. The molecule 1A-4 (docking score = –4.711 kcal/mol) was obtained as best docked (SP mode) on Ebola virus envelope glycoprotein (PDB ID-3CSY) as compared with the standards oseltamivir (docking score = –4.39 kcal/mol) and zanamivir (docking score = –3.392 kcal/mol). The obtained ZINC hit ZINC58935541 showed a good docking score of –4.892 kcal/mol. The ZINC58935541 molecule also showed a strong binding affinity towards the receptor cavity of Ebola virus envelope glycoprotein when simulated for 1.2 ns. The good QikProp parameters reflect the fact that this molecule, upon optimization into a lead, might become a good candidate for the treatment of EVD.     

Additional synthesis on thiophene-containing trisubstituted methanes (TRSMs) as inhibitors of M. tuberculosis and 3D-QSAR studies

Abstract

We earlier reported thiophene-containing trisubstituted methanes (TRSMs) as novel cores carrying anti-tubercular activity, and identified S006-830 as the phenotypic lead with potent bactericidal activity against single- and multi-drug resistant clinical isolates of Mycobacterium tuberculosis (M. tb). In this work, we carried out additional synthesis of several TRSMs. The reaction scheme essentially followed the Grignard reaction and Friedel–Crafts alkylation, followed by insertion of a dialkylaminoethyl chain. We also performed microbiological evaluations including in vitro screening against the virulent strain M. tb H37Rv, cytotoxicity assessment in the Vero C-1008 cell line, and 3D-QSAR studies with comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA). CoMFA and CoMSIA models yielded good statistical results in terms of q² and r² values, suggesting the validity of the models. It was concluded that a para-substituted benzene ring with bulkier electron-donating groups and aminoalkyl chains are required for higher inhibitory capacity against M. tuberculosis. We believe that these insights will rationally guide the design of newer, optimal, TRSMs.     

4D-QSAR analysis of a series of antifungal p450 inhibitors and 3D-pharmacophore comparisons as a function of alignment

A training set of 55 antifungal p450 analogue inhibitors was used to construct receptor-independent four-dimensional quantitative structure-activity relationship (RI 4D-QSAR) models. Ten different alignments were used to build the models, and one alignment yields a significantly better model than the other alignments. Two different methodologies were used to measure the similarity of the best 4D-QSAR models of each alignment. One method compares the residual of fit between pairs of models using the cross-correlation coefficient of their residuals of fit as a similarity measure. The other method compares the spatial distributions of the IPE types (3D-pharmacophores) of pairs of 4D-QSAR models from different alignments. Optimum models from several different alignments have nearly the same correlation coefficients, r(2), and cross-validation correlation coefficients, xv-r(2), yet the 3D-pharmacophores of these models are very different from one another. The highest 3D-pharmacophore similarity correlation coefficient between any pair of 4D-QSAR models from the 10 alignments considered is only 0.216. However, the best 4D-QSAR models of each alignment do contain some proximate common pharmacorphore sites. A test set of 10 compounds was used to validate the predictivity of the best 4D-QSAR models of each alignment. The "best" model from the 10 alignments has the highest predictivity. The inferred active sites mapped out by the 4D-QSAR models suggest that hydrogen bond interactions are not prevalent when this class of P450 analogue inhibitors binds to the receptor active site. This feature of the 4D-QSAR models is in agreement with the crystal structure results that indicate no ligand-receptor hydrogen bonds are formed.     

Drug discovery studies on quinoline-based derivatives as potential antimalarial agents

Molecular modelling studies were performed to identify the essential structural requirements of quinoline-based derivatives for improving their antimalarial activity. The developed CoMFA, CoMSIA and HQSAR models for a training set comprising 37 derivatives showed good statistical significance in terms of internal cross validation (q²) 0.70, 0.69 and 0.80 and non-cross validation (r²) 0.80, 0.79 and 0.80. Also, the predicted r² values (r²_pred) of 0.63, 0.61 and 0.72 for a test set consisting of 12 compounds suggested significant predicting ability of the models. Structural features were correlated in terms of steric, electrostatic, hydrophobic, hydrogen bond donor and hydrogen bond acceptor interactions. Furthermore, the bioactive conformation was explored and explained by docking compounds #28, 32 and 40 into the active binding site of lactate dehydrogenase of Plasmodium falciparum. The QSAR models, contour map and docking binding affinity obtained could be successfully utilized as a guiding tool for the design and discovery of novel quinoline-based derivatives with potent antimalarial activity.     

Comparison of Estrogen Receptor α and β Subtypes Based on Comparative Molecular Field Analysis (CoMFA)

Abstract

A substantial body of evidence indicates that both humans and wildlife suffer adverse health effects from exposure to environmental chemicals that are capable of interacting with the endocrine system. The recent cloning of the estrogen receptor β subtype (ER-β) suggests that the selective effects of estrogenic compounds may arise in part by the control of different subsets of estrogen-responsive promoters by the two ER subtypes, ER-α and ER-β. In order to identify the structural prerequisites for ligand-ER binding and to discriminate ER-α and ER-3 in terms of their ligand-binding specificities, Comparative Molecular Field Analysis (CoMFA) was employed to construct a three-dimensional Quantitative Structure-Activity Relationship (3D-QSAR) model on a data set of 31 structurally-diverse compounds for which competitive binding affinities have been measured against both ER-α and ER-β. Structural alignment of the molecules in CoMFA was achieved by maximizing overlap of their steric and electrostatic fields using the Steric and Electrostatic ALignment (SEAL) algorithm. The final CoMFA models, generated by correlating the calculated 3D steric and electrostatic fields with the experimentally observed binding affinities using partial least-squares (PLS) regression, exhibited excellent self-consistency (r ² > 0.99) as well as high internal predictive ability (q ² > 0.65) based on cross-validation. CoMFA-predicted values of RBA for a test set of compounds outside of the training set were consistent with experimental observations. These CoMFA models can serve as guides for the rational design of ER ligands that possess preferential binding affinities for either ER-α or ER-β. These models can also prove useful in risk assessment programs to identify real or suspected EDCs.     

Imidazole-containing farnesyltransferase inhibitors: 3D quantitative structure–activity relationships and molecular docking

One of the most promising anticancer and recent antimalarial targets is the heterodimeric zinc-containing protein farnesyltransferase (FT). In this work, we studied a highly diverse series of 192 Abbott-initiated imidazole-containing compounds and their FT inhibitory activities using 3D-QSAR and docking, in order to gain understanding of the interaction of these inhibitors with FT to aid development of a rational strategy for further lead optimization. We report several highly significant and predictive CoMFA and CoMSIA models. The best model, composed of CoMFA steric and electrostatic fields combined with CoMSIA hydrophobic and H-bond acceptor fields, had r ² = 0.878, q ² = 0.630, and r _pred² = 0.614. Docking studies on the statistical outliers revealed that some of them had a different binding mode in the FT active site based on steric bulk and available active site space, explaining why the predicted activities differed from the experimental activities. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

3D-Pharmacophore mapping of thymidine-based inhibitors of TMPK as potential antituberculosis agents

Tuberculosis (TB) is the primary cause of mortality among infectious diseases. Mycobacterium tuberculosis monophosphate kinase (TMPKmt) is essential to DNA replication. Thus, this enzyme represents a promising target for developing new drugs against TB. In the present study, the receptor-independent, RI, 4D-QSAR method has been used to develop QSAR models and corresponding 3D-pharmacophores for a set of 81 thymidine analogues, and two corresponding subsets, reported as inhibitors of TMPKmt. The resulting optimized models are not only statistically significant with r ² ranging from 0.83 to 0.92 and q ² from 0.78 to 0.88, but also are robustly predictive based on test set predictions. The most and the least potent inhibitors in their respective postulated active conformations, derived from each of the models, were docked in the active site of the TMPKmt crystal structure. There is a solid consistency between the 3D-pharmacophore sites defined by the QSAR models and interactions with binding site residues. Moreover, the QSAR models provide insights regarding a probable mechanism of action of the analogues.     

Spectrophotofluorimetric Study of the Interaction between Trazodone Hydrochloride and Bovine Serum Albumin

Abstract

The binding of trazodone hydrochloride (TZH), an antidepressant drug, to bovine serum albumin (BSA) has been investigated by fluorescence spectroscopic analysis. The fluorescence emission of BSA (λ_em=350 nm) was quenched by TZH while that of this ligand was enhanced (λ_em=443 nm). The spectral behavior was consistent with the static quenching mechanism, and the apparent binding constant, K _a (1.05×10⁴ l mol^?1) as well as binding site number, n (～1), were estimated. Thermodynamic parameters obtained from the measured data at different temperatures showed that the binding of TZH to BSA involved predominantly hydrophobic interactions as well as smaller contributions from electrostatic forces. Phenylbutazone and ibuprofen were utilized as competitive markers for sites I and II, respectively, in the interaction of TZH with BSA. This competitive displacement procedure indicated that the likely binding was site I, i.e., subdomain IIA, and this was supported by the observation that up to 50% of this site marker, phenylbutazone, could be exchanged with TZH whilst only a few percent of ibuprofen were so affected.     

In silico study directed towards identification of the key structural features of GyrB inhibitors targeting MTB DNA gyrase: HQSAR,CoMSIA and molecular dynamics simulations

ABSTRACT

Mycobacterium tuberculosis DNA gyrase subunit B (GyrB) has been identified as a promising target for rational drug design against fluoroquinolone drug-resistant tuberculosis. In this study, we attempted to identify the key structural feature for highly potent GyrB inhibitors through 2D-QSAR using HQSAR, 3D-QSAR using CoMSIA and molecular dynamics (MD) simulations approaches on a series of thiazole urea core derivatives. The best HQSAR and CoMSIA models based on IC₅₀ and MIC displayed the structural basis required for good activity against both GyrB enzyme and mycobacterial cell. MD simulations and binding free energy analysis using MM-GBSA and waterswap calculations revealed that the urea core of inhibitors has the strongest interaction with Asp79 via hydrogen bond interactions. In addition, cation-pi interaction and hydrophobic interactions of the R₂ substituent with Arg82 and Arg141 help to enhance the binding affinity in the GyrB ATPase binding site. Thus, the present study provides crucial structural features and a structural concept for rational design of novel DNA gyrase inhibitors with improved biological activities against both enzyme and mycobacterial cell, and with good pharmacokinetic properties and drug safety profiles.     

Amidocrownophanes as Anion Receptor

Crownophanes composed of 28-membered ring atoms having two hydroxy groups, two amide groups, and aromatic parts such as naphthalene rings and either pyridine or benzene ring, can bind anions with high affinity and selectivity. The anion-coordination ability of these species has been observed by ¹H NMR techniques. As anion guest molecules, we selected some halides, dihydrogenphosphate and acetate ions. It has been found that amidocrownophanes, 3 and 4, can recognize anions in the order;H₂PO ₄ ⁻ >F⁻>CH₃COO⁻>Cl⁻>>Br⁻ and I⁻, while not only 1, 2, and 5 having no hydroxy group but also 6 having 27-membered ring have no ability for anion recognition under the same conditions. In order to exhibit the recognition ability for anion receptor, plural amide groups, hydroxy groups, and m-phenylene or 1,6-pyridyl rigid part play an important role in this macrocyclic system.     

Superimposition evaluation of ecdysteroid agonist chemotypes through multidimensional QSAR

The EC₅₀ values for a training set of 66 ecdysteroids and 97 diacylhydrazines were measured in the ecdysteroid-responsive Drosophila melanogaster B_II cell line, a prototypical homologous inducible gene expression system. Each of eight superimposition hypotheses for the folded diacylhydrazine conformation was evaluated and ranked on the basis of CoMFA and 4D-QSAR Q² values for the training set and R² values for a 52-member test set comprising randomly-chosen diacylhydrazines and chronologically-chosen ecdysteroids for which data became available after model construction. Both 4D-QSAR and CoMFA rate a common superimposition as the preferred one. Two additional superimpositions, with somewhat weaker 4D-QSAR and CoMFA consensus, nonetheless share several important topological features. The resultant QSAR models address the question of relative binding orientation of the two ligand families and can be useful as a virtual screen for new chemotypes.     

Ca7.96Cu0.04Ge5O18: a new calcium germanate with GeO4 and Ge3O10 units

The title compound, octa­calcium copper penta­germanium octa­deca­oxide, represents a new inter­mediate phase between CaO and GeO₂, and has not previously been reported in the literature. The structure consists of three different Ge sites, two of them on general 8d positions, site symmetry 1, one on special position 4d, site symmetry 2. Three of the five Ca sites occur on 8d positions, site symmtery 1, one Ca is on 4b with site symmetry and one Ca is on 4c with site symmetry 2. All nine O atoms have symmetry 1 (8d position). By sharing common edges, the Ca sites form infinite bands parallel to the c axis, and these bands are inter­connected by isolated GeO₄ and Ge₃O₁₀ units. These (100) layers are stacked along a in an ABAB… sequence, with the B layer being inverted and displaced along b/2.     

The Loss of a hydroxyl group from the molecular ions of alkylnitrobenzenes

It is confirmed that the loss of HO˙ from the molecular ion of o-nitrotoluene involves exclusively a hydrogen from the methyl group. However, in higher homologues hydrogen atoms from non-benzylic sites are also implicated. With such compounds this fragmentation mode is shown not only by the ortho but, to a lesser extent, by the meta and para isomers as well. The proportion of the total ion current borne by the [M – 17]⁺ ion follows the order ortho > meta > para, which is attributed to substituent migration around the ring with a hydroxyl radical only being lost when the groups are on adjacent ring atoms. Other ions present in the spectra point to interaction between substituents to form a new heterocyclic ring.     

Identification of potential CRAC channel inhibitors: Pharmacophore mapping, 3D-QSAR modelling,and molecular docking approach

Upregulation of store-operated Ca²⁺ influx via ORAI1, an integral component of the CRAC channel, is responsible for abnormal cytokine release in active rheumatoid arthritis, and therefore ORAI1 has been proposed as an attractive molecular target. In this study, we attempted to predict the mechanical insights of ORAI1 inhibitors through pharmacophore modelling, 3D-QSAR, molecular docking and free energy analysis. Various hypotheses of pharmacophores were generated and from that, a pharmacophore hypothesis with two hydrogen bond acceptors, one hydrogen bond donor and two aromatic rings (AADRR) resulted in a statistically significant 3D-QSAR model (r² = 0.84 and q² = 0.74). We believe that the obtained statistical model is a reliable QSAR model for the diverse dataset of inhibitors against the IL-2 production assay. The visualization of contours in active and inactive compounds generated from the 3D-QSAR models and molecular docking studies revealed major interaction with GLN108, HIS113 and ASP114, and interestingly, these residues are located near the Ca²⁺ selectivity filter region. Free energy binding analysis revealed that Coulomb energy, van der Waals energy and non-polar solvation terms are more favourable for ligand binding. Thus, the present study provides the physical and chemical requirements for the development of novel ORAI1 inhibitors with improved biological activity.     

Synthetic modifications in ethyl 2-amino-4-methylthiazole-5-carboxylate: 3D QSAR analysis and antimicrobial study

Ethyl 2-amino-4-methylthiazole-5-carboxylate (1) is a very substantial derivative of the thiazole group. This derivative has been modified and has been synthesized using readily available materials. The structures of synthesized derivatives were confirmed by IR, ¹H NMR, ¹³C NMR, and mass spectral techniques. The newly prepared compounds (5a–k) were studied for their antimicrobial activities against strains of bacteria (Escherichia coli, Pseudomonas aeruginosa, Staphylococcus aureus, and Streptococcus pyogenes) and fungi (Candida albicans, Aspergillus niger, and Aspergillus clavatus) using serial dilution method. Structure-activity relationship was carried out by CoMFA and CoMSIA with the help of 3D-QSAR analysis. As per the database alignment, the CoMFA and CoMSIA models were developed. Information based on these models is useful for structure-activity relationships of the reported molecules.     

Characterization of β<Subscript>3</Subscript>-adrenergic receptor: determination of pharmacophore and 3D QSAR model for β<Subscript>3</Subscript> adrenergic receptor agonism

Summary The β₃-adrenoreceptor (β₃-AR) has been shown to mediate various pharmacological and physiological effects such as lipolysis, thermogenesis, and intestinal smooth muscle relaxation. It also plays an important role in glucose homeostasis and energy balance. Molecular modeling studies were undertaken to develop predictive pharmacophoric hypothesis and 3D-QSAR model, which may explain variations in β₃-AR agonistic activity in terms of chemical features and physicochemical properties. The two softwares, CATALYST for pharmacophoric alignment and APEX-3D for 3D-QSAR modeling were used to establish the structure activity relationships for β₃-AR agonistic activity. Among the several statistically significant models, the selection of the best pharmacophore and 3D-QSAR model was based on its ability to estimate the activity of external test sets of similar and different structural types along with the reasonable consistency of the model with the limited information of the active site of β₃-AR. The final 3D-QSAR model was derived using the pharmacophoric alignments from the hypothesis which consisted of four chemical features: basic or positive ionizable feature on the nitrogen of the aryloxypropylamino group, two ring aromatic features corresponding to the phenyl ring of the phenoxide and the benzenesulphonamido groups and a hydrogen-bond donor (HBD) in the vicinity of the nitrogen atom of the benzenesulphonamido group with the most active molecule mapping in an energetically favorable extended conformation. This hypothesis was in agreement with the site directed mutagenesis studies on human β₃-AR and correlated well the observed and estimated activity both in, training and both the external test sets. It also mapped reasonably well to six β₃-AR agonists of different structural classes under clinical development and thus this hypothesis may have a universal applicability in providing a powerful template for virtual screening and also for designing new chemical entities (NCEs) as β₃-AR agonists.^★CDRI communication number 6202. *To whom correspondence should be addressed. Fax: +91-0522-223405; E-mail: anilsak@hotmail.com