Human leukotriene A4 hydrolase (hLTA4H), which is the final and rate-limiting enzyme of arachidonic acid pathway, converts the unstable epoxide LTA4 to a proinflammatory lipid mediator LTB4 through its hydrolase function. The LTA4H is a bi-functional enzyme that also exhibits aminopeptidase activity with a preference over arginyl tripeptides. Various mutations including E271Q, R563A, and K565A have completely or partially abolished both the functions of this enzyme. The crystal structures with these mutations have not shown any structural changes to address the loss of functions. Molecular dynamics simulations of LTA4 and tripeptide complex structures with functional mutations were performed to investigate the structural and conformation changes that scripts the observed differences in catalytic functions. The observed protein-ligand hydrogen bonds and distances between the important catalytic components have correlated well with the experimental results. This study also confirms based on the structural observation that E271 is very important for both the functions as it holds the catalytic metal ion at its location for the catalysis and it also acts as N-terminal recognition residue during peptide binding. The comparison of binding modes of substrates revealed the structural changes explaining the importance of R563 and K565 residues and the required alignment of substrate at the active site. The results of this study provide valuable information to be utilized in designing potent hLTA4H inhibitors as anti-inflammatory agents.

Angiogenesis is defined as the formation of new blood vessels and is a key phenomenon manifested in a host of cancers during which tyrosine kinases play a crucial role. Vascular endothelial growth factor receptor-2 (VEGFR-2) is pivotal in cancer angiogenesis, which warrants the urgency of discovering new anti-angiogenic inhibitors that target the signalling pathways. To obtain this objective, a structure-based pharmacophore model was built from the drug target VEGFR-2 (PDB code: 4AG8), complexed with axitinib and was subsequently validated and employed as a 3D query to retrieve the candidate compounds with the key inhibitory features. The model was escalated to molecular docking studies resulting in seven candidate compounds. The molecular docking studies revealed that the seven compounds displayed a higher dock score than the reference-cocrystallised compound. The GROningen MAchine for Chemical Simulations (GROMACS) package guided molecular dynamics (MD) results determined their binding mode and affirmed stable root mean square deviation. Furthermore, these compounds have preserved their key interactions with the residues Glu885, Glu917, Cys919 and Asp1046. The obtained findings deem that the seven compounds could act as novel anti-angiogenic inhibitors and may further assist as the prototype in designing and developing new inhibitors.

Progeria is a globally noticed rare genetic disorder manifested by premature aging with no effective treatment. Under these circumstances, farnesyltransferase inhibitors (FTIs) are marked as promising drug candidates. Correspondingly, a pharmacophore model was generated exploiting the features of lonafarnib. The selected pharmacophore model was allowed to screen the InterBioScreen natural compound database to retrieve the potential lead candidates. A series of filtering steps were applied to assess the drug-likeness of the compounds. The obtained compounds were advanced to molecular docking employing the CDOCKER module available with Discovery Studio (DS). Subsequently, three compounds (Hits) have displayed a higher dock score and demonstrated key residue interactions with stable molecular dynamics simulation results compared to the reference compound. Taken together, we therefore put forth three identified Hits as FTIs that may further serve as chemical spaces in designing new compounds.

The α-glucosidase inhibitory potential of Tribulus terrestris extracts has been reported but as yet the active ingredients are unknown. This study attempted to isolate the responsible metabolites and elucidate their inhibition mechanism of α-glucosidase. By fractionating T. terristris extracts, three cinnamic acid amide derivatives (1-3) were ascertained to be active components against α-glucosidase. The lead structure, N-trans-coumaroyltyramine 1, showed significant inhibition of α-glucosidase (IC50 = 0.42 μM). Moreover, all active compounds displayed uncompetitive inhibition mechanisms that have rarely been reported for α-glucosidase inhibitors. This kinetic behavior was fully demonstrated by showing a decrease of both Km and Vmax, and Kik/Kiv ratio ranging between 1.029 and 1.053. We progressed to study how chemical modifications to the lead structure 1 may impact inhibition. An α, β-unsaturation carbonyl group and hydroxyl group in A-ring of cinnamic acid amide emerged to be critical functionalities for α-glucosidase inhibition. The molecular modeling study revealed that the inhibitory activities are tightly related to π-π interaction as well as hydrogen bond interaction between enzyme and inhibitors.

Employing iPrMgCl as an advanced base instead of lithium hexamethyldisilazane (LHMDS) resulted in dramatic improvements in aza-Claisen rearrangement. This advance is considered responsible for the increased bulkiness of the alkoxide moiety (including magnesium cation and ligands), followed by a resultant conformational change of the transition state. To support this hypothesis, various substrates of aza-Claisen rearrangement were prepared and screened. In addition, a molecular dynamic simulation study was performed to investigate and compare the structural stability of reaction intermediates.

Progeria is a rare genetic disorder characterized by premature aging that eventually leads to death and is noticed globally. Despite alarming conditions, this disease lacks effective medications; however, the farnesyltransferase inhibitors (FTIs) are a hope in the dark. Therefore, the objective of the present article is to identify new compounds from the databases employing pharmacophore based virtual screening. Utilizing nine training set compounds along with lonafarnib, a common feature pharmacophore was constructed consisting of four features. The validated Hypo1 was subsequently allowed to screen Maybridge, Chembridge, and Asinex databases to retrieve the novel lead candidates, which were then subjected to Lipinski's rule of 5 and ADMET for drug-like assessment. The obtained 3,372 compounds were forwarded to docking simulations and were manually examined for the key interactions with the crucial residues. Two compounds that have demonstrated a higher dock score than the reference compounds and showed interactions with the crucial residues were subjected to MD simulations and binding free energy calculations to assess the stability of docked conformation and to investigate the binding interactions in detail. Furthermore, this study suggests that the Hits may be more effective against progeria and further the DFT studies were executed to understand their orbital energies.

Chagas disease is one of the primary causes of heart diseases accounting to 50,000 lives annually and is listed as the neglected tropical disease. Because the currently available therapies have greater toxic effects with higher resistance, there is a dire need to develop new drugs to combat the disease. In this pursuit, the 3D QSAR ligand-pharmacophore (pharm 1) and receptor-based pharmacophore (pharm 2) search was initiated to retrieve the candidate compounds from universal natural compounds database. The validated models were allowed to map the universal natural compounds database. The obtained lead candidates were subjected to molecular docking against cysteine protease (PDB code: 1ME3) employing -Cdocker available on the discovery studio. Subsequently, two Hits have satisfied the selection criteria and were escalated to molecular dynamics simulation and binding free energy calculations. These Hits have demonstrated higher dock scores, displayed interactions with the key residues portraying an ideal binding mode complemented by mapping to all the features of pharm 1 and pharm 2. Additionally, they have rendered stable root mean square deviation (RMSD) and potential energy profiles illuminating their potentiality as the prospective antichagastic agents. The study further demonstrates the mechanism of inhibition by tetrad residues compromising of Gly23 and Asn70 holding the ligand at each ends and the residues Gly65 and Gly160 clamping the Hits at the center. The notable feature is that the Hits lie in close proximity with the residues Glu66 and Leu67, accommodating within the S1, S2 and S3 subsites. Considering these findings, the study suggests that the Hits may be regarded as effective therapeutics against Chagas disease.

Breast cancer is one of the leading causes of death noticed in women across the world. Of late the most successful treatments rendered are the use of aromatase inhibitors (AIs). In the current study, a two-way approach for the identification of novel leads has been adapted. 81 chemical compounds were assessed to understand their potentiality against aromatase along with the four known drugs. Docking was performed employing the CDOCKER protocol available on the Discovery Studio (DS v4.5). Exemestane has displayed a higher dock score among the known drug candidates and is labeled as reference. Out of 81 ligands 14 have exhibited higher dock scores than the reference. In the second approach, these 14 compounds were utilized for the generation of the pharmacophore. The validated four-featured pharmacophore was then allowed to screen Chembridge database and the potential Hits were obtained after subjecting them to Lipinski's rule of five and the ADMET properties. Subsequently, the acquired 3,050 Hits were escalated to molecular docking utilizing GOLD v5.0. Finally, the obtained Hits were consequently represented to be ideal lead candidates that were escalated to the MD simulations and binding free energy calculations. Additionally, the gene-disease association was performed to delineate the associated disease caused by CYP19A1.

Bacterial peptide deformylase (PDF) is an attractive target for developing novel inhibitors against several types of multidrug-resistant bacteria. The objective of the current study is to retrieve potential phytochemicals as prospective drugs against Staphylococcus aureus peptide deformylase (SaPDF). The current study focuses on applying ligand-based pharmacophore model (PharmL) and receptor-based pharmacophore (PharmR) approaches. Utilizing 20 known active compounds, pharmL was built and validated using Fischer's randomization, test set method and the decoy set method. PharmR was generated from the knowledge imparted by the Interaction Generation protocol implemented on the Discovery Studio (DS) v4.5 and was validated using the decoy set that was employed for pharmL. The selection of pharmR was performed based upon the selectivity score and further utilizing the Pharmacophore Comparison module available on the DS. Subsequently, the validated pharmacophore models were escalated for Taiwan Indigenous Plants (TIP) database screening and furthermore, a drug-like evaluation was performed. Molecular docking was initiated for the resultant compounds, employing CDOCKER (available on the DS) and GOLD. Eventually, the stability of the final PDF⁻hit complexes was affirmed using molecular dynamics (MD) simulation conducted by GROMACS v5.0.6. The redeemed hits demonstrated a similar binding mode and stable intermolecular interactions with the key residues, as determined by no aberrant behaviour for 50 ns. Taken together, it can be stated that the hits can act as putative scaffolds against SaPDF, with a higher therapeutic value. Furthermore, they can act as fundamental structures for designing new drug candidates.

Cervical cancer is regarded as one of the major burdens noticed in women next to breast cancer. Although, human papilloma viruses (HPVs) are regarded as the principal causative agents, they require certain other factors such as oestrogen hormone to induce cervical cancer. Aromatase is an enzyme that converts androgens into oestrogens and hindering this enzyme could subsequently hamper the formation of oestrogen thereby alleviating the disease. Accordingly, in the current investigation, a structure based pharmacophore was generated considering two proteins bearing the Protein Data Bank (PDB) codes 3EQM (pharm 1) and 3S7S (pharm 2), respectively. The two models were employed as the 3D query to screen the in-house built natural compounds database. The obtained 51 compounds were escalated to molecular docking studies to decipher on the binding affinities and to predict the quintessential binding modes which were affirmed by molecular dynamics (MD) simulations. The compound has induced dose-dependent down regulation of PP2B, Nitric oxide synthase-2 (NOS2), and Interleukin 6 (IL-6) genes in the HeLa cells and has modulated the expression of apoptotic genes such as Bax, Bcl2, and caspases-3 at different concentrations. These results guide us to comprehend that the identified aromatase inhibitor was effective against the cervical cancer cells and additionally could server as scaffolds in designing new drugs.

Bruton's tyrosine kinase (BTK) is a cytoplasmic, non-receptor tyrosine kinase which is expressed in most of the hematopoietic cells and plays an important role in many cellular signaling pathways. B cell malignancies are dependent on BCR signaling, thus making BTK an efficient therapeutic target. Over the last few years, significant efforts have been made in order to develop BTK inhibitors to treat B-cell malignancies, and autoimmunity or allergy/hypersensitivity but limited success has been achieved. Here in this study, 3D QSAR pharmacophore models were generated for Btk based on known IC50 values and experimental energy scores with extensive validations. The five features pharmacophore model, Hypo1, includes one hydrogen bond acceptor lipid, one hydrogen bond donor, and three hydrophobic features, which has the highest correlation coefficient (0.98), cost difference (112.87), and low RMS (1.68). It was further validated by the Fisher's randomization method and test set. The well validated Hypo1 was used as a 3D query to search novel Btk inhibitors with different chemical scaffold using high throughput virtual screening technique. The screened compounds were further sorted by applying ADMET properties, Lipinski's rule of five and molecular docking studies to refine the retrieved hits. Furthermore, molecular dynamic simulation was employed to study the stability of docked conformation and to investigate the binding interactions in detail. Several important hydrogen bonds with Btk were revealed, which includes the gatekeeper residues Glu475 and Met 477 at the hinge region. Overall, this study suggests that the proposed hits may be more effective inhibitors for cancer and autoimmune therapy.

Histone deacetylase 8 (HDAC8) is an enzyme involved in deacetylating the amino groups of terminal lysine residues, thereby repressing the transcription of various genes including tumor suppressor gene. The over expression of HDAC8 was observed in many cancers and thus inhibition of this enzyme has emerged as an efficient cancer therapeutic strategy. In an effort to facilitate the future discovery of HDAC8 inhibitors, we developed two pharmacophore models containing six and five pharmacophoric features, respectively, using the representative structures from two molecular dynamic (MD) simulations performed in Gromacs 4.0.5 package. Various analyses of trajectories obtained from MD simulations have displayed the changes upon inhibitor binding. Thus utilization of the dynamically-responded protein structures in pharmacophore development has the added advantage of considering the conformational flexibility of protein. The MD trajectories were clustered based on single-linkage method and representative structures were taken to be used in the pharmacophore model development. Active site complimenting structure-based pharmacophore models were developed using Discovery Studio 2.5 program and validated using a dataset of known HDAC8 inhibitors. Virtual screening of chemical database coupled with drug-like filter has identified drug-like hit compounds that match the pharmacophore models. Molecular docking of these hits reduced the false positives and identified two potential compounds to be used in future HDAC8 inhibitor design.

Porcine testicular carbonyl reductase, PTCR which is one of the short chain dehydrogenases/reductases (SDR) superfamily catalyzes the NADPH-dependent reduction of carbonyl compounds including steroids and prostaglandins. Previously we reported C-terminal tail of PTCR was deleted due to a nonsynonymous single nucleotide variation (nsSNV). Here we identified from kinetic studies that the enzymatic properties for 5α-dihydrotestosterone (5α-DHT) were different between wild-type and C-terminal-deleted PTCRs. Compared to wild-type PTCR, C-terminal-deleted PTCR has much higher reduction rate. To investigate structural difference between wild-type and C-terminal-deleted PTCRs upon 5α-DHT binding, we performed molecular dynamics simulations for two complexes. Using trajectories, molecular interactions including hydrogen bonding patterns, distance between 5α-DHT and catalytic Tyr193, and interaction energies are analyzed and compared. During the MD simulation time, the dynamic behavior of C-terminal tail in wild-type PTCR is also examined using essential dynamics analysis. The results of our simulations reveal that the binding conformation of 5α-DHT in C-terminal-deleted PTCR is more favorable for reduction reaction in PTCR, which shows strong agreement with kinetic data. These structural findings provide valuable information to understand substrate specificity of PTCR and further kinetic properties of enzymes belonging to the SDR superfamily.

Acetylcholinesterase (AChE) catalyzes the hydrolysis of neurotransmitter acetylcholine to acetate and choline in a synaptic cleft. Deficits in cholinergic neurotransmitters are linked closely with the progression of Alzheimer's disease (AD), which is a neurodegenerative disorder characterized by memory impairment, and a disordered cognitive function. Since the previously approved AChE inhibitors, donepezil (Aricept), galantamine (Reminyl), and rivastigmine (Exelon), have side effects and several studies are being carried out out to develop novel AD drugs, we have applied a three-dimensional quantitative structure-activity relationship (3D QSAR) and structure-based pharmacophore modeling methodologies to identify potential candidate inhibitors against AChE. Herein, 3D QSAR and structure-based pharmacophore models were built from known inhibitors and crystal structures of human AChE in complex with donepezil, galantamine, huperzine A, and huprine W, respectively. The generated models were used as 3D queries to screen new scaffolds from various chemical databases. The hit compounds obtained from the virtual screening were subjected to an assessment of drug-like properties, followed by molecular docking. The final hit compounds were selected based on binding modes and molecular interactions in the active site of the enzyme. Furthermore, molecular dynamics simulations for AChE in complex with the final hits were performed to evaluate that they maintained stable interactions with the active site residues. The binding free energies of the final hits were also calculated using molecular mechanics/Poisson-Boltzmann surface area method. Taken together, we proposed that these hits can be promising candidates for anti-AD drugs.

Human chymase catalyzes the hydrolysis of peptide bonds. Three chymase inhibitors with very similar chemical structures but highly different inhibitory profiles towards the hydrolase function of chymase were selected with the aim of elucidating the origin of disparities in their biological activities. As a substrate (angiotensin-I) bound crystal structure is not available, molecular docking was performed to dock the substrate into the active site. Molecular dynamics simulations of chymase complexes with inhibitors and substrate were performed to calculate the binding orientation of inhibitors and substrate as well as to characterize conformational changes in the active site. The results elucidate details of the 3D chymase structure as well as the importance of K40 in hydrolase function. Binding mode analysis showed that substitution of a heavier Cl atom at the phenyl ring of most active inhibitor produced a great deal of variation in its orientation causing the phosphinate group to interact strongly with residue K40. Dynamics simulations revealed the conformational variation in region of V36-F41 upon substrate and inhibitor binding induced a shift in the location of K40 thus changing its interactions with them. Chymase complexes with the most active compound and substrate were used for development of a hybrid pharmacophore model which was applied in databases screening. Finally, hits which bound well at the active site, exhibited key interactions and favorable electronic properties were identified as possible inhibitors for chymase. This study not only elucidates inhibitory mechanism of chymase inhibitors but also provides key structural insights which will aid in the rational design of novel potent inhibitors of the enzyme. In general, the strategy applied in the current study could be a promising computational approach and may be generally applicable to drug design for other enzymes.

Human leukotriene A4 hydrolase (hLTA4H) is a bi-functional enzyme catalyzes the hydrolase and aminopeptidase functions upon the fatty acid and peptide substrates, respectively, utilizing the same but overlapping binding site. Particularly the hydrolase function of this enzyme catalyzes the rate-limiting step of the leukotriene (LT) cascade that converts the LTA4 to LTB4. This product is a potent pro-inflammatory activator of inflammatory responses and thus blocking this conversion provides a valuable means to design anti-inflammatory agents. Four structurally very similar chemical compounds with highly different inhibitory profile towards the hydrolase function of hLTA4H were selected from the literature. Molecular dynamics (MD) simulations of the complexes of hLTA4H with these inhibitors were performed and the results have provided valuable information explaining the reasons for the differences in their biological activities. Binding mode analysis revealed that the additional thiophene moiety of most active inhibitor helps the pyrrolidine moiety to interact the most important R563 and K565 residues. The hLTA4H complexes with the most active compound and substrate were utilized in the development of hybrid pharmacophore models. These developed pharmacophore models were used in screening chemical databases in order to identify lead candidates to design potent hLTA4H inhibitors. Final evaluation based on molecular docking and electronic parameters has identified three compounds of diverse chemical scaffolds as potential leads to be used in novel and potent hLTA4H inhibitor design.

Non-small cell lung cancer (NSCLC) is a lethal non-immunogenic malignancy and proto-oncogene ROS-1 tyrosine kinase is one of its clinically relevant oncogenic markers. The ROS-1 inhibitor, crizotinib, demonstrated resistance due to the Gly2032Arg mutation. To curtail this resistance, researchers developed lorlatinib against the mutated kinase. In the present study, a receptor-ligand pharmacophore model exploiting the key features of lorlatinib binding with ROS-1 was exploited to identify inhibitors against the wild-type (WT) and the mutant (MT) kinase domain. The developed model was utilized to virtually screen the TimTec flavonoids database and the retrieved drug-like hits were subjected for docking with the WT and MT ROS-1 kinase. A total of 10 flavonoids displayed higher docking scores than lorlatinib. Subsequent molecular dynamics simulations of the acquired flavonoids with WT and MT ROS-1 revealed no steric clashes with the Arg2032 (MT ROS-1). The binding free energy calculations computed via molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) demonstrated one flavonoid (Hit) with better energy than lorlatinib in binding with WT and MT ROS-1. The Hit compound was observed to bind in the ROS-1 selectivity pocket comprised of residues from the β-3 sheet and DFG-motif. The identified Hit from this investigation could act as a potent WT and MT ROS-1 inhibitor.

Cysteine-cysteine chemokine receptor 5 (CCR5) has been discovered as a co-receptor for cellular entry of human immunodeficiency virus (HIV). Moreover, the role of CCR5 in a variety of cancers and various inflammatory responses was also discovered. Despite the fact that several CCR5 antagonists have been investigated in clinical trials, only Maraviroc has been licensed for use in the treatment of HIV patients. This indicates that there is a need for novel CCR5 antagonists. Keeping this in mind, the present study was designed. The active CCR5 inhibitors with known IC50 value were selected from the literature and utilized to develop a ligand-based common feature pharmacophore model. The validated pharmacophore model was further used for virtual screening of drug-like databases obtained from the Asinex, Specs, InterBioScreen, and Eximed chemical libraries. Utilizing computational methods such as molecular docking studies, molecular dynamics simulations, and binding free energy calculation, the binding mechanism of selected inhibitors was established. The identified Hits not only showed better binding energy when compared to Maraviroc, but also formed stable interactions with the key residues and showed stable behavior throughout the 100 ns MD simulation. Our findings suggest that Hit1 and Hit2 may be potential candidates for CCR5 inhibition, and, therefore, can be considered for further CCR5 inhibition programs.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused over 1.4 million deaths worldwide. Repurposing existing drugs offers the fastest opportunity to identify new indications for existing drugs as a stable solution against coronavirus disease 2019 (COVID-19). The SARS-CoV-2 main protease (Mpro) is a critical target for designing potent antiviral agents against COVID-19. In this study, we identify potential inhibitors against COVID-19, using an amalgam of virtual screening, molecular dynamics (MD) simulations, and binding-free energy approaches from the Korea Chemical Bank drug repurposing (KCB-DR) database. The database screening of KCB-DR resulted in 149 binders. The dynamics of protein-drug complex formation for the seven top scoring drugs were investigated through MD simulations. Six drugs showed stable binding with active site of SARS-CoV-2 Mpro indicated by steady RMSD of protein backbone atoms and potential energy profiles. Furthermore, binding free energy calculations suggested the community-acquired bacterial pneumonia drug ceftaroline fosamil and the hepatitis C virus (HCV) protease inhibitor telaprevir are potent inhibitors against Mpro. Molecular dynamics and interaction analysis revealed that ceftaroline fosamil and telaprevir form hydrogen bonds with important active site residues such as Thr24, Thr25, His41, Thr45, Gly143, Ser144, Cys145, and Glu166 that is supported by crystallographic information of known inhibitors. Telaprevir has potential side effects, but its derivatives have good pharmacokinetic properties and are suggested to bind Mpro. We suggest the telaprevir derivatives and ceftaroline fosamil bind tightly with SARS-CoV-2 Mpro and should be validated through preclinical testing.

The rapid spread of COVID-19, caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a worldwide health emergency. Unfortunately, to date, a very small number of remedies have been to be found effective against SARS-CoV-2 infection. Therefore, further research is required to achieve a lasting solution against this deadly disease. Repurposing available drugs and evaluating natural product inhibitors against target proteins of SARS-CoV-2 could be an effective approach to accelerate drug discovery and development. With this strategy in mind, we derived Marine Natural Products (MNP)-based drug-like small molecules and evaluated them against three major target proteins of the SARS-CoV-2 virus replication cycle. A drug-like database from MNP library was generated using Lipinski's rule of five and ADMET descriptors. A total of 2,033 compounds were obtained and were subsequently subjected to molecular docking with 3CLpro, PLpro, and RdRp. The docking analyses revealed that a total of 14 compounds displayed better docking scores than the reference compounds and have significant molecular interactions with the active site residues of SARS-CoV-2 virus targeted proteins. Furthermore, the stability of docking-derived complexes was analyzed using molecular dynamics simulations and binding free energy calculations. The analyses revealed two hit compounds against each targeted protein displaying stable behavior, binding affinity, and molecular interactions. Our investigation identified two hit compounds against each targeted proteins displaying stable behavior, higher binding affinity and key residual molecular interactions, with good in silico pharmacokinetic properties, therefore can be considered for further in vitro studies.