The Discovery of inhibitors of the SARS-CoV-2 S protein through computational drug repurposing

SARS-CoV-2 must bind its principal receptor, ACE2, on the target cell to initiate infection. This interaction is largely driven by the receptor binding domain (RBD) of the viral Spike (S) protein. Accordingly, antiviral compounds that can block RBD/ACE2 interactions can constitute promising antiviral agents. To identify such molecules, we performed a virtual screening of the Selleck FDA approved drugs and the Selleck database of Natural Products using a multistep computational procedure. An initial set of candidates was identified from an ensemble docking process using representative structures determined from the analysis of four 3 μ s molecular dynamics trajectories of the RBD/ACE2 complex. Two procedures were used to construct an initial set of candidates including a standard and a pharmacophore guided docking procedure. The initial set was subsequently subjected to a multistep sieving process to reduce the number of candidates to be tested experimentally, using increasingly demanding computational procedures, including the calculation of the binding free energy computed using the MMPBSA and MMGBSA methods. After the sieving process, a final list of 10 candidates was proposed, compounds which were subsequently purchased and tested ex-vivo. The results identified estradiol cypionate and telmisartan as inhibitors of SARS-CoV-2 entry into cells. Our findings demonstrate that the methodology presented here enables the discovery of inhibitors targeting viruses for which high-resolution structures are available.

A hotspot for posttranslational modifications on the androgen receptor dimer interface drives pathology and anti-androgen resistance

Mutations of the androgen receptor (AR) associated with prostate cancer and androgen insensitivity syndrome may profoundly influence its structure, protein interaction network, and binding to chromatin, resulting in altered transcription signatures and drug responses. Current structural information fails to explain the effect of pathological mutations on AR structure-function relationship. Here, we have thoroughly studied the effects of selected mutations that span the complete dimer interface of AR ligand-binding domain (AR-LBD) using x-ray crystallography in combination with in vitro, in silico, and cell-based assays. We show that these variants alter AR-dependent transcription and responses to anti-androgens by inducing a previously undescribed allosteric switch in the AR-LBD that increases exposure of a major methylation target, Arg⁷⁶¹. We also corroborate the relevance of residues Arg⁷⁶¹ and Tyr⁷⁶⁴ for AR dimerization and function. Together, our results reveal allosteric coupling of AR dimerization and posttranslational modifications as a disease mechanism with implications for precision medicine.

α4-α5 Helices on Surface of KRAS Can Accommodate Small Compounds That Increase KRAS Signaling While Inducing CRC Cell Death

KRAS is the most frequently mutated oncogene associated with the genesis and progress of pancreatic, lung and colorectal (CRC) tumors. KRAS has always been considered as a therapeutic target in cancer but until now only two compounds that inhibit one specific KRAS mutation have been approved for clinical use. In this work, by molecular dynamics and a docking process, we describe a new compound (P14B) that stably binds to a druggable pocket near the α4-α5 helices of the allosteric domain of KRAS. This region had previously been identified as the binding site for calmodulin (CaM). Using surface plasmon resonance and pulldown analyses, we prove that P14B binds directly to oncogenic KRAS thus competing with CaM. Interestingly, P14B favors oncogenic KRAS interaction with BRAF and phosphorylated C-RAF, and increases downstream Ras signaling in CRC cells expressing oncogenic KRAS. The viability of these cells, but not that of the normal cells, is impaired by P14B treatment. These data support the significance of the α4-α5 helices region of KRAS in the regulation of oncogenic KRAS signaling, and demonstrate that drugs interacting with this site may destine CRC cells to death by increasing oncogenic KRAS downstream signaling.

Shedding light on the binding mechanism of kinase inhibitors BI-2536, Volasetib and Ro-3280 with their pharmacological target PLK1

In the present work, the interactions of the novel kinase inhibitors BI-2536, Volasetib (BI-6727) and Ro-3280 with the pharmacological target PLK1 have been studied by fluorescence spectroscopy and molecular dynamics calculations. High Stern-Volmer constants were found in fluorescence experiments suggesting the formation of stable protein-ligand complexes. In addition, it was observed that the binding constant between BI-2536 and PLK1 increases about 100-fold in presence of the phosphopeptide Cdc25C-p that docks to the polo box domain of the protein and releases the kinase domain. All the determined binding constants are higher for the kinase inhibitors than for their competitor for the active center (ATP) being BI-2536 and Volasertib the inhibitors that showed more affinity for PLK1. Calculated binding free energies confirmed the higher affinity of PLK1 for BI-2536 and Volasertib than for ATP. The higher affinity of the inhibitors to PLK1 compared to ATP was mainly attributed to stronger van der Waals interactions. Results may help with the challenge of designing and developing new kinase inhibitors more effective in clinical cancer therapy.

Discovery of Diverse Natural Products as Inhibitors of SARS-CoV-2 Mpro Protease through Virtual Screening

SARS-CoV-2 is a type of coronavirus responsible for the international outbreak of respiratory illness termed COVID-19 that forced the World Health Organization to declare a pandemic infectious disease situation of international concern at the beginning of 2020. The need for a swift response against COVID-19 prompted to consider different sources to identify bioactive compounds that can be used as therapeutic agents, including available drugs and natural products. Accordingly, this work reports the results of a virtual screening process aimed at identifying antiviral natural product inhibitors of the SARS-CoV-2 M^pro viral protease. For this purpose, ca. 2000 compounds of the Selleck database of Natural Compounds were the subject of an ensemble docking process targeting the M^pro protease. Molecules that showed binding to most of the protein conformations were retained for a further step that involved the computation of the binding free energy of the ligand-M^pro complex along a molecular dynamics trajectory. The compounds that showed a smooth binding free energy behavior were selected for in vitro testing. From the resulting set of compounds, five compounds exhibited an antiviral profile, and they are disclosed in the present work.

Substituted tetrahydroisoquinolines: synthesis, characterization, antitumor activity and other biological properties

This work deals with the molecular design, synthesis and biological activity of a series of tetrahydro[1,4]dioxanisoquinolines and dimethoxyisoquinoline analogues. This study describes the synthesis strategy of these potential antitumor compounds, their multi-step synthesis and their optimization. A series of tetrahydroisoquinolines was synthesized and their cytotoxicity evaluated. Some of these tetrahydroisoquinolines showed promising KRas inhibition, antiangiogenesis activity and antiosteoporosis properties. Molecular modeling studies showed that compound 12 bind in the p1 pocket of the KRas protein making interactions with the hydrophobic residues Leu56, Tyr64, Tyr71 and Thr74 and hydrogen bonds with residues Glu37 and Asp38.

Binding of the anticancer drug BI-2536 to human serum albumin. A spectroscopic and theoretical study

BI-2536 is a potent Polo-like kinase inhibitor which induces apoptosis in diverse human cancer cell lines. The binding affinity of BI-2536 for human serum albumin (HSA) protein may define its pharmacokinetic and pharmacodynamic profile. We have studied the binding of BI-2536 to HSA by means of different spectroscopic techniques and docking calculations. We have experimentally observed that the affinity of BI-2536 for HSA is higher than that of other common HSA binding drugs. Therefore, it can be postulated that the drug dose should be increased to achieve a certain concentration of free drug in plasma, although BI-2536 could also reach tumour tissues by uptaking HSA/BI-2536 complex. Only a single binding site on HSA has been observed for BI-2536 which seems to correspond to the subdomain IIA pocket. The formation of the HSA/BI-2536 complex is a spontaneous and entropy-driven process that does not cause a significant change of the secondary structure of the protein. Its endothermic character could be related to proton release. Thermodynamic analysis showed that the main protein-drug interactions are of the van der Waals type although the presence of amide and ether groups in BI-2536 could also allow H-bonding with some residues in the subdomain IIA pocket.

Atomic level rendering of DNA-drug encounter

Computer simulations have been demonstrated to be important for unraveling atomic mechanisms in biological systems. In this study, we show how combining unbiased molecular dynamic simulations with appropriate analysis tools can successfully describe metal-based drug interactions with DNA. To elucidate the noncovalent affinity of cisplatin's family to DNA, we performed extensive all-atom molecular dynamics simulations (3.7 μs total simulation length). The results show that the parent drug, cisplatin, has less affinity to form noncovalent adducts in the major groove than its aquo complexes. Furthermore, the relative position in which the drugs enter the major groove is dependent on the compound's net charge. Based on the simulations, we estimated noncovalent binding free energies through the use of Markov state models. In addition, and to overcome the lack of experimental information, we employed two additional methods: Molecular Mechanics Poisson-Boltzmann Surface Area (MMPB-SA) and steered molecular dynamics with the Jarzynski estimator, with an overall good agreement between the three methods. All complexes show interaction energies below 3 kcal/mol with DNA but the charged hydrolysis products have slightly more favorable binding free energies than the parent drug. Moreover, this study sets the precedent for future unbiased DNA-ligand simulations of more complex binders.

Molecular dynamics analysis of the interaction between the human BCL6 BTB domain and its SMRT, NcoR and BCOR corepressors: the quest for a consensus dynamic pharmacophore

Targeting the BCL6 protein is a promising therapeutic strategy for the treatment of B cell lymphomas. One approach to treat these diseases consists of finding drug candidates able to disrupt the interactions established between BCL6 and its corepressors. Thus, this work presents a thorough comparative analysis of the interactions between the BCL6 BTB (bric-a-brac tramtrack broad complex) protein domain and its SMRT, NcoR and BCOR corepressor BBDs (BCL6 binding domain) through molecular dynamics. Moreover, a theoretical structure is presented and checked for the BCL6(BTB)-NcoR(BBD) complex. Considering the BBDs to be composed of 17 amino acids, our analyses show the region involving residues 4-15 of these 17 to play a main role in the protein-corepressor interactions. Particularly SER(11) seems to have a high relevance as it establishes specific bonds with BCL6(BTB) and is one of the only two residues sequence equivalent for the three studied corepressors. From this study, 14 pharmacophoric points have been proposed divided in two groups which coincide with residues 4-11 and 11-15, being SER(11) a hinge point. This finding suggests the possibility of searching for 2 small molecule inhibitors, mimicking 8 and 7 pharmacophoric points, respectively, which could incorporate a hydrogen donor pharmacophoric point mimicking SER(11) in any or both molecules. In short, the present work aims to contribute further knowledge in the modeling of drugs mimicking BCL6(BTB)-corepressor complexes.

Structural analysis of the inhibition of Cdk4 and Cdk6 by p16(INK4a) through molecular dynamics simulations

Cyclin-dependent kinases 4, 6 and 2 (Cdk4/6/2), are proteins that lead progression through the G1-S transition, a step strictly regulated in the process of cell proliferation. The p16(INK4a) tumor suppressor, whose expression is inhibited in a high number of cancers, binds to Cdk4/6 and inhibits phosphorylation of the retinoblastoma protein, forcing cells to remain in the G1 phase and therefore, arresting cell division. Accordingly, the design of small compounds mimicking the inhibition of p16(INK4a) appears to be a promising way to treat cancer. In order to get some insight into the key interactions governing recognition between different cyclin-dependent kinases and the p16(INK4a) tumor suppressor, the present work reports the results of molecular dynamics simulations of both, the Cdk6-p16(INK4a) complex and the Cdk4-p16(INK4a) complex, respectively at 300 K. Most of the key interactions observed, were already anticipated in the analysis of the crystal structure of Cdk6-p16(INK4a). However, a few different features found out from the analysis of these calculations provide a better understanding of the role of the T-loop conformation, a fragment of Cdks, and the way the ATP binding-site is distorted upon binding of p16(INK4a).