Discovery of Protease-Activated Receptor 4 (PAR4)-Tethered Ligand Antagonists Using Ultralarge Virtual Screening.

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.

Rosetta's Predictive Ability for Low-Affinity Ligand Binding in Fragment-Based Drug Discovery.

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.

Structure, function and pharmacology of human itch GPCRs.

The MRGPRX family of receptors (MRGPRX1-4) is a family of mas-related G-protein-coupled receptors that have evolved relatively recently1. Of these, MRGPRX2 and MRGPRX4 are key physiological and pathological mediators of itch and related mast cell-mediated hypersensitivity reactions2-5. MRGPRX2 couples to both Gi and Gq in mast cells6. Here we describe agonist-stabilized structures of MRGPRX2 coupled to Gi1 and Gq 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 Gq. Together, these findings should accelerate the structure-guided discovery of therapeutic agents for pain, itch and mast cell-mediated hypersensitivity.

Molecular basis for the evolved instability of a human G-protein coupled receptor.

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.

Computational redesign of a fluorogen activating protein with Rosetta.

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.

Viewpoints on the First Transatlantic GPCR Symposium for Early-Career Investigators.

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.

A practical guide to large-scale docking.

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.

Modeling Immunity with Rosetta: Methods for Antibody and Antigen Design.

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.

Macromolecular modeling and design in Rosetta: recent methods and frameworks.

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.

Comparative modeling and docking of chemokine-receptor interactions with Rosetta.

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.

Modeling the complete chemokine-receptor interaction.

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.

Structural basis of ligand binding modes at the neuropeptide Y Y1 receptor.

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 Y1, Y2, Y4 and Y5 receptors, with different affinity and selectivity 3 . NPY is the most powerful stimulant of food intake and this effect is primarily mediated by the Y1 receptor (Y1R) 4 . A number of peptides and small-molecule compounds have been characterized as Y1R 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 Y1R 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 Y1R to several structurally diverse antagonists and the determinants of ligand selectivity. The Y1R 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 Y1R can enable structure-based drug discovery that targets NPY receptors.

Improved in Vitro Folding of the Y2 G Protein-Coupled Receptor into Bicelles.

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 (Y2R) 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 Y2R 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.

Protocols for Molecular Modeling with Rosetta3 and RosettaScripts.

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.

Rosetta and the Design of Ligand Binding Sites.

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.

Human cytomegalovirus UL97 phosphorylates the viral nuclear egress complex.

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.

Dynamic and nucleolin-dependent localization of human cytomegalovirus UL84 to the periphery of viral replication compartments and nucleoli.

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.

A mutation deleting sequences encoding the amino terminus of human cytomegalovirus UL84 impairs interaction with UL44 and capsid localization.

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.