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.

Effect of Trace Metal Ions on the Conformational Stability of the Visual Photoreceptor Rhodopsin.

Trace metals are essential elements that play key roles in a number of biochemical processes governing human visual physiology in health and disease. Several trace metals, such as zinc, have been shown to play important roles in the visual phototransduction process. In spite of this, there has been little research conducted on the direct effect of trace metal elements on the visual photoreceptor rhodopsin. In the current study, we have determined the effect of several metal ions, such as iron, copper, chromium, manganese, and nickel, on the conformational stability of rhodopsin. To this aim, we analyzed, by means of UV-visible and fluorescence spectroscopic methods, the effects of these trace elements on the thermal stability of dark rhodopsin, the stability of its active Metarhodopsin II conformation, and its chromophore regeneration. Our results show that copper prevented rhodopsin regeneration and slowed down the retinal release process after illumination. In turn, Fe3+, but not Fe2+, increased the thermal stability of the dark inactive conformation of rhodopsin, whereas copper ions markedly decreased it. These findings stress the important role of trace metals in retinal physiology at the photoreceptor level and may be useful for the development of novel therapeutic strategies to treat retinal disease.

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.

Discovery of a Bradykinin B2 Partial Agonist Profile of Raloxifene in a Drug Repurposing Campaign.

Covid-19 urges a deeper understanding of the underlying molecular mechanisms involved in illness progression to provide a prompt therapeutical response with an adequate use of available drugs, including drug repurposing. Recently, it was suggested that a dysregulated bradykinin signaling can trigger the cytokine storm observed in patients with severe Covid-19. In the scope of a drug repurposing campaign undertaken to identify bradykinin antagonists, raloxifene was identified as prospective compound in a virtual screening process. The pharmacodynamics profile of raloxifene towards bradykinin receptors is reported in the present work, showing a weak selective partial agonist profile at the B2 receptor. In view of this new profile, its possible use as a therapeutical agent for the treatment of severe Covid-19 is discussed.

Structural Study of a New HIV-1 Entry Inhibitor and Interaction with the HIV-1 Fusion Peptide in Dodecylphosphocholine Micelles.

Previous studies support the hypothesis that the envelope GB virus C (GBV-C) E1 protein interferes the HIV-1 entry and that a peptide, derived from the region 139-156 of this protein, has been defined as a novel HIV-1 entry inhibitor. In this work, we firstly focus on the characterization of the structural features of this peptide, which are determinant for its anti-HIV-1 activity and secondly, on the study of its interaction with the proposed viral target (i.e., the HIV-1 fusion peptide). We report the structure of the peptide determined by NMR spectroscopy in dodecylphosphocholine (DPC) micelles solved by using restrained molecular dynamics calculations. The acquisition of different NMR experiments in DPC micelles (i.e., peptide-peptide titration, diffusion NMR spectroscopy, and addition of paramagnetic relaxation agents) allows a proposal of an inhibition mechanism. We conclude that a 18-mer peptide from the non-pathogenic E1 GBV-C protein, with a helix-turn-helix structure inhibits HIV-1 by binding to the HIV-1 fusion peptide at the membrane level, thereby interfering with those domains in the HIV-1, which are critical for stabilizing the six-helix bundle formation in a membranous environment.

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.

New insights into the stereochemical requirements of the bradykinin B2 receptor antagonists binding.

Bradykinin (BK) is a member of the kinin family, released in response to inflammation, trauma, burns, shock, allergy and some cardiovascular diseases, provoking vasodilatation and increased vascular permeability among other effects. Their actions are mediated through at least two G-protein coupled receptors, B1 a receptor up-regulated during inflammation episodes or tissue trauma and B2 that is constitutively expressed in a variety of cell types. The goal of the present work is to carry out a structure-activity study of BK B2 antagonism, taking into account the stereochemical features of diverse non-peptide antagonists and the way these features translate into ligand anchoring points to complementary regions of the receptor, through the analysis of the respective ligand-receptor complex. For this purpose an atomistic model of the BK B2 receptor was built by homology modeling and subsequently refined embedded in a lipid bilayer by means of a 600 ns molecular dynamics trajectory. The average structure from the last hundred nanoseconds of the molecular dynamics trajectory was energy minimized and used as model of the receptor for docking studies. For this purpose, a set of compounds with antagonistic profile, covering maximal diversity were selected from the literature. Specifically, the set of compounds include Fasitibant, FR173657, Anatibant, WIN64338, Bradyzide, CHEMBL442294, and JSM10292. Molecules were docked into the BK B2 receptor model and the corresponding complexes analyzed to understand ligand-receptor interactions. The outcome of this study is summarized in a 3D pharmacophore that explains the observed structure-activity results and provides insight into the design of novel molecules with antagonistic profile. To prove the validity of the pharmacophore hypothesized a virtual screening process was also carried out. The pharmacophore was used as query to identify new hits using diverse databases of molecules. The results of this study revealed a set of new hits with structures not connected to the molecules used for pharmacophore development. A few of these structures were purchased and tested. The results of the binding studies show about a 33% success rate with a correlation between the number of pharmacophore points fulfilled and their antagonistic potency. Some of these structures are disclosed in the present work.

Structure-Activity Studies of Novel Di-substituted [1,2,5]oxadiazolo [3,4-b]pyrazine Analogs Targeting the A-loop Regulatory Site of p38 MAP Kinase.

INTRODUCTION: In the quest for novel allosteric inhibitors of the p38 MAP kinase, we recently described the A-loop regulatory site, identified by means of molecular modeling studies together with the disclosure of a small molecule hit with a moderate inhibitory profile. Starting from this structure, we subsequently identified two additional hits with simpler molecular structures from an in silico screening study, using a substructure search in the SciFinder database. After corroboration of their inhibitory profile, analysis of their structures permitted to conclude about the suitability of the [1,2,5]oxadiazolo[3,4-b]pyrazine (furazano[ 3,4-b]pyrazine) scaffold for the development of potent A-loop regulatory site p38 MAP kinase inhibitors. Accordingly, we report the synthesis and pharmacological evaluation of a series of di-substituted analogs with a potent inhibitory profile of p38 MAP kinase, as shown by in vitro assays of their capability to inhibit IL-1ß secretion in human monocyte-derived macrophages. OBJECTIVE: To find small molecule potent inhibitors of the p38 MAP kinase A-loop regulatory site. METHODS: Starting from this structure, we subsequently identified two additional hits with simpler molecular structures from an in silico screening study, using a substructure search in the SciFinder database. After corroboration of their inhibitory profile, we carried out a hit-tolead optimization process guided by molecular modeling using a [1,2,5]oxadiazolo[3,4- b]pyrazine (furazano[3,4-b]pyrazine) scaffold. RESULTS: We report the synthesis and pharmacological evaluation of a series of di-substituted analogs with a potent inhibitory profile of p38 MAP kinase, as shown by in vitro assays of their capability to inhibit IL-1ß secretion in human monocyte-derived macrophages. CONCLUSION: We describe in the present work a series of [1,2,5]oxadiazolo[3,4-b]pyrazine (furazano[3,4-b]pyrazine), which are potent inhibitors of IL-1ß secretion in human monocytederived macrophages allosteric modulators of the p38 MAP kinase A-loop regulatory site.

Discovery of novel 2,3,5-trisubstituted pyridine analogs as potent inhibitors of IL-1ß via modulation of the p38 MAPK signaling pathway.

Interleukin-1ß is a central mediator of innate immune responses and inflammation. It plays a key role in a wide variety of pathologies, ranging from autoinflammatory diseases to metabolic syndrome and malignant tumors. It is well established that its inhibition results in a rapid and sustained reduction in disease severity, underlining the importance of having a repertoire of drugs of this class. At present, there are only three interleukin-1ß blockers approved in the clinic. All of them are biologics, requiring parenteral administration and resulting in expensive treatments. In an exercise to identify small molecule allosteric inhibitors of MAP kinases, we discovered a series of compounds that block IL-1ß release produced as a consequence of a stimulus involved in triggering an inflammatory response. The present study reports the hit-to-lead optimization process that permitted the identification of the compound 13b (AIK3-305) an orally available, potent and selective inhibitor of IL-1ß. Furthermore, the study also reports the results of an in vivo efficacy study of 13b in a LPS endotoxic shock model in male BALB/c mice, where IL-1ß inhibition is monitored in different tissues.

New Insights into the Stereochemical Requirements of the Bombesin BB1 Receptor Antagonists Binding.

Members of the family of bombesinlike peptides exert a wide range of biological activities both at the central nervous system and in peripheral tissues through at least three G-Protein Coupled Receptors: BB1, BB2 and BB3. Despite the number of peptide ligands already described, only a few small molecule binders have been disclosed so far, hampering a deeper understanding of their pharmacology. In order to have a deeper understanding of the stereochemical features characterizing binding to the BB1 receptor, we performed the molecular modeling study consisting of the construction of a 3D model of the receptor by homology modeling followed by a docking study of the peptoids PD168368 and PD176252 onto it. Analysis of the complexes permitted us to propose prospective bound conformations of the compounds, consistent with the experimental information available. Subsequently, we defined a pharmacophore describing minimal stereochemical requirements for binding to the BB1 receptor that was used in silico screening. This exercise yielded a set of small molecules that were purchased and tested, showing affinity to the BB1 but not to the BB2 receptor. These molecules exhibit scaffolds of diverse chemical families that can be used as a starting point for the development of novel BB1 antagonists.

Molecular Features of Non-Selective Small Molecule Antagonists of the Bradykinin Receptors.

Angiotensin converting enzyme 2 (ACE2) downregulation is a key negative factor for the severity of lung edema and acute lung failure observed in patients infected with SARS-CoV-2. ACE2 downregulation affects the levels of diverse peptide mediators of the renin-agiotensin-aldestosterone and kallikrein-kinin systems, compromising vascular hemostasis. Increasing evidence suggests that the inflammatory response observed in covid-19 patients is initiated by the action of kinins on the bradykinin receptors. Accordingly, the use of bradykinin antagonists should be considered as a strategy for therapeutic intervention against covid-19 illness progression. Presently, icatibant is the only bradykinin antagonist drug approved. In the present report, we investigated the molecular features characterizing non-selective antagonists targeting the bradykinin receptors and carried out a in silico screening of approved drugs, aimed at the identification of compounds with a non-selective bradykinin antagonist profile that can be evaluated for drug repurposing. The study permitted to identify eight compounds as prospective non-selective antagonists of the bradykinin receptors, including raloxifene; sildenafil; cefepime; cefpirome; imatinib; ponatinib; abemaciclib and entrectinib.

Assessment of the conformational profile of bombesin by computational methods.

In the present work, the results of a computational study aimed at assessing the conformational profile of bombesin are reported. The conformational space of the peptide was sampled by means of a 4 µs accelerated molecular dynamics simulation in water, using an explicit solvent model. The results were analyzed using Principal Component Analysis to get essential information on peptide fluctuations, along with cluster analysis to characterize different conformations in the sample. Analysis of the results suggests that the peptide adopts helical structures at the C-terminus that tend to unwind at the end of the peptide chain, since there are many structures exhibiting only two turns of a helix at the central segment of the peptide. In addition, the peptide also adopts hairpin turn structures at the N-terminus. Results of the simulation were confronted with available NMR results in a 2,2,2-trifluoroethanol/water (30% v/v) solution. Distances deduced form NOEs experiments only provide support to the presence of helical conformations that represent the most populated structures in the simulation. The absence of other conformations in the NMR experiments can be explained to be due to the α-helix enhancing nature of the solvent used in the experiments.

Importance of structure-based studies for the design of a novel HIV-1 inhibitor peptide.

Based on the structure of an HIV-1 entry inhibitor peptide two stapled- and a retro-enantio peptides have been designed to provide novel prevention interventions against HIV transmission. The three peptides show greater inhibitory potencies in cellular and mucosal tissue pre-clinical models than the parent sequence and the retro-enantio shows a strengthened proteolytic stability. Since HIV-1 fusion inhibitor peptides need to be embedded in the membrane to properly interact with their viral target, the structural features were determined by NMR spectroscopy in micelles and solved by using restrained molecular dynamics calculations. Both parent and retro-enantio peptides demonstrate a topology compatible with a shared helix-turn-helix conformation and assemble similarly in the membrane maintaining the active conformation needed for its interaction with the viral target site. This study represents a straightforward approach to design new targeted peptides as HIV-1 fusion inhibitors and lead us to define a retro-enantio peptide as a good candidate for pre-exposure prophylaxis against HIV-1.

Using molecular dynamics for the refinement of atomistic models of GPCRs by homology modeling.

Despite GPCRs sharing a common seven helix bundle, analysis of the diverse crystallographic structures available reveal specific features that might be relevant for ligand design. Despite the number of crystallographic structures of GPCRs steadily increasing, there are still challenges that hamper the availability of new structures. In the absence of a crystallographic structure, homology modeling remains one of the important techniques for constructing 3D models of proteins. In the present study we investigated the use of molecular dynamics simulations for the refinement of GPCRs models constructed by homology modeling. Specifically, we investigated the relevance of template selection, ligand inclusion as well as the length of the simulation on the quality of the GPCRs models constructed. For this purpose we chose the crystallographic structure of the rat muscarinic M3 receptor as reference and constructed diverse atomistic models by homology modeling, using different templates. Specifically, templates used in the present work include the human muscarinic M2; the more distant human histamine H1 and the even more distant bovine rhodopsin as shown in the GPCRs phylogenetic tree. We also investigated the use or not of a ligand in the refinement process. Hence, we conducted the refinement process of the M3 model using the M2 muscarinic as template with tiotropium or NMS docked in the orthosteric site and compared with the results obtained with a model refined without any ligand bound.

A multistep docking and scoring protocol for congeneric series: Implementation on kinase DFG-out type II inhibitors.

AIM: Rescoring of docking-binding poses can significantly improve molecular docking results. Our aim was to evaluate postprocessing docking protocols in order to determine the most suitable methodology for the study of the binding of congeneric compounds to protein kinases. MATERIALS & METHODS: Diverse ligand-receptor poses generated after docking were submitted to different relaxation protocols. The Molecular Mechanics Poisson-Boltzmann (Generalized Born) Surface Area approach was applied for the evaluation of the binding affinity of complexes obtained. The performance of various Molecular Mechanics Poisson-Boltzmann (Generalized Born) Surface Area methodologies was compared. RESULTS: The inclusion of a postprocessing protocol after docking enhances the quality of the results, although the best methodology is system dependent. CONCLUSION: An examination of the interactions established has allowed us to suggest useful modifications for the design of new type II inhibitors.

New insights into the stereochemical requirements of the bradykinin B1 receptor antagonists binding.

Bradykinin (BK) is a nonapeptide involved in several pathophysiological conditions including among others, septic and haemorrhagic shock, anaphylaxis, arthritis, rhinitis, asthma, inflammatory bowel disease. Accordingly, BK antagonists have long been sought after for therapeutic intervention. Action of BK is mediated through two different G-protein coupled receptors known as B1 and B2. Although there are several B1 antagonists reported in literature, their pharmacological profile is not yet optimal so that new molecules need to be discovered. In the present work we have constructed an atomistic model of the B1 receptor and docked diverse available non-peptide antagonists in order to get a deeper insight into the structure-activity relationships involving binding to this receptor. The model was constructed by homology modeling using the chemokine CXC4 and bovine rhodopsin receptors as template. The model was further refined using molecular dynamics for 600ns with the protein embedded in a POPC bilayer. From the refinement process we obtained an average structure that was used for docking studies using the Glide software. Antagonists selected for the docking studies include Compound 11, Compound 12, Chroman28, SSR240612, NPV-SAA164 and PS020990. The results of the docking study underline the role of specific receptor residues in ligand binding. The results of this study permitted to define a pharmacophore that describes the stereochemical requirements of antagonist binding, and can be used for the discovery of new compounds.

Identification of a Novel Inhibitory Allosteric Site in p38α.

In the present study, we report the discovery of a novel allosteric inhibitory site for p38α, a subclass of the mitogen-activated protein kinases (MAPK) family. The putative site was discovered after inspection of the crystallographic structure of the p38α-MK2 complex. MK2 (MAPK-activated protein kinase 2) is an interesting protein playing a dual role as modulator and substrate of p38α. This intriguing behavior is due to the ability of the two proteins to form distinctive heterodimers when p38α is phosphorylated or not. We hypothesized that the regulatory action of MK2 is due to its capability to keep p38α in an inactive conformation and consequently, we investigated the atomic structure of the p38α-MK2 complex to understand such regulatory behavior at the molecular level. After inspection of the complex structure, two peptides designed from the MK2 regulatory loop in contact with p38α with sequences Tyr1-Ser2-Asn3-His4-Gly5-Leu6 (peptide-1) and [Phe0]-peptide-1 (peptide-2) in their zwitterionic form were investigated for their phosphorylation inhibitory capability in vitro. Since both peptides exhibited inhibitory capability of the p38α kinase mediated phosphorylation of MEF2A, in a subsequent step we pursued the discovery of small molecule peptidomimetics. For this purpose we characterized in detail the peptide-p38α interaction using molecular dynamics simulations, leading to the definition of a pharmacophore for the peptide-protein interaction. This hypothesis was used as query for a in silico screening, leading to the discovery of a fused ring compound with micromolar inhibitory activity. Site-directed mutagenesis studies support that the compound binds to the putative novel allosteric site in p38α.

Blockade of the Interaction of Calcineurin with FOXO in Astrocytes Protects Against Amyloid-ß-Induced Neuronal Death.

Astrocytes actively participate in neuro-inflammatory processes associated to Alzheimer's disease (AD), and other brain pathologies. We recently showed that an astrocyte-specific intracellular signaling pathway involving an interaction of the phosphatase calcineurin with the transcription factor FOXO3 is a major driver in AD-associated pathological inflammation, suggesting a potential new druggable target for this devastating disease. We have now developed decoy molecules to interfere with calcineurin/FOXO3 interactions, and tested them in astrocytes and neuronal co-cultures exposed to amyloid-ß (Aß) toxicity. We observed that interference of calcineurin/FOXO3 interactions exerts a protective action against Aß-induced neuronal death and favors the production of a set of growth factors that we hypothesize form part of a cytoprotective pathway to resolve inflammation. Furthermore, interference of the Aß-induced interaction of calcineurin with FOXO3 by decoy compounds significantly decreased amyloid-ß protein precursor (AßPP) synthesis, reduced the AßPP amyloidogenic pathway, resulting in lower Aß levels, and blocked the expression of pro-inflammatory cytokines TNFα and IL-6 in astrocytes. Collectively, these data indicate that interrupting pro-inflammatory calcineurin/FOXO3 interactions in astrocytes triggered by Aß accumulation in brain may constitute an effective new therapeutic approach in AD. Future studies with intranasal delivery, or brain barrier permeable decoy compounds, are warranted.

Designing type II topoisomerase inhibitors: A molecular modeling approach.

Nosocomial infections are produced by pathogens with the ability to persist in hospital environments and with the propensity to develop resistance to diverse antimicrobials. In order to tackle resistance, it has been pointed as good strategy to select resilient drug targets that are evolutionally constrained to design drugs less susceptible to develop resistance. Molecular modeling can help to fulfill this goal by providing a rationalization of the observed resistance at the molecular level and, suggesting modifications on existing drugs or in the design of new ones to overcome the problem. The present report focus on type II topoisomerases, a clinical validated target for antibacterials and describe diverse modes of intervention including, inhibition of their ATPase function, stabilization of the cleavage complex or prevention of DNA strand hydrolysis. Moreover, the origin of resistance is also rationalized on the base of ligand-target interactions. Finally, efforts are described to circumvent the effect of non-susceptible strains by the design of new drugs based on existing ones, like the case of diones that act through the same mechanism as quinolones or the newly released quinole-carbonitrile derivatives that inhibit type II topoisomerases through a new mechanism.