Three-dimensional Quantitative Structure-activity Relationship Studies Research Articles

Synthesis and Biological Evaluation of Pyrazole-Pyrimidones as a New Class of Correctors of the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR).

Cystic fibrosis (CF) is caused by the functional expression defect of the cystic fibrosis transmembrane conductance regulator (CFTR) protein. Despite the recent success in CFTR modulator development, the available correctors only partially restore the F508del-CFTR channel function, and several rare CF mutations show resistance to available drugs. We previously identified compound 4172 that synergistically rescued the F508del-CFTR folding defect in combination with the existing corrector drugs VX-809 and VX-661. Here, novel CFTR correctors were designed by applying a classical medicinal chemistry approach on the 4172 scaffold. Molecular docking and three-dimensional quantitative structure-activity relationship (3D-QSAR) studies were conducted to propose a plausible binding site and design more potent and effective analogs. We identified three optimized compounds, which, in combination with VX-809 and the investigational corrector 3151, increased the plasma membrane density and function of F508del-CFTR and other rare CFTR mutants resistant to the currently approved therapies.

Three-Dimensional Quantitative Structure-Activity Relationship Study of Transient Receptor Potential Vanilloid 1 Channel Antagonists Reveals Potential for Drug Design Purposes.

Transient receptor potential vanilloid 1 (TRPV1) was reported to be a putative target for recovery from chronic pain, producing analgesic effects after its inhibition. A series of drug candidates were previously developed, without the ability to ameliorate the therapeutic outcome. Starting from previously designed compounds, derived from the hybridization of antagonist SB-705498 and partial agonist MDR-652, we performed a virtual screening on a pharmacophore model built by exploiting the Cryo-EM 3D structure of a nanomolar antagonist in complex with the human TRPV1 channel. The pharmacophore model was described by three pharmacophoric features, taking advantage of both the bioactive pose of the antagonist and the receptor exclusion spheres. The results of the screening were implemented inside a 3D-QSAR model, correlating with the negative decadic logarithm of the inhibition rate of the ligands. After the validation of the obtained 3D-QSAR model, we designed a new series of compounds by introducing key modifications on the original scaffold. Again, we determined the compounds' binding poses after alignment to the pharmacophoric model, and we predicted their inhibition rates with the validated 3D-QSAR model. The obtained values resulted in being even more promising than parent compounds, demonstrating that ongoing research still leaves much room for improvement.

Homologous Design and Three-Dimensional Quantitative Structure-Activity Relationship Study of Acaricidal 2,4-Diphenyloxazolines Containing Different Heteroatoms and Alkyl Chains.

In order to speculate the three-dimensional structure of the potential binding pocket of the chitin synthase inhibitor, a series of 2,4-diphenyloxazoline derivatives with different lengths of alkyl chains and heteroatoms were designed and synthesized by a homologous strategy. The bioassay results indicate that both the length of the alkyl chains and the type of substituents can affect the acaricidal activity against mite eggs. Compounds containing chloropropyl, alkoxyalkyl, and para-substituted phenoxyalkyl or phenylthioalkyl groups exhibit good activity, while those containing steric hindrance substituents or carbonyl substituents on the benzene ring exhibit reduced activity. Three-dimensional quantitative structure-activity relationship (3D-QSAR) study showed that there may be a narrow hydrophobic region deep in the pocket, and the steric effect plays a more important role than the electrostatic effect. The current work will provide assistance for future molecular design and target binding research.

Design, Synthesis, and Biological Activity of Novel Triketone-Containing Phenoxy Nicotinyl Inhibitors of HPPD.

4-Hydroxyphenylpyruvate dioxygenase (HPPD) is a crucial target enzyme in albino herbicides. The inhibition of HPPD activity interferes with the synthesis of carotenoids, blocking photosynthesis and resulting in bleaching and necrosis. To develop herbicides with excellent activity, a series of 3-hydroxy-2-(6-substituted phenoxynicotinoyl)-2-cyclohexen-1-one derivatives were designed via active substructure combination. The title compounds were characterized via infrared spectroscopy, 1H and 13C nuclear magnetic resonance spectroscopies, and high-resolution mass spectrometry. The structure of compound III-17 was confirmed via single-crystal X-ray diffraction. Preliminary tests demonstrated that some compounds had good herbicidal activity. Crop safety tests revealed that compound III-29 was safer than the commercial herbicide mesotrione in wheat and peanuts. Moreover, the compound exhibited the highest inhibitory activity against Arabidopsis thaliana HPPD (AtHPPD), with a half-maximal inhibitory concentration of 0.19 μM, demonstrating superior activity compared with mesotrione (0.28 μM) in vitro. A three-dimensional quantitative structure-activity relationship study revealed that the introduction of smaller groups to the 5-position of cyclohexanedione and negative charges to the 3-position of the benzene ring enhanced the herbicidal activity. A molecular structure comparison demonstrated that compound III-29 was beneficial to plant absorption and conduction. Molecular docking and molecular dynamics simulations further verified the stability of the complex formed by compound III-29 and AtHPPD. Thus, this study may provide insights into the development of green and efficient herbicides.

Novel 5-Sulfonyl-1,3,4-thiadiazole-Substituted Flavonoids as Potential Bactericides and Fungicides: Design, Synthesis, Three-Dimensional Quantitative Structure-Activity Relationship Studies.

Flavonoids, ubiquitous natural products, provide sources for drug discovery owing to their structural diversity, broad-spectrum pharmacological activity, and excellent environmental compatibility. To develop antibacterial and antifungal agents with novel mechanisms of action and innovative structures, a series of novel 5-sulfonyl-1,3,4-thiadiazole-substituted flavonoids were designed and synthesized, and their biological activities against seven agriculturally common phytopathogenic microorganisms were evaluated. The results of the antimicrobial bioassay showed that most of the target compounds displayed excellent inhibitory effects against Xanthomonas oryzae, Rhizoctonia solani, and Colletotrichum orbiculare. Compounds 1, 3, 7, 9, 13, and 14 exhibited remarkable antibacterial activity against X. oryzae pv. oryzae with EC50 values below 10 μg/mL, which were superior to bismerthiazol (70.89 μg/mL). Compound 2 (EC50 = 0.41 μg/mL) displayed the most effective inhibitory potency against R. solani in vivo, comparable protective effects with the positive control carbendizam. Preliminary mechanistic studies indicated that compound 2 induced disordered entanglement of hyphae, shrinkage of hyphal surfaces, extravasation of cellular contents, and vacuole swelling and rupture, which disrupted normal hyphal growth. Subsequently, compounds 35-53 with good antifungal activity were designed and synthesized based on reliable three-dimensional quantitative structure-activity relationship (3D-QSAR) models. Compound 49 showed high efficacy and superior antifungal activity against R. solani, with an EC50 value of 0.28 μg/mL and a half-maximal effective concentration of 0.46 μg/mL.

3D-QSAR, ADMET, and molecular docking studies of aztreonam analogs as E. colis inhibitors.

The development of multidrug resistant strains of extended-spectrum β-lactamase-producing Escherichia coli has become a global problem; therefore, the discovery of new antibacterial agents is the only available solution. To improve and propose new compounds with antibacterial activity, the three-dimensional quantitative structure-activity relationship and molecular docking studies were carried out on Aztreonam analogs as E. coli inhibitors in DNA gyrase B. This study's 3D-Quantitative structure-activity relationship model was created using on the Comparative Molecular Field Analysis and the Comparative Molecular Similarity Indices Analysis. Using the Comparative Molecular Field Analysis (Q 2 = 0.73; R 2 = 0.82), excellent predictability was achieved, and the best Comparative Molecular Similarity Indices Analysis model (Q 2 = 0.88; R 2 = 0.9). The generated model's ability to predict outcomes was assessed through external validation using a test set compound and an applicability domain technique. In this study, the steric, electrostatic, and hydrogen bond acceptor fields played a key role in antibacterial activity. The results of the molecular docking revealed that the newly generated compound A6 has the highest binding affinity with DNA gyrase B. It forms 10 hydrogen bonds with amino acid residues of Asn104, Asn274, Asn132, Ser70, Ser237, Thr105, Glu273, and 2 salt bridges with amino acid residues of Ser70 and Glu273 and one pi-pi interacting with Gys271 amino acid residue in the binding site of 5G1, and this result was validated by a new assessment method. We created some novel, highly effective DNA gyrase B inhibitors based on the earlier findings, and the most accurate model predicted their inhibitory actions. The ADMET characteristics and pharmacological similarity of these novel inhibitors were also examined. These findings would be very beneficial in guiding the optimization process for the identification of novel drugs that can address the issue of multiple drug resistance.

Investigating the common pharmacophoric points of PDK1 inhibitors as anti-cancer agents using an alignment independent 3D-QSAR, molecular docking and molecular dynamic simulation

The glycolytic enzyme 3-Phosphoinositide-dependent protein kinase-1 (PDK1), a member of the serine/threonine kinase family, is a promising target for developing novel anticancer drugs due to its vital role in growth and survival of cancer cells and tumor angiogenesis. A MIF base alignment independent, Three-Dimensional Quantitative Structure-Activity Relationship (3D-QSAR) study was conducted to highlight the correlation between PDK1 inhibitors and their structural characteristics, thus indicating the most effective molecular properties of the inhibitors as anticancer agents. Using the Kennard-Stone algorithm to separate calibration and validation subsets, Partial Least-Square (PLS) method application shows excellent agreement between experimental and predicted values (R2=0.85, Q2LOO=0.77 and R2Perd=0.93), making it a reliable model for future drug design approaches. In addition, the active site of the PDK1 receptor protein was examined via Molecular Docking with both passive and active compounds. In order to study the effect of inhibitors on PDK1 structure, 10 ns Molecular Dynamics (MD) simulation of free protein and protein-ligand complex was performed and compared. Comparing the 3D-QSAR, docking and molecular simulation outcomes, the prevailing pharmacophoric use of PDK1 inhibitors as anti-cancer agents have been highlighted. The comprehensive information of this approach reveals important structural variables and their relationship with PDK1 inhibitory activity and PDK1 cavity residues which could be valuable for designing anticancer candidates.

Preclinical Evaluation of an Imidazole-Linked Heterocycle for Alzheimer's Disease.

Humanity is facing a vast prevalence of neurodegenerative diseases, with Alzheimer's disease (AD) being the most dominant, without efficacious drugs, and with only a few therapeutic targets identified. In this scenario, we aim to find molecular entities that modulate imidazoline I2 receptors (I2-IRs) that have been pointed out as relevant targets in AD. In this work, we explored structural modifications of well-established I2-IR ligands, giving access to derivatives with an imidazole-linked heterocycle as a common key feature. We report the synthesis, the affinity in human I2-IRs, the brain penetration capabilities, the in silico ADMET studies, and the three-dimensional quantitative structure-activity relationship (3D-QSAR) studies of this new bunch of I2-IR ligands. Selected compounds showed neuroprotective properties and beneficial effects in an in vitro model of Parkinson's disease, rescued the human dopaminergic cell line SH-SY5Y from death after treatment with 6-hydroxydopamine, and showed crucial anti-inflammatory effects in a cellular model of neuroinflammation. After a preliminary pharmacokinetic study, we explored the action of our representative 2-(benzo[b]thiophen-2-yl)-1H-imidazole LSL33 in a mouse model of AD (5xFAD). Oral administration of LSL33 at 2 mg/Kg for 4 weeks ameliorated 5XFAD cognitive impairment and synaptic plasticity, as well as reduced neuroinflammation markers. In summary, this new I2-IR ligand that promoted beneficial effects in a well-established AD mouse model should be considered a promising therapeutic strategy for neurodegeneration.

Binding Studies and Lead Generation of Pteridin-7(8H)-one Derivatives Targeting FLT3.

Ligand modification by substituting chemical groups within the binding pocket is a popular strategy for kinase drug development. In this study, a series of pteridin-7(8H)-one derivatives targeting wild-type FMS-like tyrosine kinase-3 (FLT3) and its D835Y mutant (FL3D835Y) were studied using a combination of molecular modeling techniques, such as docking, molecular dynamics (MD), binding energy calculation, and three-dimensional quantitative structure-activity relationship (3D-QSAR) studies. We determined the protein–ligand binding affinity by employing molecular mechanics Poisson–Boltzmann/generalized Born surface area (MM-PB/GBSA), fast pulling ligand (FPL) simulation, linear interaction energy (LIE), umbrella sampling (US), and free energy perturbation (FEP) scoring functions. The structure–activity relationship (SAR) study was conducted using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA), and the results were emphasized as a SAR scheme. In both the CoMFA and CoMSIA models, satisfactory correlation statistics were obtained between the observed and predicted inhibitory activity. The MD and SAR models were co-utilized to design several new compounds, and their inhibitory activities were anticipated using the CoMSIA model. The designed compounds with higher predicted pIC50 values than the most active compound were carried out for binding free energy evaluation to wild-type and mutant receptors using MM-PB/GBSA, LIE, and FEP methods.

An alignment-independent three-dimensional quantitative structure-activity relationship study on ron receptor tyrosine kinase inhibitors.

Recepteur d'Origine Nantais known as RON is a member of the receptor tyrosine kinase (RTK) superfamily which has recently gained increasing attention as cancer target for therapeutic intervention. The aim of this work was to perform an alignment-independent three-dimensional quantitative structure-activity relationship (3D QSAR) study for a series of RON inhibitors. A 3D QSAR model based on GRid-INdependent Descriptors (GRIND) methodology was generated using a set of 19 compounds with RON inhibitory activities. The generated 3D QSAR model revealed the main structural features important in the potency of RON inhibitors. The results obtained from the presented study can be used in lead optimization projects for designing of novel compounds where inhibition of RON is needed.

Development of novel monoamine oxidase B (MAO-B) inhibitors by combined application of docking-based alignment, 3D-QSAR, ADMET prediction, molecular dynamics simulation, and MM_GBSA binding free energy

Unsaturated ketone derivatives are known as inhibitors of monoamine oxidase B (MAO-B), a potential drug target of Parkinson’s disease. Here, docking-based alignment, 3 D-QSAR (three-dimensional quantitative structure-activity relationship) studies, ADMET (absorption, distribution, metabolism, excretion, and toxicity) prediction, molecular dynamics (MD) simulation, and MM_GBSA binding free energy were performed on a novel series of MAO-B inhibitors. The objective is to predict new MAO-B inhibitors with high potency activity. The 3 D-QSAR models were created using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA). Molecular docking findings indicated that compounds with strong inhibitory efficacy also had a high binding affinity. 3 D-QSAR studies showed the importance of steric, electrostatic, and H-bond acceptor fields on the inhibitory activity of MAO-B. Based on the appropriate 3 D-QSAR model, a new series of MAO-B inhibitors were predicted and their pharmacokinetic characteristics were evaluated using in silico ADMET prediction. All screened compounds show good oral bioavailability without any side effects. Moreover, the dynamic behavior and stability of the most active compounds were evaluated using MD simulations. The results showed that unsaturated ketone derivatives are stable and compact during the 100 ns of MD simulation. Finally, the binding free energy of complexes was determined using the MM_GBSA method; the findings indicated that the T1 compound is more stable (ΔGbinding = −409.506 KJ/mol) than the data set's highest active compound (ΔGbinding = −31.883 KJ/mol). Communicated by Ramaswamy H. Sarma

Synthesis, Antifungal Activity, Three-Dimensional Quantitative Structure-Activity Relationship and Molecular Docking Study of 4-Acyl-3-amino-1,2,4-triazole-thioether Derivatives Containing Natural Pinene Structure

Synthesis, Antifungal Activity, Three-Dimensional Quantitative Structure-Activity Relationship and Molecular Docking Study of 4-Acyl-3-amino-1,2,4-triazole-thioether Derivatives Containing Natural Pinene Structure

3D-QSAR, molecular docking, DFT and ADMET studies on quinazoline derivatives to explore novel DHFR inhibitors

Resistance to folate antagonists is caused by mutations in the dihydrofolate reductase (DHFR) genes. These mutations affect the amino acids at positions 51, 59, 108 and 164 of DHFR, which appear to play a major role in malaria treatment failure. Therefore, the design of new drugs able to overcome the problem of antifolate drug resistance should receive urgent attention. In this study, a three-dimensional quantitative structure-activity relationship (3 D-QSAR) and molecular docking studies have been performed on antimalarial quinazoline derivatives. The CoMFA (Q2 = 0.63, R2 = 0.83 and = 0.70) and the CoMSIA (Q2 = 0.584, R2 = 0.816, and = 0.73) models show a good prediction of antimalarial activity. The reliability and robustness of the proposed models have been tested using several validation methods, which showed that the steric, electrostatic, hydrophobic and H-bond acceptor fields of the CoMSIA model play a key role in the prediction of antimalarial activity. Molecular docking studies reveal important interactions between two isomeric compounds (meta and para) and the DHFR receptor in its wild and mutant forms. The obtained outcomes of molecular docking studies have been validated using a new method based on visual inspection. The DFT study of the two isomeric compounds confirms clearly the trends of 3 D-QSAR and molecular docking for the design of new compounds. Moreover, the consistency between theoretical, 3 D-QSAR and molecular docking analysis provides guidance for the design of new drug candidates, which have been tested using ADMET properties and drug likeness analysis. Communicated by Ramaswamy H. Sarma

Human Estrogen Receptor α Antagonists. Part 1: 3-D QSAR-Driven Rational Design of Innovative Coumarin-Related Antiestrogens as Breast Cancer Suppressants through Structure-Based and Ligand-Based Studies.

The estrogen receptor α (ERα) represents a 17β-estradiol-inducible transcriptional regulator that initiates the RNA polymerase II-dependent transcriptional machinery, pointed for breast cancer (BC) development via either genomic direct or genomic indirect (i.e., tethered) pathway. To develop innovative ligands, structure-based (SB) three-dimensional (3-D) quantitative structure-activity relationship (QSAR) studies have been undertaken from structural data taken from partial agonists, mixed agonists/antagonists (selective estrogen receptor modulators (SERMs)), and full antagonists (selective ERα downregulators (SERDs)) correlated with either wild-type or mutated ERα receptors. SB and ligand-based (LB) alignments allow us to rule out guidelines for the SB/LB alignment of untested compounds. 3-D QSAR models for ERα ligands, coupled with SB/LB alignment, were revealed to be useful tools to dissect the chemical determinants for ERα-based anticancer activity as well as to predict their potency. The herein developed protocol procedure was verified through the design and potency prediction of 12 new coumarin-based SERMs, namely, 3DQ-1a to 3DQ-1e, that upon synthesis turned to be potent ERα antagonists by means of either in vitro or in vivo assays (described in the second part of this study).

Integrated 3D-QSAR, molecular docking, and molecular dynamics simulation studies on 1,2,3-triazole based derivatives for designing new acetylcholinesterase inhibitors.

Alzheimer’s disease (AD) is a multifactorial and polygenic disease. It is the most prevalent reason for dementia in the aging population. A dataset of twenty-six 1,2,3-triazole-based derivatives previously synthetized and evaluated for acetylcholinesterase inhibitory activity were subjected to the three-dimensional quantitative structure-activity relationship (3D-QSAR) study. Good predictability was achieved for comparative molecular field analysis (CoMFA) (Q2 = 0.604, R2 = 0.863, rext2 = 0.701) and comparative molecular similarity indices analysis (CoMSIA) (Q2 = 0.606, R2 = 0.854, rext2 = 0.647). The molecular features characteristics provided by the 3D-QSAR contour plots were quite useful for designing and improving the activity of acetylcholinesterase of this class. Based on these findings, a new series of 1,2,3-triazole based derivatives were designed, among which compound A1 with the highest predictive activity was subjected to detailed molecular docking and compared to the most active compound. The selected compounds were further subjected to 20 ns molecular dynamics (MD) simulations to study the comparative conformation dynamics of the protein after ligand binding, revealing promising results for the designed molecule. Therefore, this study could provide worthy guidance for further experimental analysis of highly effective acetylcholinesterase inhibitors.

Molecular Docking, 3D-QSAR, Fingerprint-Based 2D-QSAR, Analysis of Pyrimidine, and Analogs of ALK (Anaplastic Lymphoma Kinase) Inhibitors as an Anticancer Agent

Background: ALK inhibitors have become a plausible option for anticancer therapy with the availability of several FDA-approved molecules and clinical trial candidates. Hence, the design of new ALK inhibitors using computational molecular docking studies on the existing inhibitors, is an attractive approach for anticancer drug discovery. Methods: We generated six types of independent models through structural based molecular docking study, three-dimensional quantitative structure-activity relationship (3D-QSAR) study, and 2DQSAR approaches using different fingerprints, such as dendritic, linear, 2D molprint, and radial. Results: Comparison of the generated models showed that the hinge region hydrogen bond interacted with amino acids ASP1206, MET1199, and LYS1150 in docking analysis and the hydrophobic interacted with amino acids GLU1210, ARG1209, SER1206, and LYS1205 residues are responsible for the ALK inhibition. In the 3D-QSAR study, the hydrogen bond donor features of 2,4- diaryl aminopyrimidine substituents, isopropyl phenyl ring groups in hydrophobic features, and electron-withdrawing groups matched the generated contour plots. The 2D-QSAR fingerprint studies indicated that higher potency was associated with the 2-hydroxy-5-isopropyl benzamide functional group and substituted phenylamine at the second position of the pyrimidine group. Conclusion: We conclude that the incorporation of these functional groups in the design of new molecules may result in more potent ALK inhibitors.

Anthelmintic efficacy of natural saponins against Gyrodactylus kobayashii in goldfish (Carassius auratus) and their 3D-QSAR analysis.

Gyrodactylus spp. are common monogenean ectoparasites that may lead to significant fish mortality. To find effective anthelmintic agents with lower toxicity, a series of natural saponins were obtained and evaluated for their anthelmintic activity against Gyrodactylus kobayashii and acute toxicity to goldfish (Carassius auratus). Among all tested compounds, six compounds (1, 2, 3, 8, 10, and 13) shown higher anthelmintic activity and safety than widely used formaldehyde-based parasiticides, especially compound 1 having 100% anthelmintic efficacy against G. kobayashii at 0.3mg/L and a therapeutic index of 16.6. Also, the three-dimensional quantitative structure-activity relationship (3D-QSAR) studies of these saponins have been performed to explore the structural features reasonable for the anthelmintic activity against G. kobayashii. These models demonstrated that the hydroxyl group at C-17 position and the sugar moieties at C-3 position, especially the hydroxyl groups of the sugar moieties, were critical to the anthelmintic activity. The QSAR studies could provide useful information for further rational design and optimization of novel saponins for the control of gyrodactylosis.

Three-dimensional quantitative structure-activity relationship study on gas chromatographic retention index of the fragrance compounds of Liliumspp

Quantitative structure-activity relationship models on the gas chromatographic retention index of the 38 volatile fragrance compounds of Liliumspp were investigated and established by comparative molecular field analysis (CoMFA) and comparative molecular similarity index (CoMSIA) methods. The robustness and predictive performance of the developed models were assessed using external test set validation and leave-one-out cross validation. Further, the effects of the molecular structure on the gas chromatographic retention indices of these compounds were intuitively studied in light of the three-dimensional contour maps of molecular fields provided by the developed CoMSIA and CoMFA models. The validation results demonstrated that both the models could accurately predict the retention indices of the investigated components. The influence of the molecular structure on the retention indices could be reasonably explained by these models. Moreover, the prediction accuracy of the CoMSIA model was slightly higher than that of the CoMFA model. Obviously, the proposed CoMSIA model is more promising for the analysis of the volatile fragrance compounds of Lilium spp.

Molecular-Level Understanding of the Somatostatin Receptor 1 (SSTR1)-Ligand Binding: A Structural Biology Study Based on Computational Methods.

Somatostatin receptor 1 (SSTR1), a subtype of somatostatin receptors, is involved in various signaling mechanisms in different parts of the human body. Like most of the G-protein-coupled receptors (GPCRs), the available information on the structural features of SSTR1 responsible for the biological activity is scarce. In this study, we report a molecular-level understanding of SSTR1–ligand binding, which could be helpful in solving the structural complexities involved in SSTR1 functioning. Based on a three-dimensional quantitative structure–activity relationship (3D-QSAR) study using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA), we have identified that an electronegative, less-bulkier, and hydrophobic atom substitution can substantially increase the biological activity of SSTR1 ligands. A density functional theory (DFT) study has been followed to study the electron-related properties of the SSTR1 ligands and to validate the results obtained via the 3D-QSAR study. 3D models of SSTR1–ligand systems have been embedded in lipid–lipid bilayer membranes to perform molecular dynamics (MD) simulations. Analysis of the MD trajectories reveals important information about the crucial residues involved in SSTR1–ligand binding and various conformational changes in the protein that occur after ligand binding. Additionally, we have identified the probable ligand-binding site of SSTR1 and validated it using MD. We have also studied the favorable conditions that are essential for forming the most stable and lowest-energy bioactive conformation of the ligands inside the binding site. The results of the study could be useful in constructing more potent and novel SSTR1 antagonists and agonists.

Three-dimensional Quantitative Structure-activity Relationship (3DQSAR) and Molecular Docking Study of 2-((pyridin-3-yloxy)methyl) Piperazines as α7 Nicotinic Acetylcholine Receptor Modulators for the Treatment of Inflammatory Disorders.

Comparative molecular field analysis (CoMFA) of 27 analogues of 2-((pyridin-3-yloxy)methyl)piperazine derivatives was carried out using software Tripos SYBYL X. Optimal r2 (0.854) and q2 (0.541) values were obtained for the developed 3D-QSAR model. The contour plots obtained from CoMFA analysis have shown 13.84% steric contribution and 66.14% electrostatic contribution towards an anti-inflammatory activity. The homology model of the receptor protein, α7 nicotinic acetylcholine, was generated in SWISS MODELLER using auto template mode and was analysed for the quality using Procheck, QMEAN Z-score, Anolea and GROMOS plots. The QMEAN score for the model was observed to be - 3.862. The generated model of alpha 7 nicotinic acetylcholine receptor was used for docking study of 27 piperazine analogues using Auto-Dock 4.2.5.1. The dock score obtained from docking analysis was then correlated with experimental pIC50 values for in-silico validation of the developed CoMFA model and a good correlation was obtained with correlation coefficient (r2) value of -0.7378. The present investigation suggests an optimal 3D-QSAR with CoMFA model for further evaluating new chemical entities based on piperazine skeleton.