Integrative residue-intuitive machine learning and MD Approach to Unveil Allosteric Site and Mechanism for ß2AR.

Allosteric drugs offer a new avenue for modern drug design. However, the identification of cryptic allosteric sites presents a formidable challenge. Following the allostery nature of residue-driven conformation transition, we propose a state-of-the-art computational pipeline by developing a residue-intuitive hybrid machine learning (RHML) model coupled with molecular dynamics (MD) simulation, through which we can efficiently identify the allosteric site and allosteric modulator as well as reveal their regulation mechanism. For the clinical target ß2-adrenoceptor (ß2AR), we discover an additional allosteric site located around residues D792.50, F2826.44, N3187.45 and S3197.46 and one putative allosteric modulator ZINC5042. Using Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) and protein structure network (PSN), the allosteric potency and regulation mechanism are probed to further improve identification accuracy. Benefiting from sufficient computational evidence, the experimental assays then validate our predicted allosteric site, negative allosteric potency and regulation pathway, showcasing the effectiveness of the identification pipeline in practice. We expect that it will be applicable to other target proteins.

Target-specific affinity separation of the bioactive compounds from herbal extract using the spin column packed with the immobilized protein microspheres prior to LC-MS analysis.

Excellent pretreatments before instrumental analysis are critical for separation and determination of target compounds for discovery of new drugs from herb medicines. We developed a rapid and highly-selective method to separate the bioactive compounds from herbal extract using protein affinity-selection spin column, which was packed with the new sorbent materials from integrating the recombinant ß2-adrenoceptor (ß2-AR) directly out of cell lysates onto the surface of microspheres. Protein affinity-selection spin column was placed in a centrifugal tube, where after the non-specific binders were released to the filtrate under the operational centrifugation, the specific binders on the spin column were cleaned with a washing solvent for LC-MS analysis. The known agonists of ß2-AR were retained/released on protein affinity-selection spin column but not on control column, demonstrating the method with good recovery (79.4∼95.7 %) and high repeatability (RSD < 3.5 %). The adsorption features of three ligands on the spin column were described best by Prism saturation binding model, and the high-affinity binding and the large binding capacity of the spin column make it feasible to trap the trace analytes effectively. It was applied in separating bioactive compounds from Alstoniae Scholaris extract, two of which were identified as picrinine and oleanolic acid in combination with LC-MS and verified as the potential agonists towards ß2-AR though molecular docking and cell experiments. Our study demonstrated that, the spin column with the immobilized protein sorbents in the centrifugal filter device represents a promising tool, enabling rapid and target-specific affinity separation of the bioactive compounds from herbal extract.

Bitopic Ligands Support the Presence of a Metastable Binding Site at the ß2 Adrenergic Receptor.

Metastable binding sites (MBS) have been observed in a multitude of molecular dynamics simulations and can be considered low affinity allosteric binding sites (ABS) that function as stepping stones as the ligand moves toward the orthosteric binding site (OBS). Herein, we show that MBS can be utilized as ABS in ligand design, resulting in ligands with improved binding kinetics. Four homobivalent bitopic ligands (1-4) were designed by molecular docking of (S)-alprenolol ((S)-ALP) in the cocrystal structure of the ß2 adrenergic receptor (ß2AR) bound to the antagonist ALP. Ligand 4 displayed a potency and affinity similar to (S)-ALP, but with a >4-fold increase in residence time. The proposed binding mode was confirmed by X-ray crystallography of ligand 4 in complex with the ß2AR. This ligand design principle can find applications beyond the ß2AR and G protein-coupled receptors (GPCRs) as a general approach for improving the pharmacological profile of orthosteric ligands by targeting the OBS and an MBS simultaneously.

Vibrational Energy Landscapes and Energy Flow in GPCRs.

We construct and analyze disconnectivity graphs to provide the first graphical representation of the vibrational energy landscape of a protein, in this study ß2AR, a G-protein coupled receptor (GPCR), in active and inactive states. The graphs, which indicate the relative free energy of each residue and the minimum free energy barriers for energy transfer between them, reveal important composition, structural and dynamic properties that mediate the flow of energy. Prolines and glycines, which contribute to GPCR plasticity and function, are identified as bottlenecks to energy transport along the backbone from which alternative pathways for energy transport via nearby noncovalent contacts emerge, seen also in the analysis of first passage time (FPT) distributions presented here. Striking differences between the disconnectivity graphs and FPT distributions for the inactive and active states of ß2AR are found where structural and dynamic changes occur upon activation, contributing to allosteric regulation.

Improved Highly Mobile Membrane Mimetic Model for Investigating Protein-Cholesterol Interactions.

Cholesterol (CHL) plays an integral role in modulating the function and activity of various mammalian membrane proteins. Due to the slow dynamics of lipids, conventional computational studies of protein-CHL interactions rely on either long-time scale atomistic simulations or coarse-grained approximations to sample the process. A highly mobile membrane mimetic (HMMM) has been developed to enhance lipid diffusion and thus used to facilitate the investigation of lipid interactions with peripheral membrane proteins and, with customized in silico solvents to replace phospholipid tails, with integral membrane proteins. Here, we report an updated HMMM model that is able to include CHL, a nonphospholipid component of the membrane, henceforth called HMMM-CHL. To this end, we had to optimize the effect of the customized solvents on CHL behavior in the membrane. Furthermore, the new solvent is compatible with simulations using force-based switching protocols. In the HMMM-CHL, both improved CHL dynamics and accelerated lipid diffusion are integrated. To test the updated model, we have applied it to the characterization of protein-CHL interactions in two membrane protein systems, the human ß2-adrenergic receptor (ß2AR) and the mitochondrial voltage-dependent anion channel 1 (VDAC-1). Our HMMM-CHL simulations successfully identified CHL binding sites and captured detailed CHL interactions in excellent consistency with experimental data as well as other simulation results, indicating the utility of the improved model in applications where an enhanced sampling of protein-CHL interactions is desired.

Structure optimization of Cmpd-15 as negative allosteric modulators for the ß2-adrenergic receptor.

19 derivatives of 1-benzyl-3-arylpyrazole-5-carboxamides (H1-H19) and 5 derivatives of 1-benzyl-5-arylpyrazole-3-carboxamides (J1-J5) have been designed and synthesized as potential negative allosteric modulators (NAMs) for the ß2-adrenergic receptor (ß2AR). The new pyrazole derivatives were screened on the classic G-protein dependent signaling pathway at ß2AR. The majority of 1-benzyl-3-aryl-pyrazole-5-carboxamide derivatives show more potent allosteric antagonistic activity against ß2AR than Cmpd-15, the first reported ß2AR NAM. However, the 1-benzyl-5-arylpyrazole-3-carboxamide derivatives exhibit very poor or even no allosteric antagonistic activity for ß2AR. Furthermore, the active pyrazole derivatives have relative better drug-like profiles than Cmpd-15. Taken together, we discovered a series of derivatives of 1-benzyl-3-arylpyrazole-5-carboxamides as a novel scaffold of ß2AR NAM.

A Study on the Robustness and Stability of Explainable Deep Learning in an Imbalanced Setting: The Exploration of the Conformational Space of G Protein-Coupled Receptors.

G-protein coupled receptors (GPCRs) are transmembrane proteins that transmit signals from the extracellular environment to the inside of the cells. Their ability to adopt various conformational states, which influence their function, makes them crucial in pharmacoproteomic studies. While many drugs target specific GPCR states to exert their effects-thereby regulating the protein's activity-unraveling the activation pathway remains challenging due to the multitude of intermediate transformations occurring throughout this process, and intrinsically influencing the dynamics of the receptors. In this context, computational modeling, particularly molecular dynamics (MD) simulations, may offer valuable insights into the dynamics and energetics of GPCR transformations, especially when combined with machine learning (ML) methods and techniques for achieving model interpretability for knowledge generation. The current study builds upon previous work in which the layer relevance propagation (LRP) technique was employed to interpret the predictions in a multi-class classification problem concerning the conformational states of the ß2-adrenergic (ß2AR) receptor from MD simulations. Here, we address the challenges posed by class imbalance and extend previous analyses by evaluating the robustness and stability of deep learning (DL)-based predictions under different imbalance mitigation techniques. By meticulously evaluating explainability and imbalance strategies, we aim to produce reliable and robust insights.

Optical control of nanobody-mediated protein activity modulation with a photocleavable fluorescent protein.

Antibodies are crucial in various biological applications due to their specific binding to target molecules, altering protein function and structure. The advent of single-chain antibodies such as nanobodies has paved the way for broader applicability in both research and therapies due to their small size and efficient tissue penetration. Recently, several approaches have been reported to optically control the antigen-binding affinity of nanobodies. Here, we show an alternative strategy for creating photo-activatable nanobodies. By fusing the photocleavable protein PhoCl with the N-terminus of the nanobody (named optoNb60), we successfully demonstrated light-dependent restoration of the antigen-binding ability and the following modulation of the activity of a target protein, the beta-2 adrenergic receptor. Moreover, the activation of optoNb60 was monitored by the fluorescence changes upon photoconversion. The compatibility of the uncaging design with the previously reported optogenetic molecules using nanobodies will contribute to the further optimization of the response capabilities of existing optogenetic tools, thereby expanding their applicability.

Screening for promising multi-target bioactive components from Cortex Mori Radicis for the treatment of chronic cor pulmonale based on immobilized beta1-adrenergic receptor and beta2-adrenergic receptor chromatography.

Cortex Morin Radicis (CMR) is the dried root bark of Morus alba. L. It has a variety of effects such as antibacterial, anti-tumour, treatment of cardiovascular diseases or upper respiratory tract disease and so on. The pursuit for drugs selected from Cortex Mori Radicis having improved therapeutic efficacy necessitates increasing research on new assays for screening bioactive compounds with multi-targets. In this work, we applied immobilized ß1-AR and ß2-AR as the stationary phase in chromatographic column to screen bioactive compounds from Cortex Morin Radicis. Specific ligands of the two receptors (e.g. esmolol, metoprolol, atenolol, salbutamol, methoxyphenamine, tulobuterol and clorprenaline) were utilized to characterize the specificity and bioactivity of the columns. We used high performance affinity chromatography coupled with ESI-MS to screen targeted compounds of Cortex Morin Radicis. By zonal elution, we identified morin as a bioactive compound simultaneously binding to ß1-AR and ß2-AR. The compound exhibited the association constants of 3.10 × 104 and 2.60 × 104 M-1 on the ß1-AR and ß2-AR column. On these sites, the dissociation rate constants were calculated to be 0.131 and 0.097 s-1. Molecular docking indicated that the binding of morin to the two receptors occurred on Asp200, Asp121, and Val122 of ß1-AR, Asn312, Thr110, Asp113, Tyr316, Gly90, Phe193, Ile309, and Trp109 of ß2-AR. Likewise, mulberroside C was identified as the bioactive compound binding to ß2-AR. The association constants and dissociation rate constants were calculated to be 1.08 × 104 M-1 and 0.900 s-1. Molecular docking also indicated that mulberroside C could bind to ß2-AR receptor on its agonist site. Taking together, we demonstrated that the chromatographic strategy to identify bioactive natural products based on the ß1-AR and ß2-AR immobilization, has potential for screening bioactive compounds with multi-targets from complex matrices including traditional Chinese medicines.

A high-affinity, cis-on photoswitchable beta blocker to optically control ß2-adrenergic receptors in vitro and in vivo.

This study introduces (S)-Opto-prop-2, a second-generation photoswitchable ligand designed for precise modulation of ß2-adrenoceptor (ß2AR). Synthesised by incorporating an azobenzene moiety with propranolol, (S)-Opto-prop-2 exhibited a high PSScis (photostationary state for cis isomer) percentage (∼90 %) and a favourable half-life (>10 days), facilitating diverse bioassay measurements. In vitro, the cis-isomer displayed substantially higher ß2AR binding affinity than the trans-isomer (1000-fold), making (S)-Opto-prop-2 one of the best photoswitchable GPCR (G protein-coupled receptor) ligands reported so far. Molecular docking of (S)-Opto-prop-2 in the X-ray structure of propranolol-bound ß2AR followed by site-directed mutagenesis studies, identified D1133.32, N3127.39 and F2896.51 as crucial residues that contribute to ligand-receptor interactions at the molecular level. In vivo efficacy was assessed using a rabbit ocular hypertension model, revealing that the cis isomer mimicked propranolol's effects in reducing intraocular pressure, while the trans isomer was inactive. Dynamic optical modulation of ß2AR by (S)-Opto-prop-2 was demonstrated in two different cAMP bioassays and using live-cell confocal imaging, indicating reversible and dynamic control of ß2AR activity using the new photopharmacology tool. In conclusion, (S)-Opto-prop-2 emerges as a promising photoswitchable ligand for precise and reversible ß2AR modulation with light. The new tool shows superior cis-on binding affinity, one of the largest reported differences in affinity (1000-fold) between its two configurations, in vivo efficacy, and dynamic modulation. This study contributes valuable insights into the evolving field of photopharmacology, offering a potential avenue for targeted therapy in ß2AR-associated pathologies.

A universal method for surface-based binding assays by preparing immobilized ß2-adrenergic receptor stationary phase using solid binding peptide as a fusion tag.

Protein functionalized surface has the potential to develop new assays for determining the drug-like properties of potential compounds and discovering specific partners of G protein-coupled receptors (GPCRs). However, a universal method for purifying and immobilizing functional GPCRs has remained elusive. To this end, we developed a general and rapid way to purify and immobilize ß2-adrenergic receptor (ß2AR) by silicon-specific peptide. We screened CotB1p as a tag from six silica-binding peptides (minTBP-1, CotB1p, SB7, Car9, and Si4-1) by examining their affinity to macroporous silica gel. We investigated the adsorption and desorption of CotB1p-tagged ß2-adrenoceptor (ß2AR-CotB1p) under diverse conditions to propose a protocol for receptor purification and immobilization. Under optimized conditions, ß2AR immobilization were achieved by directly immersing cell lysates harboring the receptor with silica gel, and the elution of the receptor without demonstratable contaminants was realized by including l-arginine/L-lysine in the elutes. This allows purification of the receptor from Escherichia coli (E.coli) lysates with a purity of 95 %. The immobilized receptor was utilized as a stationary phase to reveal the tag impact on ligand-binding outputs by comparing the CotB1p-strategy with a typical covalent method. The KAs of salbutamol, chlorprenaline, tulobuterol, and terbutaline on ß2AR-CotB1p column were 1.26 × 106, 6.59 × 106, 7.90 × 106, and 8.97 × 105 M-1 respectively, which were two orders of magnitude higher than those on the Halo-ß2AR column. The whole immobilization was accomplished within 30 min without the requirement of any special treatment, resulting in enhanced accuracy for determining receptor-ligand binding parameters. Taken together, CotB1p-mediated strategy is simple, rapid, and universal for purification or immobilization of unstable biomolecules like GPCRs for analytical and biological applications.

Computational Modeling and Characterization of Peptides Derived from Nanobody Complementary-Determining Region 2 (CDR2) Targeting Active-State Conformation of the ß2-Adrenergic Receptor (ß2AR).

This study assessed the suitability of the complementarity-determining region 2 (CDR2) of the nanobody (Nb) as a template for the derivation of nanobody-derived peptides (NDPs) targeting active-state ß2-adrenergic receptor (ß2AR) conformation. Sequences of conformationally selective Nbs favoring the agonist-occupied ß2AR were initially analyzed by the informational spectrum method (ISM). The derived NDPs in complex with ß2AR were subjected to protein-peptide docking, molecular dynamics (MD) simulations, and metadynamics-based free-energy binding calculations. Computational analyses identified a 25-amino-acid-long CDR2-NDP of Nb71, designated P4, which exhibited the following binding free-energy for the formation of the ß2AR:P4 complex (ΔG = -6.8 ± 0.8 kcal/mol or a Ki = 16.5 µM at 310 K) and mapped the ß2AR:P4 amino acid interaction network. In vitro characterization showed that P4 (i) can cross the plasma membrane, (ii) reduces the maximum isoproterenol-induced cAMP level by approximately 40% and the isoproterenol potency by up to 20-fold at micromolar concentration, (iii) has a very low affinity to interact with unstimulated ß2AR in the cAMP assay, and (iv) cannot reduce the efficacy and potency of the isoproterenol-mediated ß2AR/ß-arrestin-2 interaction in the BRET2-based recruitment assay. In summary, the CDR2-NDP, P4, binds preferentially to agonist-activated ß2AR and disrupts Gαs-mediated signaling.

Can the heptapeptide ASSIVSF of the ß2-adrenoceptor recognize ephedrine and pseudoephedrine epimers in a complex system?

Epimer separation is crucial in the field of analytical chemistry, separation science, and the pharmaceutical industry. No reported methods could separate simultaneously epimers or even isomers and remove other unwanted, co-existing, interfering substances from complex systems like herbal extracts. Herein, we prepared a heptapeptide-modified stationary phase for the separation of 1R,2S-(-)-ephedrine [(-)-Ephe] and 1S,2S-(+)-pseudoephedrine [(+)-Pse] epimers from Ephedra sinica Stapf extract and blood samples. The heptapeptide stationary phase was comprehensively characterized by scanning electron microscopy, X-ray photoelectron spectroscopy, and Fourier transform infrared spectroscopy. The separation efficiency of the heptapeptide column was compared with an affinity column packed with full-length ß2-AR functionalized silica gel (ß2-AR column). The binding affinity of the heptapeptide with (+)-Pse was 3-fold greater than that with (-)-Ephe. Their binding mechanisms were extensively characterized by chromatographic analysis, ultraviolet spectra, circular dichroism analysis, isothermal titration calorimetry, and molecule docking. An enhanced hydrogen bonding was clearly observed in the heptapeptide-(+)-Pse complex. Such results demonstrated that the heptapeptide can recognize (+)-Pse and (-)-Ephe epimers in a complex system. This work, we believe, was the first report to simultaneously separate epimers and remove non-specific interfering substances from complex samples. The method was potentially applicable to more challenging sample separation, such as chiral separation from complex systems.

Time-resolved cryo-EM of G-protein activation by a GPCR.

G-protein-coupled receptors (GPCRs) activate heterotrimeric G proteins by stimulating guanine nucleotide exchange in the Gα subunit1. To visualize this mechanism, we developed a time-resolved cryo-EM approach that examines the progression of ensembles of pre-steady-state intermediates of a GPCR-G-protein complex. By monitoring the transitions of the stimulatory Gs protein in complex with the ß2-adrenergic receptor at short sequential time points after GTP addition, we identified the conformational trajectory underlying G-protein activation and functional dissociation from the receptor. Twenty structures generated from sequential overlapping particle subsets along this trajectory, compared to control structures, provide a high-resolution description of the order of main events driving G-protein activation in response to GTP binding. Structural changes propagate from the nucleotide-binding pocket and extend through the GTPase domain, enacting alterations to Gα switch regions and the α5 helix that weaken the G-protein-receptor interface. Molecular dynamics simulations with late structures in the cryo-EM trajectory support that enhanced ordering of GTP on closure of the α-helical domain against the nucleotide-bound Ras-homology domain correlates with α5 helix destabilization and eventual dissociation of the G protein from the GPCR. These findings also highlight the potential of time-resolved cryo-EM as a tool for mechanistic dissection of GPCR signalling events.

Probing Activation and Conformational Dynamics of the Vesicle-Reconstituted ß2 Adrenergic Receptor at the Single-Molecule Level.

G-protein-coupled receptors (GPCRs) are structurally flexible membrane proteins that mediate a host of physiological responses to extracellular ligands like hormones and neurotransmitters. Fine features of their dynamic structural behavior are hypothesized to encode the functional plasticity seen in GPCR activity, where ligands with different efficacies can direct the same receptor toward different signaling phenotypes. Although the number of GPCR crystal structures is increasing, the receptors are characterized by complex and poorly understood conformational landscapes. Therefore, we employed a fluorescence microscopy assay to monitor conformational dynamics of single ß2 adrenergic receptors (ß2ARs). To increase the biological relevance of our findings, we decided not to reconstitute the receptor in detergent micelles but rather lipid membranes as proteoliposomes. The conformational dynamics were monitored by changes in the intensity of an environmentally sensitive boron-dipyrromethene (BODIPY 493/503) fluorophore conjugated to an endogenous cysteine (located at the cytoplasmic end of the sixth transmembrane helix of the receptor). Using total internal reflection fluorescence microscopy (TIRFM) and a single small unilamellar liposome assay that we previously developed, we followed the real-time dynamic properties of hundreds of single ß2ARs reconstituted in a native-like environment─lipid membranes. Our results showed that ß2AR-BODIPY fluctuates between several states of different intensity on a time scale of seconds, compared to BODIPY-lipid conjugates that show almost entirely stable fluorescence emission in the absence and presence of the full agonist BI-167107. Agonist stimulation changes the ß2AR dynamics, increasing the population of states with higher intensities and prolonging their durations, consistent with bulk experiments. The transition density plot demonstrates that ß2AR-BODIPY, in the absence of the full agonist, interconverts between states of low and moderate intensity, while the full agonist renders transitions between moderate and high-intensity states more probable. This redistribution is consistent with a mechanism of conformational selection and is a promising first step toward characterizing the conformational dynamics of GPCRs embedded in a lipid bilayer.

Biased Signaling in Mutated Variants of ß2-Adrenergic Receptor: Insights from Molecular Dynamics Simulations.

The molecular basis of receptor bias in G protein-coupled receptors (GPCRs) caused by mutations that preferentially activate specific intracellular transducers over others remains poorly understood. Two experimentally identified biased variants of ß2-adrenergic receptors (ß2AR), a prototypical GPCR, are a triple mutant (T68F, Y132A, and Y219A) and a single mutant (Y219A); the former bias the receptor toward the ß-arrestin pathway by disfavoring G protein engagement, while the latter induces G protein signaling explicitly due to selection against GPCR kinases (GRKs) that phosphorylate the receptor as a prerequisite of ß-arrestin binding. Though rigorous characterizations have revealed functional implications of these mutations, the atomistic origin of the observed transducer selectivity is not clear. In this study, we investigated the allosteric mechanism of receptor bias in ß2AR using microseconds of all-atom Gaussian accelerated molecular dynamics (GaMD) simulations. Our observations reveal distinct rearrangements in transmembrane helices, intracellular loop 3, and critical residues R1313.50 and Y3267.53 in the conserved motifs D(E)RY and NPxxY for the mutant receptors, leading to their specific transducer interactions. Moreover, partial dissociation of G protein from the receptor core is observed in the simulations of the triple mutant in contrast to the single mutant and wild-type receptor. The reorganization of allosteric communications from the extracellular agonist BI-167107 to the intracellular receptor-transducer interfaces drives the conformational rearrangements responsible for receptor bias in the single and triple mutants. The molecular insights into receptor bias of ß2AR presented here could improve the understanding of biased signaling in GPCRs, potentially opening new avenues for designing novel therapeutics with fewer side-effects and superior efficacy.

"Photo-Adrenalines": Photoswitchable ß2 -Adrenergic Receptor Agonists as Molecular Probes for the Study of Spatiotemporal Adrenergic Signaling.

ß2 -adrenergic receptor (ß2 -AR) agonists are used for the treatment of asthma and chronic obstructive pulmonary disease, but also play a role in other complex disorders including cancer, diabetes and heart diseases. As the cellular and molecular mechanisms in various cells and tissues of the ß2 -AR remain vastly elusive, we developed tools for this investigation with high temporal and spatial resolution. Several photoswitchable ß2 -AR agonists with nanomolar activity were synthesized. The most potent agonist for ß2 -AR with reasonable switching is a one-digit nanomolar active, trans-on arylazopyrazole-based adrenaline derivative and comprises valuable photopharmacological properties for further biological studies with high structural accordance to the native ligand adrenaline.

State-dependent dynamics of extramembrane domains in the ß2 -adrenergic receptor.

G protein-coupled receptors (GPCRs) are membrane-bound signaling proteins that play an essential role in cellular signaling processes. Due to their intrinsic function of transmitting internal signals in response to external cues, these receptors are adapted to be highly dynamic in nature. The ß2 -adrenergic receptor (ß2 AR) is a representative member of the family that has been extensively analyzed in terms of its structure and activation. Although the structure of the transmembrane domain has been characterized in the different functional states of the receptor, the conformational dynamics of the extramembrane domains, especially the intrinsically disordered regions are still emerging. In this study, we analyze the state-dependent dynamics of extramembrane domains of ß2 AR using atomistic molecular dynamics simulations. We introduce a parameter, the residue excess dynamics that allows us to better quantify receptor dynamics. Using this measure, we show that the dynamics of the extramembrane domains are sensitive to the receptor state. Interestingly, the ligand-bound intermediate R ' state shows the maximal dynamics compared to either the active R*G or inactive R states. Ligand binding appears to be correlated with high residue excess dynamics that are dampened upon G protein coupling. The intracellular loop-3 (ICL3) domain has a tendency to flip towards the membrane upon ligand binding, which could contribute to receptor "priming." We highlight an important ICL1-helix-8 interplay that is broken in the ligand-bound state but is retained in the active state. Overall, our study highlights the importance of characterizing the functional dynamics of the GPCR loop domains.

Molecular determinants of ligand efficacy and potency in GPCR signaling.

Heterotrimeric guanine nucleotide-binding protein (G protein)-coupled receptors (GPCRs) bind to extracellular ligands and drugs and modulate intracellular responses through conformational changes. Despite their importance as drug targets, the molecular origins of pharmacological properties such as efficacy (maximum signaling response) and potency (the ligand concentration at half-maximal response) remain poorly understood for any ligand-receptor-signaling system. We used the prototypical adrenaline-ß2 adrenergic receptor-G protein system to reveal how specific receptor residues decode and translate the information encoded in a ligand to mediate a signaling response. We present a data science framework to integrate pharmacological and structural data to uncover structural changes and allosteric networks relevant for ligand pharmacology. These methods can be tailored to study any ligand-receptor-signaling system, and the principles open possibilities for designing orthosteric and allosteric compounds with defined signaling properties.

Coevolution-Driven Method for Efficiently Simulating Conformational Changes in Proteins Reveals Molecular Details of Ligand Effects in the ß2AR Receptor.

With the advent of AI-powered structure prediction, the scientific community is inching closer to solving protein folding. An unresolved enigma, however, is to accurately, reliably, and deterministically predict alternative conformational states that are crucial for the function of, e.g., transporters, receptors, or ion channels where conformational cycling is innately coupled to protein function. Accurately discovering and exploring all conformational states of membrane proteins has been challenging due to the need to retain atomistic detail while enhancing the sampling along interesting degrees of freedom. The challenges include but are not limited to finding which degrees of freedom are relevant, how to accelerate the sampling along them, and then quantifying the populations of each micro- and macrostate. In this work, we present a methodology that finds relevant degrees of freedom by combining evolution and physics through machine learning and apply it to the conformational sampling of the ß2 adrenergic receptor. In addition to predicting new conformations that are beyond the training set, we have computed free energy surfaces associated with the protein's conformational landscape. We then show that the methodology is able to quantitatively predict the effect of an array of ligands on the ß2 adrenergic receptor activation through the discovery of new metastable states not present in the training set. Lastly, we also stake out the structural determinants of activation and inactivation pathway signaling through different ligands and compare them to functional experiments to validate our methodology and potentially gain further insights into the activation mechanism of the ß2 adrenergic receptor.