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.

Allosteric modulation of receptors provides mechanistic safety while effectively achieving biologic endpoints otherwise difficult or impossible to obtain by other means. The theoretical case has been made for the development of a positive allosteric modulator (PAM) of the type 1 cholecystokinin receptor (CCK1R) having minimal intrinsic agonist activity to enhance meal-induced satiety for the treatment of obesity, while reducing the risk of side effects and/or toxicity. Unfortunately, such a drug does not currently exist. In this work, we have identified a PAM agonist of the CCK1R, SR146131, and determined its putative binding mode and receptor activation mechanism by combining molecular modeling, chimeric CCK1R/CCK2R constructs, and site-directed mutagenesis. We probed the structure-activity relationship of analogs of SR146131 for impact on agonism versus cooperativity of the analogs. This identified structural features that might be responsible for binding affinity and potency while retaining PAM activity. SR146131 and several of its analogs were docked into the receptor structure, which had the natural endogenous peptide agonist, cholecystokinin, already in the bound state (by docking), providing a refined structural model of the intact CCK1R holoreceptor. Both SR146131 and its analogs exhibited unique probe-dependent cooperativity with orthosteric peptide agonists and were simultaneously accommodated in this model, consistent with the derived structure-activity relationships. This provides improved understanding of the molecular basis for CCK1R-directed drug development.

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 type 1 cholecystokinin receptor (CCK1R) has multiple physiologic roles relating to nutrient homeostasis, including mediation of post-cibal satiety. This effect has been central in efforts to develop agonists of this receptor as part of a program to manage and/or prevent obesity. While a number of small molecule CCK1R agonists have been developed, none has yet been approved for clinical use, based on inadequate efficacy, side effects, or the potential for toxicity. Understanding the molecular details of docking and mechanism of action of these ligands can be helpful in the rational refinement and enhancement of small molecule drug candidates. In the current work, we have defined the mechanism of binding and activity of two triazolobenzodiazepinones, CE-326597 and PF-04756956, which are reported to be full agonist ligands. To achieve this, we utilized receptor binding with a series of allosteric and orthosteric radioligands at structurally-related CCK1R and CCK2R, as well as chimeric CCK1R/CCK2R constructs exchanging residues in the allosteric pocket, and assessment of biological activity. These triazolobenzodiazepinones docked within the intramembranous small molecule allosteric ligand pocket, with higher affinity binding to CCK2R than CCK1R, yet with biological activity exclusive to or greatly enhanced at CCK1R. These ligands exhibited cooperativity with benzodiazepine binding across the CCK1R homodimeric complex, resulting in their ability to inhibit only a fraction of the saturable binding of a benzodiazepine radioligand, unlike other small molecule antagonists and agonists of this receptor. This may contribute to the understanding of the unique short duration and reversible gallbladder contraction observed in vivo upon administration of these drugs.

Summary Mutations affecting the retinitis pigmentosa GTPase regulator-interacting protein 1 (RPGRIP1) interactome cause syndromic retinal dystrophies. RPGRIP1 interacts with the retinitis pigmentosa GTPase regulator (RPGR) through a domain homologous to RCC1 (RHD), a nucleotide exchange factor of Ran GTPase. However, functional relationships between RPGR and RPGRIP1 and their subcellular roles are lacking. We show by molecular modeling and analyses of RPGR disease-mutations that the RPGR-interacting domain (RID) of RPGRIP1 embraces multivalently the shared RHD of RPGR1–19 and RPGRORF15 isoforms and the mutations are non-overlapping with the interface found between RCC1 and Ran GTPase. RPGR disease-mutations grouped into six classes based on their structural locations and differential impairment with RPGRIP1 interaction. RPGRIP1α1 expression alone causes its profuse self-aggregation, an effect suppressed by co-expression of either RPGR isoform before and after RPGRIP1α1 self-aggregation ensue. RPGR1–19 localizes to the endoplasmic reticulum, whereas RPGRORF15 presents cytosolic distribution and they determine uniquely the subcellular co-localization of RPGRIP1α1. Disease mutations in RPGR1–19, RPGRORF15, or RID of RPGRIP1α1, singly or in combination, exert distinct effects on the subcellular targeting, co-localization or tethering of RPGRIP1α1 with RPGR1–19 or RPGRORF15 in kidney, photoreceptor and hepatocyte cell lines. Additionally, RPGRORF15, but not RPGR1–19, protects the RID of RPGRIP1α1 from limited proteolysis. These studies define RPGR- and cell-type-dependent targeting pathways with structural and functional plasticity modulating the expression of mutations in RPGR and RPGRIP1. Further, RPGR isoforms distinctively determine the subcellular targeting of RPGRIP1α1, with deficits in RPGRORF15-dependent intracellular localization of RPGRIP1α1 contributing to pathomechanisms shared by etiologically distinct syndromic retinal dystrophies.

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.

Ligand-docking-based methods are starting to play a critical role in lead discovery and optimization, thus resulting in new 'drug-candidates'. They offer the possibility to go beyond the pool of existing active compounds, and thus find novel chemotypes. A brief tutorial on ligand docking and structure-based virtual screening is presented highlighting current problems and limitations, together with the most recent methodological and algorithmic developments in the field. Recent successful applications of docking-based tools for hit discovery, lead optimization and target-biased library design are also presented. Special consideration is devoted to ongoing efforts to account for protein flexibility in structure-based virtual screening.

Computational ligand docking and screening is widely employed throughout the pharmaceutical industry to speed up the drug discovery process and identify drug candidates from very large pools of virtual compound libraries. When a ligand interacts with a receptor a number of structural changes within the ligand binding site might occur. It is therefore critical for these methods to accurately predict, or otherwise take into account the receptor flexibility upon ligand binding. This flexibility within the binding pocket explains why a diverse range of ligand sizes and shapes can sometimes bind to the same receptor pocket. This observation supersedes the notion that ligand-receptor interaction is a purely lock and key mechanism. The capability to correctly predict molecular interactions is critical for computer-aided molecular design technology. In this review, we discuss biological cases of receptor flexibility upon ligand binding that can range from large-scale movement of loops to single gate-keeper amino acid movements. In addition, we provide further evidence that rigid receptor docking alone will more than likely fail in the drug-discovery process. We then discuss computational methods, which have been developed to mimic flexibility within the binding pocket and predict ligand-receptor interactions. Early flexible receptor docking methods used soft-potential docking and rotamer libraries. More recently methods have focused on constructing an ensemble of structures generated by a variety of means including X-ray crystallography, NMR, Monte Carlo sampling, Normal Modes-based methods and Molecular Dynamics. It is evident that methods that ignore receptor flexibility can result in poorly docked solutions and therefore the challenge is to develop computational methods, which can accurately and efficiently predict this phenomenon.

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.

The development of large panels of simple-to-analyse genetic markers for tagging agronomically important genes and diversity studies in hexaploid bread wheat is an important goal in applied cereal genetic research. We have isolated and sequenced over 200 clones containing microsatellites from the wheat genome and have tested 153 primer pairs for genetic polymorphism using a panel of ten wheat varieties, including the parents of our main mapping cross. A subset comprising 49 primer pairs detects 76 loci, of which 74 can be unequivocably allocated to one of the wheat chromosomes. A relatively low frequency of the loci detected are from the D genome, and these loci show less polymorphism than those from the A and B genomes. Generally, the microsatellites show high levels of genetic polymorphism and an average of 3.5 alleles per locus with an average polymorphism information content (PIC), value of 0.51. The observed levels of polymorphism are positively correlated with the length of the microsatellite repeats. A high proportion, approximately two-thirds, of primer pairs designed to detect simple sequence repeat (SSR) variation in wheat do not generate the expected amplification products and, more significantly, often generate unresolvable PCR products. In general, our results agree closely with those obtained from other recent studies using microsatellites in plants.

The pleckstrin homology (PH) domain is a versatile fold that mediates a variety of protein–protein and protein–phosphatidylinositol lipid interactions. The Ran-binding protein 2 (RanBP2) contains four interspersed Ran GTPase-binding domains (RBD n = 1–4 ) with close structural homology to the PH domain of Bruton's tyrosine kinase. The RBD 2 , kinesin-binding domain (KBD) and RBD 3 comprise a tripartite domain (R 2 KR 3 ) of RanBP2 that causes the unfolding, microtubule binding and biphasic activation of kinesin-1, a crucial anterograde motor of mitochondrial motility. However, the interplay between Ran GTPase and R 2 KR 3 of RanBP2 in kinesin-1 activation and mitochondrial motility is elusive. We use structure–function, biochemical, kinetic and cell-based assays with time-lapse live-cell microscopy of over 260 000 mitochondrial-motility-related events to find mutually exclusive subdomains in RBD 2 and RBD 3 towards Ran GTPase binding, kinesin-1 activation and mitochondrial motility regulation. The RBD 2 and RBD 3 exhibit Ran-GTP-independent, subdomain and stereochemical-dependent discrimination on the biphasic kinetics of kinesin-1 activation or regulation of mitochondrial motility. Further, KBD alone and R 2 KR 3 stimulate and suppress, respectively, multiple biophysical parameters of mitochondrial motility. The regulation of the bidirectional transport of mitochondria by either KBD or R 2 KR 3 is highly coordinated, because their kinetic effects are accompanied always by changes in mitochondrial motile events of either transport polarity. These studies uncover novel roles in Ran GTPase-independent subdomains of RBD 2 and RBD 3 , and KBD of RanBP2, that confer antagonizing and multi-modal mechanisms of kinesin-1 activation and regulation of mitochondrial motility. These findings open new venues towards the pharmacological harnessing of cooperative and competitive mechanisms regulating kinesins, RanBP2 or mitochondrial motility in disparate human disorders.

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.

Despite remarkable advances in human genetics and other genetic model systems, the fruit fly, Drosophila melanogaster, remains a powerful experimental tool to probe with ease the inner workings of a myriad of biological and pathological processes, even when evolutionary forces impart apparent divergences to some of such processes. The understanding of such evolutionary differences provides mechanistic insights into genotype-phenotype correlations underpinning biological processes across metazoans. The pioneering work developed by the William Pak laboratory for the past four decades, and the work of others, epitomize the notion of how the Drosophila system breaks new fertile ground or complements research fields of high scientific and medical relevance. Among the three major genetic complementation groups produced by the Pak's laboratory and impairing distinct facets of photoreceptor neuronal function, the nina group (ninaA, …., ninaJ) selectively affects the biogenesis of G protein-coupled receptors (GPCRs), mediating the photoconversion and transduction of light stimuli. Among the nina genes identified, ninaA arguably assumes heightened significance for several reasons. First, it presents unique physiological selectivity toward the biogenesis of a subset of GPCRs, a standalone biological manifestation yet to be discerned for most mammalian homologues of NinaA. Second, NinaA belongs to a family of proteins, immunophilins, which are the primary targets for immunosuppressive drugs at the therapeutic forefront of a multitude of medical conditions. Third, NinaA closest homologue, cyclophilin B (CyPB/PPIB), is an immunophilin whose loss-of-function was found recently to cause osteogenesis imperfecta in the human. This report highlights advances made by studies on some members of immunophilins, the cyclophilins. Finally, it reexamines critically data and dogmas derived from past and recent genetic, structural, biological, and pathological studies on NinaA and few other cyclophilins that support some of such paradigms to be less than definite and advance our understanding of the roles of cyclophilins in cell function, disease, and therapeutic interventions.

The formation of ligand-protein complexes are critical for the correct functioning of a cell. The prediction of these interactions is important for our understanding of how the cell works and for the development of new drug molecules. Homology modeling is a method for predicting the structure of a protein based on a crystal structure template. Once a model of the protein is complete, a ligand-docking algorithm predicts the ligand-protein model interaction by searching for the best steric and energetically favorable fit. A refinement of the ligand-binding pocket improves the predicted interactions by considering the flexible nature of the ligand-binding pocket. In this chapter, we describe, from first principles, methods to identify and prepare the ligand-binding pocket in a protein model, to dock the ligand, and refine the resulting complex.

Synchrotron radiation circular dichroism (SRCD) spectroscopy is an emerging technique for sensitive determination of protein secondary structures and for monitoring of conformational changes. An important issue for its adoption as a useful technique is whether the high-intensity low-wavelength vacuum ultraviolet radiation in the SRCD chemically damages proteins. In this paper, using horse myoglobin as a test sample, it is shown that extensive irradiation in the SRCD does not produce any change in the chemical nature of the protein as detected by either SDS gel electrophoresis or mass spectrometry. In addition, no changes in the protein secondary structure are detectable from the SRCD spectra after extensive exposure to the SRCD beam.

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.

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.