Search

Here, we demonstrate a structure-based small molecule virtual screening and lead optimization pipeline using a homology model of a difficult-to-drug G-protein-coupled receptor (GPCR) target. Protease-activated receptor 4 (PAR4) is activated by thrombin cleavage, revealing a tethered ligand that activates the receptor, making PAR4 a challenging target. A virtual screen of a make-on-demand chemical library yielded a one-hit compound. From the single-hit compound, we developed a novel series of PAR4 antagonists. Subsequent lead optimization via simultaneous virtual library searches and structure-based rational design efforts led to potent antagonists of thrombin-induced activation. Interestingly, this series of antagonists was active against PAR4 activation by the native protease thrombin cleavage but not the synthetic PAR4 agonist peptide AYPGKF., Competing Interests: The authors declare no competing financial interest., (© 2024 The Authors. Published by American Chemical Society.)

Fragment-based drug discovery begins with the identification of small molecules with a molecular weight of usually less than 250 Da which weakly bind to the protein of interest. This technique is challenging for computational docking methods as binding is determined by only a few specific interactions. Inaccuracies in the energy function or slight deviations in the docking pose can lead to the prediction of incorrect binding or difficulties in ranking fragments in in silico screening. Here, we test RosettaLigand by docking a series of fragments to a cysteine-depleted variant of the TIM-barrel protein, HisF (UniProtKB Q9X0C6). We compare the computational results with experimental NMR spectroscopy screens. NMR spectroscopy gives details on binding affinities of individual ligands, which allows assessment of the ligand-ranking ability using RosettaLigand and also provides feedback on the location of the binding pocket, which serves as a reliable test of RosettaLigand's ability to identify plausible binding poses. From a library screen of 3456 fragments, we identified a set of 31 ligands with intrinsic affinities to HisF with dissociation constants as low as 400 μM. The same library of fragments was blindly screened in silico. RosettaLigand was able to rank binders before non-binders with an area under the curve of the receiver operating characteristics of 0.74. The docking poses observed for binders agreed with the binding pocket identified by NMR chemical shift perturbations for all fragments. Taken together, these results provide a baseline performance of RosettaLigand in a fragment-based drug discovery setting.

G protein-coupled receptors (GPCRs) represent the largest membrane protein family and a significant target class for therapeutics. Receptors from GPCRs' largest class, class A, influence virtually every aspect of human physiology. About 45% of the members of this family endogenously bind flexible peptides or peptides segments within larger protein ligands. While many of these peptides have been structurally characterized in their solution state, the few studies of peptides in their receptor-bound state suggest that these peptides interact with a shared set of residues and undergo significant conformational changes. For the purpose of understanding binding dynamics and the development of peptidomimetic drug compounds, further studies should investigate the peptide ligands that are complexed to their cognate receptor.

The MRGPRX family of receptors (MRGPRX1-4) is a family of mas-related G-protein-coupled receptors that have evolved relatively recently 1 . Of these, MRGPRX2 and MRGPRX4 are key physiological and pathological mediators of itch and related mast cell-mediated hypersensitivity reactions 2-5 . MRGPRX2 couples to both G i and G q in mast cells 6 . Here we describe agonist-stabilized structures of MRGPRX2 coupled to G i1 and G q in ternary complexes with the endogenous peptide cortistatin-14 and with a synthetic agonist probe, respectively, and the development of potent antagonist probes for MRGPRX2. We also describe a specific MRGPRX4 agonist and the structure of this agonist in a complex with MRGPRX4 and G q . Together, these findings should accelerate the structure-guided discovery of therapeutic agents for pain, itch and mast cell-mediated hypersensitivity., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

Membrane proteins are prone to misfolding and degradation. This is particularly true for mammalian forms of the gonadotropin-releasing hormone receptor (GnRHR). Although they function at the plasma membrane, mammalian GnRHRs accumulate within the secretory pathway. Their apparent instability is believed to have evolved through selection for attenuated GnRHR activity. Nevertheless, the molecular basis of this adaptation remains unclear. We show that adaptation coincides with a C-terminal truncation that compromises the translocon-mediated membrane integration of its seventh transmembrane domain (TM7). We also identify a series of polar residues in mammalian GnRHRs that compromise the membrane integration of TM2 and TM6. Reverting a lipid-exposed polar residue in TM6 to an ancestral hydrophobic residue restores expression with no impact on function. Evolutionary trends suggest variations in the polarity of this residue track with reproductive phenotypes. Our findings suggest that the marginal energetics of cotranslational folding can be exploited to tune membrane protein fitness., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

The use of unnatural fluorogenic molecules widely expands the pallet of available genetically encoded fluorescent imaging tools through the design of fluorogen activating proteins (FAPs). While there is already a handful of such probes available, each of them went through laborious cycles of in vitro screening and selection. Computational modeling approaches are evolving incredibly fast right now and are demonstrating great results in many applications, including de novo protein design. It suggests that the easier task of fine-tuning the fluorogen-binding properties of an already functional protein in silico should be readily achievable. To test this hypothesis, we used Rosetta for computational ligand docking followed by protein binding pocket redesign to further improve the previously described FAP DiB1 that is capable of binding to a BODIPY-like dye M739. Despite an inaccurate initial docking of the chromophore, the incorporated mutations nevertheless improved multiple photophysical parameters as well as the overall performance of the tag. The designed protein, DiB-RM, shows higher brightness, localization precision, and apparent photostability in protein-PAINT super-resolution imaging compared to its parental variant DiB1. Moreover, DiB-RM can be cleaved to obtain an efficient split system with enhanced performance compared to a parental DiB-split system. The possible reasons for the inaccurate ligand binding pose prediction and its consequence on the outcome of the design experiment are further discussed., Competing Interests: The authors have declared that no competing interests exist.

Structure-based docking screens of large compound libraries have become common in early drug and probe discovery. As computer efficiency has improved and compound libraries have grown, the ability to screen hundreds of millions, and even billions, of compounds has become feasible for modest-sized computer clusters. This allows the rapid and cost-effective exploration and categorization of vast chemical space into a subset enriched with potential hits for a given target. To accomplish this goal at speed, approximations are used that result in undersampling of possible configurations and inaccurate predictions of absolute binding energies. Accordingly, it is important to establish controls, as are common in other fields, to enhance the likelihood of success in spite of these challenges. Here we outline best practices and control docking calculations that help evaluate docking parameters for a given target prior to undertaking a large-scale prospective screen, with exemplification in one particular target, the melatonin receptor, where following this procedure led to direct docking hits with activities in the subnanomolar range. Additional controls are suggested to ensure specific activity for experimentally validated hit compounds. These guidelines should be useful regardless of the docking software used. Docking software described in the outlined protocol (DOCK3.7) is made freely available for academic research to explore new hits for a range of targets., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

In July 2021, we organized a virtual symposium aimed at early-career investigators (ECIs) in G protein-coupled receptor (GPCR) research: the first Transatlantic ECI GPCR Symposium. Here, we discuss the proceedings of this symposium and the unique networking events with GPCR leaders including the Nobel Laureates Dr. Robert Lefkowitz and Dr. Brian Kobilka., Competing Interests: The authors declare no competing financial interest., (© 2021 American Chemical Society.)

Structure-based antibody and antigen design has advanced greatly in recent years, due not only to the increasing availability of experimentally determined structures but also to improved computational methods for both prediction and design. Constant improvements in performance within the Rosetta software suite for biomolecular modeling have given rise to a greater breadth of structure prediction, including docking and design application cases for antibody and antigen modeling. Here, we present an overview of current protocols for antibody and antigen modeling using Rosetta and exemplify those by detailed tutorials originally developed for a Rosetta workshop at Vanderbilt University. These tutorials cover antibody structure prediction, docking, and design and antigen design strategies, including the addition of glycans in Rosetta. We expect that these materials will allow novice users to apply Rosetta in their own projects for modeling antibodies and antigens.

As sequencing methodologies continue to advance, the availability of protein sequences far outpaces the ability of structure determination. Homology modeling is used to bridge this gap but relies on high-identity templates for accurate model building. G-protein coupled receptors (GPCRs) represent a significant target class for pharmaceutical therapies in which homology modeling could fill the knowledge gap for structure-based drug design. To date, only about 17% of druggable GPCRs have had their structures characterized at atomic resolution. However, modeling of the remaining 83% is hindered by the low sequence identity between receptors. Here we test key inputs in the model building process using GPCRs as a focus to improve the pipeline in two critical ways: Firstly, we use a blended sequence- and structure-based alignment that accounts for structure conservation in loop regions. Secondly, by merging multiple template structures into one comparative model, the best possible template for every region of a target can be used expanding the conformational space sampled in a meaningful way. This optimization allows for accurate modeling of receptors using templates as low as 20% sequence identity, which accounts for nearly the entire druggable space of GPCRs. A model database of all non-odorant GPCRs is made available at www.rosettagpcr.org. Additionally, all protocols are made available with insights into modifications that may improve accuracy at new targets., Competing Interests: The authors have declared that no competing interests exist.

Chemokine receptors are a subset of G protein-coupled receptors defined by the distinct property of binding small protein ligands in the chemokine family. Chemokine receptors recognize their ligands by a mechanism that is distinct from other class A GPCRs that bind peptides or small molecules. For this reason, structural information on other ligand-GPCR interactions are only indirectly relevant to understanding the chemokine receptor interface. Additionally, the experimentally determined structures of chemokine-GPCR complexes represent less than 3% of the known interactions of this complex, multi-ligand/multi-receptor network. To enable predictive modeling of the remaining 97% of interactions, a general in silico protocol was designed to utilize existing chemokine receptor crystal structures, co-crystal structures, and NMR ensembles of chemokines bound to receptor fragments. This protocol was benchmarked on the ability to predict each of the three published co-crystal structures, while being blinded to the target structure. Averaging ensembles selected from the top-ranking models reproduced up to 84% of the intermolecular contacts found in the crystal structure, with the lowest Cα-RMSD of the complex at 3.3 Å. The chemokine receptor N-terminus, unresolved in crystal structures, was included in the modeling and recapitulates contacts with known sulfotyrosine binding pockets seen in structures derived from experimental NMR data. This benchmarking experiment suggests that realistic homology models of chemokine-GPCR complexes can be generated by leveraging current structural data., (Copyright © 2019 Elsevier Inc. All rights reserved.)

The Rosetta software for macromolecular modeling, docking and design is extensively used in laboratories worldwide. During two decades of development by a community of laboratories at more than 60 institutions, Rosetta has been continuously refactored and extended. Its advantages are its performance and interoperability between broad modeling capabilities. Here we review tools developed in the last 5 years, including over 80 methods. We discuss improvements to the score function, user interfaces and usability. Rosetta is available at http://www.rosettacommons.org.

The peptide ghrelin targets the growth hormone secretagogue receptor 1a (GHSR) to signal changes in cell metabolism and is a sought-after therapeutic target, although no structure is known to date. To investigate the structural basis of ghrelin binding to GHSR, we used solid-state nuclear magnetic resonance (NMR) spectroscopy, site-directed mutagenesis, and Rosetta modeling. The use of saturation transfer difference NMR identified key residues in the peptide for receptor binding beyond the known motif. This information combined with assignment of the secondary structure of ghrelin in its receptor-bound state was incorporated into Rosetta using an approach that accounts for flexible binding partners. The NMR data and models revealed an extended binding surface that was confirmed via mutagenesis. Our results agree with a growing evidence of peptides interacting via two sites at G protein-coupled receptors., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Chemokines are soluble, secreted proteins that induce chemotaxis of leukocytes and other cells. Migratory cells can sense the chemokine concentration gradient following chemokine binding and activation of chemokine receptors, a subset of the G protein-coupled receptor (GPCR) superfamily. Chemokine receptor signaling plays a central role in cell migration during inflammatory responses as well as in cancer and other diseases. Given their important role in mediating essential pathologic and physiologic processes, chemokines and their receptors are attractive targets for therapeutic development. A better understanding of the molecular basis of chemokine-GPCR interactions will aid in the understanding of the mechanistic basis for chemokine function in disease-related processes, as well as aid in the design of new therapeutics. High resolution protein structures are critical for determining these mechanisms and investigating the interactions between approximately 50 chemokines and 20 chemokine receptors. Currently, three unique structures of chemokine-GPCR complexes have been determined and have greatly broadened our knowledge of this large protein-protein interaction. While these structures represent only a small fraction of clinically relevant chemokines and receptors, they can be exploited as scaffolds for homology modeling to understand the chemokine-GPCR interactions. This chapter presents a specialized methodology to construct and validate models of chemokine-GPCR complexes using the Rosetta software suite., (© 2019 Elsevier Inc. All rights reserved.)

The poor norepinephrine innervation and high density of Gi/o-coupled α 2A - and α 2C -adrenoceptors in the striatum and the dense striatal dopamine innervation have prompted the possibility that dopamine could be an effective adrenoceptor ligand. Nevertheless, the reported adrenoceptor agonistic properties of dopamine are still inconclusive. In this study, we analyzed the binding of norepinephrine, dopamine, and several compounds reported as selective dopamine D 2 -like receptor ligands, such as the D 3 receptor agonist 7-OH-PIPAT and the D 4 receptor agonist RO-105824, to α 2 -adrenoceptors in cortical and striatal tissue, which express α 2A -adrenoceptors and both α 2A - and α 2C -adrenoceptors, respectively. The affinity of dopamine for α 2 -adrenoceptors was found to be similar to that for D 1 -like and D 2 -like receptors. Moreover, the exogenous dopamine receptor ligands also showed high affinity for α 2A - and α 2C -adrenoceptors. Their ability to activate Gi/o proteins through α 2A - and α 2C -adrenoceptors was also analyzed in transfected cells with bioluminescent resonance energy transfer techniques. The relative ligand potencies and efficacies were dependent on the Gi/o protein subtype. Furthermore, dopamine binding to α 2 -adrenoceptors was functional, inducing changes in dynamic mass redistribution, adenylyl cyclase activity, and ERK1/2 phosphorylation. Binding events were further studied with computer modeling of ligand docking. Docking of dopamine at α 2A - and α 2C -adrenoceptors was nearly identical to its binding to the crystallized D 3 receptor. Therefore, we provide conclusive evidence that α 2A - and α 2C -adrenoceptors are functional receptors for norepinephrine, dopamine, and other previously assumed selective D 2 -like receptor ligands, which calls for revisiting previous studies with those ligands.

Neuropeptide Y (NPY) receptors belong to the G-protein-coupled receptor superfamily and have important roles in food intake, anxiety and cancer biology 1,2 . The NPY-Y receptor system has emerged as one of the most complex networks with three peptide ligands (NPY, peptide YY and pancreatic polypeptide) binding to four receptors in most mammals, namely the Y 1 , Y 2 , Y 4 and Y 5 receptors, with different affinity and selectivity 3 . NPY is the most powerful stimulant of food intake and this effect is primarily mediated by the Y 1 receptor (Y 1 R) 4 . A number of peptides and small-molecule compounds have been characterized as Y 1 R antagonists and have shown clinical potential in the treatment of obesity 4 , tumour 1 and bone loss 5 . However, their clinical usage has been hampered by low potency and selectivity, poor brain penetration ability or lack of oral bioavailability 6 . Here we report crystal structures of the human Y 1 R bound to the two selective antagonists UR-MK299 and BMS-193885 at 2.7 and 3.0 Å resolution, respectively. The structures combined with mutagenesis studies reveal the binding modes of Y 1 R to several structurally diverse antagonists and the determinants of ligand selectivity. The Y 1 R structure and molecular docking of the endogenous agonist NPY, together with nuclear magnetic resonance, photo-crosslinking and functional studies, provide insights into the binding behaviour of the agonist and for the first time, to our knowledge, determine the interaction of its N terminus with the receptor. These insights into Y 1 R can enable structure-based drug discovery that targets NPY receptors.

Prerequisite for structural studies on G protein-coupled receptors is the preparation of highly concentrated, stable, and biologically active receptor samples in milligram amounts of protein. Here, we present an improved protocol for Escherichia coli expression, functional refolding, and reconstitution into bicelles of the human neuropeptide Y receptor type 2 (Y 2 R) for solution and solid-state NMR experiments. The isotopically labeled receptor is expressed in inclusion bodies and purified using SDS. We studied the details of an improved preparation protocol including the in vitro folding of the receptor, e.g., the native disulfide bridge formation, the exchange of the denaturating detergent SDS, and the functional reconstitution into bicelle environments of varying size. Full pharmacological functionality of the Y 2 R preparation was shown by a ligand affinity of 4 nM and G-protein activation. Further, simple NMR experiments are used to test sample quality in high micromolar concentration.

Previously, we published an article providing an overview of the Rosetta suite of biomacromolecular modeling software and a series of step-by-step tutorials [Kaufmann, K. W., et al. (2010) Biochemistry 49, 2987-2998]. The overwhelming positive response to this publication we received motivates us to here share the next iteration of these tutorials that feature de novo folding, comparative modeling, loop construction, protein docking, small molecule docking, and protein design. This updated and expanded set of tutorials is needed, as since 2010 Rosetta has been fully redesigned into an object-oriented protein modeling program Rosetta3. Notable improvements include a substantially improved energy function, an XML-like language termed "RosettaScripts" for flexibly specifying modeling task, new analysis tools, the addition of the TopologyBroker to control conformational sampling, and support for multiple templates in comparative modeling. Rosetta's ability to model systems with symmetric proteins, membrane proteins, noncanonical amino acids, and RNA has also been greatly expanded and improved.

Proteins that bind small molecules (ligands) can be used as biosensors, signal modulators, and sequestering agents. When naturally occurring proteins for a particular target ligand are not available, artificial proteins can be computationally designed. We present a protocol based on RosettaLigand to redesign an existing protein pocket to bind a target ligand. Starting with a protein structure and the structure of the ligand, Rosetta can optimize both the placement of the ligand in the pocket and the identity and conformation of the surrounding sidechains, yielding proteins that bind the target compound.

Unlabelled: Herpesvirus nucleocapsids exit the host cell nucleus in an unusual process known as nuclear egress. The human cytomegalovirus (HCMV) UL97 protein kinase is required for efficient nuclear egress, which can be explained by its phosphorylation of the nuclear lamina component lamin A/C, which disrupts the nuclear lamina. We found that a dominant negative lamin A/C mutant complemented the replication defect of a virus lacking UL97 in dividing cells, validating this explanation. However, as complementation was incomplete, we investigated whether the HCMV nuclear egress complex (NEC) subunits UL50 and UL53, which are required for nuclear egress and recruit UL97 to the nuclear rim, are UL97 substrates. Using mass spectrometry, we detected UL97-dependent phosphorylation of UL50 residue S216 (UL50-S216) and UL53-S19 in infected cells. Moreover, UL53-S19 was specifically phosphorylated by UL97 in vitro. Notably, treatment of infected cells with the UL97 inhibitor maribavir or infection with a UL97 mutant led to a punctate rather than a continuous distribution of the NEC at the nuclear rim. Alanine substitutions in both UL50-S216 and UL53-S19 resulted in a punctate distribution of the NEC in infected cells and also decreased virus production and nuclear egress in the absence of maribavir. These results indicate that UL97 phosphorylates the NEC and suggest that this phosphorylation modulates nuclear egress. Thus, the UL97-NEC interaction appears to recruit UL97 to the nuclear rim both for disruption of the nuclear lamina and phosphorylation of the NEC., Importance: Human cytomegalovirus (HCMV) causes birth defects and it can cause life-threatening diseases in immunocompromised patients. HCMV assembles in the nucleus and then translocates to the cytoplasm in an unusual process termed nuclear egress, an attractive target for antiviral therapy. A viral enzyme, UL97, is important for nuclear egress. It has been proposed that this is due to its role in disruption of the nuclear lamina, which would otherwise impede nuclear egress. In validating this proposal, we showed that independent disruption of the lamina can overcome a loss of UL97, but only partly, suggesting additional roles for UL97 during nuclear egress. We then found that UL97 phosphorylates the viral nuclear egress complex (NEC), which is essential for nuclear egress, and we obtained evidence that this phosphorylation modulates this process. Our results highlight a new role for UL97, the mutual dependence of the viral NEC and UL97 during nuclear egress, and differences among herpesviruses., (Copyright © 2015, American Society for Microbiology. All Rights Reserved.)

Protein-protein and protein-nucleic acid interactions within subcellular compartments are required for viral genome replication. To understand the localization of the human cytomegalovirus viral replication factor UL84 relative to other proteins involved in viral DNA synthesis and to replicating viral DNA in infected cells, we created a recombinant virus expressing a FLAG-tagged version of UL84 (UL84FLAG) and used this virus in immunofluorescence assays. UL84FLAG localization differed at early and late times of infection, transitioning from diffuse distribution throughout the nucleus to exclusion from the interior of replication compartments, with some concentration at the periphery of replication compartments with newly labeled DNA and the viral DNA polymerase subunit UL44. Early in infection, UL84FLAG colocalized with the viral single-stranded DNA binding protein UL57, but colocalization became less prominent as infection progressed. A portion of UL84FLAG also colocalized with the host nucleolar protein nucleolin at the peripheries of both replication compartments and nucleoli. Small interfering RNA (siRNA)-mediated knockdown of nucleolin resulted in a dramatic elimination of UL84FLAG from replication compartments and other parts of the nucleus and its accumulation in the cytoplasm. Reciprocal coimmunoprecipitation of viral proteins from infected cell lysates revealed association of UL84, UL44, and nucleolin. These results indicate that UL84 localization during infection is dynamic, which is likely relevant to its functions, and suggest that its nuclear and subnuclear localization is highly dependent on direct or indirect interactions with nucleolin. Importance: The protein-protein interactions among viral and cellular proteins required for replication of the human cytomegalovirus (HCMV) DNA genome are poorly understood. We sought to understand how an enigmatic HCMV protein critical for virus replication, UL84, localizes relative to other viral and cellular proteins required for HCMV genome replication and replicating viral DNA. We found that UL84 localizes with viral proteins, viral DNA, and the cellular nucleolar protein nucleolin in the subnuclear replication compartments in which viral DNA replication occurs. Unexpectedly, we also found localization of UL84 with nucleolin in nucleoli and showed that the presence of nucleolin is involved in localization of UL84 to the nucleus. These results add to previous work showing the importance of nucleolin in replication compartment architecture and viral DNA synthesis and are relevant to understanding UL84 function., (Copyright © 2014, American Society for Microbiology. All Rights Reserved.)

Protein-protein interactions are required for many biological functions. Previous work has demonstrated an interaction between the human cytomegalovirus DNA polymerase subunit UL44 and the viral replication factor UL84. In this study, glutathione S-transferase pulldown assays indicated that residues 1 to 68 of UL84 are both necessary and sufficient for efficient interaction of UL84 with UL44 in vitro. We created a mutant virus in which sequences encoding these residues were deleted. This mutant displayed decreased virus replication compared to wild-type virus. Immunoprecipitation assays showed that the mutation decreased but did not abrogate association of UL84 with UL44 in infected cell lysate, suggesting that the association in the infected cell can involve other protein-protein interactions. Further immunoprecipitation assays indicated that IRS1, TRS1, and nucleolin are candidates for such interactions in infected cells. Quantitative real-time PCR analysis of viral DNA indicated that the absence of the UL84 amino terminus does not notably affect viral DNA synthesis. Western blotting experiments and pulse labeling of infected cells with [(35)S]methionine demonstrated a rather modest downregulation of levels of multiple proteins and particularly decreased levels of the minor capsid protein UL85. Electron microscopy demonstrated that viral capsids assemble but are mislocalized in nuclei of cells infected with the mutant virus, with fewer cytoplasmic capsids detected. In sum, deletion of the sequences encoding the amino terminus of UL84 affects interaction with UL44 and virus replication unexpectedly, not viral DNA synthesis. Mislocalization of viral capsids in infected cell nuclei likely contributes to the observed decrease in virus replication.

Cholestatic jaundice and rapidly deteriorating hepatic function developed in a 19-year-old man with a lifelong history of photosensitivity. Laboratory studies revealed the characteristic increased erythrocyte and fecal protoporphyrin levels of erythropoietic protoporphyria. Progressive hepatic failure was treated by orthotopic liver transplantation six months after the first clinical indication of hepatic dysfunction. Characteristic light microscopic, fluorescence microscopic, and electron microscopic findings of erythropoietic protoporphyria were present in skin and liver. Four weeks after liver transplantation, the patient died of disseminated candidiasis. At autopsy, the donor liver had no microscopic evidence of protoporphyrin accumulation, although tissue protoporphyrin levels were mildly elevated.