Explore the vast range of titles available.

The process of ligand binding to a biological target can be represented as the equilibrium between the relevant solvated and bound states of the ligand. This which is the basis of structure-based, rigorous methods such as the estimation of relative binding affinities by free energy perturbation (FEP). Despite the growing capacity of computing power and the development of more accurate force fields, a high throughput application of FEP is currently hampered due to the need, in the current schemes, of an expert user definition of the “alchemical” transformations between molecules in the series explored. Here, we present QligFEP, a solution to this problem using an automated workflow for FEP calculations based on a dual topology approach. In this scheme, the starting poses of each of the two ligands, for which the relative affinity is to be calculated, are explicitly present in the MD simulations associated with the (dual topology) FEP transformation, making the perturbation pathway between the two ligands univocal. We show that this generalized method can be applied to accurately estimate solvation free energies for amino acid sidechain mimics, as well as the binding affinity shifts due to the chemical changes typical of lead optimization processes. This is illustrated in a number of protein systems extracted from other FEP studies in the literature: inhibitors of CDK2 kinase and a series of A 2A adenosine G protein-coupled receptor antagonists, where the results obtained with QligFEP are in excellent agreement with experimental data. In addition, our protocol allows for scaffold hopping perturbations to identify the binding affinities between different core scaffolds, which we illustrate with a series of Chk1 kinase inhibitors. QligFEP is implemented in the open-source MD package Q, and works with the most common family of force fields: OPLS, CHARMM and AMBER.

The four adenosine receptors (ARs) A 1 AR, A 2A AR, A 2B AR , and A 3 AR are G protein-coupled receptors (GPCRs) for which an exceptional amount of experimental and structural data is available. Still, limited success has been achieved in getting new chemical modulators on the market. As such, there is a clear interest in the design of novel selective chemical entities for this family of receptors. In this work, we investigate the selective recognition of ISAM-140, a recently reported A 2B AR reference antagonist. A combination of semipreparative chiral HPLC, circular dichroism and X-ray crystallography was used to separate and unequivocally assign the configuration of each enantiomer. Subsequently affinity evaluation for both A 2A and A 2B receptors demonstrate the stereospecific and selective recognition of ( S )-ISAM140 to the A 2B AR. The molecular modeling suggested that the structural determinants of this selectivity profile would be residue V250 6.51 in A 2B AR, which is a leucine in all other ARs including the closely related A 2A AR. This was herein confirmed by radioligand binding assays and rigorous free energy perturbation (FEP) calculations performed on the L249V 6.51 mutant A 2A AR receptor. Taken together, this study provides further insights in the binding mode of these A 2B AR antagonists, paving the way for future ligand optimization.

Transmembranal G Protein-Coupled Receptors (GPCRs) transduce extracellular chemical signals to the cell, via conformational change from a resting (inactive) to an active (canonically bound to a G-protein) conformation. Receptor activation is normally modulated by extracellular ligand binding, but mutations in the receptor can also shift this equilibrium by stabilizing different conformational states. In this work, we built structure-energetic relationships of receptor activation based on original thermodynamic cycles that represent the conformational equilibrium of the prototypical A 2A adenosine receptor (AR). These cycles were solved with efficient free energy perturbation (FEP) protocols, allowing to distinguish the pharmacological profile of different series of A 2A AR agonists with different efficacies. The modulatory effects of point mutations on the basal activity of the receptor or on ligand efficacies could also be detected. This methodology can guide GPCR ligand design with tailored pharmacological properties, or allow the identification of mutations that modulate receptor activation with potential clinical implications.

Generative deep learning models have emerged as a powerful approach for de novo drug design as they aid researchers in finding new molecules with desired properties. Despite continuous improvements in the field, a subset of the outputs that sequence-based de novo generators produce cannot be progressed due to errors. Here, we propose to fix these invalid outputs post hoc. In similar tasks, transformer models from the field of natural language processing have been shown to be very effective. Therefore, here this type of model was trained to translate invalid Simplified Molecular-Input Line-Entry System (SMILES) into valid representations. The performance of this SMILES corrector was evaluated on four representative methods of de novo generation: a recurrent neural network (RNN), a target-directed RNN, a generative adversarial network (GAN), and a variational autoencoder (VAE). This study has found that the percentage of invalid outputs from these specific generative models ranges between 4 and 89%, with different models having different error-type distributions. Post hoc correction of SMILES was shown to increase model validity. The SMILES corrector trained with one error per input alters 60–90% of invalid generator outputs and fixes 35–80% of them. However, a higher error detection and performance was obtained for transformer models trained with multiple errors per input. In this case, the best model was able to correct 60–95% of invalid generator outputs. Further analysis showed that these fixed molecules are comparable to the correct molecules from the de novo generators based on novelty and similarity. Additionally, the SMILES corrector can be used to expand the amount of interesting new molecules within the targeted chemical space. Introducing different errors into existing molecules yields novel analogs with a uniqueness of 39% and a novelty of approximately 20%. The results of this research demonstrate that SMILES correction is a viable post hoc extension and can enhance the search for better drug candidates. Graphical Abstract

Quantifying the effects of point mutations is of utmost interest for pharmaceutical and biotechnological applications. Reliable computational methods range from statistical and AI-based to physics-based approaches, with the optimal balance between accurate and fast predictions remaining a challenge. Free energy perturbation (FEP) simulations, a powerful physics-based approach available for decades, constitutes nowadays a method of common application in protein mutational studies. We present QresFEP-2, a novel hybrid-topology FEP protocol benchmarked on a comprehensive protein stability dataset of 10 protein systems, encompassing almost 600 mutations. QresFEP-2 combines excellent accuracy with the highest computational efficiency among available FEP protocols, and its robustness is further validated through comprehensive domain-wide mutagenesis, assessing the thermodynamic stability of over 400 mutations generated by a systematic mutation scan of the 56-residue B1 domain of streptococcal protein G (Gβ1). We also demonstrate the applicability domain of QresFEP-2 on evaluating site-directed mutagenesis effects on protein-ligand binding, tested on a GPCR, as well as on protein-protein interactions examined on the barnase/barstar complex. QresFEP-2 emerges as an open-source, physics-based alternative for advancing protein engineering, drug design, and elucidating the impact of mutations on human health. Understanding the effects of protein point mutations is crucial for pharmaceutical and biotechnological applications, yet achieving a balance between prediction accuracy and computational efficiency remains challenging. Here, the authors introduce QresFEP-2, a hybrid-topology free energy perturbation protocol that offers exceptional accuracy and efficiency, and robustly assesses effects of mutations on thermodynamic stability, protein-ligand binding and protein-protein interactions.

The four receptors that signal for adenosine, A1, A2A, A2B and A3 ARs, belong to the superfamily of G protein-coupled receptors (GPCRs). They mediate a number of (patho)physiological functions and have attracted the interest of the biopharmaceutical sector for decades as potential drug targets. The many crystal structures of the A2A, and lately the A1 ARs, allow for the use of advanced computational, structure-based ligand design methodologies. Over the last decade, we have assessed the efficient synthesis of novel ligands specifically addressed to each of the four ARs. We herein review and update the results of this program with particular focus on molecular dynamics (MD) and free energy perturbation (FEP) protocols. The first in silico mutagenesis on the A1AR here reported allows understanding the specificity and high affinity of the xanthine-antagonist 8-Cyclopentyl-1,3-dipropylxanthine (DPCPX). On the A2AAR, we demonstrate how FEP simulations can distinguish the conformational selectivity of a recent series of partial agonists. These novel results are complemented with the revision of the first series of enantiospecific antagonists on the A2BAR, and the use of FEP as a tool for bioisosteric design on the A3AR.

The adenosine A2A receptor (A2AAR) is a class A G-protein-coupled receptor (GPCR). It is an immune checkpoint in the tumor micro-environment and has become an emerging target for cancer treatment. In this study, we aimed to explore the effects of cancer-patient-derived A2AAR mutations on ligand binding and receptor functions. The wild-type A2AAR and 15 mutants identified by Genomic Data Commons (GDC) in human cancers were expressed in HEK293T cells. Firstly, we found that the binding affinity for agonist NECA was decreased in six mutants but increased for the V275A mutant. Mutations A165V and A265V decreased the binding affinity for antagonist ZM241385. Secondly, we found that the potency of NECA (EC50) in an impedance-based cell-morphology assay was mostly correlated with the binding affinity for the different mutants. Moreover, S132L and H278N were found to shift the A2AAR towards the inactive state. Importantly, we found that ZM241385 could not inhibit the activation of V275A and P285L stimulated by NECA. Taken together, the cancer-associated mutations of A2AAR modulated ligand binding and receptor functions. This study provides fundamental insights into the structure–activity relationship of the A2AAR and provides insights for A2AAR-related personalized treatment in cancer.

Overexpression of the adenosine A1 receptor (A1AR) has been detected in various cancer cell lines. However, the role of A1AR in tumor development is still unclear. Thirteen A1AR mutations were identified in the Cancer Genome Atlas from cancer patient samples. We have investigated the pharmacology of the mutations located at the 7-transmembrane domain using a yeast system. Concentration–growth curves were obtained with the full agonist CPA and compared to the wild type hA1AR. H78L3.23 and S246T6.47 showed increased constitutive activity, while only the constitutive activity of S246T6.47 could be reduced to wild type levels by the inverse agonist DPCPX. Decreased constitutive activity was observed on five mutant receptors, among which A52V2.47 and W188C5.46 showed a diminished potency for CPA. Lastly, a complete loss of activation was observed in five mutant receptors. A selection of mutations was also investigated in a mammalian system, showing comparable effects on receptor activation as in the yeast system, except for residues pointing toward the membrane. Taken together, this study will enrich the view of the receptor structure and function of A1AR, enlightening the consequences of these mutations in cancer. Ultimately, this may provide an opportunity for precision medicine for cancer patients with pathological phenotypes involving these mutations.

BackgroundAdenosine is a metabolite that suppresses antitumor immune response of T and NK cells via extracellular binding to the two subtypes of adenosine-2 receptors, A2ARs. While blockade of the A2AARs subtype effectively rescues lymphocyte activity, with four A2AAR antagonists currently in anticancer clinical trials, less is known for the therapeutic potential of the other A2BAR blockade within cancer immunotherapy. Recent studies suggest the formation of A2AAR/A2BAR dimers in tissues that coexpress the two receptor subtypes, where the A2BAR plays a dominant role, suggesting it as a promising target for cancer immunotherapy.MethodsWe report the synthesis and functional evaluation of five potent A2BAR antagonists and a dual A2AAR/A2BAR antagonist. The compounds were designed using previous pharmacological data assisted by modeling studies. Synthesis was developed using multicomponent approaches. Flow cytometry was used to evaluate the phenotype of T and NK cells on A2BAR antagonist treatment. Functional activity of T and NK cells was tested in patient-derived tumor spheroid models.ResultsWe provide data for six novel small molecules: five A2BAR selective antagonists and a dual A2AAR/A2BAR antagonist. The growth of patient-derived breast cancer spheroids is prevented when treated with A2BAR antagonists. To elucidate if this depends on increased lymphocyte activity, immune cells proliferation, and cytokine production, lymphocyte infiltration was evaluated and compared with the potent A2AAR antagonist AZD-4635. We find that A2BAR antagonists rescue T and NK cell proliferation, IFNγ and perforin production, and increase tumor infiltrating lymphocytes infiltration into tumor spheroids without altering the expression of adhesion molecules.ConclusionsOur results demonstrate that A2BAR is a promising target in immunotherapy, identifying ISAM-R56A as the most potent candidate for A2BAR blockade. Inhibition of A2BAR signaling restores T cell function and proliferation. Furthermore, A2BAR and dual A2AAR/A2BAR antagonists showed similar or better results than A2AAR antagonist AZD-4635 reinforcing the idea of dominant role of the A2BAR in the regulation of the immune system.

GPR139 is an orphan G protein-coupled receptor expressed in the brain, in particular in the habenula, hypothalamus and striatum. It has therefore been suggested that GPR139 is a possible target for metabolic disorders and Parkinson’s disease. Several surrogate agonist series have been published for GPR139. Two series published by Shi et al . and Dvorak et al . included agonists 1a and 7c respectively, with potencies in the ten-nanomolar range. Furthermore, Isberg et al . and Liu et al . have previously shown that tryptophan (Trp) and phenylalanine (Phe) can activate GPR139 in the hundred-micromolar range. In this study, we produced a mutagenesis-guided model of the GPR139 binding site to form a foundation for future structure-based ligand optimization. Receptor mutants studied in a Ca 2+ assay demonstrated that residues F109 3×33 , H187 5×43 , W241 6×48 and N271 7×38 , but not E108 3×32 , are highly important for the activation of GPR139 as predicted by the receptor model. The initial ligand-receptor complex was optimized through free energy perturbation simulations, generating a refined GPR139 model in agreement with experimental data. In summary, the GPR139 reference surrogate agonists 1a and 7c, and the endogenous amino acids l -Trp and l -Phe share a common binding site, as demonstrated by mutagenesis, ligand docking and free energy calculations.