Sundarapandian Thangapandian

3D QSAR pharmacophore based virtual screening and molecular docking for identification of potential HSP90 inhibitors.

Chemical features based 3D pharmacophore models were developed for HSP90 based on the known inhibitors using Discovery Studio V2.1. An optimal pharmacophore model was brought forth and validated using a decoy set, external test set and Fischers randomization method. The best five features pharmacophore model, Hypo1, includes two hydrogen bond acceptors, three hydrophobic features, which has the highest correlation coefficient (0.93), cost difference (73.88), low RMS (1.24), as well as it shows a high goodness of fit and enrichment factor. Hypo1 was used as a 3D query for virtual screening to retrieve potential inhibitors from Maybridge and Scaffold databases. The hit compounds were subsequently subjected to molecular docking studies and finally, 36 compounds were obtained based on consensus scoring function.

Potent bace-1 inhibitor design using pharmacophore modeling, in silico screening and molecular docking studies

BackgroundBeta-site amyloid precursor protein cleaving enzyme (BACE-1) is a single-membrane protein belongs to the aspartyl protease class of catabolic enzymes. This enzyme involved in the processing of the amyloid precursor protein (APP). The cleavage of APP by BACE-1 is the rate-limiting step in the amyloid cascade leading to the production of two peptide fragments Aβ40 and Aβ42. Among two peptide fragments Aβ42 is the primary species thought to be responsible for the neurotoxicity and amyloid plaque formation that lead to memory and cognitive defects in Alzheimer’s disease (AD). AD is a ravaging neurodegenerative disorder for which no disease-modifying treatment is currently available. Inhibition of BACE-1 is expected to stop amyloid plaque formation and emerged as an interesting and attractive therapeutic target for AD.MethodsLigand-based computational approach was used to identify the molecular chemical features required for the inhibition of BACE-1 enzyme. A training set of 20 compounds with known experimental activity was used to generate pharmacophore hypotheses using 3D QSAR Pharmacophore Generation module available in Discovery studio. The hypothesis was validated by four different methods and the best hypothesis was utilized in database screening of four chemical databases like Maybridge, Chembridge, NCI and Asinex. The retrieved hit compounds were subjected to molecular docking study using GOLD 4.1 program.ResultsAmong ten generated pharmacophore hypotheses, Hypo 1 was chosen as best pharmacophore hypothesis. Hypo 1 consists of one hydrogen bond donor, one positive ionizable, one ring aromatic and two hydrophobic features with high correlation coefficient of 0.977, highest cost difference of 121.98 bits and lowest RMSD value of 0.804. Hypo 1 was validated using Fischer randomization method, test set with a correlation coefficient of 0.917, leave-one-out method and decoy set with a goodness of hit score of 0.76. The validated Hypo 1 was used as a 3D query in database screening and retrieved 773 compounds with the estimated activity value <100 nM. These hits were docked into the active site of BACE-1 and further refined based on molecular interactions with the essential amino acids and good GOLD fitness score.ConclusionThe best pharmacophore hypothesis, Hypo 1, with high predictive ability contains chemical features required for the effective inhibition of BACE-1. Using Hypo 1, we have identified two compounds with diverse chemical scaffolds as potential virtual leads which, as such or upon further optimization, can be used in the designing of new BACE-1 inhibitors.

Ligand and structure based pharmacophore modeling to facilitate novel histone deacetylase 8 inhibitor design

Over expression of histone deacetylases (HDACs) leads to the suppression of various gene expressions including cancer suppressor gene. Thus, novel inhibitors of these enzymes can be a valid method to treat cancers. To facilitate the discovery of novel HDAC8 inhibitors, pharmacophore models were generated using ligand and receptor based approaches and validated with a database of active and inactive compounds. These validated pharmacophores have effectively been used in search of three databases and final hits were subjected to molecular docking using GOLD 4.1 program. Hit compounds that scored high GOLD fitness scores and showed interactions with catalytically important residues and metal ions were considered. Finally, three compounds have been reported as novel virtual leads to design potent HDAC8 inhibitors.

3D QSAR Pharmacophore Modeling, in Silico Screening, and Density Functional Theory (DFT) Approaches for Identification of Human Chymase Inhibitors

Human chymase is a very important target for the treatment of cardiovascular diseases. Using a series of theoretical methods like pharmacophore modeling, database screening, molecular docking and Density Functional Theory (DFT) calculations, an investigation for identification of novel chymase inhibitors, and to specify the key factors crucial for the binding and interaction between chymase and inhibitors is performed. A highly correlating (r = 0.942) pharmacophore model (Hypo1) with two hydrogen bond acceptors, and three hydrophobic aromatic features is generated. After successfully validating “Hypo1”, it is further applied in database screening. Hit compounds are subjected to various drug-like filtrations and molecular docking studies. Finally, three structurally diverse compounds with high GOLD fitness scores and interactions with key active site amino acids are identified as potent chymase hits. Moreover, DFT study is performed which confirms very clear trends between electronic properties and inhibitory activity (IC50) data thus successfully validating “Hypo1” by DFT method. Therefore, this research exertion can be helpful in the development of new potent hits for chymase. In addition, the combinational use of docking, orbital energies and molecular electrostatic potential analysis is also demonstrated as a good endeavor to gain an insight into the interaction between chymase and inhibitors.

Energy transduction and alternating access of the mammalian ABC transporter P-glycoprotein

ATP binding cassette (ABC) transporters of the exporter class harness the energy of ATP hydrolysis in the nucleotide-binding domains (NBDs) to power the energetically uphill efflux of substrates by a dedicated transmembrane domain (TMD). Although numerous investigations have described the mechanism of ATP hydrolysis and defined the architecture of ABC exporters, a detailed structural dynamic understanding of the transduction of ATP energy to the work of substrate translocation remains elusive. Here we used double electron–electron resonance and molecular dynamics simulations to describe the ATP- and substrate-coupled conformational cycle of the mouse ABC efflux transporter P-glycoprotein (Pgp; also known as ABCB1), which has a central role in the clearance of xenobiotics and in cancer resistance to chemotherapy. Pairs of spin labels were introduced at residues selected to track the putative inward-facing to outward-facing transition. Our findings illuminate how ATP energy is harnessed in the NBDs in a two-stroke cycle and elucidate the consequent conformational motion that reconfigures the TMD, two critical aspects of Pgp transport mechanism. Along with a fully atomistic model of the outward-facing conformation in membranes, the insight into Pgp conformational dynamics harmonizes mechanistic and structural data into a novel perspective on ATP-coupled transport and reveals mechanistic divergence within the efflux class of ABC transporters.

Pharmacophore based virtual screening, molecular docking studies to design potent heat shock protein 90 inhibitors

The identification of important chemical features of Heat Shock Protein 90 (HSP90) inhibitors will be helpful to discover the potent candidate to inhibit the HSP90 activity. The best hypothesis from Hip-Hop, Hypo1, one hydrogen bond donor (HBD), two hydrogen bond acceptors (HBA), and two hydrophobic (H) and structure-based hypothesis, SB_Hypo1, one HBA, one HBD and four H features, were generated using Discovery Studio and LigandScout, respectively. Test and decoy sets were used to corroborate the best hypotheses and the validated hypotheses were used to screen the chemical databases. Subsequently, the screened compounds were filtered by applying the rule of five, ADMET and molecular docking. Finally, four compounds were obtained as novel leads to inhibit the HSP90 activity.

Development, evaluation and application of 3D QSAR Pharmacophore model in the discovery of potential human renin inhibitors

BackgroundRenin has become an attractive target in controlling hypertension because of the high specificity towards its only substrate, angiotensinogen. The conversion of angiotensinogen to angiotensin I is the first and rate-limiting step of renin-angiotensin system and thus designing inhibitors to block this step is focused in this study.MethodsLigand-based quantitative pharmacophore modeling methodology was used in identifying the important molecular chemical features present in the set of already known active compounds and the missing features from the set of inactive compounds. A training set containing 18 compounds including active and inactive compounds with a substantial degree of diversity was used in developing the pharmacophore models. A test set containing 93 compounds, Fischer randomization, and leave-one-out methods were used in the validation of the pharmacophore model. Database screening was performed using the best pharmacophore model as a 3D structural query. Molecular docking and density functional theory calculations were used to select the hit compounds with strong molecular interactions and favorable electronic features.ResultsThe best quantitative pharmacophore model selected was made of one hydrophobic, one hydrogen bond donor, and two hydrogen bond acceptor features with high a correlation value of 0.944. Upon validation using an external test set of 93 compounds, Fischer randomization, and leave-one-out methods, this model was used in database screening to identify chemical compounds containing the identified pharmacophoric features. Molecular docking and density functional theory studies have confirmed that the identified hits possess the essential binding characteristics and electronic properties of potent inhibitors.ConclusionA quantitative pharmacophore model of predictive ability was developed with essential molecular features of a potent renin inhibitor. Using this pharmacophore model, two potential inhibitory leads were identified to be used in designing novel and future renin inhibitors as antihypertensive drugs.

Molecular Docking and Pharmacophore Filtering in the Discovery of Dual-Inhibitors for Human Leukotriene A4 Hydrolase and Leukotriene C4 Synthase

Combination of drugs for multiple targets has been a standard treatment in treating various diseases. A single chemical entity that acts upon multiple targets is emerging nowadays because of their predictable pharmacokinetic and pharmacodynamic properties. We have employed a computer-aided methodology combining molecular docking and pharmacophore filtering to identify chemical compounds that can simultaneously inhibit the human leukotriene hydrolase (hLTA4H) and the human leukotriene C4 synthase (hLTC4S) enzymes. These enzymes are the members of arachidonic acid pathway and act upon the same substrate, LTA4, producing different inflammatory products. A huge set of 4966 druglike compounds from the Maybridge database were docked into the active site of hLTA4H using the GOLD program. Common feature pharmacophore models were developed from the known inhibitors of both the targets using Accelrys Discovery Studio 2.5. The hits from the hLTA4H docking were filtered to match the chemical features of both the pharmacophore models. The compounds that resulted from the pharmacophore filtering were docked into the active site of hLTC4S and the hits those bind well at both the active sites and matched the pharmacophore models were identified as possible dual inhibitors for hLTA4H and hLTC4S enzymes. Reverse validation was performed to ensure the results of the study.

Potential virtual lead identification in the discovery of renin inhibitors: application of ligand and structure-based pharmacophore modeling approaches.

Renin, an enzyme by cleaving angiotensinogen to angiotensin-I, controls the first and rate-limiting step of renin-angiotensin system that is associated with blood pressure. Thus Ligand and structure-based pharmacophore models were developed in this study to identify new potential leads inhibiting this rate-limiting enzyme as an efficient way to treat blood pressure. X-ray predicted binding modes of most active compounds were used in ligand-based approach whereas the 3D structural information of renin was used in structure-based approach. Pharmacophore models were validated using various methods and utilized in database searching to identify potential hits. Drug-like filters and molecular docking studies led us identifying the final hits to be employed in designing new class of future renin inhibitors.

Construction and characterization of chimeric cellulases with enhanced catalytic activity towards insoluble cellulosic substrates

The chimeric proteins viz. CBM3-Cel9A, CBM4-Cel9A and CBM30-Cel9A, are constructed by fusion of family 3, 4, and 30 cellulose binding modules (CBMs) to N-terminus of family 9 endoglucanase (Cel9A) from Alicyclobacillus acidocaldrious. The chimeric enzymes were successfully expressed in Escherichia coli and purified to homogeneity. The chimeric enzymes showed significant increase in Avicel (8-12 folds) and filter paper (7-10 folds) degradation activities compared to Cel9A endoglucanase. Computational protein modeling and simulation on the chimeric enzymes were applied to analyze the fused CBMs effect on the increased insoluble cellulosic substrates degradation activity. Thin layer chromatography analysis of the enzymatic hydrolysis products and distribution of reducing sugars between soluble and insoluble fractions indicated processive cleavage of insoluble cellulosic substrates by the chimeras. The fused CBMs played a critical accessory role for the Cel9A catalytic domain and changed its character to facilitate the processive cleavage of insoluble cellulosic substrates.