Search

Psilocybin and its direct precursor baeocystin are indole alkaloids of psychotropic Psilocybe mushrooms. The pharmaceutical interest in psilocybin as a treatment option against depression and anxiety is currently being investigated in advanced clinical trials. Here, we report a biocatalytic route to synthesize 6-methylated psilocybin and baeocystin from 4-hydroxy-6-methyl-l-tryptophan, which was decarboxylated and phosphorylated by the Psilocybe cubensis biosynthesis enzymes PsiD and PsiK. N-Methylation was catalyzed by PsiM. We further present an in silico structural model of PsiM that revealed a well-conserved SAM-binding core along with peripheral nonconserved elements that likely govern substrate preferences.

Methane-producing archaea are among a select group of microorganisms that utilize tetrahydromethanopterin (H4MPT) as a one-carbon carrier instead of tetrahydrofolate. In H4MPT biosynthesis, β-ribofuranosylaminobenzene 5’-phosphate (RFAP) synthase catalyzes the production of RFAP, CO2, and pyrophosphate from p-aminobenzoic acid (pABA) and phosphoribosyl-pyrophosphate (PRPP). In this work, to gain insight into amino acid residues required for substrate binding, RFAP synthase from Methanothermobacter thermautotrophicus was produced in Escherichia coli, and site-directed mutagenesis was used to alter arginine 26 (R26) and aspartic acid 19 (D19), located in a conserved sequence of amino acids resembling the pABA binding site of dihydropteroate synthase. Replacement of R26 with lysine increased the KM for pABA by an order of magnitude relative to wild-type enzyme without substantially altering the KM for PRPP. Although replacement of D19 with alanine produced inactive enzyme, asparagine substitution allowed retention of some activity, and the KM for pABA increased about threefold relative to wild-type enzyme. A molecular model developed by threading RFAP synthase onto the crystal structure of homoserine kinase places R26 in the proposed active site. In the static model, D19 is located close to the active site, yet appears too far away to influence ligand binding directly. This may be indicative of the protein conformational change predicted previously in the Bi-Ter kinetic mechanism and/or formation of the active site at the interface of two subunits. Due to the vital role of RFAP synthase in H4MPT biosynthesis, insights into the mode of substrate binding and mechanism could be beneficial for developing RFAP synthase inhibitors designed to reduce the production of methane as a greenhouse gas.

Cyclophilins are peptidyl cis-trans prolyl isomerases (PPIases), whose activity is typically inhibited by cyclosporine A (CsA), a potent immunosuppressor. Cyclophilins are also chaperones. Emerging evidence supports that cyclophilins present non-overlapping PPIase and chaperone activities. The proteostasis of the disease-relevant substrates, signal transducer and activator of transcription 3 and 5 (STAT3/STAT5), heterogeneous nuclear ribonucleoprotein A2B1 (hnRNPA2B1) and M-opsin, are regulated by non-overlapping chaperone and PPIase activities of the cyclophilin domain (CY) of Ranbp2, a multifunctional and modular scaffold which controls nucleocytoplasmic shuttling and proteostasis of selective substrates. Although highly homologous, CY and the archetypal cyclophilin A (CyPA) present distinct catalytic and CsA-binding activities owing to unique structural features between these cylophilins. We explored structural idiosyncrasies between CY and CyPA to screen in silico nearly 9 million small molecules (SM) against the CY PPIase pocket and identify SMs with selective bioactivity toward STAT3, hnRNPA2B1 and/or M-opsin proteostasis. We found three classes of SMs that enhance the cytokine-stimulated transcriptional activity of STAT3 without changing latent and activated STAT3 levels, down-regulate hnRNPA2B1 or M-opsin proteostasis, or a combination of these. Further, a SM which suppresses hnRNPA2B1 proteostasis also inhibits strongly and selectively the PPIase activity of CY. This study unravels chemical probes for multimodal regulation of CY of Ranbp2 and its substrates and this regulation likely results in the allosterism stemming from the interconversion of conformational substates of cyclophilins. The results also demonstrate the feasibility of CY in drug discovery against disease-relevant substrates controlled by Ranbp2 and they open new opportunities for therapeutic interventions.

The immunophilins, cyclophilins, catalyze peptidyl cis-trans prolyl-isomerization (PPIase), a rate-limiting step in protein folding and a conformational switch in protein function. Cyclophilins are also chaperones. Noncatalytic mutations affecting the only cyclophilins with known but distinct physiological substrates, the Drosophila NinaA and its mammalian homolog, cyclophilin-B, impair opsin biogenesis and cause osteogenesis imperfecta, respectively. However, the physiological roles and substrates of most cyclophilins remain unknown. It is also unclear if PPIase and chaperone activities reflect distinct cyclophilin properties. To elucidate the physiological idiosyncrasy stemming from potential cyclophilin functions, we generated mice lacking endogenous Ran-binding protein-2 (Ranbp2) and expressing bacterial artificial chromosomes of Ranbp2 with impaired C-terminal chaperone and with (Tg-Ranbp2(WT-HA)) or without PPIase activities (Tg-Ranbp2(R2944A-HA)). The transgenic lines exhibit unique effects in proteostasis. Either line presents selective deficits in M-opsin biogenesis with its accumulation and aggregation in cone photoreceptors but without proteostatic impairment of two novel Ranbp2 cyclophilin partners, the cytokine-responsive effectors, STAT3/STAT5. Stress-induced STAT3 activation is also unaffected in Tg-Ranbp2(R2944A-HA)::Ranbp2(-/-). Conversely, proteomic analyses found that the multisystem proteinopathy/amyotrophic lateral sclerosis proteins, heterogeneous nuclear ribonucleoproteins A2/B1, are down-regulated post-transcriptionally only in Tg-Ranbp2(R2944A-HA)::Ranbp2(-/-). This is accompanied by the age- and tissue-dependent reductions of diubiquitin and ubiquitylated proteins, increased deubiquitylation activity, and accumulation of the 26 S proteasome subunits S1 and S5b. These manifestations are absent in another line, Tg-Ranbp2(CLDm-HA)::Ranbp2(-/-), harboring SUMO-1 and S1-binding mutations in the Ranbp2 cyclophilin-like domain. These results unveil distinct mechanistic and biological links between PPIase and chaperone activities of Ranbp2 cyclophilin toward proteostasis of selective substrates and with novel therapeutic potential.

Understanding the molecular basis of drug action can facilitate development of more potent and selective drugs. Here, we explore the molecular basis for action of a unique small molecule ligand that is a type 1 cholecystokinin (CCK) receptor agonist and type 2 CCK receptor antagonist, GI181771X. We characterize its binding utilizing structurally related radioiodinated ligands selective for CCK receptor subtypes that utilize the same allosteric ligand-binding pocket, using wild-type receptors and chimeric constructs exchanging the distinct residues lining this pocket. Intracellular calcium assays were performed to determine biological activity. Molecular models for docking small molecule agonists to the type 1 CCK receptor were developed using a ligand-guided refinement approach. The optimal model was distinct from the previous antagonist model for the same receptor and was mechanistically consistent with the current mutagenesis data. This study revealed a key role for Leu(7.39) that was predicted to interact with the isopropyl group in the N1 position of the benzodiazepine that acts as a "trigger" for biological activity. The molecular model was predictive of binding of other small molecule agonists, effectively distinguishing these from 1065 approved drug decoys with an area under curve value of 99%. The model also selectively enriched for agonist compounds, with 130 agonists identified by ROC analysis when seeded in 2175 non-agonist ligands of the type 1 CCK receptor (area under curve 78%). Benzodiazepine agonists in this series docked in consistent pose within this pocket, with a key role played by Leu(7.39), whereas the role of this residue was less clear for chemically distinct agonists.

Amino-terminal regions of secretin-family peptides contain key determinants for biological activity and binding specificity, although the nature of interactions with receptors is unclear. A helix N-capping motif within this region has been postulated to directly contribute to agonist activity while also stabilizing formation of a helix extending toward the peptide carboxyl terminus and docking within the receptor amino terminus. We used cysteine trapping to systematically explore spatial approximations between cysteines replacing each residue in this motif of secretin (sec), Phe(6), Thr(7), and Leu(10), and cysteines incorporated into the extracellular face of the receptor. Each peptide was a full agonist for cAMP, but had a lower binding affinity than natural hormone. These bound to COS cells expressing 61 receptor constructs incorporating cysteines in every position along each extracellular loop (ECL) and adjacent parts of transmembrane (TM) segments. Patterns of covalent labeling were distinct for each probe, with Cys(6)-sec labeling multiple residues in the carboxyl-terminal half of ECL2 and throughout ECL3, Cys(7)-sec predominantly labeling only single residues in the carboxyl-terminal end of ECL2 and the amino-terminal end of ECL3, and Cys(10)-sec not efficiently labeling any of these residues. These spatial constraints were used to refine our model of secretin bound to its receptor, now bringing ECL3 above the amino terminus of the ligand and revealing possible charge-charge interactions between this part of secretin and receptor residues in TM5, TM6, ECL2, and ECL3, which can orient and stabilize the peptide-receptor complex. This was validated by testing predicted approximations by mutagenesis and residue-residue complementation studies.

Allosteric binding pockets in peptide-binding G protein-coupled receptors create opportunities for the development of small molecule drugs with substantial benefits over orthosteric ligands. To gain insights into molecular determinants for this pocket within type 1 and 2 cholecystokinin receptors (CCK1R and CCK2R), we prepared a series of receptor constructs in which six distinct residues in TM2, -3, -6, and -7 were reversed. Two novel iodinated CCK1R- and CCK2R-selective 1,4-benzodiazepine antagonists, differing only in stereochemistry at C3, were used. When all six residues within CCK1R were mutated to corresponding CCK2R residues, benzodiazepine selectivity was reversed, yet peptide binding selectivity was unaffected. Detailed analysis, including observations of gain of function, demonstrated that residues 6.51, 6.52, and 7.39 were most important for binding the CCK1R-selective ligand, whereas residues 2.61 and 7.39 were most important for binding CCK2R-selective ligand, although the effect of substitution of residue 2.61 was likely indirect. Ligand-guided homology modeling was applied to wild type receptors and those reversing benzodiazepine binding selectivity. The models had high predictive power in enriching known receptor-selective ligands from related decoys, indicating a high degree of precision in pocket definition. The benzodiazepines docked in similar poses in both receptors, with C3 urea substituents pointing upward, whereas different stereochemistry at C3 directed the C5 phenyl rings and N1 methyl groups into opposite orientations. The geometry of the binding pockets and specific interactions predicted for ligand docking in these models provide a molecular framework for understanding ligand selectivity at these receptor subtypes. Furthermore, the strong predictive power of these models suggests their usefulness in the discovery of lead compounds and in drug development programs.

The molecular basis of ligand binding and activation of family B G protein-coupled receptors is not yet clear due to the lack of insight into the structure of intact receptors. Although NMR and crystal structures of amino-terminal domains of several family members support consistency in general structural motifs that include a peptide-binding cleft, there are variations in the details of docking of the carboxyl terminus of peptide ligands within this cleft, and there is no information about siting of the amino terminus of these peptides. There are also no empirical data to orient the receptor amino terminus relative to the core helical bundle domain. Here, we prepared a series of five new probes, incorporating photolabile moieties into positions 2, 15, 20, 24, and 25 of full agonist secretin analogues. Each bound specifically to the receptor and covalently labeled single distinct receptor residues. Peptide mapping of labeled wild-type and mutant receptors identified that the position 15, 20, and 25 probes labeled residues within the distal amino terminus of the receptor, whereas the position 24 probe labeled the amino terminus adjacent to TM1. Of note, the position 2 probe labeled a residue within the first extracellular loop of the receptor, a region not previously labeled, providing an important new constraint for docking the amino-terminal region of secretin to its receptor core. These additional experimentally derived constraints help to refine our understanding of the structure of the secretin-intact receptor complex and provide new insights into understanding the molecular mechanism for activation of family B G protein-coupled receptors.

The secretin receptor, a prototypic family B G protein-coupled receptor, forms a constitutive homodimeric complex that is stable even in the presence of hormone. Recently, a model of this agonist-bound receptor was built based on high resolution structures reported for amino-terminal domains of other family members. Although this model provided the best solution for all extant data, including 10 photoaffinity labeling constraints, a new such constraint now obtained with a position 16 photolabile probe was inconsistent with this model. As the secretin receptor forms constitutive homodimers, we explored whether secretin might dock across both protomers of the complex, an observation that could also contribute to the negative cooperativity observed. To directly explore this, we prepared six secretin analogue probes that simultaneously incorporated two photolabile benzoylphenylalanines as sites of covalent attachment, in positions known to label distinct receptor subdomains. Each bifunctional probe was a full agonist that labeled the receptor specifically and saturably, with electrophoretic migration consistent with labeling a single protomer of the homodimeric secretin receptor. No band representing radiolabeled receptor dimer was observed with any bifunctional probe. The labeled monomeric receptor bands were cleaved with cyanogen bromide to demonstrate that both of the photolabile benzoylphenylalanines within a single probe had established covalent adducts with a single receptor in the complex. These data are consistent with a model of secretin occupying a single secretin receptor protomer within the homodimeric receptor complex. A new molecular model accommodating all constraints is now proposed.

Afamin, a human plasma glycoprotein and putative transporter of hydrophobic molecules, has been shown to act as extracellular chaperone for poorly soluble, acylated Wnt proteins, forming a stable, soluble complex with functioning Wnt proteins. The 2.1-Å crystal structure of glycosylated human afamin reveals an almost exclusively hydrophobic binding cleft capable of harboring large hydrophobic moieties. Lipid analysis confirms the presence of lipids, and density in the primary binding pocket of afamin was modeled as palmitoleic acid, presenting the native O-acylation on serine 209 in human Wnt3a. The modeled complex between the experimental afamin structure and a Wnt3a homology model based on the XWnt8-Fz8-CRD fragment complex crystal structure is compelling, with favorable interactions comparable with the crystal structure complex. Afamin readily accommodates the conserved palmitoylated serine 209 of Wnt3a, providing a structural basis how afamin solubilizes hydrophobic and poorly soluble Wnt proteins.

Melanin-concentrating hormone receptor 1 (MCH-R1) is a G-protein-coupled receptor (GPCR) and a target for the development of therapeutics for obesity. The structure-based development of MCH-R1 and other GPCR antagonists is hampered by the lack of an available experimentally determined atomic structure. A ligand-steered homology modeling approach has been developed (where information about existing ligands is used explicitly to shape and optimize the binding site) followed by docking-based virtual screening. Top scoring compounds identified virtually were tested experimentally in an MCH-R1 competitive binding assay, and six novel chemotypes as low micromolar affinity antagonist "hits" were identified. This success rate is more than a 10-fold improvement over random high-throughput screening, which supports our ligand-steered method. Clearly, the ligand-steered homology modeling method reduces the uncertainty of structure modeling for difficult targets like GPCRs.

G-protein coupled receptors (GPCRs) are the largest family of cell-surface receptors involved in signal transmission. Drugs associated with GPCRs represent more than one fourth of the 100 top-selling drugs and are the targets of more than half of the current therapeutic agents on the market. Our methodology based on the internal coordinate mechanics (ICM) program can accurately identify the ligand-binding pocket in the currently available crystal structures of seven transmembrane (7TM) proteins [bacteriorhodopsin (BR) and bovine rhodopsin (bRho)]. The binding geometry of the ligand can be accurately predicted by ICM flexible docking with and without the loop regions, a useful finding for GPCR docking because the transmembrane regions are easier to model. We also demonstrate that the native ligand can be identified by flexible docking and scoring in 1.5% and 0.2% (for bRho and BR, respectively) of the best scoring compounds from two different types of compound database. The same procedure can be applied to the database of available chemicals to identify specific GPCR binders. Finally, we demonstrate that even if the sidechain positions in the bRho binding pocket are entirely wrong, their correct conformation can be fully restored with high accuracy (0.28 A) through the ICM global optimization with and without the ligand present. These binding site adjustments are critical for flexible docking of new ligands to known structures or for docking to GPCR homology models. The ICM docking method has the potential to be used to "de-orphanize" orphan GPCRs (oGPCRs) and to identify antagonists-agonists for GPCRs if an accurate model (experimentally and computationally validated) of the structure has been constructed or when future crystal structures are determined.

Circular dichroism (CD) spectroscopy is a valuable technique for the determination of protein secondary structures. Many linear and nonlinear algorithms have been developed for the empirical analysis of CD data, using reference databases derived from proteins of known structures. To date, the reference databases used by the various algorithms have all been derived from the spectra of soluble proteins. When applied to the analysis of soluble protein spectra, these methods generally produce calculated secondary structures that correspond well with crystallographic structures. In this study, however, it was shown that when applied to membrane protein spectra, the resulting calculations produce considerably poorer results. One source of this discrepancy may be the altered spectral peak positions (wavelength shifts) of membrane proteins due to the different dielectric of the membrane environment relative to that of water. These results have important consequences for studies that seek to use the existing soluble protein reference databases for the analyses of membrane proteins.

Understanding the molecular basis of ligand binding to receptors provides insights useful for rational drug design. This work describes development of a new antagonist radioligand of the type 1 cholecystokinin receptor (CCK1R), (2-fluorophenyl)-2,3-dihydro-3-[(3-isoquinolinylcarbonyl)amino]-6-methoxy-2-oxo-l-H-indole-3-propanoate (T-0632), and exploration of the molecular basis of its binding. This radioligand bound specifically with high affinity within an allosteric pocket of CCK1R. T-0632 fully inhibited binding and action of CCK at this receptor, while exhibiting no saturable binding to the closely related type 2 cholecystokinin receptor (CCK2R). Chimeric CCK1R/CCK2R constructs were used to explore the molecular basis of T-0632 binding. Exchanging exonic regions revealed the functional importance of CCK1R exon 3, extending from the bottom of transmembrane segment (TM) 3 to the top of TM5, including portions of the intramembranous pocket as well as the second extracellular loop region (ECL2). However, CCK1R mutants in which each residue facing the pocket was changed to that present in CCK2R had no negative impact on T-0632 binding. Extending the chimeric approach to ECL2 established the importance of its C-terminal region, and site-directed mutagenesis of each nonconserved residue in this region revealed the importance of Ser(208) at the top of TM5. A molecular model of T-0632-occupied CCK1R was consistent with these experimental determinants, also identifying Met(121) in TM3 and Arg(336) in TM6 as important. Although these residues are conserved in CCK2R, mutating them had a distinct impact on the two closely related receptors, suggesting differential orientation. This establishes the molecular basis of binding of a highly selective nonpeptidyl allosteric antagonist of CCK1R, illustrating differences in docking that extend beyond determinants attributable to distinct residues lining the intramembranous pocket in the two receptor subtypes.

A model of the endothelin G-protein-coupled receptor (ETA) has been constructed using a segmented approach. The model was produced using a bovine rhodopsin model as a template for the seven transmembrane α-helices. The three cytoplasmic loop regions and the C-terminal region were modeled on NMR structures of corresponding segments from bovine rhodopsin. The three extracellular loops were modeled on homologous loop regions in other proteins of known structure. The N-terminal region was modeled as a three-helix domain based on its homology with a hydrolase protein. To test the model, the FTDOCK algorithm was used to predict the ligand-binding site for the crystal structure of human endothelin. The site of docking is consistent with mutational and biochemical data. The principal sites of interaction in the endothelin ligand all lie on one face of a helix that has been implicated by structure-activity relationship studies as being essential for binding. As further support for the model, attempts to dock bigET, an inactive precursor to endothelin that does not bind to the receptor, found no sites for tight binding. The model of the receptor-ligand complex produced forms a basis for rational drug design of agonists and antagonists for this G-protein-coupled receptor.

While it is evident that the carboxyl-terminal region of natural peptide ligands bind to the amino-terminal domain of class B GPCRs, how their biologically critical amino-terminal regions dock to the receptor is unclear. We utilize cysteine trapping to systematically explore spatial approximations among residues in the first five positions of secretin and in every position within the receptor extracellular loops (ECLs). Only Cys2 and Cys5 secretin analogues exhibited full activity and retained moderate binding affinity (IC50: 92±4 and 83±1 nM, respectively). When these peptides probed 61 human secretin receptor cysteine-replacement mutants, a broad network of receptor residues could form disulfide bonds consistent with a dynamic ligand-receptor interface. Two distinct patterns of disulfide bond formation were observed: Cys2 predominantly labeled residues in the amino terminus of ECL2 and ECL3 (relative labeling intensity: Ser340, 94±7%; Pro341, 84±9%; Phe258, 73±5%; Trp274 62±8%), and Cys5 labeled those in the carboxyl terminus of ECL2 and ECL3 (Gln348, 100%; Ile347, 73±12%; Glu342, 59±10%; Phe351, 58±11%). These constraints were utilized in molecular modeling, providing improved understanding of the structure of the transmembrane bundle and interconnecting loops, the orientation between receptor domains, and the molecular basis of ligand docking. Key spatial approximations between peptide and receptor predicted by this model (H1-W274, D3-N268, G4-F258) were supported by mutagenesis and residue-residue complementation studies.—Dong, M., Xu, X., Ball, A. M., Makhoul, J. A., Lam, P. C.-H., Pinon, D. I., Orry, A., Sexton, P. M., Abagyan, R., Miller, L. J. Mapping spatial approximations between the amino terminus of secretin and each of the extracellular loops of its receptor using cysteine trapping.

SummaryThe histone acetyltransferase (HAT) p300/CBP is a transcriptional coactivator implicated in many gene regulatory pathways and protein acetylation events. Although p300 inhibitors have been reported, a potent, selective, and readily available active-site-directed small molecule inhibitor is not yet known. Here we use a structure-based, in silico screening approach to identify a commercially available pyrazolone-containing small molecule p300 HAT inhibitor, C646. C646 is a competitive p300 inhibitor with a Ki of 400 nM and is selective versus other acetyltransferases. Studies on site-directed p300 HAT mutants and synthetic modifications of C646 confirm the importance of predicted interactions in conferring potency. Inhibition of histone acetylation and cell growth by C646 in cells validate its utility as a pharmacologic probe and suggest that p300/CBP HAT is a worthy anticancer target.