Shedding Light on Dark Chemical Matter: The Discovery of a SARS-CoV-2 Mpro Main Protease Inhibitor through Intensive Virtual Screening and In Vitro Evaluation.

The development of specific antiviral therapies targeting SARS-CoV-2 remains fundamental because of the continued high incidence of COVID-19 and limited accessibility to antivirals in some countries. In this context, dark chemical matter (DCM), a set of drug-like compounds with outstanding selectivity profiles that have never shown bioactivity despite being extensively assayed, appears to be an excellent starting point for drug development. Accordingly, in this study, we performed a high-throughput screening to identify inhibitors of the SARS-CoV-2 main protease (Mpro) using DCM compounds as ligands. Multiple receptors and two different docking scoring functions were employed to identify the best molecular docking poses. The selected structures were subjected to extensive conventional and Gaussian accelerated molecular dynamics. From the results, four compounds with the best molecular behavior and binding energy were selected for experimental testing, one of which presented inhibitory activity with a Ki value of 48 ± 5 µM. Through virtual screening, we identified a significant starting point for drug development, shedding new light on DCM compounds.

A Step Forward toward Selective Activation/Inhibition of Bak, a Pro-Apoptotic Member of the Bcl-2 Protein Family: Discovery of New Prospective Allosteric Sites Using Molecular Dynamics.

Bak is a pro-apoptotic protein and a member of the Bcl-2 family that plays a key role in apoptosis, a programmed cell death mechanism of multicellular organisms. Its activation under death stimuli triggers the permeabilization of the mitochondrial outer membrane that represents a point of no return in the apoptotic pathway. This process is deregulated in many tumors where Bak is inactivated, whereas in other cases like in neurodegeneration, it exhibits an excessive response leading to disorders such as the Alzheimer disease. Members of the Bcl-2 family share a common 3D structure, exhibiting an extremely similar orthosteric binding site, a place where both pro and antiapoptotic proteins bind. This similarity raises a selectivity issue that hampers the identification of new drugs, capable of altering Bak activation in a selective manner. An alternative activation site triggered by antibodies has been recently identified, opening the opportunity to undertake new drug discovery studies. Despite this recent identification, an exhaustive study to identify cryptic pockets as prospective allosteric sites has not been yet performed. Thus, the present study aims to characterize novel hotspots in the Bak structure. For this purpose, we have carried out extensive molecular dynamics simulations using three different Bak systems including Bak in its apo form, Bak in complex with its endogen activator Bim and an intermediate form, set up by removing Bim from the previous complex. The results reported in the present work shed some light on future docking studies on Bak through the identification of new prospective allosteric sites, not previously described in this protein.

Discovery of Hit Compounds Targeting the P4 Allosteric Site of K-RAS, Identified through Ensemble-Based Virtual Screening.

Mutants of Ras are oncogenic drivers of a large number of human tumors. Despite being recognized as an attractive target for the treatment of cancer, the high affinity for its substrate tagged the protein as undruggable for a few years. The identification of cryptic pockets on the protein surface gave the opportunity to identify molecules capable of acting as allosteric modulators. Several molecules were disclosed in recent years, with sotorasib and adagrasib already approved for clinical use. The present study makes use of computational methods to characterize eight prospective allosteric pockets (P1-P8) in K-Ras, four of which (P1-P4) were previously characterized in the literature. The present study also describes the results of a virtual screening study focused on the discovery of hit compounds, binders of the P4 site that can be considered as peptidomimetics of a fragment of the SOS αI helix, a guanine exchange factor of Ras. After a detailed description of the computational procedure followed, we disclose five hit compounds, prospective binders of the P4 allosteric site that exhibit an inhibitory capability higher than 30% in a cell proliferation assay at 50 µM.

α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.

Molecular dynamics simulations of an α-synuclein NAC domain fragment with a ff14IDPSFF IDP-specific force field suggest ß-sheet intermediate states of fibrillation.

For the discovery of treatments against synucleinopathies, it is necessary to unravel and fully understand the mechanism of fibrillation of proteins involved. Among them, α-synuclein (αS) plays a key role in the development of these diseases through its aggregation into oligomers found in Lewy bodies. However, its structural disorder as an intrinsically disordered protein (IDP) makes its characterization by experimental techniques arduously difficult. Atomistic simulations aim to provide insights into this blank canvas and, fortunately, some studies have already suggested promising mechanisms. Still, it is urgent to consider the IDP features in simulations, so recently a lot of force fields designed to deal with IDPs have been developed. In this study, we have carried out a total of 12 µs simulations of an αS core fragment using a popular ff14SB AMBER force field and the ff14IDPSFF variation that includes a grid-based energy correction map (CMAP) method. The predicted chemical shifts from the simulations and those measured from the αS protein in the NMR solution indicate that ff14IDPSFF reproduces the experimental data more accurately. Moreover, structural analysis exhibits opposite trends between secondary structure propensities. The ff14SB force field preserves the α-helices found in the micelle-bound αS structure, which is used as an initial conformation, while ff14IDPSFF stands out with increased structural disorder and the formation of ß-sheets, which suggests that the IDP-specific force field can capture more suitable conformations representing the possible intermediate states of the fibrillation process.

Effects of solvents on the conformational profile of Balaram's peptide: a computational study.

The present work reports the results of a computational study aimed at characterizing the conformational profile of Balaram's peptide (Ace-Leu-Val-Val-Aib-Gly-Leu-Val-Val-NHMe) in different solvents, including chloroform, dimethyl sulfoxide, methanol and water. For this purpose, 10 µs molecular dynamics trajectories were computed in explicit solvents for each system, starting from an extended conformation. The results of the present study confirm the former NMR and CD findings and provide further insights that permit fine-tuning of the conclusions previously derived. The present results show that the peptide exhibits a helical conformation in chloroform, but a mixture of ß-hairpin and Ω-shape conformations, as the predominant structures in DMSO and MeOH. Finally, the peptide does not exhibit a preferred conformation in water, although significant populations of helical and ß-hairpin conformations are available. The present results underline the role of solvents in the conformational profile of a peptide and it is an example of the complementarity between computational methods and spectroscopy studies.

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 Mpro 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 Mpro 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-Mpro 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.

Fragment dissolved molecular dynamics: a systematic and efficient method to locate binding sites.

Diverse computational methods to support fragment-based drug discovery (FBDD) are available in the literature. Despite their demonstrated efficacy in supporting FBDD campaigns, they exhibit some drawbacks such as protein denaturation or ligand aggregation that have not yet been clearly overcome in the framework of biomolecular simulations. In the present work, we discuss a systematic semi-automatic novel computational procedure, designed to surpass these difficulties. The method, named fragment dissolved Molecular Dynamics (fdMD), utilizes simulation boxes of solvated small fragments, adding a repulsive Lennard-Jones potential term to avoid aggregation, which can be easily used to solvate the targets of interest. This method has the advantage of solvating the target with a low number of ligands, thus preventing the denaturation of the target, while simultaneously generating a database of ligand-solvated boxes that can be used in further studies. A number of scripts are made available to analyze the results and obtain the descriptors proposed as a means to trustfully discard spurious binding sites. To test our method, four test cases of different complexity have been solvated with ligand boxes and four molecular dynamics runs of 200 ns length have been run for each system, which have been extended up to 1 µs when needed. The reported results point out that the selected number of replicas are enough to identify the correct binding sites irrespective of the initial structure, even in the case of proteins having several close binding sites for the same ligand. We also propose a set of descriptors to analyze the results, among which the average MMGBSA and the average KDEEP energies have emerged as the most robust ones.

Molecular Determinants for the Activation/Inhibition of Bak Protein by BH3 Peptides.

Apoptosis is a key cell death pathway in mammalian cells. Understanding this process and its regulation has been a subject of study in the last three decades. Members of the Bcl-2 family of proteins are involved in the regulation of apoptosis through mitochondrial poration with the subsequent initiation of apoptosis. Deregulation of proapoptotic proteins contributes to the progression of many tumor processes. Understanding how these pore-forming Bcl-2 proteins Bak and Bax are activated is key to find new anticancer treatments. As no drug capable of activating Bak has been disclosed yet, the study of the structural features of BH3 peptides-known as Bak activators-relevant for binding along with its binding energy decomposition analysis, becomes essential for designing novel small-molecule mimics of BH3. Interestingly, a BH3 Bim analogue-inactivating Bak has recently been discovered, opening a question on the molecular features that determine the functions of BH3 peptides. Therefore, the present work is aimed at understanding the way BH3 peptides activate or inactivate Bak in order to identify differential structural features that can be used in drug design. For this purpose, complexes of Bak with an activator and an inhibitor have been subjected to a molecular dynamics study. Structural differences were assessed by means of the fluctuations of the corresponding principal components. Moreover, the MMPB/GBSA approach was used to compute the binding free energy of the diverse complexes to identify those residues of the BH3 peptide that exhibit the larger contributions to complex formation. The results obtained in this work show differences between activators and inhibitors, both in structural and energetic terms, which can be used in the design of new molecules that can activate or inactivate proapoptotic Bak.

Modeling and subtleties of K-Ras and Calmodulin interaction.

K-Ras, one of the most common small GTPases of the cell, still presents many riddles, despite the intense efforts to unveil its mysteries. Such is the case of its interaction with Calmodulin, a small acidic protein known for its role as a calcium ion sensor. Although the interaction between these two proteins and its biological implications have been widely studied, a model of their interaction has not been performed. In the present work we analyse this intriguing interaction by computational means. To do so, both conventional molecular dynamics and scaled molecular dynamics have been used. Our simulations suggest a model in which Calmodulin would interact with both the hypervariable region and the globular domain of K-Ras, using a lobe to interact with each of them. According to the presented model, the interface of helixes α4 and α5 of the globular domain of K-Ras would be relevant for the interaction with a lobe of Calmodulin. These results were also obtained when bringing the proteins together in a step wise manner with the umbrella sampling methodology. The computational results have been validated using SPR to determine the relevance of certain residues. Our results demonstrate that, when mutating residues of the α4-α5 interface described to be relevant for the interaction with Calmodulin, the interaction of the globular domain of K-Ras with Calmodulin diminishes. However, it is to be considered that our simulations indicate that the bulk of the interaction would fall on the hypervariable region of K-Ras, as many more interactions are identified in said region. All in all our simulations present a suitable model in which K-Ras could interact with Calmodulin at membrane level using both its globular domain and its hypervariable region to stablish an interaction that leads to an altered signalling.

Identification of Potential Small Molecule Binding Pockets in p38α MAP Kinase.

Given the essential role played by protein kinases in regulating cellular pathways, their dysregulation can result in the onset and/or progression of various human diseases. Structural analysis of diverse protein kinases suggests that these proteins exhibit a remarkable plasticity that allows them to adopt distinct conformations in response to interactions with other proteins, providing an opportunity for designing allosteric modulators. The present work reports the results of an in silico screening study aimed at identifying novel prospective allosteric binding sites in the paradigmatic p38α MAP kinase. The process was carried out using a protein ensemble generated from a 6 µs accelerated molecular dynamics simulation. The results of this calculation were first used to study the flexibility of the protein using Principal Component Analysis, followed by a Cluster Analysis aimed at producing an ensemble of conformations representative of the sampling process. Representative structures of the diverse clusters were subsequently screened for hot spots using FTMAP. The procedure permitted the identification of diverse allosteric sites of p38α already described in the literature including the DFG pocket, the lipid binding pocket, the DEF site, the docking groove, the CD and ED sites, and the backside site as well as a novel site recently reported: the A-loop regulatory site. Furthermore, the study also permitted the identification of ten novel prospective allosteric sites named NP1 to NP10, involving in most of the cases protein structural elements that control kinase activation including the activation loop, the catalytic loop, the αC helix, the L16 loop, and the glycine-rich loop.

Assessment of the Sampling Performance of Multiple-Copy Dynamics versus a Unique Trajectory.

The goal of the present study was to ascertain the differential performance of a long molecular dynamics trajectory versus several shorter ones starting from different points in the phase space and covering the same sampling time. For this purpose, we selected the 16-mer peptide Bak16BH3 as a model for study and carried out several samplings in explicit solvent. These samplings included an 8 µs trajectory (sampling S1); two 4 µs trajectories (sampling S2); four 2 µs trajectories (sampling S3); eight 1 µs trajectories (sampling S4); 16 0.5 µs trajectories (sampling S5), and 80 0.1 µs trajectories (sampling S6). Moreover, the 8 µs trajectory was further extended to 16 µs to have reference values of the diverse properties measured. The diverse samplings were compared qualitatively and quantitatively. Among the former, we carried out a comparison of the conformational profiles of the peptide using cluster analysis. Moreover, we also gained insight into the interchange among these structures along the sampling process. Among the latter, we computed the number of new conformational patterns sampled with time using strings defined from the conformations attained by each of the residues in the peptide. We also compared the locations and depths of the obtained minima on the free energy surface using principal component analysis. Finally, we also compared the helical profiles per residue at the end of the sampling process. The results suggest that a few short molecular dynamics trajectories may provide better sampling than one unique trajectory. Moreover, this procedure can also be advantageous to avoid getting trapped in a local minimum. However, caution should be exercised since short trajectories need to be long enough to overcome local barriers surrounding the starting point and the required sampling time depends on the number of degrees of freedom of the system under study. An effective way to gain insight into the minimum MD trajectory length is to monitor the convergence of different structural features, as shown in the present work.

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 topology applied to the discovery of 1-benzyl-2-(3-fluorophenyl)-4-hydroxy-3-(3-phenylpropanoyl)-2H-pyrrole-5-one as a non-ligand-binding-pocket antiandrogen.

We report the discovery of 1-benzyl-2-(3-fluorophenyl)-4-hydroxy-3-(3-phenylpropanoyl)-2H-pyrrole-5-one as a novel non-ligand binding pocket (non-LBP) antagonist of the androgen receptor (AR) through the application of molecular topology techniques. This compound, validated through time-resolved fluorescence resonance energy transfer and fluorescence polarization biological assays, provides the basis for lead optimization and structure-activity relationship analysis of a new series of non-LBP AR antagonists. Induced-fit docking and molecular dynamics studies have been performed to establish a consistent hypothesis for the interaction of the new active molecule on the AR surface.

Assessment of the bound conformation of bombesin to the BB1 and BB2 receptors.

Bombesin is an endogenous peptide involved in a wide spectrum of physiological activities ranging from satiety, control of circadian rhythm and thermoregulation in the central nervous system, to stimulation of gastrointestinal hormone release, activation of macrophages and effects on development in peripheral tissues. Actions of the peptide are mediated through the two high affinity G-protein coupled receptors BB1R and BB2R. Under pathophysiological conditions, these receptors are overexpressed in many different types of tumors, such as prostate cancer, breast cancer, small and non-small cell lung cancer and pancreatic cancer. This observation has been used for designing cell markers, but it has not been yet exploited for therapeutical purposes. Despite the enormous biological interest of the peptide, little is known about the stereochemical features that contribute to their activity. On the one hand, mutagenesis studies identified a few receptor residues important for high bombesin affinity and on the other, a few studies focused on the relevance of diverse residues of the peptide for receptor activation. Models of the peptide bound to BB1R and BB2R can be helpful to improve our understanding of the stereochemical features granting bombesin activity. Accordingly, the present study describes the computational process followed to construct such models by means of Steered Molecular Dynamics, using models of the peptide and its receptors. Present results provide new insights into the structure-activity relationships of bombesin and its receptors, as well as render an explanation for the differential binding affinity observed towards BB1R and BB2R. Finally, these models can be further exploited to help for designing novel small molecule peptidomimetics with improved pharmacokinetics profile.

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.

Autochthonous Peruvian Natural Plants as Potential SARS-CoV-2 Mpro Main Protease Inhibitors.

Over 750 million cases of COVID-19, caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), have been reported since the onset of the global outbreak. The need for effective treatments has spurred intensive research for therapeutic agents based on pharmaceutical repositioning or natural products. In light of prior studies asserting the bioactivity of natural compounds of the autochthonous Peruvian flora, the present study focuses on the identification SARS-CoV-2 Mpro main protease dimer inhibitors. To this end, a target-based virtual screening was performed over a representative set of Peruvian flora-derived natural compounds. The best poses obtained from the ensemble molecular docking process were selected. These structures were subjected to extensive molecular dynamics steps for the computation of binding free energies along the trajectory and evaluation of the stability of the complexes. The compounds exhibiting the best free energy behaviors were selected for in vitro testing, confirming the inhibitory activity of Hyperoside against Mpro, with a Ki value lower than 20 µM, presumably through allosteric modulation.

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, Arg761. We also corroborate the relevance of residues Arg761 and Tyr764 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.

Molecular determinants of Bim(BH3) peptide binding to pro-survival proteins.

Proteins of the Bcl-2 family regulate apoptosis through the formation of heterodimers between antiapoptotic or pro-survival proteins and proapoptotic or pro-death proteins. Overexpression of antiapoptotic proteins not only contributes to the progression of many cancers, but also confers resistance to the chemo- and radiotherapeutic treatments. It has been demonstrated that peptides containing the BH3 domain of proapoptotic Bcl-2 family members are able to bind and inhibit antiapoptotic proteins. For this reason, the design of small molecules mimicking the BH3 domain of proapoptotic proteins has emerged as a promising therapeutic strategy for cancer treatment during the last years. However, BH3 domains exhibit different affinities for binding to antiapoptotic proteins; whereas Bim(BH3) and Puma(BH3) are able to bind all antiapoptotic proteins, others like Bad(BH3) and Bmf(BH3) show preference for some proteins over others. Consequently, the ability of a BH3-mimetic to kill tumor cells will depend on the BH3 peptide used as template and thus will have a selective or pan-inhibition effect. Recently, it has been suggested that this last approach could be interesting. Therefore, the present work is aimed to elucidate how the nonselective peptide Bim(BH3) is able to bind to all of the Bcl-2 family antiapoptotic proteins. To unravel the molecular determinants of this pan-inhibition, we used the MM-PB/GBSA approaches to calculate the binding free energy of the different complexes studied and to determine which residues of the peptide have the largest contribution to complex formation. Results obtained in the present work show that the binding of Bim(BH3) to pro-survival proteins is mainly hydrophobic and that specific interactions are fully distributed along the peptide sequence.

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.