3D-QSAR Analysis Research Articles

Combining ligand-based and structure-based drug design approaches to study the structure-activity relationships of a β-carboline derivative series

In this study, we investigated the structure-activity relationships of a series of β-carboline alkaloid derivatives using the 2D-QSAR and molecular docking, in order to identify the mode of interaction between β-carboline derivatives and the PLK1 kinase, and determine their key substituents responsible for the cytotoxic activity. The obtained QSAR models using multiple linear regression (MLR) and partial least squares (PLS) methods showed a high correlation between the experimental activity and the predicted one by PLS (R2PLS = 0.82, q2 = 0.72) and MLR (R2MLR = 0.82, q2 = 0.72). An external dataset was used to test the extrapolation power of the models which resulted in an R2PLS (EV) = 0.76; RMSE = 0.39. The 2D-QSAR analysis reveals that lipophilicity plays an important role in the cytotoxic activity of this group of β-carboline derivatives. Indeed, the molecular docking study into the active site of the polo-like kinase (PLK1) revealed that the most active ligand 57 shows higher binding energy and interacts, especially by H-bonds and hydrophobic interactions, with the active site of the PLK1 kinase. Consequently, the results obtained from the 2D-QSAR and docking studies provided a useful tool to design new and potent β-carboline derivatives as cytotoxic agents.

Synthesis, bioactivity, action mode and 3D-QSAR of novel anthranilic diamide derivatives

To study the pesticide effect, action mode, structure-activity relationships (SARs) of anthranilic diamide insecticide and screen highly active pesticides, novel anthranilic diamide derivatives were synthesized. Bioassays indicated that all of the title compounds displayed 100% mortality against diamondback moth and oriental armyworm at 100mg/L, among which 12v and 12w showed 100% insecticidal acitvity at 5mg/L. Surprisingly compound 12w exhibited better insecticidal acitvity than commercialized chlorantraniliprole against Pyrausta nubilalis (0.1mg/L) and Cnaphalocrocis Medinalis (2mg/L). 3D-QSAR and SARs statistical analysis revealed that title compounds with R2 fixed as methoxy had the highest probability possessing high activity. The calcium fluorescence measurements on neurons revealed that E series compounds containing pyrazinyl may have a molecular target different from caffeine on ryanodine receptors rather than the voltage-gated calcium channel present on cytomembran.

Eco-friendly synthesis, in vitro anti-proliferative evaluation, and 3D-QSAR analysis of a novel series of monocationic 2-aryl/heteroaryl-substituted 6-(2-imidazolinyl)benzothiazole mesylates

Herein, we describe the synthesis of twenty-one novel water-soluble monocationic 2-aryl/heteroaryl-substituted 6-(2-imidazolinyl)benzothiazole mesylates 3a-3u and present the results of their anti-proliferative assays. Efficient syntheses were achieved by three complementary simple two-step synthetic protocols based on the condensation reaction of aryl/heteroaryl carbaldehydes or carboxylic acid. We developed an eco-friendly synthetic protocol using glycerol as green solvent, particularly appropriate for the condensation of thermally and acid-sensitive heterocycles such as furan, benzofuran, pyrrole, and indole. Screening of anti-proliferative activity was performed on four human tumour cell lines in vitro including pancreatic cancer (CFPAC-1), metastatic colon cancer (SW620), hepatocellular carcinoma (HepG2), and cervical cancer (HeLa), as well as in normal human fibroblast cell lines. All tested compounds showed strong to moderate anti-proliferative activity on tested cell lines depending on the structure containing aryl/heteroaryl moiety coupled to 6-(2-imidazolinyl)benzothiazole moiety. The most potent cytostatic effects on all tested cell lines with [Formula: see text] values ranging from 0.1 to 3.70[Formula: see text] were observed for benzothiazoles substituted with naphthalene-2-yl 3c, benzofuran-2-yl 3e, indole-3-yl 3j, indole-2-yl 3k, quinoline-2-yl 3s, and quinoline-3-yl 3t and derivatives substituted with phenyl 3a, naphthalene-1-yl 3b, benzothiazole-2-yl 3g, benzothiazole-6-yl 3h, N-methylindole-3-yl 3l, benzimidazole-2-yl 3n, benzimidazole-5(6)-yl 3o, and quinolone-4-yl 3u with [Formula: see text] values ranging from 1.1 to 29.1[Formula: see text]. Based on obtained anti-proliferative activities, 3D-QSAR models for five cell lines were derived. Molecular volume, molecular surface, the sum of hydrophobic surface areas, molecular mass, and possibility of making dispersion forces were identified by QSAR analyses as molecular properties that are positively correlated with anti-proliferative activity, while compound's capability to accept H-bond was identified as a negatively correlated property. Comparison of molecular properties identified for different cell lines enabled assumptions about similarity of mode of action through which anti-proliferative activities against different cell lines are accomplished. Novel compounds that are predicted to have enhanced activities in comparison with herein presented ones were designed using 3D-QSAR analysis as guideline.

Insight into the interaction mechanism of human SGLT2 with its inhibitors: 3D-QSAR studies, homology modeling, and molecular docking and molecular dynamics simulations.

Human sodium-dependent glucose co-transporter 2 (hSGLT2) is a crucial therapeutic target in the treatment of type 2 diabetes. In this study, both comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were applied to generate three-dimensional quantitative structure-activity relationship (3D-QSAR) models. In the most accurate CoMFA-based and CoMSIA-based QSAR models, the cross-validated coefficients (r2cv) were 0.646 and 0.577, respectively, while the non-cross-validated coefficients (r2) were 0.997 and 0.991, respectively, indicating that both models were reliable. In addition, we constructed a homology model of hSGLT2 in the absence of a crystal structure. Molecular docking was performed to explore the bonding mode of inhibitors to the active site of hSGLT2. Molecular dynamics (MD) simulations and binding free energy calculations using MM-PBSA and MM-GBSA were carried out to further elucidate the interaction mechanism. With regards to binding affinity, we found that hydrogen-bond interactions of Asn51 and Glu75, located in the active site of hSGLT2, with compound 40 were critical. Hydrophobic and electrostatic interactions were shown to enhance activity, in agreement with the results obtained from docking and 3D-QSAR analysis. Our study results shed light on the interaction mode between inhibitors and hSGLT2 and may aid in the development of C-aryl glucoside SGLT2 inhibitors.

3D-QSAR study of steroidal and azaheterocyclic human aromatase inhibitors using quantitative profile of protein\u2013ligand interactions

Aromatase is a member of the cytochrome P450 superfamily responsible for a key step in the biosynthesis of estrogens. As estrogens are involved in the control of important reproduction-related processes, including sexual differentiation and maturation, aromatase is a potential target for endocrine disrupting chemicals as well as breast cancer therapy. In this work, 3D-QSAR combined with quantitative profile of protein–ligand interactions was employed in the identification and characterization of critical steric and electronic features of aromatase-inhibitor complexes and the estimation of their quantitative contribution to inhibition potency. Bioactivity data on pIC50 values of 175 steroidal and 124 azaheterocyclic human aromatase inhibitors (AIs) were used for the 3D-QSAR analysis. For the quantitative description of the effects of the hydrophobic contact and nitrogen–heme–iron coordination on aromatase inhibition, the hydrophobicity density field model and the smallest dual descriptor Δf(r)S were introduced, respectively. The model revealed that hydrophobic contact and nitrogen–heme–iron coordination primarily determines inhibition potency of steroidal and azaheterocyclic AIs, respectively. Moreover, hydrogen bonds with key amino acid residues, in particular Asp309 and Met375, and interaction with the heme–iron are required for potent inhibition. Phe221 and Thr310 appear to be quite flexible and adopt different conformations according to a substituent at 4- or 6-position of steroids. Flexible docking results indicate that proper representation of the residues’ flexibility is critical for reasonable description of binding of the structurally diverse inhibitors. Our results provide a quantitative and mechanistic understanding of inhibitory activity of steroidal and azaheterocyclic AIs of relevance to adverse outcome pathway development and rational drug design.

Study of Radical Scavenging Activities of a Series of Flavonoids through 3D-QSAR Analysis

Study of Radical Scavenging Activities of a Series of Flavonoids through 3D-QSAR Analysis

Design of Novel N-Myristoyltransferase Inhibitors of Leishmania donovani Using Four-Dimensional Quantitative Structure-Activity Relationship Analysis

N-Myristoylation protein is catalyzed by N-myristoyltransferase (NMT), an essential target in Leishmania donovani, the causative agent of kala-azar. Four-dimensional quantitative structure-activity relationship (4D-QSAR) analysis was applied to a series of 77 Leishmania donovani NMT inhibitors. Then, three new compounds were proposed using QSAR models. In addition, molecular docking was performed to predict the binding affinities and interaction modes among the proposed compounds and the NMT active site. In silico absorption, distribution, metabolism and excretion (ADME) evaluation was performed and potential inhibitors demonstrated satisfactory pharmacokinetic properties.

Multiple molecular modelling studies on some derivatives and analogues of glutamic acid as matrix metalloproteinase-2 inhibitors

Matrix metalloproteinase-2 (MMP-2) is a potential target in anticancer drug discovery due to its association with angiogenesis, metastasis and tumour progression. In this study, 67 glutamic acid derivatives, synthesized and evaluated as MMP-2 inhibitors, were taken into account for multi-QSAR modelling study (regression-based 2D-QSAR, classification-based LDA-QSAR, Bayesian classification QSAR, HQSAR, 3D-QSAR CoMFA and CoMSIA as well as Open3DQSAR). All these QSAR studies were statistically validated individually. Regarding the 3D-QSAR analysis, the Open3DQSAR results were better than CoMFA and CoMSIA, although all these 3D-QSAR models supported each other. The importance of biphenylsulphonyl moiety over phenylacetyl/naphthylacetyl moieties was established due to its association with favourable steric and hydrophobic characters. HQSAR, LDA-QSAR and Bayesian classification QSAR studies also suggested that the biphenylsulphonamido group was better than the phenylacetylcarboxamido function. Additionally, glutamines were proven to be far better inhibitors than isoglutamines. Observations obtained from the current study were revalidated and supported by the earlier reported molecular modelling studies. Depending on these observations, newer glutamic acid-based compounds may be designed further in future for potent MMP-2 inhibitory activity.

Computational Analysis of Artimisinin Derivatives on the Antitumor Activities

The study on antitumor activities of artemisinin and its derivatives has been closely focused on in recent years. Herein, 2D and 3D QSAR analysis was performed on the basis of a series of artemisinin derivatives with known bioactivities against the non-small-cell lung adenocarcinoma A549 cells. Four QSAR models were successfully established by CoMSIA, CoMFA, topomer CoMFA and HQSAR approaches with respective characteristic values q2 = 0.567, R2 = 0.968, ONC = 5; q2 = 0.547, R2 = 0.980, ONC = 7; q2 = 0.559, R2 = 0.921, ONC = 7 and q2 = 0.527, R2 = 0.921, ONC = 6. The predictive ability of CoMSIA with r2 = 0.991 is the best one compared with the other three approaches, such as CoMFA (r2 = 0.787), topomer CoMFA (r2 = 0.819) and HQSAR (r2 = 0.743). The final QSAR models can provide guidance in structural modification of artemisinin derivatives to improve their anticancer activities.

Combined molecular modelling and 3D-QSAR study for understanding the inhibition of NQO1 by heterocyclic quinone derivatives.

A combination of three-dimensional quantitative structure-activity relationship (3D-QSAR), and molecular modelling methods were used to understand the potent inhibitory NAD(P)H:quinone oxidoreductase 1 (NQO1) activity of a set of 52 heterocyclic quinones. Molecular docking results indicated that some favourable interactions of key amino acid residues at the binding site of NQO1 with these quinones would be responsible for an improvement of the NQO1 activity of these compounds. The main interactions involved are hydrogen bond of the amino group of residue Tyr128, π-stacking interactions with Phe106 and Phe178, and electrostatic interactions with flavin adenine dinucleotide (FADH) cofactor. Three models were prepared by 3D-QSAR analysis. The models derived from Model I and Model III, shown leave-one-out cross-validation correlation coefficients (q2LOO ) of .75 and .73 as well as conventional correlation coefficients (R2 ) of .93 and .95, respectively. In addition, the external predictive abilities of these models were evaluated using a test set, producing the predicted correlation coefficients (r2pred ) of .76 and .74, respectively. The good concordance between the docking results and 3D-QSAR contour maps provides helpful information about a rational modification of new molecules based in quinone scaffold, in order to design more potent NQO1 inhibitors, which would exhibit highly potent antitumor activity.

Design and synthesis of 2,6-di(substituted phenyl)thiazolo[3,2-b]-1,2,4-triazoles as α-glucosidase and α-amylase inhibitors, co-relative Pharmacokinetics and 3D QSAR and risk analysis

Ten fused heterocyclic derivatives bearing the 2,6-di(subsituted phenyl)thiazolo[3,2-b]-1,2,4-triazoles as central rings were synthesized and structures of the compounds were established by analytical and spectral data using FTIR, EI-MS, 1H NMR and 13C NMR techniques. In vitro inhibitory activities of synthesized compounds on α-amylase, α-glucosidase and α-burylcholinesterase (α-BuChE) were evaluated using a purified enzyme assays. Compound 5c demonstrated strong and selective α-amylase inhibitory activity (IC50=1.1μmol/g). 5g exhibited excellent inhibition against α-glucosidase (IC50=1.2μmol/g) when compared with acarbose (IC50=4.7μmol/g) as a positive reference. Compound 5i was found to be most potent derivative against α-BuChE with the IC50 of 1.5μmol/g which was comparable to the value obtained for (4.7μmol/g) positive control (i.e. galantamine hydrobromide). Molecular dockings of synthesized compounds into the binding sites of human pancreatic α-amylase, intestinal maltase-glucoamylase and neuronal α-butrylcholinesterase allowed to shed light on the affinity and binding mode of these novel inhibitors. Preliminary structure–activity relationship (SAR) studies were carried out to understand the relationship between molecular structural features and inhibition activities of synthesized derivatives. These data suggested that compounds 5c, 5g and 5i are promising candidates for hitto- lead follow-up in the drug-discovery process for the treatment of Alzheimer’s disease and hyperinsulinamia.

2D-QSAR analysis on some 8-methoxy quinoline derivatives as H37RV (MTB) inhibitors with comparison of different model

2D-QSAR analysis on some 8-methoxy quinoline derivatives as H37RV (MTB) inhibitors with comparison of different model

Synthesis and antiproliferative activity of amino-substituted benzimidazo[1,2- $${\varvec{a}}$$ a ]quinolines as mesylate salts designed by 3D-QSAR analysis

An experimental search for new benzimidazole derivatives with enhanced antiproliferative activity was successfully guided by QSAR modelling. Robust 3D-QSAR models were derived on an available database of compounds with previously measured activities. Our QSAR analysis revealed that an increase of the antiproliferative activities towards H460, HCT 116, MCF-7 and SW 620 cells will be obtained if new compounds are charged at a pH range from 5 to 7 and if their hydrophobicity is increased compared to the dataset compounds. Novel benzimidazo[1,2-a]quinolines bearing quarternary amino groups with corresponding aliphatic chains were designed, and their antiproliferative activities were computationally predicted. Using uncatalysed microwave-assisted amination reactions, 14 novel compounds were obtained to assess their antiproliferative activities towards H460, HCT 116, MCF-7, and SW 620 tumour cell lines in vitro. Novel compounds showed antiproliferative activities at micromolar and submicromolar inhibition concentrations. Experimental measurements of antiproliferative activities validation the QSAR models showing very good agreement between experimentally measured activities and computational predictions. In an attempt to elucidate the mode of action through which benzimidazole derivatives accomplish their antiproliferative activities, thermal denaturation experiments were performed to test their DNA-binding properties.

Structure-Activity Relationships of Pentacyclic Triterpenoids as Potent and Selective Inhibitors against Human Carboxylesterase 1.

Human carboxylesterase 1 (hCE1), one of the most important serine hydrolases distributed in liver and adipocytes, plays key roles in endobiotic homeostasis and xenobiotic metabolism. This study aimed to find potent and selective inhibitors against hCE1 from phytochemicals and their derivatives. To this end, a series of natural triterpenoids were collected and their inhibitory effects against human carboxylesterases (hCEs) were assayed using D-Luciferin methyl ester (DME) and 6,8-dichloro-9,9-dimethyl-7-oxo-7,9-dihydroacridin-2-yl benzoate (DDAB) as specific optical substrate for hCE1, and hCE2, respectively. Following screening of a series of natural triterpenoids, oleanolic acid (OA), and ursolic acid (UA) were found with strong inhibitory effects on hCE1 and relative high selectivity over hCE2. In order to get the highly selective and potent inhibitors of hCE1, a series of OA and UA derivatives were synthesized from OA and UA by chemical modifications including oxidation, reduction, esterification, and amidation. The inhibitory effects of these derivatives on hCEs were assayed and the structure-activity relationships of tested triterpenoids as hCE1 inhibitors were carefully investigated. The results demonstrated that the carbonyl group at the C-28 site is essential for hCE1 inhibition, the modifications of OA or UA at this site including esters, amides and alcohols are unbeneficial for hCE1 inhibition. In contrast, the structural modifications on OA and UA at other sites, such as converting the C-3 hydroxy group to 3-O-β-carboxypropionyl (compounds 20 and 22), led to a dramatically increase of the inhibitory effects against hCE1 and very high selectivity over hCE2. 3D-QSAR analysis of all tested triterpenoids including OA and UA derivatives provide new insights into the fine relationships linking between the inhibitory effects on hCE1 and the steric-electrostatic properties of triterpenoids. Furthermore, both inhibition kinetic analyses and docking simulations demonstrated that compound 22 was a potent competitive inhibitor against hCE1-mediated DME hydrolysis. All these findings are very helpful for medicinal chemists to design and develop highly selective and more potent hCE1 inhibitors for biomedical applications.

Synthesis, biological screening, POM, and 3D-QSAR analyses of some novel pyrazolic compounds

A series of new pyrazolic heterocyclic compounds were prepared in good and excellent yields and characterized by proton and carbon nuclear magnetic resonance, infrared, and mass spectroscopy studies. These products were screened in vitro against three bacterial pathogens, namely Bacillus subtilis, Micrococcus luteus, and Escherichia coli and antifungal potential, against Fusarium oxysporum f.sp.albedinis. A considerable and excellent activity was recorded with respect to the two studied microorganisms. A good correlation was obtained between the experimental results and the theoretical predictions of bioavailability using Petra/Osiris/Molinspiration suite (Petra/Osiris/Molinspiration containing Lipinski’s rule-of-five). The quantitative structure activity relationship approach has been analyzed to support the Petra/Osiris/Molinspiration results and composite indexes of some quantum chemical parameters were constructed in order to characterize the inhibition performance of the tested molecules.

Design, synthesis, cytotoxicity and 3D-QSAR analysis of new 3,6-disubstituted-1,2,4,5-tetrazine derivatives as potential antitumor agents

We synthesized two new series of 3-substituted-6-(2,5-dimethylpyrazol-1-yl)-1,2,4,5-tetrazines and analysed them for a potential role as antitumor agents. Twenty-two compounds were obtained, and four molecular structures were determined by X-ray diffraction analysis. Using flow cytometry and MTT assay, potential action on cell toxicity was determined for each of the compounds for four cancer cell lines. The potency and selectivity demonstrated by these compounds are dependent on the cancer cell line, where the following compounds were found the most promising agents against certain cell lines: compounds 1i and 1j for HL-60 cells, 1a and 1b on HCT116 cells, 1f on Hela cells and 2h on H1975 cells. The action exerted by these compounds is comparable to the well-known cancer treatment drug etoposide and higher than vatalanib. To arrive at the structural requirements for activity on each cell line, a SAR and 3D-QSAR analysis was carried out. From the 3D-QSAR models, steric and electronic features were identified in the aromatic centres, and were key components for cytotoxic activity on HL-60 cell lines. The cytometry results suggest that some tetrazine derivatives induce apoptosis on HCT116 cells.

Pharmacophore generation, atom-based 3D-QSAR and molecular dynamics simulation analyses of pyridine-3-carboxamide-6-yl-urea analogues as potential gyrase B inhibitors

DNA gyrase subunit B (GyrB) is an attractive drug target for the development of antibacterial agents with therapeutic potential. In the present study, computational studies based on pharmacophore modelling, atom-based QSAR, molecular docking, free binding energy calculation and dynamics simulation were performed on a series of pyridine-3-carboxamide-6-yl-urea derivatives. A pharmacophore model using 49 molecules revealed structural and chemical features necessary for these molecules to inhibit GyrB. The best fitted model AADDR.13 was generated with a coefficient of determination (r²) of 0.918. This model was validated using test set molecules and had a good r² of 0.78. 3D contour maps generated by the 3D atom-based QSAR revealed the key structural features responsible for the GyrB inhibitory activity. Extra precision molecular docking showed hydrogen bond interactions with key amino acid residues of ATP-binding pocket, important for inhibitor binding. Further, binding free energy was calculated by the MM-GBSA rescoring approach to validate the binding affinity. A 10 ns MD simulation of inhibitor #47 showed the stability of the predicted binding conformations. We identified 10 virtual hits by in silico high-throughput screening. A few new molecules were also designed as potent GyrB inhibitors. The information obtained from these methodologies may be helpful to design novel inhibitors of GyrB.

2D-SAR and 3D-QSAR analyses for acetylcholinesterase inhibitors.

Alzheimer's disease (AD) accounts for almost three quarters of dementia patients and interferes people's normal life. Great progress has been made recently in the study of Acetylcholinesterase (AChE), known as one of AD's biomarkers. In this study, acetylcholinesterase inhibitors (AChEI) were collected to build a two-dimensional structure-activity relationship (2D-SAR) model and three-dimensional quantitative structure-activity relationship (3D-QSAR) model based on feature selection method combined with random forest. After calculation, the prediction accuracy of the 2D-SAR model was 89.63% by using the tenfold cross-validation test and 87.27% for the independent test set. Three cutting ways were employed to build 3D-QSAR models. A model with the highest [Formula: see text] (cross-validated correlation coefficient) and [Formula: see text](non-cross-validated correlation coefficient) was obtained to predict AChEI activity. The mean absolute error (MAE) of the training set and the test set was 0.0689 and 0.5273, respectively. In addition, molecular docking was also employed to reveal that the ionization state of the compounds had an impact upon their interaction with AChE. Molecular docking results indicate that Ser124 might be one of the active site residues.

Insight into the structural requirements of pyrimidine-based phosphodiesterase 10A (PDE10A) inhibitors by multiple validated 3D QSAR approaches

Schizophrenia is a complex disorder of thinking and behaviour (0.3−0.7% of the population is affected). The over-expression of phosphodiesterase 10A (PDE10A) enzyme may be a potential target for schizophrenia and Huntington’s disease. Because 3D QSAR analysis is one of the most frequently used modelling techniques, in the present study, five different 3D QSAR tools, namely CoMFA, CoMSIA, kNN-MFA, Open3DQSAR and topomer CoMFA methods, were used on a dataset of pyrimidine-based PDE10A inhibitors. All developed models were validated internally and externally. The non-commercial Open3DQSAR produced the best statistical results amongst 3D QSAR tools. The structural interpretations obtained from different methods were thoroughly analysed and were justified on the basis of information obtained from the crystal structure. Information from one method was mostly validated by the results of other methods and vice versa. In the current work, the use of multiple tools in the same analysis revealed more complete information about the structural requirements of these compounds. On the basis of the observations of the 3D QSAR studies, 12 new compounds were designed for better PDE10A inhibitory activity. The current investigation may help in further designing new PDE10A inhibitors with promising activity.

Molecular docking, 3D-QSAR and structural optimization on imidazo-pyridine derivatives dually targeting AT1 and PPARg.

Telmisartan, a bifunctional agent of blood pressure lowering and glycemia reduction, was previously reported to antagonize angiotensin II type 1 (AT1) receptor and partially activate peroxisome proliferator-activated receptor γ (PPARγ) simultaneously. Through the modification to telmisartan, researchers designed and obtained imidazo-\\pyridine derivatives with the IC50s of 0.49∼94.1 nM against AT1 and EC50s of 20∼3640 nM towards PPARγ partial activation. For minutely inquiring the interaction modes with the relevant receptor and analyzing the structure-activity relationships, molecular docking and 3D-QSAR (Quantitative structure-activity relationships) analysis of these imidazo-\\pyridines on dual targets were conducted in this work. Docking approaches of these derivatives with both receptors provided explicit interaction behaviors and excellent matching degree with the binding pockets. The best CoMFA (Comparative Molecular Field Analysis) models exhibited predictive results of q2=0.553, r2=0.954, SEE=0.127, r2pred=0.779 for AT1 and q2=0.503, r2=1.00, SEE=0.019, r2pred=0.604 for PPARγ, respectively. The contour maps from the optimal model showed detailed information of structural features (steric and electrostatic fields) towards the biological activity. Combining the bioisosterism with the valuable information from above studies, we designed six molecules with better predicted activities towards AT1 and PPARγ partial activation. Overall, these results could be useful for designing potential dual AT1 antagonists and partial PPARγ agonists.