Harnessing the power of an X-ray laser for serial crystallography of membrane proteins crystallized in lipidic cubic phase

G Protein-Coupled Receptor-Ligand Pose and Functional Class Prediction

G protein-coupled receptor (GPCR) transmembrane protein family members play essential roles in physiology. Numerous pharmaceuticals target GPCRs, and many drug discovery programs utilize virtual screening (VS) against GPCR targets. Improvements in the accuracy of predicting new molecules that bind to and either activate or inhibit GPCR function would accelerate such drug discovery programs. This work addresses two significant research questions. First, do ligand interaction fingerprints provide a substantial advantage over automated methods of binding site selection for classical docking? Second, can the functional status of prospective screening candidates be predicted from ligand interaction fingerprints using a random forest classifier? Ligand interaction fingerprints were found to offer modest advantages in sampling accurate poses, but no substantial advantage in the final set of top-ranked poses after scoring, and, thus, were not used in the generation of the ligand-receptor complexes used to train and test the random forest classifier. A binary classifier which treated agonists, antagonists, and inverse agonists as active and all other ligands as inactive proved highly effective in ligand function prediction in an external test set of GPR31 and TAAR2 candidate ligands with a hit rate of 82.6% actual actives within the set of predicted actives.

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

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

Systematic analyses of the sequence conservation and ligand interaction patterns of purinergic P1 and P2Y receptors provide a structural basis for receptor selectivity

Purinergic receptors are membrane proteins that regulate numerous cellular functions by catalyzing reactions involving purine nucleotides or nucleosides. Among the three receptor families, i.e., P1, P2X, and P2Y, the P1 and P2Y receptors share common structural features of class A GPCR. Comprehensive sequence and structural analysis revealed that the P1 and P2Y receptors belong to two distinct groups. They exhibit different ligand-binding site features that can distinguish between specific activators. These specific amino acid residues in the binding cavity may be involved in the selectivity and unique pharmacological behavior of each subtype. In this study, we conducted a structure-based analysis of purinergic P1 and P2Y receptors to identify their evolutionary signature and obtain structural insights into ligand recognition and selectivity. The structural features of the P1 and P2Y receptor classes were compared based on sequence conservation and ligand interaction patterns. Orthologous protein sequences were collected for the P1 and P2Y receptors, and sequence conservation was calculated based on Shannon entropy to identify highly conserved residues. To analyze the ligand interaction patterns, we performed docking studies on the P1 and P2Y receptors using known ligand information extracted from the ChEMBL database. We analyzed how the conserved residues are related to ligand-binding sites and how the key interacting residues differ between P1 and P2Y receptors, or between agonists and antagonists. We extracted new similarities and differences between the receptor subtypes, and the results can be used for designing new ligands by predicting hotspot residues that are important for functional selectivity.

Electrically stimulated droplet injector for reduced sample consumption in serial crystallography

With advances in X-ray free-electron lasers (XFELs), serial femtosecond crystallography (SFX) has enabled the static and dynamic structure determination for challenging proteins such as membrane protein complexes. In SFX with XFELs, the crystals are typically destroyed after interacting with a single XFEL pulse. Therefore, thousands of new crystals must be sequentially introduced into the X-ray beam to collect full data sets. Because of the serial nature of any SFX experiment, up to 99% of the sample delivered to the X-ray beam during its "off-time" between X-ray pulses is wasted due to the intrinsic pulsed nature of all current XFELs. To solve this major problem of large and often limiting sample consumption, we report on improvements of a revolutionary sample-saving method that is compatible with all current XFELs. We previously reported 3D-printed injection devices coupled with gas dynamic virtual nozzles (GDVNs) capable of generating samples containing droplets segmented by an immiscible oil phase for jetting crystal-laden droplets into the path of an XFEL. Here, we have further improved the device design by including metal electrodes inducing electrowetting effects for improved control over droplet generation frequency to stimulate the droplet release to matching the XFEL repetition rate by employing an electrical feedback mechanism. We report the improvements in this electrically triggered segmented flow approach for sample conservation in comparison with a continuous GDVN injection using the microcrystals of lysozyme and 3-deoxy-D-manno-octulosonate 8-phosphate synthase and report the segmented flow approach for sample injection applied at the Macromolecular Femtosecond Crystallography instrument at the Linear Coherent Light Source for the first time.

GPCR Agonist-to-Antagonist Conversion: Enabling the Design of Nucleoside Functional Switches for the A2A Adenosine Receptor

Modulators of the G protein-coupled A2A adenosine receptor (A2AAR) have been considered promising agents to treat Parkinson's disease, inflammation, cancer, and central nervous system disorders. Herein, we demonstrate that a thiophene modification at the C8 position in the common adenine scaffold converted an A2AAR agonist into an antagonist. We synthesized and characterized a novel A2AAR antagonist, 2 (LJ-4517), with Ki = 18.3 nM. X-ray crystallographic structures of 2 in complex with two thermostabilized A2AAR constructs were solved at 2.05 and 2.80 Å resolutions. In contrast to A2AAR agonists, which simultaneously interact with both Ser2777.42 and His2787.43, 2 only transiently contacts His2787.43, which can be direct or water-mediated. The n-hexynyl group of 2 extends into an A2AAR exosite. Structural analysis revealed that the introduced thiophene modification restricted receptor conformational rearrangements required for subsequent activation. This approach can expand the repertoire of adenosine receptor antagonists that can be designed based on available agonist scaffolds.

Stability, flow alignment and a phase transition of the lipidic cubic phase during continuous flow injection

Continuous flow injection is a key technology for serial crystallography measurements of protein crystals suspended in the lipidic cubic phase (LCP). To date, there has been little discussion in the literature regarding the impact of the injection process itself on the structure of the lipidic phase. This is despite the fact that the phase of the injection matrix is critical for the flow properties of the stream and potentially for sample stability. Here we report small-angle X-ray scattering measurements of a monoolein:water mixture during continuous delivery using a high viscosity injector. We observe both an alignment and modification of the LCP as a direct result of the injection process. The orientation of the cubic lattice with respect to the beam was estimated based on the anisotropy of the diffraction pattern and does not correspond to a single low order zone axis. The solvent fraction was also observed to impact the stability of the cubic phase during injection. In addition, depending on the distance traveled by the lipid after exiting the needle, the phase is observed to transition from a pure diamond phase (Pn3m) to a mixture containing both gyriod (Ia3d) and lamellar (L_α) phases. Finite element modelling of the observed phase behaviour during injection indicates that the pressure exerted on the lipid stream during extrusion accounts for the variations in the phase composition of the monoolein:water mixture.

Molecular Simulations and Drug Discovery of Adenosine Receptors

G protein-coupled receptors (GPCRs) represent the largest family of human membrane proteins. Four subtypes of adenosine receptors (ARs), the A1AR, A2AAR, A2BAR and A3AR, each with a unique pharmacological profile and distribution within the tissues in the human body, mediate many physiological functions and serve as critical drug targets for treating numerous human diseases including cancer, neuropathic pain, cardiac ischemia, stroke and diabetes. The A1AR and A3AR preferentially couple to the Gi/o proteins, while the A2AAR and A2BAR prefer coupling to the Gs proteins. Adenosine receptors were the first subclass of GPCRs that had experimental structures determined in complex with distinct G proteins. Here, we will review recent studies in molecular simulations and computer-aided drug discovery of the adenosine receptors and also highlight their future research opportunities.

Pink-beam serial femtosecond crystallography for accurate structure-factor determination at an X-ray free-electron laser

Serial femtosecond crystallography (SFX) at X-ray free-electron lasers (XFELs) enables essentially radiation-damage-free macromolecular structure determination using microcrystals that are too small for synchrotron studies. However, SFX experiments often require large amounts of sample in order to collect highly redundant data where some of the many stochastic errors can be averaged out to determine accurate structure-factor amplitudes. In this work, the capability of the Swiss X-ray free-electron laser (SwissFEL) was used to generate large-bandwidth X-ray pulses [Δλ/λ = 2.2% full width at half-maximum (FWHM)], which were applied in SFX with the aim of improving the partiality of Bragg spots and thus decreasing sample consumption while maintaining the data quality. Sensitive data-quality indicators such as anomalous signal from native thaumatin micro-crystals and de novo phasing results were used to quantify the benefits of using pink X-ray pulses to obtain accurate structure-factor amplitudes. Compared with data measured using the same setup but using X-ray pulses with typical quasi-monochromatic XFEL bandwidth (Δλ/λ = 0.17% FWHM), up to fourfold reduction in the number of indexed diffraction patterns required to obtain similar data quality was achieved. This novel approach, pink-beam SFX, facilitates the yet underutilized de novo structure determination of challenging proteins at XFELs, thereby opening the door to more scientific breakthroughs.

Purinergic GPCR transmembrane residues involved in ligand recognition and dimerization

We compare the GPCR-ligand interactions and highlight important residues for recognition in purinergic receptors-from both X-ray crystallographic and cryo-EM structures. These include A1 and A2A adenosine receptors, and P2Y1 and P2Y12 receptors that respond to ADP and other nucleotides. These receptors are important drug discovery targets for immune, metabolic and nervous system disorders. In most cases, orthosteric ligands are represented, except for one allosteric P2Y1 antagonist. This review catalogs the residues and regions that engage in contacts with ligands or with other GPCR protomers in dimeric forms. Residues that are in proximity to bound ligands within purinergic GPCR families are correlated. There is extensive conservation of recognition motifs between adenosine receptors, but the P2Y1 and P2Y12 receptors are each structurally distinct in their ligand recognition. Identifying common interaction features for ligand recognition within a receptor class that has multiple structures available can aid in the drug discovery process.

Pathways and Mechanism of Caffeine Binding to Human Adenosine A2A Receptor

Caffeine (CFF) is a common antagonist to the four subtypes of adenosine G-protein-coupled receptors (GPCRs), which are critical drug targets for treating heart failure, cancer, and neurological diseases. However, the pathways and mechanism of CFF binding to the target receptors remain unclear. In this study, we have performed all-atom-enhanced sampling simulations using a robust Gaussian-accelerated molecular dynamics (GaMD) method to elucidate the binding mechanism of CFF to human adenosine A_2A receptor (A_2AAR). Multiple 500–1,000 ns GaMD simulations captured both binding and dissociation of CFF in the A_2AAR. The GaMD-predicted binding poses of CFF were highly consistent with the x-ray crystal conformations with a characteristic hydrogen bond formed between CFF and residue N6.55 in the receptor. In addition, a low-energy intermediate binding conformation was revealed for CFF at the receptor extracellular mouth between ECL2 and TM1. While the ligand-binding pathways of the A_2AAR were found similar to those of other class A GPCRs identified from previous studies, the ECL2 with high sequence divergence serves as an attractive target site for designing allosteric modulators as selective drugs of the A_2AAR.

DNA-Mediated Stack Formation of Nanodiscs

Membrane-scaffolding proteins (MSPs) derived from apolipoprotein A-1 have become a versatile tool in generating nano-sized discoidal membrane mimetics (nanodiscs) for membrane protein research. Recent efforts have aimed at exploiting their controlled lipid protein ratio and size distribution to arrange membrane proteins in regular supramolecular structures for diffraction studies. Thereby, direct membrane protein crystallization, which has remained the limiting factor in structure determination of membrane proteins, would be circumvented. We describe here the formation of multimers of membrane-scaffolding protein MSP1D1-bounded nanodiscs using the thiol reactivity of engineered cysteines. The mutated positions N42 and K163 in MSP1D1 were chosen to support chemical modification as evidenced by fluorescent labeling with pyrene. Minimal interference with the nanodisc formation and structure was demonstrated by circular dichroism spectroscopy, differential light scattering and size exclusion chromatography. The direct disulphide bond formation of nanodiscs formed by the MSP1D1_N42C variant led to dimers and trimers with low yield. In contrast, transmission electron microscopy revealed that the attachment of oligonucleotides to the engineered cysteines of MSP1D1 allowed the growth of submicron-sized tracts of stacked nanodiscs through the hybridization of nanodisc populations carrying complementary strands and a flexible spacer.

The Specificity of Downstream Signaling for A1 and A2AR Does Not Depend on the C-Terminus, Despite the Importance of This Domain in Downstream Signaling Strength

Recent efforts to determine the high-resolution crystal structures for the adenosine receptors (A₁R and A_2AR) have utilized modifications to the native receptors in order to facilitate receptor crystallization and structure determination. One common modification is a truncation of the unstructured C-terminus, which has been utilized for all the adenosine receptor crystal structures obtained to date. Ligand binding for this truncated receptor has been shown to be similar to full-length receptor for A_2AR. However, the C-terminus has been identified as a location for protein-protein interactions that may be critical for the physiological function of these important drug targets. We show that variants with A_2AR C-terminal truncations lacked cAMP-linked signaling compared to the full-length receptor constructs transfected into mammalian cells (HEK-293). In addition, we show that in a humanized yeast system, the absence of the full-length C-terminus affected downstream signaling using a yeast MAPK response-based fluorescence assay, though full-length receptors showed native-like G-protein coupling. To further study the G protein coupling, we used this humanized yeast platform to explore coupling to human-yeast G-protein chimeras in a cellular context. Although the C-terminus was essential for Gα protein-associated signaling, chimeras of A₁R with a C-terminus of A_2AR coupled to the A₁R-specific Gα (i.e., Gαi1 versus Gαs). This surprising result suggests that the C-terminus is important in the signaling strength, but not specificity, of the Gα protein interaction. This result has further implications in drug discovery, both in enabling the experimental use of chimeras for ligand design, and in the cautious interpretation of structure-based drug design using truncated receptors.