The Role of Cholesterol in M2 Clustering and Viral Budding Explained.

The influenza A M2 homotetrameric channel consists of four transmembrane (TM) and four amphipathic helices (AHs). This viral proton channel is suggested to form clusters in the catenoid budding neck areas in raft-like domains of the plasma membrane, resulting in cell membrane scission and viral release. The channel clustering environment is rich in cholesterol. Previous experiments have shown that cholesterol significantly contributes to lipid bilayer undulations in viral buds. However, a clear explanation of membrane curvature from the distribution of cholesterol around the M2TM-AH clusters is lacking. Using coarse-grained molecular dynamics simulations of M2TM-AH in bilayers, we observed that M2 channels form specific, C2-symmetric, clusters with conical shapes driven by the attraction of their AHs. We showed that cholesterol stabilized the formation of M2 channel clusters by filling and bridging the conical gap between M2 channels at specific sites in the N-termini of adjacent channels or via the C-terminal region of TM and AHs, with the latter sites displaying a longer interaction time and higher stability. The potential of mean force calculations showed that when cholesterols occupy the identified interfacial binding sites between two M2 channels, the dimer is stabilized by 11 kJ/mol. This translates to the cholesterol-bound dimer being populated by almost 2 orders of magnitude compared to a dimer lacking cholesterol. We demonstrated that the cholesterol-bridged M2 channels can exert a lateral force on the surrounding membrane to induce the necessary negative Gaussian curvature profile, which permits spontaneous scission of the catenoid membrane neck and leads to viral buds and scission.

Conformational heterogeneity and structural features for function of the prototype viroporin influenza AM2.

The 97-residue influenza A matrix 2 (ΑM2) protein, a prototype for viroporins, transports protons through water molecules and His37. We discuss structural biology and molecular biophysics experiments and some functional assays that have transformed over 40 years our understanding of the structure and function of AM2. The structural studies on ΑM2 have been performed with different conditions (pH, temperature, lipid, constructs) and using various protein constructs, e.g., AM2 transmembrane (AM2TM) domain, AM2 conductance domain (AM2CD), ectodomain-containing or ectodomain-truncated, AM2 full length (AM2FL) and aimed to describe the different conformations and structural details that are necessary for the stability and function of AM2. However, the conclusions from these experiments appeared sometimes ambiguous and caused exciting debates. This was not due to inaccurate measurements, but instead because of the different membrane mimetic environment used, e.g., detergent, micelles or phospholipid bilayer, the method (e.g., X-ray crystallography, solid state NMR, solution NMR, native mass spectrometry), the used protein construct (e.g., AM2TM or AM2CD), or the amino acids residues to follow observables (e.g., NMR chemical shifts). We present these results according to the different used biophysical methods, the research groups and often by keeping a chronological order for presenting the progress in the research. We discuss ideas for additional research on structural details of AM2 and how the present findings can be useful to explore new routes of influenza A inhibition. The AM2 research can provide inspiration to study other viroporins as drug targets.

The Tuberculosis Drug Candidate SQ109 and Its Analogs Have Multistage Activity against Plasmodium falciparum.

Toward repositioning the antitubercular clinical candidate SQ109 as an antimalarial, analogs were investigated for structure-activity relationships for activity against asexual blood stages of the human malaria parasite Plasmodium falciparum pathogenic forms, as well as transmissible, sexual stage gametocytes. We show that equipotent activity (IC50) in the 100-300 nM range could be attained for both asexual and sexual stages, with the activity of most compounds retained against a multidrug-resistant strain. The multistage activity profile relies on high lipophilicity ascribed to the adamantane headgroup, and antiplasmodial activity is critically dependent on the diamine linker. Frontrunner compounds showed conserved activity against genetically diverse southern African clinical isolates. We additionally validated that this series could block transmission to mosquitoes, marking these compounds as novel chemotypes with multistage antiplasmodial activity.

Hit-to-Lead Optimization of Heterocyclic Carbonyloxycarboximidamides as Selective Antagonists at Human Adenosine A3 Receptor.

Antagonism of the human adenosine A3 receptor (hA3R) has potential therapeutic application. Alchemical relative binding free energy calculations of K18 and K32 suggested that the combination of a 3-(2,6-dichlorophenyl)-isoxazolyl group with 2-pyridinyl at the ends of a carbonyloxycarboximidamide group should improve hA3R affinity. Of the 25 new analogues synthesized, 37 and 74 showed improved hA3R affinity compared to K18 (and K32). This was further improved through the addition of a bromine group to the 2-pyridinyl at the 5-position, generating compound 39. Alchemical relative binding free energy calculations, mutagenesis studies and MD simulations supported the compounds' binding pattern while suggesting that the bromine of 39 inserts deep into the hA3R orthosteric pocket, so highlighting the importance of rigidification of the carbonyloxycarboximidamide moiety. MD simulations highlighted the importance of rigidification of the carbonyloxycarboximidamide, while suggesting that the bromine of 39 inserts deep into the hA3R orthosteric pocket, which was supported through mutagenesis studies 39 also selectively antagonized endogenously expressed hA3R in nonsmall cell lung carcinoma cells, while pharmacokinetic studies indicated low toxicity enabling in vivo evaluation. We therefore suggest that 39 has potential for further development as a high-affinity hA3R antagonist.

Copper-Catalyzed One-Pot Synthesis of Thiazolidin-2-imines.

The synthesis of thiazolines, thiazolidines, and thiazolidinones has been extensively studied, due to their biological activity related to neurodegenerative diseases, such as Parkinson's and Alzheimer's, as well as their antiparasitic and antihypertensive properties. The closely related thiazolidin-2-imines have been studied less, and efficient strategies for synthesizing them, mainly based on the reaction of propargylamines with isothiocyanates, have been explored less. The use of one-pot approaches, providing modular, straightforward, and sustainable access to these compounds, has also received very little attention. Herein, we report a novel, one-pot, multicomponent, copper-catalyzed reaction among primary amines, ketones, terminal alkynes, and isothiocyanates, toward thiazolidin-2-imines bearing quaternary carbon centers on the five-membered ring, in good to excellent yields. Density functional theory calculations, combined with experimental mechanistic findings, suggest that the copper(I)-catalyzed reaction between the in situ-formed propargylamines and isothiocyanates proceeds with a lower energy barrier in the pathway leading to the S-cyclized product, compared to that of the N-cyclized one, toward the chemo- and regioselective formation of 5-exo-dig S-cyclized thiazolidin-2-imines.

Computational Workflow for Refining AlphaFold Models in Drug Design Using Kinetic and Thermodynamic Binding Calculations: A Case Study for the Unresolved Inactive Human Adenosine A3 Receptor.

A structure-based drug design pipeline that considers both thermodynamic and kinetic binding data of ligands against a receptor will enable the computational design of improved drug molecules. For unresolved GPCR-ligand complexes, a workflow that can apply both thermodynamic and kinetic binding data in combination with alpha-fold (AF)-derived or other homology models and experimentally resolved binding modes of relevant ligands in GPCR-homologs needs to be tested. Here, as test case, we studied a congeneric set of ligands that bind to a structurally unresolved G protein-coupled receptor (GPCR), the inactive human adenosine A3 receptor (hA3R). We tested three available homology models from which two have been generated from experimental structures of hA1R or hA2AR and one model was a multistate alphafold 2 (AF2)-derived model. We applied alchemical calculations with thermodynamic integration coupled with molecular dynamics (TI/MD) simulations to calculate the experimental relative binding free energies and residence time (τ)-random accelerated MD (τ-RAMD) simulations to calculate the relative residence times (RTs) for antagonists. While the TI/MD calculations produced, for the three homology models, good Pearson correlation coefficients, correspondingly, r = 0.74, 0.62, and 0.67 and mean unsigned error (mue) values of 0.94, 1.31, and 0.81 kcal mol-1, the τ-RAMD method showed r = 0.92 and 0.52 for the first two models but failed to produce accurate results for the multistate AF2-derived model. With subsequent optimization of the AF2-derived model by reorientation of the side chain of R1735.34 located in the extracellular loop 2 (EL2) that blocked ligand's unbinding, the computational model showed r = 0.84 for kinetic data and improved performance for thermodynamic data (r = 0.81, mue = 0.56 kcal mol-1). Overall, after refining the multistate AF2 model with physics-based tools, we were able to show a strong correlation between predicted and experimental ligand relative residence times and affinities, achieving a level of accuracy comparable to an experimental structure. The computational workflow used can be applied to other receptors, helping to rank candidate drugs in a congeneric series and enabling the prioritization of leads with stronger binding affinities and longer residence times.

Accurate calculation of affinity changes to the close state of influenza A M2 transmembrane domain in response to subtle structural changes of adamantyl amines using free energy perturbation methods in different lipid bilayers.

Experimental binding free energies of 27 adamantyl amines against the influenza M2(22-46) WT tetramer, in its closed form at pH 8, were measured by ITC in DPC micelles. The measured Kd's range is ~44 while the antiviral potencies (IC50) range is ~750 with a good correlation between binding free energies computed with Kd and IC50 values (r = 0.76). We explored with MD simulations (ff19sb, CHARMM36m) the binding profile of complexes with strong, moderate and weak binders embedded in DMPC, DPPC, POPC or a viral mimetic membrane and using different experimental starting structures of M2. To predict accurately differences in binding free energy in response to subtle changes in the structure of the ligands, we performed 18 alchemical perturbative single topology FEP/MD NPT simulations (OPLS2005) using the BAR estimator (Desmond software) and 20 dual topology calculations TI/MD NVT simulations (ff19sb) using the MBAR estimator (Amber software) for adamantyl amines in complex with M2(22-46) WT in DMPC, DPPC, POPC. We observed that both methods with all lipids show a very good correlation between the experimental and calculated relative binding free energies (r = 0.77-0.87, mue = 0.36-0.92 kcal mol-1) with the highest performance achieved with TI/MBAR and lowest performance with FEP/BAR in DMPC bilayers. When antiviral potencies are used instead of the Kd values for computing the experimental binding free energies we obtained also good performance with both FEP/BAR (r = 0.83, mue = 0.75 kcal mol-1) and TI/MBAR (r = 0.69, mue = 0.77 kcal mol-1).

Structure-Based Lead Optimization of Enterovirus D68 2A Protease Inhibitors.

Enterovirus D68 (EV-D68) virus is a nonpolio enterovirus that typically causes respiratory illness and, in severe cases, can lead to paralysis and death in children. There is currently no vaccine or antiviral for EV-D68. We previously discovered the viral 2A protease (2Apro) as a viable antiviral drug target and identified telaprevir as a 2Apro inhibitor. 2Apro is a viral cysteine protease that cleaves the viral VP1-2A polyprotein junction. In this study, we report the X-ray crystal structures of EV-D68 2Apro, wild-type, and the C107A mutant and the structure-based lead optimization of telaprevir. Guided by the X-ray crystal structure, we predicted the binding pose of telaprevir in 2Apro using molecular dynamics simulations. We then utilized this model to inform structure-based optimization of the telaprevir's reactive warhead and P1-P4 substitutions. These efforts led to the discovery of 2Apro inhibitors with improved antiviral activity than telaprevir. These compounds represent promising lead compounds for further development as EV-D68 antivirals.

Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory.

We selected 145 reference organic molecules that include model fragments used in computer-aided drug design. We calculated 158 conformational energies and barriers using force fields, with wide applicability in commercial and free softwares and extensive application on the calculation of conformational energies of organic molecules, e.g. the UFF and DREIDING force fields, the Allinger's force fields MM3-96, MM3-00, MM4-8, the MM2-91 clones MMX and MM+, the MMFF94 force field, MM4, ab initio Hartree-Fock (HF) theory with different basis sets, the standard density functional theory B3LYP, the second-order post-HF MP2 theory and the Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory, with the latter used for accurate reference values. The data set of the organic molecules includes hydrocarbons, haloalkanes, conjugated compounds, and oxygen-, nitrogen-, phosphorus- and sulphur-containing compounds. We reviewed in detail the conformational aspects of these model organic molecules providing the current understanding of the steric and electronic factors that determine the stability of low energy conformers and the literature including previous experimental observations and calculated findings. While progress on the computer hardware allows the calculations of thousands of conformations for later use in drug design projects, this study is an update from previous classical studies that used, as reference values, experimental ones using a variety of methods and different environments. The lowest mean error against the DLPNO-CCSD(T) reference was calculated for MP2 (0.35 kcal mol-1), followed by B3LYP (0.69 kcal mol-1) and the HF theories (0.81-1.0 kcal mol-1). As regards the force fields, the lowest errors were observed for the Allinger's force fields MM3-00 (1.28 kcal mol-1), ΜΜ3-96 (1.40 kcal mol-1) and the Halgren's MMFF94 force field (1.30 kcal mol-1) and then for the MM2-91 clones MMX (1.77 kcal mol-1) and MM+ (2.01 kcal mol-1) and MM4 (2.05 kcal mol-1). The DREIDING (3.63 kcal mol-1) and UFF (3.77 kcal mol-1) force fields have the lowest performance. These model organic molecules we used are often present as fragments in drug-like molecules. The values calculated using DLPNO-CCSD(T) make up a valuable data set for further comparisons and for improved force field parameterization.

Molecular Biophysics of Class A G Protein Coupled Receptors-Lipids Interactome at a Glance-Highlights from the A2A Adenosine Receptor.

G protein-coupled receptors (GPCRs) are embedded in phospholipid membrane bilayers with cholesterol representing 34% of the total lipid content in mammalian plasma membranes. Membrane lipids interact with GPCRs structures and modulate their function and drug-stimulated signaling through conformational selection. It has been shown that anionic phospholipids form strong interactions between positively charged residues in the G protein and the TM5-TM6-TM 7 cytoplasmic interface of class A GPCRs stabilizing the signaling GPCR-G complex. Cholesterol with a high content in plasma membranes can be identified in more specific sites in the transmembrane region of GPCRs, such as the Cholesterol Consensus Motif (CCM) and Cholesterol Recognition Amino Acid Consensus (CRAC) motifs and other receptor dependent and receptor state dependent sites. Experimental biophysical methods, atomistic (AA) MD simulations and coarse-grained (CG) molecular dynamics simulations have been applied to investigate these interactions. We emphasized here the impact of phosphatidyl inositol-4,5-bisphosphate (PtdIns(4,5)P2 or PIP2), a minor phospholipid component and of cholesterol on the function-related conformational equilibria of the human A2A adenosine receptor (A2AR), a representative receptor in class A GPCR. Several GPCRs of class A interacted with PIP2 and cholesterol and in many cases the mechanism of the modulation of their function remains unknown. This review provides a helpful comprehensive overview for biophysics that enter the field of GPCRs-lipid systems.

A Study of the Activity of Adamantyl Amines against Mutant Influenza A M2 Channels Identified a Polycyclic Cage Amine Triple Blocker, Explored by Molecular Dynamics Simulations and Solid-State NMR.

We compared the anti-influenza potencies of 57 adamantyl amines and analogs against influenza A virus with serine-31 M2 proton channel, usually termed as WT M2 channel, which is amantadine sensitive. We also tested a subset of these compounds against viruses with the amantadine-resistant L26F, V27A, A30T, G34E M2 mutant channels. Four compounds inhibited WT M2 virus in vitro with mid-nanomolar potency, with 27 compounds showing sub-micromolar to low micromolar potency. Several compounds inhibited L26F M2 virus in vitro with sub-micromolar to low micromolar potency, but only three compounds blocked L26F M2-mediated proton current as determined by electrophysiology (EP). One compound was found to be a triple blocker of WT, L26F, V27A M2 channels by EP assays, but did not inhibit V27A M2 virus in vitro, and one compound inhibited WT, L26F, V27A M2 in vitro without blocking V27A M2 channel. One compound blocked only L26F M2 channel by EP, but did not inhibit virus replication. The triple blocker compound is as long as rimantadine, but could bind and block V27A M2 channel due to its larger girth as revealed by molecular dynamics simulations, while MAS NMR informed on the interaction of the compound with M2(18-60) WT or L26F or V27A.

Study of SQ109 analogs binding to mycobacterium MmpL3 transporter using MD simulations and alchemical relative binding free energy calculations.

N-geranyl-N΄-(2-adamantyl)ethane-1,2-diamine (SQ109) is a tuberculosis drug that has high potency against Mycobacterium tuberculosis (Mtb) and may function by blocking cell wall biosynthesis. After the crystal structure of MmpL3 from Mycobacterium smegmatis in complex with SQ109 became available, it was suggested that SQ109 inhibits Mmpl3 mycolic acid transporter. Here, we showed using molecular dynamics (MD) simulations that the binding profile of nine SQ109 analogs with inhibitory potency against Mtb and alkyl or aryl adducts at C-2 or C-1 adamantyl carbon to MmpL3 was consistent with the X-ray structure of MmpL3 - SQ109 complex. We showed that rotation of SQ109 around carbon-carbon bond in the monoprotonated ethylenediamine unit favors two gauche conformations as minima in water and lipophilic solvent using DFT calculations as well as inside the transporter's binding area using MD simulations. The binding assays in micelles suggested that the binding affinity of the SQ109 analogs was increased for the larger, more hydrophobic adducts, which was consistent with our results from MD simulations of the SQ109 analogues suggesting that sizeable C-2 adamantyl adducts of SQ109 can fill a lipophilic region between Y257, Y646, F260 and F649 in MmpL3. This was confirmed quantitatively by our calculations of the relative binding free energies using the thermodynamic integration coupled with MD simulations method with a mean assigned error of 0.74 kcal mol-1 compared to the experimental values.

Comparative Study of Receptor-, Receptor State-, and Membrane-Dependent Cholesterol Binding Sites in A2A and A1 Adenosine Receptors Using Coarse-Grained Molecular Dynamics Simulations.

We used coarse-grained molecular dynamics (CG MD) simulations to study protein-cholesterol interactions for different activation states of the A2A adenosine receptor (A2AR) and the A1 adenosine receptor (A1R) and predict new cholesterol binding sites indicating amino acid residues with a high residence time in three biologically relevant membranes. Compared to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)-cholesterol and POPC-phosphatidylinositol-bisphosphate (PIP2)-cholesterol, the plasma mimetic membrane best described the cholesterol binding sites previously detected for the inactive state of A2AR and revealed the binding sites with long-lasting amino acid residues. We observed that using the plasma mimetic membrane and plotting residues with cholesterol residence time ≥2 µs, our CG MD simulations captured most obviously the cholesterol-protein interactions. For the inactive A2AR, we identified one more binding site in which cholesterol is bound to residues with a long residence time compared to the previously detected, for the active A1R, three binding sites, and for the inactive A1R, two binding sites. We calculated that for the active states, cholesterol binds to residues with a much longer residence time compared to the inactive state for both A2AR and A1R. The stability of the identified binding sites to A1R or A2AR with CG MD simulations was additionally investigated with potential of mean force calculations using umbrella sampling. We observed that the binding sites with residues to which cholesterol has a long residence time in A2AR have shallow binding free energy minima compared to the related binding sites in A1R, suggesting a stronger binding for cholesterol to A1R. The differences in binding sites in which cholesterol is stabilized and interacts with residues with a long residence time between active and inactive states of A1R and A2AR can be important for differences in functional activity and orthosteric agonist or antagonist affinity and can be used for the design of allosteric modulators, which can bind through lipid pathways. We observed a stronger binding for cholesterol to A1R (i.e., generally higher association rates) compared to A2AR, which remains to be demonstrated. For the active states, cholesterol binds to residues with much longer residence times compared to the inactive state for both A2AR and A1R. Taken together, binding sites of active A1R may be considered as promising allosteric targets.

Synthesis and Testing of Analogs of the Tuberculosis Drug Candidate SQ109 against Bacteria and Protozoa: Identification of Lead Compounds against Mycobacterium abscessus and Malaria Parasites.

SQ109 is a tuberculosis drug candidate that has high potency against Mycobacterium tuberculosis and is thought to function at least in part by blocking cell wall biosynthesis by inhibiting the MmpL3 transporter. It also has activity against bacteria and protozoan parasites that lack MmpL3, where it can act as an uncoupler, targeting lipid membranes and Ca2+ homeostasis. Here, we synthesized 18 analogs of SQ109 and tested them against M. smegmatis, M. tuberculosis, M. abscessus, Bacillus subtilis, and Escherichia coli, as well as against the protozoan parasites Trypanosoma brucei, T. cruzi, Leishmania donovani, L. mexicana, and Plasmodium falciparum. Activity against the mycobacteria was generally less than with SQ109 and was reduced by increasing the size of the alkyl adduct, but two analogs were ∼4-8-fold more active than SQ109 against M. abscessus, including a highly drug-resistant strain harboring an A309P mutation in MmpL3. There was also better activity than found with SQ109 with other bacteria and protozoa. Of particular interest, we found that the adamantyl C-2 ethyl, butyl, phenyl, and benzyl analogs had 4-10× increased activity against P. falciparum asexual blood stages, together with low toxicity to a human HepG2 cell line, making them of interest as new antimalarial drug leads. We also used surface plasmon resonance to investigate the binding of inhibitors to MmpL3 and differential scanning calorimetry to investigate binding to lipid membranes. There was no correlation between MmpL3 binding and M. tuberculosis or M. smegmatis cell activity, suggesting that MmpL3 is not a major target in mycobacteria. However, some of the more active species decreased lipid phase transition temperatures, indicating increased accumulation in membranes, which is expected to lead to enhanced uncoupler activity.

A proof-of-concept study of the secondary structure of influenza A, B M2 and MERS- and SARS-CoV E transmembrane peptides using folding molecular dynamics simulations in a membrane mimetic solvent.

Here, we have carried out a proof-of-concept molecular dynamics (MD) simulation with adaptive tempering in a membrane mimetic environment to study the folding of single-pass membrane peptides. We tested the influenza A M2 viroporin, influenza B M2 viroporin, and protein E from coronaviruses MERS-Cov-2 and SARS-CoV-2 peptides with known experimental secondary structures in membrane bilayers. The two influenza-derived peptides are significantly different in the peptide sequence and secondary structure and more polar than the two coronavirus-derived peptides. Through a total of more than 50 µs of simulation time that could be accomplished in trifluoroethanol (TFE), as a membrane model, we characterized comparatively the folding behavior, helical stability, and helical propensity of these transmembrane peptides that match perfectly their experimental secondary structures, and we identified common motifs that reflect their quaternary organization and known (or not) biochemical function. We showed that BM2 is organized into two structurally distinct parts: a significantly more stable N-terminal half, and a fast-converting C-terminal half that continuously folds and unfolds between α-helical structures and non-canonical structures, which are mostly turns. In AM2, both the N-terminal half and C-terminal half are very flexible. In contrast, the two coronavirus-derived transmembrane peptides are much more stable and fast helix-formers when compared with the influenza ones. In particular, the SARS-derived peptide E appears to be the fastest and most stable helix-former of all the four viral peptides studied, with a helical structure that persists almost without disruption for the whole of its 10 µs simulation. By comparing the results with experimental observations, we benchmarked TFE in studying the conformation of membrane and hydrophobic peptides. This work provided accurate results suggesting a methodology to run long MD simulations and predict structural properties of biologically important membrane peptides.

Dual A1/A3 Adenosine Receptor Antagonists: Binding Kinetics and Structure-Activity Relationship Studies Using Mutagenesis and Alchemical Binding Free Energy Calculations.

Drugs targeting adenosine receptors (AR) can provide treatment for diseases. We report the identification of 7-(phenylamino)-pyrazolo[3,4-c]pyridines L2-L10, A15, and A17 as low-micromolar to low-nanomolar A1R/A3R dual antagonists, with 3-phenyl-5-cyano-7-(trimethoxyphenylamino)-pyrazolo[3,4-c]pyridine (A17) displaying the highest affinity at both receptors with a long residence time of binding, as determined using a NanoBRET-based assay. Two binding orientations of A17 produce stable complexes inside the orthosteric binding area of A1R in molecular dynamics (MD) simulations, and we selected the most plausible orientation based on the agreement with alanine mutagenesis supported by affinity experiments. Interestingly, for drug design purposes, the mutation of L2506.51 to alanine increased the binding affinity of A17 at A1R. We explored the structure-activity relationships against A1R using alchemical binding free energy calculations with the thermodynamic integration coupled with the MD simulation (TI/MD) method, applied on the whole G-protein-coupled receptor-membrane system, which showed a good agreement (r = 0.73) between calculated and experimental relative binding free energies.

Synthesis and in vitro proof-of-concept studies on bispecific iron oxide magnetic nanoparticles targeting PSMA and GRP receptors for PET/MR imaging of prostate cancer.

Prostate cancer (PCa) is the most common malignancy worldwide in men. This is a proof-of-concept study describing the development of 68Ga-magnetic iron oxide nanoparticles (mNP) targeting prostate specific membrane antigen (PSMA) and gastrin releasing peptide (GRPR) receptors as potential tools for diagnosis of PCa with PET/MRI. Two pharmacophores targeting PSMA, 1, and GRPR, 2, were coupled to mNPs carrying -SH (mNP-S1/2) or -NH2 (mNP-N1/2) groups. The mNP-S1/2 and mNP-N1/2 were characterized for their size, zeta potential, structure, and efficiency of functionalization using dynamic light scattering (DLS), FT-IR and RP-HPLC. A direct 68Ga-labelling procedure was followed, where 68Ga-mNP-N1/2 proved superior to 68Ga-mNP-S1/2 regarding radiolabelling efficiency, and thus were further evaluated in vitro. Toxicity studies in PCa cells (LNCaP, PC-3) showed low toxicity, and minimal hemolysis of red blood cells. In vitro assays in cells expressing PSMA (LNCaP), and GRPR (PC-3), showed specific time-dependent binding (40 min to plateau), high avidity (PC-3: Kd = 28.27 nM, LNCaP: Kd = 11.49 nM) and high internalization rates for 68Ga-mNP-N1/2 in both cell lines.

Discovery of a High Affinity Adenosine A1/A3 Receptor Antagonist with a Novel 7-Amino-pyrazolo[3,4-d]pyridazine Scaffold.

Here we describe the design and synthesis of pyrazolo[3,4-d]pyridazines as adenosine receptor (AR) ligands. We demonstrate that the introduction of a 3-phenyl group, together with a 7-benzylamino and 1-methyl group at the pyrazolopyridazine scaffold, generated the antagonist compound 10b, which displayed 21 nM affinity and a residence time of ∼60 min, for the human A1R, 55 nM affinity and a residence time of ∼73 min, for the human A3R and 1.7 µΜ affinity for the human A2BR while not being toxic. Strikingly, the 2-methyl analog of 10b, 15b, had no significant affinity. Docking calculations and molecular dynamics simulations of the ligands inside the orthosteric binding area suggested that the 2-methyl group in 15b hinders the formation of hydrogen bonding interactions with N6.55 which are considered critical for the stabilization inside the orthosteric binding cavity. We, therefore, demonstrate that 10a is a novel scaffold for the development of high affinity AR ligands. From the mutagenesis experiments the biggest effect was observed for the Y2717.46A mutation which caused an ∼10-fold reduction in the binding affinity of 10b.

Investigation of Tumor Cells and Receptor-Ligand Simulation Models for the Development of PET Imaging Probes Targeting PSMA and GRPR and a Possible Crosstalk between the Two Receptors.

Prostate-specific membrane antigen (PSMA) and gastrin-releasing peptide receptor (GRPR) have both been used in nuclear medicine as targets for molecular imaging and therapy of prostate (PCa) and breast cancer (BCa). Three bioconjugate probes, the PSMA specific: [68Ga]Ga-1, ((HBED-CC)-Ahx-Lys-NH-CO-NH Glu or PSMA-11), the GRPR specific: [68Ga]Ga-2, ((HBED-CC)-4-amino-1-carboxymethyl piperidine-[D-Phe6, Sta13]BN(6-14), a bombesin (BN) analogue), and 3 (the BN analogue: 4-amino-1-carboxymethyl piperidine-[(R)-Phe6, Sta13]BN(6-14) connected with the fluorescent dye, BDP-FL), were synthesized and tested in vitro with PCa and BCa cell lines, more specifically, with PCa cells, PC-3 and LNCaP, with BCa cells, T47D, MDA-MB-231, and with the in-house created PSMA-overexpressing PC-3(PSMA), T47D(PSMA), and MDA-MB-231(PSMA). In addition, biomolecular simulations were conducted on the association of 1 and 2 with PSMA and GRPR. The PSMA overexpression resulted in an increase of cell-bound radioligand [68Ga]Ga-1 (PSMA) for PCa and BCa cells and also of [68Ga]Ga-2 (GRPR), especially in those cell lines already expressing GRPR. The results were confirmed by fluorescence-activated cell sorting with a PE-labeled PSMA-specific antibody and the fluorescence tracer 3. The docking calculations and molecular dynamics simulations showed how 1 enters the PSMA funnel region and how pharmacophore Glu-urea-Lys interacts with the arginine patch, the S1', and S1 subpockets by forming hydrogen and van der Waals bonds. The chelating moiety of 1, that is, HBED-CC, forms additional stabilizing hydrogen bonding and van der Waals interactions in the arene-binding site. Ligand 2 is diving into the GRPR transmembrane (TM) helical cavity, thereby forming hydrogen bonds through its amidated end, water-mediated hydrogen bonds, and π-π interactions. Our results provide valuable information regarding the molecular mechanisms involved in the interactions of 1 and 2 with PSMA and GRPR, which might be useful for the diagnostic imaging and therapy of PCa and BCa.