Uncoupling the Structure-Activity Relationships of β2 Adrenergic Receptor Ligands from Membrane Binding

High-throughput kinetics in drug discovery

The importance of a drug's kinetic profile and interplay of structure-kinetic activity with PK/PD has long been appreciated in drug discovery. However, technical challenges have often limited detailed kinetic characterization of compounds to the latter stages of projects. This review highlights the advances that have been made in recent years in techniques, instrumentation, and data analysis to increase the throughput of detailed kinetic and mechanistic characterization, enabling its application earlier in the drug discovery process.

The contribution of NaV1.6 to the efficacy of voltage-gated sodium channel inhibitors in wild type and NaV1.6 gain-of-function (GOF) mouse seizure control

BACKGROUND AND PURPOSE

Inhibitors of voltage-gated sodium channels (NaVs) are important anti-epileptic drugs, but the contribution of specific channel isoforms is unknown since available inhibitors are non-selective. We aimed to create novel, isoform selective inhibitors of Nav channels as a means of informing the development of improved antiseizure drugs.

EXPERIMENTAL APPROACH

We created a series of compounds with diverse selectivity profiles enabling block of NaV1.6 alone or together with NaV1.2. These novel NaV inhibitors were evaluated for their ability to inhibit electrically evoked seizures in mice with a heterozygous gain-of-function mutation (N1768D/+) in Scn8a (encoding NaV1.6) and in wild-type mice.

KEY RESULTS

Pharmacologic inhibition of NaV1.6 in Scn8aN1768D/+ mice prevented seizures evoked by a 6-Hz shock. Inhibitors were also effective in a direct current maximal electroshock seizure assay in wild-type mice. NaV1.6 inhibition correlated with efficacy in both models, even without inhibition of other CNS NaV isoforms.

CONCLUSIONS AND IMPLICATIONS

Our data suggest NaV1.6 inhibition is a driver of efficacy for NaV inhibitor anti-seizure medicines. Sparing the NaV1.1 channels of inhibitory interneurons did not compromise efficacy. Selective NaV1.6 inhibitors may provide targeted therapies for human Scn8a developmental and epileptic encephalopathies and improved treatments for idiopathic epilepsies.

Allosteric modulation of α1β3γ2 GABAA receptors by farnesol through the neurosteroid sites

In mammalian cells, all-trans farnesol, a 15-carbon isoprenol, is a product of the mevalonate pathway. It is the natural substrate of alcohol dehydrogenase and a substrate for CYP2E1, two enzymes implicated in ethanol metabolism. Studies have shown that farnesol is present in the human brain and inhibits voltage-gated Ca2+ channels at much lower concentrations than ethanol. Here we show that farnesol modulates the activity of γ-aminobutyric acid type A receptors (GABAARs), some of which also mediate the sedative activity of ethanol. Electrophysiology experiments performed in HEK cells expressing human α1β3γ2 or α6β3γ2 GABAARs revealed that farnesol increased chloride currents through positive allosteric modulation of these receptors and showed dependence on both the alcoholic functional group of farnesol and the length of the alkyl chain for activity. In silico studies using long-timescale unbiased all-atom molecular dynamics (MD) simulations of the human α1β3γ2 GABAA receptors revealed that farnesol modulates the channel by directly binding to the transmembrane neurosteroid-binding site, after partitioning into the surrounding membrane and reaching the receptor by lateral diffusion. Channel activation by farnesol was further characterized by several structural and dynamic variables, such as global twisting of the receptor's extracellular domain, tilting of the transmembrane M2 helices, radius, cross-sectional area, hydration status, and electrostatic potential of the channel pore. Our results expand the pharmacological activities of farnesol to yet another class of ion channels implicated in neurotransmission, thus providing a novel path for understanding and treatment of diseases involving GABAA receptor dysfunction.

Rational design, synthesis, and pharmacological evaluation of a cohort of novel beta-adrenergic receptors ligands enables an assessment of structure-activity relationships

Biomedical applications of molecules that are able to modulate β-adrenergic signaling have become increasingly attractive over the last decade, revealing that β-adrenergic receptors (β-ARs) are key targets for a plethora of therapeutic interventions, including cancer. Despite successes in β-AR drug discovery, identification of β-AR ligands that are useful as selective chemical tools in pharmacological studies of the three β-AR subtypes, or lead compounds for drug development is still a highly challenging task. This is mainly due to the intrinsic plasticity of β-ARs as G protein-coupled receptors in conjunction with the requirement for functional receptor subtype selectivity, tissue specificity and minimal off-target effects. With the aim to provide insight into structure-activity relationships for the three β-AR subtypes, we have synthesized and obtained the pharmacological profile of a series of structurally diverse compounds (named MC) that were designed based on the aryloxy-propanolamine scaffold of SR59230A. Comparative analysis of their predicted binding mode within the active and inactive states of the receptors in combination with their pharmacological profile revealed key structural elements that control their activity as agonists or antagonists, in addition to clues about substituents that mediate selectivity for one receptor subtype over the others. We anticipate that these results will facilitate selective β-AR drug development efforts.

Membrane Lipids Are an Integral Part of Transmembrane Allosteric Sites in GPCRs: A Case Study of Cannabinoid CB1 Receptor Bound to a Negative Allosteric Modulator, ORG27569, and Analogs

A growing number of G-protein-coupled receptor (GPCR) structures reveal novel transmembrane lipid-exposed allosteric sites. Ligands must first partition into the surrounding membrane and take lipid paths to these sites. Remarkably, a significant part of the bound ligands appears exposed to the membrane lipids. The experimental structures do not usually account for the surrounding lipids, and their apparent contribution to ligand access and binding is often overlooked and poorly understood. Using classical and enhanced molecular dynamics simulations, we show that membrane lipids are critical in the access and binding of ORG27569 and its analogs at the transmembrane site of cannabinoid CB1 receptor. The observed differences in the binding affinity and cooperativity arise from the functional groups that interact primarily with lipids. Our results demonstrate the significance of incorporating membrane lipids as an integral component of transmembrane sites for accurate characterization, binding-affinity calculations, and lead optimization in drug discovery.

L-DOPA and Droxidopa: From Force Field Development to Molecular Docking into Human β2-Adrenergic Receptor

The increasing interest in the molecular mechanism of the binding of different agonists and antagonists to β₂-adrenergic receptor (β₂AR) inactive and active states has led us to investigate protein–ligand interactions using molecular docking calculations. To perform this study, the 3.2 Å X-ray crystal structure of the active conformation of human β₂AR in the complex with the endogenous agonist adrenaline has been used as a template for investigating the binding of two exogenous catecholamines to this adrenergic receptor. Here, we show the derivation of L-DOPA and Droxidopa OPLS all atom (AA) force field (FF) parameters via quantum mechanical (QM) calculations, molecular dynamics (MD) simulations in aqueous solutions of the two catecholamines and the molecular docking of both ligands into rigid and flexible β₂AR models. We observe that both ligands share with adrenaline similar experimentally observed binding anchor sites, which are constituted by Asp113/Asn312 and Ser203/Ser204/Ser207 side chains. Moreover, both L-DOPA and Droxidopa molecules exhibit binding affinities comparable to that predicted for adrenaline, which is in good agreement with previous experimental and computational results. L-DOPA and Droxidopa OPLS AA FFs have also been tested by performing MD simulations of these ligands docked into β₂AR proteins embedded in lipid membranes. Both hydrogen bonds and hydrophobic interaction networks observed over the 1 μs MD simulation are comparable with those derived from molecular docking calculations and MD simulations performed with the CHARMM FF.

Applications of Molecular Dynamics Simulation in Protein Study

Molecular Dynamics (MD) Simulations is increasingly used as a powerful tool to study protein structure-related questions. Starting from the early simulation study on the photoisomerization in rhodopsin in 1976, MD Simulations has been used to study protein function, protein stability, protein-protein interaction, enzymatic reactions and drug-protein interactions, and membrane proteins. In this review, we provide a brief review for the history of MD Simulations application and the current status of MD Simulations applications in protein studies.

Identifying a selective oligopeptide clamp in the gas phase

A new experimental scheme that can both identify selective oligopeptides and resolve their folding conformations is reported, which reveals that a “V” type pentapeptide can provide binding sites that precisely clamp four β2-blockers.

Relative Binding Free-Energy Calculations at Lipid-Exposed Sites: Deciphering Hot Spots

Relative binding free-energy (RBFE) calculations are experiencing resurgence in the computer-aided drug design of novel small molecules due to performance gains allowed by cutting-edge molecular mechanic force fields and computer hardware. Application of RBFE to soluble proteins is becoming a routine, while recent studies outline necessary steps to successfully apply RBFE at the orthosteric site of membrane-embedded G-protein-coupled receptors (GPCRs). In this work, we apply RBFE to a congeneric series of antagonists that bind to a lipid-exposed, extra-helical site of the P2Y1 receptor. We find promising performance of RBFE, such that it may be applied in a predictive manner on drug discovery programs targeting lipid-exposed sites. Further, by the application of the microkinetic model, binding at a lipid-exposed site can be split into (1) membrane partitioning of the drug molecule followed by (2) binding at the extra-helical site. We find that RBFE can be applied to calculate the free energy of each step, allowing the uncoupling of observed binding free energy from the influence of membrane affinity. This protocol may be used to identify binding hot spots at extra-helical sites and guide drug discovery programs toward optimizing intrinsic activity at the target.

Analysis of L-DOPA and droxidopa binding to human β2-adrenergic receptor

Over the last two decades, an increasing number of studies has been devoted to a deeper understanding of the molecular process involved in the binding of various agonists and antagonists to active and inactive conformations of β2-adrenergic receptor (β2AR). The 3.2 Å x-ray crystal structure of human β2AR active state in combination with the endogenous low affinity agonist adrenaline offers an ideal starting structure for studying the binding of various catecholamines to adrenergic receptors. We show that molecular docking of levodopa (L-DOPA) and droxidopa into rigid and flexible β2AR models leads for both ligands to binding anchor sites comparable to those experimentally reported for adrenaline, namely D113/N312 and S203/S204/S207 side chains. Both ligands have a hydrogen bond network that is extremely similar to those of noradrenaline and dopamine. Interestingly, redocking neutral and protonated versions of adrenaline to rigid and flexible β2AR models results in binding poses that are more energetically stable and distinct from the x-ray crystal structure. Similarly, lowest energy conformations of noradrenaline and dopamine generated by docking into flexible β2AR models had binding free energies lower than those of best poses in rigid receptor models. Furthermore, our findings show that L-DOPA and droxidopa molecules have binding affinities comparable to those predicted for adrenaline, noradrenaline, and dopamine, which are consistent with previous experimental and computational findings and supported by the molecular dynamics simulations of β2AR-ligand complexes performed here.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

A Retrospective Look at the Impact of Binding Site Environment on the Optimization of TRPA1 Antagonists

Transient receptor potential ankyrin 1 (TRPA1) antagonists have generated broad interest in the pharmaceutical industry for the treatment of both pain and asthma. Over the past decade, multiple antagonist classes have been reported in the literature with a wide range of structural diversity. Our own work has focused on the development of proline sulfonamide and hypoxanthine-based antagonists, two antagonist classes with distinct physicochemical properties and pharmacokinetic (PK) trends. Late in our discovery program, cryogenic electron microscopy (cryoEM) studies revealed two different antagonist binding sites: a membrane-exposed proline sulfonamide transmembrane site and an intracellular hypoxanthine site near the membrane interface. A retrospective look at the discovery program reveals how the different binding sites, and their location relative to the cell membrane, influenced the optimization trajectories and overall drug profiles of each antagonist class.

Membrane-facilitated receptor access and binding mechanisms of long-acting β2-adrenergic receptor (β2-AR) agonists

The drugs salmeterol, formoterol, and salbutamol constitute the frontline treatment for asthma and other chronic pulmonary diseases. These drugs activate the β2-adrenergic receptors (β2-AR), a class A G-protein-coupled receptor (GPCR) and differ significantly in their clinical onset and duration of actions. According to the "microkinetic model," the long duration of action of salmeterol and formoterol compared to salbutamol were attributed, at least in part, to their high lipophilicity and increased local concentrations in the membrane near the receptor. However, the structural and molecular bases of how the lipophilic drugs reach the binding site of the receptor from the surrounding membrane remain unknown. Using a variety of classical and enhanced molecular dynamics simulation techniques, we investigated the membrane partitioning characteristics, binding, and unbinding mechanisms of the ligands. The obtained results offer remarkable insight into the functional role of membrane lipids in the ligand association process. Strikingly, salmeterol entered the binding site from the bilayer through transmembrane helices 1 and 7. The entry was preceded by membrane-facilitated rearrangement and presentation of its phenyl-alkoxy-alkyl tail as a passkey to an access route gated by F193, a residue known critical for salmeterol's affinity. Formoterol's access is through the aqueous path shared by other β2-AR agents. We observed a novel secondary path for salbutamol that is distinct from its primary route. Our study offers a mechanistic description for the membrane-facilitated access and binding of ligands to β2-AR and establishes a groundwork for recognizing membrane lipids as an integral component in the molecular recognition process. Significance Statement The cell membrane's functional role behind the duration of action of long-acting β2-adrenergic receptor (β2-AR) agonists such as salmeterol has been a subject of debate for a long time. We investigated the binding and unbinding mechanisms of the three commonly used β2-AR agonists, salmeterol, formoterol, and salbutamol, using advanced simulation techniques. The obtained results offer unprecedented insights into the active role of membrane lipids in facilitating access and binding of the ligands, affecting the molecular recognition process and their pharmacology.

Ligand binding at the protein-lipid interface: strategic considerations for drug design

Many drug targets are embedded within the phospholipid bilayer of cellular membranes, including G protein-coupled receptors, ion channels, transporters and membrane-bound enzymes. Increasing evidence from biophysical and structural studies suggests that many small-molecule drugs commonly associate with these targets at binding sites at the protein-phospholipid interface. Without a direct path from bulk solvent to a binding site, a drug must first partition in the phospholipid membrane before interacting with the protein target. This membrane access mechanism necessarily affects the interpretation of potency data, structure-activity relationships, pharmacokinetics and physicochemical properties for drugs that target these sites. With an increasing number of small-molecule intramembrane binding sites revealed through X-ray crystallography and cryogenic electron microscopy, we suggest that ligand-lipid interactions likely play a larger role in small-molecule drug action than commonly appreciated. This Perspective introduces key concepts and drug design considerations to aid discovery teams operating within this target space, and discusses challenges and future opportunities in the field.

Insights into adenosine A2A receptor activation through cooperative modulation of agonist and allosteric lipid interactions

The activation process of G protein-coupled receptors (GPCRs) has been extensively studied, both experimentally and computationally. In particular, Molecular Dynamics (MD) simulations have proven useful in exploring GPCR conformational space. The typical behaviour of class A GPCRs, when subjected to unbiased MD simulations from their crystallized inactive state, is to fluctuate between inactive and intermediate(s) conformations, even with bound agonist. Fully active conformation(s) are rarely stabilized unless a G protein is also bound. Despite several crystal structures of the adenosine A2a receptor (A2aR) having been resolved in complex with co-crystallized agonists and Gs protein, its agonist-mediated activation process is still not completely understood. In order to thoroughly examine the conformational landscape of A2aR activation, we performed unbiased microsecond-length MD simulations in quadruplicate, starting from the inactive conformation either in apo or with bound agonists: endogenous adenosine or synthetic NECA, embedded in two homogeneous phospholipid membranes: 1,2-dioleoyl-sn-glycerol-3-phosphoglycerol (DOPG) or 1,2-dioleoyl-sn-glycerol-3-phosphocholine (DOPC). In DOPC with bound adenosine or NECA, we observe transition to an intermediate receptor conformation consistent with the known adenosine-bound crystal state. In apo state in DOPG, two different intermediate conformations are obtained. One is similar to that observed with bound adenosine in DOPC, while the other is closer to the active state but not yet fully active. Exclusively, in DOPG with bound adenosine or NECA, we reproducibly identify receptor conformations with fully active features, which are able to dock Gs protein. These different receptor conformations can be attributed to the action/absence of agonist and phospholipid-mediated allosteric effects on the intracellular side of the receptor.

Identification of CNS-Penetrant Aryl Sulfonamides as Isoform-Selective NaV1.6 Inhibitors with Efficacy in Mouse Models of Epilepsy

Nonselective antagonists of voltage-gated sodium (Na_V) channels have been long used for the treatment of epilepsies. The efficacy of these drugs is thought to be due to the block of sodium channels on excitatory neurons, primarily Na_V1.6 and Na_V1.2. However, these currently marketed drugs require high drug exposure and suffer from narrow therapeutic indices. Selective inhibition of Na_V1.6, while sparing Na_V1.1, is anticipated to provide a more effective and better tolerated treatment for epilepsies. In addition, block of Na_V1.2 may complement the anticonvulsant activity of Na_V1.6 inhibition. We discovered a novel series of aryl sulfonamides as CNS-penetrant, isoform-selective Na_V1.6 inhibitors, which also displayed potent block of Na_V1.2. Optimization focused on increasing selectivity over Na_V1.1, improving metabolic stability, reducing active efflux, and addressing a pregnane X-receptor liability. We obtained compounds 30-32, which produced potent anticonvulsant activity in mouse seizure models, including a direct current maximal electroshock seizure assay.

Drug Concentration Asymmetry in Tissues and Plasma for Small Molecule-Related Therapeutic Modalities

The well accepted "free drug hypothesis" for small-molecule drugs assumes that only the free (unbound) drug concentration at the therapeutic target can elicit a pharmacologic effect. Unbound (free) drug concentrations in plasma are readily measurable and are often used as surrogates for the drug concentrations at the site of pharmacologic action in pharmacokinetic-pharmacodynamic analysis and clinical dose projection in drug discovery. Furthermore, for permeable compounds at pharmacokinetic steady state, the free drug concentration in tissue is likely a close approximation of that in plasma; however, several factors can create and maintain disequilibrium between the free drug concentration in plasma and tissue, leading to free drug concentration asymmetry. These factors include drug uptake and extrusion mechanisms involving the uptake and efflux drug transporters, intracellular biotransformation of prodrugs, membrane receptor-mediated uptake of antibody-drug conjugates, pH gradients, unique distribution properties (covalent binders, nanoparticles), and local drug delivery (e.g., inhalation). The impact of these factors on the free drug concentrations in tissues can be represented by K _p,uu, the ratio of free drug concentration between tissue and plasma at steady state. This review focuses on situations in which free drug concentrations in tissues may differ from those in plasma (e.g., K _p,uu > or <1) and discusses the limitations of the surrogate approach of using plasma-free drug concentration to predict free drug concentrations in tissue. This is an important consideration for novel therapeutic modalities since systemic exposure as a driver of pharmacologic effects may provide limited value in guiding compound optimization, selection, and advancement. Ultimately, a deeper understanding of the relationship between free drug concentrations in plasma and tissues is needed.

Does the Lipid Bilayer Orchestrate Access and Binding of Ligands to Transmembrane Orthosteric/Allosteric Sites of G Protein-Coupled Receptors?

The ligand-binding sites of many G protein-coupled receptors (GPCRs) are situated around and deeply embedded within the central pocket formed by their seven transmembrane-spanning α-helical domains. Generally, these binding sites are assumed accessible to endogenous ligands from the aqueous phase. Recent advances in the structural biology of GPCRs, along with biophysical and computational studies, suggest that amphiphilic and lipophilic molecules may gain access to these receptors by first partitioning into the membrane and then reaching the binding site via lateral diffusion through the lipid bilayer. In addition, several crystal structures of class A and class B GPCRs bound to their ligands offer unprecedented details on the existence of lipid-facing allosteric binding sites outside the transmembrane helices that can only be reached via lipid pathways. The highly organized structure of the lipid bilayer may direct lipophilic or amphiphilic drugs to a specific depth within the bilayer, changing local concentration of the drug near the binding site and affecting its binding kinetics. Additionally, the constraints of the lipid bilayer, including its composition and biophysical properties, may play a critical role in "pre-organizing" ligand molecules in an optimal orientation and conformation to facilitate receptor binding. Despite its clear involvement in molecular recognition processes, the critical role of the membrane in binding ligands to lipid-exposed transmembrane binding sites remains poorly understood and warrants comprehensive investigation. Understanding the mechanistic basis of the structure-membrane interaction relationship of drugs will not only provide useful insights about receptor binding kinetics but will also enhance our ability to take advantage of the apparent membrane contributions when designing drugs that target transmembrane proteins with improved efficacy and safety. In this minireview, we summarize recent structural and computational studies on membrane contributions to binding processes, elucidating both lipid pathways of ligand access and binding mechanisms for several orthosteric and allosteric ligands of class A and class B GPCRs.

Cell-based assays and molecular dynamics analysis of a boron-containing agonist with different profiles of binding to human and guinea pig beta2 adrenoceptors

The design of beta2 adrenoceptor (β₂AR) agonists is attractive because of their wide-ranging applications in medicine, and the details of agonist interactions with β₂AR are interesting because it is considered a prototype for G-protein coupled receptors. Preclinical studies for agonist development have involved biological assays with guinea pigs due to a similar physiology to humans. Boron-containing Albuterol derivatives (BCADs) designed as bronchodilators have improved potency and efficacy compared with their boron-free precursor on guinea pig β₂ARs (gpβ₂ARs), and two of the BCADs (BR-AEA and boronterol) conserve these features on cells expressing human β₂ARs (hβ₂ARs). The aim of this study was to test the BCAD Politerol on gpβ₂ARs and hβ₂ARs in vitro and in silico. Politerol displayed higher potency and efficacy on gpβ₂AR than on hβ₂AR in experimental assays, possible explanations are provided based on molecular modeling, and molecular dynamics simulations of about 0.25 µs were performed for the free and bound states adding up to 2 µs in total. There were slight differences, particularly in the role of the boron atom, in the interactions of Politerol with gpβ₂ARs and hβ₂ARs, affecting movements of transmembrane domains 5-7, known to be pivotal in receptor activation. These findings could be instrumental in the design of compounds selective for hβ₂ARs.

Link between a high k on for drug binding and a fast clinical action: to be or not to be?

Review articles on binding kinetics essentially focus on drugs that dissociate slowly from their target since this is required for the successful treatment of many pathophysiological conditions. Recently, the therapeutic benefit of a high k on (i.e. the second order association rate constant) has also been linked to fast association and to a fast clinical action. Other studies, however, called this assertion into question since additional factors, like the dosing paradigm and the binding mechanism, are important as well. The still ongoing reticence about integrating binding kinetics in lead optimization programs motivated us to critically review the link between the drug's kinetic rate constants and their in vitro and in vivo target occupancy profile, with special focus on k on. The presented simulations tally with a positive link between a drug's effective/observed association rate (which is quite easy to determine in vitro) and the swiftness of its clinical action. On the other hand, the simulations show that the k on-concept should not be confounded with the effective association process since increasing this parameter only enhances the drug's in vitro and in vivo association under certain conditions: the binding mechanism should be suitable, rebinding (and thus the factors within the target's micro-environment that favour this mechanism) should not be too prominent and the dosage should not be kept in par with the drug's affinity. Otherwise, increasing k on could be ineffective or even be counter-productive.

Recent Trends and Applications of Molecular Modeling in GPCR⁻Ligand Recognition and Structure-Based Drug Design

G protein-coupled receptors represent the largest family of human membrane proteins and are modulated by a variety of drugs and endogenous ligands. Molecular modeling techniques, especially enhanced sampling methods, have provided significant insight into the mechanism of GPCR–ligand recognition. Notably, the crucial role of the membrane in the ligand-receptor association process has earned much attention. Additionally, docking, together with more accurate free energy calculation methods, is playing an important role in the design of novel compounds targeting GPCRs. Here, we summarize the recent progress in the computational studies focusing on the above issues. In the future, with continuous improvement in both computational hardware and algorithms, molecular modeling would serve as an indispensable tool in a wider scope of the research concerning GPCR–ligand recognition as well as drug design targeting GPCRs.

The Molecular Mechanism Underlying Ligand Binding to the Membrane-Embedded Site of a G-Protein-Coupled Receptor

The crystal structure of P2Y₁ receptor (P2Y₁R), a class A GPCR, revealed a special extra-helical site for its antagonist, BPTU, which locates in-between the membrane and the protein. However, due to the limitation of crystallization experiments, the membrane was mimicked by use of detergents, and the information related to the binding of BPTU to the receptor in the membrane environment is rather limited. In the present work, we conducted a total of ∼7.5 μs all-atom simulations in explicit solvent using conventional molecular dynamics and multiple enhanced sampling methods, with models of BPTU and a POPC bilayer, both in the absence and presence of P2Y₁R. Our simulations revealed that BPTU prefers partitioning into the interface of polar/lipophilic region of the lipid bilayer before associating with the receptor. Then, it interacts with the second extracellular loop of the receptor and reaches the binding site through the lipid-receptor interface. In addition, by use of funnel-metadynamics simulations which efficiently enhance the sampling of bound and unbound states, we provide a statistically accurate description of the underlying binding free energy landscape. The calculated absolute ligand-receptor binding affinity is in excellent agreement with the experimental data (Δ G_{b0_theo} = -11.5 kcal mol^-1, Δ G_{b0_exp}= -11.7 kcal mol^-1). Our study broadens the view of the current experimental/theoretical models and our understanding of the protein-ligand recognition mechanism in the lipid environment. The strategy used in this work is potentially applicable to investigate ligands association/dissociation with other membrane-embedded sites, allowing identification of compounds targeting membrane receptors of pharmacological interest.

Structural insights into positive and negative allosteric regulation of a G protein-coupled receptor through protein-lipid interactions

Lipids are becoming known as essential allosteric modulators of G protein-coupled receptor (GPCRs). However, how they exert their effects on GPCR conformation at the atomic level is still unclear. In light of recent experimental data, we have performed several long-timescale molecular dynamics (MD) simulations, totalling 24 μs, to rigorously map allosteric modulation and conformational changes in the β₂ adrenergic receptor (β2AR) that occur as a result of interactions with three different phospholipids. In particular, we identify different sequential mechanisms behind receptor activation and deactivation, respectively, mediated by specific lipid interactions with key receptor regions. We show that net negatively charged lipids stabilize an active-like state of β2AR that is able to dock G_sα protein. Clustering of anionic lipids around the receptor with local distortion of membrane thickness is also apparent. On the other hand, net-neutral zwitterionic lipids inactivate the receptor, generating either fully inactive or intermediate states, with kinetics depending on lipid headgroup charge distribution and hydrophobicity. These chemical differences alter membrane thickness and density, which differentially destabilize the β2AR active state through lateral compression effects.

GPCRs: What Can We Learn from Molecular Dynamics Simulations?

Advances in the structural biology of G-protein Coupled Receptors have resulted in a significant step forward in our understanding of how this important class of drug targets function at the molecular level. However, it has also become apparent that they are very dynamic molecules, and moreover, that the underlying dynamics is crucial in shaping the response to different ligands. Molecular dynamics simulations can provide unique insight into the dynamic properties of GPCRs in a way that is complementary to many experimental approaches. In this chapter, we describe progress in three distinct areas that are particularly difficult to study with other techniques: atomic level investigation of the conformational changes that occur when moving between the various states that GPCRs can exist in, the pathways that ligands adopt during binding/unbinding events and finally, the influence of lipids on the conformational dynamics of GPCRs.

Interactions between β2-Adrenoceptor Ligands and Membrane: Atomic-Level Insights from Magic-Angle Spinning NMR

To understand the relationship between structural properties of the β2-adrenoceptor ligands and their interactions with membranes, we have investigated the location and distribution of five β2 agonists with distinct clinical durations and onsets of action (indacaterol, two indacaterol analogues, salmeterol and formoterol) in monounsaturated model membranes using magic angle spinning NMR to measure these interactions through both ¹H nuclear Overhauser enhancement (NOE) and paramagnetic relaxation enhancement (PRE) techniques. The hydrophilic aromatic groups of all five β2 agonists show maximum distribution in the lipid/water interface, but distinct location and dynamic behavior were observed for the lipophilic aromatic rings. Our study elucidates at atomic level that the hydrophobicity and substitution geometry of lipophilic groups play important roles in compound-lipid interactions.

Pharmacological property optimization for allosteric ligands: A medicinal chemistry perspective

New strategies to potentially improve drug safety and efficacy emerge with allosteric programs. Biased allosteric modulators can be designed with high subtype selectivity and defined receptor signaling endpoints, however, selecting the most meaningful parameters for optimization can be perplexing. Historically, "potency hunting" at the expense of physicochemical and pharmacokinetic optimization has led to numerous tool compounds with excellent pharmacological properties but no path to drug development. Conversely, extensive physicochemical and pharmacokinetic screening with only post hoc bias and allosteric characterization has led to inefficacious compounds or compounds with on-target toxicities. This field is rapidly evolving with new mechanistic understanding, changes in terminology, and novel opportunities. The intent of this digest is to summarize current understanding and debates within the field. We aim to discuss, from a medicinal chemistry perspective, the parameter choices available to drive SAR.