Mycobacterium tuberculosis serine/threonine protein kinases: structural information for the design of their specific ATP-competitive inhibitors

Advanced modeling of salt-inducible kinase (SIK) inhibitors incorporating protein flexibility through molecular dynamics and cross-docking

Salt-inducible kinases (SIK1, SIK2, and SIK3) regulate metabolism and immune responses, making them promising targets for inflammatory and autoimmune diseases. Understanding inhibitor selectivity among isoforms is crucial for therapeutic development. In this study, 44 compounds were investigated as SIK inhibitors using molecular modeling. A flexible treatment of the kinases via molecular dynamics (MD) simulations captured binding site conformational changes, followed by molecular docking to generate protein kinase (PK)-ligand complex models. Ligand orientations were validated against crystallographic data using LigRMSD and interaction fingerprints (IFPs). A genetic algorithm was applied to select conformations that maximize correlation between docking energies and biological activities, yielding R² values of 0.821, 0.646, and 0.620 for SIK1, SIK2, and SIK3, respectively. Our results highlight the importance of protein flexibility in achieving accurate correlations between docking energies and experimental pIC50 values, enhancing inhibitor selectivity predictions.

Molecular docking, free energy calculations, ADMETox studies, DFT analysis, and dynamic simulations highlighting a chromene glycoside as a potential inhibitor of PknG in Mycobacterium tuberculosis

Introduction

Tuberculosis (TB), caused by the Mycobacterium tuberculosis (M.tb), remains a serious medical concern globally. Resistant M.tb strains are emerging, partly because M.tb can survive within alveolar macrophages, resulting in persistent infection. Protein kinase G (PknG) is a mycobacterial virulence factor that promotes the survival of M.tb in macrophages. Targeting PknG could offer an opportunity to suppress the resistant M.tb strains.

Methods

In the present study, multiple computational tools were adopted to screen a library of 460,000 molecules for potential inhibitors of PknG of M.tb.

Results and discussions

Seven Hits (1-7) were identified with binding affinities exceeding that of the reference compound (AX20017) towards the PknG catalytic domain. Next, the ADMETox studies were performed to identify the best hit with appropriate drug-like properties. The chromene glycoside (Hit 1) was identified as a potential PknG inhibitor with better pharmacokinetic and toxicity profiles rendering it a potential drug candidate. Furthermore, quantum computational analysis was conducted to assess the mechanical and electronic properties of Hit 1, providing guidance for further studies. Molecular dynamics simulations were also performed for Hit 1 against PknG, confirming the stability of its complex. In sum, the findings in the current study highlight Hit 1 as a lead with potential for development of drugs capable of treating resistant TB.

Rosmarinic acid turned α-syn oligomers into non-toxic species preserving microtubules in Raw 264.7 cells

Parkinson's disease (PD) is the second most prevalent neurodegenerative disorder worldwide, and the therapeutic is focused on several approaches including the inhibition of fibril formation by small compounds, avoiding the formation of cytotoxic oligomers. Thus, we decided to explore the capacity of compounds carrying catechol moieties to inhibit the progression of α-synuclein. Overall, the compounds rosmarinic acid (1), carnosic acid (2), carnosol (3), epiisorosmanol (4), and rosmanol (5) avoid the progression of fibril formation assessed by Thiofavine T (ThT), and atomic force microscopy images showed that morphology is influenced for the actions of compounds over fibrillization. Moreover, ITC experiments showed a Kd varying from 28 to 51 µM, the ΔG showed that the reaction between compounds and α-syn is spontaneous, and ΔH is associated with an exothermic reaction, suggesting the interactions of hydrogen bonds among compounds and α-syn. Docking experiments reinforce this idea showing the intermolecular interactions are mostly hydrogen bonding within the sites 2, 9, and 3/13 of α-synuclein, and compounds 1 and 5. Thus, compound 1, rosmarinic acid, interestingly interacts better with site 9 through catechol and Lysines. In cultured Raw 264. 7 cells, the presence of compounds showed that most of them can promote cell differentiation, especially rosmarinic acid, and rosmanol, both preserving tubulin cytoskeleton. However, once we evaluated whether or not the aggregates pre-treated with compounds could prevent the disruption of microtubules of Raw 264.7 cells, only pre-treated aggregates with rosmarinic acid prevented the disruption of the cytoskeleton. Altogether, we showed that especially rosmarinic acid not only inhibits α-syn but stabilizes the remaining aggregates turning them into not-toxic to Raw 264.7 cells suggesting a main role in cell survival and antigen processing in response to external α-syn aggregates.

Virtual Screening Identifies Novel and Potent Inhibitors of Mycobacterium tuberculosis PknB with Antibacterial Activity

Mycobacterium tuberculosis protein kinase B (PknB) is essential to mycobacterial growth and has received considerable attention as an attractive target for novel anti-tuberculosis drug development. Here, virtual screening, validated by biological assays, was applied to select candidate inhibitors of M. tuberculosis PknB from the Specs compound library (www.specs.net). Fifteen compounds were identified as hits and selected for in vitro biological assays, of which three indoles (2, AE-848/42799159; 4, AH-262/34335013; 10, AP-124/40904362) inhibited growth of M. tuberculosis H37Rv with minimal inhibitory concentrations of 6.2, 12.5, and 6.2 μg/mL, respectively. Two compounds, 2 and 10, inhibited M. tuberculosis PknB activity in vitro, with IC₅₀ values of 14.4 and 12.1 μM, respectively, suggesting this to be the likely basis of their anti-tubercular activity. In contrast, compound 4 displayed anti-tuberculosis activity against M. tuberculosis H37Rv but showed no inhibition of PknB activity (IC₅₀ > 128 μM). We hypothesize that hydrolysis of its ethyl ester to a carboxylate moiety generates an active species that inhibits other M. tuberculosis enzymes. Molecular dynamics simulations of modeled complexes of compounds 2, 4, and 10 bound to M. tuberculosis PknB indicated that compound 4 has a lower affinity for M. tuberculosis PknB than compounds 2 and 10, as evidenced by higher calculated binding free energies, consistent with experiment. Compounds 2 and 10 therefore represent candidate inhibitors of M. tuberculosis PknB that provide attractive starting templates for optimization as anti-tubercular agents.

Flavanones: A potential natural inhibitor of the ATP binding site of PknG of Mycobacterium tuberculosis

Over the years, Mycobacterium tuberculosis has been one of the major causes of death worldwide. As several clinical isolates of the bacteria have developed drug resistance against the target sites of the current therapeutic agents, the development of a novel drug is the pressing priority. According to recent studies on Mycobacterium tuberculosis, ATP binding sites of Mycobacterium tuberculosis serine/threonine protein kinases (MTPKs) have been identified as the new promising drug target. Among the several other protein kinases (PKs), Protein kinase G (PknG) was selected for the study because of its crucial role in modulating bacterium's metabolism to survive in host macrophages. In this work, we have focused on the H37Rv strain of Mycobacterium tuberculosis. A list of 477 flavanones obtained from the PubChem database was docked one by one against the crystallized and refined structure of PknG by in-silico techniques. Initially, potential inhibitors were narrowed down by preliminary docking. Flavanones were then selected using binding energies ranging from -7.9 kcal.mol^-1 to -10.8 kcal.mol^-1. This was followed by drug-likeness prediction, redocking analysis, and molecular dynamics simulations. Here, we have used experimentally confirmed drug AX20017 as a reference to determine candidate compounds that can act as potential inhibitors for PknG. PubChem165506, PubChem242065, PubChem688859, PubChem101367767, PubChem3534982, and PubChem42607933 were identified as possible target site inhibitors for PknG with a desirable negative binding energy of -8.1, -8.3, -8.4, -8.8, -8.6 and -7.9 kcal.mol^-1 respectively. Communicated by Ramaswamy H. Sarma.

Identifying RO9021 as a Potential Inhibitor of PknG from Mycobacterium tuberculosis: Combinative Computational and In Vitro Studies

Tuberculosis (TB) is an infectious disease caused by Mycobacterium tuberculosis (Mtb). Despite being considered curable and preventable, the increase of antibiotic resistance is becoming a serious public health problem. Mtb is a pathogen capable of surviving in macrophages, causing long-term latent infection where the mycobacterial serine/threonine protein kinase G (PknG) plays a protective role. Therefore, PknG is an important inhibitory target to prevent Mtb from entering the latency stage. In this study, we use a pharmacophore-based virtual screening and biochemical assays to identify the compound RO9021 (CHEMBL3237561) as a PknG inhibitor. In detail, 1.5 million molecules were screened using a scalable cloud-based setup, identifying 689 candidates, which were further subjected to additional screening employing molecular docking. Molecular docking spotted 62 compounds with estimated binding affinities of -7.54 kcal/mol (s.d. = 0.77 kcal/mol). Finally, 14 compounds were selected for in vitro experiments considering previously reported biological activities and commercial availability. In vitro assays of PknG activity showed that RO9021 inhibits the kinase activity similarly to AX20017, a known inhibitor. The inhibitory effect was found to be dose dependent with a relative IC₅₀ value of 4.4 ± 1.1 μM. Molecular dynamics simulations predicted that the PknG-RO9021 complex is stable along the tested timescale. Altogether, our study indicates that RO9021 is a noteworthy drug candidate for further developing new anti-TB drugs that hold excellent reported pharmacokinetic parameters.

Specificity and Reactivity of Mycobacterium tuberculosis Serine/Threonine Kinases PknG and PknB

Mycobacterium tuberculosis (Mtb), the causative agent of Tuberculosis, has 11 eukaryotic-like serine/threonine protein kinases, which play essential roles in cell growth, signal transduction, and pathogenesis. Protein kinase G (PknG) regulates the carbon and nitrogen metabolism by phosphorylation of the glycogen accumulation regulator (GarA) protein at Thr21. Protein kinase B (PknB) is involved in cell wall synthesis and cell shape, as well as phosphorylates GarA but at Thr22. While PknG seems to be constitutively activated and recognition of GarA requires phosphorylation in its unstructured tail, PknB activation is triggered by phosphorylation of its activation loop, which allows binding of the forkhead-associated domain of GarA. In the present work, we used molecular dynamics and quantum-mechanics/molecular mechanics simulations of the catalytically competent complex and kinase activity assays to understand PknG/PknB specificity and reactivity toward GarA. Two hydrophobic residues in GarA, Val24 and Phe25, seem essential for PknG binding and allow specificity for Thr21 phosphorylation. On the other hand, phosphorylated residues in PknB bind Arg26 in GarA and regulate its specificity for Thr22. We also provide a detailed analysis of the free energy profile for the phospho-transfer reaction and show why PknG has a constitutively active conformation not requiring priming phosphorylation in contrast to PknB. Our results provide new insights into these two key enzymes relevant for Mtb and the mechanisms of serine/threonine phosphorylation in bacteria.

Molecular Modeling Studies of N-phenylpyrimidine-4-amine Derivatives for Inhibiting FMS-like Tyrosine Kinase-3

Overexpression and frequent mutations in FMS-like tyrosine kinase-3 (FLT3) are considered risk factors for severe acute myeloid leukemia (AML). Hyperactive FLT3 induces premature activation of multiple intracellular signaling pathways, resulting in cell proliferation and anti-apoptosis. We conducted the computational modeling studies of 40 pyrimidine-4,6-diamine-based compounds by integrating docking, molecular dynamics, and three-dimensional structure-activity relationship (3D-QSAR). Molecular docking showed that K644, C694, F691, E692, N701, D829, and F830 are critical residues for the binding of ligands at the hydrophobic active site. Molecular dynamics (MD), together with Molecular Mechanics Poison-Boltzmann/Generalized Born Surface Area, i.e., MM-PB(GB)SA, and linear interaction energy (LIE) estimation, provided critical information on the stability and binding affinity of the selected docked compounds. The MD study suggested that the mutation in the gatekeeper residue F691 exhibited a lower binding affinity to the ligand. Although, the mutation in D835 in the activation loop did not exhibit any significant change in the binding energy to the most active compound. We developed the ligand-based comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) models. CoMFA (q² = 0.802, r² = 0.983, and QF32 = 0.698) and CoMSIA (q² = 0.725, r² = 0.965 and QF32 = 0.668) established the structure-activity relationship (SAR) and showed a reasonable external predictive power. The contour maps from the CoMFA and CoMSIA models could explain valuable information about the favorable and unfavorable positions for chemical group substitution, which can increase or decrease the inhibitory activity of the compounds. In addition, we designed 30 novel compounds, and their predicted pIC₅₀ values were assessed with the CoMSIA model, followed by the assessment of their physicochemical properties, bioavailability, and free energy calculation. The overall outcome could provide valuable information for designing and synthesizing more potent FLT3 inhibitors.

Natural products from Brazilian biodiversity identified as potential inhibitors of PknA and PknB of M. tuberculosis using molecular modeling tools

Mycobacterium tuberculosis was discovered in 1882 by Robert Koch but, since its discovery, the tuberculosis (TB) epidemic has endured, being one of the top 10 causes of death worldwide. Drug-resistant TB continues to be a public health threat and bioactive compounds with a new mode of action (MoA) are needed to overcome this. Since natural products are described as important sources for the development of new drugs, the objective of this work was to identify potential ligands from Brazilian natural products (NPs) for M. tuberculosis targets using molecular modeling tools. Using chemogenomics we identified the Serine/Threonine Protein Kinase PknB as a putative target for 13 NPs from a database from Brazilian biodiversity (NuBBE). Literature data supported further investigation of NuBBE105, NuBBE598, NuBBE936, NuBBE964, NuBBE1045, and NuBBE1180 by molecular docking and dynamics. Key interactions were observed with PknB and simulations confirmed stability and favorable binding energies. Considering structural similarity with PknB, we further explored binding of the NPs to PknA, critical for M. tuberculosis survival, and all of them resembled important interactions with the enzyme, showing stable and favorable binding energies, whilst van der Waals interactions seem to play a key role for binding to PknA and PknB. NuBBE936 and NuBBE1180 have already had their antimycobacterial activity reported and our results can provide a basis for their MoA. Finally, the other NPs which have not been tested against M. tuberculosis deserve further investigation, aiming at the discovery of antimycobacterial drug candidates with innovative MoA.

Computational Modeling to Explain Why 5,5-Diarylpentadienamides are TRPV1 Antagonists

Several years ago, the crystallographic structures of the transient receptor potential vanilloid 1 (TRPV1) in the presence of agonists and antagonists were reported, providing structural information about its chemical activation and inactivation. TRPV1’s activation increases the transport of calcium and sodium ions, leading to the excitation of sensory neurons and the perception of pain. On the other hand, its antagonistic inactivation has been explored to design analgesic drugs. The interactions between the antagonists 5,5-diarylpentadienamides (DPDAs) and TRPV1 were studied here to explain why they inactivate TRPV1. The present work identified the structural features of TRPV1–DPDA complexes, starting with a consideration of the orientations of the ligands inside the TRPV1 binding site by using molecular docking. After this, a chemometrics analysis was performed (i) to compare the orientations of the antagonists (by using LigRMSD), (ii) to describe the recurrent interactions between the protein residues and ligand groups in the complexes (by using interaction fingerprints), and (iii) to describe the relationship between topological features of the ligands and their differential antagonistic activities (by using a quantitative structure–activity relationship (QSAR) with 2D autocorrelation descriptors). The interactions between the DPDA groups and the residues Y511, S512, T550, R557, and E570 (with a recognized role in the binding of classic ligands), and the occupancy of isoquinoline or 3-hydroxy-3,4-dihydroquinolin-2(1H)-one groups of the DPDAs in the vanilloid pocket of TRPV1 were clearly described. Based on the results, the structural features that explain why DPDAs inactivate TRPV1 were clearly exposed. These features can be considered for the design of novel TRPV1 antagonists.

Sequential Dihedral Angles (SDAs): A Method for Evaluating the 3D Structure of Proteins

One of the most important steps in modeling three-dimensional (3D) structures of proteins is the evaluation of the constructed models. The present study suggests that the correctness of a structure may be tested by using the characteristics of sequential dihedral angles (SDAs) between adjacent alpha-carbons (C_α) in the main chains of proteins. From our studies on protein structures in the protein data bank (PDB), the SDAs between the C_α in the main chains are limited in their values. In addition, the sum of the absolute values of the three sequential dihedral angles (SDAs) can never be 0 degree. Moreover, 48 degrees is the lowest value existing for the sum of the absolute values of three sequential dihedral angles (SDAs). Thus, the SDAs between the alpha-carbons along the main chains of proteins may be a useful parameter for evaluating anomalies in protein structures.

The latest automated docking technologies for novel drug discovery

INTRODUCTION

Molecular docking has been consolidated as one of the most important methods in the molecular modeling field. It has been recognized as a prominent tool in the study of protein-ligand complexes, to describe intermolecular interactions, to accurately predict poses of multiple ligands, to discover novel promising bioactive compounds. Molecular docking methods have evolved in terms of their accuracy and reliability; but there are pending issues to solve for improving the connection between the docking results and the experimental evidence.

AREAS COVERED

In this article, the author reviews very recent innovative molecular docking applications with special emphasis on reverse docking, treatment of protein flexibility, the use of experimental data to guide the selection of docking poses, the application of Quantum mechanics(QM) in docking, and covalent docking.

EXPERT OPINION

There are several issues being worked on in recent years that will lead to important breakthroughs in molecular docking methods in the near future These developments are related to more efficient exploration of large datasets and receptor conformations, advances in electronic description, and the use of structural information for guiding the selection of results.

Considerations for Docking of Selective Angiotensin-Converting Enzyme Inhibitors

The angiotensin-converting enzyme (ACE) is a two-domain dipeptidylcarboxypeptidase, which has a direct involvement in the control of blood pressure by performing the hydrolysis of angiotensin I to produce angiotensin II. At the same time, ACE hydrolyzes other substrates such as the vasodilator peptide bradykinin and the anti-inflammatory peptide N-acetyl-SDKP. In this sense, ACE inhibitors are bioactive substances with potential use as medicinal products for treatment or prevention of hypertension, heart failures, myocardial infarction, and other important diseases. This review examined the most recent literature reporting ACE inhibitors with the help of molecular modeling. The examples exposed here demonstrate that molecular modeling methods, including docking, molecular dynamics (MD) simulations, quantitative structure-activity relationship (QSAR), etc, are essential for a complete structural picture of the mode of action of ACE inhibitors, where molecular docking has a key role. Examples show that too many works identified ACE inhibitory activities of natural peptides and peptides obtained from hydrolysates. In addition, other works report non-peptide compounds extracted from natural sources and synthetic compounds. In all these cases, molecular docking was used to provide explanation of the chemical interactions between inhibitors and the ACE binding sites. For docking applications, most of the examples exposed here do not consider that: (i) ACE has two domains (nACE and cACE) with available X-ray structures, which are relevant for the design of selective inhibitors, and (ii) nACE and cACE binding sites have large dimensions, which leads to non-reliable solutions during docking calculations. In support of the solution of these problems, the structural information found in Protein Data Bank (PDB) was used to perform an interaction fingerprints (IFPs) analysis applied on both nACE and cACE domains. This analysis provides plots that identify the chemical interactions between ligands and both ACE binding sites, which can be used to guide docking experiments in the search of selective natural components or novel drugs. In addition, the use of hydrogen bond constraints in the S₂ and S₂′ subsites of nACE and cACE are suggested to guarantee that docking solutions are reliable.

Molecular modeling of a series of dehydroquinate dehydratase type II inhibitors of Mycobacterium tuberculosis and design of new binders

Tuberculosis, caused by Mycobacterium tuberculosis (M. tuberculosis), is still responsible for a large number of fatal cases, especially in developing countries with alarming rates of incidence and prevalence worldwide. Mycobacterium tuberculosis has a remarkable ability to develop new resistance mechanisms to the conventional antimicrobials treatment. Because of this, there is an urgent need for novel bioactive compounds for its treatment. The dehydroquinate dehydratase II (DHQase II) is considered a key enzyme of shikimate pathway, and it can be used as a promising target for the design of new bioactive compounds with antibacterial action. The aim of this work was the construction of QSAR models to aid the design of new potential DHQase II inhibitors. For that purpose, various molecular modeling approaches, such as activity cliff, QSAR models and computer-aided ligand design were utilized. A predictive in silico 4D-QSAR model was built using a database comprising 86 inhibitors of DHQase II, and the model was used to predict the activity of the designed ligands. The obtained model proved to predict well the DHQase II inhibition for an external validation dataset ([Formula: see text] = 0.72). Also, the Activity Cliff analysis shed light on important structural features applied to the ligand design.

Structural Requirements of N-alpha-Mercaptoacetyl Dipeptide (NAMdP) Inhibitors of Pseudomonas Aeruginosa Virulence Factor LasB: 3D-QSAR, Molecular Docking, and Interaction Fingerprint Studies

The zinc metallopeptidase Pseudomonas elastase (LasB) is a virulence factor of Pseudomonas aeruginosa (P. aeruginosa), a pathogenic bacterium that can cause nosocomial infections. The present study relates the structural analysis of 118 N-alpha-mercaptoacetyl dipeptides (NAMdPs) as LasB inhibitors. Field-based 3D-QSAR and molecular docking methods were employed to describe the essential interactions between NAMdPs and LasB binding sites, and the chemical features that determine their differential activities. We report a predictive 3D-QSAR model that was developed according to the internal and external validation tests. The best model, including steric, electrostatic, hydrogen bond donor, hydrogen bond acceptor, and hydrophobic fields, was found to depict a three-dimensional map with the local positive and negative effects of these chemotypes on the LasB inhibitory activities. Furthermore, molecular docking experiments yielded bioactive conformations of NAMdPs inside the LasB binding site. The series of NAMdPs adopted a similar orientation with respect to phosphoramidon within the LasB binding site (crystallographic reference), where the backbone atoms of NAMdPs are hydrogen-bonded to the LasB residues N112, A113, and R198, similarly to phosphoramidon. Our study also included a deep description of the residues involved in the protein-ligand interaction patterns for the whole set of NAMdPs, through the use of interaction fingerprints (IFPs).

Predictive Binding Affinity of Plant-Derived Natural Products Towards the Protein Kinase G Enzyme of Mycobacterium tuberculosis (MtPknG)

Tuberculosis (TB), caused by Mycobacterium tuberculosis, is a growing public health concern worldwide, especially with the emerging challenge of drug resistance to the current drugs. Efforts to discover and develop novel, more effective, and safer anti-TB drugs are urgently needed. Products from natural sources, such as medicinal plants, have played an important role in traditional medicine and continue to provide some inspiring templates for the design of new drugs. Protein kinase G, produced by M. tuberculosis (MtPKnG), is a serine/threonine kinase, that has been reported to prevent phagosome-lysosome fusion and help prolong M. tuberculosis survival within the host’s macrophages. Here, we used an in silico, target-based approach (docking) to predict the interactions between MtPknG and 84 chemical constituents from two medicinal plants (Pelargonium reniforme and Pelargonium sidoides) that have a well-documented historical use as natural remedies for TB. Docking scores for ligands towards the target protein were calculated using AutoDock Vina as the predicted binding free energies. Ten flavonoids present in the aerial parts of P. reniforme and/or P.sidoides showed docking scores ranging from −11.1 to −13.2 kcal/mol. Upon calculation of all ligand efficiency indices, we observed that the (−ΔG/MW) ligand efficiency index for flavonoids (4), (5) and (7) was similar to the one obtained for the AX20017 control. When taking all compounds into account, we observed that the best (−ΔG/MW) efficiency index was obtained for coumaric acid, coumaraldehyde, p-hydroxyphenyl acetic acid and p-hydroxybenzyl alcohol. We found that methyl gallate and myricetin had ligand efficiency indices superior and equal to the AX20017 control efficiency, respectively. It remains to be seen if any of the compounds screened in this study exert an effect in M. tuberculosis-infected macrophages.

Docking, Interaction Fingerprint, and Three-Dimensional Quantitative Structure-Activity Relationship (3D-QSAR) of Sigma1 Receptor Ligands, Analogs of the Neuroprotective Agent RC-33

The human Sigma1 receptor (S1R), which has been identified as a target with an important role in neuropsychological disorders, was first crystallized 3 years ago. Since S1R structure has no relation with another previous crystallized structures, the presence of the new crystal is an important hallmark for the design of agonists and antagonists against this important target. Some years ago, our group identified RC-33, a potent and selective S1R agonist, endowed with neuroprotective properties. In this work, drawing on new structural information, we studied the interactions of RC-33 and its analogs with the S1R binding site by using computational methods such as docking, interaction fingerprints, and receptor-guided alignment three dimensional quantitative structure–activity relationship (3D-QSAR). We found that RC-33 and its analogs adopted similar orientations within S1R binding site, with high similitude with orientations of the crystallized ligands; such information was used for identifying the residues involved in chemical interactions with ligands. Furthermore, the structure-activity relationship of the studied ligands was adequately described considering classical QSAR tests. All relevant aspects of the interactions between the studied compounds and S1R were covered here, through descriptions of orientations, binding interactions, and features that influence differential affinities. In this sense, the present results could be useful in the future design of novel S1R modulators.

Rational Design of Novel Phosphoinositide 3-Kinase Gamma (PI3Kγ) Selective Inhibitors: A Computational Investigation Integrating 3D-QSAR, Molecular Docking and Molecular Dynamics Simulation

Phosphoinositide 3-kinase gamma (PI3Kγ) draws an increasing attention due to its link with deadly cancer, chronic inflammation and allergy. But the development of PI3Kγ selective inhibitors is still a challenging endeavor because of the high sequence homology with the other PI3K isoforms. In order to acquire valuable information about the interaction mechanism between potent inhibitors and PI3Kγ, a series of PI3Kγ isoform-selective inhibitors were analyzed by a systematic computational method, combining 3D-QSAR, molecular docking, molecular dynamic (MD) simulations, free energy calculations and decomposition. The general structure-activity relationships were revealed and some key residues relating to selectivity and high activity were highlighted. It provides precious guidance for rational virtual screening, modification and design of selective PI3Kγ inhibitors. Finally, ten novel inhibitors were optimized and P10 showed satisfactory predicted bioactivity, demonstrating the feasibility to develop potent PI3Kγ inhibitors through this computational modeling and optimization.

Substituent effects on the stability, physicochemical properties and chemical reactivity of nitroimidazole derivatives with potential antiparasitic effect: a computational study

Neglected tropical diseases caused by parasitic pathogens have caused an increase in research interest in drug discovery.