Molecular docking and dynamics simulation, receptor-based hypothesis: application to identify novel sirtuin 2 inhibitors

Virtual Screening in the Identification of Sirtuins’ Activity Modulators

Sirtuins are NAD⁺-dependent deac(et)ylases with different subcellular localization. The sirtuins’ family is composed of seven members, named SIRT-1 to SIRT-7. Their substrates include histones and also an increasing number of different proteins. Sirtuins regulate a wide range of different processes, ranging from transcription to metabolism to genome stability. Thus, their dysregulation has been related to the pathogenesis of different diseases. In this review, we discussed the pharmacological approaches based on sirtuins’ modulators (both inhibitors and activators) that have been attempted in in vitro and/or in in vivo experimental settings, to highlight the therapeutic potential of targeting one/more specific sirtuin isoform(s) in cancer, neurodegenerative disorders and type 2 diabetes. Extensive research has already been performed to identify SIRT-1 and -2 modulators, while compounds targeting the other sirtuins have been less studied so far. Beside sections dedicated to each sirtuin, in the present review we also included sections dedicated to pan-sirtuins’ and to parasitic sirtuins’ modulators. A special focus is dedicated to the sirtuins’ modulators identified by the use of virtual screening.

Repurposing of FDA-approved drugs as potential inhibitors of the SARS-CoV-2 main protease: Molecular insights into improved therapeutic discovery

With numerous infections and fatalities, COVID-19 has wreaked havoc around the globe. The main protease (Mpro), which cleaves the polyprotein to form non-structural proteins, thereby helping in the replication of SARS-CoV-2, appears as an attractive target for antiviral therapeutics. As FDA-approved drugs have shown effectiveness in targeting Mpro in previous SARS-CoV(s), molecular docking and virtual screening of existing antiviral, antimalarial, and protease inhibitor drugs were carried out against SARS-CoV-2 Mpro. Among 53 shortlisted drugs with binding energies lower than that of the crystal-bound inhibitor α-ketoamide 13 b (-6.7 kcal/mol), velpatasvir, glecaprevir, grazoprevir, baloxavir marboxil, danoprevir, nelfinavir, and indinavir (-9.1 to -7.5 kcal/mol) were the most significant on the list (hereafter referred to as the 53-list). Molecular dynamics (MD) simulations confirmed the stability of their Mpro complexes, with the MMPBSA binding free energy (ΔGbind) ranging between -124 kJ/mol (glecaprevir) and -28.2 kJ/mol (velpatasvir). Despite having the lowest initial binding energy, velpatasvir exhibited the highest ΔGbind value for escaping the catalytic site during the MD simulations, indicating its reduced efficacy, as observed experimentally. Available inhibition assay data adequately substantiated the computational forecast. Glecaprevir and nelfinavir (ΔGbind = -95.4 kJ/mol) appear to be the most effective antiviral drugs against Mpro. Furthermore, the remaining FDA drugs on the 53-list can be worth considering, since some have already demonstrated antiviral activity against SARS-CoV-2. Hence, theoretical pKi (Ki = inhibitor constant) values for all 53 drugs were provided. Notably, ΔGbind directly correlates with the average distance of the drugs from the His41-Cys145 catalytic dyad of Mpro, providing a roadmap for rapid screening and improving the inhibitor design against SARS-CoV-2 Mpro.

Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development

Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.

Experimental and theoretical investigations on the interaction of l-methionine molecules with α-chymotrypsin in the aqueous solution using various methods

l-Methionine (l-Met) is one of the necessary amino acids that play unparalleled roles, influencing both the protein structure and metabolism. Understanding the interactions between proteins and small molecules can be realized by various perspectives, and this is significant for the progression of basic sciences and drug development. In this study, the variations in the stability, function, and structure of α-Chymotrypsin (α-Chy) in the presence of l-Met were investigated using spectroscopic and computational approaches. The results of the UV-vis absorption demonstrated that α-Chy had a maximum peak at 280 nm due to the Trp residue. Hyperchromism shift was seen in the presence of l-Met. Ground state system was formed in the presence of l-Met, as confirmed by the fluorescence studies. Fluorescence variations also revealed static quenching. The CD spectra also represented the alteration of the enzyme with an increase in the α-helix and a decrease in the β-sheet. The activity of α-Chy was incremented in the presence of l-Met. Therefore, l-Met served as an activator. Molecular docking results also indicated a negative amount for the Gibbs free energy of the binding of l-Met to the enzyme. α-Chy became more stable in the presence of l-Met, based on the molecular dynamics simulation.

Synthesis and Evaluation of Novel Benzofuran Derivatives as Selective SIRT2 Inhibitors

A series of benzofuran derivatives were designed and synthesized, and their inhibitory activites were measured against the SIRT1–3. The enzymatic assay showed that all the compounds showed certain anti-SIRT2 activity and selective over SIRT1 and SIRT3 with IC₅₀ (half maximal inhibitory concentration) values at the micromolar level. The preliminary structure–activity relationships were analyzed and the binding features of compound 7e (IC₅₀ 3.81 µM) was predicted using the CDOCKER program. The results of this research could provide informative guidance for further optimizing benzofuran derivatives as potent SIRT2 inhibitors.

Computational Studies on Sirtuins from Trypanosoma cruzi: Structures, Conformations and Interactions with Phytochemicals

Background

The silent-information regulator 2 proteins, otherwise called sirtuins, are currently considered as emerging anti-parasitic targets. Nicotinamide, a pan-sirtuin inhibitor, is known to cause kinetoplast alterations and the arrested growth of T. cruzi, the protozoan responsible for Chagas disease. These observations suggested that sirtuins from this parasite (TcSir2rp1 and TcSir2rp3) could play an important role in the regulation of the parasitic cell cycle. Thus, their inhibition could be exploited for the development of novel anti-trypanosomal compounds.

Methods

Homology modeling was used to determine the three-dimensional features of the sirtuin TcSir2rp1 from T. cruzi. The apo-form of human SIRT2 and the same structure solved in complex with its co-substrate NAD⁺ allowed the modeling of TcSir2rp1 in the open and closed conformational states. Molecular docking studies were then carried out. A library composed of fifty natural and diverse compounds that are known to be active against this parasite, was established based on the literature and virtually screened against TcSir2rp1 and TcSir2rp3, which was previously modeled by our group.

Results

In this study, two conformational states of TcSir2rp1 were described for the first time. The molecular docking results of compounds capable of binding sirtuins proved to be meaningful when the closed conformation of the protein was taken into account for calculations. This specific conformation was then used for the virtual screening of antritrypanosomal phytochemicals against TcSir2rp1 and TcSir2rp3. The calculations identified a limited number of scaffolds extracted from Vismia orientalis, Cussonia zimmermannii, Amomum aculeatum and Anacardium occidentale that potentially interact with both proteins.

Conclusions

The study provided reliable models for future structure-based drug design projects concerning sirtuins from T. cruzi. Molecular docking studies highlighted not only the advantages of performing in silico interaction studies on their closed conformations but they also suggested the potential mechanism of action of four phytochemicals known for their anti-trypanosomal activity in vitro.

T. cruzi is a protozoan pathogen responsible for Chagas disease. Current therapies rely only on a very small number of drugs, most of which are inadequate because of their severe host toxicity or because of their susceptibility to drug-resistance mechanisms. To determine efficient therapeutic alternatives, the identification of new biotargets and detailed knowledge of their structures are essential. Sirtuins from T. cruzi have been recently considered as promising targets for the development of new treatments for Chagas disease. Inhibition of their activity has been shown to significantly interfere with the life cycle of the parasite. T. cruzi possesses genes encoding two sirtuin-like proteins, TcSIR2rp1 and TcSIR2rp3. The structures of these enzymes were theoretically elucidated in this work, which also focused on the impact of their possible conformational states on computational interaction studies. A small library of phytochemicals that are active against the parasite was built and screened against the most meaningful conformations, identifying a restricted number of scaffolds that potentially interact with the modeled proteins. For these hits, a mechanism of action related to interactions with sirtuins was proposed.

Identification of inhibitor binding site in human sirtuin 2 using molecular docking and dynamics simulations

The ability to identify the site of a protein that can bind with high affinity to small, drug-like compounds has been an important goal in drug design. Sirtuin 2 (SIRT2), histone deacetylase protein family, plays a central role in the regulation of various pathways. Hence, identification of drug for SIRT2 has attracted great interest in the drug discovery community. To elucidate the molecular basis of the small molecules interactions to inhibit the SIRT2 function we employed the molecular docking, molecular dynamics simulations, and the molecular mechanism Poisson-Boltzmann/surface area (MM-PBSA) calculations. Five well know inhibitors such as suramin, mol-6, sirtinol, 67, and nf675 were selected to establish the nature of the binding mode of the inhibitors in the SIRT2 active site. The molecular docking and dynamics simulations results revealed that the hydrogen bonds between Arg97 and Gln167 are crucial to inhibit the function of SIRT2. In addition, the MM-PBSA calculations revealed that binding of inhibitors to SIRT2 is mainly driven by van der Waals/non-polar interactions. Although the five inhibitors are very different in structure, shape, and electrostatic potential, they are able to fit in the same binding pocket. These findings from this study provide insights to elucidate the binding pattern of SIRT2 inhibitors and help in the rational structure-based design of novel SIRT2 inhibitors with improved potency and better resistance profile.

Insight the C-site pocket conformational changes responsible for sirtuin 2 activity using molecular dynamics simulations

Sirtuin belongs to a family of typical histone deacetylase which regulates the fundamental cellular biological processes including gene expression, genome stability, mitosis, nutrient metabolism, aging, mitochondrial function, and cell motility. Michael et. al. reported that B-site mutation (Q167A and H187A) decreased the SIRT2 activity but still the structural changes were not reported. Hence, we performed 5 ns molecular dynamics (MD) simulation on SIRT2 Apo-form and complexes with substrate/NAD(+) and inhibitor of wild type (WT), Q167A, and H187A. The results revealed that the assembly and disassembly of C-site induced by presence of substrate/NAD(+) and inhibitor, respectively. This assembly and disassembly was mainly due to the interaction between the substrate/NAD(+) and inhibitor and F96 and the distance between F96 and H187 which are present at the neck of the C-site. MD simulations suggest that the conformational change of L3 plays a major role in assembly and disassembly of C-site. Our current results strongly suggest that the distinct conformational change of L3 as well as the assembly and disassembly of C-site plays an important role in SIRT2 deacetylation function. Our study unveiled the structural changes of SIRT2 in presence of NAD(+) and inhibitor which should be helpful to improve the inhibitory potency of SIRT2.

Investigation of silent information regulator 1 (Sirt1) agonists from Traditional Chinese Medicine

Silent information regulator 1 (Sirt1), a class III nicotinamide adenine dinucleotide dependent histone deacetylases, is important in cardioprotection, neuroprotection, metabolic disease, calorie restriction, and diseases associated with aging. Traditional Chinese Medicine (TCM) compounds from TCM Database@Taiwan ( http://tcm.cmu.edu.tw/ ) were employed for screening potent Sirt1 agonists, and molecular dynamics (MD) simulation was implemented to simulate ligand optimum docking poses and protein structure under dynamic conditions. TCM compounds such as (S)-tryptophan-betaxanthin, 5-O-feruloylquinic acid, and RosA exhibited good binding affinity across different computational methods, and their drug-like potential were validated by MD simulation. Docking poses indicate that the carboxylic group of the three candidates generated H-bonds with residues in the protein chain from Ser441 to Lys444 and formed H-bond, π-cation interactions, or hydrophobic contacts with Phe297 and key active residue, His363. During MD, stable π-cation interactions with residues Phe273 or Arg274 were formed by (S)-tryptophan-betaxanthin and RosA. All candidates were anchored to His363 by stable π- or H-bonds. Hence, we propose (S)-tryptophan-betaxanthin, 5-O-feruloylquinic acid, and RosA as potential lead compounds that can be further tested in drug development process for diseases associated with aging An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:28.